Solar and Stellar Astrophysics (astro-ph.SR)

Abstract: We investigate the influence of an accretion disk on the angular momentum (AM) evolution of young M dwarfs, which parameters govern the AM distribution after the disk phase, and whether this leads to a mass-independent distribution of SAM. We find that above an initial rate $\Dot{M}_\mathrm{crit} \sim 10^{-8}~\Msolpyr$ accretion "erases" the initial SAM of M dwarfs during the disk lifetime, and stellar rotation converges to values of SAM that are largely independent of initial conditions. For stellar masses $> 0.3~\mathrm{M_\odot}$, we find that observed initial accretion rates $\Dot{M}_\mathrm{init}$ are comparable to or exceed $\Dot{M}_\mathrm{crit}$. Furthermore, stellar SAM after the disk phase scales with the stellar magnetic field strength as a power-law with an exponent of $-1.1$. For lower stellar masses, $\Dot{M}_\mathrm{init}$ is predicted to be smaller than $\Dot{M}_\mathrm{crit}$ and the initial conditions are imprinted in the stellar SAM after the disk phase. To explain the observed mass-independent distribution of SAM, the stellar magnetic field strength has to range between 20~G and 500~G (700~G and 1500~G) for a 0.1~$\mathrm{M_\odot}$ (0.6~$\mathrm{M_\odot}$) star. These values match observed large-scale magnetic field measurements of young M~dwarfs and the positive relation between stellar mass and magnetic field strength agrees with a theoretically-motivated scaling relation. The scaling law between stellar SAM, mass, and the magnetic field strength is consistent for young stars, where these parameters are constrained by observations. Due to the very limited number of available data, we advocate for efforts to obtain more such measurements. Our results provide new constraints on the relation between stellar mass and magnetic field strength and can be used as initial conditions for future stellar spin models, starting after the disk phase. (shortened)

Abstract: We present a distant T$-$type brown dwarf candidate at $\approx2.55$ kpc discovered in the Cosmic Evolution Early Release Science (CEERS) fields by James Webb Space Telescope (JWST) NIRCam. In addition to the superb sensitivity, we utilised 7 filters from JWST in near-IR and thus is advantageous in finding faint, previously unseen brown dwarfs. From the model spectra in new JWST/NIRCam filter wavelengths, the selection criteria of F115W-F277W$<$-0.8 and F277W-F444W$>$1.1 were chosen to target the spectrum features of brown dwarfs having temperatures from 500K to 1300K. Searching through the data from Early Release Observations (ERO) and Early Release Science (ERS), we find 1 promising candidate in the CEERS field. The result of SED fitting suggested an early T spectral type with a low effective temperature of T$_\text{eff}\approx$1300K, the surface gravity of $\log{g}\approx5.25\text{cm s}^{-2}$, and an eddy diffusion parameter of logK$_{zz}\approx7\text{cm}^2 \text{s}^{-1}$, which indicates an age of $\approx$1.8Gyr and a mass of $\approx0.05$M$_{\odot}$. In contrast to typically found T$-$dwarf within several hundred parsecs, the estimated distance of the source is $\approx2.55$kpc, showing the JWST's power to extend the search to a much larger distance.

Abstract: We present an image classification algorithm using deep learning convolutional neural network architecture, which classifies the morphologies of eclipsing binary systems based on their light curves. The algorithm trains the machine with light curve images generated from the observational data of eclipsing binary stars in contact, detached and semi-detached morphologies, whose light curves are provided by Kepler, ASAS and CALEB catalogues. The structure of the architecture is explained, the parameters of the network layers and the resulting metrics are discussed. Our results show that the algorithm, which is selected among 132 neural network architectures, estimates the morphological classes of an independent validation dataset, 705 real data, with an accuracy of 92%.

Abstract: Nanoflare heating through small-scale magnetic reconnection events is one of the prime candidates to explain heating of the solar corona. However, direct signatures of nanoflares are difficult to determine, and unambiguous observational evidence is still lacking. Numerical models that include accelerated electrons, and can reproduce flaring conditions, are essential in understanding how low-energetic events act as a heating mechanism of the corona, and how such events are able to produce signatures in the spectral lines that can be detected through observations. We investigate the effects of accelerated electrons in synthetic spectra from a 3D radiative magnetohydrodynamics simulation to better understand small-scale heating events and their impact on the solar atmosphere. We synthesised the chromospheric Ca II and Mg II lines and the transition region Si IV resonance lines from a quiet Sun numerical simulation that includes accelerated electrons. We calculated the contribution function to the intensity to better understand how the lines are formed, and what factors are contributing to the detailed shape of the spectral profiles. The synthetic spectra are highly affected by variations in temperature and vertical velocity. Beam heating exceeds conductive heating at the heights where the spectral lines form, indicating that the electrons should contribute to the heating of the lower atmosphere and hence affect the line profiles. However, we find that it is difficult to determine specific signatures from the non-thermal electrons due to the complexity of the atmospheric response to the heating in combination with the relatively low energy output (~1e21 erg/s). Still, our results contribute to the understanding of small-scale heating events in the solar atmosphere, and give further guidance to future observations.

Abstract: Oscillations are abundant in the solar corona. Coronal loop oscillations are typically studied using highly idealised models of magnetic flux tubes. In order to improve our understanding of coronal oscillations, it is necessary to consider the effect of realistic magnetic field topology and density structuring. We analyse the damping of coronal oscillations using a self-consistent 3D radiation-MHD simulation of the solar atmosphere spanning from the convection zone into the corona, the associated oscillation dissipation and heating, and finally the physical processes responsible for the damping and dissipation. The simulated corona formed in such a model does not depend on any prior assumptions about the shape of the coronal loops. We find that the bundle of magnetic loops shows damped transverse oscillations in response to perturbations in two separate instances with oscillation periods of 177 s and 191 s, velocity amplitudes of 10 km/s and 16 km/s and damping times of 176 s and 198 s, respectively. The coronal oscillations lead to the development of velocity shear in the simulated corona resulting in the formation of vortices seen in the velocity field caused by the Kelvin-Helmholtz instability, contributing to the damping and dissipation of the transverse oscillations. The oscillation parameters and evolution observed are in line with the values typically seen in observations of coronal loop oscillations. The dynamic evolution of the coronal loop bundle suggests the models of monolithic and static coronal loops with constant lengths might need to be re-evaluated by relaxing the assumption of highly idealised waveguides.

Abstract: We present the software package binary_c-python which provides a convenient and easy-to-use interface to the binary_c framework, allowing the user to rapidly evolve individual systems and populations of stars. binary_c-python is available on Pip and on GitLab. binary_c-python contains many useful features to control and process the output of binary_c, like by providing binary_c-python with logging statements that are dynamically compiled and loaded into binary_c. Moreover, we have recently added standardised output of events like Roche-lobe overflow or double compact-object formation to binary_c, and automatic parsing and managing of that output in binary_c-python. binary_c-python uses multiprocessing to utilise all the cores on a particular machine, and can run populations with HPC cluster workload managers like HTCondor and Slurm, allowing the user to run simulations on large computing clusters. We provide documentation that is automatically generated based on docstrings and a suite of Jupyter notebooks. These notebooks consist of technical tutorials on how to use binary_c-python and use-case scenarios aimed at doing science. Much of binary_c-python is covered by unit tests to ensure reliability and correctness, and the test coverage is continually increased as the package is improved.

Abstract: Aims. The main goal of this work is to date young open clusters using $\delta$ Sct stars. Seismic indices such as the large separation and the frequency at maximum power can help to constrain the models to better characterise the stars. We propose a reliable method to identify some radial modes, which gives us greater confidence in the constrained models. Methods. We extract the frequency content of a sample of $\delta$ Sct stars belonging to the same open cluster. We estimate the low-order large separation by means of different techniques and the frequency at maximum power for each member of the sample. We use a grid of models built with the typical parameters of $\delta$ Sct stars, including mass, metallicity and rotation as independent variables, and determine the oscillation modes. We select the observed frequencies whose ratios match those of the models. Once we find a range of radial modes matching the observed frequencies, mainly the fundamental mode, we add it to the other seismic parameters to derive the stellar age. Assuming star groups have similar chemistry and age, we estimate their mean age by computing a weighted probability density function fit to the age distribution of the seismically constrained models. Results. We estimate the age of Trumpler 10 to be $30_{-20}{+30}$ Myr, and that of Praesepe to be $580 \pm 230$ Myr. In this latter case, we find two apparent populations of $\delta$ Sct stars in the same cluster, one at $510 \pm 140$ Myr and another at $890 \pm 140$ Myr. This may be due to two different formation events, different rotational velocities of the members in our sample of stars (as rapid rotation may modify the observed large separation), or to membership of unresolved binary systems.

Abstract: Early dust grain growth in protostellar envelopes infalling on young discs has been suggested in recent studies, supporting the hypothesis that dust particles start to agglomerate already during the Class 0/I phase of young stellar objects (YSOs). If this early evolution were confirmed, it would impact the usually assumed initial conditions of planet formation, where only particles with sizes $\lesssim 0.25 \mu$m are usually considered for protostellar envelopes. We aim to determine the maximum grain size of the dust population in the envelope of the Class 0/I protostar L1527 IRS, located in the Taurus star-forming region (140 pc). We use Atacama Large millimetre/sub-millimetre Array (ALMA) and Atacama Compact Array (ACA) archival data and present new observations, in an effort to both enhance the signal-to-noise ratio of the faint extended continuum emission and properly account for the compact emission from the inner disc. Using observations performed in four wavelength bands and extending the spatial range of previous studies, we aim to place tight constraints on the spectral ($\alpha$) and dust emissivity ($\beta$) indices in the envelope of L1527 IRS. We find a rather flat $\alpha \sim$ 3.0 profile in the range 50-2000 au. Accounting for the envelope temperature profile, we derive values for the dust emissivity index, 0.9 < $\beta$ < 1.6, and reveal a tentative, positive outward gradient. This could be interpreted as a distribution of mainly ISM-like grains at 2000 au, gradually progressing to (sub-)millimetre-sized dust grains in the inner envelope, where at R=300 au, $\beta$ = 1.1 +/- 0.1. Our study supports a variation of the dust properties in the envelope of L1527 IRS. We discuss how this can be the result of in-situ grain growth, dust differential collapse from the parent core, or upward transport of disc large grains.

Abstract: The short-period AP Dor eclipsing binary's first in-depth and multiband photometric solutions are presented. We made use of our eight nights of ground-based at a southern hemisphere observatory, and twelve sectors of TESS observations. We extracted eight and 1322 minima from our observations and TESS, respectively. We suggested a new linear ephemeris based on the trend of orbital period variations using the Markov chain Monte Carlo (MCMC) approach. The PHysics Of Eclipsing BinariEs (PHOEBE) Python code and the MCMC approach were used for the light curve analysis. This system did not require a starspot for the light curve solutions. We calculated the absolute parameters of the system using Gaia DR3 parallax method. The orbital angular momentum (J_0) of the AP Dor indicates that this system is located in a region of contact binaries. According to our results, this system is an overcontact binary system with a mass ratio of 0.584, a fillout factor of 48%, and an inclination of 53deg. The positions of AP Dor stars on the Hertzsprung-Russell (HR) diagram are represented.

Abstract: New instrumental and telescopes covering the optical and ultra-violet spectral regions have revealed a range of small-scale dynamic features, many which may be related. For example, the range of spicule-like features hints towards a spectrum of features and not just two types; however, direct observational evidence in terms of tracking spicules across multiple wavelengths are needed in order to provide further insight into the dynamics of the Sun's outer atmosphere. This paper uses H $\alpha$ data obtained with the CRisp Imaging SpectroPolarimeter instrument on the Swedish 1-m Solar Telescope, and in the transition region using the Interface Region Imaging Spectrograph with the SJI 1400 {\AA} channel plus spectral data via the Si IV 1394 {\AA} line to track spicules termed Rapid Blue-shifted Excursions (RBEs). The RBEs as seen in the H $\alpha$ blue-wing images presented here can be sub-divided into two categories; a single or multi-threaded feature. Based on the H $\alpha$ spectra, the features can be divided into events showing broadening and line core absorption, events showing broadening and line core emission, events with a pure blue shifted H $\alpha$ profile without any absorption in the red wing, broadened line profile with the absorption in the blue stronger compared to the red wing. From the RBE-like events which have a Si IV 1394 {\AA} line profile, 78% of them show a Si IV line flux increase. Most of these features show a second broadened Si IV component which is slightly blue-shifted.

Abstract: We present a detailed model atmosphere analysis of 14001 DA white dwarfs from the Montreal White Dwarf Database with ultraviolet photometry from the GALEX mission. We use the 100 pc sample, where the extinction is negligible, to demonstrate that there are no major systematic differences between the best-fit parameters derived from optical only data and the optical + UV photometry. GALEX FUV and NUV data improve the statistical errors in the model fits, especially for the hotter white dwarfs with spectral energy distributions that peak in the UV. Fitting the UV to optical spectral energy distributions also reveals UV-excess or UV-deficit objects. We use two different methods to identify outliers in our model fits. Known outliers include objects with unusual atmospheric compositions, strongly magnetic white dwarfs, and binary white dwarfs, including double degenerates and white dwarf + main-sequence systems. We present a list of 89 newly identified outliers based on GALEX UV data; follow-up observations of these objects will be required to constrain their nature. Several current and upcoming large scale spectroscopic surveys are targeting $>10^5$ white dwarfs. In addition, the ULTRASAT mission is planning an all-sky survey in the NUV band. A combination of the UV data from GALEX and ULTRASAT and optical data on these large samples of spectroscopically confirmed DA white dwarfs will provide an excellent opportunity to identify unusual white dwarfs in the solar neighborhood.