Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn improved parametric method for evaluating radiative accelerations in stellar interiors
58
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The single-valued parameter (SVP) method is a parametric method that offers the possibility of computing radiative accelerations in stellar interiors much faster than other methods. It has been implemented in a few stellar evolution numerical codes for about a decade. In the present paper, we describe improvements we have recently brought in the process of preparing, from atomic/opacity databases, the SVP tables that are needed to use the method, and their extension to a larger stellar mass domain (from 1 to 10 solar mass) on the main-sequence. We discuss the validity domain of the method. We also present the website from where new tables and codes can be freely accessed and implemented in stellar evolution codes.
rate research
Read More
A brief review of various methods to calculate radiative accelerations for stellar evolution and an analysis of their limitations are followed by applications to Pop I and Pop II stars. Recent applications to Horizontal Branch (HB) star evolution are also described. It is shown that models including atomic diffusion satisfy Schwarzschilds criterion on the interior side of the core boundary on the HB without the introduction of overshooting. Using stellar evolution models starting on the Main Sequence and calculated throughout evolution with atomic diffusion, radiative accelerations are shown to lead to abundance anomalies similar to those observed on the HB of M15.
Atomic diffusion has been recognized as an important process that has to be considered in any computations of stellar models. In solar-type and cooler stars, this process is dominated by gravitational settling, which is now included in most stellar evolution codes. In hotter stars, radiative accelerations compete with gravity and become the dominant ingredient in the diffusion flux for most heavy elements. Introducing radiative accelerations into the computations of stellar models modifies the internal element distribution and may have major consequences on the stellar structure. Coupling these processes with hydrodynamical stellar motions has important consequences that need to be investigated in detail. We aim to include the computations of radiative accelerations in a stellar evolution code (here the TGEC code) using a simplified method (SVP) so that it may be coupled with sophisticated macroscopic motions. We also compare the results with those of the Montreal code in specific cases for validation and study the consequences of these coupled processes on accurate models of A- and early-type stars. We implemented radiative accelerations computations into the Toulouse-Geneva stellar evolution code following the semi-analytical prescription proposed by Alecian and LeBlanc. This allows more rapid computations than the full description used in the Montreal code. We present results for A-type stellar models computed with this updated version of TGEC and compare them with similar published models obtained with the Montreal evolution code. We discuss the consequences for the coupling with macroscopic motions, including thermohaline convection.
We study the effects of rotation on the growth and saturation of the double-diffusive fingering (thermohaline) instability at low Prandtl number. Using direct numerical simulations, we estimate the compositional transport rates as a function of the relevant non-dimensional parameters - the Rossby number, inversely proportional to the rotation rate, and the density ratio which measures the relative thermal and compositional stratifications. Within our explored range of parameters, we generally find rotation to have little effect on vertical transport. However, we also present one exceptional case where a cyclonic large scale vortex (LSV) is observed at low density ratio and fairly low Rossby number. The LSV leads to significant enhancement in the fingering transport rates by concentrating compositionally dense downflows at its core. We argue that the formation of such LSVs could be relevant to solving the missing mixing problem in RGB stars.
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.
Context: There is a wide discrepancy in current estimates of the strength of convection flows in the solar interior obtained using different helioseismic methods applied to observations from SDO/HMI. The cause for these disparities is not known. Aims: As one step in the effort to resolve this discrepancy, we aim to characterize the multi-ridge fitting code for ring-diagram helioseismic analysis that is used to obtain flow estimates from local power spectra of solar oscillations. Methods: We updated the multi-ridge fitting code developed by Greer et al.(2014) to solve several problems we identified through our inspection of the code. In particular, we changed the merit function to account for the smoothing of the power spectra, model for the power spectrum, and noise estimates. We used Monte Carlo simulations to generate synthetic data and to characterize the noise and bias of the updated code by fitting these synthetic data. Results: The bias in the output fit parameters, apart from the parameter describing the amplitude of the p-mode resonances in the power spectrum, is below what can be measured from the Monte-Carlo simulations. The amplitude parameters are underestimated; this is a consequence of choosing to fit the logarithm of the averaged power. We defer fixing this problem as it is well understood and not significant for measuring flows in the solar interior. The scatter in the fit parameters from the Monte-Carlo simulations is well-modeled by the formal error estimates from the code. Conclusions: We document and demonstrate a reliable multi-ridge fitting method for ring-diagram analysis. The differences between the updated fitting results and the original results are less than one order of magnitude and therefore we suspect that the changes will not eliminate the aforementioned orders-of-magnitude discrepancy in the amplitude of convective flows in the solar interior.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
