No Arabic abstract
... [C]urrent stellar models predict oscillation frequencies that are systematically affected by simplified modelling of the near-surface layers. We use three-dimensional radiation hydrodynamics simulations to better model the near-surface equilibrium structure of dwarfs with spectral types F3, G2, K0 and K5, and examine the differences between oscillation mode frequencies. ... We precisely match stellar models to the simulations gravities and effective temperatures at the surface, and to the temporally- and horizontally-averaged densities and pressures at their deepest points. We then replace the near-surface structure with that of the averaged simulation and compute the change in the oscillation mode frequencies. We also fit the differences using several parametric models currently available in the literature. The surface effect in the stars of solar-type and later is qualitatively similar and changes steadily with decreasing effective temperature. In particular, the point of greatest frequency difference decreases slightly as a fraction of the acoustic cut-off frequency and the overall scale of the surface effect decreases. The surface effect in the hot, F3-type star follows the same trend in scale (i.e. it is larger in magnitude) but shows a different overall variation with mode frequency. We find that the two-term fit by Ball & Gizon (2014) is best able to reproduce the surface terms across all four spectral types, although the scaled solar term and a modified Lorentzian function also match the three cooler simulations reasonably well. ... Our simplified results suggest that the surface effect is generally larger in hotter stars (and correspondingly smaller in cooler stars) and of similar shape in stars of solar type and cooler. However, we cannot presently predict whether this will remain so when other components of the surface effect are included.
The space-borne missions have provided us with a wealth of high-quality observational data that allows for seismic inferences of stellar interiors. This requires the computation of precise and accurate theoretical frequencies, but imperfect modeling of the uppermost stellar layers introduces systematic errors. To overcome this problem, an empirical correction has been introduced by Kjeldsen et al. (2008, ApJ, 683, L175) and is now commonly used for seismic inferences. Nevertheless, we still lack a physical justification allowing for the quantification of the surface-effect corrections. We used a grid of these simulations computed with the CO$^5$BOLD code to model the outer layers of solar-like stars. Upper layers of the corresponding 1D standard models were then replaced by the layers obtained from the horizontally averaged 3D models. The frequency differences between these patched models and the 1D standard models were then calculated using the adiabatic approximation and allowed us to constrain the Kjeldsen et al. power law, as well as a Lorentzian formulation. We find that the surface effects on modal frequencies depend significantly on both the effective temperature and the surface gravity. We further provide the variation in the parameters related to the surface-effect corrections using their power law as well as a Lorentzian formulation. Scaling relations between these parameters and the elevation (related to the Mach number) is also provided. The Lorentzian formulation is shown to be more robust for the whole frequency spectrum, while the power law is not suitable for the frequency shifts in the frequency range above $ u_{rm max}$.
Chemical element transport processes are among the crucial physical processes needed for precise stellar modelling. Atomic diffusion by gravitational settling nowadays is usually taken into account, and is essential for helioseismic studies. On the other hand, radiative accelerations are rarely accounted for, act differently on the various chemical elements, and can strongly counteract gravity in some stellar mass domains. In this study we aim at determining whether radiative accelerations impact the structure of solar-like oscillating main-sequence stars observed by asteroseismic space missions. We implemented the calculation of radiative accelerations in the CESTAM code using the Single-Valued Parameter method. We built and compared several grids of stellar models including gravitational settling, but some with and others without radiative accelerations. We found that radiative accelerations may not be neglected for stellar masses larger than 1.1~M$_{odot}$ at solar metallicity. The difference in age due to their inclusion in models can reach 9% for the more massive stars of our grids. We estimated that the percentage of the PLATO core program stars whose modelling would require radiative accelerations ranges between 33 and 58% depending on the precision of the seismic data. We conclude that, in the context of Kepler, TESS, and PLATO missions, which provide (or will provide) high quality seismic data, radiative accelerations can have a significant effect when inferring the properties of solar-like oscillators properly. This is particularly important for age inferences. However, the net effect for each individual star results from the competition between atomic diffusion including radiative accelerations and other internal transport processes. This will be investigated in a forthcoming companion paper.
Over the last decades, realistic 3D radiative-MHD simulations have become the dominant theoretical tool for understanding the complex interactions between the plasma and the magnetic field on the Sun. Most of such simulations are based on approximations of magnetohydrodynamics, without directly considering the consequences of the very low degree of ionization of the solar plasma in the photosphere and bottom chromosphere. The presence of large amount of neutrals leads to a partial decoupling of the plasma and the magnetic field. As a consequence of that, a series of non-ideal effects (ambipolar diffusion, Hall effect and battery effect) arises. The ambipolar effect is the dominant one in the solar chromosphere. Here we report on the first three-dimensional realistic simulations of magneto-convection including ambipolar diffusion and battery effects. The simulations are done using the newly developed Mancha3D code. Our results reveal that ambipolar diffusion causes measurable effects on the amplitudes of waves excited by convection in the simulations, on the absorption of Poynting flux and heating and on the formation of chromospheric structures. We provide a low limit on the chromospheric temperature increase due to the ambipolar effect using the simulations with battery-excited dynamo fields.
The relevance of the standing accretion shock instability (SASI) compared to neutrino-driven convection in three-dimensional (3D) supernova-core environments is still highly controversial. Studying a 27 Msun progenitor, we demonstrate, for the first time, that violent SASI activity can develop in 3D simulations with detailed neutrino transport despite the presence of convection. This result was obtained with the Prometheus-Vertex code with the same sophisticated neutrino treatment so far used only in 1D and 2D models. While buoyant plumes initially determine the nonradial mass motions in the postshock layer, bipolar shock sloshing with growing amplitude sets in during a phase of shock retraction and turns into a violent spiral mode whose growth is only quenched when the infall of the Si/SiO interface leads to strong shock expansion in response to a dramatic decrease of the mass accretion rate. In the phase of large-amplitude SASI sloshing and spiral motions, the postshock layer exhibits nonradial deformation dominated by the lowest-order spherical harmonics (l=1, m=0,-1,+1) in distinct contrast to the higher multipole structures associated with neutrino-driven convection. We find that the SASI amplitudes, shock asymmetry, and nonradial kinetic energy in 3D can exceed those of the corresponding 2D case during extended periods of the evolution. We also perform parametrized 3D simulations of a 25 Msun progenitor, using a simplified, gray neutrino transport scheme, an axis-free Yin-Yang grid, and different amplitudes of random seed perturbations. They confirm the importance of the SASI for another progenitor, its independence of the choice of spherical grid, and its preferred growth for fast accretion flows connected to small shock radii and compact proto-neutron stars as previously found in 2D setups.
We explore the possibility of the formation of globular clusters under ultraviolet (UV) background radiation. One-dimensional spherical symmetric radiation hydrodynamics (RHD) simulations by Hasegawa et al. have demonstrated that the collapse of low-mass (10^6-10^7 solar masses) gas clouds exposed to intense UV radiation can lead to the formation of compact star clusters like globular clusters (GCs) if gas clouds contract with supersonic infall velocities. However, three-dimensional effects, such as the anisotropy of background radiation and the inhomogeneity in gas clouds, have not been studied so far. In this paper, we perform three-dimensional RHD simulations in a semi-cosmological context, and reconsider the formation of compact star clusters in strong UV radiation fields. As a result, we find that although anisotropic radiation fields bring an elongated shadow of neutral gas, almost spherical compact star clusters can be procreated from a supersonic infall cloud, since photo-dissociating radiation suppresses the formation of hydrogen molecules in the shadowed regions and the regions are compressed by UV heated ambient gas. The properties of resultant star clusters match those of GCs. On the other hand, in weak UV radiation fields, dark matter-dominated star clusters with low stellar density form due to the self-shielding effect as well as the positive feedback by ionizing photons. Thus, we conclude that the supersonic infall under a strong UV background is a potential mechanism to form GCs.