Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEquation of state SAHA-S meets stellar evolution code CESAM2k
62
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present an example of an interpolation code of the SAHA-S equation of state that has been adapted for use in the stellar evolution code CESAM2k. The aim is to provide the necessary data and numerical procedures for its implementation in a stellar code. A technical problem is the discrepancy between the sets of thermodynamic quantities provided by the SAHA-S equation of state and those necessary in the CESAM2k computations. Moreover, the independent variables in a practical equation of state (like SAHA-S) are temperature and density, whereas for modelling calculations the variables temperature and pressure are preferable. Specifically for the CESAM2k code, some additional quantities and their derivatives must be provided. To provide the bridge between the equation of state and stellar modelling, we prepare auxiliary tables of the quantities that are demanded in CESAM2k. Then we use cubic spline interpolation to provide both smoothness and a good approximation of the necessary derivatives. Using the B-form of spline representation provides us with an efficient algorithm for three-dimensional interpolation. The table of B-spline coefficients provided can be directly used during stellar model calculations together with the module of cubic spline interpolation. This implementation of the SAHA-S equation of state in the CESAM2k stellar structure and evolution code has been tested on a solar model evolved to the present. A comparison with other equations of state is briefly discussed. The choice of a regular net of mesh points for specific primary quantities in the SAHA-S equation of state, together with accurate and consistently smooth tabulated values, provides an effective algorithm of interpolation in modelling calculations. The proposed module of interpolation procedures can be easily adopted in other evolution codes.
rate research
Read More
We use the multi-epoch radial velocities acquired by the APOGEE survey to perform a large scale statistical study of stellar multiplicity for field stars in the Milky Way, spanning the evolutionary phases between the main sequence and the red clump. We show that the distribution of maximum radial velocity shifts (drvm) for APOGEE targets is a strong function of logg, with main sequence stars showing drvm as high as $sim$300 kms, and steadily dropping down to $sim$30 kms for logg$sim$0, as stars climb up the Red Giant Branch (RGB). Red clump stars show a distribution of drvm values comparable to that of stars at the tip of the RGB, implying they have similar multiplicity characteristics. The observed attrition of high drvm systems in the RGB is consistent with a lognormal period distribution in the main sequence and a multiplicity fraction of 0.35, which is truncated at an increasing period as stars become physically larger and undergo mass transfer after Roche Lobe Overflow during H shell burning. The drvm distributions also show that the multiplicity characteristics of field stars are metallicity dependent, with metal-poor ([Fe/H]$lesssim-0.5$) stars having a multiplicity fraction a factor 2-3 higher than metal-rich ([Fe/H]$gtrsim0.0$) stars. This has profound implications for the formation rates of interacting binaries observed by astronomical transient surveys and gravitational wave detectors, as well as the habitability of circumbinary planets.
The Saha equation describes the relative number density of consecutive ionization levels of a given atomic species under conditions of thermodynamic equilibrium in an ionized gas. Because the number density in the denominator may be very small, special steps must be taken to ensure numerical stability. In this paper we recast the equation into a form in which each ionization fraction is normalized by the total number density of the atomic species, analogous to the Boltzmann equation describing the distribution of excitation states for a given ion.
We present the updated version of the code used to compute stellar evolutionary tracks in Padova. It is the result of a thorough revision of the major input physics, together with the inclusion of the pre-main sequence phase, not present in our previous releases of stellar models. Another innovative aspect is the possibility of promptly generating accurate opacity tables fully consistent with any selected initial chemical composition, by coupling the OPAL opacity data at high temperatures to the molecular opacities computed with our AESOPUS code (Marigo & Aringer 2009). In this work we present extended sets of stellar evolutionary models for various initial chemical compositions, while other sets with different metallicities and/or different distributions of heavy elements are being computed. For the present release of models we adopt the solar distribution of heavy elements from the recent revision by Caffau et al. (2011), corresponding to a Suns metallicity Z=0.0152. From all computed sets of stellar tracks, we also derive isochrones in several photometric systems. The aim is to provide the community with the basic tools to model star clusters and galaxies by means of population synthesis techniques.
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.
The first gravitational wave detections of mergers between black holes and neutron stars represent a remarkable new regime of high-energy transient astrophysics. The signals observed with LIGO-Virgo detectors come from mergers of extreme physical objects which are the end products of stellar evolution in close binary systems. To better understand their origin and merger rates, we have performed binary population syntheses at different metallicities using the new grid-based binary population synthesis code ComBinE. Starting from newborn pairs of stars, we follow their evolution including mass loss, mass transfer and accretion, common envelopes and supernova explosions. We apply the binding energies of common envelopes based on dense grids of detailed stellar structure models, make use of improved investigations of the subsequent Case BB Roche-lobe overflow and scale supernova kicks according to the stripping of the exploding stars. We demonstrate that all the double black hole mergers, GW150914, LVT151012, GW151226, GW170104, GW170608 and GW170814, as well as the double neutron star merger GW170817, are accounted for in our models in the appropriate metallicity regime. Our binary interaction parameters are calibrated to match the accurately determined properties of Galactic double neutron star systems, and we discuss their masses and types of supernova origin. Using our default values for the input physics parameters, we find a double neutron star merger rate of about 3.0 Myr^-1 for Milky-Way equivalent galaxies. Our upper limit to the merger-rate density of double neutron stars is R=400 yr^-1 Gpc^-3 in the local Universe (z=0).
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
