Solar and Stellar Astrophysics (astro-ph.SR)

Abstract: Our recent multi-band photometric study of the colour width of the lower main sequence of the open cluster M37 has revealed the presence of a sizeable initial chemical composition spread in the cluster. If initial chemical composition spreads are common amongst open clusters, this would have major implications for cluster formation models and the foundation of the chemical tagging technique. Here we present a study of the unevolved main sequence of the open cluster M38, employing Gaia DR3 photometry and astrometry, together with newly acquired Sloan photometry. We have analysed the distribution of the cluster's lower main sequence stars with a differential colour-colour diagram made of combinations of Gaia and Sloan magnitudes, like in the study of M37. We employed synthetic stellar populations to reproduce the observed trend of M38 stars in this diagram, and found that the observed colour spreads can be explained simply by the combined effect of differential reddening across the face of the cluster and the presence of unresolved binaries. There is no need to include in the synthetic sample a spread of initial chemical composition as instead necessary to explain the main sequence of M37. Further photometric investigations like ours, as well as accurate differential spectroscopic analyses on large samples of open clusters, are necessary to understand whether chemical abundance spreads are common among the open cluster population.

Abstract: The shapes of Stokes profiles contain much information about the atmospheric conditions that produced them. However, a variety of different atmospheric structures can produce very similar profiles. Thus, it is important for proper interpretation of observations to have a good understanding of how the shapes of Stokes profiles depend on the underlying atmosphere. An excellent tool in this regard is forward modeling, i.e. computing and studying synthetic spectra from realistic simulations of the solar atmosphere. Modern simulations routinely produce several hundred thousand spectral profiles per snapshot. With such numbers, it becomes necessary to use automated procedures in order to organize the profiles according to their shape. Here we illustrate the use of two complementary methods, k-means and k-Shape, to cluster similarly shaped profiles, and demonstrate how the resulting clusters can be combined with knowledge of the simulation's atmosphere to interpret spectral shapes. We generate synthetic Stokes profiles for the Ca II 854.2 nm line using the Multi3D code from a Bifrost simulation snapshot. We then apply the k-means and k-Shape clustering techniques to group the profiles together according to their shape. We show and compare the classes of profile shapes we retrieve from applying both k-means and k-Shape to our synthetic intensity spectra. We then show the structure of the underlying atmosphere for two particular classes of profile shapes retrieved by the clustering, and demonstrate how this leads to an interpretation for the formation of those profile shapes. Furthermore, we apply both methods to the subset of our profiles containing the strongest Stokes V signals, and demonstrate how k-Shape can be qualitatively better than k-means at retrieving complex profile shapes when using a small number of clusters.

Abstract: We infer initial masses of the synthesized radioactive nickel-56 in a sample of recent Type Ia supernovae applying a new formalism introduced recently by Khatami & Kasen (2019). It is shown that the nickel masses we derive do not differ significantly from previous estimates based on the traditional Arnett-model. We derive the $\beta$ parameter for our sample SNe and show that these are consistent with the fiducial value of $\sim 1.6$ given by Khatami & Kasen (2019) from SN Ia hydrodynamical simulations.

Abstract: The recent pre-print by Saio et al. 2023 argues that the supergiant Betelgeuse is already undergoing carbon burning, based on the assumption that all of its light variations are caused by radial pulsations. However, the angular diameter measurements of the star are in conflict with the stellar radius required by their models, as we show in this note. We discuss the feasibility that the Great Dimming was caused by constructive mode interference using long-term brightness measurements and comment on differences in modeling frameworks adopted in Saio et al. 2023 vs Joyce et al. 2020.

Abstract: We analyse the relation between the activity cycle length and the Rossby number and collected a sample of 44 main sequence stars with well-known activity cycle periods and rotation periods. We find a linear behaviour in the double-logarithmic relation between the Rossby number and cycle period. The bifurcation into a long and a short period branch is clearly real but it depends, empirically, on the colour index B-V, indicating a physical dependence on effective temperature and position on the main sequence. Furthermore, there is also a correlation between cycle length and convective turnover time with the relative depth of the convection zone. Based on this, we derive empirical relations between cycle period and Rossby number, and for the short period cycle branch relations, we estimate a scatter of the relative deviation between 14% and 28% on the long-period cycle branch. With these relations, we obtain a good match with the 10.3 yr period for the well known 11-year solar Schwabe cycle and a long-period branch value of 104 yr for the Gleissberg cycle of the Sun. Finally, we suggest that the cycles on the short-period branch appear to be generated in the deeper layers of the convective zone, while long-period branch cycles seem to be related to fewer deep layers in that zone. We show that for a broader B-V range, the Rossby number is a more suitable parameter for universal relation with cycle-rotation than just the rotation period alone. As proof, we demonstrate that our empirical stellar relations are consistent with the 11-year solar Schwabe cycle, in contrast to earlier studies using just the rotation period in their relations. Previous studies have tried to explain the cycle position of the Sun in the cycle-rotation presentation via other kinds of dynamo, however, in our study, no evidence is found that would suggest another type of dynamo for the Sun and other stars.

Abstract: We simulate the hydrodynamics of the wind flow in the B[e] supergiant binary system GG~Carinae and obtain the mass accretion rate onto the secondary and the observed lightcurve. We find an inhomogeneous Bondi-Hoyle-Lyttleton accretion into a curved accretion tail, and confirm that the accretion rate is modulated along the orbit, with a maximum close to periastron. We show that the accretion itself cannot account for the periodical variation in brightness. Instead, we explain the observed variation in the light curve with absorption by the accretion tail. Our results are in general agreement with previously derived stellar masses, orbital parameters, and the system orientation, but imply that the B[e] supergiant is more luminous. We find an effect related to the orbital motion of the two stars, in which the accretion tail is cut by the primary and we term it the Lizard Autotomy Effect. As part of the effect, the primary is self accreting wind that it ejected earlier. The Lizard Autotomy Effect creates an outwardly expanding spiral shell made up of broken segments. We suggest that such a tail exists in other B[e] supergiant systems and can be the source of the circumstellar material observed in such systems. The accretion also forms a disc around the secondary near periastron that later vanishes. We suggest that the formation of such a disc can launch jets that account for the bipolar structure observed around some B[e] supergiants.

Abstract: Magnetic fields are extremely rare in close, hot binaries, with only 1.5\% of such systems known to contain a magnetic star. The eccentric $\epsilon$ Lupi system stands out in this population as the only close binary in which both stars are known to be magnetic. We report the discovery of strong, variable radio emission from $\epsilon$ Lupi using the upgraded Giant Metrewave Radio Telescope (uGMRT) and the MeerKAT radio telescope.The light curve exhibits striking, unique characteristics including sharp, high-amplitude pulses that repeat with the orbital period, with the brightest enhancement occurring near periastron. The characteristics of the light curve point to variable levels of magnetic reconnection throughout the orbital cycle, making $\epsilon$ Lupi the first known high-mass, main sequence binary embedded in an interacting magnetosphere. We also present a previously unreported enhancement in the X-ray light curve obtained from archival XMM-Newton data. The stability of the components' fossil magnetic fields, the firm characterization of their relatively simple configurations, and the short orbital period of the system make $\epsilon$ Lupi an ideal target to study the physics of magnetospheric interactions. This system may thus help us to illuminate the exotic plasma physics of other magnetically interacting systems such as moon-planet, planet-star, and star-star systems including T Tauri binaries, RS CVn systems, and neutron star binaries.

Abstract: A possible key element for large-scale energy release in the solar corona is an MHD kink instability in a single twisted magnetic flux tube. An initial helical current sheet fragments in a turbulent way into smaller-scale sheets, similarly to a nanoflare storm. As the loop expands in the radial direction during the relaxation process, an unstable loop can disrupt nearby stable loops and trigger an MHD avalanche. Exploratory investigations have been conducted in previous works with relatively simplified loop configurations. Here, we address a more realistic environment that comprehensively accounts for most of the physical effects involved in a stratified atmosphere, typical of an active region. The question is whether the avalanche process will be triggered, with what timescales, and how it will develop, as compared with the original, simpler approach. Three-dimensional MHD simulations describe the interaction of magnetic flux tubes, which have a stratified atmosphere, including chromospheric layers, the thin transition region to the corona, and the related transition from high-beta to low-beta regions. The model also includes the effects of thermal conduction and of optically thin radiation. Our simulations address the case where one flux tube among a few is twisted at the footpoints faster than its neighbours. We show that this flux tube becomes kink unstable first, in conditions in agreement with those predicted by analytical models. It rapidly involves nearby stable tubes, instigating significant magnetic reconnection and dissipation of energy as heat. The heating determines the development of chromospheric evaporation, while the temperature rises up to about 10 MK, close to microflares observations. This work confirms that avalanches are a viable mechanism for the storing and release of magnetic energy in plasma confined in closed coronal loops, as a result of photospheric motions.

Abstract: Linear theory is a well developed framework for characterizing instabilities in weakly collisional plasmas, such as the solar wind. In the previous instalment of this series, we analyzed ~1.5M proton and alpha particle Velocity Distribution Functions (VDFs) observed by Helios I and II to determine the statistical properties of the standard instability parameters such as the growth rate, frequency, the direction of wave propagation, and the power emitted or absorbed by each component, as well as to characterize their behavior with respect to the distance from the Sun and collisional processing. In this work, we use this comprehensive set of instability calculations to train a Machine Learning algorithm consisting of three interlaced components that: 1) predict if an interval is unstable from observed VDF parameters; 2) predict the instability properties for a given unstable VDF; and 3) classify the type of the unstable mode. We use these methods to map the properties in multi-dimensional phase space to find that the parallel-propagating, proton-core-induced Ion Cyclotron mode dominates the young solar wind, while the oblique Fast Magnetosonic mode regulates the proton beam drift in the collisionally old plasma.