Solar and Stellar Astrophysics (astro-ph.SR)

Abstract: The solar inner corona is a region that plays a critical role in energizing the solar wind and propelling it to supersonic and supra-Alfvenic velocities. Despite its importance, this region remains poorly understood because of being least explored due to observational limitations. The coronal radio sounding technique in this context becomes useful as it helps in providing information in parts of this least explored region. To shed light on the dynamics of the solar wind in the inner corona, we conducted a study using data obtained from coronal radio-sounding experiments carried out by the Akatsuki spacecraft during the 2021 Venus-solar conjunction event. By analyzing X-band radio signals recorded at two ground stations (IDSN in Bangalore and UDSC in Japan), we investigated plasma turbulence characteristics and estimated flow speed measurements based on isotropic quasi-static turbulence models. Our analysis revealed that the speed of the solar wind in the inner corona (at heliocentric distances from 5 to 13 solar radii), ranging from 220-550 km/sec, was higher than the expected average flow speeds in this region. By integrating our radio-sounding results with EUV images of the solar disk, we gained a unique perspective on the properties and energization of high-velocity plasma streams originating from coronal holes. We tracked the evolution of fast solar wind streams emanating from an extended coronal hole as they propagated to increasing heliocentric distances. Our study provides unique insights into the least-explored inner coronal region by corroborating radio sounding results with EUV observations of the corona.

Abstract: Optical spectroscopy offers the most direct view of the stellar properties and the accretion indicators. Standard accretion tracers, such as $H\beta$, $H\alpha$, and, Ca II triplet lines, and most photospheric features, fall in the optical wavelengths. However, these tracers are not readily observable from deeply embedded protostars because of the large line of sight extinction (Av $\sim$ 50-100 mag) toward them. In some cases, however, it is possible to observe protostars at optical wavelengths if the outflow cavity is aligned along the line-of-sight that allows observations of the photosphere, or the envelope is very tenuous and thin such that the extinction is low. In such cases, we can not only detect these protostars at optical wavelengths but also follow up spectroscopically. We have used the HOPS catalog (Furlan et al. 2016) of protostars in Orion to search for optical counterparts for protostars in the Gaia DR3 survey. Out of the 330 protostars in the HOPS sample, an optical counterpart within 2" is detected for 62 of the protostars. For 17 out of 62 optically detected protostars, we obtained optical spectra { (between 5500 to 8900 $\AA$) using the Aries-Devasthal Faint Object Spectrograph \& Camera (ADFOSC) on the 3.6-m Devasthal Optical Telescope (DOT) and Hanle Faint Object Spectrograph Camera (HFOSC) on 2-m Himalayan Chandra Telescope (HCT)}. We detect strong photospheric features, such as the TiO bands in the spectra {(of 4 protostars)}, hinting that photospheres can form early on in the star formation process. We further determined the spectral types of protostars, which show photospheres similar to a late M-type. Mass accretion rates derived for the protostars are similar to those found for T-Tauri stars, in the range of 10$^{-7}$ to 10$^{-8}$ $M_\odot$/yr.

Abstract: Asteroseismology of solar-like oscillations in giant stars allow the derivation of their masses and radii. For members of open clusters this allows an age estimate of the cluster which should be identical to the age estimate from the colour-magnitude diagram, but independent of the uncertainties that are present for that type of analysis. Thus, a more precise and accurate age estimate can be obtained. We aim to measure asteroseismic properties of oscillating giant members of the open cluster NGC 6866 and utilise these for a cluster age estimate. Model comparisons allow constraints on the stellar physics, and here we investigate the efficiency of convective-core overshoot and effects of rotation during the main-sequence, which has a significant influence on the age for these relatively massive giants. We identify six giant members of NGC 6866 and derive asteroseismic measurements for five of them. This constrains the convective-core overshoot and enables a more precise and accurate age estimate than previously possible. Asteroseismology establishes the helium-core burning evolutionary phase for the giants, which have a mean mass of 2.8 $M_{\odot}$. Their radii are significantly smaller than predicted by current 1D stellar models unless the amount of convective-core overshoot on the main sequence is reduced to $\alpha_{ov} \leq 0.1 H_p$ in the step-overshoot description. Our measurements also suggest that rotation has affected the evolution of the stars in NGC 6866 in a way that is consistent with 3D simulations but not with current 1D stellar models. The cluster age is estimated to be 0.43 $\pm$ 0.05 Gyr, significantly younger and more precise than most previous estimates. We derive a precise cluster age while constraining convective-core overshooting and effects of rotation in the models. We uncover potential biases for automated age estimates of helium-core burning stars.

Abstract: We aim to combine asteroseismology, spectroscopy, and evolutionary models to establish a comprehensive picture of the evolution of Galactic blue supergiant stars (BSG). To start such an investigation, we selected three BSG candidates for our analysis: HD 42087 (PU Gem), HD 52089 ($\epsilon$ CMa) and HD 58350 ($\eta$ CMa). These stars show pulsations and were suspected to be in an evolutionary stage either preceding or succeding the red supergiant (RSG) stage. For our analysis, we utilized the 2-min cadence TESS data to study the photometric variability and obtained new spectroscopic observations at the CASLEO observatory. We calculated CMFGEN non-LTE radiative transfer models and derived stellar and wind parameters using the iterative spectral analysis pipeline XTGRID. The spectral modeling was limited to changing only the effective temperature, surface gravity, CNO abundances, and mass-loss rates. Finally, we compared the derived metal abundances with predictions from Geneva stellar evolution models. The frequency spectra of all three stars show either stochastic oscillations, nonradial strange modes, or a rotational splitting. We conclude that the rather short sectoral observing windows of TESS prevent establishing a reliable mode identification of low frequencies connected to mass-loss variabilities. The spectral analysis confirmed gradual changes in the mass-loss rates and the derived CNO abundances comply with the values reported in the literature. We were able to achieve a quantitative match with stellar evolution models for the stellar masses and luminosities. However, the spectroscopic surface abundances turned out to be inconsistent with theoretical predictions. The stars show N enrichment, typical for CNO cycle processed material, but the abundance ratios do not reflect the associated levels of C and O depletion.

Abstract: Inertial wave modes in the Sun are of interest owing to their potential to reveal new insight into the solar interior. These predominantly retrograde-propagating modes in the solar subsurface appear to deviate from the thin-shell Rossby-Haurwitz model at high azimuthal orders. We present new measurements of sectoral equatorial inertial modes at $m>15$ where the modes appear to become progressively less retrograde compared to the canonical Rossby-Haurwitz dispersion relation in a co-rotating frame. We use a spectral eigenvalue solver to compute the spectrum of solar inertial modes in the presence of differential rotation. Focussing specifically on equatorial Rossby modes, we find that the numerically obtained mode frequencies lie along distinct ridges, one of which lies strikingly close to the observed mode frequencies in the Sun. We also find that the $n=0$ ridge is deflected strongly in the retrograde direction. This suggests that the solar measurements may not correspond to the fundamental $n=0$ Rossby-Haurwitz solutions as was initially suspected, but to a those for a higher $n$. The numerically obtained eigenfunctions also appear to sit deep within the convection zone -- unlike those for the $n=0$ modes -- which differs substantially from solar measurements and complicates inference.

Abstract: Explosive phenomena are known to trigger a wealth of shocks in warm plasma environments, including the solar chromosphere and molecular clouds where the medium consists of both ionised and neutral species. Partial ionisation is critical in determining the behaviour of shocks, since the ions and neutrals locally decouple, allowing for substructure to exist within the shock. Accurately modelling partially ionised shocks requires careful treatment of the ionised and neutral species, and their interactions. Here we study a partially-ionised switch-off slow-mode shock using a multi-level hydrogen model with both collisional and radiative ionisation and recombination rates that are implemented into the two-fluid (P\underline{I}P) code, and study physical parameters that are typical of the solar chromosphere. The multi-level hydrogen model differs significantly from MHD solutions due to the macroscopic thermal energy loss during collisional ionisation. In particular, the plasma temperature both post-shock and within the finite-width is significantly cooler that the post-shock MHD temperature. Furthermore, in the mid to lower chromosphere, shocks feature far greater compression then their single-fluid MHD analogues. The decreased temperature and increased compression reveal the importance of non-equilibrium ionised in the thermal evolution of shocks in partially ionised media. Since partially ionised shocks are not accurately described by the Rankine-Hugoniot shock jump conditions, it may be incorrect to use these to infer properties of lower atmospheric shocks.

Abstract: Radial substructures have now been observed in a wide range of protoplanetary discs (PPDs), from young to old systems, however their formation is still an area of vigorous debate. Recent magnetohydrodynamic (MHD) simulations have shown that rings and gaps can form naturally in PPDs when non-ideal MHD effects are included. However these simulations employ ad-hoc approximations to the magnitudes of the magnetic diffusivities in order to facilitate ring growth. We replace the parametrisation of these terms with a simple chemical network and grain distribution model to calculate the non-ideal effects in a more self-consistent way. We use a range of grain distributions to simulate grain formation for different disc conditions. Including ambipolar diffusion, we find that large grain populations (> 1{\mu}m), and those including a population of very small polyaromatic hydrocarbons (PAHs) facilitate the growth of periodic, stable rings, while intermediate sized grains suppress ring formation. Including Ohmic diffusion removes the positive influence of PAHs, with only large grain populations still producing periodic ring and gap structures. These results relate closely to the degree of coupling between the magnetic field and the neutral disc material, quantified by the non-dimensional Elsasser number {\Lambda} (the ratio of magnetic forces to Coriolis force). For both the ambipolar-only and ambipolar-ohmic cases, if the total Elsasser number is initially of order unity along the disc mid-plane, ring and gap structures may develop.