Search

We study local magnetohydrodynamical instabilities of differential rotation in magnetized, stably stratified regions of stars and planets using a Cartesian Boussinesq model. We consider arbitrary latitudes and general shears (with gravity direction misaligned from this by an angle |$\phi$|⁠), to model radial (⁠|$\phi =0$|⁠), latitudinal (⁠|$\phi =\pm 90^\circ$|⁠), and mixed differential rotations, and study both non-diffusive [including magnetorotational instability (MRI) and Solberg–Høiland instability] and diffusive instabilities [including Goldreich–Schubert–Fricke (GSF) and MRI with diffusion]. These instabilities could drive turbulent transport and mixing in radiative regions, including the solar tachocline and the cores of red giant stars, but their dynamics are incompletely understood. We revisit linear axisymmetric instabilities with and without diffusion and analyse their properties in the presence of magnetic fields, including deriving stability criteria and computing growth rates, wave vectors, and energetics, both analytically and numerically. We present a more comprehensive analysis of axisymmetric local instabilities than prior work, exploring arbitrary differential rotations and diffusive processes. The presence of a magnetic field leads to stability criteria depending upon angular velocity rather than angular momentum gradients. We find MRI operates for much weaker differential rotations than the hydrodynamic GSF instability, and that it typically prefers much larger length-scales, while the GSF instability is impeded by realistic strength magnetic fields. We anticipate MRI to be more important for turbulent transport in the solar tachocline than the GSF instability when |$\phi \gt 0$| in the Northern (and vice versa in the Southern) hemisphere, though the latter could operate just below the convection zone when MRI is absent for |$\phi \lt 0$|⁠. [ABSTRACT FROM AUTHOR]

Stochastic and secular variations in the spin frequency |$\nu$| of a rotation-powered pulsar complicate the interpretation of the measured braking index, n , in terms of a power-law spin-down torque |$\propto \nu ^{n_{\rm pl}}$|⁠. Both categories of variation can lead to anomalous braking indices, with |$\vert n \vert = \vert \nu \ddot{\nu } / \dot{\nu }^2 \vert \gg 1$|⁠ , where the overdot symbolizes a derivative with respect to time. Here, we quantify the combined effect of stochastic and secular deviations from pure power-law spin-down on measurements of n. Through analytical calculations, Monte Carlo simulations involving synthetic data, and modern Bayesian timing techniques, it is shown that the variance of n satisfies the predictive, falsifiable formula |$\langle n^{2} \rangle = (n_{\rm pl}+\dot{K}_{\rm dim})^{2}+\sigma _{\rm dim}^{2}$|⁠ , where |$\dot{K}_{\rm dim}$| is inversely proportional to the time-scale |$\tau _K$| over which the proportionality constant of the power-law spin-down torque varies, |$\sigma _{\rm dim}$| is proportional to the timing noise amplitude and inversely proportional to the square root of the total observing time, and the average is over an ensemble of random realizations of the timing noise process. The anomalous regime |$\langle n^2 \rangle \gg 1$| occurs for |$\dot{K}_{\rm dim} \gg 1$|⁠ , |$\sigma _{\rm dim} \gg 1$|⁠ , or both. The sign of n depends in part on the sign of |$\dot{K}_{\rm dim}$|⁠ , so it is possible to measure unequal numbers of positive and negative n values in a large sample of pulsars. The distinguishable impact of stochastic and secular anomalies on phase residuals is quantified to prepare for extending the analysis of synthetic data to real pulsars. [ABSTRACT FROM AUTHOR]

The mechanism of X-ray outbursts in Be X-ray binaries remains a mystery, and understanding their circumstellar discs is crucial for a solution of the mass-transfer problem. In particular, it is important to identify the Be star activities (e.g. pulsations) that cause mass ejection and, hence, disc formation. Therefore, we investigated the relationship between optical flux oscillations and the infrared (IR) excess in a sample of five Be X-ray binaries. Applying the Lomb–Scargle technique to high-cadence optical light curves from the Transiting Exoplanet Survey Satellite (⁠|${\it TESS}$|⁠), we detected several significant oscillation modes in the 3–24 h period range for each source. We also measured the IR excess (a proxy for disc growth) of those five sources, using J -band light curves from Palomar Gattini-IR. In four of the five sources, we found anticorrelations between the IR excess and the amplitude of the main flux oscillation modes. This result is inconsistent with the conventional idea that non-radial pulsations drive mass ejections. We propose an alternative scenario where internal temperature variations in the Be star cause transitions between pulsation-active and mass-ejection-active states. [ABSTRACT FROM AUTHOR]

Red giant stars host solar-like oscillations which have mixed character, being sensitive to conditions both in the outer convection zone and deep within the interior. The properties of these modes are sensitive to both core rotation and magnetic fields. While asteroseismic studies of the former have been done on a large scale, studies of the latter are currently limited to tens of stars. We aim to produce the first large catalogue of both magnetic and rotational perturbations. We jointly constrain these parameters by devising an automated method for fitting the power spectra directly. We successfully apply the method to 302 low-luminosity red giants. We find a clear bimodality in core rotation rate. The primary peak is at |$\delta \nu _{\mathrm{rot}}$| = 0.32  |$\mu$| Hz, and the secondary at |$\delta \nu _{\mathrm{rot}}$| = 0.47  |$\mu$| Hz. Combining our results with literature values, we find that the percentage of stars rotating much more rapidly than the population average increases with evolutionary state. We measure magnetic splittings of 2 |$\sigma$| significance in 23 stars. While the most extreme magnetic splitting values appear in stars with masses |$\gt $| 1.1 M |$_{\odot }$|⁠ , implying they formerly hosted a convective core, a small but statistically significant magnetic splitting is measured at lower masses. Asymmetry between the frequencies of a rotationally split multiplet has previously been used to diagnose the presence of a magnetic perturbation. We find that of the stars with a significant detection of magnetic perturbation, 43 per cent do not show strong asymmetry. We find no strong evidence of correlation between the rotation and magnetic parameters. [ABSTRACT FROM AUTHOR]

The viscous decretion disc (VDD) scenario has been used to model the observables of the Be stars. Its capability to predict individual observables has been confirmed for several Be stars. Here, we simultaneously analyse the spectroscopic and |$BVI$| -band photometric data for the Be star MT91-213 with the Monte Carlo radiative-transfer code hdust to determine the stellar parameters, geometry, and physical conditions for its circumstellar disc. MT91-213 is the primary component of a binary system whose companion is the pulsar PSR J2032+4127. We find that the VDD model can simultaneously reproduce the multiple observables qualitatively, but not quantitatively. We determine the mass of the primary star to be 13.1  |${\rm M}_{\odot }$| which is smaller than reported in the literature. We present a dynamical scenario for the evolving disc density from a diffuse to a dense phase. Also, we determine that the inclination of the disc is about 40 |$^{\circ }$| which means it is 20 |$^{\circ }$| tilted from the orbit of the secondary star. Our results indicate that the mass loss rate for MT91-213 is |$\sim 10^{-7}$| to |$10^{-6} {\rm M}_{\odot }\,\mathrm{yr}^{-1}$| which is in agreement with the suggested values in the literature, required to explain the observed X-ray synchrotron luminosity, |$L_\mathrm{x}$|⁠ , for PSR J2032+4127. [ABSTRACT FROM AUTHOR]

At metallicities lower than that of the Small Magellanic Cloud, it remains essentially unexplored how fossil magnetic fields, forming large-scale magnetospheres, could affect the evolution of massive stars, thereby impacting the fundamental building blocks of the early Universe. We extend our stellar evolution model grid with representative calculations of main-sequence, single-star models with initial masses of 20 and 60 M |$_\odot$|⁠ , including appropriate changes for low-metallicity environments (⁠|$Z = 10^{-3}$| – |$10^{-6}$|⁠). We scrutinize the magnetic, rotational, and chemical properties of the models. When lowering the metallicity, the rotational velocities can become higher and the tendency towards quasi-chemically homogeneous evolution increases. While magnetic fields aim to prevent the development of this evolutionary channel, the weakening stellar winds lead to less efficient magnetic braking in our models. Since the stellar radius is almost constant during a blueward evolution caused by efficient chemical mixing, the surface magnetic field strength remains unchanged in some models. We find core masses at the terminal-age main sequence between 22 and 52 M |$_\odot$| for initially 60 M |$_\odot$| models. This large difference is due to the vastly different chemical and rotational evolution. We conclude that in order to explain chemical species and, in particular, high nitrogen abundances in the early Universe, the adopted stellar models need to be under scrutiny. The assumptions regarding wind physics, chemical mixing, and magnetic fields will strongly impact the model predictions. [ABSTRACT FROM AUTHOR]

Context. The evolution and fate of massive stars are thought to be affected by rotationally induced internal mixing. The surface boron abundance is a sensitive tracer of this in early B-type main sequence stars. Aims. We test current stellar evolution models of massive main sequence stars which include rotational mixing through a systematic study of their predicted surface boron depletion. Methods. We construct a dense grid of rotating single star models using MESA, for which we employ a new nuclear network which follows all the stable isotopes up to silicon, including lithium, beryllium, boron, as well as the radioactive isotope aluminium-26. We also compile the measured physical parameters of the 90 Galactic early B-type stars with boron abundance information. We then compare each observed stars with our models through a Bayesian analysis, which yields the mixing efficiency parameter with which the star is reproduced the best, and the probability that it is represented by the stellar models. Results. We find that about two-thirds of the sample stars are well represented by the stellar models, with the best agreement achieved for a rotational mixing efficiency of ∼50% compared to the widely adopted value. The remaining one third of the stars, of which many are strongly boron depleted slow rotators, are largely incompatible with our models, for any rotational mixing efficiency. We investigate the observational incidence of binary companions and surface magnetic fields, and discuss their evolutionary implications. Conclusions. Our results confirm the concept of rotational mixing in radiative stellar envelopes. On the other hand, we find that a different boron depletion mechanism, and likely a different formation path, is required to explain about one-third of the sample stars. The large spread in the surface boron abundances of these stars may hold a clue to understanding their origin. [ABSTRACT FROM AUTHOR]

Context. In the partially ionized material of stellar interiors, the strongest forces acting on electrons and ions are the Coulomb interactions between charges. The dynamics of the plasma as a whole depend on the magnitudes of the average electrostatic interactions and the average kinetic energies of the particles that constitute the stellar material. An important question is how these interactions of real gases are related to the observable stellar properties. Specifically, the relationships between rotation, magnetic activity, and the thermodynamic properties of stellar interiors are still not well understood. These connections are crucial for understanding and interpreting the abundant observational data provided by space-based missions, such as Kepler/K2 and TESS, and the future data from the PLATO mission. Aims. In this study, we investigate the electrostatic effects within the interiors of low-mass main sequence (MS) stars. Specifically, we introduce a global quantity, a global plasma parameter, which allows us to compare the importance of electrostatic interactions across a range of low-mass theoretical models (0.7 − 1.4 M⊙) with varying ages and metallicities. We then correlate the electrostatic properties of the theoretical models with the observable rotational trends on the MS. Methods. We use the open-source 1D stellar evolution code MESA to compute a grid of main-sequence stellar models. Our models span the log g − Teff space of a set of 66 Kepler main-sequence stars. Results. We identify a correlation between the prominence of electrostatic effects in stellar interiors and stellar rotation rates. The variations in the magnitude of electrostatic interactions with age and metallicity further suggest that understanding the underlying physics of the collective effects of plasma can clarify key observational trends related to the rotation of low-mass stars on the MS. These results may also advance our understanding of the physics behind the observed weakened magnetic braking in stars. [ABSTRACT FROM AUTHOR]

Context. Many fast-rotating stars (rotation periods of < 2 days) are found to be unresolved binaries with separations of tens of AU. This correlation between fast rotators and binarity leads to the question of whether the formation of binary stars inherently produces fast rotators. Aims. Our goal is to understand the spin evolution of protostars and whether the formation of companions plays a role in spinning up stars. Methods. We used magneto-hydrodynamic simulations to study the formation of multiple star systems from turbulent and non-turbulent protostellar cores. We tracked the angular momentum accreted by individual star and inner disc systems by using a sink (star) particle technique. We ran a resolution study to extrapolate protostellar properties. Results. We find in all simulations that the primary star can experience a spin-up event correlated with the formation of companions, namely fragmentation into binaries or higher-order systems. The primary star can spin up by up to 84% of its pre-fragmentation angular momentum and by up to 18% of its pre-fragmentation mass-specific angular momentum. The mechanism for the spin-up is gravitational disc instabilities in the circumstellar disc around the primary star, which leads to the accretion of material with high specific angular momentum. The simulations that experience the strongest disc instabilities fragment to form companions. Simulations with weaker spin-up events experience disc instabilities triggered by a companion flyby, and the disc instability in these cases typically does not produce further fragments (i.e. they remain binary systems). Conclusions. The primary star in multiple star systems can end up with a higher spin than single stars. This is because gravitational instabilities in the circumstellar disc around the primary star can trigger a spin-up event. In the strongest spin-up events, the instability is likely to cause disc fragmentation and the formation of companions. This spin-up mechanism, coupled with shorter disc lifetimes due to truncated circumstellar discs (and thus short spin-down times), may help produce fast rotators. [ABSTRACT FROM AUTHOR]

Context. Grids of stellar evolution models with rotation using the Geneva stellar evolution code (GENEC) have been published for a wide range of metallicities. Aims. We introduce the last remaining grid of GENEC models, with a metallicity of Z = 10−5. We study the impact of this extremely metal-poor initial composition on various aspects of stellar evolution, and compare it to the results from previous grids at other metallicities. We provide electronic tables that can be used to interpolate between stellar evolution tracks and for population synthesis. Methods. Using the same physics as in the previous papers of this series, we computed a grid of stellar evolution models with GENEC spanning masses between 1.7 and 500 M⊙, with and without rotation, at a metallicity of Z = 10−5. Results. Due to the extremely low metallicity of the models, mass-loss processes are negligible for all except the most massive stars. For most properties (such as evolutionary tracks in the Hertzsprung-Russell diagram, lifetimes, and final fates), the present models fit neatly between those previously computed at surrounding metallicities. However, specific to this metallicity is the very large production of primary nitrogen in moderately rotating stars, which is linked to the interplay between the hydrogen- and helium-burning regions. Conclusions. The stars in the present grid are interesting candidates as sources of nitrogen-enrichment in the early Universe. Indeed, they may have formed very early on from material previously enriched by the massive short-lived Population III stars, and as such constitute a very important piece in the puzzle that is the history of the Universe. [ABSTRACT FROM AUTHOR]

Context. The Transiting Exoplanet Survey Satellite (TESS) performs high-precision photometry over almost the whole sky primarily in search of exoplanet transits. It also provides exquisite data to study stellar variability, in particular for pulsating hot B subdwarf (sdB) stars. Aims. We present a detailed analysis of a new hybrid (p- and g-mode) sdB pulsator, TIC 441725813 (TYC 4427-1021-1), discovered and monitored by TESS for 670 days. Methods. The TESS light curves available for this star were analysed using prewhitening techniques to extract mode frequencies accurately. The pulsation spectrum was then interpreted through methods that include asymptotic period spacing relationships and the identification of rotational multiplets. We also exploited a high signal-to-noise ratio (S/N), low-resolution spectrum of TIC 441725813 using grids of non-local thermodynamic equilibrium (NLTE) model atmospheres to derive its atmospheric parameters. Results. The light curve analysis reveals that frequencies are mostly found in the g-mode region, but several p-modes are also detected, indicating that TIC 441725813 is a hybrid sdB pulsator. We identify 25 frequencies that can be associated with ℓ = 1 g-modes, 15 frequencies corresponding to ℓ = 2 g-modes, and six frequencies characteristic of p-modes. Interestingly, several frequency multiplets interpreted as rotational splittings of deep-probing g-modes indicate a slow rotation period of at least 85.3 ± 3.6 day, while splittings of mostly envelope-probing p-modes suggest a significantly shorter rotation period of 17.9 ± 0.7 day, which implies the core (mainly the helium mantle with possibly the deeper partially mixed helium-burning core that it surrounds) rotates at least 4.7 times slower than the envelope. The radial velocity curves indicate that TIC 441725813 is in a close binary system with a low-luminosity companion, possibly a white dwarf. While elusive in the available TESS photometry, a low-frequency signal that would correspond to a period of ∼6.7 h is found, albeit at a low S/N. Furthermore, we estimate the inclination angle to be ∼60° by two independent means. Conclusions. TIC 441725813 is a particularly interesting sdB star whose envelope rotates faster than the core. We hypothesise that this might be caused by the effects of a tidal interaction with a companion, although in the present case, the presence of such a companion will have to be further investigated. This analysis paves the way towards a more detailed seismic probing of TIC 441725813 using optimisation techniques, which will be presented in a second paper. [ABSTRACT FROM AUTHOR]

Context. Constraining the coplanarity status of very low mass (VLM) tight binary systems provides valuable information towards understanding the dominant mechanisms in star and planetary formation at the lower end of the initial mass function. Aims. We sought to constrain the v sin i degeneracies of two nearby tight VLM binary systems, 2MASS J0746+20AB and 2MASS J1314+13AB, by independently determining each companion's rotational period and so characterize each system's spin-orbit alignment. Methods. Long observational baseline high cadence I band photometry data were obtained for both systems using the Galway Ultra Fast Imager (GUFI) on the Mount Graham International Observatory's Vatican Advanced Technology Telescope. Previously known rotational periods determined in other passbands were used as a basis for recovering additional periodic modulations in the time-series data. Results. Using the known rotational period of 3.32 hours for 2MASS J0746425+200032A, we recovered an underlying periodic modulation of 2.14 ± 0.11 hours, which we associate with 2MASS J0746425+200032B. Breaking each components' v sin i corresponds to equatorial inclination angles of 32 ± 4 degrees and 37 ± 4 degrees for components A & B respectively. We recover a weaker 2.06 ± 0.05 hours modulation separate from the known 3.79 hour signature for J1314203+132001B, which we associate with J1314203+132001A. We place a lower limit on J1314203+132001A's equatorial inclination angle to be in the range of 24.5−3+3.5 degrees, which deviates from the system's orbital plane previously determined to be 49.34−0.23+0.28 degrees. Conclusions. We confirm long term, consistent periodic modulations from both binary systems and report the first definitive rotational period for J1314203+132001A of 2.06 ± 0.05 hours, in addition to coplanarity to within 10 degrees in the spin-orbit alignment of the 2MASS J0746425+200032AB system. The lack of separate v sin i values for the 2MASS J1314203+132001AB system limits a definitive assessment but best estimates suggest coplanarity to be unlikely. [ABSTRACT FROM AUTHOR]

Context. Tidal deformation breaks the axisymmetric structure of stars, and this may affect stellar activity. This effect has been demonstrated in theoretical analyses and simulations, but it lacks observational support. In this paper, we use spot-modulated detached binaries to study the effect of tides on stellar activity. We show this effect by analyzing the properties of the spot lifetime, the harmonic decay timescale, and the orbital parameters. Aims. We aim to explore the differences in spot lifetimes between binaries and single stars, the main mechanisms of spot decay in binaries, and the correlation between orbital parameters and spot lifetimes. These differences will provide clues to the effect of tides on stellar activity. Methods. We collected data of 311 spot-modulated detached binaries and 3272 single stars. The relative orbital parameters of the binaries were derived by combining Kepler photometry, stellar atmospheric parameters from LAMOST DR9 and Gaia DR3, and 2MASS photometry. We then used the ACF method to obtain the rotational periods, lifetimes, and harmonic decay timescales. Finally, we analyzed the correlation between the lifetime of spots and orbital parameters, explored the dominant decay mechanism of spots, and examined the differences in spots for binaries and single stars. Results. The relative lifetime of a starspot is correlated with the sum of the fractional radii, the orbital eccentricity, and the synchronization ratio. Longer lifetimes are observed in close, circular, and synchronous binaries than in the other binaries. The main mechanism for the decay of star spots in binaries is large-scale convective motion. However, on close, cool, and fast-rotating binaries, horizontal diffusion or subphotospheric diffusion are dominant. Compared to single stars, the median lifetime of a starspot on binaries was found to be longer. Moreover, this difference decreases with rotation period. Additionally, it should be noted that spots on binaries experience increased horizontal or subphotospheric diffusion at the same rotation period and effective temperature. Conclusions. According to the observation results, we conclude that the lifetime of starspots on detached close binaries is affected by tidal interactions. [ABSTRACT FROM AUTHOR]

Context. The estimation of rotation periods of stars is a key challenge in stellar astrophysics. Given the large amount of data available from ground-based and space-based telescopes, there is a growing interest in finding reliable methods to quickly and automatically estimate stellar rotation periods with a high level of accuracy and precision. Aims. This work aims to develop a computationally inexpensive approach, based on machine learning techniques, to accurately predict thousands of stellar rotation periods. Methods. The innovation in our approach is the use of the XGBoost algorithm to predict the rotation periods of Kepler targets by means of regression analysis. Therefore, we focused on building a robust supervised machine learning model to predict surface stellar rotation periods from structured data sets built from the Kepler catalogue of K and M stars. We analysed the set of independent variables extracted from Kepler light curves and investigated the relationships between them and the ground truth. Results. Using the extreme gradient boosting (GB) method, we obtained a minimal set of variables that can be used to build machine learning models for predicting stellar rotation periods. Our models have been validated by predicting the rotation periods of about 2900 stars. The results are compatible with those obtained by classical techniques and comparable to those obtained by other recent machine learning approaches, with the advantage of using fewer predictors. When restricting the analysis to stars with rotation periods of less than 45 d, our models are on average wrong less than 5% of the time. Conclusions. We have developed an innovative approach based on a machine learning method to accurately fit the rotation periods of stars. Based on the results of this study, we conclude that the best models generated by the proposed methodology can compete with the latest state-of-the-art approaches, while offering the advantage of being computationally cheaper, easy to train, and reliant only on small sets of predictors. [ABSTRACT FROM AUTHOR]

Context. Grids of stellar evolution models with rotation using the Geneva stellar evolution code (GENEC) have been published for a wide range of metallicities. Aims. We introduce the last remaining grid of GENEC models, with a metallicity of Z = 10−5. We study the impact of this extremely metal-poor initial composition on various aspects of stellar evolution, and compare it to the results from previous grids at other metallicities. We provide electronic tables that can be used to interpolate between stellar evolution tracks and for population synthesis. Methods. Using the same physics as in the previous papers of this series, we computed a grid of stellar evolution models with GENEC spanning masses between 1.7 and 500 M⊙, with and without rotation, at a metallicity of Z = 10−5. Results. Due to the extremely low metallicity of the models, mass-loss processes are negligible for all except the most massive stars. For most properties (such as evolutionary tracks in the Hertzsprung-Russell diagram, lifetimes, and final fates), the present models fit neatly between those previously computed at surrounding metallicities. However, specific to this metallicity is the very large production of primary nitrogen in moderately rotating stars, which is linked to the interplay between the hydrogen- and helium-burning regions. Conclusions. The stars in the present grid are interesting candidates as sources of nitrogen-enrichment in the early Universe. Indeed, they may have formed very early on from material previously enriched by the massive short-lived Population III stars, and as such constitute a very important piece in the puzzle that is the history of the Universe. [ABSTRACT FROM AUTHOR]

With the release of a large amount of astronomical data, an increasing number of close-in hot Jupiters have been discovered. Calculating their evolutionary curves using star–planet interaction models presents a challenge. To expedite the generation of evolutionary curves for these close-in hot Jupiter systems, we utilized tidal interaction models established on mesa to create 15 745 samples of star–planet systems and 7500 samples of stars. Additionally, we employed a neural network (Multilayer Perceptron – MLP) to predict the evolutionary curves of the systems, including stellar effective temperature, radius, stellar rotation period, and planetary orbital period. The median relative errors of the predicted evolutionary curves were found to be 0.15 per cent, 0.43 per cent, 2.61 per cent, and 0.57 per cent, respectively. Furthermore, the speed at which we generate evolutionary curves exceeds that of model-generated curves by more than four orders of magnitude. We also extracted features of planetary migration states and utilized lightgbm to classify the samples into six categories for prediction. We found that by combining three types that undergo long-term double synchronization into one label, the classifier effectively recognized these features. Apart from systems experiencing long-term double synchronization, the median relative errors of the predicted evolutionary curves were all below 4 per cent. Our work provides an efficient method to save significant computational resources and time with minimal loss in accuracy. This research also lays the foundation for analysing the evolutionary characteristics of systems under different migration states, aiding in the understanding of the underlying physical mechanisms of such systems. Finally, to a large extent, our approach could replace the calculations of theoretical models. [ABSTRACT FROM AUTHOR]

We present a software package designed to produce photometric light curves and measure rotation periods from full-frame images taken by the Transiting Exoplanet Survey Satellite (TESS), which we name ' tessilator '. tessilator is the only publicly available code that will run a full light curve and rotation period (⁠|$P_{\rm rot}$|⁠) analysis based on just a (list of) target identifier(s) or sky position(s) via a simple command-line prompt. This paper sets out to introduce the rationale for developing tessilator , and then describes the methods, considerations, and assumptions for: extracting photometry; dealing with potential contamination; accounting for natural and instrumental systematic effects; light-curve normalization and detrending; removing outliers and unreliable data; and finally, measuring the |$P_{\rm rot}$| value and several periodogram attributes. Our methods have been tuned specifically to optimize TESS light curves and are independent from the pipelines developed by the TESS Science Processing Operations Center (SPOC), meaning tessilator can, in principle, analyse any target across the entire celestial sphere. We compare tessilator |$P_{\rm rot}$| measurements with TESS -SPOC-derived light curves of 1560 (mainly solar-type and low-mass) stars across four benchmark open clusters (Pisces–Eridanus, the Pleiades, the Hyades, and Praesepe) and a sample of nearby field M-dwarfs. From a vetted subsample of 864 targets we find an excellent return of |$P_{\rm rot}$| matches for the first three open clusters (⁠|$\gt 85$|  per cent) and a moderate (⁠|$\sim 60$|  per cent) match for the 700 Myr Praesepe and MEarth sample, which validates tessilator as a tool for measuring |$P_{\rm rot}$|⁠. The tessilator code is available at https://github.com/alexbinks/tessilator. [ABSTRACT FROM AUTHOR]

In recent years, there has been a push to understand how chemical composition affects the magnetic activity levels of main sequence low-mass stars. Results indicate that more metal-rich stars are more magnetically active for a given stellar mass and rotation period. This metallicity dependence has implications for how the rotation periods and activity levels of low-mass stars evolve over their lifetimes. Numerical modelling suggests that at late ages more metal-rich stars should be rotating more slowly and be more magnetically active. In this work, we study the rotation and activity evolution of low-mass stars using a sample of Kepler field stars. We use the gyro-kinematic age dating technique to estimate ages for our sample and use the photometric activity index as our proxy for magnetic activity. We find clear evidence that, at late ages, more metal-rich stars have spun down to slower rotation in agreement with the theoretical modelling. However, further investigation is required to definitively determine whether the magnetic activity evolution occurs in a metallicity dependent way. [ABSTRACT FROM AUTHOR]

We extend our model of magnetic braking (MB), driven by an α–Ω dynamo mechanism, from fully convective M-dwarfs (FCMDs) to explain the surface and internal spin |$P_\mathrm{spin}$| evolution of partly convective dwarfs (PCDs) starting from the disc-dispersal stage to the main-sequence turn-off. In our model, the spin of the core is governed by shear at the core-envelope boundary while the spin of the envelope is governed by MB and shear. We show that (1) the most massive FCMDs experience a stronger spin-down than PCDs and less massive FCMDs, (2) the stalled spin-down and enhanced activity of K-dwarfs and the pile-up of G-dwarfs older than a few Gyr are stellar-structure and MB-dependent, and weakly dependent on core-envelope coupling effects, (3) our expression of the core-envelope convergence time-scale |$\tau _{\text{converge}}(M_\ast , P_{\text{spin}})$| between a few 10 to 100 Myr strongly depends on stellar structure but weakly on MB strength and shear, such that fast and massive rotators achieve corotation earlier, (4) our estimates of the surface magnetic fields are in general agreement with observations and our wind mass loss evolution explains the weak winds from the solar analog |$\pi ^1$| UMa, and (5) with our model the massive young Sun hypothesis as a solution to the faint young Sun problem can likely be ruled out, because the maximum mass lost by winds from our Sun with our model is about an order of magnitude smaller than required to solve the problem. [ABSTRACT FROM AUTHOR]

When a neutron star is spun-up or spun-down, the changing strains in its solid elastic crust can give rise to sudden fractures known as starquakes. Early interest in starquakes focused on their possible connection to pulsar glitches. While modern glitch models rely on pinned superfluid vorticity rather than crustal fracture, starquakes may nevertheless play a role in the glitch mechanism. Recently, there has been interest in the issue of starquakes resulting in non-axisymmetric shape changes, potentially linking the quake phenomenon to the building of neutron star mountains, which would then produce continuous gravitational waves. Motivated by this issue, we present a simple model that extends the energy minimization-based calculations, originally developed to model axisymmetric glitches, to also include non-axisymmetric shape changes. We show that the creation of a mountain in a quake necessarily requires a change in the axisymmetric shape too. We apply our model to the specific problem of the spin-up of an initially non-rotating star, and estimate the maximum mountain that can be built in such a process, subject only to the constraints of energy and angular momentum conservation. [ABSTRACT FROM AUTHOR]

The origin of extended main-sequence turn-offs (eMSTO) in star clusters younger than 2 Gyr still challenges our current understanding of stellar evolution. Exploiting data from Gaia Data Release 3 (DR3), we investigate eMSTOs in a large sample of 32 Galactic open clusters younger than 2.4 Gyr. We first validate Gaia rotational velocities from radial velocity spectrometer spectra by comparing them with literature values and assessing their correlation with magnetic activity measurements from LAMOST spectra. We detect a general positive correlation between turn-off colour and projected stellar rotation, with slow-rotating stars predominantly found on the bluer side of the turn-off. Comparing our observations with theoretical models, we find that the eMSTO morphology is well reproduced by a single population formed with a high rotation rate, and observed with rotation axis inclination ranging between 0° (pole-on) and 90° (edge-on). This contrasts with observations of Magellanic Cloud clusters, where a population of non-rotating stars appears to be ubiquitous in clusters younger than 700 Myr. However, we note that our interpretation, while successfully explaining the overall eMSTO morphology, cannot fully explain the observed projected rotational velocities. Additionally, two young clusters, NGC 3532 and NGC 2287, exhibit moderate evidence of a split-main sequence in colour and rotation, suggesting a possible small spread in the initial rotation rate. Finally, we advise caution in determining the ages of young clusters from non-rotating isochrones, as neglecting the effects of stellar rotation can impact the isochrone dating by up to factors of 5–20 per cent. [ABSTRACT FROM AUTHOR]

The presence of extended Main-Sequence Turn-Off (eMSTO) in the open clusters has been attributed to various factors, such as spread in rotation rates, binary stars, and dust-like extinction from stellar excretion discs. We present a comprehensive analysis of the eMSTO in the open cluster NGC 2355. Using spectra from the Gaia –European Southern Observatory (ESO) archives, we find that the stars in the red part of the eMSTO have a higher mean v sin i value of 135.3 |$\pm$| 4.6 km s |$^{-1}$| compared to the stars in the blue part that have an average v sin i equal to 81.3 |$\pm$| 5.6 km s |$^{-1}$|⁠. This suggests that the eMSTO in NGC 2355 is possibly caused by the spread in rotation rates of stars. We do not find any substantial evidence of the dust-like extinction from the eMSTO stars using ultraviolet data from the Swift survey. The estimated synchronization time for low-mass ratio close binaries in the blue part of the eMSTO suggests that they would be mostly slow-rotating if present. However, the stars in the blue part of the eMSTO are preferentially located in the outer region of the cluster indicating that they may lack low-mass ratio close binaries. The spread in rotation rates of eMSTO stars in NGC 2355 is most likely caused by the star-disc interaction mechanism. The stars in the lower main sequence beyond the eMSTO region of NGC 2355 are slow rotating (mean v sin i = 26.5 |$\pm$| 1.3 km s |$^{-1}$|⁠) possibly due to the magnetic braking of their rotations. [ABSTRACT FROM AUTHOR]

General-relativistic equilibria of differentially rotating stars are expected in a number of astrophysical scenarios, from core-collapse supernovae to the remnant of binary neutron-star mergers. The latter, in particular, have been the subject of extensive studies where they were modelled with a variety of laws of differential rotation with varying degree of realism. Starting from accurate and fully general-relativistic simulations of binary neutron-star mergers with various equations of state and mass ratios, we establish the time when the merger remnant has reached a quasi-stationary equilibrium and extract in this way realistic profiles of differential rotation. This allows us to explore how well traditional laws reproduce such differential-rotation properties and to derive new laws of differential rotation that better match the numerical data in the low-density Keplerian regions of the remnant. In this way, we have obtained a novel and somewhat surprising result: the dynamical stability line to quasi-radial oscillations computed from the turning-point criterion can have a slope that is not necessarily negative with respect to the central rest-mass density, as previously found with traditional differential-rotation laws. Indeed, for stellar models reproducing well the properties of the merger remnants, the slope is actually positive, thus reflecting remnants with angular momentum at large distances from the rotation axis, and hence with cores having higher central rest-mass densities and slower rotation rates. [ABSTRACT FROM AUTHOR]

We present the comprehensive analysis of the high-amplitude |$\delta$| Sct star V2367 Cygni. First, we perform the frequency analysis for the whole available set of the Kepler and TESS photometry. Most of the frequency peaks are harmonics or combinations of the three known independent frequencies with the highest amplitudes, i.e. |$\nu _1=5.661\,06$|  d |$^{-1}$|⁠ , |$\nu _2=7.14898$|  d |$^{-1}$|⁠ , and |$\nu _3=7.77557$|  d |$^{-1}$|⁠. The total number of independent frequencies is 26 and 25 from the Kepler and TESS light curve, respectively. Then, using the |${\it UBVRI}$| time-series photometry, we unambiguously identify the dominant frequency |$\nu _1$| as the radial mode, whereas in the case of frequencies |$\nu _2$| and |$\nu _3$| the most probable mode degrees are |$\ell =0$| or |$\ell =2$|⁠. However, only the frequency |$\nu _2$| can be associated with a radial mode, and only if higher order effects of rotation are taken into account. Including the rotational mode coupling, we constructed complex seismic models of V2367 Cyg, which fit |$\nu _1$| and |$\nu _2$| as radial modes, and reproduce the amplitude of bolometric flux variations (the parameter f) for the dominant mode. The empirical values of f are derived from the |${\it UBVRI}$| amplitudes and phases. We rely on the Bayesian analysis based on Monte Carlo simulations to derive constraints on evolutionary stage, mass, rotation, overshooting from the convective core, and efficiency of convective transport in the envelope. Our seismic analysis clearly indicates that V2367 Cyg is in a post-main sequence phase of evolution. This is the first extensive seismic modelling that takes into account the effect of rotational coupling between pulsation modes. [ABSTRACT FROM AUTHOR]

In our Galactic Centre, about |$10\,000$| to |$100\,000$| stars are estimated to have survived tidal disruption events, resulting in partially disrupted remnants. These events occur when a supermassive black hole (SMBH) tidally interacts with a star, but not enough to completely disrupt the star. We use the 1D stellar evolution code Kepler and the 3D smoothed particle hydrodynamics code Phantom to model the tidal disruption of 1, 3, and |$10\, \mathrm{M}_\odot$| stars at zero-age main sequence (ZAMS), middle-age main sequence (MAMS), and terminal-age main sequence (TAMS). We map the disruption remnants into Kepler in order to understand their post-distribution evolution. We find distinct characteristics in the remnants, including increased radius, rapid core rotation, and differential rotation in the envelope. The remnants undergo composition mixing that affects their stellar evolution. Although the remnants formed by disruption of ZAMS models evolve similarly to unperturbed models of the same mass, for MAMS and TAMS stars, the remnants have higher luminosity and effective temperature. Potential observational signatures include peculiarities in nitrogen and carbon abundances, higher luminosity, rapid rotation, faster evolution, and unique tracks in the Hertzsprung–Russell diagram. [ABSTRACT FROM AUTHOR]

The properties of blue supergiants are key for constraining the end of the main sequence (MS) of massive stars. Whether the observed drop in the relative number of fast-rotating stars below ≈21 kK is due to enhanced mass-loss rates at the location of the bistability jump, or the result of the end of the MS is still debated. Here, we combine newly derived estimates of photospheric and wind parameters with Gaia distances and wind terminal velocities from the literature to obtain upper limits on the mass-loss rates for a sample of 116 Galactic luminous blue supergiants. The parameter space covered by the sample ranges between 35–15 kK in Teff and 4.8–5.8 dex in log(L/L⊙). Our results show no increase in the mass-loss rates over the bistability jump. Therefore, we argue that the drop in rotational velocities cannot be explained by enhanced mass loss. Since a large jump in the mass-loss rates is commonly included in evolutionary models, we suggest an urgent revision of the default prescriptions currently in use. [ABSTRACT FROM AUTHOR]

Context. Lithium (Li) abundance is an age indicator for G, K, and M stellar types, as its abundance decreases over time for these spectral types. However, despite all of the observational efforts made over the past few decades, the role of rotation, stellar activity, and metallicity in the depletion of Li is still unclear. Aims. Our purpose is to investigate how Li depletion is affected by rotation and metallicity in G and K members of the roughly Pleiades-aged open cluster M35. Methods. We have collected an initial sample of 165 candidate members observed with the WIYN/Hydra spectrograph. In addition, we have taken advantage of three previous spectroscopic studies of Li in M35. As a result, we have collected a final sample of 396 stars observed with the same instrument, which we have classified as non-members, possible non-members, possible members, and probable members of the cluster. We have measured iron abundances, Li equivalent widths, and Li abundances for the 110 M35 members added to the existing sample by this study. Finally, rotation periods for cluster members have been obtained from the literature or derived from Zwicky Transient Facility light curves. Results. We have confirmed that fast G and K rotators are Li-rich in comparison with slow rotators of similar effective temperature. This trend, which is also seen in previous studies, is more evident when binaries are not taken into account. Furthermore, while we derived an average metallicity of [Fe/H] = −0.26 ± 0.09 from our spectra, the distribution of Li in M35 is similar to those observed for the Pleiades and M34 open clusters, which have solar metallicity and slightly different ages. In addition, we have shown that an empirical relationship proposed to remove the contribution of the Fe I line at 670.75 nm to the blended feature at 670.78 nm overestimates by 5–15 mÅ the contribution of this iron line for M35 members. Conclusions. M35 fast G and K rotators have depleted less Li than their slower counterparts. Furthermore, a 0.2−0.3 dex difference in metallicity appears to make little difference in the Li distributions of open clusters with ages between 100 and 250 Myr. [ABSTRACT FROM AUTHOR]

Studies of the rotation and activity of M type stars are essential in order to enhance our understanding of stellar dynamos and angular momentum evolution. Using the outstanding photometric capabilities of space telescopes, even those rotation signals with low amplitudes can now be investigated in previously unrivaled detail. By combining data of K2 and the TESS prime mission, the star spot activity of M dwarfs can be monitored over a timescale of half a decade. In the framework of our study on the rotation–activity relation for bright and nearby M dwarfs, we also aim to investigate their long-term activity. While K2 was observing fields distributed around the ecliptic plane, the TESS prime mission was oriented along a line of ecliptic longitude, with one camera centered on an ecliptic pole. Due to these different observing strategies, the overlap between K2 and the TESS prime mission is marginal. However, 45 stars from our sample were observed with both missions, and of these, two early M-type stars that fulfill our selection criteria, EPIC 202059229 and EPIC 245919787, were analyzed in more detail. We find that, for both stars, the rotation period did not change, while the rotational phase did change for EPIC 245919787 by ∼0.2. The amplitude of the spot-induced variability changed for both stars but more significantly for EPIC 245919787. By comparing the cumulative flare frequency distributions, we find that the flare activity for EPIC 202059229 is unchanged, while it slightly changes for EPIC 245919787 between the K2 and TESS epochs. Using a combination of light curves from K2 and TESS that span a baseline of up to 4.5 yr, we are able to measure significant differential rotation for EPIC 245919787. Furthermore, we show that combining missions like K2 and TESS is a promising method for detecting stellar activity cycles. [ABSTRACT FROM AUTHOR]

Context. Understanding the effects of rotation in stellar evolution is key to modelling early-type stars, half of which have equatorial velocities over 100 km s−1. The nearby star Altair is an example of such fast-rotating stars, and furthermore, it has the privilege of being modelled by a detailed 2D concordance model that reproduces most of its observables. Aims. The aim of this paper is to include new asteroseismic frequencies to improve our knowledge of Altair, especially its age. Methods. We processed images of Altair obtained during July 2022 by the Transiting Exoplanet Survey Satellite using the halo photometry technique to obtain its light curve over this observation period. Results. By analysing the light curve, we derived a set of 22 new frequencies in the oscillation spectrum of Altair and confirmed 12 previously known frequencies. Compared with model predictions, we could associate ten frequencies with ten axisymmetric modes. This identification is based on the modelled visibility of the modes. Moreover, nine of the modelled frequencies can be adjusted to simultaneously match their corresponding observed frequencies, once the core hydrogen mass fraction of the concordance model is set to Xcore/Xini ≃ 0.972, with Xini = 0.739. Using the combined results of a 1D MESA model computing the pre-main sequence and a 2D time-dependent ESTER model computing the main sequence, we find that this core hydrogen abundance sets the age of Altair to 88 ± 10 Myr, which is slightly younger than previous estimates. [ABSTRACT FROM AUTHOR]

Aims. We aim to measure the impact of our mass-loss recipe in the evolution of massive stars at the metallicity of the Small Magellanic Cloud (SMC). Methods. We used the Geneva-evolution code (GENEC) to run evolutionary tracks for stellar masses ranging from 20 to 85 M⊙ at SMC metallicity (ZSMC = 0.002). We upgraded the recipe for stellar winds by replacing Vink's formula with our self-consistent m-CAK prescription, which reduces the value of the mass-loss rate, Ṁ, by a factor of between two and six depending on the mass range. Results. The impact of our new [weaker] winds is wide, and it can be divided between direct and indirect impact. For the most massive models (60 and 85 M⊙) with Ṁ ≳ 2 × 10−7M⊙ yr−1, the impact is direct because lower mass loss make stars remove less envelope, and therefore they remain more massive and less chemically enriched at their surface at the end of their main sequence (MS) phase. For the less massive models (20 and 25 M⊙) with Ṁ ≲ 2 × 10−8M⊙ yr−1, the impact is indirect because lower mass loss means the stars keep high rotational velocities for a longer period of time, thus extending the H-core burning lifetime and subsequently reaching the end of the MS with higher surface enrichment. In either case, given that the conditions at the end of the H-core burning change, the stars will lose more mass during their He-core burning stages anyway. For the case of Mzams = 20–40 M⊙, our models predict stars will evolve through the Hertzsprung gap, from O-type supergiants to blue supergiants (BSGs), and finally red supergiants (RSGs), with larger mass fractions of helium compared to old evolution models. New models also sets the minimal initial mass required for a single star to become a Wolf-Rayet (WR) at metallicity Z = 0.002 at Mzams = 85 M⊙. Conclusions. These results reinforce the importance of upgrading mass-loss prescriptions in evolution models, in particular for the earlier stages of stellar lifetime, even for Z ≪ Z⊙. New values for Ṁ need to be complemented with upgrades in additional features such as convective-core overshooting and distribution of rotational velocities, besides more detailed spectroscopical observations from projects such as XShootU, in order to provide a robust framework for the study of massive stars at low-metallicity environments. [ABSTRACT FROM AUTHOR]

Context. Numerical computations of stellar oscillations for models representative of B-type stars predict fewer modes to be excited than observations reveal from modern space-based photometric data. One shortcoming of state-of-the-art evolution models of B-type stars that may cause a lack of excited modes is the absence of microscopic diffusion in most such models. Aims. We investigate whether the inclusion of microscopic diffusion in stellar models of B-type stars, notably radiative levitation experienced by isotopes, leads to extra mode driving by the opacity mechanism compared to the case of models that do not include microscopic diffusion. Methods. We consider the case of slowly to moderately rotating stars and use non-rotating equilibrium models, while we account for (uniform) rotation in the computations of the pulsation frequencies. We calculate 1D stellar structure and evolution models with and without microscopic diffusion and examine the effect of radiative levitation on mode excitation, for both low-radial order pressure and gravity modes and for high-radial order gravity modes. As is common practice in asteroseismology, rotation is included in the pulsation computations according to the mode's frequency regime. For modes having frequencies below twice the rotation frequency, that is, modes in the sub-inertial regime, we adopt the traditional approximation of rotation. For modes in the super-inertial regime with frequency above twice the rotation frequency, rotation is treated perturbatively up to first order in the rotation. We consider macroscopic envelope mixing induced by internal gravity waves to compute the modes and study its effect on the surface abundances. Results. We find systematically more modes to be excited for the stellar models including microscopic diffusion compared to those without it, in agreement with observational findings of pulsating B-type dwarfs. Furthermore, the models with microscopic diffusion predict that excited modes occur earlier on in the evolution compared to modes without it. In order to maintain realistic surface abundances during the main sequence, we include macroscopic envelope mixing by internal gravity waves. Along with microscopic diffusion, such macroscopic envelope mixing ensures both more excited modes and surface abundances consistent with spectroscopic studies of B-type stars. Conclusions. While radiative levitation has so far largely been neglected in stellar evolution computations of B-type stars for computational convenience, it impacts mode excitation predictions for stellar models of such stars. We conclude that the process of radiative levitation is able to reduce the discrepancy between predicted and observed excited pulsation modes in B-type stars. [ABSTRACT FROM AUTHOR]

Context. In 2020, the Transiting Exoplanet Survey Satellite (TESS) observed a rapidly rotating M7 dwarf, TIC 277539431, producing a flare at 81° latitude, the highest latitude flare located to date. This is in stark contrast to solar flares that occur much closer to the equator, typically below 30°. The mechanisms that allow flares at high latitudes to occur are poorly understood. Aims. We studied five sectors of TESS monitoring, and obtained 36 ks of XMM-Newton observations to investigate the coronal and flaring activity of TIC 277539431. Methods. From the observations, we infer the optical flare frequency distribution; flare loop sizes and magnetic field strengths; the soft X-ray flux, luminosity, and coronal temperatures; as well as the energy, loop size, and field strength of a large flare in the XMM-Newton observations. Results. We find that the corona of TIC 277539431 does not differ significantly from other low-mass stars on the canonical saturated activity branch with respect to coronal temperatures and flaring activity, but shows lower luminosity in soft X-ray emission by about an order of magnitude, consistent with other late M dwarfs. Conclusions. The lack of X-ray flux, the high-latitude flare, the star's viewing geometry, and the otherwise typical stellar corona taken together can be explained by the migration of flux emergence to the poles in rapid rotators like TIC 277539431 that drain the star's equatorial regions of magnetic flux, but preserve its ability to produce powerful flares. [ABSTRACT FROM AUTHOR]

Aims. The most reliable sources for determining absolute stellar parameters are the double-lined eclipsing binary systems. Some of these systems also show apsidal motion, characterized by the variable log k2. This point grants us the possibility of investigating the stellar interior, specifically the degree of stellar mass concentration. The first studies carried out about four decades ago on this topic showed appreciable discrepancies not only with respect to the comparison between the observed absolute dimensions and their theoretical counterparts, but mainly concerning the degree of mass concentration through the analysis of their apsidal motions. Fortunately, this scenario has been gradually improving with the advances in the quality of observational techniques and advances in the input physics of the evolutionary stellar models (e.g. opacities, thermonuclear reactions, equations of state, numerical techniques). These new developments in the input physics has improved the comparison between observations and the values predicted by theory, including the apsidal motion rates. This progress has lead us to investigate second-order effects such as rotation and dynamic tides. In this paper we deal with the effects of rotation on the degree of mass concentration. Methods. The stellar models were calculated using the MESA package. The mass range studied here was 0.40–25.0 M⊙ for a generic chemical composition characterized by X = 0.70 and Z = 0.02. The present models were computed without taking into account core overshooting in order to highlight the effects of rotation. Each model was followed from the pre-main sequence until the central hydrogen content is of the order of or less than 0.03, covering the range of masses and evolutionary status of the double-lined eclipsing binaries showing apsidal motion. Regarding the calculation of the internal structure coefficients, we integrated the Radau equation using the fifth-order Runge-Kutta method and a tolerance level of 10−7. As an auxiliary tool, homology transformations have been used to explain, qualitatively, how the equation of state, thermonuclear reactions, and convection impact the degree of mass concentration of the models. Results. In our past studies on the effects of rotation on log k2 an average correction was used for all models together. For this a small range of masses (2.0, 7.0, 15.0 M⊙) has been used. In the present paper the corrections due to the effects of rotation in log k2 are presented for each mass individually and for three evolutionary phases. This point is particularly important, given that these corrections show a clear dependence on the mass and on the evolutionary status. Such corrections can be easily introduced into the theoretical calculation of apsidal motion rates through a linear equation. A typical correction due to the rotation effects is of the order of −0.03 in log k2. [ABSTRACT FROM AUTHOR]

Context. The first stars might have been fast rotators. This would have important consequences for their radiative, mechanical, and chemical feedback. Aims. We discuss the impact of fast initial rotation on the evolution of massive Population III models and on their nitrogen and oxygen stellar yields. Methods. We explore the evolution of Population III stars with initial masses in the range of 9 M⊙ ≤ Mini ≤ 120 M⊙, starting with an initial rotation on the zero-age main sequence equal to 70% of the critical one. Results. We find that with the physics of rotation considered here, our rapidly rotating Population III stellar models do not follow a homogeneous evolution. They lose very little mass in the case in which mechanical winds are switched on when the surface rotation becomes equal to or larger than the critical velocity. The impact on the ionising flux appears to be modest when compared to moderately rotating models. Fast rotation favours, in models with initial masses above ∼20 M⊙, the appearance of a very extended intermediate convective zone around the H-burning shell during the core He-burning phase. This shell has important consequences for the sizes of the He- and CO-cores, and thus impacts the final fate of stars. Moreover, it has a strong impact on nucleosynthesis, boosting the production of primary 14N. Conclusions. Fast initial rotation significantly impacts the chemical feedback of Population III stars. Observations of extremely metal-poor stars and/or starbursting regions are essential to provide constraints on the properties of the first stars. [ABSTRACT FROM AUTHOR]

Context. The final stages of a protoplanetary disk are essential for our understanding of the formation and evolution of planets. Photoevaporation is an important mechanism that contributes to the dispersal of an accretion disk and has significant consequences for the disk's lifetime. However, the combined effects of photoevaporation and star-disk interaction have not been investigated in previous studies. Aims. A photoevaporative disk evolution model including a detailed formulation of the inner star-disk interaction region will improve the understanding of the final stages of disk evolution. Methods. We combined an implicit disk evolution model with a photoevaporative mass-loss profile. By including the innermost disk regions down to 0.01 AU, we could calculate the star-disk interaction, the stellar spin evolution, and the transition from an accreting disk to the propeller regime self-consistently. Starting from an early Class II star-disk system (with an age of 1 Myr), we calculated the long-term evolution of the system until the disk becomes almost completely dissolved. Results. Photoevaporation has a significant effect on disk structure and evolution. The radial extent of the dead zone decreases, and the number of episodic accretion events (outbursts) is reduced by high stellar X-ray luminosities. Reasonable accretion rates (>10−8M⊙yr) in combination with photoevaporative gaps are possible for a dead zone that is still massive enough to develop episodic accretion events. Furthermore, the stellar spin evolution during the Class II evolution is less affected by the star-disk interaction in the case of high X-ray luminosities. Conclusions. Our results suggest that the formation of planets, especially habitable planets, in the dead zone is strongly impaired in the case of strong X-ray luminosities. Additionally, the importance of the star-disk interaction during the Class II phase with respect to the stellar spin evolution is reduced. [ABSTRACT FROM AUTHOR]

Context. The manifestation of surface spots on magnetic chemically peculiar (mCP) stars is most commonly explained by the atomic diffusion theory, which requires a calm stellar atmosphere and only moderate rotation. While very successful and well described, this theory still needs to be revised and fine-tuned to the observations. Aims. Our study aims to enlarge the sample of known photometrically variable mCP stars (ACV variables) to pave the way for more robust and significant statistical studies. We derive accurate physical parameters for these objects and discuss our results in the framework of the atomic diffusion theory. Methods. We studied 1314 candidate ACV variables that were selected from the Zwicky Transient Factory catalogue of periodic variables based on light curve characteristics. We investigated these objects using photometric criteria, a colour-magnitude diagram, and spectroscopic data from the LAMOST and Gaia missions to confirm their status as ACV variables. Results. We present a sample of 1232 new ACV variables, including information on distance from the Sun, mass, fractional age on the main sequence, fraction of the radius between the zero-age and terminal-age main sequence, and the equatorial velocity and its ratio to the critical velocity. Conclusions. Our results confirm that the employed selection process is highly efficient for detecting ACV variables. We have identified 38 stars with vequ in excess of 150 km s−1 (with extreme values up to 260 km s−1). This challenges current theories that cannot explain the occurrence of such fast-rotating mCP stars. [ABSTRACT FROM AUTHOR]

We present the discovery of three super-Earth candidates orbiting HD 48948, a bright K-dwarf star with an apparent magnitude of |$m_V$|  = 8.58 mag. As part of the HARPS-N Rocky Planet Search programme, we collect 189 high-precision radial velocity measurements using the HARPS-N spectrograph from 2013 October 6, to 2023 April 16. Various methodologies are applied to extract the radial velocities from the spectra, and we conduct a comprehensive comparative analysis of the outcomes obtained through these diverse extraction techniques. To ensure the robustness of our findings, we employ several methods to address stellar variability, with a focus on Gaussian Process regression. To account for the impact of stellar variability and correlated noise in the radial velocity data set, we include activity indicators, such as |$\log R^{^{\prime }}_{\mathrm{HK}}$| and bisector span, in the multidimensional Gaussian Process regression. Our analysis reveals three planetary candidates with orbital periods of 7.3, 38, and 151 d, and minimum masses estimated at |$4.88 \pm 0.21$| M |$_{\oplus }$|⁠ , |$7.27 \pm 0.70$| M |$_{\oplus }$|⁠ , and |$10.59 \pm 1.00$| M |$_{\oplus }$|⁠ , respectively. The outermost planet resides within the (temperate) habitable zone, positioned at a projected distance of |$0.029\,{\rm arcsec}$| from its star. Given the close proximity of this planetary system, situated at a distance of 16.8 parsecs, HD 48498 emerges as a promising target (closest super-Earth around FGK stars) for future high-contrast direct imaging and high-resolution spectroscopic studies. [ABSTRACT FROM AUTHOR]

Aims. The most reliable sources for determining absolute stellar parameters are the double-lined eclipsing binary systems. Some of these systems also show apsidal motion, characterized by the variable log k2. This point grants us the possibility of investigating the stellar interior, specifically the degree of stellar mass concentration. The first studies carried out about four decades ago on this topic showed appreciable discrepancies not only with respect to the comparison between the observed absolute dimensions and their theoretical counterparts, but mainly concerning the degree of mass concentration through the analysis of their apsidal motions. Fortunately, this scenario has been gradually improving with the advances in the quality of observational techniques and advances in the input physics of the evolutionary stellar models (e.g. opacities, thermonuclear reactions, equations of state, numerical techniques). These new developments in the input physics has improved the comparison between observations and the values predicted by theory, including the apsidal motion rates. This progress has lead us to investigate second-order effects such as rotation and dynamic tides. In this paper we deal with the effects of rotation on the degree of mass concentration. Methods. The stellar models were calculated using the MESA package. The mass range studied here was 0.40–25.0 M⊙ for a generic chemical composition characterized by X = 0.70 and Z = 0.02. The present models were computed without taking into account core overshooting in order to highlight the effects of rotation. Each model was followed from the pre-main sequence until the central hydrogen content is of the order of or less than 0.03, covering the range of masses and evolutionary status of the double-lined eclipsing binaries showing apsidal motion. Regarding the calculation of the internal structure coefficients, we integrated the Radau equation using the fifth-order Runge-Kutta method and a tolerance level of 10−7. As an auxiliary tool, homology transformations have been used to explain, qualitatively, how the equation of state, thermonuclear reactions, and convection impact the degree of mass concentration of the models. Results. In our past studies on the effects of rotation on log k2 an average correction was used for all models together. For this a small range of masses (2.0, 7.0, 15.0 M⊙) has been used. In the present paper the corrections due to the effects of rotation in log k2 are presented for each mass individually and for three evolutionary phases. This point is particularly important, given that these corrections show a clear dependence on the mass and on the evolutionary status. Such corrections can be easily introduced into the theoretical calculation of apsidal motion rates through a linear equation. A typical correction due to the rotation effects is of the order of −0.03 in log k2. [ABSTRACT FROM AUTHOR]

In the standard two-component crust-superfluid model of a neutron star, timing noise can arise when the two components are perturbed by stochastic torques. Here it is demonstrated how to analyse fluctuations in radio pulse times of arrival with a Kalman filter to measure physical properties of the two-component model, including the crust-superfluid coupling time-scale and the variances of the crust and superfluid torques. The analysis technique, validated previously on synthetic data, is applied to observations with the Molonglo Observatory Synthesis Telescope of the representative pulsar PSR J1359−6038. It is shown that the two-component model is preferred to a one-component model, with log Bayes factor 6.81 ± 0.02. The coupling time-scale and the torque variances on the crust and superfluid are measured with 90 per cent confidence to be |$10^{7.1^{+0.8}_{-0.5}}$| |$\rm {s}$| and |$10^{-24.0^{+0.4}_{-5.6}}$| |$\rm {rad^2~s^{-3}}$| and |$10^{-21.7^{+3.5}_{-0.9}}$| |$\rm {rad^2~s^{-3}}$|⁠ , respectively. [ABSTRACT FROM AUTHOR]

Fast radio bursts (FRBs) are cosmological radio transients with millisecond durations and extremely high brightness temperatures. One FRB repeater, FRB 180916.J0158+65 (FRB 180916B), was confirmed to appear 16.35-day periodic activities with 5-day activity window. Another FRB repeater, FRB 121102, and two soft gamma-ray repeaters (SGRs), SGR 1935+2154 and SGR 1806−20, also show possible periodic activities. These periodicities might originate from the precession process of young magnetars due to the anisotropic pressure from the inner magnetic fields as proposed in the literature. In this work, we analyse a self-consistent model for the rotation evolution of magnetars and obtain the evolutions of magnetar precession and obliquity. We find that if the FRB repeaters and the SGRs with (possible) periodic activities originate from the magnetar precession, their ages would be constrained to be hundreds to tens of thousands of years, which is consistent with the typical ages of magnetars. Assuming that the FRB emission is beaming in the magnetosphere as proposed in the literature, we calculate the evolution of the observable probability and the duty cycle of the active window period. We find that for a given magnetar the observable probability increases with the magnetar age in the early stage and decreases with the magnetar age in the later stage; meanwhile, there are one or two active windows in one precession period if the emission is not perfectly axisymmetric with respect to the deformation axis of a magnetar, which could be tested by the future observation for repeating FRB sources. [ABSTRACT FROM AUTHOR]

ABSTRACT: In many rotating fluids, the lowest order force balance is between gravity, pressure, and rotational acceleration (‘GPR’ balance). Terrestrial GPR balance takes the form of geostrophy and hydrostasy, which together yield the terrestrial thermal wind equation. By contrast, stellar GPR balance is an oblateness equation, which determines the departures of the thermal variables from spherical symmetry; its curl yields the ‘stellar thermal wind equation’. In this sense, the stellar thermal wind should be viewed not as a consequence of geostrophy, but of baroclinicity in the oblateness. Here, we treat the full stellar oblateness, including the thermal wind, using pressure coordinates. We derive the generalized stellar thermal wind equation and identify the parameter regime for which it holds. In the case of the Sun, not considering the full oblateness has resulted in conflicting calculations of the theoretical aspherical temperature anomaly. We provide new calculation here and find that the baroclinic anomaly is ∼3–60 times smaller than the barotropic anomaly. Thus, the anomaly from the thermal wind may not be measurable helioseismically, but if measurement were possible, this would potentially yield a new way to bracket the depth of the solar tachocline.

Young α-rich (YAR) stars within the old Galactic thick disc exhibit a dual characteristic of relative youth determined with asteroseismology and abundance enhancement in α elements measured from high-resolution spectroscopy. The youth origin of YAR stars has been proposed to be binary evolution via mass transfer or stellar mergers. If that is the case, YAR stars should spin rapidly and thus be magnetically active, because they are mass and angular momentum gainers. In this study, to seek this binary footprint, we select YAR stars on the main-sequence turn-off or the subgiant branch (MSTO-SGB) from APOGEE DR17, whose ages and projected rotation velocities (v sin  i) can be precisely measured. With APOGEE v sin  i and LAMOST spectra, we find that YAR stars are indeed fast rotators and magnetically active. In addition, we observe low [C/N] ratios and high Gaia RUWE in some YAR stars, suggesting that these MSTO-SGB stars probably have experienced mass transfer from red-giant companions. Our findings underscore that magnetic activity can serve as a valuable tool for probing the binary evolution for other chemically peculiar stars, such as red giants with lithium anomalies and carbon-enhanced metal-poor stars. [ABSTRACT FROM AUTHOR]

We use the magnetic field components measured by Zeeman Doppler imaging (ZDI) to calculate the stellar surface force and torque due to magnetic stresses for the fast rotators σ Ori E, 36 Lyn, and CU Vir, and the slow rotator τ Sco. If we assume the stars have spherical photospheres, the estimated torques give spin-down time-scales no larger than 7 × 105 yr. For σ Ori E, the predicted spin-down time-scale, ≃ 6000 yr, is much less than the observationally measured time-scale of ≃ 106 yr. However, for CU Vir, we find that the spin-down time-scale from its ZDI map is 7 × 105 yr in good agreement with its average rate of spin-down from 1960 to 2010. With the exception of τ Sco, the net force due to magnetic stresses at the stellar surface are large compared to the surface-integrated pressure. We discuss possible reasons for the large values of the forces (and torques), and suggest that the likely explanation is that rotation and the magnetic stresses create significant departures from spherical symmetry. [ABSTRACT FROM AUTHOR]

This study investigates the influence of different braking mechanisms on the formation of three long-period radio pulsars (PSRs J0250+5854, J2251−3711, and J0901−4046): plasma-filled magnetosphere in combination with magnetic field decay, fall-back disc, and r-mode instability. These braking mechanisms can also affect the thermal evolution of pulsars. By comparing the model-predicted values with observational data such as spin periods, period derivatives, and upper limits of the bolometric luminosity of these pulsars, we find that these three braking mechanisms can reasonably explain the spin period and the period derivative within a certain range of parameters for these sources. The model-predicted bolometric luminosity associated with magnetic field dissipation exceeds the upper limit for PSR J0901−4046 but falls below the upper limits for PSR J0250+5854 and PSR J2251−3711. The model-predicted bolometric luminosity within the fall-back disc model exceeds the observed results, whereas the luminosity within the r-mode instability falls below the observed upper limit for these three pulsars. However, a conflict arises when we consider the pulsar radio activity and the accretion phases within the fall-back disc model simultaneously. By combining data from radio and X-ray observations, along with precise measurements of surface thermal emissions, we can enhance our understanding of the braking mechanisms involved in the formation of long-period radio pulsars or constrain key model parameters. Finding more long-period pulsars in the future and conducting multiband observations will enhance our understanding of the formation and nature of long-period radio pulsars. [ABSTRACT FROM AUTHOR]

Context. Binary stars cannot be fully understood without assessing the interaction effects between the two components and the impact of these effects on observational diagnostics. Periastron brightening events, also known as the heartbeat phenomenon, are a clear manifestation of this type of interaction. Aims. We aim to explore the role of tidal shear energy dissipation in stars undergoing periastron brightening events. Methods. We performed a computation from first principles that uses a quasi-hydrodynamic Lagrangian scheme to simultaneously solve the orbital motion and the equations of motion of a 3D grid of volume elements covering the inner, rigidly rotating "core" of a tidally perturbed star. The equations of motion include the gravitational acceleration of both stars, the centrifugal, Coriolis, gas pressure accelerations, and viscous coupling between volume elements. The method is illustrated for a grid of model binary systems with a 10 M⊙ primary that is perturbed by a 6.97 M⊙ companion in eccentric orbits (e = 0 − 0.7). The model is then applied to the heartbeat star MACHO 80.7443.1718. Results. We find an increase by factors ∼10−6–10−3 in tidal shear energy dissipation at periastron, consistent with the majority of observed heartbeat stars. The magnitude of the periastron effect correlates with the degree of departure from synchronicity: stars rotating much faster or much slower than the synchronous rate at periastron present the strongest effect. We confirm that for eccentricities ≤0.3, pseudo-synchronization occurs for 0.8 < ω/Ωave < 1, where Ωave is the average orbital angular velocity. The minimum energy rotation rate (pseudo-synchronism) for e = 0.5 and 0.7 occurs for 1.0 < ω/Ωave < 1.15. The tidal shear energy dissipation model reproduces from first principles the ∼23% maximum brightness enhancement at periastron of MACHO 80.7443.1718. Conclusions. Our results suggest that the magnitude and shape of the heartbeat signal may serve as diagnostics for the internal stellar rotation and turbulent viscosity values. [ABSTRACT FROM AUTHOR]

Context. The theory of rotational and chemical evolution is incomplete, thereby limiting the accuracy of model-dependent stellar mass and age determinations. The γ Doradus (γ Dor) pulsators are excellent points of calibration for the current state-of-the-art stellar evolution models, as their gravity modes probe the physical conditions in the deep stellar interior. Yet, individual asteroseismic modelling of these stars is not always possible because of insufficient observed oscillation modes. Aims. This paper presents a novel method to derive distributions of the stellar mass, age, core-boundary mixing efficiency, and initial rotation rates for γ Dor stars. Methods. We computed a grid of rotating stellar evolution models covering the entire γ Dor instability strip. We then used the observed distributions of the luminosity, effective temperature, buoyancy travel time, and near-core rotation frequency of a sample of 539 stars to assign a statistical weight to each of our models. This weight is a measure of how likely the combination of a specific model is. We then computed weighted histograms to derive the most likely distributions of the fundamental stellar properties. Results. We find that the rotation frequency at zero-age main sequence follows a normal distribution, peaking at around 25% of the critical Keplerian rotation frequency. The probability-density function for extent of the core-boundary mixing zone, given by a factor of fCBM times the local pressure scale height (assuming an exponentially decaying parameterisation), decreases linearly with increasing fCBM. Conclusions. Converting the distribution of fractions of critical rotation at the zero-age main sequence to units of d−1, we find most F-type stars start the main sequence with a rotation frequency between 0.5 d−1 and 2 d−1. Regarding the core-boundary mixing efficiency, we find that it is generally weak in this mass regime. [ABSTRACT FROM AUTHOR]

Stellar population synthesis (SPS) is essential for understanding galaxy formation and evolution. However, the recent discovery of rotation-driven phenomena in star clusters warrants a review of uncertainties in SPS models caused by overlooked factors, including stellar rotation. In this study, we investigated the impact of rotation on SPS specifically using the PARSEC V2.0 rotation model and its implications for high-redshift galaxies with the JWST. Rotation enhances the ultraviolet (UV) flux for up to ∼400 Myr after the starburst, with the slope of UV increasing as the population becomes faster rotating and metal-poorer. Using the Prospector tool, we constructed simulated galaxies and deduce their properties associated with dust and star formation. Our results suggest that rapid rotation models result in a gradual UV slope up to 0.1 dex higher and an approximately 50% increase in dust attenuation for identical wide-band spectral energy distributions. Furthermore, we investigated biases if the stellar population was characterized by rapid rotation and demonstrate that accurate estimation can be achieved for rotation rates up to ωi = 0.6. Accounting for the bias in the case of rapid rotation aligns specific star formation rates more closely with predictions from theoretical models. Notably, this also implies a slightly higher level of dust attenuation than previously anticipated, while still allowing for a "dust-free" interpretation of the galaxy. The impact of rapid rotation SPS models on the rest-UV luminosity function is found to be minimal. Overall, our findings have potentially important implications for comprehending dust attenuation and mass assembly history in the high-redshift Universe. [ABSTRACT FROM AUTHOR]

Context. Estimating the physical states of the surfaces of fast-rotating stars is challenging due to several intrinsic processes, which include radiative flux inhomogeneities on the photosphere induced by rotation and circumstellar signatures in their spectra. The analysis of their spectra ultimately requires the use of synthetic grids of spectra accounting for all these physical processes. Aims. In this paper, we present the 'von ZeiPEl's code for gravity darKening specTRal synthesis' (ZPEKTR) code, which is designed to perform the spectral synthesis of fast-rotating stars, accounting for gravity darkening, limb-darkening effects in the continuum and geometrical deformation induced by fast rotation. Methods. We consider colatitudinal temperature and surface-gravity variations, assuming both the classical prescription developed by von Zeipel and the new formulation by Espinosa-Lara. The code runs either with a rectangular or a triangular mesh on the stellar surface. We compare the temperature and gravitational distribution as a function of the stellar latitude arising from both models. Results. The line profiles of He I 4388, 4471, 4922, and 6678 Å produced with both formalisms are compared at three different rotation rates and illustrate differences in shape and central intensity. We also illustrate the fittings of 31 line spectra of classical Be stars averaged from the Be Stars Observation Survey (BeSOS) database and make a comparison among their apparent physical parameters and ages determined from plane-parallel non-local thermodynamical equilibrium (non-LTE) models and parameters determined from classical von Zeipel models, finding a displacement of more evolved objects towards the zero-age main sequence. We also compare the distributions of projected rotation velocities of these objects obtained with and without the inclusion of gravity-darkening effects with ZPEKTR. Conclusions. We observe a shift of the histogram of rotation velocities calculated accounting for effects of gravity darkening concerning rotation velocities obtained through the fittings with classical plane-parallel non-LTE models. We show that models that do not account for gravity darkening can underestimate the rotation velocity, because the stellar latitudes that contribute the higher velocities are those in the equator with the least radiative flux. We envisage near-future improvements to the code, such as the inclusion of differential rotation and treatment of tidal forces in binary stellar systems. [ABSTRACT FROM AUTHOR]

Context. The Li abundance observed in pre-main sequence and main sequence late-type stars is strongly age-dependent, but also shows a complex pattern depending on several parameters, such as rotation, chromospheric activity, and metallicity. The best way to calibrate these effects, and with the aim of studying Li as an age indicator for FGK stars, is to calibrate coeval groups of stars, such as open clusters (OCs) and associations. Aims. We present a considerable target sample of 42 OCs and associations – with an age range from 1 Myr to 5 Gyr – observed within the Gaia-ESO survey (GES), and using the latest data provided by GES iDR6 and the most recent release of Gaia that was then available, EDR3. As part of this study, we update and improve the membership analysis for all 20 OCs presented in our previous article. Methods. We perform detailed membership analyses for all target clusters to identify likely candidates, using all available parameters provided by GES, complemented with detailed bibliographical searches, and based on numerous criteria: from radial velocity distributions, to the astrometry (proper motions and parallaxes) and photometry provided by Gaia, to gravity indicators (log g and the γ index), [Fe/H] metallicity, and Li content in diagrams of (Li equivalent widths) EW(Li) versus Teff. Results. We obtain updated lists of cluster members for the whole target sample, as well as a selection of Li-rich giant contaminants obtained as an additional result of the membership process. Each selection of cluster candidates was thoroughly contrasted with numerous existing membership studies using data from Gaia to ensure the most robust results. Conclusions. These final cluster selections will be used in the third and last paper of this series, which reports the results of a comparative study characterising the observable Li dispersion in each cluster and analysing its dependence on several parameters, allowing us to calibrate a Li–age relation and obtain a series of empirical Li envelopes for key ages in our sample. [ABSTRACT FROM AUTHOR]

Spectropolarimetic campaigns have established that large-scale magnetic fields are present at the surfaces of approximately 10% of massive dwarf stars. However, there is a dearth of magnetic field measurements for their deep interiors. Asteroseismology of gravity-mode pulsations combined with rotating magneto-hydrodynamical calculations of the early-B main-sequence star HD 43317 constrain its magnetic field strength to be approximately 5 × 105 G just outside its convective core. This proof-of-concept study for magneto-asteroseismology opens a new window into the observational characterisation of magnetic fields inside massive stars. [ABSTRACT FROM AUTHOR]