arxiv: v1 [astro-ph.ep] 11 Jun 2016
|
|
- Anabel Robinson
- 6 years ago
- Views:
Transcription
1 A Simple Analytical Model for Rocky Planet Interior Li Zeng 1,a and Stein B. Jacobsen 1,b 1 Department of Earth and Planetary Sciences, Harvard University, Cambridge, MA arxiv: v1 [astro-ph.ep] 11 Jun 2016 a astrozeng@gmail.com b jacobsen@neodymium.harvard.edu ABSTRACT Based on the PREM (Dziewonski & Anderson 1981), the internal gravity of the Earth increases approximately linearly in the core and stays approximately constant in the mantle. Tested with numerical calculations in Zeng et al. (2016), this behavior of internal gravity can be generalized to any two-layer (iron-silicate) rocky exoplanet with core mass fraction (CMF) in between and mass in between 0.1 and 30 M. Based on this fact, (1) CMF CRF 2 (core radius fractionsquared), (2)P c gs 2 (interiorpressurescalesassurfacegravitysquared), (3) energy released of core formation is 1 the total gravitational energy, (4) effective heat capacity of the mantle ( 10 ) Mp M J/K, (5) moment of inertia 1 M 3 p Rp 2, and some other results can be derived. These results, though approximate, are handy to use owing to their simplicity and lucidity. They lead us deeper into the understanding of the interior structures of those planets. Subject headings: Extrasolar planets, interior structure, core mass fraction, core radius fraction, internal pressure, internal gravity, energy of differentiation, moment of inertia 1. Introduction Numerical calculations using realistic Equation of State(EOS)(c.f. Manipulate Planet at (Zeng et al. 2016) indicate that for a two-layer (core+mantle) rocky planet of M and 0.2 CMF 0.35, the internal gravity profile is approximately (Fig. 1): In the core, the gravity g increases linearly with radius from 0 at ( the) center to g s r (surface value) at the core-mantle boundary (CMB): g core (r) = g s R core r
2 2 In the mantle, g stays constant: g mantle (r) = g s In this approximation, core density is constant, then followed by a density drop of 2 3 at the CMB, and mantle density decease as 1/r in radius. From it, some important relations regarding rocky planet interior can be derived. The two first-order differential equations (Zeng & Sasselov 2013; Zeng & Seager 2008; Seager et al. 2007) governing rocky planet interior are: (1) hydrostatic equilibrium (force balance) equation: (2) mass conservation equation: dp dr = Gmρ r 2 = g ρ (1) dm dr = 4πr2 ρ (2) Eq. 1 and Eq. 2 can be combined to give a relation between the internal pressure P and the mass m (mass contained within radius r, now used as the independent variable instead): dp dm = Gm 4πr = 1 4 4πG g2 m = 1 dln(m) g2 4πG dm (3) Integrating Eq. 3, we get (ln(m) stands for natural logarithmic of m): dp = 1 4πG g 2 dln(m) (4) Therefore, the internal pressure of a planet is on the order of: P g2 4πG (5) where g 2 is an average of g 2. Here let s define planet mass as M p, planet radius as R p, and planet mean density as ρ p. Surface gravity g s and characteristic interior pressure P c are: g s GM p R 2 p (6)
3 3 P c g s 2 4πG = GM2 p 4πR 4 p (7) Later on, we will show that P c also approximates the CMB pressure P CMB. If g s is in S.I. unit (m/s 2 ) and P c in GPa, the following approximate relation holds: P c g 2 s (8) For example, Earth s surface gravity is 10 m/s 2, so 10 2 = 100 GPa is close to the actual CMB pressure of 136 GPa. 2. Internal Density Profile { Based on the assumption of gravity profile, the density profile can be easily derived: ρcore (r) = 3gs 1 4πGR core = ρ p = constant, so m CRF core(r) r 3 ρ mantle (r) = gs 1, so m 2πG r r mantle(r) r 2 Fig. 2 compares the PREM density profile (Dziewonski & Anderson 1981) to the approximation. Notice that CRF stands for core radius fraction, g is continuous in r, and ρ may have discontinuities. The 1 -dependence approximates the compression of mantle material towards depth, while the relatively small core (CMF 0.35) allows core density to be r approximated as constant. Anywhere in the mantle, In particular, at the core-mantle boundary (CMB) : ( ) 2 m r = (9) M p R p CMF M core M p = ( Rcore R p ) 2 = CRF 2 (10) In reality: CMF CRF 2 (11) Forexample, applyingeq.11toearthandkepler-93b(dressing et al.2015;ballard et al. 2014):
4 4 Fig. 1. gravity (black) and density (red) profile of the approximation. Fig. 2. black: PREM density profile. red: approximation.
5 5 For Earth, CMF = 0.325, and CRF = 3480km 6371km = CRF2 = = 0.298, which is about 9% smaller than For Kepler-93b, calculated from the model in (Zeng et al. 2016), CMF K93b = 0.278, and CRF K93b = CRF 2 K93b = 0.243, which is about 12% smaller than Inreality, CRF 2 tendstounderestimatecmfby 10%. Nevertheless, itprovidesaquick { CMF 0.3 way to estimate CMF from CRF or vice versa. For Earth-like rocky planets: CRF 0.5 Eq. 11 can be generalized to a rocky planet which has volatile envelope, as long as the formula is only applied to the solid portion of that planet. 3. Internal Pressure Profile 3.1. Pressure in the Mantle Since gravity in mantle is constant, Eq 4 can be integrated straightforward to get: P mantle (m) = g2 s 4πG ln ( ) Mp = P c ln m ( ) Mp = 2P c ln m ( Rp r ) (12) From Eq. 12, at the bottom of the mantle, the CMB pressure (P CMB ) is: ( ) ( ) 1 P CMB = P c ln = g2 s 1 CMF 4πG CMF (13) Considering its applicability to 0.2 CMF 0.35, then P CMB ( )P c. And P CMB determines the state of the mantle material in contact with the core. Surface gravity g s can be measured independent of the stellar parameters (Southworth et al. 2007), thus, for rocky planets it is possible to estimate P CMB even without accurately knowing the mass and radius in some cases: g s = 2π (1 e 2 ) 1/2 A RV P orb (R/a) 2 Sin[i] (14)
6 6 where semi-amplitude A RV andorbital eccentricity ecan beconstrained fromthe radialvelocity curve. R/a is the radius over semi-major axis ratio which could be constrained directly from the transit light curve. Orbital period P orb can be constrained from both Range of Applicability At any pressure P, the bulk modulus K P and K = ln(p). In the mantle: ln(ρ) P ln(ρ) K mantle = P mantle ln(ρ mantle ) = g2 s dln(m) 4πG dln(r) = 2P c (15) Thus, in the approximation, the mantle bulk modulus K mantle is treated as a constant equal to twice the characteristic internal pressure P c. For Earth, P c, = g 2 = 115 GPa, so in the approximation, K 4πG mantle, = 230 GPa. ComparingittotheisentropicbulkmodulusK s ofearthmantleaccordingtoprem(dziewonski & Anderson K,LID = 130 GPa 1981): K,670km = GPa, they are on the same order of magnitude. And K,D = 640 GPa this approximation is further justified by invoking the Birch-Murnaghan 2nd-order (BM2) EOS (Birch 1947, 1952), which when fitted to PREM gives K 0 200GPa for both core and mantle: P = 3 2 K 0 [ ( ρ ρ 0 )7 ( ] )5 3 ρ 3 ρ 0 (16) The bulk modulus in BM2 is obtained by the differentiating Eq. 16: K [ P ln(ρ/ρ 0 ) = 3 ( )7 ( 2 K 7 ρ ρ ρ 0 3 ρ 0 ] )5 3 (17) { when P K0,K K 0 Eq. 17 suggests: when P K 0,K/P 7 2. Since P c is the characteristic pressure in the mantle, we have K 2P 2P c, which is the same as Eq. 15. Therefore, this 3 approximation holds for rocky planets from 1M up to higher mass as long as BM2 EOS holds for that mass.
7 Pressure in the Core Since ρ core = const in the approximation, from Eq. 1 we have: dp core (r) dr = G r 2 ( ) ( ) gs r 3 3gs = 3P c r (18) GR core 4πGR core Rcore 2 Integrating it gives the pressure-dependence on r as a parabolic function: ( r ) 2 (19) P core (r) = P P c R core Central pressure P 0 can be determined by connecting it to P CMB at CMB: P 0 3 ( ) 1 2 P c = P CMB = P c ln CMF (20) So, P 0 = P c [ ( )] ln P 0 3P c (21) CMF Thus, in our approximation, the pressure dependence on radius is piecewise: parabolic in core and logarithmic in mantle. For Earth, it can be closely matched to the pressure profile calculated from PREM shown in Fig Energy of Core Formation According to Virial Theorem (Haswell 2010), the total gravitational energy is: surface P E grav = 3 center ρ dm (22) With the analytic forms of P and ρ in our approximation, Eq. 22 can be integrated: E grav = GM2 p 3R p (2 15 CRF3 ) (23)
8 8 Fig. 3. Piecewise approximation to pressure profile inside Earth (generally applicable to two-layer rocky planets of 0.2 CMF 0.35): Black curve: pressure profile from PREM. Red dashed curve: parabolic function in r as an approximation for pressure in core. Light-green dashed curve: logarithmic function in r as an approximation for pressure in mantle.
9 9 Comparingittothegravitationalenergyofauniform-density spheree grav, uniform sphere = 3 GMp 2 5 R p, the difference of the two can be taken as the energy released during core formation (release of gravitational energy due to concentration of denser material towards the center): E diff = E grav, uniform sphere E grav = GM2 p R p ( ) 1 15 (1 CRF3 ) = GM2 p R p ( ) 1 15 (1 CMF3/2 ) (24) Since 0.2 CMF 0.35, the term CMF 3/2 is small enough to be dropped to give: E diff 1 15 GM 2 p R p 1 10 E grav (25) Therefore, the energy released during core formation is 10% of the total gravitational energy of such an Earth-like rocky planet. For Earth, E diff, J. 5. Heat Capacity of the Rocky Mantle Since the temperatures inside the mantle of such a planet are likely above the Debye temperature, the heat capacity per mole of atoms can be approximated as 3R (gas constant R=8.314 J K 1 mol 1 ). The specific heat capacity (heat capacity per unit mass) is 3R/µ where µ is the average atomic weight of the composition, which for Mg-silicates (MgO, SiO 2, or any proportion of them combined, including MgSiO 3 or Mg 2 SiO 4 ) is 0.02 kg/mol. The specific thermal energy u th of the mantle material is thus: u th = 3RT µ (26) where T is temperature. The total thermal energy of mantle is calculated by integration: E th,mantle = Mp M core u th dm = M p 1 CMF u th dx (27) where x m M p. In the approximation, mantle density ρ mantle 1 r 1 m. On the other hand, with the assumption of adiabatic temperature gradient in the mantle and the introduction of the Grüneisen parameter γ ln(t) ln(ρ) adiabat,
10 10 u th = 3R T mp µ ( ) γ ρ = 3R T mp µ ρ 0 ( m M p ) γ 2 (28) where T mp (mantle potential temperature) is defined as the temperature when mantle adiabat is extrapolated to zero pressure. Here, γ = 1 is taken for simplicity. Then, E th,mantle = 2 M p 3R T mp (1 CMF) (29) µ whereas CMF 0.3 for the planets concerned, so CMF CRF 0.5 = 1, and total 2 thermal energy of mantle (mostly contributed by vibrations of atoms in crystal lattice, while electron contribution is ignored here) is: E th,mantle 3R T mp µ M p ( ) ( ) Tmp Mp J (30) 1000K M One can define an effective heat capacity C th,mantle of the mantle with regard to T mp : C th,mantle 3R µ M p ( Mp M ) J/K (31) Eq. 30 suggests that the thermal energy can be calculated by treating the mantle as uncompressed att mp ofmassm p. ForEarth, T mp 1700KandM p = kgandµ = 0.02 kg/mol (Stacey & Davis 2008), thus, E th,mantle, J and C th,mantle, J/K. Detailed calculation in Stacey & Davis (2008) shows the effective heat capacity of Earth mantle is J/K, very close to our estimate. Since the core is small, they can be regarded as approximations for thermal energy and heat capacity for the entire planet, and generally applied to other rocky planets having known their masses and mantle potential temperatures. Mantle dictates the cooling of the core (Stacey & Davis 2008). 6. Moment of Inertia Moment of inertia is calculated from the following formula, where x is the distance of the mass element to the rotational axis: I = ρ( r) x 2 dv (32) V
11 11 Let s consider two simple cases first: { thin spherical shell with radius Rp, I shell = 2 3 MR2 p uniform solid sphere with radius R p, I solid sphere = 2 5 MR2 p Defining C I/MR 2 p asthe moment of inertia factor, forshell C=2 3, andfor sphere C=2 5. Smaller C means mass is more concentrated towards the center. For the approximation, the moment of inertia of the core and the mantle can each be calculated separately first, then combined to give the full moment of inertia of the planet: I core = 2 5 M pr 2 p CMF2 (33) I mantle = 1 3 M pr 2 p (1 CMF 2 ) (34) I total = I core +I mantle = 1 3 M prp ( ) CMF2 (35) Considering only 0.2 CMF 0.35, the term 1 5 CMF2 can be ignored, and C 1. We 3 already know that in the solar system, C for Mercury, Venus, and Earth are indeed very close to 1 (Rubie et al. 2007). Here we generalize this fact of C 1 to all Earth-like rocky 3 3 planets: I planet 1 3 M p R 2 p (36) For Earth, M R 2 = kg m 2, and I 1 3 M R 2 = kg m 2. Therefore, angular momentum of Earth s rotation is L = I Ω = kg m 2 /s. And the total rotationalenergyofearthise rot = 2 I 1 Ω 2 = L2 2 Ω 2 is the angular frequency of Earth s rotation. = J, whereω = rad/sec 7. More Applications The aforementioned 5 results are just examples of what can be derived from this model, while it can be used to derive more results as you wish. The model itself only depends on 3 parameters: (1) surface gravity g s (2) planet radius R p (3) core radius fraction CRF. The procedure is as follows:
12 12 1. Given the mass M p and radius R p of a rocky planet, the surface gravity g s = G Mp R 2 p be calculated immediately. can 2. Then, CMF can be determined from the mass-radius relations described in (Zeng et al. 2016; Zeng & Jacobsen 2016b). 3. CRF can be estimated as CMF. Then all 3 parameters are obtained, and one may calculate the approximate properties of this planet as one wishes. 8. Conclusion Here we present a simple approximate model for Earth-like rocky planetary interiors. It predicts (1) CMF CRF 2 (core radius fraction squared), (2) P c gs 2 (interior pressure scales as surface gravity squared), (3) energy released of core formation ( is 1 the total 10 gravitational energy, (4) effective heat capacity of the mantle Mp J/K, and (5) moment of inertia 1 M 3 p Rp 2. Other results may be derived following the procedure. Although being approximate, these relations are straightforward to apply, as in many cases mass and radius are only measured approximately anyway. When they are combined with mass-radius relations of (Zeng et al. 2016; Zeng & Jacobsen 2016b), they can provide us more insights and intuitions towards rocky planetary interiors, and complements the numerical approach. M ) 9. Acknowledgement This work was supported by a grant from the Simons Foundation(SCOL[award#337090] to L.Z.). Part of this research was conducted under the Sandia Z Fundamental Science Program and supported by the Department of Energy National Nuclear Security Administration under Award Numbers de-na and de-na to S. B. Jacobsen (PI) with Harvard University. This research is the authors views and not those of the DOE. The authors would like to thank Dimitar D. Sasselov for insightful suggestions and helpful comments on this paper.
13 13 REFERENCES Ballard, S., Huber, D., Chaplin, W., et al. 2014, ApJ Birch, F. 1947, Physical Review, 71, , J. Geophys. Res., 57, 227 Dressing, C. D., Charbonneau, D., Dumusque, X., et al. 2015, ApJ, 800, 135 Dziewonski, A. M., & Anderson, D. L. 1981, Physics of the Earth and Planetary Interiors, 25, 297 Haswell, C. A. 2010, Transiting Exoplanets Rubie, D., Nimmo, F., & Melosh, H. 2007, in Treatise on Geophysics, ed. G. Schubert (Amsterdam: Elsevier), Seager, S., Kuchner, M., Hier-Majumder, C. A., & Militzer, B. 2007, The Astrophysical Journal, 669, 1279 Southworth, J., Wheatley, P. J., & Sams, G. 2007, MNRAS, 379, L11 Stacey, F. D., & Davis, P. M Zeng, L., & Jacobsen, S. B. 2016b, submitted to ApJ Zeng, L., & Sasselov, D. 2013, Publications of the Astronomical Society of the Pacific, 125, pp. 227 Zeng, L., Sasselov, D. D., & Jacobsen, S. B. 2016, ApJ, 819, 127 Zeng, L., & Seager, S. 2008, Publications of the Astronomical Society of the Pacific, 120, 983 This preprint was prepared with the AAS L A TEX macros v5.2.
A Simple Analytical Model for Rocky Planet Interiors
A Simple Analytical Model for Rocky Planet Interiors Li Zeng 1,a and Stein B. Jacobsen 1,b 1 Department of Earth and Planetary Sciences, Harvard University, Cambridge, MA 02138 a astrozeng@gmail.com b
More informationCOMPOSITIONS OF HOT SUPER-EARTH ATMOSPHERES: EXPLORING KEPLER CANDIDATES
Draft version April 29, 2013 Preprint typeset using L A TEX style emulateapj v. 11/10/09 COMPOSITIONS OF HOT SUPER-EARTH ATMOSPHERES: EXPLORING KEPLER CANDIDATES Y. Miguel Max Planck Institut fuer Astronomie,
More informationA Computational Tool to Interpret the Bulk Composition of Solid Exoplanets based on Mass and Radius Measurements
PUBLICATIONS OF THE ASTRONOMICAL SOCIETY OF THE PACIFIC, 120:983 991, 2008 September 2008. The Astronomical Society of the Pacific. All rights reserved. Printed in U.S.A. A Computational Tool to Interpret
More informationThe structure of rocky planets
he structure of rocky planets Christophe Sotin Iron planets (Mercury) errestrial planets Ocean / Icy planets Icy Moons Uranus and Neptune Giant planets References: Sotin et al. (21) errestrial Planet Interiors;
More informationMass-Radius Relationships for Solid Exoplanets
Mass-Radius Relationships for Solid Exoplanets Sara Seager Massachusetts Institute of Technology See Seager, Kuchner, HierMichelson Summer Workshop 2007 Image credit: NASA/JPL-Caltech/R. Hurt (SSC) Majumder,
More informationFundamental Stellar Parameters. Radiative Transfer. Stellar Atmospheres
Fundamental Stellar Parameters Radiative Transfer Stellar Atmospheres Equations of Stellar Structure Basic Principles Equations of Hydrostatic Equilibrium and Mass Conservation Central Pressure, Virial
More informationExamination, course FY2450 Astrophysics Wednesday 23 rd May, 2012 Time:
Page 1 of 18 The Norwegian University of Science and Technology Department of Physics Contact person Name: Robert Hibbins Tel: 93551, mobile: 94 82 08 34 Examination, course FY2450 Astrophysics Wednesday
More informationMars Internal Structure, Activity, and Composition. Tilman Spohn DLR Institute of Planetary Research, Berlin
Mars Internal Structure, Activity, and Composition Tilman Spohn DLR Institute of Planetary Research, Berlin Interior Structure Interior Structure models aim at 2 the bulk chemistry of the planet the masses
More informationPressure Volume Temperature Equation of State
Pressure Volume Temperature Equation of State S.-H. Dan Shim ( ) Acknowledgement: NSF-CSEDI, NSF-FESD, NSF-EAR, NASA-NExSS, Keck Equations relating state variables (pressure, temperature, volume, or energy).
More informationAstronomy 6570 Physics of the Planets
Astronomy 6570 Physics of the Planets Planetary Rotation, Figures, and Gravity Fields Topics to be covered: 1. Rotational distortion & oblateness 2. Gravity field of an oblate planet 3. Free & forced planetary
More informationInformation on internal structure from shape, gravity field and rotation
Information on internal structure from shape, gravity field and rotation Seismological information is available only for the Earth and in limited amounts for the Moon. Various geodetic data put constraints
More informationarxiv: v1 [astro-ph.ep] 13 Nov 2013
Astronomy & Astrophysics manuscript no. ms c ESO 2018 September 16, 2018 On the radius of habitable planets Y. Alibert 1,2 1 Physikalisches Institut & Center for Space and Habitability, Universitaet Bern,
More informationTransits of planets: mean densities
Chapter 3 Transits of planets: mean densities Close-in (short period) planets have a relatively high chance to transit in front of the star. A transit introduces a small periodic dimming of the star which
More information### dv where ρ is density, R is distance from rotation axis, and dv is
Comments This is one of my favorite problem assignments in our Solid Earth Geophysics class, typically taken by junior and senior concentrators and by first-year graduate students. I encourage students
More informationPlanetary Interiors. Ulrich Christensen
Planetary Interiors Ulrich Christensen Earth as a prototype planet Informations from shape, gravity and rotation Internal structure of terrestrial planets and icy moons The interior of gas planets Introduction
More informationEquations of Stellar Structure
Equations of Stellar Structure Stellar structure and evolution can be calculated via a series of differential equations involving mass, pressure, temperature, and density. For simplicity, we will assume
More informationStructure and evolution of (giant) exoplanets: some news from the theoretical front. I. Baraffe University of Exeter
Structure and evolution of (giant) exoplanets: some news from the theoretical front I. Baraffe University of Exeter I) Structure of Jupiter and Saturn II) Exoplanets: Interior structure and evolutionary
More informationarxiv: v2 [astro-ph.ep] 22 Jun 2015
Preprint 23 June 2015 Compiled using MNRAS LATEX style file v3.0 A Variable Polytrope Index Applied to Planet and Material Models S. P. Weppner, 1 J. P. McKelvey, 2 K. D. Thielen 1 A. K. Zielinski 1 1
More informationEART164: PLANETARY ATMOSPHERES
EART16: PLANETARY ATMOSPHERES Francis Nimmo Last Week How do planets form? They accrete from the solar nebula (dust+gas) They may subsequently migrate Where do atmospheres come from? Primary, secondary,
More informationORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points)
ORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points) PROBLEM 1 A) Summarize the content of the three Kepler s Laws. B) Derive any two of the Kepler
More informationEART162: PLANETARY INTERIORS
EART162: PLANETARY INTERIORS Francis Nimmo Last Week Global gravity variations arise due to MoI difference (J 2 ) We can also determine C, the moment of inertia, either by observation (precession) or by
More informationMinimum Radii of Super-Earths: Constraints from Giant Impacts
Minimum Radii of Super-Earths: Constraints from Giant Impacts Robert A. Marcus 1,a, Dimitar Sasselov 1, Lars Hernquist 1, Sarah T. Stewart 2 1 Astronomy Department, Harvard University, Cambridge, MA 02138
More informationAstro Instructors: Jim Cordes & Shami Chatterjee.
Astro 2299 The Search for Life in the Universe Lecture 8 Last time: Formation and function of stars This time (and probably next): The Sun, hydrogen fusion Virial theorem and internal temperatures of stars
More informationAy123 Set 1 solutions
Ay13 Set 1 solutions Mia de los Reyes October 18 1. The scale of the Sun a Using the angular radius of the Sun and the radiant flux received at the top of the Earth s atmosphere, calculate the effective
More informationarxiv: v1 [astro-ph.ep] 12 Dec 2016
Astronomy & Astrophysics manuscript no. manuscript_arxiv c ESO 2018 December 20, 2018 Frequency-dependent tidal dissipation in a viscoelastic Saturnian core and expansion of Mimas semi-major axis arxiv:1612.03664v1
More informationPrecise Stellar Parameters are Crucial for Constraining Transiting Exoplanet Interiors
Precise Stellar Parameters are Crucial for Constraining Transiting Exoplanet Interiors Andrew Vanderburg! NASA Sagan Fellow The University of Texas at Austin Juliette C. Becker, Lars A. Buchhave, Annelies
More informationGS388 Handout: Radial density distribution via the Adams-Williamson equation 1
GS388 Handout: Radial density distribution via the Adams-Williamson equation 1 TABLE OF CONTENTS ADIABATIC COMPRESSION: THE ADAMS WILLIAMSON EQUATION...1 EFFECT OF NON-ADIABATIC TEMPERATURE GRADIENT...3
More informationarxiv: v1 [astro-ph.ep] 29 Aug 2018
Draft version August, 8 Typeset using L A TEX preprint style in AASTeX6 THRESHOLD RADII OF VOLATILE-RICH PLANETS M. Lozovsky, R. Helled, C. Dorn, and J. Venturini arxiv:88.987v [astro-ph.ep] 9 Aug 8 Center
More informationEART162: PLANETARY INTERIORS
EART162: PLANETARY INTERIORS Francis Nimmo Last Week Applications of fluid dynamics to geophysical problems Navier-Stokes equation describes fluid flow: Convection requires solving the coupled equations
More informationInterior Structure of Rocky and Vapor Worlds
Interior Structure of Rocky and Vapor Worlds Diana Valencia, 4 April 2011 NASA Sagan Fellow, MIT Exploring Strange New Worlds: From Giant Planets to Super-Earths Super-Earths in the context of exoplanets
More informationω = 0 a = 0 = α P = constant L = constant dt = 0 = d Equilibrium when: τ i = 0 τ net τ i Static Equilibrium when: F z = 0 F net = F i = ma = d P
Equilibrium when: F net = F i τ net = τ i a = 0 = α dp = 0 = d L = ma = d P = 0 = I α = d L = 0 P = constant L = constant F x = 0 τ i = 0 F y = 0 F z = 0 Static Equilibrium when: P = 0 L = 0 v com = 0
More informationASTR 610 Theory of Galaxy Formation Lecture 15: Galaxy Interactions
ASTR 610 Theory of Galaxy Formation Lecture 15: Galaxy Interactions Frank van den Bosch Yale University, spring 2017 Galaxy Interactions & Transformations In this lecture we discuss galaxy interactions
More informationThermal Evolution of Terrestrial Exoplanets. System Earth Modelling Masters Thesis
Thermal Evolution of Terrestrial Exoplanets System Earth Modelling Masters Thesis Supervisor: Dr. A. P. van den Berg Institute of Earth Sciences Utrecht University Yue Zhao August 30, 2011 Abstract The
More informationStatic Stellar Structure
Most of the Life of A Star is Spent in Equilibrium Evolutionary Changes are generally slow and can usually be handled in a quasistationary manner We generally assume: Hydrostatic Equilibrium Thermodynamic
More information3 Hydrostatic Equilibrium
3 Hydrostatic Equilibrium Reading: Shu, ch 5, ch 8 31 Timescales and Quasi-Hydrostatic Equilibrium Consider a gas obeying the Euler equations: Dρ Dt = ρ u, D u Dt = g 1 ρ P, Dɛ Dt = P ρ u + Γ Λ ρ Suppose
More informationINTRODUCTION TO THE PHYSICS OF THE EARTH S INTERIOR
INTRODUCTION TO THE PHYSICS OF THE EARTH S INTERIOR SECOND EDITION JEAN-PAULPOIRIER Institut de Physique du Globe de Paris PUBLISHED BY THE PRESS SYNDICATE OF THE UNIVERSITY OF CAMBRIDGE The Pitt Building,
More informationFree-Fall Timescale of Sun
Free-Fall Timescale of Sun Free-fall timescale: The time it would take a star (or cloud) to collapse to a point if there was no outward pressure to counteract gravity. We can calculate the free-fall timescale
More informationPressure Volume Temperature (P-V-T) Relationships and Thermo elastic Properties of Geophysical Minerals
Pressure Volume Temperature (P-V-T) Relationships and Thermo elastic Properties of Geophysical Minerals A PROPOSAL FOR Ph.D PROGRAMME BY MONIKA PANWAR UNDER THE SUPERVISION OF DR SANJAY PANWAR ASSISTANT
More informationThermal and compositional structure of the Mantle and Lithosphere
Chapter 1 Thermal and compositional structure of the Mantle and Lithosphere 1.1 Primordial heat of the Earth The most widely accepted planetary formation theory says that the solar system accreted from
More informationB ν (T) = 2hν3 c 3 1. e hν/kt 1. (4) For the solar radiation λ = 20µm photons are in the Rayleigh-Jean region, e hν/kt 1+hν/kT.
Name: Astronomy 18 - Problem Set 8 1. Fundamental Planetary Science problem 14.4 a) Calculate the ratio of the light reflected by Earth at 0.5 µm to that emitted by the Sun at the same wavelength. The
More informationELASTICITY AND CONSTITUTION OF THE EARTH'S INTERIOR*
JOURNAL OF GEOPHYSICAL RESEARCH VOLUME 57, NO. 2 JUNE, 1952 ELASTICITY AND CONSTITUTION OF THE EARTH'S INTERIOR* BY FRANCIS BIRCtt Harvard University, Cambridge, Massachusetts (Received January 18, 1952)
More informationAST1100 Lecture Notes
AST1100 Lecture Notes 5 The virial theorem 1 The virial theorem We have seen that we can solve the equation of motion for the two-body problem analytically and thus obtain expressions describing the future
More informationTerrestrial Exoplanet Radii, Structure and Tectonics. Diana Valencia (OCA), KITP, 31 March 2010
Terrestrial Exoplanet Radii, Structure and Tectonics Diana Valencia (OCA), KITP, 31 March 2010 Earth s Structure Olivine Structure Romanowicz, 2008 Earth s Structure Upper mantle (Mg,Fe)2SiO4 α +(Mg,Fe)2Si2O6
More informationThe structure and evolution of stars. Introduction and recap
The structure and evolution of stars Lecture 3: The equations of stellar structure 1 Introduction and recap For our stars which are isolated, static, and spherically symmetric there are four basic equations
More informationComputational support for a pyrolitic lower mantle containing ferric iron
SUPPLEMENTARY INFORMATION DOI: 1.138/NGEO2458 Computational support for a pyrolitic lower mantle containing ferric iron Xianlong Wang, Taku Tsuchiya and Atsushi Hase NATURE GEOSCIENCE www.nature.com/naturegeoscience
More informationExplaining Planetary-Rotation Periods Using an Inductive Method
Explaining Planetary-Rotation Periods Using an Inductive Method Gizachew Tiruneh, Ph. D., Department of Political Science, University of Central Arkansas June 19, 2009 This paper uses an inductive method
More informationAdvanced Newtonian gravity
Foundations of Newtonian gravity Solutions Motion of extended bodies, University of Guelph h treatment of Newtonian gravity, the book develops approximation methods to obtain weak-field solutions es the
More informationNewton s Gravitational Law
1 Newton s Gravitational Law Gravity exists because bodies have masses. Newton s Gravitational Law states that the force of attraction between two point masses is directly proportional to the product of
More informationBinary star formation
Binary star formation So far we have ignored binary stars. But, most stars are part of binary systems: Solar mass stars: about 2 / 3 are part of binaries Separations from: < 0.1 au > 10 3 au Wide range
More informationCHAPTER X. Second Half Review 2017
CHAPTER X Second Half Review 217 Here is a quick overview of what we covered in the second half of the class. Remember that the final covers the whole course but there will naturally be a bias towards
More informationPhysics H7A, Fall 2011 Homework 6 Solutions
Physics H7A, Fall 2011 Homework 6 Solutions 1. (Hooke s law) (a) What are the dimensions (in terms of M, L, and T ) of the spring constant k in Hooke s law? (b) If it takes 4.00 J of work to stretch a
More informationUniform Circular Motion
Circular Motion Uniform Circular Motion Uniform Circular Motion Traveling with a constant speed in a circular path Even though the speed is constant, the acceleration is non-zero The acceleration responsible
More informationChapter 4 Making Sense of the Universe: Understanding Motion, Energy, and Gravity. Copyright 2009 Pearson Education, Inc.
Chapter 4 Making Sense of the Universe: Understanding Motion, Energy, and Gravity How do we describe motion? Precise definitions to describe motion: Speed: Rate at which object moves speed = distance time
More information2. Equations of Stellar Structure
2. Equations of Stellar Structure We already discussed that the structure of stars is basically governed by three simple laws, namely hyostatic equilibrium, energy transport and energy generation. In this
More informationAY202a Galaxies & Dynamics Lecture 7: Jeans Law, Virial Theorem Structure of E Galaxies
AY202a Galaxies & Dynamics Lecture 7: Jeans Law, Virial Theorem Structure of E Galaxies Jean s Law Star/Galaxy Formation is most simply defined as the process of going from hydrostatic equilibrium to gravitational
More information[2 marks] Show that derivative of the angular velocity. What is the specific angular momentum j as a function of M and R in this Keplerian case?
k!! Queen Mary University of London M. Sc. EXAM I N AT1 0 N ASTMOOS Angular Momentum and Accretion in Astrophysics Fkiday, 26th May, 2006 18:15-19:45 Time Allowed: lh 30m This paper has two Sections and
More informationII Planet Finding.
II Planet Finding http://sgoodwin.staff.shef.ac.uk/phy229.html 1.0 Introduction There are a lot of slides in this lecture. Much of this should be familiar from PHY104 (Introduction to Astrophysics) and
More informationDepartment of Physics PRELIMINARY EXAMINATION 2016 Part I. Short Questions
Department of Physics PRELIMINARY EXAMINATION 2016 Part I. Short Questions Thursday May 19th, 2016, 14-17h Examiners: Prof. J. Cline, Prof. H. Guo, Prof. G. Gervais (Chair), and Prof. D. Hanna INSTRUCTIONS
More informationDETAILED MODELS OF SUPER-EARTHS: HOW WELL CAN WE INFER BULK PROPERTIES?
The Astrophysical Journal, 665:1413 1420, 2007 August 20 # 2007. The American Astronomical Society. All rights reserved. Printed in U.S.A. A DETAILED MODELS OF SUPER-EARTHS: HOW WELL CAN WE INFER BULK
More informationAy101 Set 1 solutions
Ay11 Set 1 solutions Ge Chen Jan. 1 19 1. The scale of the Sun a 3 points Venus has an orbital period of 5 days. Using Kepler s laws, what is its semi-major axis in units of AU Start with Kepler s third
More informationPHYSICS 231 INTRODUCTORY PHYSICS I
PHYSICS 231 INTRODUCTORY PHYSICS I Lecture 12 Last Lecture Newton s Law of gravitation F grav = GMm R 2 G = 6.67"10 #11 Nm 2 /kg 2 Kepler s Laws of Planetary motion 1. Ellipses with sun at focus 2. Sweep
More informationProblem set: solar irradiance and solar wind
Problem set: solar irradiance and solar wind Karel Schrijver July 3, 203 Stratification of a static atmosphere within a force-free magnetic field Problem: Write down the general MHD force-balance equation
More informationThe Interiors of the Stars
The Interiors of the Stars Hydrostatic Equilibrium Stellar interiors, to a good first approximation, may be understood using basic physics. The fundamental operating assumption here is that the star is
More informationAssignment 1. Due Feb. 11, 2019
Assignment 1 Due Feb. 11, 2019 Show all work and turn in answers on separate pages, not on these pages. Circle your final answers for clarity. Be sure to show/explain all of your reasoning and that your
More informationread 9.4-end 9.8(HW#6), 9.9(HW#7), 9.11(HW#8) We are proceding to Chap 10 stellar old age
HW PREVIEW read 9.4-end Questions 9.9(HW#4), 9(HW#4) 9.14(HW#5), 9.8(HW#6), 9.9(HW#7), 9.11(HW#8) We are proceding to Chap 10 stellar old age Chap 11 The death of high h mass stars Contraction of Giant
More informationData from: The Extrasolar Planet Encyclopaedia.
Data from: The Extrasolar Planet Encyclopaedia http://exoplanet.eu/ 2009->10 Status of Exoplanet Searches Direct Detection: 5->9 planets detected Sensitive to large planets in large orbits around faint
More informationSilicate Atmospheres, Clouds, and Fractional Vaporization of Hot Earth-like Exoplanets
Silicate Atmospheres, Clouds, and Fractional Vaporization of Hot Earth-like Exoplanets Laura Schaefer and Bruce Fegley, Jr. Planetary Chemistry Laboratory Department of Earth and Planetary Sciences Washington
More informationGrade XI. Physics Exam Preparation Booklet. Chapter-wise Important Questions. #GrowWithGreen
Grade XI Physics Exam Preparation Booklet Chapter-wise Important Questions #GrowWithGreen Units and Measurements Q1. After reading the physics book, Anamika recalled and noted down the expression for the
More informationarxiv: v1 [physics.pop-ph] 18 May 2015
A Guided Tour of Planetary Interiors arxiv:1505.05894v1 [physics.pop-ph] 18 May 2015 Alexander. Klotz Department of Physics, McGill University 3600 rue University, Montreal, Quebec H3A 2T8, Canada klotza@physics.mcgill.ca
More informationSolved examples of Gravitation
Solved examples of Gravitation Example 1 The time period of Moon around the Earth is n times that of Earth around the Sun. If the ratio of the distance of the Earth from the Sun to that of the distance
More informationPHYSICS OF ASTROPHSYICS - Energy.
PHYSICS OF ASTROPHSYICS - Energy http://apod.nasa.gov/apod/ ENERGY Result of a force acting through a distance. units = erg = dyne cm i.e., force x distance = gm cm 2 /sec 2 Two types: kinetic - energy
More informationLecture 3. Basic Physics of Astrophysics - Force and Energy.
Lecture 3 Basic Physics of Astrophysics - Force and Energy http://apod.nasa.gov/apod/ Tue Wed Thur Thur ISB 165 Wed 5 Thur 4 Momentum (generally m F = Forces is the product of mass and velocity - a vector
More informationGravitation. Luis Anchordoqui
Gravitation Kepler's law and Newton's Synthesis The nighttime sky with its myriad stars and shinning planets has always fascinated people on Earth. Towards the end of the XVI century the astronomer Tycho
More informationElliptical-like orbits on a warped spandex fabric
Elliptical-like orbits on a warped spandex fabric Prof. Chad A. Middleton CMU Physics Seminar September 3, 2015 Submitted to the American Journal of Physics Kepler s 3 Laws of planetary motion Kepler s
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Earth, Atmospheric, and Planetary Sciences Department. Final Exam
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Earth, Atmospheric, and Planetary Sciences Department Physics 8.282J EAPS 12.402J May 20, 2005 Final Exam Name Last First (please print) 1. Do any
More informationChapter 9 Circular Motion Dynamics
Chapter 9 Circular Motion Dynamics Chapter 9 Circular Motion Dynamics... 9. Introduction Newton s Second Law and Circular Motion... 9. Universal Law of Gravitation and the Circular Orbit of the Moon...
More informationL = I ω = const. I = 2 3 MR2. When the balloon shrinks (because the air inside it cools down), the moment of inertia decreases, R = 1. L = I ω = I ω.
PHY 30 K. Solutions for mid-term test #3. Problem 1: Out in space, there are no forces acting on the balloon except gravity and hence no torques (with respect to an axis through the center of mass). Consequently,
More informationSPECIFIC ANGULAR MOMENTUM DISTRIBUTION FOR SOLAR ANALOGS AND TWINS: WHERE IS THE SUN MISSING HIS ANGULAR MOMENTUM?
5th International Workshop on Astronomy and Relativistic Astrophysics (IWARA2011) International Journal of Modern Physics: Conference Series Vol. 18 (2012) 58 62 c World Scientific Publishing Company DOI:
More informationarxiv: v1 [astro-ph.ep] 23 Mar 2010
Formation of Terrestrial Planets from Protoplanets under a Realistic Accretion Condition arxiv:1003.4384v1 [astro-ph.ep] 23 Mar 2010 Eiichiro Kokubo Division of Theoretical Astronomy, National Astronomical
More informationComputational Modeling of Volcanism on Earth-like
PHYSICS Computational Modeling of Volcanism on Earth-like Planets EDWARD KIM MATTHEW J.S. BEACH ETHAN BUCHMAN ABSTRACT The aim of this paper is to model the level of volcanic activity on Earth-like planets
More informationTHIRD-YEAR ASTROPHYSICS
THIRD-YEAR ASTROPHYSICS Problem Set: Stellar Structure and Evolution (Dr Ph Podsiadlowski, Michaelmas Term 2006) 1 Measuring Stellar Parameters Sirius is a visual binary with a period of 4994 yr Its measured
More informationarxiv: v1 [astro-ph] 19 Jul 2007
Draft version February 1, 2008 Preprint typeset using L A TEX style emulateapj v. 11/12/01 MASS-RADIUS RELATIONSHIPS FOR SOLID EXOPLANETS S. Seager 1,2, M. Kuchner 3, C. A. Hier-Majumder 1, B. Militzer
More informationEOS-FIT V6.0 R.J. ANGEL
EOS-FIT V6. R.J. AGEL Crystallography Laboratory, Dept. Geological Sciences, Virginia Tech, Blacksburg, VA46, USA http://www.geol.vt.edu/profs/rja/ ITRODUCTIO EosFit started as a program to fit equations
More informationChapter 13. Gravitation
Chapter 13 Gravitation e = c/a A note about eccentricity For a circle c = 0 à e = 0 a Orbit Examples Mercury has the highest eccentricity of any planet (a) e Mercury = 0.21 Halley s comet has an orbit
More informationarxiv:astro-ph/ v2 31 Jan 2007
THE PRECESSION OF MERCURY AND THE DEFLECTION OF STARLIGHT FROM SPECIAL RELATIVITY ALONE Robert L. Kurucz Harvard-Smithsonian Center for Astrophysics arxiv:astro-ph/0608434v2 31 Jan 2007 60 Garden St, Cambridge,
More informationAstr 1050 Fri., Feb. 24, 2017
Astr 1050 Fri., Feb. 24, 2017 Chapter 7 & 8: Overview & Formation of the Solar System Reading: Chapters 7 on Solar System Chapter 8: Earth & Terrestrial Planets Reminders: New homework on MA up this afternoon,
More informationCHAPTER 10 GRAVITATION
CHAPTER 10 GRAVITATION Earth attracts everything towards it by an unseen force of attraction. This force of attraction is known as gravitation or gravitation pull. Universal Law of Gravitation:- Every
More informationarxiv: v1 [astro-ph.sr] 8 Aug 2014
Magnetic Modulation of Stellar Angular Momentum Loss Cecilia Garraffo 1, Jeremy J. Drake 1, and Ofer Cohen 1 ABSTRACT arxiv:1408.1965v1 [astro-ph.sr] 8 Aug 2014 Angular Momentum Loss is important for understanding
More information43. A person sits on a freely spinning lab stool that has no friction in its axle. When this person extends her arms,
43. A person sits on a freely spinning lab stool that has no friction in its axle. When this person extends her arms, A) her moment of inertia increases and her rotational kinetic energy remains the same.
More informationSelf-Similar Gas Dynamics of Voids and Supernova Ejecta
Self-Similar Gas Dynamics of Voids and Supernova Ejecta arxiv:1109.2682v2 [astro-ph.sr] Li-Le Wang Director: Prof. Yu-Qing Lou Department of Physics and Tsinghua Center for Astrophysics (THCA) Tsinghua
More informationStellar structure Conservation of mass Conservation of energy Equation of hydrostatic equilibrium Equation of energy transport Equation of state
Stellar structure For an isolated, static, spherically symmetric star, four basic laws / equations needed to describe structure: Conservation of mass Conservation of energy (at each radius, the change
More informationCanadian Journal of Physics. Higher Order derivatives of bulk modulus of materials at infinite pressure
Higher Order derivatives of bulk modulus of materials at infinite pressure Journal: Manuscript ID cjp-015-0699.r Manuscript Type: Article Date Submitted by the Author: 0-Apr-016 Complete List of Authors:
More informationAST1100 Lecture Notes
AST1100 Lecture Notes 20: Stellar evolution: The giant stage 1 Energy transport in stars and the life time on the main sequence How long does the star remain on the main sequence? It will depend on the
More informationVISUAL PHYSICS ONLINE
VISUAL PHYSICS ONLINE EXCEL SIMULATION MOTION OF SATELLITES DOWNLOAD the MS EXCEL program PA50satellite.xlsx and view the worksheet Display as shown in the figure below. One of the most important questions
More informationExtrasolar planets and their hosts: Why exoplanet science needs X-ray observations
Extrasolar planets and their hosts: Why exoplanet science needs X-ray observations Dr. Katja Poppenhaeger Sagan Fellow Harvard-Smithsonian Center for Astrophysics Exoplanets Exoplanets in 2005 (from the
More informationequals the chemical potential µ at T = 0. All the lowest energy states are occupied. Highest occupied state has energy µ. For particles in a box:
5 The Ideal Fermi Gas at Low Temperatures M 5, BS 3-4, KK p83-84) Applications: - Electrons in metal and semi-conductors - Liquid helium 3 - Gas of Potassium 4 atoms at T = 3µk - Electrons in a White Dwarf
More information7 - GRAVITATION Page 1 ( Answers at the end of all questions )
7 - GRAVITATION Page 1 1 ) The change in the value of g at a height h above the surface of the earth is the same as at a depth d below the surface of earth. When both d and h are much smaller than the
More informationAST Homework V - Solutions
AST 341 - Homework V - Solutions TA: Marina von Steinkirch, steinkirch@gmail.com State University of New York at Stony Brook November, 010 1 (1 point) Derive the homologous form of the luminosity equation
More informationConstraining the Radius of Neutron Stars Through the Moment of Inertia
Constraining the Radius of Neutron Stars Through the Moment of Inertia Neutron star mergers: From gravitational waves to nucleosynthesis International Workshop XLV on Gross Properties of Nuclei and Nuclear
More informationE = K + U. p mv. p i = p f. F dt = p. J t 1. a r = v2. F c = m v2. s = rθ. a t = rα. r 2 dm i. m i r 2 i. I ring = MR 2.
v = v i + at x = x i + v i t + 1 2 at2 E = K + U p mv p i = p f L r p = Iω τ r F = rf sin θ v 2 = v 2 i + 2a x F = ma = dp dt = U v dx dt a dv dt = d2 x dt 2 A circle = πr 2 A sphere = 4πr 2 V sphere =
More information