Astronomy & Astrophysics manuscript no. main ©ESO 2021 September 8, 2021 Corona and XUV emission modelling of the Sun and Sun-like stars Munehito Shoda1 and Shinsuke Takasao2 1 National Astronomical Observatory of Japan, National Institutes of Natural Sciences, 2-21-1 Osawa, Mitaka, Tokyo, 181-8588, Japan e-mail: [email protected] 2 Department of Earth and Space Science, Graduate School of Science, Osaka University, Toyonaka, Osaka 560-0043, Japan Received month dd, yyyy; accepted month dd, yyyy ABSTRACT The X-ray and extreme-ultraviolet (EUV) emissions from the low-mass stars significantly affect the evolution of the planetary at- mosphere. However, it is, observationally difficult to constrain the stellar high-energy emission because of the strong interstellar extinction of EUV photons. In this study, we simulate the XUV (X-ray+EUV) emission from the Sun-like stars by extending the solar coronal heating model that self-consistently solves, with sufficiently high resolution, the surface-to-coronal energy transport, turbulent coronal heating, and coronal thermal response by conduction and radiation. The simulations are performed with a range of loop lengths and magnetic filling factors at the stellar surface. With the solar parameters, the model reproduces the observed so- lar XUV spectrum below the Lyman edge, thus validating its capability of predicting the XUV spectra of other Sun-like stars. The model also reproduces the observed nearly-linear relation between the unsigned magnetic flux and the X-ray luminosity. From the simulation runs with various loop lengths and filling factors, we also find a scaling relation, namely log LEUV = 9:93 + 0:67 log LX, where LEUV and LX are the luminosity in the EUV (100 Å < λ ≤ 912 Å) and X-ray (5 Å < λ ≤ 100 Å) range, respectively, in cgs. By assuming a power–law relation between the Rossby number and the magnetic filling factor, we reproduce the renowned relation be- tween the Rossby number and the X-ray luminosity. We also propose an analytical description of the energy injected into the corona, which, in combination with the conventional Rosner–Tucker–Vaiana scaling law, semi-analytically explains the simulation results. This study refines the concepts of solar and stellar coronal heating and derives a theoretical relation for estimating the hidden stellar EUV luminosity from X-ray observations. Key words. Sun: corona – Stars: coronae – Ultraviolet: stars – X-rays: stars 1. Introduction netic field of the star (Pevtsov et al. 2003; Vidotto et al. 2014; Kochukhov et al. 2020), stellar EUV emissions are poorly char- The Sun has an aura of hot plasma called the corona, which has acterised because they are difficult to observe. EUV photons suf- a temperature of a few million Kelvin (Edlén 1943). The image fer from strong absorption by the interstellar medium (Rumph of the corona has been captured by the Atmospheric Imaging et al. 1994), for which stellar EUV spectra are observable only Assembly (Lemen et al. 2012) of the NASA’s Solar Dynamics for nearby stars in the limited range of wavelength (≤ 360 Å, Observatory (Pesnell et al. 2012). The corona is not unique to the Ribas et al. 2005; Johnstone et al. 2021). One thus needs to in- Sun and has been observed to surround low-mass main-sequence directly estimate or reconstruct the EUV spectrum of a star from stars in general (Güdel et al. 1997). Due to its high temperature, other observables. The proposed reconstruction methods include the corona is the principal source of stellar XUV (X-ray + EUV: the inversion of the differential emission measure from UV extreme ultraviolet) emissions. and/or X-ray observations (Sanz-Forcada et al. 2011; Duvvuri Stellar XUV emissions drive the expansion and thermal et al. 2021), the empirical correlation between other observable escape of planetary atmospheres (Vidal-Madjar et al. 2003; lines and XUV emission (Linsky et al. 2014; Youngblood et al. Lecavelier des Etangs et al. 2012; Owen & Wu 2013; Ehren- 2017; France et al. 2018; Sreejith et al. 2020), and/or a combi- reich et al. 2015; Airapetian et al. 2017). Therefore, describing nation of them (Diamond-Lowe et al. 2021). However, the em- the stellar XUV emission in terms of the stellar fundamental pa- pirical estimations of stellar EUV emission are based on obser- rameters (luminosity, mass, radius, etc.) is essential in under- vations with large uncertainty and a limited number of samples, arXiv:2106.08915v2 [astro-ph.SR] 7 Sep 2021 standing the evolution of a planet and its habitability. Since stel- and therefore require further validation from different perspec- lar XUV emissions decrease over time (Güdel et al. 1997; Ribas tives. et al. 2005; Telleschi et al. 2005; Claire et al. 2012; Guinan et al. To circumvent the intrinsic difficulty of stellar EUV obser- 2016) presumably in response to the stellar spin-down (Kraft vation, this study uses the solar-stellar connection to theoreti- 1967; Skumanich 1972; Barnes 2003; Irwin & Bouvier 2009; cally estimate the stellar XUV emission. The solar-stellar theo- Matt et al. 2015) caused by magnetised stellar wind (Weber & retical connection is a natural strategy to predict the stellar prop- Davis 1967; Kawaler 1988; Shoda et al. 2020), the long-term erties. For example, by extending the solar coronal theory, Shi- evolution of the XUV activity of the host star needs to be de- bata & Yokoyama(2002) modelled the X-ray characteristics of scribed as well (Tu et al. 2015; Johnstone et al. 2021). the stellar coronae, which was later extended by Takasao et al. While the observational characteristics of stellar X-ray emis- (2020) considering the size distribution of active regions. In this sions have been established, including their correlations with the study, by extending the solar atmospheric model, the structure rotation (Pallavicini et al. 1981; Wright et al. 2011) and mag- of the upper stellar atmosphere and the XUV emission are pre- Article number, page 1 of 25 A&A proofs: manuscript no. main dicted. To this end, the coronal heating problem must be explic- notation meaning itly solved. To solve the coronal heating problem, the following three r radial distance from the stellar centre issues must be addressed (Klimchuk 2006). 1. Energy generation and transfer: the source of the coronal s coordinate along the flux-tube axis thermal energy probably comes from the magneto- G gravitational constant convection on the surface (Steiner et al. 1998). The energy generation and transfer to the corona, possibly in the form k Boltzmann constant of Alfvén waves (Alfvén 1947; Osterbrock 1961; Kudoh & B Shibata 1999; De Pontieu et al. 2007; McIntosh et al. 2011; mH hydrogen mass Srivastava et al. 2017), needs to be solved. 2. Energy dissipation in the corona: the magnetic (and kinetic) me electron mass energy in the corona needs to dissipate to sustain the high-temperature corona. The field-braiding process must be h Planck constant considered as a promising mechanism of magnetic-energy dissipation (Parker 1972, 1988). M solar mass 3. Thermal response to the heating: the density and temper- R solar radius ature of a coronal loop are determined by the energy bal- ance among heating, conduction and radiation. The Rosner- T solar effective temperature Tucker-Vaiana (RTV) scaling law (Rosner et al. 1978) origi- nates from the thermal response to coronal heating. Thus, the Ω solar angular rotation rate RTV scaling law or its generalised form (Serio et al. 1981; Zhuleku et al. 2020) should be reproduced by the model (An- Table 1. Notations of the coordinates (r and s) and constant parameters. tolin & Shibata 2010). Note that subscript denotes the solar fundamental parameter. For these issues, a model of the coronal heating should 1. in- clude the photosphere (stellar surface) and chromosphere, 2. appropriately consider the field-braiding process in the corona, XUV emission that includes a significant contribution from the and 3. implement the thermal conduction and radiative cooling. transition region. Classically, numerical models of a corona have often focused Considering the difficulty of the TR problem, we perform on the thermal responses to heating events by one-dimensional a series of 1D, high-resolution magnetohydrodynamic (MHD) (1D) (expanding) flux-tube models (Antiochos & Sturrock 1978; simulations for a wide range of parameters. This 1D model fa- Peres et al. 1982; Antiochos et al. 1999). A zero-dimensional cilitates a coronal loop simulation with sufficiently high numeri- model of the coronal thermal evolution is also proposed (Klim- cal resolution. The field-braiding process (or turbulence), which chuk et al. 2008). The advancement of numerical techniques and is essentially 3D, must be appropriately modelled when solving increase in computational power have made models highly so- the coronal heating problem by 1D simulation. In this study, an phisticated. Several models deal with the realistic energy gen- approximated formulation of turbulent dissipation, developed in eration by explicitly solving the magneto-convection (Hansteen previous studies (Dmitruk et al. 2002; Shoda et al. 2018), is em- et al. 2015; Rempel 2017). Other models have focused more ployed as it is likely to reproduce the average heating rate of the on energy dissipation with simpler numerical settings (Moriyasu coronal loop (van Ballegooijen et al. 2011). For simplicity, we et al. 2004; Rappazzo et al. 2008; van Ballegooijen et al. 2011; focus on the Sun-like stars that exhibit solar mass, luminosity, Dahlburg et al. 2016). These models have predicted that the radius, and metallicity. signature of coronal heating could be explained by convection- The rest of the manuscript is structured as follows. In Sec- driven energy injection. By extending these studies, we aimed tion2, the coronal model and numerical methods are detailed. to construct a stellar coronal model that would satisfy the three Section3 produces the numerical results, focusing on the depen- requirements.