Source author record

D. R. Reese

D. R. Reese appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2020arXiv

Oscillations of 2D ESTER models. I. The adiabatic case

Recent numerical and theoretical considerations have shown that low-degree acoustic modes in rapidly rotating stars follow an asymptotic formula and recent observations of pulsations in rapidly rotating delta Scuti stars seem to match these expectations. However, a key question is whether strong gradients or discontinuities can adversely affect this pattern to the point of hindering its identification. Other important questions are how rotational splittings are affected by the 2D rotation profiles expected from baroclinic effects and whether it is possible to probe the rotation profile using these splittings. Accordingly, we numerically calculate pulsation modes in continuous and discontinuous rapidly rotating models produced by the 2D ESTER (Evolution STEllaire en Rotation) code. This spectral multi-domain code self-consistently calculates the rotation profile based on baroclinic effects and allows us to introduce discontinuities without loss of numerical accuracy. Pulsations are calculated using an adiabatic version of the Two-dimensional Oscillation Program (TOP) code. The variational principle is used to confirm the high accuracy of the pulsation frequencies and to derive an integral formula that closely matches the generalised rotational splittings, except when modes are involved in avoided crossings. This potentially allows us to probe the the rotation profile using inverse theory. Acoustic glitch theory, applied along the island mode orbit deduced from ray dynamics, can correctly predict the periodicity of the glitch frequency pattern produced by a discontinuity or the Gamma1 dip related to the He II ionisation zone in some of the models. The asymptotic frequency pattern remains sufficiently well preserved to potentially allow its detection in observed stars.

preprint2019arXiv

A realistic two-dimensional model of Altair

Fast rotation is responsible for important changes in the structure and evolution of stars. Optical long baseline interferometry now permits the study of its effects on the stellar surface, mainly gravity darkening and flattening. We aim to determine the fundamental parameters of the fast-rotating star Altair, in particular its evolutionary stage, mass, and differential rotation, using state-of-the-art stellar interior and atmosphere models together with interferometric, spectroscopic, and asteroseismic observations. We use ESTER 2D stellar models to produce the relevant surface parameters needed to create intensity maps from atmosphere models. Interferometric and spectroscopic observables are computed from these intensity maps and several stellar parameters are then adjusted using the MCMC algorithm Emcee. We determined Altair's equatorial radius to be 2.008 +/- 0.006 Rsun, the position angle 301.1 +/- 0.3 degrees, the inclination 50.7 +/- 1.2 degrees, and the equatorial angular velocity 0.74 +/- 0.01 times the Keplerian angular velocity. This angular velocity leads to a flattening of 0.220 +/- 0.003. We also deduce from the spectroscopically derived vsini ~ 243 km/s, a true equatorial velocity of ~314 km/s corresponding to a rotation period of 7h46m (~3 c/d). The data also impose a strong correlation between mass, metallicity, hydrogen abundance, and core evolution. Thanks to asteroseismic data, we constrain the mass of Altair to 1.86 +/- 0.03 Msun and further deduce its metallicity Z = 0.019 and its core hydrogen mass fraction Xc = 0.71, assuming an initial solar hydrogen mass fraction X = 0.739. These values suggest that Altair is ~100 Myrs old. Finally, the 2D ESTER model also gives the internal differential rotation of Altair, showing that its core rotates approximately 50% faster than the envelope, while the surface differential rotation does not exceed 6%.

preprint2016arXiv

In-depth study of 16CygB using inversion techniques

The 16Cyg binary system hosts the solar-like Kepler targets with the most stringent observational constraints. Moreover, this system is particularly interesting since both stars are very similar in mass but the A component is orbited by a red dwarf, whereas the B component is orbited by a Jovian planet and thus could have formed a more complex planetary system. In our previous study, we showed that seismic inversions of integrated quantities could be used to constrain microscopic diffusion in the A component. In this study, we analyse the B component in the light of a more regularised inversion. We wish to analyse independently the B component of the 16Cyg binary system using the inversion of an indicator dedicated to analyse core conditions, denoted tu. Using this independent determination, we wish to analyse any differences between both stars due to the potential influence of planetary formation on stellar structure and/or their respective evolution. First, we recall the observational constraints for 16CygB and the method we used to generate reference stellar models. The inversion results were then used to analyse the differences between the A and B components. The tu indicator for 16CygB shows a disagreement with models including microscopic diffusion and sharing the chemical composition previously derived for 16CygA. Small changes in chemical composition are insufficient to solve the problem but that extra mixing can account for the differences seen between both stars. We use a parametric approach to analyse the impact of extra mixing in the form of turbulent diffusion on the behaviour of the tu values. We conclude on the necessity of further investigations using models with a physically motivated implementation of extra mixing including additional constraints to further improve the accuracy with which the fundamental parameters of this system are determined.

preprint2016arXiv

SpaceInn hare-and-hounds exercise: Estimation of stellar properties using space-based asteroseismic data

Context: Detailed oscillation spectra comprising individual frequencies for numerous solar-type stars and red giants are or will become available. These data can lead to a precise characterisation of stars. Aims: Our goal is to test and compare different methods for obtaining stellar properties from oscillation frequencies and spectroscopic constraints, in order to evaluate their accuracy and the reliability of the error bars. Methods: In the context of the SpaceInn network, we carried out a hare-and-hounds exercise in which one group produced "observed" oscillation spectra for 10 artificial solar-type stars, and various groups characterised these stars using either forward modelling or acoustic glitch signatures. Results: Results based on the forward modelling approach were accurate to 1.5 % (radius), 3.9 % (mass), 23 % (age), 1.5 % (surface gravity), and 1.8 % (mean density). For the two 1 Msun stellar targets, the accuracy on the age is better than 10 % thereby satisfying PLATO 2.0 requirements. The average accuracies for the acoustic radii of the base of the convection zone, the He II ionisation, and the Gamma_1 peak were 17 %, 2.4 %, and 1.9 %, respectively. Glitch fitting analysis seemed to be affected by aliasing problems for some of the targets. Conclusions: Forward modelling is the most accurate approach, but needs to be complemented by model-independent results from, e.g., glitch analysis. Furthermore, global optimisation algorithms provide more robust error bars.

preprint2015arXiv

A test of the asteroseismic numax scaling relation for solar-like oscillations in main-sequence and sub-giant stars

Large-scale analyses of stellar samples comprised of cool, solar-like oscillators now commonly utilize the so-called asteroseismic scaling relations to estimate fundamental stellar properties. In this paper we present a test of the scaling relation for the global asteroseismic parameter $ν_{\rm max}$, the frequency at which a solar-like oscillator presents its strongest observed pulsation amplitude. The classic relation assumes that this characteristic frequency scales with a particular combination of surface gravity and effective temperature that also describes the dependence of the cut-off frequency for acoustic waves in an isothermal atmosphere, i.e., $ν_{\rm max} \propto gT_{\rm eff}^{-1/2}$. We test how well the oscillations of cool main-sequence and sub-giant stars adhere to this relation, using a sample of asteroseismic targets observed by the NASA \emph{Kepler} Mission. Our results, which come from a grid-based analysis, rule out departures from the classic $gT_{\rm eff}^{-1/2}$ scaling dependence at the level of $\simeq 1.5\,\rm per cent$ over the full $\simeq 1560\,\rm K$ range in $T_{\rm eff}$ that we tested. There is some uncertainty over the absolute calibration of the scaling. However, any variation with $T_{\rm eff}$ is evidently small, with limits similar to those above.

preprint2015arXiv

Constraints on the structure of 16 Cyg A and 16 Cyg B using inversion techniques

Constraining mixing processes and chemical composition is a central problem in stellar physics as their impact on stellar age determinations leads to biases in our studies of stellar evolution, galactic history and exoplanetary systems. In two previous papers, we showed how seismic inversion techniques could offer strong constraints on such processes by pointing out weaknesses in theoretical models. We now apply our technique to the solar analogues 16CygA and 16CygB, being amongst the best targets in the Kepler field to test the diagnostic potential of seismic inversions. The combination of various seismic indicators helps to provide more constrained and accurate fundamendal parameters for these stars. We use the latest seismic, spectroscopic and interferometric observational constraints in the litterature for this system to determine reference models independently for both stars. We carry out seismic inversions of the acoustic radius, the mean density and a core conditions indicator. We note that a degeneracy exists for the reference models. Namely, changing the diffusion coefficient or the chemical composition within the observational values leads to 5% changes in mass, 3% changes in radius and up to 8% changes in age. We use acoustic radius and mean density inversions to improve our reference models then carry out inversions for a core conditions indicator. Thanks to its sensitivity to microscopic diffusion and chemical composition mismatches, we are able to reduce the mass dispersion to 2%, namely [0.96, 1.0] M_sun, the radius dispersion to 1%, namely [1.188, 1.200] R_sun and the age dispersion to 3%, namely [7.0, 7.4] Gy, for 16CygA. For 16CygB, we can check the consistency of the models but not reduce independently the age dispersion. Nonetheless, assuming consistency with the age of 16CygA helps to further constrain its mass and radius.

preprint2015arXiv

Inequalities on stellar rotational splittings derived from assumptions on the rotation profile

Context: A number of pulsating stars with rotational splittings have been observed thanks to the CoRoT and Kepler missions. This is particularly true of evolved (sub-giant and giant) stars, and has led various groups to investigate their rotation profiles via different methods. Aims: We would like to set up some criteria which will help us to know whether a decreasing rotation profile, or one which satisfies Rayleigh's stability criterion, is compatible with a set of observed rotational splittings for a given reference model. Methods: We derive inequalities on the rotational splittings using a reformulated version of the equation which relates the splittings to the rotation profile and kernels. Results: These inequalities are tested out on some simple examples. The first examples show how they are able to reveal when a rotation profile is increasing somewhere or inconsistent with Rayleigh's criterion in a main sequence star, depending on the profile and the $\ell$ values of the splittings. The next example illustrates how a slight mismatch between an observed evolved star and a reference model can lead to erroneous conclusions about the rotation profile. We also show how frequency differences between the star and the model, which should normally reveal this mismatch, can be masked by frequency corrections for near-surface effects.

preprint2015arXiv

Using seismic inversions to obtain an internal mixing processes indicator for main-sequence solar-like stars

Determining accurate and precise stellar ages is a major problem in astrophysics. These determinations are either obtained through empirical relations or model-dependent approaches. Currently, seismic modelling is one of the best ways of providing accurate ages. However, current methods are affected by simplifying assumptions concerning mixing processes. In this context, providing new structural indicators which are less model-dependent and more sensitive to such processes is crucial. We build a new indicator for core conditions on the main sequence, which should be more sensitive to structural differences and applicable to older stars than the indicator t presented in a previous paper. We also wish to analyse the importance of the number and type of modes for the inversion, as well as the impact of various constraints and levels of accuracy in the forward modelling process that is used to obtain reference models for the inversion. First, we present a method to obtain new structural kernels and use them to build an indicator of central conditions in stars and test it for various effects including atomic diffusion, various initial helium abundances and metallicities, following the seismic inversion method presented in our previous paper. We then study its accuracy for 7 different pulsation spectra including those of 16CygA and 16CygB and analyse its dependence on the reference model by using different constraints and levels of accuracy for its selection We observe that the inversion of the new indicator using the SOLA method provides a good diagnostic for additional mixing processes in central regions of stars. Its sensitivity allows us to test for diffusive processes and chemical composition mismatch. We also observe that octupole modes can improve the accuracy of the results, as well as modes of low radial order.

preprint2014arXiv

Are the stars of a new class of variability detected in NGC~3766 fast rotating SPB stars?

A recent photometric survey in the NGC~3766 cluster led to the detection of stars presenting an unexpected variability. They lie in a region of the Hertzsprung-Russell (HR) diagram where no pulsation are theoretically expected, in between the $δ$ Scuti and slowly pulsating B (SPB) star instability domains. Their variability periods, between $\sim$0.1--0.7~d, are outside the expected domains of these well-known pulsators. The NCG~3766 cluster is known to host fast rotating stars. Rotation can significantly affect the pulsation properties of stars and alter their apparent luminosity through gravity darkening. Therefore we inspect if the new variable stars could correspond to fast rotating SPB stars. We carry out instability and visibility analysis of SPB pulsation modes within the frame of the traditional approximation. The effects of gravity darkening on typical SPB models are next studied. We find that at the red border of the SPB instability strip, prograde sectoral (PS) modes are preferentially excited, with periods shifted in the 0.2--0.5~d range due to the Coriolis effect. These modes are best seen when the star is seen equator-on. For such inclinations, low-mass SPB models can appear fainter due to gravity darkening and as if they were located between the $δ$~Scuti and SPB instability strips.

preprint2014arXiv

Seismic constraints on the radial dependence of the internal rotation profiles of six Kepler subgiants and young red giants

Context : We still do not know which mechanisms are responsible for the transport of angular momentum inside stars. The recent detection of mixed modes that contain the signature of rotation in the spectra of Kepler subgiants and red giants gives us the opportunity to make progress on this issue. Aims: Our aim is to probe the radial dependance of the rotation profiles for a sample of Kepler targets. For this purpose, subgiants and early red giants are particularly interesting targets because their rotational splittings are more sensitive to the rotation outside the deeper core than is the case for their more evolved counterparts. Methods: We first extract the rotational splittings and frequencies of the modes for six young Kepler red giants. We then perform a seismic modeling of these stars using the evolutionary codes CESAM2k and ASTEC. By using the observed splittings and the rotational kernels of the optimal models, we perform inversions of the internal rotation profiles of the six stars. Results: We obtain estimates of the mean rotation rate in the core and in the convective envelope of these stars. We show that the rotation contrast between the core and the envelope increases during the subgiant branch. Our results also suggest that the core of subgiants spins up with time, contrary to the RGB stars whose core has been shown to spin down. For two of the stars, we show that a discontinuous rotation profile with a deep discontinuity reproduces the observed splittings significantly better than a smooth rotation profile. Interestingly, the depths that are found most probable for the discontinuities roughly coincide with the location of the H-burning shell, which separates the layers that contract from those that expand. These results will bring observational constraints to the scenarios of angular momentum transport in stars.

preprint2014arXiv

The PLATO 2.0 Mission

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

preprint2014arXiv

The puzzling new class of variable stars in NGC 3766 : old friend pulsators?

The recent variability survey of the NGC 3766 cluster revealed a considerable number of periodic variable stars in a region of the H-R diagram where no pulsation is expected. This region lies between the instability strips of the delta Scuti and SPB stars. Moreover the periods of the new phenomenon, P~0.1-0.7 d, do not allow to associate it a priori to either of these two types of pulsations. Stars in the NGC 3766 cluster are known as fast rotators with rotational velocities typically larger than half of their critical velocity. Rotation can affect both the geometrical properties and period domain of pulsations. It also alters the apparent stellar luminosity through gravity darkening, effect seldom taken considered in theoretical studies of the rotation-pulsation interaction. We explore if both of these effects are able to deliver a consistent interpretation for the observed properties of the "new variables" in NGC 3766: explaining their presence outside the known instability strips and their variability periods. We carry out an instability analysis of SPB models within the framework of the Traditional Approximation of Rotation and study the visibility of modes according to the angle of view and rotation. We also check how gravity darkening affects the effective temperature and luminosity of stellar models for different angles of view and rotation velocities. At the red (cold) border of the instability strip, prograde sectoral modes are preferentially excited and their visibilities are maximum when seen equator-on. Furthermore low-mass SPB models seen equator-on can appear in the gap between non-rotating SPB and delta Scuti stars due to gravity darkening. In that case, periods of these most visible modes are shifted to the 0.2-0.5 d range due to the effects of the Coriolis force. We hence suggest that the new variable stars observed in NGC 3766 are actually fast rotating SPB pulsators.

preprint2013arXiv

2D non-perturbative modeling of oscillations in rapidly rotating stars

We present and discuss results of a recently developped two dimensional non-perturbative method to compute accurate adiabatic oscillation modes of rapidly rotating stars . The 2D calculations fully take into account the centrifugal distorsion of the star while the non-perturbative method includes the full influence of the Coriolis acceleration. These characteristics allows us to compute oscillation modes of rapid rotators - from high order p-modes in $δ$Scuti stars, to low order p- and g-modes in $β$ Cephei or Be stars.

preprint2013arXiv

Stable higher order finite-difference schemes for stellar pulsation calculations

Context: Calculating stellar pulsations requires a sufficient accuracy to match the quality of the observations. Many current pulsation codes apply a second order finite-difference scheme, combined with Richardson extrapolation to reach fourth order accuracy on eigenfunctions. Although this is a simple and robust approach, a number of drawbacks exist thus making fourth order schemes desirable. A robust and simple finite-difference scheme, which can easily be implemented in either 1D or 2D stellar pulsation codes is therefore required. Aims: One of the difficulties in setting up higher order finite-difference schemes for stellar pulsations is the so-called mesh-drift instability. Current ways of dealing with this defect include introducing artificial viscosity or applying a staggered grids approach. However these remedies are not well-suited to eigenvalue problems, especially those involving non-dissipative systems, because they unduly change the spectrum of the operator, introduce supplementary free parameters, or lead to complications when applying boundary conditions. Methods: We propose here a new method, inspired from the staggered grids strategy, which removes this instability while bypassing the above difficulties. Furthermore, this approach lends itself to superconvergence, a process in which the accuracy of the finite differences is boosted by one order. Results: This new approach is shown to be accurate, flexible with respect to the underlying grid, and able to remove mesh-drift.

preprint2013arXiv

Study of KIC 8561221 observed by Kepler: an early red giant showing depressed dipolar modes

The continuous high-precision photometric observations provided by the CoRoT and Kepler space missions have allowed us to better understand the structure and dynamics of red giants using asteroseismic techniques. A small fraction of these stars shows dipole modes with unexpectedly low amplitudes. The reduction in amplitude is more pronounced for stars with higher frequency of maximum power. In this work we want to characterize KIC 8561221 in order to confirm that it is currently the least evolved star among this peculiar subset and to discuss several hypotheses that could help explain the reduction of the dipole mode amplitudes. We used Kepler short- and long-cadence data combined with spectroscopic observations to infer the stellar structure and dynamics of KIC 8561221. We then discussed different scenarios that could contribute to the reduction of the dipole amplitudes such as a fast rotating interior or the effect of a magnetic field on the properties of the modes. We also performed a detailed study of the inertia and damping of the modes. We have been able to characterize 37 oscillations modes, in particular, a few dipole modes above nu_max that exhibit nearly normal amplitudes. We have inferred a surface rotation period of around 91 days and uncovered the existence of a variation in the surface magnetic activity during the last 4 years. As expected, the internal regions of the star probed by the l = 2 and 3 modes spin 4 to 8 times faster than the surface. With our grid of standard models we are able to properly fit the observed frequencies. Our model calculation of mode inertia and damping give no explanation for the depressed dipole modes. A fast rotating core is also ruled out as a possible explanation. Finally, we do not have any observational evidence of the presence of a strong deep magnetic field inside the star.

preprint2012arXiv

Estimating stellar mean density through seismic inversions

Determining the mass of stars is crucial both to improving stellar evolution theory and to characterising exoplanetary systems. Asteroseismology offers a promising way to estimate stellar mean density. When combined with accurate radii determinations, such as is expected from GAIA, this yields accurate stellar masses. The main difficulty is finding the best way to extract the mean density from a set of observed frequencies. We seek to establish a new method for estimating stellar mean density, which combines the simplicity of a scaling law while providing the accuracy of an inversion technique. We provide a framework in which to construct and evaluate kernel-based linear inversions which yield directly the mean density of a star. We then describe three different inversion techniques (SOLA and two scaling laws) and apply them to the sun, several test cases and three stars. The SOLA approach and the scaling law based on the surface correcting technique described by Kjeldsen et al. (2008) yield comparable results which can reach an accuracy of 0.5 % and are better than scaling the large frequency separation. The reason for this is that the averaging kernels from the two first methods are comparable in quality and are better than what is obtained with the large frequency separation. It is also shown that scaling the large frequency separation is more sensitive to near-surface effects, but is much less affected by an incorrect mode identification. As a result, one can identify pulsation modes by looking for an l and n assignment which provides the best agreement between the results from the large frequency separation and those from one of the two other methods. Non-linear effects are also discussed as is the effects of mixed modes. In particular, it is shown that mixed modes bring little improvement as a result of their poorly adapted kernels.

preprint2012arXiv

Multisite spectroscopic seismic study of the beta Cep star V2052 Oph: inhibition of mixing by its magnetic field

We used extensive ground-based multisite and archival spectroscopy to derive observational constraints for a seismic modelling of the magnetic beta Cep star V2052 Ophiuchi. The line-profile variability is dominated by a radial mode (f_1=7.14846 d^{-1}) and by rotational modulation (P_rot=3.638833 d). Two non-radial low-amplitude modes (f_2=7.75603 d^{-1} and f_3=6.82308 d^{-1}) are also detected. The four periodicities that we found are the same as the ones discovered from a companion multisite photometric campaign (Handler et al. 2012) and known in the literature. Using the photometric constraints on the degrees l of the pulsation modes, we show that both f_2 and f_3 are prograde modes with (l,m)=(4,2) or (4,3). These results allowed us to deduce ranges for the mass (M \in [8.2,9.6] M_o) and central hydrogen abundance (X_c \in [0.25,0.32]) of V2052 Oph, to identify the radial orders n_1=1, n_2=-3 and n_3=-2, and to derive an equatorial rotation velocity v_eq \in [71,75] km s^{-1}. The model parameters are in full agreement with the effective temperature and surface gravity deduced from spectroscopy. Only models with no or mild core overshooting (alpha_ov \in [0,0.15] local pressure scale heights) can account for the observed properties. Such a low overshooting is opposite to our previous modelling results for the non-magnetic beta Cep star theta Oph having very similar parameters, except for a slower surface rotation rate. We discuss whether this result can be explained by the presence of a magnetic field in V2052 Oph that inhibits mixing in its interior.

preprint2012arXiv

Seismic evidence for a rapidly rotating core in a lower-giant-branch star observed with Kepler

Rotation is expected to have an important influence on the structure and the evolution of stars. However, the mechanisms of angular momentum transport in stars remain theoretically uncertain and very complex to take into account in stellar models. To achieve a better understanding of these processes, we desperately need observational constraints on the internal rotation of stars, which until very recently were restricted to the Sun. In this paper, we report the detection of mixed modes - i.e. modes that behave both as g modes in the core and as p modes in the envelope - in the spectrum of the early red giant KIC7341231, which was observed during one year with the Kepler spacecraft. By performing an analysis of the oscillation spectrum of the star, we show that its non-radial modes are clearly split by stellar rotation and we are able to determine precisely the rotational splittings of 18 modes. We then find a stellar model that reproduces very well the observed atmospheric and seismic properties of the star. We use this model to perform inversions of the internal rotation profile of the star, which enables us to show that the core of the star is rotating at least five times faster than the envelope. This will shed new light on the processes of transport of angular momentum in stars. In particular, this result can be used to place constraints on the angular momentum coupling between the core and the envelope of early red giants, which could help us discriminate between the theories that have been proposed over the last decades.

preprint2011arXiv

2D computations of g modes

We present complete 2D computations of g modes in distorted polytropic models of stars performed with the Two-dimensional Oscillation Program (TOP). We computed low-degree modes (l=1 modes with radial order n=-1...-14, and l=2,3 modes with n=-1...-5 and -16...-20) of a nonrotating model and followed them by slowly increasing the rotation rate up to 70 % of the Keplerian break-up velocity. We use these computations to determine the domain of validity of perturbative methods up to the 3rd order. We study the evolution of the regularities of the spectrum and show quantitative agreement with the traditional approximation for not too large values of the ratio of the rotation rate to the pulsation frequency. We also show the appearance of new types of modes, called "rosette" modes due to their spatial structure. Thanks to the ray theory for gravito-inertial waves that we developed, we can associate these modes with stable periodic rays.

preprint2011arXiv

On the effects of rotation on acoustic stellar pulsations: validity domains of perturbative methods and close frequency pairs

Pulsation frequencies of acoustic modes are calculated for realistic rotating stellar models using both a perturbative and a two-dimensional approach. A comparison between the two yields validity domains which are similar to those previously obtained in Reese, Lignieres and Rieutord (2006) for polytropic models. One can also construct validity domains based on polynomial fits to the frequencies from the two-dimensional approach, and these also turn out to be similar, thus further confirming the agreement between the perturbative and two-dimensional approach at low rotation rates. Furthermore, as was previously shown in Espinosa, Perez Hernandez and Roca Cortes (2004), adjacent frequencies in multiplets come close together, thus forming pairs. This phenomena, exclusive to two-dimensional calculations, is shown to be an unlikely explanation of the close frequency pairs observed in delta Scuti stars. A systematic search for all close frequency pairs in the calculated spectrum is also carried out. The number of close frequency pairs is shown to agree with what is expected based on a Poisson distribution, but does not match the number or distribution of close pairs in stars like FG Vir. Furthermore, a lack of close frequency pairs appears at low rotation rates, where frequency multiplets do not overlap. Delta Scuti stars currently reported as having close frequency pairs do not fall in this interval.

preprint2011arXiv

Regular modes in rotating stars

Despite more and more observational data, stellar acoustic oscillation modes are not well understood as soon as rotation cannot be treated perturbatively. In a way similar to semiclassical theory in quantum physics, we use acoustic ray dynamics to build an asymptotic theory for the subset of regular modes which are the easiest to observe and identify. Comparisons with 2D numerical simulations of oscillations in polytropic stars show that both the frequency and amplitude distributions of these modes can accurately be described by an asymptotic theory for almost all rotation rates. The spectra are mainly characterized by two quantum numbers; their extraction from observed spectra should enable one to obtain information about stellar interiors.

preprint2010arXiv

Gravity modes in rapidly rotating stars. Limits of perturbative methods

CoRoT and Kepler missions are now providing high-quality asteroseismic data for a large number of stars. Among intermediate-mass and massive stars, fast rotators are common objects. Taking the rotation effects into account is needed to correctly understand, identify, and interpret the observed oscillation frequencies of these stars. A classical approach is to consider the rotation as a perturbation. In this paper, we focus on gravity modes, such as those occurring in gamma Doradus, slowly pulsating B (SPB), or Be stars. We aim to define the suitability of perturbative methods. With the two-dimensional oscillation program (TOP), we performed complete computations of gravity modes -including the Coriolis force, the centrifugal distortion, and compressible effects- in 2-D distorted polytropic models of stars. We started with the modes l=1, n=1-14, and l=2-3, n=1-5,16-20 of a nonrotating star, and followed these modes by increasing the rotation rate up to 70% of the break-up rotation rate. We then derived perturbative coefficients and determined the domains of validity of the perturbative methods. Second-order perturbative methods are suited to computing low-order, low-degree mode frequencies up to rotation speeds ~100 km/s for typical gamma Dor stars or ~150 km/s for B stars. The domains of validity can be extended by a few tens of km/s thanks to the third-order terms. For higher order modes, the domains of validity are noticeably reduced. Moreover, perturbative methods are inefficient for modes with frequencies lower than the Coriolis frequency 2Omega. We interpret this failure as a consequence of a modification in the shape of the resonant cavity that is not taken into account in the perturbative approach.

preprint2010arXiv

On the interpretation of echelle diagrams for solar-like oscillations. Effect of centrifugal distortion

This work aims at determining the impact of slow to moderate rotation on the regular patterns often present in solar-like oscillation spectra. We focus on the well-known asteroseismic diagnostic echelle diagrams, examining how rotation may modify the estimates of the large and small spacings, as well as the identification of modes. We illustrate the work with a real case: the solar-like star $η$ Bootis. The modeling takes into account rotation effects on the equilibrium models through an effective gravity and on the oscillation frequencies through both perturbative and non-perturbative calculations. We compare the results of both type of calculations in the context of the regular spacings (like the small spacings and the scaled small spacings) and echelle diagrams. We show that for echelle diagrams the perturbative approach remains valid for rotational velocities up to 40-50 km/s. We show that for the rotational velocities measured in solar-like stars, theoretical oscillation frequencies must be corrected up to the second-order in terms of rotation rate, including near degeneracy effects. For rotational velocities of about 16 km/S and higher, diagnostics on large spacings and on modal identification through echelle diagrams can be significantly altered by the presence of the $m\neq0$ components of the rotationally split modes. We found these effects to be detectable in the observed frequency range. Analysis of the effects of rotation on small spacings and scaled small spacings reveals that these can be of the order of, or even larger than surface effects, typically turbulence, microscopic diffusion, etc. Furthermore, we show that scaled spacings are significantly affected by stellar distortion even for small stellar rotational velocities (from 10-15 km/s) and therefore some care must be taken when using them as indicators for probing deep stellar interiors.

preprint2009arXiv

Mode identification in rapidly rotating stars

Context: Recent calculations of pulsation modes in rapidly rotating polytropic models and models based on the Self-Consistent Field method (MacGregor et al. 2007) have shown that the frequency spectrum of low degree pulsation modes can be described by an empirical formula similar to Tassoul's asymptotic formula (Tassoul 1980), provided that the underlying rotation profile is not too differential (Lignieres & Georgeot 2008, Reese et al. 2009). Aims: Given the simplicity of this asymptotic formula, we investigate whether it can provide a means by which to identify pulsation modes in rapidly rotating stars. Methods: We develop a new mode identification scheme which consists in scanning a multidimensional parameter space for the formula coefficients which yield the best-fitting asymptotic spectra. This mode identification scheme is then tested on artificial spectra based on the asymptotic formula, on random frequencies and on spectra based on full numerical eigenmode calculations for which the mode identification is known beforehand. We also investigate the effects of adding random frequencies to mimic the effects of chaotic modes which are also expected to show up in such stars (Lignieres & Georgeot 2008).

preprint2009arXiv

Pulsation modes in rapidly rotating stellar models based on the Self-Consistent Field method

Context: New observational means such as the space missions CoRoT and Kepler and ground-based networks are and will be collecting stellar pulsation data with unprecedented accuracy. A significant fraction of the stars in which pulsations are observed are rotating rapidly. Aims: Our aim is to characterise pulsation modes in rapidly rotating stellar models so as to be able to interpret asteroseismic data from such stars. Methods: The pulsation code developed in Lignières et al. (2006) and Reese et al. (2006) is applied to stellar models based on the self-consistent field (SCF) method (Jackson et al. 2004, 2005, MacGregor et al. 2007). Results: Pulsation modes in SCF models follow a similar behaviour to those in uniformly rotating polytropic models, provided that the rotation profile is not too differential. Pulsation modes fall into different categories, the three main ones being island, chaotic, and whispering gallery modes, which are rotating counterparts to modes with low, medium, and high l-|m| values, respectively. The frequencies of the island modes follow an asymptotic pattern quite similar to what was found for polytropic models. Extending this asymptotic formula to higher azimuthal orders reveals more subtle behaviour as a function of m and provides a first estimate of the average advection of pulsation modes by rotation. Further calculations based on a variational principle confirm this estimate and provide rotation kernels that could be used in inversion methods. When the rotation profile becomes highly differential, it becomes more and more difficult to find island and whispering gallery modes at low azimuthal orders. At high azimuthal orders, whispering gallery modes, and in some cases island modes, reappear.