Self-consistent 2D models of fast rotating early-type star

This work aims at presenting the first two-dimensional models of an isolated rapidly rotating star that include the derivation of the differential rotation and meridional circulation in a self-consistent way.We use spectral methods in multidomains, together with a Newton algorithm to determine the steady state solutions including differential rotation and meridional circulation for an isolated non-magnetic, rapidly rotating early-type star. In particular we devise an asymptotic method for small Ekman numbers (small viscosities) that removes the Ekman boundary layer and lifts the degeneracy of the inviscid baroclinic solutions.For the first time, realistic two-dimensional models of fast-rotating stars are computed with the actual baroclinic flows that predict the differential rotation and the meridional circulation for intermediate-mass and massive stars. These models nicely compare with available data of some nearby fast-rotating early-type stars like Ras Alhague ($\alpha$ Oph), Regulus ($\alpha$ Leo), and Vega ($\alpha$ Lyr). It is shown that baroclinicity drives a differential rotation with a slow pole, a fast equator, a fast core, and a slow envelope. The differential rotation is found to increase with mass, with evolution (here measured by the hydrogen mass fraction in the core), and with metallicity. The core-envelope interface is found to be a place of strong shear where mixing will be efficient.Two-dimensional models offer a new view of fast-rotating stars, especially of their differential rotation, which turns out to be strong at the core-envelope interface. They also offer more accurate models for interpreting the interferometric and spectroscopic data of early-type stars.Comment: 16 pages, 17 figures, to appear in Astronomy and Astrophysic

An algorithm for computing the 2D structure of fast rotating stars

Stars may be understood as self-gravitating masses of a compressible fluid whose radiative cooling is compensated by nuclear reactions or gravitational contraction. The understanding of their time evolution requires the use of detailed models that account for a complex microphysics including that of opacities, equation of state and nuclear reactions. The present stellar models are essentially one-dimensional, namely spherically symmetric. However, the interpretation of recent data like the surface abundances of elements or the distribution of internal rotation have reached the limits of validity of one-dimensional models because of their very simplified representation of large-scale fluid flows. In this article, we describe the ESTER code, which is the first code able to compute in a consistent way a two-dimensional model of a fast rotating star including its large-scale flows. Compared to classical 1D stellar evolution codes, many numerical innovations have been introduced to deal with this complex problem. First, the spectral discretization based on spherical harmonics and Chebyshev polynomials is used to represent the 2D axisymmetric fields. A nonlinear mapping maps the spheroidal star and allows a smooth spectral representation of the fields. The properties of Picard and Newton iterations for solving the nonlinear partial differential equations of the problem are discussed. It turns out that the Picard scheme is efficient on the computation of the simple polytropic stars, but Newton algorithm is unsurpassed when stellar models include complex microphysics. Finally, we discuss the numerical efficiency of our solver of Newton iterations. This linear solver combines the iterative Conjugate Gradient Squared algorithm together with an LU-factorization serving as a preconditionner of the Jacobian matrix.Comment: 40 pages, 12 figures, accepted in J. Comput. Physic

Physical processes leading to surface inhomogeneities: the case of rotation

In this lecture I discuss the bulk surface heterogeneity of rotating stars, namely gravity darkening. I especially detail the derivation of the omega-model of Espinosa Lara & Rieutord (2011), which gives the gravity darkening in early-type stars. I also discuss the problem of deriving gravity darkening in stars owning a convective envelope and in those that are members of a binary system.Comment: 23 pages, 11 figure, Lecture given to the school on the cartography of the Sun and the stars (May 2014 in Besan\c{c}on), to appear in LNP, Neiner and Rozelot edts V2: typos correcte

Non-classical light state transfer in su(2)su(2) resonator networks

We use a normal mode approach to show full and partial state transfer in a class of coupled resonator networks with underlying $su(2)$ symmetry that includes the so-called $J_{x}$ photonic lattice. Our approach defines an auxiliary Hermitian coupling matrix describing the network that yields the normal modes of the system and its time evolution in terms of orthogonal polynomials. These results provide insight on the full quantum state reconstruction time in a general $su(2)$ network of any size and the full quantum transfer time in the $J_{x}$ network of size $4 n + 1$ with $n=1,2,3,\ldots$ In the latter, our approach shows that the Fock state probability distribution of the initial state is conserved but the amplitudes suffer a phase shift proportional to $\pi/2$ that results in partial quantum state transfer for any other network size.Comment: 13 pages, 1 figur

Solar energetic electron events measured by MESSENGER and Solar Orbiter. Peak intensity and energy spectrum radial dependences: statistical analysis

Context/Aims: We present a list of 61 solar energetic electron (SEE) events measured by the MESSENGER mission and the radial dependences of the electron peak intensity and the peak-intensity energy spectrum. The analysis comprises the period from 2010 to 2015, when MESSENGER heliocentric distance varied between 0.31 and 0.47 au. We also show the radial dependencies for a shorter list of 12 SEE events measured in February and March 2022 by spacecraft near 1 au and by Solar Orbiter around its first close perihelion at 0.32 au. Results: Due to the elevated background intensity level of the particle instrument on board MESSENGER, the SEE events measured by this mission are necessarily large and intense; most of them accompanied by a CME-driven shock, being widespread in heliolongitude, and displaying relativistic ($\sim$1 MeV) electron intensity enhancements. The two main conclusions derived from the analysis of the large SEE events measured by MESSENGER, which are generally supported by Solar Orbiter's data results, are: (1) There is a wide variability in the radial dependence of the electron peak intensity between $\sim$0.3 au and $\sim$1 au, but the peak intensities of the energetic electrons decrease with radial distance from the Sun in 27 out of 28 events. On average and within the uncertainties, we find a radial dependence consistent with $R^{-3}$. (2) The electron spectral index found in the energy range around 200 keV ($\delta$200) of the backward-scattered population near 0.3 au measured by MESSENGER is harder in 19 out of 20 (15 out of 18) events by a median factor of $\sim$20% ($\sim$10%) when comparing to the anti-sunward propagating beam (backward-scattered population) near 1 au.Comment: 20 pages, 13 figure

Polarization due to rotational distortion in the bright star Regulus

This is the full published article (retrieved from the 6 months post-publication posting on arXiv) including the Methods and Supplementary Information sections: 33 pages, 10 figures, 8 tablesPolarization in stars was first predicted by Chandrasekhar [1] who calculated a substantial linear polarization at the stellar limb for a pure electron-scattering atmosphere. This polarization will average to zero when integrated over a spherical star but could be detected if the symmetry is broken, for example by the eclipse of a binary companion. Nearly 50 years ago, Harrington and Collins [2] modeled another way of breaking the symmetry and producing net polarization - the distortion of a rapidly rotating hot star. Here we report the first detection of this effect. Observations of the linear polarization of Regulus, with two different high-precision polarimeters, range from +42 parts-per-million (ppm) at a wavelength of 741 nm to -22 ppm at 395 nm. The reversal from red to blue is a distinctive feature of rotation-induced polarization. Using a new set of models for the polarization of rapidly rotating stars we find that Regulus is rotating at 96.5(+0.6/-0.8)% of its critical angular velocity for breakup, and has an inclination greater than 76.5 degrees. The rotation axis of the star is at a position angle of 79.5+/-0.7 degrees. The conclusions are independent of, but in good agreement with, the results of previously published interferometric observations of Regulus [3]. The accurate measurement of rotation in early-type stars is important for understanding their stellar environments [4], and course of their evolution [5].Peer reviewedFinal Accepted Versio

First near-relativistic solar electron events observed by EPD onboard Solar Orbiter

Context. Solar Orbiter, launched in February 2020, started its cruise phase in June 2020, in coincidence with its first perihelion at 0.51 au from the Sun. The in situ instruments onboard, including the Energetic Particle Detector (EPD), operate continuously during the cruise phase enabling the observation of solar energetic particles. Aims. In situ measurements of the first near-relativistic solar electron events observed in July 2020 by EPD are analyzed and the solar origins and the conditions for the interplanetary transport of these particles investigated. Methods. Electron observations from keV energies to the near-relativistic range were combined with the detection of type III radio bursts and extreme ultraviolet (EUV) observations from multiple spacecraft in order to identify the solar origin of the electron events. Electron anisotropies and timing as well as the plasma and magnetic field environment were evaluated to characterize the interplanetary transport conditions. Results. All electron events were clearly associated with type III radio bursts. EUV jets were also found in association with all of them except one. A diversity of time profiles and pitch-angle distributions was observed. Different source locations and different magnetic connectivity and transport conditions were likely involved. The July 11 event was also detected by Wind, separated 107 degrees in longitude from Solar Orbiter. For the July 22 event, the Suprathermal Electron and Proton (STEP) sensor of EPD allowed for us to not only resolve multiple electron injections at low energies, but it also provided an exceptionally high pitch-angle resolution of a very anisotropic beam. This, together with radio observations of local Langmuir waves suggest a very good magnetic connection during the July 22 event. This scenario is challenged by a high-frequency occultation of the type III radio burst and a nominally non-direct connection to the source; therefore, magnetic connectivity requires further investigation