Solar and Stellar Astrophysics (astro-ph.SR)

Abstract: Solar energetic particle (SEP) events are one of the most crucial aspects of space weather. Their prediction depends on various factors including the source solar eruptions such as flares and coronal mass ejections (CMEs). The Geostationary Solar Energetic Particle (GSEP) Events catalog was developed as an extensive data set towards this effort for solar cycles 22, 23 and 24. In the present work, we review and extend the GSEP data set by; (1) adding "weak" SEP events that have proton enhancements from 0.5 to 10 pfu in the E>10 MeV channel, and (2) improving the associated solar source eruptions information. We analyze and discuss spatio-temporal properties such as flare magnitudes, locations, rise times, and speed and width of CMEs. We check for the correlation of these parameters with peak proton fluxes and event fluences. Our study also focuses on understanding feature importance towards the optimal performance of machine learning (ML) models for SEP event forecasting. We implement random forest (RF), extreme gradient boosting (XGBoost), logistic regression (LR) and support vector machines (SVM) classifiers in a binary classification schema. Based on the evaluation of our best models, we find both the flare and CME parameters are requisites to predict the occurrence of an SEP event. This work is a foundation for our further efforts on SEP event forecasting using robust ML methods.

Abstract: The $\gamma$-process nucleosynthesis in core-collapse supernovae is generally accepted as a feasible process for the synthesis of neutron-deficient isotopes beyond iron. However, crucial discrepancies between theory and observations still exist: the average production of $\gamma$-process yields from massive stars are too low to reproduce the solar distribution in galactic chemical evolution calculations, and the yields of the Mo and Ru isotopes are by a further factor of 10 lower than the yields of the other $\gamma$-process nuclei. We investigate the $\gamma$-process in 5 sets of core-collapse supernova models published in literature with initial masses 15, 20, and 25 M$_{\odot}$ at solar metallicity. We compared the $\gamma$-process overproduction factors from the different models. To highlight the possible effect of nuclear physics input, we also considered 23 ratios of two isotopes close to each other in mass, relative to their solar values. Further, we investigated the contribution of C-O shell mergers in the supernova progenitors as an additional site of the $\gamma$-process. Our analysis shows that a large scatter among the different models exists for both the $\gamma$-process integrated yields and the isotopic ratios. We found only 10 ratios that agree with their solar values, all the others differ by at least a factor of 3 from the solar values in all the considered sets of models. The $\gamma$-process within C-O shell mergers mostly influence the isotopic ratios that involve intermediate and heavy proton-rich isotopes with $\rm A>100$.

Abstract: Space asteroseismology is revolutionizing our knowledge of the internal structure and dynamics of stars. A breakthrough is ongoing with the recent discoveries of signatures of strong magnetic fields in the core of red giant stars. The key signature for such a detection is the asymmetry these fields induce in the frequency splittings of observed dipolar mixed gravito-acoustic modes. We investigate the ability of the observed asymmetries of the frequency splittings of dipolar mixed modes to constrain the geometrical properties of deep magnetic fields. We use the powerful analytical Racah-Wigner algebra used in Quantum Mechanics to characterize the geometrical couplings of dipolar mixed oscillation modes with various possible realistic fossil magnetic fields' topologies and compute the induced perturbation of their frequencies. First, in the case of an oblique magnetic dipole, we provide the exact analytical expression of the asymmetry as a function of the angle between the rotation and magnetic axes. Its value provides a direct measure of this angle. Second, considering a combination of axisymmetric dipolar and quadrupolar fields, we show how the asymmetry is blind to unravel the relative strength and sign of each component. Finally, in the case of a given multipole, we show that a negative asymmetry is a signature of non-axisymmetric topologies. Therefore, asymmetries of dipolar mixed modes provide key but only partial information on the geometrical topology of deep fossil magnetic fields. Asteroseismic constraints should therefore be combined with spectropolarimetric observations and numerical simulations, which aim to predict the more probable stable large-scale geometries.

Abstract: We present a multi-band photometric analysis of CRTS J163819.6+03485, the first low mass ratio (LMR) contact binary system with a period under the contact binary (CB) period limit. The unprecedented combination of mass ratio and period makes this system unique for eclipsing binary (EB) research. Using new multi-band photometric observations, we explored the parameter space of this unique total EB system through a detailed scan in the mass ratio - inclination plane and using the PIKAIA genetic algorithm optimizer. The best set of relative physical parameters and corresponding uncertainties was adopted through Markov Chain Monte Carlo sampling of the parameter space. The resulting mass ratio of the system is $q = 0.16 \pm 0.01$. The absolute parameters were derived by adopting an empirical mass-luminosity relation. Period changes are also investigated by using new observations and archival photometric light curves from massive astronomical surveys, which revealed in a preliminary solution the presence of a possible low-mass tertiary companion. The origin and evolutionary status of the system are investigated through the detached-binary formation scenario.

Abstract: The ages of solar-like stars have been at the center of many studies such as exoplanet characterization or Galactic-archaeology. While ages are usually computed from stellar evolution models, relations linking ages to other stellar properties, such as rotation and magnetic activity, have been investigated. With the large catalog of 55,232 rotation periods, $P_{\rm rot}$, and photometric magnetic activity index, $S_{\rm ph}$ from Kepler data, we have the opportunity to look for such magneto-gyro-chronology relations. Stellar ages are obtained with two stellar evolution codes that include treatment of angular momentum evolution, hence using $P_{\rm rot}$ as input in addition to classical atmospheric parameters. We explore two different ways of predicting stellar ages on three subsamples with spectroscopic observations: solar analogs, late-F and G dwarfs, and K dwarfs. We first perform a Bayesian analysis to derive relations between $S_{\rm ph}$ and ages between 1 and 5 Gyr, and other stellar properties. For late-F and G dwarfs, and K dwarfs, the multivariate regression favors the model with $P_{\rm rot}$ and $S_{\rm ph}$ with median differences of 0.1%.and 0.2% respectively. We also apply Machine Learning techniques with a Random Forest algorithm to predict ages up to 14 Gyr with the same set of input parameters. For late-F, G and K dwarfs together, predicted ages are on average within 5.3% of the model ages and improve to 3.1% when including $P_{\rm rot}$. These are very promising results for a quick age estimation for solar-like stars with photometric observations, especially with current and future space missions.

Abstract: WD0810-353 is a white dwarf within the 20pc volume around the Sun. Using Gaia astrometric distance and proper motions, and a radial velocity derived from Gaia spectroscopy, it has been predicted that this star will pass within 1pc of the Solar System in about 30kyr. However, WD0810-353 has been also shown to host a magnetic field with strength of the order of 30MG. Its spectrum is therefore not like those of normal DA stars of similar effective temperature. We have obtained and analysed new polarised spectra of the star around Halpha. Our analysis suggests that the visible surface of the star shows two regions of different field strength (~30 and ~45MG, respectively), and opposite polarity. The spectra do not change over a 4 year time span, meaning that either the stellar rotation period is no shorter than several decades, or that the field is symmetric about the rotation axis. Taking into account magnetic shift and splitting, we obtain an estimate of the radial velocity of the star (+83+/- 140km/s); we reject both the value an the claimed precision deduced from the Gaia DR3 spectroscopy (-373.7+/- 8.2km/s), and we conclude that there will probably be no close encounter between the Solar System and WD0810-353. We also reject the suggestion that the star is a hypervelocity runaway star, a survivor of a Type Ia Supernova explosion. It is just a stellar remnant in the Solar neighborhood with a very strong and complex magnetic field.

Abstract: The Tayler-Spruit dynamo mechanism has been proposed two decades ago as a plausible mechanism to transport angular momentum in radiative stellar layers. Direct numerical simulations are still needed to understand its trigger conditions and the saturation mechanisms. The present study follows up on (Petitdemange et al. 2023), where we reported the first numerical simulations of a Tayler-Spruit dynamo cycle. Here we extend the explored parameter space to assess in particular the influence of stratification on the dynamo solutions. We also present numerical verification of theoretical assumptions made in (Spruit 2002), which are instrumental in deriving the classical prescription for angular momentum transport implemented in stellar evolution codes. A simplified radiative layer is modeled numerically by considering the dynamics of a stably-stratified, differentially rotating, magnetized fluid in a spherical shell. Our simulations display a diversity of magnetic field topologies and amplitudes depending on the flow parameters, including hemispherical solutions. The Tayler-Spruit dynamos reported here are found to satisfy magnetostrophic equilibrium and achieve efficient turbulent transport of angular momentum, following Spruit's heuristic prediction.

Abstract: Very massive stars (VMS) up to 200-300 $M_\odot$ have been found in the Local Universe. If they would lose little mass they produce intermediate-mass black holes or pair-instability supernovae (PISNe). Until now, VMS modellers have extrapolated mass-loss vs. metallicity ($Z$) exponents from optically-thin winds, resulting in a range of PISN thresholds that might be unrealistically high in $Z$, as VMS develop optically-thick winds. We utilize the transition mass-loss rate of Vink and Gr\"afener (2012) that accurately predicts mass-loss rates of Of/WNh ("slash") stars that characterize the morphological transition from absorption-dominated O-type spectra to emission-dominated WNh spectra. We develop a wind efficiency framework, where optically thin winds transition to enhanced winds, enabling us to study VMS evolution at high redshift where individual stars cannot be resolved. We present a MESA grid covering $Z_\odot/2$ to $Z_\odot/100$. VMS above the transition evolve towards lower luminosity, skipping the cool supergiant phase but directly forming pure He stars at the end of hydrogen burning. Below the transition, VMS evolve as cooler luminous blue variables (LBVs) or yellow hypergiants (YHGs), naturally approaching the Eddington limit. Strong winds in this YHG/LBV regime -- combined with a degeneracy in luminosity -- result in a mass-loss runaway where a decrease in mass increases wind mass loss. Our models indicate an order-of-magnitude lower threshold than usually assumed, at $Z_\odot/20$ due to our mass-loss runaway. While future work on LBV mass loss could affect the PISN threshold, our framework will be critical for establishing definitive answers on the PISN threshold and galactic chemical evolution modelling.

Abstract: Massive stars are progenitors of supernovae, neutron stars and black holes. During the hydrogen-core burning phase their convective cores are the prime drivers of their evolution, but inferences of core masses are subject to unconstrained boundary mixing processes. Moreover, uncalibrated transport mechanisms can lead to strong envelope mixing and differential radial rotation. Ascertaining the efficiency of the transport mechanisms is challenging because of a lack of observational constraints. Here we deduce the convective core mass and robustly demonstrate non-rigid radial rotation in a supernova progenitor, the $12.0^{+1.5}_{-1.5}$ solar-mass hydrogen-burning star HD 192575, using asteroseismology, TESS photometry, high-resolution spectroscopy, and Gaia astrometry. We infer a convective core mass ($M_{\rm cc} = 2.9^{+0.5}_{-0.8}$ solar masses), and find the core to be rotating between 1.4 and 6.3 times faster than the stellar envelope depending on the location of the rotational shear layer. Our results deliver a robust inferred core mass of a massive star using asteroseismology from space-based photometry. HD 192575 is a unique anchor point for studying interior rotation and mixing processes, and thus also angular momentum transport mechanisms inside massive stars.

Abstract: White dwarfs stars are known to be polluted by their active planetary systems, but little attention has been paid to the accretion of wind from low-mass companions. The capture of stellar or substellar wind by white dwarfs is one of few methods available to astronomers which can assess mass-loss rates from unevolved stars and brown dwarfs, and the only known method to extract their chemical compositions. In this work, four white dwarfs with closely-orbiting, L-type brown dwarf companions are studied to place limits on the accretion of a substellar wind, with one case of a detection, and at an extremely non-solar abundance $m_{\rm Na}/m_{\rm Ca}>900$. The mass-loss rates and upper limits are tied to accretion in the white dwarfs, based on limiting cases for how the wind is captured, and compared with known cases of wind pollution from close M dwarf companions, which manifest in solar proportions between all elements detected. For wind captured in a Bondi-Hoyle flow, mass-loss limits $\dot M\lesssim 5\times10^{-17}$ M$_\odot$ yr$^{-1}$ are established for three L dwarfs, while for M dwarfs polluting their hosts, winds in the range $10^{-13} - 10^{-16}$ M$_\odot$ yr$^{-1}$ are found. The latter compares well with the $\dot M\sim 10^{-13} - 10^{-15}$ M$_\odot$ yr$^{-1}$ estimates obtained for nearby, isolated M dwarfs using Ly$\alpha$ to probe their astropsheres. These results demonstrate that white dwarfs are highly-sensitive stellar and substellar wind detectors, where further work on the actual captured wind flow is needed.

Abstract: We present the first JWST spectral energy distribution of a Y dwarf. This spectral energy distribution of the Y0 dwarf WISE J035934.06$-$540154.6 consists of low-resolution ($\lambda$/$\Delta\lambda$ $\sim$ 100) spectroscopy from 1$-$12 $\mu$m and three photometric points at 15, 18, and 21 $\mu$m. The spectrum exhibits numerous fundamental, overtone, and combination rotational-vibrational bands of H$_2$O, CH$_4$, CO, CO$_2$, and NH$_3$, including the previously unidentified $\nu_3$ band of NH$_3$ at 3 $\mu$m. Using a Rayleigh-Jeans tail to account for the flux emerging at wavelengths greater than 21 $\mu$m, we measure a bolometric luminosity of $1.523\pm0.090\times10^{20}$ W. We determine a semi-empirical effective temperature estimate of $467^{+16}_{-18}$ K using the bolometric luminosity and evolutionary models to estimate a radius. Finally, we compare the spectrum and photometry to a grid of atmospheric models and find reasonably good agreement with a model having $T_{\mathrm{eff}}$=450 K, log $g$=3.25 [cm s$^{-2}$], [M/H]=$-0.3$. However, the low surface gravity implies an extremely low mass of 1 $M_{\rm{Jup}}$ and a very young age of 20 Myr, the latter of which is inconsistent with simulations of volume-limited samples of cool brown dwarfs.

Abstract: Recent observations demonstrate that misalignments and other out-of-plane structures are common in protoplanetary discs. Many of these have been linked to a central host binary with an orbit that is inclined with respect to the disc. We present simulations of misaligned circumbinary discs with a range of parameters to gain a better understanding of the link between those parameters and the disc morphology in the wave-like regime of warp propagation that is appropriate to protoplanetary discs. The simulations confirm that disc tearing is possible in protoplanetary discs as long as the mass ratio, $\mu$, and disc-binary inclination angle, $i$, are not too small. For the simulations presented here this corresponds to $\mu > 0.1$ and $i \gtrsim 40^\circ$. For highly eccentric binaries, tearing can occur for discs with smaller misalignment. Existing theoretical predictions provide an estimate of the radial extent of the disc in which we can expect breaking to occur. However, there does not seem to be a simple relationship between the disc properties and the radius within the circumbinary disc at which the breaks appear, and furthermore the radius at which the disc breaks can change as a function of time in each case. We discuss the implications of our results for interpreting observations and suggest some considerations for modelling misaligned discs in the future.

Abstract: We present optical and infrared imaging and spectroscopy of the R Coronae Borealis-type (R Cor Bor) star IRAS 00450+7401. Optical spectra further confirm its classification as a cool R Cor Bor system, having a hydrogen-deficient carbon star spectral sub-class of HdC5 or later. Mid-infrared spectroscopy reveals the typical ~8 um ``hump'' seen in other R Cor Bor stars and no other features. A modern-epoch spectral energy distribution shows bright emission from hot dust having Tdust>600 K. Historical infrared data reveal generally cooler dust color temperatures combined with long-term fading trends, but provide no discernible correlation between flux level and temperature. Investigating the most mid-infrared variable R Cor Bor stars found in IRAS, AKARI, and WISE data reveals similar fading trends, bursts that can show a factor of up to 10 change in flux density between epochs, and blackbody-fit dust color temperatures that span 400-1300 K. While some R Cor Bor stars such as IRAS 00450+7401 appear to undergo fade/burst cycles in the mid-infrared, significant gaps in temporal coverage prevent conclusively identifying any preferred timescale for their mid-infrared variability and circumstellar dust temperature changes.