Solar and Stellar Astrophysics (astro-ph.SR)

Abstract: COCONUT is a global coronal magnetohydrodynamic model recently developed. In order to achieve robustness and fast convergence to steady-state, several assumptions have been made during its development, such as prescribing filtered photospheric magnetic maps for representing the magnetic field in the lower corona. This filtering leads to smoothing and lower magnetic field values at the inner boundary, resulting in an unrealistically high plasma beta.In this paper, we examine the effects of prescribing such filtered magnetograms and formulate a method for achieving more realistic plasma beta values without losing computational performance. We demonstrate the effects of the highly pre-processed magnetic maps and the resulting high plasma beta on the features in the domain. Then, in our new approach, we shift the inner boundary to 2 Rs and preserve the prescribed highly filtered magnetic map. This effectively reduces the prescribed plasma beta and leads to a more realistic setup. The method is applied on a magnetic dipole, a minimum (2008) and a maximum (2012) solar activity case, to demonstrate its effects. The results obtained with the proposed approach show significant improvements in the resolved density and radial velocity profiles, and far more realistic values of the plasma \{beta} at the boundary and inside the computational domain. This is also demonstrated via synthetic white light imaging and with the validation against tomography. The computational performance comparison shows similar convergence to a limit residual on the same grid when compared to the original setup. Considering that the grid can be further coarsened with this new setup, the operational performance can be additionally increased if needed. The newly developed method is thus deemed as a good potential replacement of the original setup for operational purposes.

Abstract: The powerful, high-energy magnetic activities of young stars play important roles in the magnetohydrodynamics in the innermost parts of the protoplanetary disks. In addition, the associated UV and X-ray emission dictates the photochemistry; moreover, the corona activities can affect the atmosphere of a newborn extra-solar planet. How the UV and X-ray photons are generated, and how they illuminate the disks, are not well understood. Here we report the analyses of the optical and infrared (OIR) photometric monitoring observations and the high angular-resolution centimeter band images of the low-mass (M1 type) pre-main sequence star, DM Tau. We found that the OIR photometric light curves present periodic variations, which is consistent with that the host young star is rotating in the same direction as the natal disk and is hosting at least one giant cold spot. In addition, we resolved that the ionized gas in the DM Tau disk is localized, and its spatial distribution is varying with time. All the present observations can be coherently interpreted, if the giant cold spot is the dominant anisotropic UV and/or X-ray source that illuminates the ambient cone-like region. These results indicate that a detailed theoretical model of the high-energy protostellar emission is essential in the understanding of the space weather around the extra-solar planets and the origin of life.

Abstract: The importance of the existence of a radiative core in generating a solar-like magnetic dynamo is still unclear. Analytic models and magnetohydrodynamic simulations of stars suggest the thin layer between a star's radiative core and its convective zone can produce shearing that reproduces key characteristics of a solar-like dynamo. However, recent studies suggest fully and partially convective stars exhibit very similar period-activity relations, hinting that dynamos generated by stars with and without radiative cores hold similar properties. Here, using kinematic ages, we discover an abrupt change in the stellar spin-down law across the fully convective boundary. We found that fully convective stars exhibit a higher angular momentum loss rate, corresponding to a torque that is $\sim$ 2.25 times higher for a given angular velocity than partially convective stars around the fully convective boundary. This requires a dipole field strength that is larger by a factor of $\sim$2.5, a mass loss rate that is $\sim$4.2 times larger, or some combination of both of those factors. Since stellar-wind torques depend primarily on large-scale magnetic fields and mass loss rates, both of which derive from magnetic activity, the observed abrupt change in spin-down law suggests that the dynamos of partially and fully convective stars may be fundamentally different

Abstract: Interplanetary solar radio type III bursts provide the means for remotely studying and tracking energetic electrons propagating in the interplanetary medium. Due to the lack of direct radio source imaging, several methods have been developed to determine the source positions from space-based observations. Moreover, none of the methods consider the propagation effects of anisotropic radio-wave scattering, which would strongly distort the trajectory of radio waves, delay their arrival times, and affect their apparent characteristics. We investigate the source positions and directivity of an interplanetary type III burst simultaneously observed by Parker Solar Probe, Solar Orbiter, STEREO, and Wind and compare the results of applying the intensity fit and timing methods with ray-tracing simulations of radio-wave propagation with anisotropic density fluctuations. The simulation calculates the trajectories of the rays, their time profiles at different viewing sites, and the apparent characteristics for various density fluctuation parameters. The results indicate that the observed source positions are displaced away from the locations where emission is produced, and their deduced radial distances are larger than expected from density models. This suggests that the apparent position is affected by anisotropic radio-wave scattering, which leads to an apparent position at a larger heliocentric distance from the Sun. The methods to determine the source positions may underestimate the apparent positions if they do not consider the path of radio-wave propagation and incomplete scattering at a viewing site close to the intrinsic source position.

Abstract: NSV 14264 and NSV 14172 are suspected to be variable stars of RR Lyr type (Brun, 1964). They were observed during three nights in October 2018 with a 25cm diameter telescope. These observations completed by ASAS-SN survey data bring to the conclusion that these two stars are not RR Lyraes but constant stars in the limit of the precision of the present photometry. The analysis of GAIA data allows to say that NSV 14264 is a main sequence dwarf similar to the Sun but that NSV 14172 is a yellow giant star located in the HR diagram at the limit between RR Lyraes and CW cepheids; however, it does not pulsate with significant amplitude.

Abstract: YZ Ret is the first X-ray flash detected classical nova, and is also well observed in optical, X-ray, and gamma-ray. We propose a comprehensive model that explains the observational properties. The white dwarf mass is determined to be $\sim 1.33 ~M_\odot$ that reproduces multiwavelength light curves of YZ Ret, from optical, X-ray, and to gamma-ray. We show that a shock is naturally generated far outside the photosphere because winds collide with themselves. The derived lifetime of shock explains some of the temporal variations of emission lines. The shocked shell significantly contributes to the optical flux in the nebular phase. The decline trend of shell emission in the nebular phase is close to $\propto t^{-1.75}$ and the same as the universal decline law of classical novae, where $t$ is the time from the outburst.

Abstract: Aims. We tested the new atomic model using atmospheric models of stars of different spectral types: the Sun (dG2), HD22049 (dK2, Epsilon Eridani), GJ 832 (dM2), and GJ 581 (dM3). Methods. We used new Breit-Pauli distorted-wave (DW) multiconfiguration calculations, which proved to be relevant for many transitions in the mid-infrared (MIR) range. The new atomic model of Mg I includes the following: i) recomputed ECS data through the DW method, including the superlevels. ii) For the nonlocal thermodynamic equilibrium (NLTE) population calculations, 5676 theoretical transitions were added (3001 term-to-term). iii) All of these improvements were studied in the Sun and the stars listed above. Results. The Mg distribution between ionization states for stars with different effective temperatures was compared. For the Sun and Epsilon Eridani, Mg II predominates with more than 95 %, while for GJ 832 and GJ 581, Mg I represents more than 72 % of the population. Moreover, in the latter stars, the amount of Mg forming molecules in their atmosphere is at least two orders of magnitude higher. Regarding the NLTE population, a noticeable lower variability in the departure coefficients was found, indicating a better population coupling for the new model. Comparing the synthetic spectrum calculated with the older and new Mg I atomic model, these results show minimal differences in the visible range but they are stronger in the IR for all of the stars. This aspect should be considered when using lines from this region as indicators. Nevertheless, some changes in the spectral type were found, also emphasizing the need to test the atomic models in different atmospheric conditions. The most noticeable changes occurred in the FUV and NUV, obtaining a higher flux for the new atomic model regardless of the spectral type.

Abstract: We present our findings on the spectral analysis of seven magnetic white dwarfs that were presumed to be double degenerates. We obtained time-resolved spectroscopy at the Gemini Observatory to look for evidence of binarity or fast rotation. We find three of our targets have rotation periods of less than an hour based on the shifting positions of the Zeeman-split H$\alpha$ components: 13, 35, and 39 min, and we find one more target with a ~hour long period that is currently unconstrained. We use offset dipole models to determine the inclination, magnetic field strength, and dipole offset of each target. The average surface field strengths of our fast rotators vary by 1-2 MG between different spectra. In all cases, the observed absorption features are too shallow compared to our models. This could be due to extra flux from a companion for our three low-mass targets, but the majority of our sample likely requires an inhomogeneous surface composition. Including an additional magnetic white dwarf with similar properties presented in the literature, we find that 5 of the 8 targets in this sample show field variations on minute/hour timescales. A crystallization driven dynamo can potentially explain the magnetic fields in three of our targets with masses above $0.7~M_{\odot}$ but another mechanism is still needed to explain their rapid rotation. We suggest that rapid rotation or low-masses point to binary evolution as the likely source of magnetism in 7 of these 8 targets.

Abstract: White dwarf stars are the most common stellar fossils. When in binaries, they make up the dominant form of compact object binary within the Galaxy and can offer insight into different aspects of binary formation and evolution. One of the most remarkable white dwarf binary systems identified to date is AR Scorpii (henceforth AR Sco). AR Sco is composed of an M-dwarf star and a rapidly-spinning white dwarf in a 3.56-hour orbit. It shows pulsed emission with a period of 1.97 minutes over a broad range of wavelengths, which led to it being known as a white dwarf pulsar. Both the pulse mechanism and the evolutionary origin of AR Sco provide challenges to theoretical models. Here we report the discovery of the first sibling of AR Sco, J191213.72-441045.1 (henceforth J1912-4410), which harbours a white dwarf in a 4.03-hour orbit with an M-dwarf and exhibits pulsed emission with a period of 5.30 minutes. This discovery establishes binary white dwarf pulsars as a class and provides support for proposed formation models for white dwarf pulsars.