Transient and asymmetric dust structures in the TeV-bright nova RS Oph revealed by spectropolarimetry

A long-standing question related to nova eruptions is how these eruptions might lead to dust formation, despite the ostensibly in-hospitable environment for dust within the hot, irradiated ejecta. In the novae of systems such as the symbiotic binary RS Ophiuchi (RS Oph), ejecta from the white dwarf collide with pre-existing circumstellar material fed by the wind from the red-giant companion, o ﬀ ering a particularly clear view of some nova shocks and any associated dust production. In this work, we use the spectropolari-metric monitoring of the recurrent nova RS Oph starting two days after its eruption in August 2021 to show that: 1) dust was present in the RS Oph system as early as two days into the 2021 eruption; 2) the spatial distribution of this early dust was asymmetric, with components both aligned with and perpendicular to the orbital plane of the binary; 3) between two and nine days after the start of the eruption, this early dust was gradually destroyed; and 4) dust was again created, aligned roughly with the orbital plane of the binary more than 80 days after the start of the outburst, most likely as a result of shocks that arose as the ejecta interacted with circumbinary material concentrated in the orbital plane. The modeling of X-rays and very-high-energy (GeV and TeV) emission from RS Oph days to months into the 2021 eruption suggests that collisions between the ejecta and the circumbinary material may have led to shock formation in two distinct regions: the polar regions perpendicular to the orbital plane, where collimated outﬂows have been observed after prior eruptions, and a circumbinary torus in the orbital plane. The observations described here indicate that dust formed in approximately the same two regions, supporting the connection between shocks and dust in novae and revealing a very early onset of asymmetry. The spectropolarimetric signatures of RS Oph in the ﬁrst week into the 2021 outburst indicate: 1) polarized ﬂux across the H α emission line and 2) the position angle orientation relative to the radio axis is similar to what is seen from the spectropolarimetric signatures of active galactic nuclei (AGNs).


Introduction
With its seven recorded outbursts, RS Oph is perhaps the beststudied recurrent nova.It belongs to a small group of objects in which a white dwarf (WD) accretes from the wind of its red giant companion and ignites under degenerate conditions with a sudden and major increase in luminosity when significant material has accumulated on its surface.The RS Oph binary system consists of a massive carbon-oxygen WD with M WD =1.2-1.4M ⊙ (where M WD is the mass of the WD, e.g., Mikołajewska & Shara 2017) and a M2 III red giant (Zamanov et al. 2018) in a 453.6 day orbit with an inclination of i = 49 • -52 • (Brandi et al. 2009).
Detailed observations during the previous outburst allowed us to identify asymmetric structures on the expanding material (e.g., O'Brien et al. 2006;Bode et al. 2007;Ribeiro et al. 2009;Rupen et al. 2008;Sokoloski et al. 2008;Orlando et al. 2009;Montez et al. 2022) whose origin is still a matter of debate: Are the asymmetries generated once the expanding ejecta interacts with the surrounding circumstellar material?;Are they the direct result of radio jets anchored in the accretion disk?
⋆ e-mail: ynikolov@nao-rozhen.org After the 2006 outburst, infrared observations taken as early as 3.8 days after the start of the outburst indicated the presence of dust, mostly the silicate type, at about 17 AU from the WD, suggesting that this dust was already present before the outburst and survived the harsh environment during it (Barry et al. 2008).Later in the outburst, dust was detected on days 208 after the 2006 eruption by Evans et al. (2007).Derdzinski et al. (2017) conducted a theoretical study of the possibility that dust can be created behind radiative shocks during nova outbursts and concluded that under certain conditions, silicate-type dust can be created as early as one or two days into the outburst.This hypothesis can be tested during the 2021 outburst of RS Oph.
The most recent outburst of RS Oph, on 2021 August 8.93 UT, was first reported by K. Geary 1 .Hereafter we adopt the time of eruption determined by Munari & Valisa (2021) as T 0 =2459435.00JD.We present spectropolarimetric observations of the recurrent nova RS Oph from day 2 to day 352 after the start of the recent 2021 outburst, aimed at detecting variable intrinsic linear polarization and the presence of dust.In Section 2 we present details of the observations and analysis.Sections

Observations and Data Analysis
We obtained spectropolarimetric observations of RS Oph on eleven nights, with nine of them after the 2021 outburst and two observations taken 26 and 16 months before, during quiescence.The spectropolarimetric observations were secured with the 2-Channel-Focal-Reducer Rozhen (FoReRo2)2 , similar to that described by Jockers et al. (2000), attached to the Cassegrain focus of the 2.0m RCC telescope of the Bulgarian Rozhen National Astronomical Observatory.A rotatable Super-Achromatic True Zero-Order Waveplate 53 retarder (APSAW-5) was added to FoReRo2.The observations consist of sets of eight polarized spectra.The resolving power is R ≈ 1100.A beam-swapping technique was used (Bagnulo et al. 2009) to minimize the instrumental polarization.The instrumental polarization was corrected using a standard star with zero degrees of polarization.The offset between the position angle in the celestial and instrumental polarization was corrected using strongly polarized standard stars.For the spectral observations we used two instruments -the FoReRo2 for the low-resolution observations and the ESPERO spectrograph for the high-resolution observations.The ESPERO spectrograph is a fibre-fed Echelle spectrograph mounted on the 2m telescope at Rozhen NAO (Bonev et al. 2017).It has a resolution of ∼30 000 and covers the range 3900Å -9000Å.The journal of observations is given in Table 1.The table contains the date (in format YYYY-MM-DD), UT, of the beginning of the observation, the number of days after the start of the 2021 outburst (the sign "-" indicates the days before the outburst), and the total exposure time.

Spectral evolution
On Figure 1 we present the spectral evolution of RS Oph in the period August -October 2021 based on the low-resolution data.The most prominent feature on the spectra is the Hα emission line.After the outburst, the equivalent width of the Hα line is increased, reaching ∼ −3000 Å on 2021 Sep, followed by a decrease.The complex structure of the line is not well defined due to the low resolution of the spectra.On 2021 Aug 10 -2 days after the outburst, a P Cyg profile with velocity up to ∼ -4260 km s −1 is detected.The measured velocity is in agreement with the value, reported by Mikolajewska et al. (2021).HeI emission lines at 5876 Å, 6678 Å and 7065 Å rapidly evolve after the outburst, becoming prominent on the spectra after 2 days.The OI lines at 7477 Å and 7773 Å respectively are welldefined, especially on the spectrum obtained on 2021 Aug 10.The most prominent Fe feature in the spectra is the FeII multiplet at 5316 Å, being more evident on the earliest spectrum.In addition, some coronal lines of [FeX]   On Figure 2 we present the spectral evolution of selected emission lines in the spectra of RS Oph based on high-resolution data.The observations were obtained on 2021 September 17 (40 days after the outburst) and on 2022 May 20 (285 days after the outburst).For comparison, we present also a spectrum obtained before the outburst, on 2019 July 17.The Balmer Hβ emission line observed on 2021 September 17 is complex, revealing a structure similar to that observed in Hα line (Zamanov et al. 2022).A sharp P Cyg profile is visible at the top of a strong and broad emission.The triple-peaked He I line has radial velocities of -235 km/s, -30 km/s and 130 km/s.The radial velocities of the peaks of the OI 8446 Å line are -190 km/s, -25 km/s and 137 km/s.The He II lines are produced by shocked ejected material during its cooling phase, so the presence of the He II 4686 Å line 40 days after the outburst confirms the shock wave (Azzollini et al. 2023).The spectral observations reveal a typical spectral behaviour of a nova with a complex structure dominated by broad Balmer, He, O, and Fe lines with some P Cyg profiles, evolving towards a supersoft source phase.

Interstellar polarization toward the Recurrent Nova RS Oph
The observed polarization is a vectorial sum of intrinsic polarization and interstellar polarization.The recurrent nova phenomenon gives us the opportunity to study an object before and after the nova outburst.A 3D map of the interstellar polarization in the field around RS Oph is presented in Figure 3.The data are taken from Heiles (2000), Cropper (1990) and Nikolov et al. (2019).The observed degree of polarization and position angle of RS Oph during quiescence is similar to the degree of polarization and position angle of the stars of the vicinity of RS Oph published by Cropper (1990).Moreover, the spectropolarimetric observations during quiescence do not show the characteristics of an object with an intrinsic polarization, such as a depolarization effect across the H α emission line, a variable degree of polarization and a variable position angle (Nikolov et al. 2019).
On Figure 4 we show the observed degree of polarization and position angle of RS Oph taken 26 and 16 months before the 2021 outburst, during quiescence, and position angle and degree of polarization of 2MASS 17500625-0643499, a nearby source at an angular distance of ≈ 2.2 arcmin.The wavelength dependence of the observed degree of polarization of 2MASS 17500625-0643499 is the same as the wavelength dependence of RS Oph during quiescence at the same position angle.The observations of RS Oph during quiescence thus indicate that the observed pre-eruption polarization is interstellar, as no intrinsic polarization was detected (Nikolov et al. 2019;Somero et al. 2017).One of the methods to separate intrinsic and interstellar polarization is based on the depolarization effect of the emission lines.Using this method, Kawabata et al. (2006) determined the intrinsic polarization of Nova V475 Scuti, in which a clear depolarization effect was well visible on day 25 after the outburst.Our observations demonstrate that in the case of RS Oph, this method is applicable shortly after the outburst when the scattering region is relatively smaller than the emission line region and the degree of polarization in the Hα emission line is approximately 0. The He I 5876 and Hα lines have different degrees of polarization in the core of their profiles.Using them for the determination of the interstellar polarization would give a different (smaller) power law index.This method gives accurate results when the scattering region is much smaller than the emission-line region.In the case of RS Oph, the scattering region gradually expanded into the emission-line region, and the depolarization effect decreased over time.

Intrinsic polarization of RS Oph after 2021 outburst
To calculate the Stokes parameters of the intrinsic polarization Q int (λ) and U int (λ) of RS Oph, we derived these values from our observations and corrected them for the interstellar polarization using the equations: where Q obs (λ) and U obs (λ) are the Stokes parameters of the observed polarization after the 2021 outburst.Q IS (λ) and U IS (λ) represent the Stokes parameters of the interstellar polarization.
As interstellar Stokes parameters, we used the parameters Q IS (λ) and U IS (λ) obtained on 2019 July 2. The degree of intrinsic polarization (P int (λ)) and intrinsic position angle (θ int (λ)) were calculated using the following equations: where Θ 0 depends on the sign of Q int (λ) and U int (λ) (Bagnulo et al. 2009).
In Figure 5 we present the normalized intensity (normalized Stokes I) parameter and the degree of intrinsic polarization as a function of wavelength for each of the spectra taken during our campaign.The polarization spectra on days 34.8 (2021 September 12) and 35.8 (2021 September 13) are noisy because of the low continuum signal.The spectra obtained on days 34.8, 35.8 and 81.8 after the 2021 outburst are thus binned with a dispersion of 20 Å/pix.The maximum degree of intrinsic polarization, of about 1%, is observed on day 5.9.
Strong depolarization effects are observed in FeII 5169 Å, HeI 5876 Å, Hα, HeI 6678 Å, HeI 7065 Å and OI 7771 Å emission lines until 6.7 days after the outburst.The small polarization of the Hα and HeI 5876 Å emission lines on days 2, 5.9, and 6.7 is noticeable.The polarization in the Hα emission line is approximately 0.05% on day 2 and 0.3% on days 5.9 and 6.7.The HeI 5876 Å emission line polarization is approximately 0.35% on day 2 and 0.5%-0.4% on days 5.9 and 6.7, respectively.The typical error on the degree of polarization in the emission lines is 0.05%.The depolarization in the emission lines suggests that the scattering material is located within the emission line region.On the other hand, the polarization in the HeI 5876 Å emission line is higher than the Hα emission line.This is because the Hα region is more extended than HeI 5876 Å region.It is worth noting that the depolarization effect decreased with time.This can be explained by the scattering region expanding into the emissionline region.
Figure 6 shows the position angle of the intrinsic polarization as a function of wavelength from the observations taken after the 2021 outburst (days 2, 5.9, 6.7, 81.8 and 352).We omit data on days 34.8, 35.8, 300 and 301 due to their low signal-to-noise ratios.The data presented in Figure 6 were smoothed using a moving average of over seven points to better illustrate the position angle structures across the H α emission line.The position angle of the intrinsic polarization is approximately 90 • throughout the wavelength range explored, with the exception of the region around Hα on day 2, which shows a position angle of 180 • (see discussion in Section 3.3.3).The 90 • position angle, which implies a scattering structure that extends predominantly in the N-S direction on the sky, is aligned with the highly collimated outflows detected during the outbursts in 2006 (e.g.Rupen et al. 2008) and 2021 (Munari et al. 2022).

Wavelength dependence of the intrinsic polarization
The degree of intrinsic polarization after the 2021 outburst increases toward shorter wavelengths.We use a monotonic function (power-law) as a description of the variation of the degree of intrinsic polarization with the wavelength.We fit the continuum data with a power law as: where c is a scaling constant and β is a power law index.This index indicates the nature of the scattering process: i) Mie scattering by dust grains or, ii) Thompson scattering by electrons.The first process is wavelength dependent whereas Thompson scattering is not.Figure 7 shows the power-law index as a function of time after the outburst.The index increases from β = −2.2 on day 2 and reaches β = 0, which corresponds to Thompson scattering, on day 9.After that, the index decreased again until day 352, reaching values of β = −1.6.

Long-term variability of the degree of intrinsic polarization and position angle in broad-bands
We used the IRAF SBAND procedure to calculate broadband linear polarization.Two sets of bands were chosen -standard V band and filter C. We selected a region with a weak depolarization effect in the wavelength range between λ 1 = 5200Å and λ 2 = 5600Å to define the filter C. The transmission function of the C filter is 100 %.
Long-term variability of the degree of intrinsic polarization and position angle are presented in Figure 8. P int (V) and P int (C) represent the degree of intrinsic polarization in the V band and in the C band, respectively.The intrinsic polarization increased from day 2 to day 5.9 after the outburst.After day 5.9 the degree of intrinsic polarization decreased, and we obtain P int (V) = 0.34% ± 0.05% and P int (C) = 0.37% ± 0.05% on day 352 after the outburst.On day 5.9, we detected the greatest degree of intrinsic polarization, with P int (V) = 1.2% ± 0.05% and P int (C) = 1.39% ± 0.05%.The position angle of the intrinsic  polarization (P.A. int (V) and P.A. int (C)) was practically constant during the observations and had a value of approximately 90 • .

Evolution of the polarized flux.
In Figure 9, we present the position angle, the evolution of the polarized flux, and the degree of the intrinsic polarization during the first week of the outburst.The position angle (P.A.) of the intrinsic polarization on the 2nd day into the outburst (left column, top) shows a complex structure, with P.A. ≈ 90 • at wavelengths between λ = 6200 Å and λ = 6600 Å and P.A. ≈ 180 • for wavelengths between λ = 6600 Å and λ = 6650 Å. Radio interferometric imaging of RS Oph obtained 34 days after the start of the 2021 nova outburst shows a bipolar structure (Munari et al. 2022) in the east-west direction; the orbital plane is perpendicular to these radio lobes.When the P.A. was 90 • , it was thus aligned with the observed bipolar structure.When the P.A. was 180 • , it was aligned with the orbital plane.The wavelength range with P.A. = 180 • corresponds to the small peak in the polarized flux spectrum near a velocity of v = 3000 km s −1 (see the left column in Figure 9).The two main observed values of the position angle correspond to the mutually perpendicular planes.

Quick dust formation after the outburst
Our observations show the presence of dust as early as 2 days after the outburst.The wavelength dependence of the intrinsic polarization gives information about the nature of the scattering process.If the scattering is dominated by dust, we expect an increase in the intrinsic polarization toward shorter wavelengths.By modelling the continuum of the polarized spectra, we found that it was strongly wavelength-dependent from day 2 until day 9 after the outburst.This finding provides the most direct observational support of the presence of dust in a nova at such early times.
Although other novae have shown optical polarization within a few days of the start of the eruption, the case for early dust associated with the nova has been less definitive.As early as 1.36 days after the outburst of T Pyx, Pavana et al. (2019) detected polarization that could be attributed to asymmetries during the outburst or the presence of silicate dust that survived the eruption; however, given the unknown amount of interstellar polarization toward T Pyx, it was not possible for those authors to determine the wavelength dependence of the intrinsic polarization and how the intrinsic polarization and position an- gle varied with time, precluding a definitive conclusion about its origin.Kawabata et al. (2000) observed Nova V4444 Sgr from days 2 to 10 after the start of its 1999 outburst and found a wavelength-dependent intrinsic polarization with a power law index in the range: −2 ≲ β ≲ −1.Given that this dependence did not change with time, however, the authors concluded that the most likely origin was scattering by pre-existing dust located at several dozen AU from the central star.Nova V1494 Aql was observed between days 2.39 and 3.39 after the start of its outburst by Kawabata et al. (2001).In their Fig. 4, the visual negative slope of the degree of intrinsic polarization with the wavelength increases from day 2.39 to day 3.39 and decreases since then.
Evidence for dust formation in novae has been found through observations in the IR and optical at a much later phase of the outburst.Infrared observations of RS Oph with the Spitzer telescope from 208 to 430 days after the start of the 2006 eruption of RS Oph revealed the presence of silicate dust in the environs of the system (see Evans et al. 2007;Woodward et al. 2008).Evans et al. (2002) found that the wavelength-dependence of polarization in the case of V4362 Sgr suggests scattering by small dust grains (80 days after the start of the outburst).In the case of V339 Del, Kawakita et al. (2019)  A. was practically constant throughout our observing campaign, while the degree of intrinsic polarization reached a maximum at day 5.9 and has been decreasing since then.
neighbouring emission-free regions as values for the interstellar polarization.Nova V339 Del produced dust grains around 35 days into the outburst (Evans et al. 2017).
Our detection of wavelength-dependent intrinsic polarization indicates that, most likely, dust particles were created in dense regions shielded from the strong radiation field of the eruption.Such wavelength-dependent polarization was not observed during quiescence (down to our detection error of 0.05%), which strongly suggests that dust was not present during the quiescence phase.Moreover, Cropper (1990) found that the observed polarization of RS Oph a year after the start of the 1985 outburst was comparable with the polarization of nearby stars.
The possibility that dust can be created in radiative shock regions was studied by Derdzinski et al. (2017), who concluded that dust can be formed in, dense, shielded shock regions.These dense shocks, however, have been modelled by Metzger et al. (2014Metzger et al. ( , 2015) ) in the context of novae in CVs; similar studies tailored for embedded nova such as RS Oph are necessary to establish if dense regions behind shocks in those systems are also expected to be conducive to dust formation.Our data show that dust was present between days 2 and 9, and at about this same time, the γ-ray emission detected by Fermi peaked.The GeV-TeV emission has been modelled as due to multiple shocks, with the "slow" component first expanding into a high-density environment (Diesing et al. 2022).Our finding thus provides strong observational support that the scenario proposed by Derdzinski et al. (2017) is plausible for embedded novae as well as classical novae (see, e.g., Chomiuk et al. 2021;Babul et al. 2022, for a discussion of dust formation in classical novae).
The changes in the wavelength dependence of the intrinsic polarization with time, parameterized by the power law index β, indicate that the process of scattering switched from a mix of dust with different-sized grains to Thompson scattering after day 9.This switch indicates that any dust that was formed within the first week or so of the eruption quickly evaporated.Conditions then became favourable for dust formation again much later in the outburst.In their infrared observation on day 12 after the 2021 outburst, Woodward et al. (2021) did not find evidence of dust.The spectropolarimetric observations on days 300 and 352 after the outburst do show the presence of dust grains.Further infrared observations can confirm the presence of dust months to years after the 2021 outburst.

Early development of asymmetries in the ejecta
After the previous outburst in 2006, observations at a variety of wavelengths revealed evidence of the presence of asymmetric structures aligned preferentially in the East-West direction.Chesneau et al. (2007) found a structure moving rapidly in the E-W direction 5.5 days into the outburst.Bode et al. (2007) and Ribeiro et al. (2009)  Our detection of linear polarization on day 2 constitutes the earliest evidence of an asymmetric structure in RS Oph reported to date.With a position angle of 90 • aligned with the previously known E-W structures, our finding indicates that some structure was present in the equatorial plane before day 2, consistent with the scenario in which the ejecta is shaped after its interaction with an equatorial density enhancement (Munari et al. 2022).From the orientation of the radio lobes, which are perpendicular to the orbital plane, we can conclude that P.A. = 90 • corresponds to the scattering of light by material in the orbital plane.On the other hand, P.A. = 180 • corresponds to the scattering of light by material in the polar regions.
The polarized flux represents the spectrum of the scattered light.We observed a noticeable asymmetry of the polarized flux on day 2 -a strong blue-shifted peak with a radial velocity of v r ≈ −1200 km/s at a position angle P.A. ≈ 90 • and a small red-shifted peak with a radial velocity of v r ≈ +3000km/s at a position angle P.A. ≈ 180 • (Figure 9).The interaction between the ejecta and an equatorial density enhancement is responsible for the formation of a ring-like structure in the orbital plane (see Munari et al. 2022).As depicted in Figure 10, the west lobe is moving toward us and the portion of the torus that is moving away from us (red-shifted) is hidden by the west lobe.The material within the west lobe absorbs the scattered light from the red-shifted part of the torus.This half-torus structure is similar to the structure observed in radio wavelengths by O'Brien et al. (2006) after the 2006 outburst.It is possible to separate the equatorial from the polar scattering regions only in the case when the red-shifted component of the equatorial torus is hidden.Such asymmetry of the polarized flux and position angle is not observed during the next observations between days 5.9 and 6.7.The P.A. on days 5.9 and 6.7 shows that matter in the orbital plane is now responsible for the scattering.Further proof of the torus-like structure in the orbital plane is a triple-peaked profile of the permitted lines (Munari et al. 2022).The velocities of the double peak of the polarized flux across the H α emission line are highlighted in Figure 9 and represent the radial velocities of the slowly expanding torus.Our values of the radial velocities of the torus and the east arc are consistent with the observations, Similar torus-like structures have been observed in classical novae, but in this case, they are thought to arise from the propeller effect of the orbital motion on the ejecta (Chomiuk et al. 2014).Given the long orbital periods of symbiotic recurrent novae, it is not entirely clear if this effect could be at work or if the torus-like structure had been present in the system before the outburst.

Different types of objects -similar spectropolarimetric signatures
The similar spectropolarimetric signatures between RS Oph and AGNs relate to the following characteristics: -Polarized flux structure across the H α emission line.
-Position angle orientation relative to the radio axis.
Three-quarters of the AGNs with detected intrinsic polarization presented by Smith et al. (2002) exhibit a decrease in the degree of polarization across the core of the H α emission line.In two of them, Mrk 6 and Mrk 335, the red-shifted peak of the polarized flux is higher than the blue-shifted peak, similar to the polarized flux in RS Oph observed on day 5.9 and day 6.9 into the 2021 outburst (Fig. 9).In the case of RS Oph, the depolarization in the emission lines is due to the scattering material being located within the emission line region.A similar explanation was suggested by Cohen et al. (1999).Those authors observed weakly polarized narrow lines in a sample of 13 FR II radio galaxies.The narrow lines arise outside of the scattering region and the net effect would be a decrease of the degree of polarization in the core of the emission lines.A more complex explanation of the decrease in the degree of polarization for AGNs has been proposed by Smith et al. (2002), where the broad-line region takes the form of a rotating disc surrounded by a co-planar ring of scatterers -this geometry will produce a high degree of polarization in the wings of the H α emission line (see Smith et al. 2002, Sec 4.3).A similar explanation may also be valid for RS Oph if the accretion disc survives the outburst and lies in the same plane as the scattering region.
A position angle parallel to the radio axis (with P.A. = 90 • in the case of RS Oph) is similar to that observed in type 1 Seyfert galaxies, where the optical polarization P.A. is often parallel to the radio axis (Smith et al. 2004), whereas P.A. = 180 • is most similar to type 2 Seyfert galaxies, in which the optical polarization position angle is perpendicular to the projected radio source axis (e.g., Antonucci 1983;Brindle et al. 1990).Scattering in a plane perpendicular to the radio axis produces a position angle parallel to the radio axis.On the other hand, scattering in a plane aligned with the radio axis would produce a position angle perpendicular to the radio axis (Smith et al. 2002).We suggest that similar geometry, torus and polar scattering region can reproduce similar spectropolarimetric signatures.
Theoretical results based on hydrodynamical simulations for classical novae predict that much of the outflowing mass is focused into the equatorial plane and forms spiral arms (Shen & Quataert 2022).Moreover, 3D hydrodynamical simulations of RS Oph in the quiescent phase show that the circumstellar mass distribution is highly structured with a mass enhancement in the orbital plane (Walder et al. 2008;Booth et al. 2016).For RS Oph, the interaction between the nova ejecta with the accretion disc and focused mass in the orbital plane forms an equatorial ringlike structure and limits the ejecta within a bipolar direction perpendicular to the orbital plane (Booth et al. 2016).
For those novae that form an equatorial-like structure (or torus) during the outburst, the scattering geometry limits the orientation of the position angle, parallel or perpendicular to the torus, depending on the relative contribution of the equatorial and polar scattering region.

Conclusions
The optical linear spectropolarimetric observations in the range from 5000 Å to 8200 Å of the recurrent nova RS Oph after the 2021 outburst indicate that: -The variable degree of intrinsic polarization in RS Oph with time and wavelength suggests the presence of dust as early as day 2.This dust could have been created in the strongly shielded, dense shock regions (Evans et al. 2007).-The continuum intrinsic polarization can be fitted with a power law.The power law index β increased from β = −2.2 on day 2 to β = 0 on day 9.The flat polarization continuum suggests that dust was destroyed, and the scattering was then due to Thompson electron scattering.In the last three observations, after day 300 into the 2021 outburst, the power index β decreases to β=-1.6.This is an indicator of dust formation and the position angle P.A. = 90 • means that the dust is located in the orbital plane.-The polarized light two days after the outburst shows two position angles at 90 • and 180 • degrees.The position angle P.A. int ≈ 90 • is aligned with the radio/optical/X-rays elongated structures detected at much later times after the 2006 outburst.The P.A. int ≈ 180 • degrees indicates that during these early times, there is a region responsible for the light scattering that is perpendicular to the ejecta in the polar directions.These provide the first evidence of the development of asymmetries in the ejecta.-We suggest that similar geometry (torus and polar scattering region) presented in AGNs and RS Oph, can reproduce similar spectropolarimetric signatures.

Fig. 1 .
Fig. 1.Spectral evolution of RS Oph after the 2021 outburst.For comparison, we also include a spectrum taken in quiescence, on July 2019.The most prominent lines are highlighted.

Fig. 2 .
Fig. 2. Selected high-resolution emission lines in the spectra in quiescence (2019) and 40 and 285 days after the outburst.The line profiles are very complex, with at least three components at different radial velocities.

Fig. 3 .Fig. 4 .
Fig. 3.The interstellar polarization of the field stars around RS Oph (Heiles 2000; Cropper 1990).The degree of polarization is proportional to the length of the orange straight lines.The horizontal bar at the top left presents 4.0% polarization.The P.A. of the stars in the direction of RS Oph (with blue bar) is similar to that of RS Oph during quiescence (red line; Nikolov et al. 2019).The color of every star corresponds to its distance.The background image represents 100 µm dust emission maps (Schlegel et al. 1998).

Fig. 5 .
Fig. 5. Intrinsic polarization of RS Oph after the 2021 outburst.The degree of polarization is presented in blue.Strong depolarization effects are well visible in the first week after the start of the 2021 outburst (left column).The normalized intensity is presented in red.The orange line represents the power law fit of the degree of polarization data.The index of the law increases in the first week after the start of the 2021 outburst (left column).Later, from day 300 to day 352, the power law index decreases (right column).

Fig. 6 .
Fig.6.Intrinsic position angle of RS Oph after the 2021 outburst.Aside from the 180 • position angle around the Hα region on day 2, the position angle of the intrinsic polarization is on average 90 • , implying that the scattering region that produced the observed polarization was oriented perpendicular to the elongated structures observed at other wavelengths during the last two outbursts.The data were smoothed with a moving average over seven points to better illustrate the position angle structures.

Fig. 8 .
Fig.8.Long-term variability of the position angle (top) and degree of intrinsic polarization (bottom) of RS Oph after the start of the 2021 outburst.The P.A. was practically constant throughout our observing campaign, while the degree of intrinsic polarization reached a maximum at day 5.9 and has been decreasing since then.

Fig. 9 .
Fig.9.Evolution of the polarized flux across the H α emission line.The position angles (P.A.) are presented in red in the first row from day 2 (left) to day 6.7 (right) into the 2021 outburst.The P.A. on the 2nd day into the outburst (left column, top) shows a complex structure.The evolution of the polarized flux (in blue) and 0.01*flux (in orange) are presented in the second row.Asymmetries in the blue-and red-shifted peaks of the polarized flux are well visible.In the third row, the degree of polarization is presented.

Fig. 10 .
Fig. 10.Schematic view of the polar and equatorial scattering regions that are detected at different position angles and their contribution to the polarized flux.

Table 1 .
Journal of observations.
and [FeXI] at 6375 Å and 7892 Å respectively are visible on the later spectra, when RS Oph entered a supersoft source phase (Page 2021; Munari et al. 2021).
Nikolov et al.: RS Oph -asymmetry of the ejected material after the 2021 outburst.