Normalized to: Sonoi, T.
[1]
oai:arXiv.org:1811.05229 [pdf] - 1815045
Calibration of the mixing-length parameter $\alpha$ for the MLT and FST
models by matching with CO$^5$BOLD models
Submitted: 2018-11-13
The CoRoT and Kepler missions provided a wealth of high-quality data for
solar-like oscillations. To make the best of such data for seismic inferences,
we need theoretical models with precise near-surface structure, which has
significant influence on solar-like oscillation frequencies. The mixing-length
parameter, $\alpha$, is a key factor for the near-surface structure. In the
convection formulations used in evolution codes, the $\alpha$ is a free
parameter that needs to be properly specified. We calibrated $\alpha$ values by
matching entropy profiles of 1D envelope models with those of 3D CO$^5$BOLD
models. For such calibration, previous works concentrated on the classical
mixing-length theory (MLT). Here we also analyzed the full spectrum turbulence
(FST) models. For the atmosphere part in the 1D models, we use the Eddington
grey $T(\tau)$ relation and the one with the solar-calibrated Hopf-like
function. For both the MLT and FST models with a mixing length $l=\alpha H_p$,
calibrated $\alpha$ values increase with increasing $g$ or decreasing $T_{\rm
eff}$. For the FST models, we also calibrated values of $\alpha^*$ defined as
$l=r_{\rm top}-r+\alpha^*H_{p,{\rm top}}$. $\alpha^*$ is found to increase with
$T_{\rm eff}$ and $g$. As for the correspondence to the 3D models, the solar
Hopf-like function gives a photospheric-minimum entropy closer to a 3D model
than the Eddington $T(\tau)$. The structure below the photosphere depends on
the convection model. However, not a single convection model gives the best
correspondence since the averaged 3D quantities are not necessarily related via
an EOS. Although the FST models with $l=r_{\rm top}-r+\alpha^*H_{p,{\rm top}}$
are found to give the frequencies closest to the solar observed ones, a more
appropriate treatment of the top part of the 1D convective envelope is
necessary.
[2]
oai:arXiv.org:1809.08904 [pdf] - 1791054
Influence of metallicity on the near-surface effect on oscillation
frequencies
Submitted: 2018-09-24, last modified: 2018-09-25
The CoRoT and Kepler missions have provided high-quality measurements of the
frequency spectra of solar-like pulsators, enabling us to probe stellar
interiors with a very high degree of accuracy by comparing the observed and
modeled frequencies. However, the frequencies computed with 1D models suffer
from systematic errors related to the poor modeling of the uppermost layers of
stars. These biases are what is commonly named the near surface effect. The
dominant effect is related to the turbulent pressure that modifies the
hydrostatic equilibrium and thus the frequencies. This has already been
investigated using grids of 3D RMHD simulations, which also were used to
constrain the parameters of the empirical correction models. However, the
effect of metallicity has not been considered so far. We study the impact of
metallicity on the surface effect across the HR diagram, and provide a method
for accounting for it when using the empirical correction models. We computed a
grid of patched 1D stellar models with the stellar evolution code CESTAM in
which poorly modeled surface layers have been replaced by averaged
stratification computed with the 3D RMHD code CO5BOLD. We found that
metallicity has a strong impact on the surface effect: keeping T_eff and log g
constant, the frequency residuals can vary by up to a factor two. Therefore,
the influence of metallicity cannot be neglected. We found that a correct way
of accounting for it is to consider the surface Rosseland mean opacity. It
allowed us to give a physically-grounded justification as well as a scaling
relation for the frequency differences at nu_max as a function of T_eff, log g
and kappa. Finally, we provide prescriptions for the fitting parameters of the
correction functions. We show that the impact of metallicity through the
Rosseland mean opacity must be taken into account when studying and correcting
the surface effect.
[3]
oai:arXiv.org:1808.08391 [pdf] - 1743894
Mean density inversions for red giants and red clump stars
Buldgen, Gaël;
Rendle, B.;
Sonoi, T.;
Davies, G. R.;
Miglio, A.;
Salmon, S. J. A. J.;
Reese, D. R.;
Bossini, D.;
Eggenberger, P.;
Noels, A.;
Scuflaire, R.
Submitted: 2018-08-25
Since the CoRoT and Kepler missions, the availability of high quality seismic
spectra for red giants has made them the standard clocks and rulers for
Galactic Archeology. With the expected excellent data from the TESS and PLATO
missions, red giants will again play a key role in Galactic studies and stellar
physics, thanks to the precise masses and radii determined by asteroseismology.
The determination of these quantities is often based on so-called scaling laws,
which have been used extensively for main-sequence stars. We show how the SOLA
inversion technique can provide robust determinations of the mean density of
red giants within 1 per cent of the real value, using only radial oscillations.
Combined with radii determinations from Gaia of around 2 per cent precision,
this approach provides robust, less model-dependent masses with an error lower
than 10 per cent. It will improve age determinations, helping to accurately
dissect the Galactic structure and history. We present results on artificial
data of standard models, models including an extended atmosphere from averaged
3D simulations and non-adiabatic frequency calculations to test surface
effects, and on eclipsing binaries. We show that the inversions provide very
robust mean density estimates, using at best seismic information. However, we
also show that a distinction between red-giant branch and red-clump stars is
required to determine a reliable estimate of the mean density. The stability of
the inversion enables an implementation in automated pipelines, making it
suitable for large samples of stars.
[4]
oai:arXiv.org:1701.07244 [pdf] - 1554003
On the computation of eigenfrequencies for equilibrium models including
turbulent pressure
Submitted: 2017-01-25
The space-borne missions have provided a wealth of highly accurate data.
However, our inability to properly model the upper-most region of solar-like
stars prevents us from making the best of these observations. This problem is
called "surface effect" and a key ingredient to solve it is turbulent pressure
for the computation of both the equilibrium models and the oscillations. While
3D hydrodynamic simulations help to include properly the turbulent pressure in
the equilibrium models, the way this surface effect is included in the
computation of stellar oscillations is still subject to uncertainties. We aim
at determining how to properly include the effect of turbulent pressure and its
Lagrangian perturbation in the adiabatic computation of the oscillations. We
also discuss the validity of the gas-gamma model (GGM) and reduced gamma model
(RGM) approximations, which have been used to compute adiabatic oscillations of
equilibrium models including turbulent pressure. We use a patched model of the
Sun with an inner part constructed by a 1D stellar evolution code (CESTAM) and
an outer part by the 3D hydrodynamical code (CO$^5$BOLD). Then, the adiabatic
oscillations are computed using the ADIPLS code for the GGM and RGM and with
the MAD code imposing the adiabatic condition on an existing time-dependent
convection (TDC) formalism. We show that the computation of the oscillations
using the TDC formalism in the adiabatic limit improves significantly the
agreement with the observed frequencies compared to the GGM and RGM. Of the
components of the turbulent pressure perturbation, the perturbation of the
density and advection term is found to contribute most to the frequency shift.
We propose a formalism to evaluate the frequency shift due to the inclusion of
the term with the turbulent pressure perturbation in the variational principle
in order to extrapolate our result to other stars.
[5]
oai:arXiv.org:1510.01151 [pdf] - 1358893
Stellar oscillations. II The non-adiabatic case
Submitted: 2015-10-05
A leap forward has been performed due to the space-borne missions, MOST,
CoRoT and Kepler. They provided a wealth of observational data, and more
precisely oscillation spectra, which have been (and are still) exploited to
infer the internal structure of stars. While an adiabatic approach is often
sufficient to get information on the stellar equilibrium structures it is not
sufficient to get a full understanding of the physics of the oscillation.
Indeed, it does not permit one to answer some fundamental questions about the
oscillations, such as: What are the physical mechanisms responsible for the
pulsations inside stars? What determines the amplitudes? To what extent the
adiabatic approximation is valid? All these questions can only be addressed by
considering the energy exchanges between the oscillations and the surrounding
medium.
This lecture therefore aims at considering the energetical aspects of stellar
pulsations with particular emphasis on the driving and damping mechanisms. To
this end, the full non-adiabatic equations are introduced and thoroughly
discussed. Two types of pulsation are distinguished, namely the self-excited
oscillations that result from an instability and the solar-like oscillations
that result from a balance between driving and damping by turbulent convection.
For each type, the main physical principles are presented and illustrated using
recent observations obtained with the ultra-high precision photometry
space-borne missions (MOST, CoRoT and Kepler). Finally, we consider in detail
the physics of scaling relations, which relates the seismic global indices with
the global stellar parameters and gave birth to the development of statistical
(or ensemble) asteroseismology. Indeed, several of these relations rely on the
same cause: the physics of non-adiabatic oscillations.
[6]
oai:arXiv.org:1510.00300 [pdf] - 1304191
Surface-effect corrections for solar-like oscillations using 3D
hydrodynamical simulations
Submitted: 2015-10-01
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 $\nu_{\rm max}$.
[7]
oai:arXiv.org:1505.05452 [pdf] - 1232230
Angular momentum redistribution by mixed modes in evolved low-mass
stars. II. Spin-down of the core of red giants induced by mixed modes
Submitted: 2015-05-20
The detection of mixed modes in subgiants and red giants by the CoRoT and
\emph{Kepler} space-borne missions allows us to investigate the internal
structure of evolved low-mass stars. In particular, the measurement of the mean
core rotation rate as a function of the evolution places stringent constraints
on the physical mechanisms responsible for the angular momentum redistribution
in stars. It showed that the current stellar evolution codes including the
modelling of rotation fail to reproduce the observations. An additional
physical process that efficiently extracts angular momentum from the core is
thus necessary.
Our aim is to assess the ability of mixed modes to do this. To this end, we
developed a formalism that provides a modelling of the wave fluxes in both the
mean angular momentum and the mean energy equations in a companion paper. In
this article, mode amplitudes are modelled based on recent asteroseismic
observations, and a quantitative estimate of the angular momentum transfer is
obtained. This is performed for a benchmark model of 1.3 $M_{\odot}$ at three
evolutionary stages, representative of the evolved pulsating stars observed by
CoRoT and Kepler.
We show that mixed modes extract angular momentum from the innermost regions
of subgiants and red giants. However, this transport of angular momentum from
the core is unlikely to counterbalance the effect of the core contraction in
subgiants and early red giants. In contrast, for more evolved red giants, mixed
modes are found efficient enough to balance and exceed the effect of the core
contraction, in particular in the hydrogen-burning shell. Our results thus
indicate that mixed modes are a promising candidate to explain the observed
spin-down of the core of evolved red giants, but that an other mechanism is to
be invoked for subgiants and early red giants.
[8]
oai:arXiv.org:1505.05447 [pdf] - 1232229
Angular momentum redistribution by mixed modes in evolved low-mass
stars. I. Theoretical formalism
Submitted: 2015-05-20
Seismic observations by the space-borne mission \emph{Kepler} have shown that
the core of red giant stars slows down while evolving, requiring an efficient
physical mechanism to extract angular momentum from the inner layers. Current
stellar evolution codes fail to reproduce the observed rotation rates by
several orders of magnitude, and predict a drastic spin-up of red giant cores
instead. New efficient mechanisms of angular momentum transport are thus
required.
In this framework, our aim is to investigate the possibility that mixed modes
extract angular momentum from the inner radiative regions of evolved low-mass
stars. To this end, we consider the Transformed Eulerian Mean (TEM) formalism,
introduced by Andrews \& McIntyre (1978), that allows us to consider the
combined effect of both the wave momentum flux in the mean angular momentum
equation and the wave heat flux in the mean entropy equation as well as their
interplay with the meridional circulation.
In radiative layers of evolved low-mass stars, the quasi-adiabatic
approximation, the limit of slow rotation, and the asymptotic regime can be
applied for mixed modes and enable us to establish a prescription for the wave
fluxes in the mean equations. The formalism is finally applied to a $1.3
M_\odot$ benchmark model, representative of observed CoRoT and \emph{Kepler}
oscillating evolved stars.
We show that the influence of the wave heat flux on the mean angular momentum
is not negligible and that the overall effect of mixed modes is to extract
angular momentum from the innermost region of the star. A quantitative and
accurate estimate requires realistic values of mode amplitudes. This is
provided in a companion paper.
[9]
oai:arXiv.org:1404.5264 [pdf] - 1209071
Stability Analysis of Strange-Modes in Hot Massive Stars with
Time-Dependent Convection
Submitted: 2014-04-21
We carry out a nonadiabatic analysis of strange-modes in hot massive stars
with time-dependent convection (TDC). In envelopes of such stars, convective
luminosity is not so dominant as that in envelopes of stars in the redder side
of the classical instability strip. Around the Fe opacity bump, however,
convection non-negligibly contributes to energy transfer. Indeed, instability
of modes excited at the Fe bump is likely to be suppressed with TDC compared
with the case of adopting the frozen-in convection approximation. But we make
sure that unstable strange-modes certainly appear in hot massive stars even by
taking into account TDC. We also examine properties of the strange-mode
instability, which is related to destabilization of strange-modes without
adiabatic counterparts. In this type of instability, the phase lag between
density and pressure varies from 0 to $180^{\circ}$ in an excitation zone
unlike the case of the $\kappa$-mechanism. In addition, we confirm by comparing
models with $Z=0$ and $Z=0.02$ that dominance of radiation pressure is
important for this type of instability.
[10]
oai:arXiv.org:1311.2819 [pdf] - 1180657
Analysis of strange-mode instability with time-dependent convection in
hot massive stars
Submitted: 2013-11-12
We carry out nonadiabatic analysis of strange-modes in hot massive stars with
time-dependent convection (TDC) for the first time. Although convective
luminosity in envelopes of hot massive stars is not as dominative as in stars
near the red edge of the classical Cepheid instability strip in the
Hertzsprung-Russell (H-R) diagram, we have found that the strange-mode
instability can be affected by the treatment of convection. However, existence
of the instability around and over the Humphreys-Davidson (H-D) limit is
independent of the treatment. This implies that the strange-mode instability
could be responsible for the lack of observed stars over the H-D limit
regardless of uncertainties on convection theories.
[11]
oai:arXiv.org:1202.6432 [pdf] - 1116972
Dipole low-order g mode instability of metal-poor low-mass main-sequence
stars due to the $\varepsilon$-mechanism
Submitted: 2012-02-28
We analyzed vibrational stability of metal-poor low-mass main-sequence stars
due to the $\varepsilon$-mechanism. Since outer convection zones of the
metal-poor stars is limited only to the very outer layers, the uncertainty in
treatment of convection does not affect the result significantly. We found that
the dipole g$_1$- and g$_2$-modes certainly become unstable due to the
$\varepsilon$-mechanism for $Z\la 6\times 10^{-4}$. Besides that, we found that
as the metallicity decreases the mass range of the $\varepsilon$-mechanism
instability extends toward higher mass.
[12]
oai:arXiv.org:1202.5121 [pdf] - 1116845
Vibrational Instability of Metal-Poor Low-Mass Main-Sequence Stars
Submitted: 2012-02-23
We find that low-degree low-order g-modes become unstable in metal-poor
low-mass stars due to the $\varepsilon$-mechanism of the pp-chain. Since the
outer convection zone of these stars is limited only to the very outer layers,
the uncertainty in the treatment of convection does not affect the result
significantly. The decrease in metallicity leads to decrease in opacity and
hence increase in luminosity of a star. This makes the star compact and results
in decrease in the density contrast, which is favorable to the
$\varepsilon$-mechanism instability. We find also instability for high order
g-modes of metal-poor low-mass stars by the convective blocking mechanism.
Since the effective temperature and the luminosity of metal-poor stars are
significantly higher than those of Pop I stars, the stars showing $\gamma$
Dor-type pulsation are substantially less massive than in the case of Pop I
stars. We demonstrate that those modes are unstable for about $1\,M_\odot$
stars in the metal-poor case.
[13]
oai:arXiv.org:1202.0980 [pdf] - 1116395
Fully Nonadiabatic Analysis of Vibrational Instability of Population III
Stars due to the $\varepsilon$-Mechanism
Submitted: 2012-02-05
A linear nonadiabatic analysis of the vibrational stability of population III
main-sequence stars was carried out. It was demonstrated that, in the case of
massive stars with $M \gtrsim 5\Mo$, helium burning (triple alpha reaction)
starts during the main-sequence stage and produces $^{12}{\rm C}$, leading to
the activation of a part of the CNO-cycle. It was found that, despite of that,
those stars with $M\lesssim 13\Mo$ become unstable against the dipole g$_1$-
and g$_2$-modes due to the $\varepsilon$-mechanism, during the early
evolutionary phase at which the pp-chain is still the dominant nuclear energy
source. The instability due to the $\varepsilon$-mechanism occurs against
g-modes having a large amplitude in the off-centered $^3$He accumulation shell
in the deep interior, and the growth time is much shorter than the evolutionary
timescale. This instability is therefore likely to induce mixing in the stellar
interior, and have a significant influence on the evolution of these stars.
[14]
oai:arXiv.org:1201.0884 [pdf] - 459632
Vibrational instability of Population III very massive main-sequence
stars due to the $\varepsilon$-mechanism
Submitted: 2012-01-04
Very massive stars are thought to be formed in the early Universe because of
a lack of cooling process by heavy elements, and might have been responsible
for the later evolution of the Universe. We had an interest in vibrational
stability of their evolution and carried out the linear nonadiabatic analysis
of radial and nonradial oscillations for population III very massive
main-sequence stars with $500-3000M_{\sun}$. We found that only the radial
fundamental mode becomes unstable due to the $\varepsilon$-mechanism for these
stars. The instability appears just after the CNO cycle is activated and the
nuclear energy generation rate becomes large enough to stop the
pre--main-sequence contraction, and continues during the early stage of the
core hydrogen burning. Besides, we roughly estimated amount of mass loss due to
the instability to know its significance.