Solar and Stellar Astrophysics (astro-ph.SR)

Abstract: We present the results of optical spectroscopy of stellar companions to three binary millisecond pulsars, PSRs J0621$+$2514, J2317$+$1439 and J2302$+$4442, obtained with the Gran Telescopio Canarias. The spectrum of the J0621$+$2514 companion shows a blue continuum and prominent Balmer absorption lines. The latter are also resolved in the spectrum of the J2317$+$1439 companion, showing that both are DA-type white dwarfs. No spectral features are detected for the J2302$+$4442 companion, however, its broadband magnitudes and the spectral shape of the continuum emission imply that this is also a DA-type white dwarf. Based on the spectral analyses, we conclude that the companions of J0621$+$2514 and J2317$+$1439 are relatively hot, with effective temperatures $T_{\rm eff}$$=$8600$\pm$200 and 9600$\pm$2000~K, respectively, while the J2302$+$4442 companion is significantly cooler, $T_{\rm eff}$$<$6000~K. We also estimated the distance to J0621$+$2514 of 1.1$\pm$0.3 kpc and argue that its companion and the companion of J2317$+$1439 are He-core white dwarfs providing constraints on their cooling ages of $\lesssim$2 Gyr.

Abstract: We investigate the origin of mesoscale structures in the solar wind called microstreams defined as enhancements in solar wind speed and temperature that last several hours. They were first clearly detected in Helios and Ulysses solar wind data and are now omnipresent in the "young" solar wind measured by Parker Solar Probe and Solar Orbiter. These recent data reveal that microstreams transport a profusion of Alfv\'enic perturbations in the form of velocity spikes and magnetic switchbacks. In this study we use a very high-resolution 2.5 MHD model of the corona and the solar wind to simulate the emergence of magnetic bipoles interacting with the pre-existing ambient corona and the creation of jets that become microstreams propagating in the solar wind. Our high-resolution simulations reach sufficiently high Lundquist numbers to capture the tearing mode instability that develops in the reconnection region and produces plasmo\"ids released with the jet into the solar wind. Our domain runs from the lower corona to 20 Rs, this allows us to track the formation process of plasmo\"ids and their evolution into Alfv\'enic velocity spikes. We obtain perturbed solar wind flows lasting several hours with velocity spikes occurring at characteristic periodicities of about 19 minutes. We retrieve several properties of microstreams measured in the pristine solar wind by Parker Solar Probe, namely an increase in wind velocity of about 100 km/s during the streams passage together with superposed velocity spikes of also about 100 km/s released into the solar wind.

Abstract: The Solar Orbiter mission completed its first remote-sensing observation windows in the spring of 2022. On 2/4/2022, an M-class flare followed by a filament eruption was seen both by the instruments on board the mission and from several observatories in Earth's orbit. The complexity of the observed features is compared with the predictions given by the standard flare model in 3D. We use the observations from a multi-view dataset, which includes EUV imaging to spectroscopy and magnetic field measurements. These data come from IRIS, SDO, Hinode, as well as several instruments on Solar Orbiter. Information given by SDO/HMI and Solar Orbiter PHI/HRT shows that a parasitic polarity emerging underneath the filament is responsible for bringing the flux rope to an unstable state. As the flux rope erupts, Hinode/EIS captures blue-shifted emission in the transition region and coronal lines in the northern leg of the flux rope prior to the flare peak. Solar Orbiter SPICE captures the whole region, complementing the Doppler diagnostics of the filament eruption. Analyses of the formation and evolution of a complex set of flare ribbons and loops show that the parasitic emerging bipole plays an important role in the evolution of the flaring region. While the analysed data are overall consistent with the standard flare model, the present particular magnetic configuration shows that surrounding magnetic activity such as nearby emergence needs to be taken into account to fully understand the processes at work. This filament eruption is the first to be covered from different angles by spectroscopic instruments, and provides an unprecedented diagnostic of the multi-thermal structures present before and during the flare. This dataset of an eruptive event showcases the capabilities of coordinated observations with the Solar Orbiter mission.

Abstract: Atmospheric gravity waves (AGWs) are low-frequency, buoyancy-driven waves that are generated by turbulent convection and propagate obliquely throughout the solar atmosphere. Their proposed energy contribution to the lower solar atmosphere and sensitivity to atmospheric parameters (e.g. magnetic fields and radiative damping) highlight their diagnostic potential. We investigate AGWs near a quiet Sun disk center region using multi-wavelength data from the Interferometric BIdimensional Spectrometer (IBIS) and the Solar Dynamics Observatory (SDO). These observations showcase the complex wave behavior present in the entire acoustic-gravity wave spectrum. Using Fourier spectral analysis and local helioseismology techniques on simultaneously observed line core Doppler velocity and intensity fluctuations, we study both the vertical and horizontal properties of AGWs.Propagating AGWs with perpendicular group and phase velocities are detected at the expected temporal and spatial scales throughout the lower solar atmosphere. We also find previously unobserved, varied phase difference distributions among our velocity and intensity diagnostic combinations. Time-distance analysis indicates that AGWs travel with an average group speed of 4.5 kms$^{-1}$, which is only partially described by a simple simulation suggesting that high-frequency AGWs dominate the signal. Analysis of the median magnetic field (4.2 G) suggests that propagating AGWs are not significantly affected by quiet Sun photospheric magnetic fields. Our results illustrate the importance of multi-height observations and the necessity of future work to properly characterize this observed behavior.

Abstract: We compute the evidence in favour of two models, one based on field-aligned thermal conduction alone and another that includes thermal misbalance as well, in explaining the damping of slow magneto-acoustic waves in hot coronal loops. Our analysis is based on the computation of the marginal likelihood and the Bayes factor for the two damping models. We quantify their merit in explaining the apparent relationship between slow mode periods and damping times, measured with SOHO/SUMER in a set of hot coronal loops. The results indicate evidence in favour of the model with thermal misbalance in the majority of the sample, with a small population of loops for which thermal conduction alone is more plausible. The apparent possibility of two different regimes of slow-wave damping, if due to differences between the loops of host active regions and/or the photospheric dynamics, may help with revealing the coronal heating mechanism.

Abstract: Poynting flux is the flux of magnetic energy, which is responsible for chromospheric and coronal heating in the solar atmosphere. It is defined as a cross product of electric and magnetic fields, and in ideal MHD conditions it can be expressed in terms of magnetic field and plasma velocity. Poynting flux has been computed for active regions and plages, but estimating it in the quiet Sun (QS) remains challenging due to resolution effects and polarimetric noise. However, with upcoming DKIST capabilities, these estimates will become more feasible than ever before. Here, we study QS Poynting flux in Sunrise/IMaX observations and MURaM simulations. We explore two methods for inferring transverse velocities from observations - FLCT and a neural network based method DeepVel - and show DeepVel to be the more suitable method in the context of small-scale QS flows. We investigate the effect of azimuthal ambiguity on Poynting flux estimates, and we describe a new method for azimuth disambiguation. Finally, we use two methods for obtaining the electric field. The first method relies on idealized Ohm's law, whereas the second is a state-of-the-art inductive electric field inversion method PDFI SS. We compare the resulting Poynting flux values with theoretical estimates for chromospheric and coronal energy losses and find that some of Poynting flux estimates are sufficient to match the losses. Using MURaM simulations, we show that photospheric Poynting fluxes vary significantly with optical depth, and that there is an observational bias that results in underestimated Poynting fluxes due to unaccounted shear term contribution.

Abstract: We examine the East-West asymmetry of the peak intensity in energetic storm particle (ESP) events using the improved Particle Acceleration and Transport in the Heliosphere (iPATH) model. We find that injection efficiency peaks east of the nose of coronal mass ejection shock where the shock exhibits a quasi-parallel geometry. We show that the peak intensity at the eastern flank is generally larger than that at the western flank and it positively correlates with the injection efficiency. We also examine this asymmetry for heavy ions, which depends sensitively on the ion energy. Comparison between the modelling results with the measurements of ESP events at 1 au shows a reasonable agreement. We suggest that the injection efficiency can be a primary factor leading to the East-West asymmetry of the peak intensity in ESP events. Additionally, the charge-to-mass (Q/A) dependence of the maximum particle energy affects this asymmetry for heavy ions.

Abstract: This paper reports the ULTRACAM discovery of dipolar surface spots in two cool magnetic white dwarfs with Balmer emission lines, while a third system exhibits a single spot, similar to the prototype GD 356. The light curves are modeled with simple, circular, isothermal dark spots, yielding relatively large regions with minimum angular radii of 20 deg. For those stars with two light curve minima, the dual spots are likely observed at high inclination (or colatitude), however, identical and antipodal spots cannot simultaneously reproduce both the distinct minima depths and the phases of the light curve maxima. The amplitudes of the multi-band photometric variability reported here are all several times larger than that observed in the prototype GD 356; nevertheless, all DAHe stars with available data appear to have light curve amplitudes that increase toward the blue in correlated ratios. This behavior is consistent with cool spots that produce higher contrasts at shorter wavelengths, with remarkably similar spectral properties given the diversity of magnetic field strengths and rotation rates. These findings support the interpretation that some magnetic white dwarfs generate intrinsic chromospheres as they cool, and that no external source is responsible for the observed temperature inversion. Spectroscopic time-series data for DAHe stars is paramount for further characterization, where it is important to obtain well-sampled data, and consider wavelength shifts, equivalent widths, and spectropolarimetry.

Abstract: Binary stars evolve into chemically-peculiar objects and are a major driver of the Galactic enrichment of heavy elements. During their evolution they undergo interactions, including tides, that circularize orbits and synchronize stellar spins, impacting both individual systems and stellar populations. Using Zahn's tidal theory and MESA main-sequence model grids, we derive the governing parameters $\lambda_{lm}$ and $E_2$, and implement them in the new MINT library of the stellar population code BINARY_C. Our MINT equilibrium tides are 2 to 5 times more efficient than the ubiquitous BSE prescriptions while the radiative-tide efficiency drops sharply with increasing age. We also implement precise initial distributions based on bias-corrected observations. We assess the impact of tides and initial orbital-parameter distributions on circularization and synchronization in eight open clusters, comparing synthetic populations and observations through a bootstrapping method. We find that changing the tidal prescription yields no statistically-significant improvement as both calculations typically lie within 0.5$\sigma$. The initial distribution, especially the primordial concentration of systems at $\log_{10}(P/{\rm d}) \approx 0.8, e\approx 0.05$ dominates the statistics even when artificially increasing tidal strength. This confirms the inefficiency of tides on the main sequence and shows that constraining tidal-efficiency parameters using the $e-\log_{10}(P/{\rm d})$ distribution alone is difficult or impossible. Orbital synchronization carries a more striking age-dependent signature of tidal interactions. In M35 we find twice as many synchronized rotators in our MINT calculation as with BSE. This measure of tidal efficiency is verifiable with combined measurements of orbital parameters and stellar spins.

Abstract: In this paper we revisit metal-enriched pair instability supernovae (PISNe) models which undergo chemically homogeneous evolution (CHE). By calculating multiple models, we intend to clarify mass ranges for the PISNe, $^{56}$Ni masses from the PISNe, and mass loss histories of CHE-PISNe models for metallicities consistent with the Small Magellanic Cloud (SMC) and with the Large Magellanic Cloud (LMC). We show that for an initial velocity of $v_{\rm i}/ v_{\rm k}$ = 0.1, these models undergo CHE and He-rich (Type Ib) PISNe occur in a lower mass range ($M_{\rm i} \sim 110-170 M_\odot$) than for more slowly rotating models. Interestingly, bright PISNe which have $^{56}$Ni masses larger than 10 $M_\odot$ occur in a relatively small mass range, $M_{\rm i} \sim 140-170 M_\odot$. Another notable characteristic of CHE-PISNe is the large late time mass loss rates; consequently, CSM interaction may be observable in their light curves. We also show some examples of O-rich (Type Ic) CHE-PISNe produced by $v_{\rm i}/ v_{\rm k}$ = 0.2 models. We expect these models to exhibit interaction with O-rich CSM, behavior which is consistent with the observed properties of the recently discovered PISN candidate, SN2018ibb. Finally, we present a collapsing $v_{\rm i}/ v_{\rm k}$ = 0.2 model which has sufficient angular momentum to be regarded as a candidate for a Super-Kilonova.