A More Realistic Representation of Overshoot at the Base of the Solar Convective Envelope As Seen by Helioseismology

Total Page:16

File Type:pdf, Size:1020Kb

A More Realistic Representation of Overshoot at the Base of the Solar Convective Envelope As Seen by Helioseismology Mon. Not. R. Astron. Soc. 414, 1158–1174 (2011) doi:10.1111/j.1365-2966.2011.18460.x A more realistic representation of overshoot at the base of the solar convective envelope as seen by helioseismology J. Christensen-Dalsgaard,1,2 M. J. P. F. G. Monteiro,3,4 M. Rempel2 and M. J. Thompson2,5 1Danish AsteroSeismology Centre, and Department of Physics and Astronomy, Aarhus University, Ny Munkegade, DK-8000 Aarhus C, Denmark 2High Altitude Observatory, NCAR†, Boulder, CO 80307-3000, USA 3Centro de Astrof´ısica da Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 4Departamento de F´ısica e Astronomia, Faculdade de Cienciasˆ da Universidade do Porto, Portugal 5School of Mathematics & Statistics, University of Sheffield, Hounsfield Road, Sheffield S3 7RH Accepted 2011 February 1. Received 2011 January 31; in original form 2010 October 28 ABSTRACT The stratification near the base of the Sun’s convective envelope is governed by processes of convective overshooting and element diffusion, and the region is widely believed to play a key role in the solar dynamo. The stratification in that region gives rise to a characteristic signal in the frequencies of solar p modes, which has been used to determine the depth of the solar convection zone and to investigate the extent of convective overshoot. Previous helioseismic investigations have shown that the Sun’s spherically symmetric stratification in this region is smoother than that in a standard solar model without overshooting, and have ruled out simple models incorporating overshooting, which extend the region of adiabatic stratification and have a more-or-less abrupt transition to subadiabatic stratification at the edge of the overshoot region. In this paper we consider physically motivated models which have a smooth transition in stratification bridging the region from the lower convection zone to the radiative interior beneath. We find that such a model is in better agreement with the helioseismic data than a standard solar model. Key words: asteroseismology – convection – Sun: helioseismology – Sun: interior – stars: interiors. age estimates of stars, and so improved constraints on theories of 1 INTRODUCTION overshooting obtained from a study of how overshooting works An understanding of the overshoot region at the bottom of the Sun’s in the solar case can be important for understanding stars more convective envelope is important for a number of reasons. The widely. overshoot region approximately coincides with the solar tachocline, Helioseismology provides a means of probing directly the con- a region of rotational shear which is generally believed to play a key ditions inside the Sun, because the frequencies of resonant modes role in the solar dynamo: overshooting is likely to be important for set up by acoustic waves propagating in the solar interior depend in helping to store the magnetic flux below the convection zone during particular on the adiabatic sound speed c which is given by the solar cycle. Bulk motion in the overshoot region also affects the p k T 2 = 1 1 B thermal stratification and it may contribute to significant mixing of c . (1) ρ muμ chemical elements, for example transporting fragile elements such Here, is the logarithmic derivative of pressure p with respect to as lithium to hotter regions where they are destroyed more easily 1 density ρ at constant specific entropy, T is temperature, μ is the than in the convection zone. More generally, convective overshoot mean molecular weight, k is Boltzmann’s constant and m is the in stars (particularly those with convective cores) is likely to be B u atomic mass unit. Hence, the sound-speed gradient with respect to an important and as yet imperfectly understood process affecting depth depends on the temperature gradient, which itself depends on the mechanism by which heat is transported. The transition between E-mail: [email protected] (JCD); [email protected] (MJPFGM); fully radiative heat transport beneath the convection zone and con- [email protected] (MR); [email protected] (MJT) vective heat transport within the convection zone is manifested in †The National Center for Atmospheric Research is operated by the Uni- the temperature gradient and hence too in the sound-speed gradient. versity Corporation for Atmospheric Research under sponsorship of the If the transition in sound-speed gradient takes place over a distance National Science Foundation. that is small compared with the vertical wavelength of the acoustic C 2011 The Authors Monthly Notices of the Royal Astronomical Society C 2011 RAS Overshoot in the Sun 1159 waves near the base of the convection zone, then the transition ap- When the observed and model frequencies are compared, it is pears to the waves to be more-or-less sharp and this gives rise to an found that the amplitude of δωp in the Sun is comparable with or oscillatory signal in the mode frequencies ω: the form of the signal smaller than that in models without overshooting, implying that gives information about the location and nature of this ‘acoustic the amount of overshooting of the kind predicted by these mixing- glitch’ (Gough 2002a). length models is very small. Monteiro et al. (1994) and Christensen- Monteiro, Christensen-Dalsgaard & Thompson (1994) repre- Dalsgaard et al. (1995) concluded that the extent of any such over- sented the effect of the base of the convection zone in terms of shoot at the base of the convection zone was less than one tenth an additive contribution δωp to the frequencies, relative to those of of a pressure scaleheight. A similar limit was found by Basu et al. a corresponding model in which the transition had been smoothed (1994), while Roxburgh & Vorontsov (1994) obtained a somewhat out. For low-degree modes (see Monteiro, Christensen-Dalsgaard weaker limit. Basu & Antia (1994) noted that the composition gra- & Thompson 2000) a break in the first derivative of the sound speed dient caused by the inclusion of helium settling produced a sharper gives rise to a contribution of the form transition in the sound-speed gradient, even in models without over- shoot, and hence a larger oscillatory signal. From these analyses, δω = A(ω)cos(2ωτ¯ + 2φ ) , (2) p d 0 it would appear that the transition in sound-speed gradient at the while a break in the second derivative gives rise to a similar signal base of the solar convection zone is if anything smoother than in the but with a sine term instead of a cosine, and with a different fre- non-overshoot models. A caveat is that what we purport to measure quency dependence of A(ω). Here,τ ¯d is essentially the value of τ at in the above studies is the spherically symmetric component of the the location of the acoustic glitch, where structure: departures from sphericity, such as a latitudinal depen- dence to the shape of the base of the convection zone, could make R dr τ = (3) the transition appear smoother than it is locally. The helioseismic r c evidence, however, is that the location of the base of the convection is the acoustic depth beneath the surface, r being the corresponding zone is independent of latitude (Monteiro & Thompson 1998; Basu distance to the centre and R the surface radius of the Sun. Also, φ0 & Antia 2001). Changes of the base of the convection zone on time- is a phase introduced by the reflection of the mode at the turning scales shorter than the observation interval could also have a similar points, depending in particular on the near-surface structure. The effect by introducing a time-averaging effect on the mode frequen- function A(ω) is an amplitude which depends on the sharpness and cies that would mimic a smoother transition. The importance of this nature of the convection-zone base: the smoother the transition, the effect is difficult to estimate, however, and depends strongly on the smaller in general will be the amplitude. However, if moderate- 3D nature of convection at the base of the envelope. degree data are used, as in the case of Sun where we have accurate Overshoot has been addressed in the last two decades by a va- data for modes whose degree is above 3, the above expression needs riety of 2D and 3D numerical simulations (Roxburgh & Simmons to include additional terms, both in the amplitude (Monteiro et al. 1993; Hurlburt et al. 1994; Singh, Roxburgh & Chan 1995, 1998; 1994) and in the argument of the signal (Christensen-Dalsgaard, Saikia et al. 2000; Brummell, Clune & Toomre 2002; Rogers & Monteiro & Thompson 1995), to account for the first-order effect Glatzmaier 2005a,b). Whilst the non-local mixing-length models of the mode degree on the signal. have clearly predicted an adiabatic overshoot region of a sizeable Overshoot at the base of the solar convection has traditionally fraction of a pressure scaleheight and a rather sharp transition to the been modelled using non-local mixing-length theory (e.g., Zahn radiative zone beneath, the numerical simulations show a greater 1991). Such models mostly predict an overshoot region that is nearly variety of possible behaviours. The work by Brummell et al. (2002) adiabatically stratified; in terms of the logarithmic temperature gra- is currently one of the best resolved and most turbulent investiga- −6 dient ∇=dlnT/dlnp one finds that δ ≡∇−∇ad ∼−10 ,where tions: it shows strongly subadiabatic overshoot with very smooth ∇ad = (∂ ln T/∂ ln p)s is the adiabatic temperature gradient, s be- transition towards the radiative temperature gradient. Most of the ing specific entropy. The depth of the overshoot region is typically earlier 2D and 3D simulation were more in the laminar regime and between 0.2Hp and 0.4Hp,whereHp is the pressure scaleheight at found, depending on their parameters (mainly the stiffness of the the base of the convection zone, with a very steep transition towards subadiabatic layer), both nearly adiabatic overshoot and extended the radiative temperature gradient.
Recommended publications
  • Future Prospects for Helioseismology
    Future Prospects for Helioseismology Michael J. Thompson High Altitude Observatory National Center for Atmospheric Research Boulder, Colorado, USA Abstract Helioseismology has been spectacularly successful at revealing many aspects of the Sun s internal structure and dynamics. "hat then are the new challenges and opportunities for the $eld? &n this talk & #ill briefly review the current state of helioseismology, and #ill then discuss science opportunities and future directions for observations from the ground and space. Solar oscillations: sound waves in the Sun Resonant cavities are of Refraction of #aves provides lower di*ering e+tent boundary of resonant cavity Resonant modes the sun l=-,m,. l=-,m,/ l=-,m,- 1+amples of spherical harmonics Spherical harmonic degree l, order m Schematic of a resonant mode of the sun l=/,m,0 Helioseismology Measure mode fre3uencies and other properties 1igenfunctions 4 spherical harmonics 5re3uencies νnlm(t7 depend on conditions in solar interior determining #ave propagation &nversion provides maps such as of sound speed, density and rotation Spherical harmonics Ring-Analysis Helioseismology · Observe line-of9sight velocity · 8ata used in this talk are 68opplergrams7! from SOHO/MD& · Make plane-wave decomposition of patches 6tiles) of the solar disk y t + :omography · 2ake surface observations at end-points of ray paths of waves propagating through the solar interior · Cross9correlating the observational signal at the two end-points reveals the wave travel time between the two points · &nversions of those travel times then reveal the subsurface structure and dynamics 6and magnetic field) Temperature Sound-speed inferences Convection zone Build-up of helium* Gravitational settling 2 ∝ Γ µ Opacity c 1T / Helium ionization zone Base of convection zone (Sound speed is in units of hundreds of km/s) ;&ndirect inference from small fre3uency separations Solar Rotation and the Tachocline Solar rotation Solar tachocline Near-surface Tachocline shear layer Tachocline lies substantially<.
    [Show full text]
  • Arxiv:1908.06042V4 [Astro-Ph.SR] 21 Jan 2020 Which Are Born on the Interface Between the Tachocline and the Overshoot Layer, Are Developed
    Thermomagnetic Ettingshausen-Nernst effect in tachocline, magnetic reconnection phenomenon in lower layers, axion mechanism of solar luminosity variations, coronal heating problem solution and mechanism of asymmetric dark matter variations around black hole V.D. Rusov1,∗ M.V. Eingorn2, I.V. Sharph1, V.P. Smolyar1, M.E. Beglaryan3 1Department of Theoretical and Experimental Nuclear Physics, Odessa National Polytechnic University, Odessa, Ukraine 2CREST and NASA Research Centers, North Carolina Central University, Durham, North Carolina, U.S.A. 3Department of Computer Technology and Applied Mathematics, Kuban State University, Krasnodar, Russia Abstract It is shown that the holographic principle of quantum gravity (in the hologram of the Uni- verse, and therefore in our Galaxy, and of course on the Sun!), in which the conflict between the theory of gravitation and quantum mechanics disappears, gives rise to the Babcock-Leighton holographic mechanism. Unlike the solar dynamo models, it generates a strong toroidal mag- netic field by means of the thermomagnetic Ettingshausen-Nernst (EN) effect in the tachocline. Hence, it can be shown that with the help of the thermomagnetic EN effect, a simple estimate of the magnetic pressure of an ideal gas in the tachocline of e.g. the Sun can indirectly prove that by using the holographic principle of quantum gravity, the repulsive toroidal magnetic field Sun 7 Sun of the tachocline (Btacho = 4:1 · 10 G = −Bcore) precisely \neutralizes" the magnetic field in the Sun core, since the projections of the magnetic fields in the tachocline and the core have equal values but opposite directions. The basic problem is a generalized problem of the antidy- namo model of magnetic flux tubes (MFTs), where the nature of both holographic effects (the thermomagnetic EN effect and Babcock-Leighton holographic mechanism), including magnetic cycles, manifests itself in the modulation of asymmetric dark matter (ADM) and, consequently, the solar axion in the Sun interior.
    [Show full text]
  • Perspective Helioseismology
    Proc. Natl. Acad. Sci. USA Vol. 96, pp. 5356–5359, May 1999 Perspective Helioseismology: Probing the interior of a star P. Demarque* and D. B. Guenther† *Department of Astronomy, Yale University, New Haven, CT 06520-8101; and †Department of Astronomy and Physics, Saint Mary’s University, Halifax, Nova Scotia, Canada B3H 3C3 Helioseismology offers, for the first time, an opportunity to probe in detail the deep interior of a star (our Sun). The results will have a profound impact on our understanding not only of the solar interior, but also neutrino physics, stellar evolution theory, and stellar population studies in astrophysics. In 1962, Leighton, Noyes and Simon (1) discovered patches of in the radial eigenfunction. In general low-l modes penetrate the Sun’s surface moving up and down, with a velocity of the more deeply inside the Sun: that is, have deeper inner turning order of 15 cmzs21 (in a background noise of 330 mzs21!), with radii than higher l-valued p-modes. It is this property that gives periods near 5 minutes. Termed the ‘‘5-minute oscillation,’’ the p-modes remarkable diagnostic power for probing layers of motions were originally believed to be local in character and different depth in the solar interior. somehow related to turbulent convection in the solar atmo- The cut-away illustration of the solar interior (Fig. 2 Left) sphere. A few years later, Ulrich (2) and, independently, shows the region in which p-modes propagate. Linearized Leibacher and Stein (3) suggested that the phenomenon is theory predicts (2, 4) a characteristic pattern of the depen- global and that the observed oscillations are the manifestation dence of the eigenfrequencies on horizontal wavelength.
    [Show full text]
  • Enhanced Rotational Mixing in the Radiative Zones of Massive Stars
    MNRAS 000,1{21 (2017) Preprint 26 October 2018 Compiled using MNRAS LATEX style file v3.0 Enhanced Rotational Mixing in the Radiative Zones of Massive Stars Adam S. Jermyn,1? Christopher A. Tout1 and Shashikumar M. Chitre1;2 1Institute of Astronomy, University of Cambridge, Madingley Rd, Cambridge CB3 0HA, UK 2Centre for Basic Sciences, University of Mumbai, India Accepted XXX. Received YYY; in original form ZZZ ABSTRACT Convection in the cores of massive stars becomes anisotropic when they rotate. This anisotropy leads to a misalignment of the thermal gradient and the thermal flux, which in turn results in baroclinicity and circulation currents in the upper radiative zone. We show that this induces a much stronger meridional flow in the radiative zone than previously thought. This drives significantly enhanced mixing, though this mixing does not necessarily reach the surface. The extra mixing takes on a similar form to convec- tive overshooting, and is relatively insensitive to the rotation rate above a threshold, and may help explain the large overshoot distances inferred from observations. This has significant consequences for the evolution of these stars by enhancing core-envelope mixing. Key words: Stars - evolution; Stars - interiors; Stars - rotation 1 INTRODUCTION motions, which set a baroclinic boundary condition at the base of the radiative zone. This boundary condition is of- The problem of meridional mixing in the radiative zones of ten more significant than the locally-generated baroclinicity rotating stars has been studied for nearly a century (Edding- and so results in enhanced circulation. This scenario has not ton 1929; Maeder & Zahn 1998).
    [Show full text]
  • Solar Turbulent Dynamo Near Convective Overshoot Layer and Radiative Tachocline
    SOLAR TURBULENT DYNAMO NEAR CONVECTIVE OVERSHOOT LAYER AND RADIATIVE TACHOCLINE V. N. Kryvodubskyj Astronomical Observatory, National Taras Shevchenko University of Kyiv 3 Observatorna Str., 04053 Kyiv, Ukraine e-mail: [email protected] In order to extend the abilities of the αΩ-dynamo model to explain the observed regularities and anomalies of the solar magnetic activity, the magnetic quenching of the α-effect were included in the model, and newest helioseismically determined inner rotation of the Sun has been used. Allowance for the radial inhomogeneity of turbulent velocity in derivations of the helicity parameter resulted in a change of sign of the α-effect from positive to negative in the northern hemisphere near the bottom of the solar convection zone (SCZ). The change of sign is very important for direction of Parker’s dynamo-waves propagation and for parity of excited magnetic fields. The period of the dynamo-wave calculated with taking into account the magnetic α-quenching is about seven years, that agrees by order of magnitude with the observed mean duration of the sunspot cycles. Using the modern helioseismology data to define dynamo-parameters, we conclude that north– south asymmetry should exist in the meridional field. The calculated configuration of the net meridional field is likely to explain the magnetic anomaly of polar fields (the apparent magnetic “monopole”) observed near the maxima of solar cycles. INTRODUCTION One of the most important problems in the solar physics is to understand the origin of the Sun’s magnetic cycle. Currently, it is widely believed that the solar cycle results from a turbulent hydromagnetic dynamo.
    [Show full text]
  • ASTR 565 Stellar Structure and Evolution Fall 2019 New Mexico State University
    1 ASTR 565 Stellar Structure and Evolution Fall 2019 New Mexico State University Jason Jackiewicz Department of Astronomy December 5, 2019 2 Contents I Stellar Structure 7 1 Equilibrium and time scales 9 1.1 Energy equilibrium .......................................... 9 1.2 Local thermodynamic equilibrium .................................. 10 1.3 Are stars a one-fluid plasma? .................................... 10 1.4 Simple time scales of stars ...................................... 11 1.4.1 Dynamical timescale ..................................... 11 1.4.2 Thermal timescale ...................................... 12 1.4.3 Nuclear timescale ....................................... 12 1.5 Problems ............................................... 12 2 Nuclear interactions 15 2.1 Coulomb barrier ........................................... 15 2.2 Nuclear reaction rates ........................................ 15 2.3 Gamow peak ............................................. 17 2.4 Problems ............................................... 19 3 Energy release in nuclear reactions 21 3.1 Mass excess .............................................. 21 3.2 Binding energy ............................................ 22 3.3 Hydrogen burning .......................................... 23 3.3.1 PP-I chain .......................................... 24 3.3.2 PP-II and PP-III chains ................................... 25 3.3.3 CNO cycle .......................................... 26 3.4 Problems ............................................... 27 4 Distribution
    [Show full text]
  • Arxiv:1512.07298V1
    Dynamics of Turbulent Convection and Convective Overshoot in a Moderate Mass Star I. N. Kitiashvili1, A. G. Kosovichev2, N. N. Mansour1, A. A. Wray1 1NASA Ames Research Center, Moffett Field, Mountain View, CA 94035, USA 2New Jersey Institute of Technology, Newark, NJ 07102, USA ABSTRACT Continued progress in observational stellar astrophysics requires a deep under- standing of the underlying convection dynamics. We present results of realistic 3D radiative hydrodynamic simulations of the outer layers of a moderate mass star (1.47 M⊙), including the full convection zone, the overshoot region, and the top layers of the radiative zone. The simulation results show that the surface granulation has a broad range of scales, from 2 to 12 Mm, and that large granules are organized in well-defined clusters, consisting of several granules. Comparison of the mean structure profiles from 3D simulations with the corresponding 1D standard stellar model shows an increase of the stellar radius by ∼ 800 km, as well as significant changes in the thermodynamic structure and turbulent properties of the ionization zones. Convective downdrafts in the intergranular lanes between granulation clusters reach speeds of more than 20 km s−1, penetrate through the whole convection zone, hit the radiative zone, and form a 8 Mm thick overshoot layer. Contrary to semi-empirical overshooting models, our results show that the 3D dynamic overshoot region consists of two layers: a nearly adiabatic extension of the convection zone and a deeper layer of enhanced subadiabatic stratification. This layer is formed because of heating caused by the braking of the overshooting convective plumes.
    [Show full text]
  • Convective Hydration in the Tropical Tropopause Layer During the Stratoclim Aircraft Campaign
    1 2 Convective hydration in the tropical tropopause layer during the StratoClim aircraft campaign: 3 Pathway of an observed hydration patch 4 5 Keun-Ok Lee1, Thibaut Dauhut1, Jean-Pierre Chaboureau1, Sergey Khaykin2, Martina Krämer3 and 6 Christian Rolf3 7 1Laboratoire d’Aérologie, Université de Toulouse, CNRS, UPS, Toulouse, France 8 2 LATMOS-IPSL, Universié Versailles St-Quentin; Sorbonne Université, CNRS/INSU, Guyancourt, France 9 3Institute for Energy and Climate Research – Stratosphere (IEK-7), Forschungzentrum Jülich, Jülich, Germany 10 11 ABSTRACT 12 The source and pathway of the hydration patch in the TTL (Tropical Tropopause Layer) that was measured 13 during the StratoClim field campaign during the Asian summer monsoon in 2017, and its connection to 14 convective overshoots are investigated. During flight #7, two remarkable layers are measured in the TTL 15 namely, (1) moist layer (ML) with water vapour content of 4.8−5.7 ppmv in altitudes of 18−19 km altitudes in 16 the lower stratosphere, and (2) ice layer (IL) with ice content up to 1.9 eq. ppmv in altitudes of 17−18 km in 17 the upper troposphere around 06:30 UTC on 8 August to the south of Kathmandu (Nepal). A Meso-NH 18 convection-permitting simulation succeeds in reproducing the characteristics of ML and IL. Through analysis, 19 we show that ML and IL are generated by convective overshoots that occurred over the Sichuan basin about 20 1.5 day before. Overshooting clouds develop up to 19 km, hydrating the lower stratosphere of up to 20 km 21 with 6401 t of water vapour by a strong-to-moderate mixing of the updraughts with the stratospheric air.
    [Show full text]
  • Modeling the Dynamical Coupling of Solar Convection with the Radiative Interior
    The Astrophysical Journal, 742:79 (20pp), 2011 December 1 doi:10.1088/0004-637X/742/2/79 C 2011. The American Astronomical Society. All rights reserved. Printed in the U.S.A. MODELING THE DYNAMICAL COUPLING OF SOLAR CONVECTION WITH THE RADIATIVE INTERIOR Allan Sacha Brun1, Mark S. Miesch2, and Juri Toomre3 1 Laboratoire AIM Paris-Saclay, CEA/Irfu Universite´ Paris-Diderot CNRS/INSU, 91191 Gif-sur-Yvette, France 2 High Altitude Observatory, NCAR, Boulder, CO 80307-3000, USA 3 JILA & Department of Astrophysical and Planetary Sciences, University of Colorado, Boulder, CO 80309-0440, USA Received 2011 February 11; accepted 2011 August 27; published 2011 November 8 ABSTRACT The global dynamics of a rotating star like the Sun involves the coupling of a highly turbulent convective envelope overlying a seemingly benign radiative interior. We use the anelastic spherical harmonic code to develop a new class of three-dimensional models that nonlinearly couple the convective envelope to a deep stable radiative interior. The numerical simulation assumes a realistic solar stratification from r = 0.07 up to 0.97R (with R the solar radius), thus encompassing part of the nuclear core up through most of the convection zone. We find that a tachocline naturally establishes itself between the differentially rotating convective envelope and the solid body rotation of the interior, with a slow spreading that is here diffusively controlled. The rapid angular momentum redistribution in the convective envelope leads to a fast equator and slow poles, with a conical differential rotation achieved at mid-latitudes, much as has been deduced by helioseismology.
    [Show full text]
  • Arxiv:2009.06972V2 [Astro-Ph.SR] 1 Oct 2020
    MNRAS 000,1{19 (2020) Preprint 2 October 2020 Compiled using MNRAS LATEX style file v3.0 Asteroseismic Inference of Subgiant Evolutionary Parameters with Deep Learning Marc Hon1?, Earl P. Bellinger2;1;3, Saskia Hekker3;2, Dennis Stello1;2;4, and James S. Kuszlewicz3;2 1School of Physics, The University of New South Wales, Sydney NSW 2052, Australia 2Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, Denmark 3Max-Planck-Institut fur¨ Sonnensystemforschung, Justus-von-Liebig-Weg 3, 37077 G¨ottingen, Germany 4Sydney Institute for Astronomy (SIfA), School of Physics, University of Sydney, NSW 2006, Australia Accepted XXX. Received YYY; in original form ZZZ ABSTRACT With the observations of an unprecedented number of oscillating subgiant stars ex- pected from NASA's TESS mission, the asteroseismic characterization of subgiant stars will be a vital task for stellar population studies and for testing our theories of stellar evolution. To determine the fundamental properties of a large sample of subgiant stars efficiently, we developed a deep learning method that estimates distributions of fundamental parameters like age and mass over a wide range of input physics by learn- ing from a grid of stellar models varied in eight physical parameters. We applied our method to four Kepler subgiant stars and compare our results with previously deter- mined estimates. Our results show good agreement with previous estimates for three of them (KIC 11026764, KIC 10920273, KIC 11395018). With the ability to explore a vast range of stellar parameters, we determine that the remaining star, KIC 10005473, is likely to have an age 1 Gyr younger than its previously determined estimate.
    [Show full text]
  • Three-Dimensional Simulations of the Upper Radiation–Convection Transition Layer in Subgiant Stars   F
    Mon. Not. R. Astron. Soc. 347, 1208–1216 (2004) Three-dimensional simulations of the upper radiation–convection transition layer in subgiant stars F. J. Robinson,1 P. Demarque,1 L. H. Li,1 S. Sofia,1 Y.-C. Kim,2 K. L. Chan3 and D. B. Guenther4 1Astronomy Department, Yale University, Box 208101, New Haven, CT 06520-8101, USA 2Astronomy Department, Yonsei University, Seoul, Korea 3Hong Kong University of Science & Technology, Hong Kong, China 4Department of Astronomy and Physics, Saint Mary’s University, Halifax, Nova Scotia B3A 4R2, Canada Accepted 2003 October 13. Received 2003 October 8; in original form 2003 May 19 ABSTRACT This paper describes three-dimensional (3D) large eddy simulations of stellar surface con- vection using realistic model physics. The simulations include the present Sun, a subgiant of one solar mass and a lower-gravity subgiant, also of one solar mass. We examine the thermal structure (superadiabaticity) after modification by 3D turbulence, the overshoot of convec- tive motions into the radiative atmosphere and the range of convection cell sizes. Differences between models based on the mixing length theory (MLT) and the simulations are found to increase significantly in the more evolved stages as the surface gravity decreases. We find that the full width at half maximum (FWHM) of the turbulent vertical velocity correlation provides a good objective measure of the vertical size of the convective cells. Just below the convection surface, the FWHM is close to the mean vertical size of the granules and 2 × FWHM is close to the mean horizontal diameter of the granules. For the Sun, 2 × FWHM = 1200 km, a value close to the observed mean granule size.
    [Show full text]
  • Numerical Simulations of Penetration and Overshoot in the Sun
    Numerical Simulations of Penetration and Overshoot in the Sun Tamara M. Rogers Astronomy and Astrophysics Department, University of California, Santa Cruz, CA 95064 [email protected] Gary A. Glatzmaier Earth Sciences Department, University of California, Santa Cruz, CA 95064 C.A. Jones Department of Applied Math, University of Leeds, Leeds LS2 9JT, UK ABSTRACT We present numerical simulations of convective overshoot in a two- dimensional model of the solar equatorial plane. The model equations are solved in the anelastic approximation with enhanced thermal conductivity and viscosity for numerical stability. The simulated domain extends from 0.001 R to 0.93 R, spanning both convective and radiative regions. We show that convective pen- etration leads to a slightly extended, mildly subadiabatic temperature gradient beneath the convection zone, spanning approximately 0.05 Hp,belowwhichthere is a rapid transition to a strongly subadiabatic region. A slightly higher tempera- ture is maintained in the overshoot region by adiabatic heating from overshooting plumes. This enhanced temperature may partially account for the sound speed discrepancy between the standard solar model and helioseismology. Simulations conducted with tracer particles suggest that a fully mixed region exists down to at least 0.684 R. Subject headings: convection: overshoot,mixing 1. Introduction One of the main unsolved problems in all of stellar evolution theory is the treatment of convective-radiative boundaries. In the Sun, understanding convective overshoot is cru- cial in determining properties of the solar tachocline and the solar dynamo. Helioseismic –2– observations have shown that the differential rotation observed at the solar surface persists throughout the convection zone (Thompson et al.
    [Show full text]