PHYS 633: Introduction to Stellar Astrophysics Spring Semester 2006 Rich Townsend ([email protected])

Total Page:16

File Type:pdf, Size:1020Kb

PHYS 633: Introduction to Stellar Astrophysics Spring Semester 2006 Rich Townsend (Rhdt@Bartol.Udel.Edu) PHYS 633: Introduction to Stellar Astrophysics Spring Semester 2006 Rich Townsend ([email protected]) Governing Equations In the foregoing analysis, we have seen how a star usually maintains a state of hydrostatic equilibrium. We have seen how the energy created in the star through nuclear reactions, or released through contraction, makes its way out in the form of the star’s luminosity. We have seen how the physical processes responsible for transporting this energy can either be radiation or convection. And, finally, we have seen how nuclear reactions within the star, and transport processes such as convective mixing, can alter the stars chemical composition. These four key pieces of physics lead to 3 + I of the differential equations governing stellar structure (here, I is the number of elements whose evolution we are following). Let’s quickly review these equations. First, we have the most general (spherically-symmetric) form for the equation governing momen- tum conservation, ∂P Gm 1 ∂2r = − − . (1) ∂m 4πr2 4πr2 ∂t2 If the second term on the right-hand side vanishes, then we recover the equation of hydrostatic equilibrium, expressed in Lagrangian form (i.e., derivatives with respect to the mass coordinate m, rather than the radial coordinate r). If this term does not vanish at some point in the star, then the mass shells at that point will being to expand or contract, with an acceleration given by ∂2r/∂t2. Second, we have the most general form for the energy conservation equation, ∂l ∂T δ ∂P = − − c + . (2) ∂m ν P ∂t ρ ∂t The first and second terms on the right-hand side come from nuclear energy gen- eration and neutrino losses, respectively. The third and fourth terms together compose what we have previously termed the gravothermal energy release rate, gi, P ∂ρ ∂u ∂T δ ∂P ≡ − = −c + ; (3) gi ρ2 ∂t ∂t P ∂t ρ ∂t they give a net positive contribution to the luminosity gradient when the star contracts. Third, we have the equation(s) of energy transport. These allow us to de- termine the physical temperature gradient ∇ ≡ ∂ ln T/∂ ln P at each point in the star, from which we can obtain the ‘real’ temperature gradient as ∂T GmT = − ∇. (4) ∂m 4πr4P 1 In radiative regions, the physical gradient is equal to the radiative gradient ∇rad, 3 κρH 3 κlP ∇ ≡ ∇ ≡ P F = . (5) rad rad 4ac T 4 16πacG mT 4 If ∇rad exceeds the adiabatic temperature gradient ∇rad, convection will set in, and we need instead to calculate ∇ from a convective model — usually, mixing-length theory. Finally, we have the I − 1 equations describing the change in the mass fractions describing the star’s chemical composition, ∂Xi mi X i X i X ji X ik = pjirji − pjirik = − . (6) ∂t ρ qji qik j k j k The final equation governing the chemical composition is the closure relation X ∂Xi = 0. (7) ∂t i To describe mixing in convective regions, the above equations are replaced by Z mb dX¯i ∂Xi dma dmb = dm − Xi,a(Xi,a − X¯i) + (Xi,b − X¯i). (8) dt ma ∂t dt dt There is one final differential equation we must add to the ones above — one based not so much on deep physics, and more on the definition of the mass coordinate m. This is of course the mass-radius relation ∂r 1 = (9) ∂m 4πr2ρ This equation, plus eqns. (1), (2), (4), (6) and (5), together form the 4 + I differential equations governing stellar structure. To solve them, we’re going to need informationr regarding the stellar material — for instance, its equation of state, thermodynamic coefficients, opacity, nuclear rates, and so on. We’re also going to need boundary conditions, to ensure that we obtain only those solutions that are physically-realistic. First, however, let’s step back and have a look over the equations. As Kippen- hahn & Weigert discussion in their §9, certain subsystems can be distinguished in the equations. For instance, the mechanical equations, (1) and (9), can be solved if the density can be expressed solely as a function of the pressure. One possible pressure-density relation is the polytropic relation, P = Kρ1+1/n; sub- stituting this into the mechanical equations leads to the Lane-Emden equation. Another pair of subsystems can be formed by splitting the chemical compo- sition equations from the others. This split is highlighted in the title of Kippen- hahn & Weigert: ‘Stellar structure and Evolution’ indicates two seperate topics, and it is not far from the truth to view the chemical composition equations as relating to evolution, while the remaining equations relate to stellar structure. 2 In a sence, the chemical composition is a star’s memory; it records where a star has been during its past evolution, and dictates where it will go in the future. As we shall discuss later, armed solely with knowledge of the total mass of a star and its internal chemical composition distribution, it is in principle possible to reconstruct its full internal structure. Let’s now look at the additional information we need to solve the differential equations. This comprises: • The equation of state ρ = ρ(P, T, Xi), (10) which relates the density at each point to the temperature and pressure, • Three thermodynamic coefficient equations, for instance cP = cP (P, T, Xi), (11) Γ1 = Γ1(P, T, Xi), (12) ∇rad = ∇rad(P, T, Xi), (13) from which the 21 other possible thermodynamic coefficients can be de- rived, • The equation for the Rosseland mean opacity (with conduction included), κ = κ(P, T, Xi), (14) which tells us how easily radiation can flow through stellar material, • The nuclear rate equations rik = rik(P, T, Xi), (15) = (P, T, Xi), (16) 3 = ν (P, T, Xi) (17) Combining these with the differential equations, we are able to write the right- hand sides of the differential equations purely in terms of r, P , T , l, and Xi — that is, 4+I unknowns. Since there are 4+I differential equations, we therefore have sufficient information to calculate a solution — at last! Equations. (1) and (2) both contain time derivatives as well as space deriva- tives, marking them out as mixed differential equations. Often, mixed equations can introduce quite a bit of difficulty in numerical solution techniques. However, in the present case there are a couple of simplifications we can make. Let’s first consider the momentum equation, (1). Suppose the pressure gradient term on the left-hand side could be magically set to zero (corresponding to ‘turning off’ the pressure supporting the star against gravity). The approximate timescale τdyn for the surface layers to collapse down to the center is given by the solution of GM R 0 = − 4 + 2 , (18) 4πR τdyn 3 which is r R τ ∼ . (19) dyn GM This ‘dynamical timescale’, which turns out to be around 20 minutes for the Sun, tells us how quickly a star will readjust itself into hydrostatic equilibrium, in response to a disturbance. It also tells us how important the time derivative term in eqn. (1) is. If the internal structure of a star is changing on a timescale much larger than τdyn, then we can neglect the time derivative, to give ∂P Gm = − . (20) ∂m 4πr2 Let’s do a similar thing for the energy equation. Suppose that the nuclear energy sources and neutrino losses went away, so that the star had to provide its luminosity from its gravitational and internal energy, ∂l ∂T δ ∂P = −c + . (21) ∂m P ∂t ρ ∂t On the left-hand side, we have approximately the energy being lost from the star per unit mass; and on the right hand side, we have approximately the rate- of-change of thermal and gravitational energy per unit mass supplying this loss. Hence, the star will run out of energy over a timescale L L τKH ∼ ∼ , (22) Eg Ei where we have used Eg and Ei before to denote the gravitational and ther- mal energy, respectively. This timescale is of course is the Kelvin-Helmholtz timescale. As with the dynamical timescale, the Kelvin-Helmholtz timescale tells us how quickly a star will respond to departures from thermal equilibrium (these departures coming, e.g., from losses through the surface that require a contraction of the star). If the internal structure of a star is changing on a timescale much larger than τKH, then we can neglect the time derivative in the energy equation, to give ∂l = − (23) ∂m ν Typically, a star evolves during nuclear burning over a nuclear timescale L L τn ∼ ∼ 2 , (24) En fMc where the mass-energy conversion fraction f is ∼ 10−4. This timescale is much longer than the dynamical and Kelvin-Helmholtz timescales, so the star re- mains in hydrostatic and thermal equilibrium during nuclear burning. Of course, when nuclear reactions cease, the star then evolves over the Kelvin-Helmholtz timescale (i.e., much more rapidly), until conditions are reached where fresh nu- clear burning can occur. During this rapid evolution, the star is not in thermal equilibrium, and we cannot neglect the time derivative in the energy equa- tion (2). 4.
Recommended publications
  • Used Time Scales GEORGE E
    PROCEEDINGS OF THE IEEE, VOL. 55, NO. 6, JIJNE 1967 815 Reprinted from the PROCEEDINGS OF THE IEEI< VOL. 55, NO. 6, JUNE, 1967 pp. 815-821 COPYRIGHT @ 1967-THE INSTITCJTE OF kECTRITAT. AND ELECTRONICSEN~INEERS. INr I’KINTED IN THE lT.S.A. Some Characteristics of Commonly Used Time Scales GEORGE E. HUDSON A bstract-Various examples of ideally defined time scales are given. Bureau of Standards to realize one international unit of Realizations of these scales occur with the construction and maintenance of time [2]. As noted in the next section, it realizes the atomic various clocks, and in the broadcast dissemination of the scale information. Atomic and universal time scales disseminated via standard frequency and time scale, AT (or A), with a definite initial epoch. This clock time-signal broadcasts are compared. There is a discussion of some studies is based on the NBS frequency standard, a cesium beam of the associated problems suggested by the International Radio Consultative device [3]. This is the atomic standard to which the non- Committee (CCIR). offset carrier frequency signals and time intervals emitted from NBS radio station WWVB are referenced ; neverthe- I. INTRODUCTION less, the time scale SA (stepped atomic), used in these emis- PECIFIC PROBLEMS noted in this paper range from sions, is only piecewise uniform with respect to AT, and mathematical investigations of the properties of in- piecewise continuous in order that it may approximate to dividual time scales and the formation of a composite s the slightly variable scale known as UT2. SA is described scale from many independent ones, through statistical in Section 11-A-2).
    [Show full text]
  • 14 Timescales in Stellar Interiors
    14Timescales inStellarInteriors Having dealt with the stellar photosphere and the radiation transport so rel- evant to our observations of this region, we’re now ready to journey deeper into the inner layers of our stellar onion. Fundamentally, the aim we will de- velop in the coming chapters is to develop a connection betweenM,R,L, and T in stars (see Table 14 for some relevant scales). More specifically, our goal will be to develop equilibrium models that describe stellar structure:P (r),ρ (r), andT (r). We will have to model grav- ity, pressure balance, energy transport, and energy generation to get every- thing right. We will follow a fairly simple path, assuming spherical symmetric throughout and ignoring effects due to rotation, magneticfields, etc. Before laying out the equations, let’sfirst think about some key timescales. By quantifying these timescales and assuming stars are in at least short-term equilibrium, we will be better-equipped to understand the relevant processes and to identify just what stellar equilibrium means. 14.1 Photon collisions with matter This sets the timescale for radiation and matter to reach equilibrium. It de- pends on the mean free path of photons through the gas, 1 (227) �= nσ So by dimensional analysis, � (228)τ γ ≈ c If we use numbers roughly appropriate for the average Sun (assuming full Table3: Relevant stellar quantities. Quantity Value in Sun Range in other stars M2 1033 g 0.08 �( M/M )� 100 × � R7 1010 cm 0.08 �( R/R )� 1000 × 33 1 3 � 6 L4 10 erg s− 10− �( L/L )� 10 × � Teff 5777K 3000K �( Teff/mathrmK)� 50,000K 3 3 ρc 150 g cm− 10 �( ρc/g cm− )� 1000 T 1.5 107 K 106 ( T /K) 108 c × � c � 83 14.Timescales inStellarInteriors P dA dr ρ r P+dP g Mr Figure 28: The state of hydrostatic equilibrium in an object like a star occurs when the inward force of gravity is balanced by an outward pressure gradient.
    [Show full text]
  • Stellar Structure and Evolution
    Lecture Notes on Stellar Structure and Evolution Jørgen Christensen-Dalsgaard Institut for Fysik og Astronomi, Aarhus Universitet Sixth Edition Fourth Printing March 2008 ii Preface The present notes grew out of an introductory course in stellar evolution which I have given for several years to third-year undergraduate students in physics at the University of Aarhus. The goal of the course and the notes is to show how many aspects of stellar evolution can be understood relatively simply in terms of basic physics. Apart from the intrinsic interest of the topic, the value of such a course is that it provides an illustration (within the syllabus in Aarhus, almost the first illustration) of the application of physics to “the real world” outside the laboratory. I am grateful to the students who have followed the course over the years, and to my colleague J. Madsen who has taken part in giving it, for their comments and advice; indeed, their insistent urging that I replace by a more coherent set of notes the textbook, supplemented by extensive commentary and additional notes, which was originally used in the course, is directly responsible for the existence of these notes. Additional input was provided by the students who suffered through the first edition of the notes in the Autumn of 1990. I hope that this will be a continuing process; further comments, corrections and suggestions for improvements are most welcome. I thank N. Grevesse for providing the data in Figure 14.1, and P. E. Nissen for helpful suggestions for other figures, as well as for reading and commenting on an early version of the manuscript.
    [Show full text]
  • Evolution of the Sun, Stars, and Habitable Zones
    Evolution of the Sun, Stars, and Habitable Zones AST 248 fs The fraction of suitable stars N = N* fs fp nh fl fi fc L/T Hertzsprung-Russell Diagram Parts of the H-R Diagram •Supergiants •Giants •Main Sequence (dwarfs) •White Dwarfs Making Sense of the H-R Diagram •The Main Sequence is a sequence in mass •Stars on the main sequence undergo stable H fusion •All other stars are evolved •Evolved stars have used up all their core H •Main Sequence ® Giants ® Supergiants •Subsequent evolution depends on mass Hertzsprung-Russell Diagram Evolutionary Timescales Pre-main sequence: Set by gravitational contraction •The gravitational potential energy E is ~GM2/R •The luminosity is L •The timescale is ~E/L We know L, M, R from observations For the Sun, L ~ 30 million years Evolutionary Timescales Main sequence: •Energy source: nuclear reactions, at ~10-5 erg/reaction •Luminosity: 4x1033 erg/s This requires 4x1038 reactions/second Each reaction converts 4 H ® He 56 The solar core contains 0.1 M¤, or ~10 H atoms 1056 atoms / 4x1038 reactions/second -> 3x1017 sec, or 1010 years. This is the nuclear timescale. Stellar Lifetimes τ ~ M/L On the main sequence, L~M3 Therefore, τ~M-2 10 τ¤ = 10 years τ~1010/M2 years Lower mass stars live longer than the Sun Post-Main Sequence Timescales Timescale τ ~ E/L L >> Lms τ << τms Habitable Zones Refer back to our discussion of the Greenhouse Effect. 2 0.25 Tp ~ (L*/D ) The habitable zone is the region where the temperature is between 0 and 100 C (273 and 373 K), where liquid water can exist.
    [Show full text]
  • Useful Constants
    Appendix A Useful Constants A.1 Physical Constants Table A.1 Physical constants in SI units Symbol Constant Value c Speed of light 2.997925 × 108 m/s −19 e Elementary charge 1.602191 × 10 C −12 2 2 3 ε0 Permittivity 8.854 × 10 C s / kgm −7 2 μ0 Permeability 4π × 10 kgm/C −27 mH Atomic mass unit 1.660531 × 10 kg −31 me Electron mass 9.109558 × 10 kg −27 mp Proton mass 1.672614 × 10 kg −27 mn Neutron mass 1.674920 × 10 kg h Planck constant 6.626196 × 10−34 Js h¯ Planck constant 1.054591 × 10−34 Js R Gas constant 8.314510 × 103 J/(kgK) −23 k Boltzmann constant 1.380622 × 10 J/K −8 2 4 σ Stefan–Boltzmann constant 5.66961 × 10 W/ m K G Gravitational constant 6.6732 × 10−11 m3/ kgs2 M. Benacquista, An Introduction to the Evolution of Single and Binary Stars, 223 Undergraduate Lecture Notes in Physics, DOI 10.1007/978-1-4419-9991-7, © Springer Science+Business Media New York 2013 224 A Useful Constants Table A.2 Useful combinations and alternate units Symbol Constant Value 2 mHc Atomic mass unit 931.50MeV 2 mec Electron rest mass energy 511.00keV 2 mpc Proton rest mass energy 938.28MeV 2 mnc Neutron rest mass energy 939.57MeV h Planck constant 4.136 × 10−15 eVs h¯ Planck constant 6.582 × 10−16 eVs k Boltzmann constant 8.617 × 10−5 eV/K hc 1,240eVnm hc¯ 197.3eVnm 2 e /(4πε0) 1.440eVnm A.2 Astronomical Constants Table A.3 Astronomical units Symbol Constant Value AU Astronomical unit 1.4959787066 × 1011 m ly Light year 9.460730472 × 1015 m pc Parsec 2.0624806 × 105 AU 3.2615638ly 3.0856776 × 1016 m d Sidereal day 23h 56m 04.0905309s 8.61640905309
    [Show full text]
  • 7A: Introduction to Astrophysics 2016
    7A: Introduction to Astrophysics 2016 General Information: • Instructor: Mariska Kriek ([email protected]) • Graduate Student Instructors: o Nick Kern ([email protected]) o Michael Medford ([email protected]) • Lectures: Tuesday & Thursday 11:00 am - 12:30 pm, 131 Campbell Hall • Discussion sections: - 101: Monday 1-2 pm, 121 Campbell - Nick Kern - 102: Monday 4-5 pm, 121 Campbell - Michael Medford • The Astronomy Learning Center (TALC); Thursday 5-7 pm, 121 Campbell This is a large, collaborative "office hour" where students work on their homework assignments in an informal group setting. TALC is staffed by GSIs who serve as guides, rather than tutors, in helping students with their homework problems. In addition to supervised group work, students may discuss difficulties in their conceptual understanding of lecture and reading topics with their peers and the GSIs. There is no TALC on 08/25, 09/29, and 11/03. • Midterms (during lecture hours): o Thursday 09/29 in 131 Campbell Hall o Thursday 11/03 in 131 Campbell Hall • Final exam: Wednesday 12/14, 8-11 am, Campbell 131 (last names starting with J-Z) and Campbell 121 (last names starting with A-H). If you miss a midterm or the final exam without any notice you will receive zero credit for that portion of the course grade. Exams can be rescheduled in exceptional cases. If you miss the final exam for a good reason, your grade will be an Incomplete. • Office hours: - Nick Kern: Wednesday 2:30-3:30 pm, Campbell Hall 233 - Mariska Kriek: Tuesday 4:30-5:30 pm, Campbell Hall 355 - Michael Medford: Monday 12:00-1:00 pm, Campbell Hall 233 • Book: "An Introduction to Modern Astrophysics (2nd edition)" by Carroll & Ostlie (required) • Prerequisites: Physics 7A & 7B (7B may be taken concurrently).
    [Show full text]
  • Evolving Pulsation of the Slowly Rotating Magnetic $\Beta $ Cep Star
    MNRAS 000,1{14 (2019) Preprint 19 December 2019 Compiled using MNRAS LATEX style file v3.0 Evolving pulsation of the slowly rotating magnetic β Cep star ξ1 CMa G.A. Wade?1, A. Pigulski2, S. Begy1, M. Shultz3, G. Handler4, J. Sikora5, H. Neilson6, H. Cugier2, C. Erba3, A.F.J. Moffat7, B. Pablo8, A. Popowicz9, W. Weiss10, K. Zwintz11 1Dept. of Physics & Space Science, Royal Military College of Canada, PO Box 17000 Station Forces, Kingston, ON, Canada K7K 0C6 2Instytut Astronomiczny, Uniwersytet Wroc lawski, Kopernika 11, 51-622 Wroc law, Poland 3Dept. of Physics and Astronomy, University of Delaware, 217 Sharp Lab, Newark, DE 19716, USA 4Nicolaus Copernicus Astronomical Center, Bartycka 18, 00{716 Warszawa, Poland 5Dept. of Physics, Engineering Physics and Astronomy, Queen's University, 99 University Avenue, Kingston, ON K7L 3N6, Canada 6Dept. of Astronomy & Astrophysics, University of Toronto, 50 St. George Street, Toronto, ON M5S 3H4, Canada 7D´ept. de physique and Centre de Recherche en Astrophysique du Qu´ebec (CRAQ), Universit´ede Montr´eal, C.P. 6128, Succ. Centre-Ville, Montr´eal, QC H3C 3J7, Canada 8American Association of Variable Star Observers, 49 Bay State Road, Cambridge, MA 02138, USA 9Silesian University of Technology, Institute of Automatic Control, Akademicka 16, Gliwice, Poland 10Institut fur¨ Astrophysik, Universit¨at Wien, Turkenschanzstrasse¨ 17, A-1180 Wien, Austria 11Institut fur¨ Astro- und Teilchenphysik, Universit¨at Innsbruck, Technikerstrasse 25/8, A-6020 Innsbruck, Austria Accepted . Received ; in original form ABSTRACT Recent BRITE-Constellation space photometry of the slowly rotating, magnetic β Cep pulsator ξ1 CMa permits a new analysis of its pulsation properties.
    [Show full text]
  • Stellar Structure and Evolution
    STELLAR STRUCTURE AND EVOLUTION O.R. Pols Astronomical Institute Utrecht September 2011 Preface These lecture notes are intended for an advanced astrophysics course on Stellar Structure and Evolu- tion given at Utrecht University (NS-AP434M). Their goal is to provide an overview of the physics of stellar interiors and its application to the theory of stellar structure and evolution, at a level appro- priate for a third-year Bachelor student or beginning Master student in astronomy. To a large extent these notes draw on the classical textbook by Kippenhahn & Weigert (1990; see below), but leaving out unnecessary detail while incorporating recent astrophysical insights and up-to-date results. At the same time I have aimed to concentrate on physical insight rather than rigorous derivations, and to present the material in a logical order, following in part the very lucid but somewhat more basic textbook by Prialnik (2000). Finally, I have borrowed some ideas from the textbooks by Hansen, Kawaler & Trimble (2004), Salaris & Cassissi (2005) and the recent book by Maeder (2009). These lecture notes are evolving and I try to keep them up to date. If you find any errors or incon- sistencies, I would be grateful if you could notify me by email ([email protected]). Onno Pols Utrecht, September 2011 Literature C.J. Hansen, S.D. Kawaler & V. Trimble, Stellar Interiors, 2004, Springer-Verlag, ISBN 0-387- • 20089-4 (Hansen) R. Kippenhahn & A. Weigert, Stellar Structure and Evolution, 1990, Springer-Verlag, ISBN • 3-540-50211-4 (Kippenhahn; K&W) A. Maeder, Physics, Formation and Evolution of Rotating Stars, 2009, Springer-Verlag, ISBN • 978-3-540-76948-4 (Maeder) D.
    [Show full text]
  • Kasliwal Phd Thesis
    Bridging the Gap: Elusive Explosions in the Local Universe Thesis by Mansi M. Kasliwal Advisor Professor Shri R. Kulkarni In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy California Institute of Technology Pasadena, California 2011 (Defended April 26, 2011) ii c 2011 Mansi M. Kasliwal All rights Reserved iii Acknowledgements The first word that comes to my mind to describe my learning experience at Caltech is exhilarating. I have no words to thank my “Guru”, Professor Shri Kulkarni. Shri taught me the “Ps” necessary to Pursue the Profession of a Professor in astroPhysics. In addition to Passion & Perseverance, I am now Prepared for Papers, Proposals, Physics, Presentations, Politics, Priorities, oPPortunity, People-skills, Patience and PJs. Thanks to the awesomely fantastic Palomar Transient Factory team, especially Peter Nugent, Robert Quimby, Eran Ofek and Nick Law for sharing the pains and joys of getting a factory off-the-ground. Thanks to Brad Cenko and Richard Walters for the many hours spent taming the robot on mimir2:9. Thanks to Avishay Gal-Yam and Lars Bildsten for illuminating discussions on different observational and theoretical aspects of transients in the gap. The journey from idea to first light to a factory churning out thousands of transients has been so much fun that I would not trade this experience for any other. Thanks to Marten van Kerkwijk for supporting my “fishing in new waters” project with the Canada France Hawaii Telescope. Thanks to Dale Frail for supporting a “kissing frogs” radio program. Thanks to Sterl Phinney, Linqing Wen and Samaya Nissanke for great discussions on the challenge of finding the light in the gravitational sound wave.
    [Show full text]
  • ASTR 127 LECTURE NOTES JAN 13, 2015 1. Administrative Topics 1.1
    ASTR 127 LECTURE NOTES JAN 13, 2015 MICHAEL TOPPING 1. Administrative Topics 1.1. Overview of the discussion. This is the discussion for ASTR 127: Stellar Atmospheres, Interiors, and Evolution. 1.2. What to do during the discussion section? • Review of the lectures • Talk about relevant, but not required topics • Do example problems from the textbook, etc. • Go over previous homework • Do activities { Worksheets { Problem solving { Group work 1.3. Email Policy. Give me 24 hours to respond to emails, although I may respond earlier. After 24 hours feel free to email me again. Don't email me questions about what we did in lecture, because I don't go to the lectures. 1 2 MICHAEL TOPPING 1.4. Office Hours. • Thursdays 3pm - 4pm • Fridays 2pm-3pm • Or by appointment • Location: PAB 1-704A ASTR 127 LECTURE NOTES JAN 13, 2015 3 2. Timescales Timescales are important to know because they can help us determine the causes for different aspects of stellar evolution/structure. For example, If you calculate the lifetime of a star as if it were a ball of hydrogen on fire, you would get a lifetime much shorter than the age of the sun currently. Here are three different timescales in order of shortest to longest: 2.1. Dynamical Timescale. Start with the acceleration due to gravity: d2r GM (1) a = = − dt2 r2 From the chain rule, we know that: d2r d dr d 1 d (2) = v(r) = v = v2 dt2 dt dt dr 2 dr Therefore: d 2GM (3) v2 = − dr r2 Integrate both sides: Z Z 2GM (4) dv2 = − dr r2 With initial conditions of r = R, v = 0, the velocity becomes: 1 dr 1 1 2 (5) v = = 2GM( − ) dt r R 4 MICHAEL TOPPING Solving for the time interval dt we have: −dr (6) dt = 1 2 1 1 2GM( r − R ) Rearranging this and integrating we then get: 0 R 1 Z dr (7) t = ( ) 2 ff R 1 2GM 2 R (1 − r ) r dr With the substitution x = R and dx = R we obtain: 3 0 R 1 Z dx (8) t = ( ) 2 ff 1 1 2GM 2 1 (1 − x ) π The integral evaluates to 2 , then we have a final answer of: 1 3π 2 (9) t = ff 32Gρ This result (Equation 9) is the free fall timescale, sometimes called the dynamical timescale.
    [Show full text]
  • Numerical Simulations of Mass Transfer in Close and Contact Binaries Using Bipolytropes
    Louisiana State University LSU Digital Commons LSU Doctoral Dissertations Graduate School 2017 Numerical Simulations of Mass Transfer in Close and Contact Binaries Using Bipolytropes Kundan Vaman Kadam Louisiana State University and Agricultural and Mechanical College, [email protected] Follow this and additional works at: https://digitalcommons.lsu.edu/gradschool_dissertations Part of the Physical Sciences and Mathematics Commons Recommended Citation Kadam, Kundan Vaman, "Numerical Simulations of Mass Transfer in Close and Contact Binaries Using Bipolytropes" (2017). LSU Doctoral Dissertations. 4325. https://digitalcommons.lsu.edu/gradschool_dissertations/4325 This Dissertation is brought to you for free and open access by the Graduate School at LSU Digital Commons. It has been accepted for inclusion in LSU Doctoral Dissertations by an authorized graduate school editor of LSU Digital Commons. For more information, please [email protected]. NUMERICAL SIMULATIONS OF MASS TRANSFER IN CLOSE AND CONTACT BINARIES USING BIPOLYTROPES A Dissertation Submitted to the Graduate Faculty of the Louisiana State University and Agricultural and Mechanical College in partial fulfillment of the requirements for the degree of Doctor of Philosophy in The Department of Physics and Astronomy by Kundan Vaman Kadam B.S. in Physics, University of Mumbai, 2005 M.S., University of Mumbai, 2007 May 2017 Acknowledgments I am grateful to my advisor, Geoff Clayton, for supporting me wholeheartedly throughout the graduate school and keeping me motivated even when things didn't go as planned. I'd like to thank Hartmut Kaiser for allowing me to work on my dissertation through the NSF CREATIV grant. I am thankful to Juhan Frank for insightful discussions on the physics of binary systems and their simulations.
    [Show full text]
  • 8.901 Lecture Notes Astrophysics I, Spring 2019
    8.901 Lecture Notes Astrophysics I, Spring 2019 I.J.M. Crossfield (with S. Hughes and E. Mills)* MIT 6th February, 2019 – 15th May, 2019 Contents 1 Introduction to Astronomy and Astrophysics 6 2 The Two-Body Problem and Kepler’s Laws 10 3 The Two-Body Problem, Continued 14 4 Binary Systems 21 4.1 Empirical Facts about binaries................... 21 4.2 Parameterization of Binary Orbits................. 21 4.3 Binary Observations......................... 22 5 Gravitational Waves 25 5.1 Gravitational Radiation........................ 27 5.2 Practical Effects............................ 28 6 Radiation 30 6.1 Radiation from Space......................... 30 6.2 Conservation of Specific Intensity................. 33 6.3 Blackbody Radiation......................... 36 6.4 Radiation, Luminosity, and Temperature............. 37 7 Radiative Transfer 38 7.1 The Equation of Radiative Transfer................. 38 7.2 Solutions to the Radiative Transfer Equation........... 40 7.3 Kirchhoff’s Laws........................... 41 8 Stellar Classification, Spectra, and Some Thermodynamics 44 8.1 Classification.............................. 44 8.2 Thermodynamic Equilibrium.................... 46 8.3 Local Thermodynamic Equilibrium................ 47 8.4 Stellar Lines and Atomic Populations............... 48 *[email protected] 1 Contents 8.5 The Saha Equation.......................... 48 9 Stellar Atmospheres 54 9.1 The Plane-parallel Approximation................. 54 9.2 Gray Atmosphere........................... 56 9.3 The Eddington Approximation................... 59 9.4 Frequency-Dependent Quantities.................. 61 9.5 Opacities................................ 62 10 Timescales in Stellar Interiors 67 10.1 Photon collisions with matter.................... 67 10.2 Gravity and the free-fall timescale................. 68 10.3 The sound-crossing time....................... 71 10.4 Radiation transport.......................... 72 10.5 Thermal (Kelvin-Helmholtz) timescale............... 72 10.6 Nuclear timescale..........................
    [Show full text]