Search

The $\sim$5 Myr PDS 70 is the only known system with protoplanets residing in the cavity of the circumstellar disk from which they formed, ideal for studying exoplanet formation and evolution within its natal environment. Here we report the first spin constraint and C/O measurement of PDS 70b from Keck/KPIC high-resolution spectroscopy. We detected CO (3.8 $\sigma$) and H$_2$O (3.5 $\sigma$) molecules in the PDS 70b atmosphere via cross-correlation, with a combined CO and H$_2$O template detection significance of 4.2 $\sigma$. Our forward model fits, using BT-Settl model grids, provide an upper limit for the spin-rate of PDS 70b ($<$29 km s$^{-1}$). The atmospheric retrievals constrain the PDS 70b C/O ratio to ${0.28}^{+0.20}_{-0.12}$ ($<$0.63 under 95$\%$ confidence level) and a metallicity [C/H] of ${-0.2}^{+0.8}_{-0.5}$ dex, consistent with that of its host star. The following scenarios can explain our measured C/O of PDS 70b in contrast with that of the gas-rich outer disk (for which C/O $\gtrsim$ 1). First, the bulk composition of PDS 70b might be dominated by dust+ice aggregates rather than disk gas. Another possible explanation is that the disk became carbon-enriched $\textit{after}$ PDS 70b was formed, as predicted in models of disk chemical evolution and as observed in both very low mass star and older disk systems with $\textit{JWST}$/MIRI. Because PDS 70b continues to accrete and its chemical evolution is not yet complete, more sophisticated modeling of the planet and the disk, and higher quality observations of PDS 70b (and possibly PDS 70c), are necessary to validate these scenarios., Comment: Accepted to ApJ Letters; 15 pages, 3 figures

The ground-based measurement of gravitational waves (GW) from merging binary black holes (BBH) uniquely allows determination of spins of stellar-remnant black holes (BH), thereby offering insights into their formation mechanisms. The observed population of BBH mergers exhibits two intriguing peculiarities related to BH spins, namely, a positively biased distribution of effective spin parameter, $\chi_{\rm eff}$, and an apparent anti-correlation between merger mass ratio, $q$, and $\chi_{\rm eff}$. Here we investigate the potential mechanisms for such observed properties, in BBH mergers via isolated binary evolution. We synthesize BBH mergers with the fast binary evolution code BSE. The role of various physical assumptions is explored, including tidal spin-up, compact remnant mass, and mass transfer physics. We compare the properties of WR-BH binaries that formed through stable mass transfer (SMT) and common envelope evolution (CE). We find that both the asymmetry in $\chi_{\rm eff}$ distribution and the $\chi_{\rm eff}-q$ anti-correlation are natural outcomes of isolated-binary BBH formation with the standard assumptions on mass transfer physics. The anti-correlation is particularly pronounced for SMT-channel BBH mergers that experience a mass-ratio-reversal, i.e., those where the second-born BH is the more massive member. The anti-correlation arises from the dependence of orbital shrinking during mass transfer and the Roche lobe size on the system's mass ratio. This characteristic $\chi_{\rm eff}-q$ trend diminishes with increasing metallicity and when the isolated-binary BBH merger population is mixed with a significant contribution of dynamically formed BBH mergers. Our results demonstrate that isolated massive binary evolution via the SMT sub-channel can reproduce the observable BBH merger population, with characteristic signatures in mass, mass ratio, and spin distributions.[abr.], Comment: 25 pages, 18 figures, 3 appendices. Submitted to Astronomy and Astrophysics. Comments are welcome

In this paper we provide a detailed investigation of the energisation processes in two-dimensional, two and a half-dimensional and three-dimensional collapsing magnetic trap models. Using kinematic magnetohydrodynamic models of collapsing magnetic traps, we examine the importance of Fermi acceleration in comparison with betatron acceleration in these models. We extend previous work by investigating particle orbits in two-dimensional models without and with a guide field component and from full three-dimensional models. We compare the outcomes for the different models and how they depend on the chosen initial conditions. While in the literature betatron acceleration has been emphasised as the major mechanism for particle energisation in collapsing magnetic traps, we find that Fermi acceleration can play a significant role as well for particle orbits with suitable initial conditions., Comment: Accepted for publication in MNRAS, 15 pages, 15 figures

Internal gravity waves (hereafter IGWs) are one of the mechanisms that can play a key role to redistribute efficiently angular momentum in stars along their evolution. The study of IGWs is thus of major importance since space-based asteroseismology reveals a transport of angular momentum in stars, which is stronger by two orders of magnitude than the one predicted by stellar models ignoring their action or those of magnetic fields. IGWs trigger angular momentum transport when they are damped by heat or viscous diffusion, when they meet a critical layer or when they break. Theoretical prescriptions have been derived for the transport of angular momentum induced by IGWs because of their radiative and viscous dampings and of the critical layers they encounter along their propagation. However, none has been proposed for the transport triggered by their nonlinear breaking. In this work, we aim to derive such a physical and robust prescription, which can be implemented in stellar structure and evolution codes. We adapt an analytical saturation model, which has been developed for IGWs nonlinear convective breaking in the Earth atmosphere and has been successfully compared to in-situ measurements in the stratosphere, to the case of deep spherical stellar interiors. In a first step, we neglect the modification of IGWs by the Coriolis acceleration and the Lorentz force, which are discussed and taken into account in a second step. We derive a complete semi-analytical prescription for the transport of angular momentum by IGWs, which takes into account both their radiative damping and their potential nonlinear breaking because of their convective and vertical shear instabilities. This allows us to bring the physical prescription for the interactions between IGWs and the differential rotation to the same level of realism that the one used in global circulation models for the atmosphere., Comment: 9 pages, 1 figure, accepted for publication in Astronomy & Astrophysics

BabyIAXO is the intermediate stage of the International Axion Observatory (IAXO) to be hosted at DESY. Its primary goal is the detection of solar axions following the axion helioscope technique. Axions are converted into photons in a large magnet that is pointing to the sun. The resulting X-rays are focused by appropriate X-ray optics and detected by sensitive low-background detectors placed at the focal spot. The aim of this article is to provide an accurate quantitative description of the different components (such as the magnet, optics, and X-ray detectors) involved in the detection of axions. Our efforts have focused on developing robust and integrated software tools to model these helioscope components, enabling future assessments of modifications or upgrades to any part of the IAXO axion helioscope and evaluating the potential impact on the experiment's sensitivity. In this manuscript, we demonstrate the application of these tools by presenting a precise signal calculation and response analysis of BabyIAXO's sensitivity to the axion-photon coupling. Though focusing on the Primakoff solar flux component, our virtual helioscope model can be used to test different production mechanisms, allowing for direct comparisons within a unified framework., Comment: 36 pages, 18 figures, 4 tables, Submitted to JHEP

We utilize a hybrid approach that integrates the traditional cross-correlation function (CCF) and machine learning to detect spectroscopic multi-systems, specifically focusing on double-line spectroscopic binary (SB2). Based on the ninth data release (DR9) of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), which includes a medium-resolution survey (MRS) containing 29,920,588 spectra, we identify 27,164 double-line and 3124 triple-line spectra, corresponding to 7096 SB2 candidates and 1903 triple-line spectroscopic binary (SB3) candidates, respectively, representing about 1% of the selection dataset from LAMOST-MRS DR9. Notably, 70.1% of the SB2 candidates and 89.6% of the SB3 candidates are newly identified. Compared to using only the traditional CCF technique, our method significantly improves the efficiency of detecting SB2, saves time on visual inspections by a factor of four., Comment: 18 pages, 11 figures, accepted by ApJS, Data available via China-VO PaperData repository

We present starkiller, an open-source Python package for forward-modeling flux retrieval from integral field unit spectrograph (IFU) datacubes. Starkiller simultaneously provides stellar spectral classification, relative velocity, and line-of-sight extinction for all sources in a catalog, alongside a source-subtracted datacube. It performs synthetic difference imaging by simulating all catalog sources in the field of view, using the catalog for positions and fluxes to scale stellar models, independent of the datacube. This differencing method is particularly powerful for subtracting both point-sources and trailed or even streaked sources from extended astronomical objects. We demonstrate starkiller's effectiveness in improving observations of extended sources in dense stellar fields for VLT/MUSE observations of comets, asteroids and nebulae. We also show that starkiller can treat satellite-impacted VLT/MUSE observations. The package could be applied to tasks as varied as dust extinction in clusters and stellar variability; the stellar modeling using Gaia fluxes is provided as a standalone function. The techniques can be expanded to imagers and to other IFUs., Comment: 21 pages, 16 figures, submitted to AAS journals, open source software

Asteroseismology has emerged as a powerful tool to unravel the intricate relationships between evolved stars and their planetary systems. In this study, we leverage this technique to investigate the evolutionary stages of five exoplanet host stars, each exhibiting solar-like oscillations. Building on our previous work that identified two host stars as red clump and red giant branch (RGB) stars, this study focuses on a new and broader sample. By precisely measuring asteroseismic parameters such as the period spacing of dipole gravity modes ($\Delta\Pi_{1}$), we provide definitive confirmation of these stars' evolutionary states as subgiants or RGB stars. These results are not only crucial for understanding the internal structures of evolved stars but also for predicting the eventual fate of their planetary companions, which may face engulfment as their host stars expand. This research highlights the profound role of asteroseismology in advancing our knowledge of planetary system evolution and opens new pathways for exploring how stellar evolution impacts planetary survival. Our findings set the stage for future studies on the dynamic fates of exoplanets, providing key insights into the intricate processes of stellar and planetary evolution.

The progenitors of Type II-P supernovae (SN) have been confirmed to be red supergiants. However, the upper mass limit of the directly probed progenitors is much lower than that predicted by current theories, and the accurate determination of the progenitor masses is key to understand the final fate of massive stars. Located at a distance of only 6.72 Mpc, the Type II-P SN 2024ggi is one of the closest SN in the last decade. Previous studies have analyzed its progenitor by direct detection, but the derived progenitor mass may be influenced by the very uncertain circumstellar extinction and pulsational brightness variability. In this work, we try to constrain the progenitor mass with an environmental analysis based on images from the Hubble Space Telescope. We found that stars in the progenitor environment have a uniform spatial distribution without significant clumpiness, and we derived the star formation history of the environment with a hierarchical Bayesian method. The progenitor is associated with the youngest population in the SN environment with an age of log($t$/yr) = 7.41 (i.e. 25.7 Myr), which corresponds to an initial mass of $10.2^{+0.06}_{-0.09}$ $M_\odot$. Our work provides an independent measurement of the progenitor mass, which is not affected by circumstellar extinction and pulsational brightness variability., Comment: 8 pages, 3 figures; submitted to ApJL

While mass transfer in binary systems is a crucial aspect of binary evolution models, it remains far from understood. HD 352 is a spectroscopic binary exhibiting ellipsoidal variability, likely due to a tidally deformed giant donor filling its Roche lobe and transferring matter to a faint companion. Here, we analyse VLTI/PIONIER interferometric observations of the system, obtained between 2010 to 2020. We demonstrate that observations near the system's quadrature cannot be explained by simple symmetric disk models, but are consistent with the shape of a Roche-lobe-filling star. We think that this is the first case of tidal deformation of a red giant being observed directly, thanks to the interferometric technique. By combining our interferometric modeling results with the analysis of the optical spectrum, multi-frequency spectral energy distribution, and published radial velocities and light curves, we constrain the system parameters and show that HD 352 will likely soon enter the common envelope phase, although we cannot reject the hypothesis that it is undergoing stable mass transfer against theoretical predictions. This has important consequences for modeling a large class of binary systems. Additionally, our observations confirm that Roche-lobe-filling giants can be resolved with interferometry under favorable conditions. Such observations may help resolve the mass transfer dichotomy in systems like symbiotic binaries, where the predominant mass transfer mode remains unclear., Comment: 10 pages, 6 figures, 2 tables; accepted in Astronomy & Astrophysics

The coldest Y spectral type brown dwarfs are similar in mass and temperature to cool and warm ($\sim$200 -- 400 K) giant exoplanets. We can therefore use their atmospheres as proxies for planetary atmospheres, testing our understanding of physics and chemistry for these complex, cool worlds. At these cold temperatures, their atmospheres are cold enough for water clouds to form, and chemical timescales increase, increasing the likelihood of disequilibrium chemistry compared to warmer classes of planets. JWST observations are revolutionizing the characterization of these worlds with high signal-to-noise, moderate resolution near- and mid-infrared spectra. The spectra have been used to measure the abundances of prominent species like water, methane, and ammonia; species that trace chemical reactions like carbon monoxide; and even isotopologues of carbon monoxide and ammonia. Here, we present atmospheric retrieval results using both published fixed-slit (GTO program 1230) and new averaged time series observations (GO program 2327) of the coldest known Y dwarf, WISE 0855-0714 (using NIRSpec G395M spectra), which has an effective temperature of $\sim$ 264 K. We present a detection of deuterium in an atmosphere outside of the solar system via a relative measurement of deuterated methane (CH$_{3}$D) and standard methane. From this, we infer the D/H ratio of a substellar object outside the solar system for the first time. We also present a well-constrained part-per-billion abundance of phosphine (PH$_{3}$). We discuss our interpretation of these results and the implications for brown dwarf and giant exoplanet formation and evolution., Comment: 17 pages, 7 figures, submitted to the Astrophysical Journal Letters

In this article we study the nature of the recently identified populations of hot companions to red supergiant stars (RSGs). To this end, we compile the literature on the most well characterised systems with the aim of better understanding the hot companions identified with ultra-violet photometry and confirmed with Hubble Space Telescope spectra in the Local Group. We identify 9 systems with current masses greater than around 8 solar masses that have constraints on their orbital periods, which are in the range 3 to 75a. Antares (Alpha Sco) is the obvious outlier in this distribution, having an estimated orbital period of around 2000 a. Mass-ratio (q=M2/MRSG) estimates are available only for 5 of the compiled systems and range between 0.16 < q < 1. Ongoing efforts to identify and characterise hot companions to RSGs in the SMC have revealed 88 hot companions that have observational constraints free of contamination from the RSG component. We present a summary of a recently conducted HST UV spectroscopic survey that aims to characterise a subset of these companions. These companions show a flat-q distribution in the range 0.3 < q < 1., Comment: To appear in the proceedings of the conference LIAC41: The eventful life of massive star multiples

Deuterium fractionation is a well-established evolutionary tracer in low-mass star formation, but its applicability to the high-mass regime remains an open question. The abundances and ratios of deuterated species have often been proposed as reliable evolutionary indicators for different stages of the high-mass star formation. We investigate the role of N$_2$H$^+$ and key deuterated molecules as tracers of the different stages of the high-mass star formation, and test whether their abundance ratios can serve as reliable evolutionary indicators. We conducted APEX observations of o-H$_2$D$^+$ (1$_{10}$-1$_{11}$), N$_2$H$^+$ (4-3), and N$_2$d$^+$ (3-2) in 40 high-mass clumps at different evolutionary stages, selected from the ATLASGAL survey. Molecular column densities ($N$) and abundances ($X$), were derived through spectral line modelling, both under local thermodynamic equilibrium (LTE) and non-LTE conditions. The $N$(o-H$_2$D$^+$) show the smallest deviation from LTE results when derived under non-LTE assumptions. In contrast, N$_2$D$^+$ shows the largest discrepancy between the $N$ derived from LTE and non-LTE. In all the cases discussed, we found that $X$(o-H$_2$D$^+$) decreases more significantly with time than in the case of $X$(N$_2$D$^+$); whereas $X$(N$_2$H$^+$) increases slightly. Therefore, the validity of the recently proposed $X$(o-H$_2$D$^+$)/$X$(N$_2$D$^+$) ratio as a reliable evolutionary indicator was not observed for this sample. While the deuteration fraction derived from N$_2$D$^+$ and N$_2$H$^+$ clearly decreases with clump evolution, the interpretation of this trend is complex, given the different distribution of the two tracers. Our results suggest that a careful consideration of the observational biases and beam-dilution effects are crucial for an accurate interpretation of the evolution of the deuteration process during the high-mass star formation process., Comment: 13 pages, 5 figures, 3 tables. Accepted for publication in A&A

A considerable number of asymptotic giant branch (AGB) stars exhibit UV excess and/or X-ray emission indicative of the presence of a binary companion. AGB stars are so bright that they easily outshine their companions, making almost impossible their identification. Y Gem has been known to be a far-UV- and X-ray-bright AGB star for some decades now, but the nature of its companion, whether it is a main-sequence star or a white dwarf (WD) in a symbiotic system (SySt), is disputed. Our goal is to uncover the true nature of Y Gem to help us peer into the possible misidentified population of SySts. Multi-wavelength IR, optical, UV and X-ray observations have been analyzed to investigate the properties of the stellar components and accretion process in Y Gem. In particular, an optical spectrum of Y Gem is presented here for the first time, while X-ray data are interpreted by means of reflection models produced by an accretion disk and material in its vicinity. The optical spectrum exhibits the typical saw-shaped features of molecular absorptions in addition to narrow recombination and forbidden emission lines. The presence of the emission lines and the analysis of the extinction-corrected UV spectrum suggest a hot component with $T_\mathrm{eff}\approx$60,000 K, $L$=140 L$_{\odot}$, and $R$=0.11 R$_{\odot}$, very likely an accreting WD. The late component is found to be an 1.1 M$_\odot$ AGB star with $T_\mathrm{eff}$=3350 K and $R$=240 R$_\odot$. Using IR, optical, UV, and X-ray data, we found that Y Gem is a S-type SySt whose compact component is accreting at an estimated mass accretion rate of $\dot{M}_\mathrm{acc}=2.3\times10^{-7}$ M$_\odot$ yr$^{-1}$. At such accretion rate, the accreting WD has reached the stable and steady burning phase where no recurrent events are expected., Comment: 13 pages, 9 Figures, 4 Tables; Accepted to A&A

An established method measuring the hydrogen ionisation fraction in shock excited gas is the BE99 method, which utilises six bright forbidden emission lines of [SII]6716, 6731, [NII]6548, 6583, and [OI]6300, 6363. We aim to extent the BE99 method by including more emission lines in the blue and near-infrared part of the spectrum ($\lambda$ = 3500-11000A), and considering higher hydrogen ionisation fractions ($x_e > 0.3$). In addition, we investigate how a non-equilibrium state of the gas and the presence of extinction influence the BE99 technique. We find that plenty additional emission line ratios can in principle be exploited as extended curves (or stripes) in the ($x_e, T_e$)-diagram. If the BE99 equilibrium is reached and extinction is corrected for, all stripes overlap in one location in the ($x_e, T_e$)-diagram indicating the existing gas parameters. The application to the Par Lup 3-4 outflow shows that the classical BE99 lines together with the [NI]5198+5200 lines do not meet in one locationin the ($x_e, T_e$)-diagram. This indicates that the gas parameters derived from the classical BE99 method are not fully consistent with other observed line ratios. A multi-line approach is necessary to determine the gas parameters. From our analysis we derive $n_e \sim$ 45 000 cm^-3 - 53000 cm^-3 , $T_e$ = 7600K - 8000K, and $x_e \sim$ 0.027 - 0.036 for the Par Lup 3-4 outflow. For the 244-440 Proplyd we were able to use the line ratios of [SII]6716+6731, [OI]6300+6363, and [OII]7320, 7330 in the BE99 diagram to estimate the ionisation fraction at knot E3 ($x_e = 0.58 \pm 0.05$). In conclusion, exploiting new line ratios reveals more insights on the state of the gas. Our analysis indicates, however, that a multi-line approach is more robust in deriving gas parameters, especially for high density gas., Comment: 12 pages, 10 figures, accepted to A&A

Magnetohydrodynamic (MHD) waves, playing a crucial role in transporting energy through the solar atmosphere, manifest in various chromospheric structures. Here, we investigated MHD waves in a long-lasting dark fibril using high-temporal-resolution (2~s cadence) Atacama Large Millimeter/submillimeter Array (ALMA) observations in Band 6 (centered at 1.25~mm). We detected oscillations in brightness temperature, horizontal displacement, and width at multiple locations along the fibril, with median periods and standard deviations of $240\pm114$~s, $225\pm102$~s, and $272\pm118$~s, respectively. Wavelet analysis revealed a combination of standing and propagating waves, suggesting the presence of both MHD kink and sausage modes. Less dominant than standing waves, oppositely propagating waves exhibit phase speeds (median and standard deviation of distributions) of $74\pm204$~km/s, $52\pm197$~km/s, and $28\pm254$~km/s for the three observables, respectively. This work demonstrates ALMA's capability to effectively sample dynamic fibrillar structures, despite previous doubts, and provides valuable insights into wave dynamics in the upper chromosphere.

Massive stars lead the ionization and mechanical feedback within young star-forming regions. The Large Magellanic Cloud (LMC) is an ideal galaxy for studying individual massive stars and quantifying their feedback contribution to the environment. We analyze eight exemplary targets in LMC N11 B from the Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) program, using novel spectra from HST (COS and STIS) in the UV, and from VLT (X-shooter) in the optical. We model the spectra of early to late O-type stars by using state-of-the-art PoWR atmosphere models. We determine the stellar and wind parameters (e.g., $T_\star$, $\log g$, $L_{\star}$, $\dot{M}$, $v_\infty$) of the analyzed objects, chemical abundances (C, N, O), ionizing and mechanical feedback ($Q_\mathrm{H}$, $Q_\mathrm{He{\small{I}}}$, $Q_\mathrm{He{\small{II}}}$, $L_\mathrm{mec}$) and X-rays. We report ages of $2-4.5$ Myr and masses of $30-60$ $M_\odot$ for the analyzed stars in N11 B, consistent with a scenario of sequential star formation. We note that the observed wind-momentum luminosity relation is consistent with theoretical predictions. We detect nitrogen enrichment in most of the stars, up to a factor of seven. However, we do not find a correlation between nitrogen enrichment and projected rotational velocity. Finally, based on their spectral type, we estimate the total ionizing photons injected from the O-type stars in N11 B into its environment. We report $\log$ ($\sum$ $Q_\mathrm{H}$)$=50.5$ ph s$^{-1}$, $\log$ ($\sum$ $Q_\mathrm{He{\small{I}}}$)$=49.6$ ph s$^{-1}$ and $\log$ ($\sum$ $Q_\mathrm{He{\small{II}}}$)$=44.4$ ph s$^{-1}$, consistent with the total ionizing budget in N11., Comment: Accepted in Astronomy and Astrophysics

Accurate calculations of starspot spectra are essential for multiple applications in astronomy. The current standard is to represent starspot spectra by spectra of stars that are cooler than the quiet star regions. This implies approximating a starspot as a non-magnetic 1D structure in radiative-convective equilibrium, parametrizing convective energy transport by mixing length theory. It is the inhibition of convection by the starspot magnetic field that is emulated by using a lower spot temperature relative to the quiet stellar regions. Here, we take a different approach avoiding the approximate treatment of convection and instead self-consistently accounting for the interaction between matter, radiation, and the magnetic field. We simulate spots on G2V, K0V, M0V stars with the 3D radiative magnetohydrodynamics code MURaM and calculate spectra ($R \approx 500$ from 250~nm to 6000~nm) using ray-by-ray radiative transfer with the MPS-ATLAS code. We find that the 1D models fail to return accurate umbral and penumbral spectra on K0V and M0V stars where convective and radiative transfer of energy is simultaneously important over a broad range of atmospheric heights rendering mixing length theory inaccurate. However, 1D models work well for G2V stars, where both radiation and convection significantly contribute to energy transfer only in a narrow region near the stellar surface. Quantitatively, the 1D approximation leads to errors longward of 500 nm of about 50\% for both umbral and penumbral flux contrast relative to quiet star regions on M0V stars, and less than 2\% (for umbrae) and 10\% (for penumbrae) for G2V stars., Comment: Accepted for publication in the Astrophysical Journal Letters

New emerging flux (NEF) has long been considered a mechanism for solar eruptions, but detailed process remains an open question. In this work, we explore how NEF drives a coronal magnetic configuration to erupt. This configuration is created by two magnetic sources of strengths $M$ and $S$ embedded in the photosphere, one electric-current-carrying flux rope (FR) floating in the corona, and an electric current induced on the photospheric surface by the FR. The source $M$ is fixed accounting for the initial background field, and $S$ changes playing the role of NEF. We introduce the channel function $C$ to forecast the overall evolutionary behavior of the configuration. Location, polarity, and strength of NEF governs the evolutionary behavior of FR before eruption. In the case of $|S/M|<1$ with reconnection occur between new and old fields, the configuration in equilibrium evolves to the critical state, invoking the catastrophe. In this case, if polarities of the new and old fields are opposite, reconnection occurs as NEF is close to FR; and if polarities are the same, reconnection happens as NEF appears far from FR. With different combinations of the relative polarity and the location, the evolutionary behavior of the system gets complex, and the catastrophe may not occur. If $|S/M|>1$ and the two fields have opposite polarity, the catastrophe always takes place; but if the polarities are the same, catastrophe occurs only as NEF is located far from FR; otherwise, the evolution ends up either with failed eruption or without catastrophe at all., Comment: 18 pages, 7 figures, accepted by ApJL

The Sun is depleted in refractory elements compared to nearby solar twins, which may be linked to the formation of giant or terrestrial planets. Here we present high-resolution, high signal-to-noise spectroscopic data for 17 solar-like stars hosting planets, obtained with Magellan II/MIKE, to investigate whether this depletion is related to planet formation. We derive stellar parameters, including stellar atmosphere, age, radius, mass, and chemical abundances for 22 elements from carbon to europium through line-by-line differential analysis. Our uncertainties range from 0.01 dex for Fe and Si to 0.08 dex for Sr, Y, and Eu. By comparing the solar abundances to those of the 17 stars, we investigate the differential abundance ([X/Fe]$_{\rm solar}$ - [X/Fe]$_{\rm star}$) versus condensation temperature ($T_c$) trend. In particular, we apply Galactic chemical evolution corrections to five solar twins within the full sample. Our results conform to previous studies that the Sun is relatively depleted in refractory compared to volatile elements. For both five solar twins and the rest of solar-like stars, we find that all stars hosting known gas giant planets exhibit negative $T_c$ trend slopes, suggesting that the Sun is relatively depleted in refractory elements compared to similar giant-planet-host stars. Additionally, we find no correlation between $T_c$ trend slopes and the total mass of detected terrestrial planets in each system, suggesting that terrestrial planet formation may not be the cause of refractory element depletion in the Sun., Comment: 26 pages, 10 figures, 7 tables; accepted for publication in ApJ

The recent discovery of planetesimals orbiting white dwarfs has renewed interest in the final chapters of the evolution of planetary systems. Although observational and theoretical studies have examined the dynamical evolution of these systems, studies of their magnetic star-planet interactions, as powered by unipolar induction, have thus far only been assessed theoretically. This fifth paper of the ROME (Radio Observations of Magnetized Exoplanets) series presents the results of a targeted mini-survey of nine white dwarfs within 25 pc without known stellar mass companions in search of radio emissions generated by magnetic interactions between white dwarfs and their planetary remnants. This $\sim$5 GHz Arecibo radio telescope survey achieved mJy-level sensitivity over $<$1 s integration times. Although no exoplanet-induced stellar radio flares were detected, this is the first survey to search for magnetic star-planet interactions between white dwarfs and planetary companions, cores, or disrupted planetesimals. It is also the most extensive and sensitive radio survey for intrinsic coronal emissions from apparently isolated white dwarfs. The study of radio emissions from white dwarf systems may present a new means to detect and measure DC and DQ magnetic fields, search for white dwarf coronae, characterize the density and spatial distribution of white dwarf magnetospheric plasma, characterize the dynamical and electrical properties of planetary cores, and offer new constraints on the modeling of double degenerate merger events., Comment: 14 pages, 3 tables, 2 figures. Summary presented at the 56th Annual Meeting of the AAS Division for Planetary Sciences, abstract 407.04. In press at The Astrophysical Journal

We investigate observations of circumbinary disks (CBD), to find evidence for an equilibrium eccentricity predicted by current binary accretion theory. Although stellar binary demographics in the Milky Way show no evidence for a preferred eccentricity for binary systems, we show that actively accreting systems lie on a predicted equilibrium eccentricity curve. We constrain our sample to only systems that have well defined orbital parameters (e.g,. eccentricity, mass-ratio, inclination angle). We find observations are consistent with theory for stellar binaries that are aligned with the disk and that are separated enough that tidal circularization is negligible. This suggests that eccentricity in these systems evolves after the dissipation of the CBD, given the flat eccentricity distribution of binary systems in the Milky Way., Comment: 5 pages, 2 figures

Throughout a planetary system's formation evolution, some of the planetary material may end up falling into the host star and be engulfed by it, leading to a potential variation of the stellar composition. The present study explores how planet engulfment may impact the chemical composition of the stellar surface and discusses what would be the rate of events with an observable imprint, for Sun-like stars. We use data from the NGPPS calculations by the Generation III Bern model to analyse the conditions under which planet engulfment may occur. Additionally, we use stellar models computed with Cesam2k20 to account for how the stellar internal structure and its processes may affect the dilution of the signal caused by planet engulfment. Our results show that there are three different phases associated to different mechanisms under which engulfment events may happen. Moreover, systems that undergo planet engulfment are more likely to come from protoplanetary disks that are more massive and more metal-rich than non-engulfing systems. Engulfment events leading to an observable signal happen after the dissipation of the protoplanetary disk when the convective envelope of the stars becomes thinner. With the stellar convective layer shrinking as the star evolves in the main sequence, they display a higher variation of chemical composition, which also correlates with the amount of engulfed material. By accounting for the physical processes happening in the stellar interior and in the optimistic case of being able to detect variations above 0.02 dex in the stellar composition, we find an engulfment rate no higher than $20\%$ for Sun-like stars that may reveal detectable traces of planet engulfment. Engulfment events that lead to an observable variation of the stellar composition are rare due to the specific conditions required to result in such signatures., Comment: Accepted to A&A journal

Context: L$_2$ Puppis (L$_2$ Pup) is a nearby red giant star and an important object in late-type star research because it has a dust disc and potentially a companion. Aims: L$_2$ Pup is often called the second-closest Asymptotic Giant Branch (AGB) star to the sun, second only to R Doradus. However, whether the star is indeed on the AGB or the Red Giant Branch (RGB) is questionable. We review its evolutionary state. Methods: We analyse high-resolution optical archive spectra to search for absorption lines of the third dredge-up indicator technetium (Tc) in L$_2$~Pup. We also compare the star to a sample of well-known AGB stars in terms of luminosity and pulsation properties and place it in a Gaia-2MASS diagram. Results: L$_2$ Pup is found to be Tc-poor. Thus, it is not undergoing third dredge-up events. The star is fainter than the RGB tip and fainter than all Tc-rich stars in the comparison sample. L$_2$ Pup pulsates in the fundamental mode, similar to Mira variables, but its pulsation properties do not allow us to distinguish between the RGB and AGB stages. Conclusions: In conclusion, L$_2$ Pup could be an RGB or early AGB star, but it is more likely to be an RGB than an AGB star. Our findings are important for a better understanding of the L$_2$ Pup system and its past and future evolution, Comment: 8 pages, 4 figures, accepted for publication in Astronomy & Astrophysics

The Kepler and K2 missions enabled robust calculations of planet occurrence rates around FGKM-type stars. However, these missions observed too few stars with earlier spectral types to tightly constrain the occurrence rates of planets orbiting hotter stars. Using TESS, we calculate the occurrence rate of small ($1 \, R_\oplus < R_{\rm p} < 8 \, R_\oplus$), close-in ($P_{\rm orb} < 10$~days) planets orbiting A-type stars for the first time. We search a sample of 20,257 bright ($6 < T < 10$) A-type stars for transiting planets using a custom pipeline and vet the detected signals, finding no reliable small planets. We characterize the pipeline completeness using injection/recovery tests and determine the $3\sigma$ upper limits of the occurrence rates of close-in sub-Saturns ($4 \, R_\oplus < R_{\rm p} < 8 \, R_\oplus$), sub-Neptunes ($2 \, R_\oplus < R_{\rm p} < 4 \, R_\oplus$), and super-Earths ($1 \, R_\oplus < R_{\rm p} < 2 \, R_\oplus$). We find upper limits of $2.2 \pm 0.4$ sub-Saturns and $9.1 \pm 1.8$ sub-Neptunes per 1000 A-type stars, which may be more than $3\times$ and $6\times$ lower than \textit{Kepler}-era estimates for Sun-like stars. We calculate an upper limit of $186 \pm 34$ super-Earths per 1000 A-type stars, which may be more than $1.5\times$ lower than that for M dwarfs. Our results hint that small, close-in planets become rarer around early-type stars and that their occurrence rates decrease faster than that of hot Jupiters with increasing host star temperature. We discuss plausible explanations for these trends, including star-disk interactions and enhanced photoevaporation of planet atmospheres., Comment: 35 pages, 11 figures, 1 table

While WISE is the largest, best quality infrared all-sky survey to date, a smaller coverage mission, Spitzer, was designed to have better sensitivity and spatial resolution at similar wavelengths. Confusion and contamination in WISE data result in discrepancies between them. We present a novel approach to work with WISE measurements with the goal of maintaining both its high coverage and vast amount of data while taking full advantage of the higher sensitivity and spatial resolution of Spitzer. We have applied machine learning (ML) techniques to a complete WISE data sample of open cluster members, using a training set of paired data from high-quality Spitzer Enhanced Imaging Products (SEIP), MIPS and IRAC, and allWISE catalogs, W1 (3.4 {\mu}m) to W4 (22 {\mu}m) bands. We have tested several ML regression models with the aim of predicting mid-infrared fluxes at MIPS1 (24 {\mu}m) and IRAC4 (8 {\mu}m) bands from WISE fluxes and quality flags. In addition, to improve the prediction quality, we have implemented feature selection techniques to remove irrelevant WISE variables. We have notably enhanced WISE detection capabilities, mostly at lowest magnitudes, which previously showed the largest discrepancies with Spitzer. In our particular case, extremely randomized trees was found to be the best algorithm to predict mid-infrared fluxes from WISE variables. We have tested our results in the SED of members of IC 348. We show discrepancies in the measurements of Spitzer and WISE and demonstrate the good concordance of our predicted fluxes with the real ones. ML is a fast and powerful tool that can be used to find hidden relationships between datasets, as the ones that exist between WISE and Spitzer fluxes. We believe this approach could be employed for other samples from the allWISE catalog with SEIP positional counterparts, and in other astrophysical studies with analogous discrepancies., Comment: 13 pages, 10 figures

Data access -- or the availability of new and archival data for use by the larger community -- is key for scientific advancement. How data is presented, searched, and formatted determines accessibility and it can be difficult to find a solution that fits the needs of a given subdiscipline. We present a generalized roadmap for developing a specialty astronomy database with web application based on the development of the SpExoDisks (Spectra of Exoplanet forming Disks) database (spexodisks.com), which provides infrared spectra of protoplanetary disks. Expertise in an astronomy subdiscipline can provide two necessary components for creating a database: access to a large volume of specialized data and knowledge of how that data should be presented to the community. However, there are a variety of steps and decisions for database development that can fall outside astronomy expertise. Here we offer generalized discussions on design and process that are accompanied by real-world examples from the SpExoDisks developer team and website. Starting from the database portal design and data organization, we demonstrate on-demand data distribution and query using publicly accessible database software. These systems support interactive visualizations such that users can explore spectra directly from their browsers. We also offer details that show how the technical concepts in SpExoDisks are implemented, particularly emphasizing sustainability and long-term management of the codebase and processes. Finally, we illustrate the utility that a specialty website can offer to the community by providing a specific example of how the combined spectra from SpExoDisks can enhance our understanding of protoplanetary disks., Comment: 25 pages, 2 appendices, 6 figures, 1 table, accepted to PASP

Understanding the interstellar and potentially circumstellar extinction in the sight lines of classical T Tauri stars is an important ingredient for constructing reliable spectral energy distributions, which catalyze protoplanetary disk chemistry, for example. Therefore, some attempts of measuring $A_{V}$ toward individual stars have been made using partly different wavelength regimes and different underlying assumptions. We used strong lines of Ly{\alpha} fluorescent H2 and derived the extinction based on the assumption of optically thin transitions. We investigated a sample of 72 classical T Tauri stars observed with the Hubble Space Telescope in the framework of the ULLYSES program. We computed $A_{V}$ and $R_{V}$ values for the 34 objects with sufficient data quality and an additionally $A_{V}$ value for the canonical $R_{V}$ = 3.1 value. Our results agree largely with values obtained from optical data. Moreover, we confirm the degeneracy between $A_{V}$ and $R_{V}$ and present possibilities to break this. Finally, we discuss whether the assumption of optical thin lines is valid., Comment: 12 pages, 6 figures, accepted to A&A

Energetic Neutral Atom (ENA) observations provide valuable insights into the plasma conditions in the heliosphere and the surrounding interstellar medium. Unlike plasma detectors, which measure charged particles tied to the magnetic fields at their location, ENA detectors capture former ions that were neutralized in distant regions and traverse the heliosphere in straight trajectories. ENA fluxes near the Sun represent line-of-sight integrals of parent ion fluxes multiplied by neutralization (production) rates and reduced by the probability of ENA reionization (loss) processes. So far, most ENA analyses have focused on charge exchange between hydrogen atoms and protons as the primary source of ENAs. Here, we examine various ENA production and loss processes throughout the heliosphere in the broad energy range (5 eV to 500 keV) encompassing the next-generation ENA instruments aboard the Interstellar Mapping and Acceleration Probe (IMAP) mission. Our study considers binary collisions involving the most abundant species: protons, electrons, {\alpha}-particles, He+ ions, photons, as well as hydrogen and helium atoms. Our findings indicate that, in addition to ENAs produced by charge exchange of energetic protons with hydrogen atoms, a significant portion of high-energy ENAs originate from the charge exchange with helium atoms. Below 10 keV, the dominant ENA loss processes are charge exchange collisions with protons and photoionization. However, stripping ionization processes, e.g., from collisions with ambient interstellar neutral hydrogen, become the main loss mechanism for higher energies because the charge exchange rate rapidly decreases., Comment: 22 pages, 7 figures, 3 tables, submitted to ApJS

The knowledge of beta decay transitional probabilities and GamowTeller (GT) strength functions from highly excited states of nuclides is of particular importance for applications to astrophysical network calculations of nucleosynthesis in explosive stellar events. These quantities are challenging to achieve from measurements or computations using various nuclear models. Due to unavailability of feasible alternatives, many theoretical studies often rely on the Brink Axel (BA) hypothesis, that is, the response of strength functions depends merely on the transition energy of the parent nuclear ground state and is independent of the underlying details of the parent state, for the calculation of stellar rates. BA hypothesis has been used in many applications from nuclear structure determination to nucleosynthesis yield in the astrophysical matter., Comment: 29 Pages, 8Figures 5 Tables

The typically large distances, extinction, and crowding of Galactic supermassive star clusters have so far hampered the identification of their very low mass members, required to extend our understanding of star and planet formation, and early stellar evolution, to starburst. This situation has now evolved thanks to the James Webb Space Telescope (JWST), and its unmatched resolution and sensitivity in the infrared. In this paper, the third of the series of the Extended Westerlund 1 and 2 Open Clusters Survey (EWOCS), we present JWST/NIRCam and JWST/MIRI observations of the supermassive star cluster Westerlund 1. These observations are specifically designed to unveil the cluster members down to the BD mass regime, and to allow us to select and study the protoplanetary disks and to study the mutual feedback between the cluster members and the surrounding environment. Westerlund 1 was observed as part of JWST GO-1905 for 23.6 hours. The data have been reduced using the JWST calibration pipeline, together with specific tools necessary to remove artifacts. Source identification and photometry were performed with DOLPHOT. The MIRI images show a plethora of different features. Diffuse nebular emission is observed around the cluster, which is typically composed of myriads of droplet-like features pointing toward the cluster center or the group of massive stars surrounding the WR star W72/A. A long pillar is also observed in the NW. The MIRI images also show resolved shells and outflows surrounding the M-type RSG W20, W26, W75, and W237, the sgB[e] star W9 and the YHG W4. The color-magnitude diagrams built using the NIRCam photometry show a clear cluster sequence, which is marked in its upper part by the 1828 NIRCam stars with X-ray counterparts. NIRCam observations using the F115W filter have reached the 23.8 mag limit with 50\% completeness (roughly corresponding to a 0.06 Msol brown dwarf)., Comment: Accepted for publication in Astronomy & Astrophysics

The disk mass and substructure in young stellar objects suggest that planet formation may start at the protostellar stage through the growth of dust grains. To accurately estimate the grain size at the protostellar stage, we have observed the Class I protostar TMC-1A using the Jansky Very Large Array (VLA) at the Q (7 mm) and Ka (9 mm) bands at a resolution of ~0.2" and analyzed archival data of Atacama Large Millimeter/submillimeter Array (ALMA) at Band 6 (1.3 mm) and 7 (0.9 mm) that cover the same spatial scale. The VLA images show a compact structure with a size of ~25 au and a spectral index of ~2.5. The ALMA images show compact and extended structures with a spectral index of ~2 at the central ~40 au region and another index of ~3.3 in the outer region. Our SED analysis using the observed fluxes at the four bands suggests one branch with a small grain size of ~0.12 mm and another with a grown grain size of ~4 mm. We also model polarized dust continuum emission adopting the two grain sizes and compare them with an observational result of TMC-1A, suggesting that the small grain size is preferable to the grown grain size. The small grain size implies gravitational instability in the TMC-1A disk, which is consistent with a spiral-like component recently identified., Comment: 18 pages, 9 figures, accepted by the Astrophysical Journal

This study investigates the temporal and spatial variations in lithium abundance within the Milky Way using a sample of 22,034 main-sequence turn-off (MSTO) stars and subgiants, characterised by precise stellar ages, 3D NLTE (non-local thermodynamic equilibrium) lithium abundances, and birth radii. Our results reveal a complex variation in lithium abundance with stellar age: a gradual increase from 14 Gyr to 6 Gyr, followed by a decline between 6 Gyr and 4.5 Gyr, and a rapid increase thereafter. We find that young Li-rich stars (ages $<$ 4 Gyr, A(Li) $>$ 2.7 dex) predominantly originate from the outer disc. By binning the sample according to guiding center radius and z$_{\rm max}$, we observe that these young Li-rich stars migrate radially to the local and inner discs. In addition, the stars originating from the inner disc experienced a rapid Li enrichment process between 8 Gyr and 6 Gyr. Our analysis suggests that the age range of Li-dip stars is 4-5 Gyr, encompassing evolution stages from MSTO stars to subgiants. The Galactic radial profile of A(Li) (with respect to birth radius), as a function of age, reveals three distinct periods: 14-6 Gyr ago, 6-4 Gyr ago, and 4-1 Gyr ago. Initially, the lithium abundance gradient is positive, indicating increasing Li abundance with birth radius. During the second period, it transitions to a negative and broken gradient, mainly affected by Li-dip stars. In the final period, the gradient reverts to a positive trend., Comment: Accepted for publication in MNRAS. 9 pages, 10 figures

The magnetic dynamo mechanism of giant stars remains an open question, which can be explored by investigating their activity-rotation relations with multiple proxies. By using the data from the LAMOST and \emph{GALEX} surveys, we carried out a comprehensive study of activity-rotation relations of evolved stars based on \cahk lines, $\rm{H\alpha}$ lines and near ultraviolet (NUV) emissions. Our results show that evolved stars and dwarfs obey a similar power-law in the unsaturated region of the activity-rotation relation, indicating a common dynamo mechanism in both giant and dwarfs. There is no clear difference in the activity levels between red giant branch stars and red clump stars, nor between single giants and those in binaries. Additionally, our results show that the NUV activity levels of giants are comparable to those of G- and K-type dwarfs and are higher than those of M dwarfs., Comment: APJ accepted

Recent literature reports a color deviation between observed Gaia color-magnitude diagrams (CMDs) and theoretical model isochrone predictions, particularly in the very low-mass regime. To assess its impact on cluster age determination via isochrone fitting, we quantified the color deviations for three benchmark clusters, Hyades, Pleiades, and Praesepe, both for the Gaia color (BP-RP) and (G-RP). In general, the (G-RP) color deviations are smaller than the (BP-RP) ones. Empirical color correction functions based on these benchmarks are derived for the currently available MIST and PARSEC 1.2S isochrone models. Applying the correction functions to 31 additional open clusters and 3 moving groups results in a significantly improved alignment between the isochrones and observed CMDs. With our empirical corrections, isochrones provide age estimates consistent with literature values obtained through the spectral Lithium Depletion Boundary method, validating the effectiveness of our approach. The corresponding metallicities with PARSEC 1.2S also show a good agreement with the spectroscopic results. The empirical color correction function we present in this work offers a tool for a consistent age determination within the full mass range of stellar clusters using the isochrone fitting method., Comment: 27 pages. Accepted for publication in ApJ

Star-forming regions host a large and evolving suite of molecular species. Molecular transition lines, particularly of complex molecules, can reveal the physical and dynamical environment of star formation. We aim to study the large-scale structure and environment of high-mass star formation through single-dish observations of CH$_3$CCH, CH$_3$OH, and H$_2$CO. We have conducted a wide-band spectral survey with the IRAM 30-m telescope and the 100-m GBT towards the high-mass star-forming region DR21(OH)/N44. We use a multi-component local thermodynamic equilibrium model to determine the large-scale physical environment near DR21(OH) and the surrounding dense clumps. We follow up with a radiative transfer code for CH$_3$OH to look at non-LTE behaviour. We then use a gas-grain chemical model to understand the formation routes of these molecules in their observed environments. We disentangle multiple components of DR21(OH) in each of the three molecules. We find a warm and cold component each towards the dusty condensations MM1 and MM2, and a fifth broad, outflow component. We also reveal warm and cold components towards other dense clumps in our maps: N40, N36, N41, N38, and N48. We find thermal mechanisms are adequate to produce the observed abundances of H$_2$CO and CH$_3$CCH while non-thermal mechanisms are needed to produce CH$_3$OH. Through a combination of wide-band mapping observations, LTE and non-LTE model analysis, and chemical modelling, we disentangle the different velocity and temperature components within our clump-scale beam, a scale that links a star-forming core to its parent cloud. We find numerous warm, 20-80 K components corresponding to known cores and outflows in the region. We determine the production routes of these species to be dominated by grain chemistry., Comment: Submitted to A&A, 21 pages, 23 figures

Planets lose mass to atmospheric outflows, and this mass loss is thought to be central in shaping the bimodal population of gaseous giant and rocky terrestrial exoplanets in close orbits. We model the escape of planetary atmospheres in three dimensional gas dynamic simulations in order to study their emergent morphology. Planetary outflows show a range of shapes from fast, isotropic outflows bounded by bow shocks to slower motion confined to thin streams. We show that a crucial factor is the role of the tidal gravity and orbiting reference frame in which planets lose mass. Flows can be characterized by the dimensionless Rossby number evaluated at the scale of the Hill sphere. Flows with a low Rossby number are significantly deviated and shaped by the stellar gravity, while those with a high Rossby number are comparatively unaffected. Rossby number alone is sufficient to predict outflow morphology as well as kinematic gradients across transit. The known exoplanet population should span a range of outflow Rossby numbers and thus shapes. We can use this information to constrain outflow physics and to inform observing strategies., Comment: Submitted to AAS journals, comments welcome! Reproduction materials at: https://github.com/morganemacleod/PWMorphology

The Medium-Resolution Spectrometer on the Mid-Infrared Instrument on JWST obtained spectra of three carbon stars in the Large Magellanic Cloud. Two of the spectra differ significantly from spectra obtained ~16-19 years earlier with the Infrared Spectrograph on the Spitzer Space Telescope. The one semi-regular variable among the three has changed little. The long-period Mira variable in the sample shows changes consistent with its pulsation cycle. The short-period Mira shows dramatic changes in the strength of its molecular absorption bands, with some bands growing weaker and some stronger. Whether these variations result from its pulsation cycle or its evolution is not clear., Comment: 7 pages, 3 figures

Shock-generated transients, such as hot flow anomalies (HFAs), upstream of planetary bow shocks, play a critical role in electron acceleration. Using multi-mission data from NASA's Magnetospheric Multiscale (MMS) and ESA's Cluster missions, we demonstrate the transmission of HFAs through Earth's quasi-parallel bow shock, associated with acceleration of electrons up to relativistic energies. Energetic electrons, initially accelerated upstream, are shown to remain broadly confined within the transmitted transient structures downstream, where betatron acceleration further boosts their energy due to elevated compression levels. Additionally, high-speed jets form at the compressive edges of HFAs, exhibiting a significant increase in dynamic pressure and potentially contributing to driving further localized compression. Our findings emphasize the efficiency of quasi-parallel shocks in driving particle acceleration far beyond the immediate shock transition region, expanding the acceleration region to a larger spatial domain. Finally, this study underscores the importance of multi-scale observational approach in understanding the convoluted processes behind collisionless shock physics and their broader implications.

Empirical exoplanet mass-radius relations have been used to study the demographics and compositions of small exoplanets for many years. However, the heterogeneous nature of these measurements hinders robust statistical analysis of this population, particularly with regard to the masses of planets. For this reason, we perform a homogeneous and consistent re-analysis of the radial velocity (RV) observations of 85 small exoplanets using publicly available HARPS RV data and the fitting toolkit Pyaneti. For the entire sample, we run 12 different models to investigate the impact of modelling choices, including the use of multi-dimensional Gaussian Processes (GPs) to mitigate stellar activity. We find that the way orbital eccentricity is modelled can significantly impact the RV amplitude found in some cases. We also find that the addition of a GP to mitigate stellar activity does impact the RV amplitude found - though if the GP is modelled on activity indicators as well as the RVs the results are more robust. The RV amplitude found for every planet in our sample using all the models is made available for other groups to perform demographics studies. Finally, we provide a list of recommendations for the RV community moving forward., Comment: Accepted for publication in A&A

An accurate estimation of the continuum excess emission from accretion spots and inner circumstellar disk regions is crucial for a proper derivation of fundamental stellar parameters in accreting systems. However, the presence of starspots can make disentangling the complicated multi-component emission in these systems challenging. Subtraction of a single-temperature spectral template is insufficient to account for the composite stellar emission, as we demonstrated in a recent campaign involving Weak-Lined T Tauri Stars. Here, we model the moderate resolution near-infrared spectra of Classical T Tauri Stars, presenting new spectral models that incorporate spotted stars plus emission from accretion hot-spots and a warm inner disk, allowing us to simultaneously reconstruct the entire 0.8-2.4 micrometer spectrum of our sixteen targets. Using these models, we re-derive the continuum excess emission. Our results indicate that accounting for starspots resolves the need to include a previously proposed intermediate temperature component in the IYJ excess, and highlights the importance of a proper treatment of starspots in studies of accreting low-mass stars., Comment: Accepted for publication in ApJ

Understanding and predicting the structure and evolution of coronal mass ejections (CMEs) in the heliosphere remains one of the most sought-after goals in heliophysics and space weather research. A powerful tool for improving current knowledge and capabilities consists of multi-spacecraft observations of the same event, which take place when two or more spacecraft fortuitously find themselves in the path of a single CME. Multi-probe events can not only supply useful data to evaluate the large-scale of CMEs from 1D in-situ trajectories, but also provide additional constraints and validation opportunities for CME propagation models. In this work, we analyse and simulate the coronal and heliospheric evolution of a slow, streamer-blowout CME that erupted on 23 September 2021 and was encountered in situ by four spacecraft approximately equally distributed in heliocentric distance between 0.4 and 1 au. We employ the Open Solar Physics Rapid Ensemble Information (OSPREI) modelling suite in ensemble mode to predict the CME arrival and structure in a hindcast fashion and to compute the "best-fit" solutions at the different spacecraft individually and together. We find that the spread in the predicted quantities increases with heliocentric distance, suggesting that there may be a maximum (angular and radial) separation between an inner and an outer probe beyond which estimates of the in-situ magnetic field orientation (parameterised by flux rope model geometry) increasingly diverge. We discuss the importance of these exceptional observations and the results of our investigation in the context of advancing our understanding of CME structure and evolution as well as improving space weather forecasts., Comment: 20 pages, 12 figures, 3 tables; accepted for publication in MNRAS

Thermal changes in coronal loops are well-studied, both in quiescent active regions and in flaring scenarios. However, relatively little attention has been paid to loop emission in the hours before the onset of a solar flare; here, we present the findings of a study of over 50 off-limb flares of GOES class C5.0 and above. We investigated the integrated emission variability for Solar Dynamics Observatory Atmospheric Imaging Assembly channels 131, 171, 193, and 304 \r{A}ngstroms for 6 hours before each flare, and compared these quantities to the same time range and channels above active regions without proximal flaring. We find significantly increased emission variability in the 2-3 hours before flare onset, particularly for the 131 and 304 channels. This finding suggests a potential new flare prediction methodology. The emission trends between the channels are not consistently well-correlated, suggesting a somewhat chaotic thermal environment within the coronal portion of the loops that disturbs the commonly-observed heating and cooling cycles of quiescent active region loops. We present our approach, the resulting statistics, and discuss the implications for heating sources in these pre-flaring active regions., Comment: Accepted for publication to the Astrophysical Journal Letters

It is widely assumed that dust opacities in molecular clouds follow a power-law profile with an index, $\beta$. Recent studies of the Orion Molecular Cloud (OMC) 2/3 complex, however, show a flattening in the spectral energy distribution (SED) at $ \lambda > 2$ mm implying non-constant indices on scales $\gtrsim$ 0.08 pc. The origin of this flattening is not yet known but it may be due to the intrinsic properties of the dust grains or contamination from other sources of emission. We investigate the SED slopes in OMC 2/3 further using observations of six protostellar cores with NOEMA from 2.9 mm to 3.6 mm and ALMA-ACA in Band 4 (1.9 -- 2.1 mm) and Band 5 (1.6 -- 1.8 mm) on core and envelope scales of $\sim 0.02 - 0.08$ pc. We confirm flattened opacity indices between 2.9 mm and 3.6 mm for the six cores with $\beta \approx -0.16 - 1.45$, which are notably lower than the $\beta$ values of $> 1.3$ measured for these sources on $0.08$ pc scales from single-dish data. Four sources have consistent SED slopes between the ALMA data and the NOEMA data. We propose that these sources may have a significant fraction of emission coming from large dust grains in embedded disks, which biases the emission more at longer wavelengths. Two sources, however, had inconsistent slopes between the ALMA and NOEMA data, indicating different origins of emission. These results highlight how care is needed when combining multi-scale observations or extrapolating single-band observations to other wavelengths., Comment: Accepted to ApJ, 18 pages, 4 figures, 6 Tables

We present a single-fluid approach for the simulation of partially-ionized plasmas (PIPs) which is designed to capture the non-ideal effects introduced by neutrals while remaining close in computational efficiency to single-fluid MHD. This is achieved using a model which treats the entire partially-ionized plasma as a single mixture, which renders internal ionization/recombination source terms unnecessary as both the charged and neutral species are part of the mixture's conservative system. Instead, the effects of ionization and the differing physics of the species are encapsulated as material properties of the mixture. Furthermore, the differing dynamics between the charged and neutral species is captured using a relative-velocity quantity, which impacts the bulk behavior of the mixture in a manner similar to the treatment of the ion-electron relative-velocity as current in MHD. Unlike fully-ionized plasmas, the species composition of a PIP changes rapidly with its thermodynamic state. This is captured through a look-up table referred to as the tabulated equation of state (TabEoS), which is constructed prior to runtime using empirical physicochemical databases and efficiently provides the ionization fraction and other material properties of the PIP specific to the thermodynamic state of each computational cell. Crucially, the use of TabEoS also allows our approach to self-consistently capture the non-linear feedback cycle between the PIP's macroscopic behavior and the microscopic physics of its internal particles, which is neglected in many fluid simulations of plasmas today.

Consistent stellar evolution and nonlinear radial stellar pulsation calculations were carried out for models of asymptotic giant branch stars with initial masses $1.5M_\odot\le M_\mathrm{ZAMS}\le 3M_\odot$ and initial metal abundance $Z=0.006$. All the models are shown to be either the fundamental mode or the first overtone pulsators. The lower limit of the first overtone period increases with increasing mass of the Mira model from $\Pi_{1,\min}\approx 80$ days for $M=1.3M_\odot$ to $\Pi_{1,\min}\approx 120$ days for $M=2.6M_\odot$. The upper limit of the first overtone period and lower limit of the fundamental mode period depend on the stellar structure during mode switching and range from $\Pi_{1,\max}=130$, $\Pi_{0,\min}=190$ days for $M=0.96M_\odot$ to $\Pi_{1,\max}=210$, $\Pi_{0,\min}=430$ days for $M=2.2M_\odot$. The slope of the theoretical period--luminosity relation of Mira variables perceptibly increases with decreasing $Z$. Fourier spectra of the kinetic energy of twelve hydrodynamic models show a split of the fundamental mode maximum into several equidistant components. Frequency intervals between split components fall within the range $0.03 \le \Delta\nu/\nu_0 \le 0.1$. The superposition of radial oscillations with the fundamental mode splitting leads to the long-term amplitude variations with the cycle length from 10 to 30 times longer than the fundamental mode period. A more thorough analysis of hydrodynamic models is required for understanding the origin of the principal pulsation mode splitting., Comment: 16 pages, 8 figures, accepted to Astron. Letters

The gamma-ray binary HESS J0632+057 consists of a Be star and an undetected compact object in a $\sim$317 day orbit. The interpretation of the emission from this system is complicated by the lack of a clear orbital solution, as two different and incompatible orbital solutions were obtained by previous radial velocity studies of this source. In order to address this, we report on 24 new observations, covering $\sim$60 per cent of the orbit which we have undertaken with the Southern African Large Telescope (SALT). We obtained new radial velocity measurements from cross-correlation of the narrower spectral features, and by fitting Voigt profiles to the wings of the Balmer emission lines. Additionally, we find an indication of orbital variability in the equivalent widths and V/R of the Balmer lines. Using the combined data from this study, as well as archival data from the earlier radial velocity studies, we have derived updated orbital solutions. Using reported H $\alpha$ emission radial velocities - previously not considered for the orbital solution - along with the new SALT data, a solution is obtained where the brighter peak in the X-ray and gamma-ray light curves is closer to periastron. However, continuing sparse coverage in the data around the expected phases of periastron indicates that the orbital solution could be improved with further observation., Comment: 8 pages, 6 figures

As mixed with real pulsations, the reflection of super-Nyquist frequencies (SNFs) pose a threat to asteroseismic properties. Although SNFs have been studied in several pulsating stars, a systematic survey remains scarcely explored. Here we propose a method to identify SNFs from Kepler and TESS photometry by characterizing their periodic frequency modulations using a sliding Fourier transform. After analyzing long cadence photometry in the Kepler legacy, we have identified 304 SNFs in 56 stars from 45607 frequencies in $\sim600$ $\gamma$ Doradus stars, corresponding to a fraction of approximately $0.67\%$ and $9.2\%$, respectively. Most SNFs are detected in the frequency range of pressure mode over 120 $\mu$Hz and the fraction of SNF detection increases as frequency up to $\sim7\%$. We barely found two potential SNFs mixed with gravity modes in two $\gamma$ Doradus stars. These findings indicate that SNFs have a negligible impact on the global seismic properties, such as those derived from period spacing in $\gamma$ Doradus stars. However, we stress that SNFs must be carefully and systematically examined by this method in other pulsating stars, particularly $\delta$ Scuti and hot B subdwarf stars, to establish a solid foundation for precise asteroseismolgy of various types of pulsators., Comment: 7 pages, 6 figures, accepted for publication in A&A

In the astrophysics community it is common practice to model collisionless dust, entrained in a gas flow, as a pressureless fluid. However a pressureless fluid is fundamentally different from a collisionless fluid - the latter of which generically possess a non-zero anisotropic pressure or stress tensor. In this paper we derive a fluid model for collisionless dust, entrained in a turbulent gas, starting from the equations describing the motion of individual dust grains. We adopt a covariant formulation of our model to allow for the geometry and coordinate systems prevalent in astrophysics, and provide a closure valid for the accretion disc context. We show that the continuum mechanics properties of a dust fluid corresponds to a higher-dimensional anisotropic Maxwell fluid, after the extra dimensions are averaged out, with a dynamically important rheological stress tensor. This higher-dimensional treatment has the advantage of keeping the dust velocity and velocity of the fluid seen, and their respective moments, on the same footing. This results in a simplification of the constitutive relation describing the evolution of the dust Rheological stress., Comment: 51 pages, 7 figures, accepted for publication in JfM

Inferring the coupling of different atmospheric layers requires observing spectral lines sensitive to the atmospheric parameters, particularly the magnetic field vector, at various heights. The best way to tackle this goal is to perform multi-line observations simultaneously. For instance, the new version of the Gregor Infrared Spectrograph instrument offers the possibility to observe the spectral lines at 8542 and 10830 A simultaneously for the first time. The first spectral window contains the Ca II 8542 A spectral line, while the Si I 10827 A transition and He I 10830 A triplet infrared lines can be found in the second spectral window. As the sensitivity to the atmospheric parameters and the height of formation of those transitions is different, combining them can help understand the properties of the solar photosphere and chromosphere and how they are magnetically coupled. Traditionally, the analysis of the Si I 10827 A transition assumes local thermodynamic equilibrium (LTE), which is not the best approximation to model this transition. Hence, in this work, we examine the potential of performing non-LTE (NLTE) inversions of the full Stokes vector of the Si I 10827 A spectral line. The results indicate that we properly infer the atmospheric parameters through an extended range of atmospheric layers in comparison with the LTE case (only valid for the spectral line wings, i.e., the low photosphere), with no impact on the robustness of the solution and just a minor increase in computational time. Thus, the NLTE assumption will help to accurately constrain the photospheric physical parameters when performing combined inversions with, e.g., the Ca II 8542 A spectral line., Comment: 7 pages, 8 figures, regular paper