Things I don't understand about Vlasov-Poisson numerics
|
|
- Isaac Chambers
- 5 years ago
- Views:
Transcription
1 Things I don't understand about Vlasov-Poisson numerics
2 1. Modeling stellar systems 2. Interpolation in simulations of cold dark matter 3. Orbit integration in N-body systems 4. How accurately should we integrate?
3 1. Modeling stellar systems
4 1. Modeling stellar systems the problem: model a stellar system such as an E0 galaxy the most we can hope to determine is ν(r,vlos), the surface brightness at projected position R in a unit interval of line-ofsight velocity vlos given a spherical, stationary, non-rotating stellar system with known ν(r,vlos), find the gravitational potential Φ(r) and the distribution function (df) f(r,vr,v ) Jeans s theorem tells us that the df can only depend on energy E = v 2 /2 + Φ(r) and angular momentum L=rv, that is, f(r,vr,v ) = f(e,l)
5 given a spherical, stationary, non-rotating stellar system with known ν(r,vlos), find the gravitational potential Φ(r) and the df f(e,l) why this is difficult: given Φ(r) it s straightforward to determine f(e,l) but in general you can t determine a function of two variables and a function of one variable [ν(r,vlos) and Φ(r)] from a single function of two variables f(e,l)
6 given a spherical, stationary, non-rotating stellar system with known ν(r,vlos), find the gravitational potential Φ(r) and the df f(e,l) in general you can t determine a function of two variables and a function of one variable [ν(r,vlos) and Φ(r)] from a single function of two variables f(e,l) Dejonghe & Merritt (1992): assume Φ(r) and determine f(e,l) if Φ(r) is wrong, it s likely to give f(e,l) < 0 in some places, which is unphysical therefore Φ(r) is determined through a nonlinear, implicit, integral inequality other difficulties: the integral equation for f(e,l) given Φ(r) is ill-conditioned so noise is amplified determination of ν(r,vlos) is subject to systematic errors because of template mismatches and continuum subtraction
7 given a spherical, stationary, non-rotating stellar system with known ν(r,vlos), find the gravitational potential Φ(r) and the df f(e,l) how this is done in practice: the gold standard is Schwarzschild s (1978) method: assume some potential Φ(r) divide phase space into bins of volume Δwj = ΔxjΔvj integrate N orbits from a wide variety of initial conditions and track the fraction of time orbit i spends in phase-space bin j, fij assign luminosity λi to orbit i and compute the expected surface brightness νλ(r,vlos). Calculate χ 2 = [ν(rk,vlos,k) - νλ(rk,vlos,k)] 2 and adjust the luminosities λi to minimize χ 2 df is then f(xj,vj)= i λifij/δwj repeat for other potentials and find the potential that minimizes χ 2
8 given a spherical, stationary, non-rotating stellar system with known ν(r,vlos), find the gravitational potential Φ(r) and the df f(e,l) how this is done in practice: the gold standard is Schwarzschild s (1978) method problems with Schwarzschild s method: the problem is under-determined and ill-conditioned so you must regularize, e.g., find maximum-entropy solution, apply smoothness constraints to the df, etc. the number of degrees of freedom is impossible to estimate the method is noisy because orbits have square-root singularities at their turning points in estimating the parameters of the potential one should marginalize over the df (the λi) but all existing methods just find the maximumlikelihood λi (Magorrian 2007)
9 given a spherical, stationary, non-rotating stellar system with known ν(r,vlos), find the gravitational potential Φ(r) and the df f(e,l) how this is done in practice: the gold standard is Schwarzschild s (1978) method there are many problems with Schwarzschild s method despite all these concerns, Schwarzschild s method, and simpler methods validated by Schwarzschild s method, is the basis for almost all we know about the dynamics of hot stellar systems existence and masses of black holes in the centers of galaxies distribution of dark matter in elliptical galaxies mass-to-light ratios of elliptical galaxies, tilt of the fundamental plane cores vs. cusps in dwarf spheroidal galaxies
10 N test particles orbit a point mass M in a steady state. We have three-dimensional position and velocity vectors for all N particles. What is the best estimate of M? 1. Virial theorem: A simple toy problem drawbacks: can t be generalized to potentials depending on several parameters inefficient use of data if there is a range of radii no error estimates
11 A simple toy problem 2. Projected mass estimator: drawbacks: can t be generalized to case where only line-of-sight velocities are available, without introducing dependence on eccentricity distribution can t be generalized to potentials depending on several parameters
12 A simple toy problem 3. Orbital roulette (Beloborodov & Levin 2004) drawbacks: can t be generalized to case where only line-of-sight velocities are available fails for zero eccentricity orbits result depends on the choice of statistical test
13 4. What would Bayes do? A simple toy problem
14 A simple toy problem
15 4. What would Bayes do? A simple toy problem
16 A simple toy problem Test on 1000 realizations of a 10- particle system; uniform distribution in e 2 and central mass M=1
17 A simple toy problem Test on 1000 realizations of a 10- particle system with central mass M=1 uniform in e 2
18 now apply the same method to a general integrable potential...
19
20 Bayesian approach to determining potential parameters fails, because there is a unique orbit consistent with every star so all potentials have the same likelihood
21 Bayesian approach to estimating potential parameters the Bayesian approach works well for power-law potentials Φ(r) = cr -k if it is applied to variables with the same dimension as c, i.e. variables of the form v 2 r k and vr 2 r k for a general integrable potential and a non-smooth distribution function f(j), a Bayesian approach to determining the gravitational potential fails bigger observational errors give better determinations of the potential parameters! (Magorrian 2013) the results from Bayesian analyses depend on arbitrary assumptions about the smoothness of the distribution function
22 Bayesian approach to estimating potential parameters the Bayesian approach works well for power-law potentials Φ(r) = cr -k if it is applied to variables with the same dimension as c, what i.e. variables are we of the really form doing v 2 r k and when vr 2 r k we estimate gravitational for a general integrable parameters potential of and galaxies a non-smooth using distribution function f(j), a Bayesian approach to determining the Schwarzschild s method? gravitational potential fails bigger observational errors give better determinations of the potential parameters! (Magorrian 2013) the results from Bayesian analyses depend on arbitrary assumptions about the smoothness of the distribution function
23 2. N-body simulations of cold dark matter
24 standard cosmological model involves cold dark matter: at high redshift, dark matter occupies a 3-dimensional manifold or sheet in 6-dimensional phase space, with near-zero comoving velocity as structure develops this sheet is stretched and folded but never torn, and remains 3-dimensional N-body simulations follow a Monte Carlo sampling of the mass distribution in this sheet particles in a simulation are vertices in a 3-dimensional mesh that describes the sheet how can we best characterize the properties of the sheet, i.e. how should we interpolate between particles? Schulz et al. (2013)
25 N-body simulations follow a Monte Carlo sampling of the mass distribution in this sheet particles in a simulation are vertices in a 3-dimensional mesh that describes the sheet how can we best characterize the properties of the sheet, i.e. how should we interpolate between particles? Simplest interpolation schemes use only position information: traditionally these were developed to generate gravitational potential from Poisson s equation. Two approaches: replace each particle with a smooth density kernel (as in softening and SPH) divide space into a grid and assign mass of particle to nearby grid points (NGP, CIC, etc.) Voroni and Delaunay tesselations (Bernardeau & van de Weygaert 1996) k nearest neighbors recursive binary splitting (Ascasibar & Binney 2005) More sophisticated schemes use both position and velocity information (interpolation in phase space)
26 Phase-space interpolation schemes for cold dark matter interpolate onto a grid in phase space extremely memory-intensive: grid with N elements per dimension needs N 6 grid points Yoshikawa et al. (2013): direct solution of Vlasov-Poisson equations using 646 phase-space grid (run on 64 nodes, 32 GB RAM per node, 180 s per timestep) kernel density estimation needs a metric for phase space Delaunay tesselation (Arad et al. 2004) expensive in time and memory k nearest neighbors needs a metric recursive binary splitting (Maciejewski et al. 2009)
27 Phase-space interpolation schemes for cold dark matter interpolate onto a grid in phase space kernel density estimation Delaunay tesselation k nearest neighbors recursive binary splitting Lagrangian tesselation (Abel et al. 2012, Shandarin et al. 2012) tesselate the phase space at t=0 and use an N-body simulation to follow the nodes of the tesselation
28 Phase-space interpolation schemes for cold dark matter particle-based Lagrangian tesselation Shandarin et al. (2012)
29 Phase-space interpolation schemes for cold dark matter Abel et al. (2012) grid-based density estimate Lagrangian tesselation density estimate
30 Phase-space interpolation schemes for cold dark matter interpolate onto a grid in phase space kernel density estimation Delaunay tesselation k nearest neighbors recursive binary splitting Lagrangian tesselation (Shandarin et al. 2012) geodesic deviation (Vogelsberger et al. 2008)
31 Phase-space interpolation schemes for cold dark matter
32 Phase-space interpolation schemes for cold dark matter interpolate onto a grid in phase space kernel density estimation Delaunay tesselation k nearest neighbors recursive binary splitting Lagrangian tesselation (Shandarin et al. 2012) geodesic deviation (Vogelsberger et al. 2008) symplectic interpolation
33 Symplectic interpolation
34 Symplectic interpolation assume we tesselate 3-space with tetrahedra and use linear or quadratic interpolation within each tetrahedron: linear interpolation quadratic interpolation number of available interpolation constants number of constraints from symplectic condition 3 30 number of constraints from the 4 vertices of the tetrahedron number of constraints from adjacent tetrahedra -?
35 Symplectic interpolation assume we tesselate 3-space with tetrahedra and use linear or quadratic interpolation within each tetrahedron: is symplectic interpolation a useful tool for improving the resolution linear interpolation and power quadratic of large interpolation cosmological N-body simulations? number of available interpolation constants number of constraints from symplectic condition 3 30 number of constraints from the 4 vertices of the tetrahedron number of constraints from adjacent tetrahedra -?
36 3. Orbit integration in N-body systems
37 Orbit integration in N-body systems
38 Orbit integration in N-body systems standard approach is Leapfrog or Verlet method advantages of leapfrog time-reversible symplectic no memory needed for intermediate steps disadvantages of leapfrog only second-order accurate not suitable for collisional simulations not symplectic if variable timestep Consider an N-body system with stars of mass m and radius R. Typical velocity is v 2 =GNm/R. Crossing or dynamical time is tcr = R/v. Typical inter-particle separation is r = R/N 1/3. Leapfrog fails if the timestep h exceeds time required to cross r, h > r/v ~ tcr/n 1/3 as N increases and the potential becomes smoother, the maximum timestep must shrink
39 Impulse integrator with Hanno Rein (IAS and Toronto)
40 Impulse-approximation integrator symplectic second-order accurate (like leapfrog) no memory needed for intermediate steps no restriction to timesteps less than inter-particle travel time disadvantages: algebraic formulas are more complicated, so slower to evaluate implicit with Hanno Rein
41 rel. energy error rel. energy error e-05 1e-06 1e-07 1e-08 1e e-05 1e-06 1e-07 1e-08 1e-09 N = 100,000 Impulse-approximation integrator time [crossing times] N = 10,000,000 2nd order impulse integrator leapfrog integratorr time [crossing times] 2nd order impulse integrator leapfrog integratorr with Hanno Rein ape velocity neareast star N fixed particles in a uniform sphere, plus one test particle on circular orbit 1000 energy of test particle should be conserved 100 errors in leapfrog get bigger as N increases errors in impulse integrator get smaller as N increases N = 10,000,000 2nd order impulse integrator leapfrog integrator Test particle velocity
42 rel. energy error rel. energy error e-05 1e-06 1e-07 1e-08 1e e-05 1e-06 1e-07 1e-08 1e-09 N = 100,000 Impulse-approximation integrator time [crossing times] N = 10,000,000 2nd order impulse integrator leapfrog integratorr time [crossing times] 2nd order impulse integrator leapfrog integratorr with Hanno Rein ape velocity neareast star why does the error saturate? N = 10,000,000 why does it saturate at the same timestep for both methods? 1000 what determines the saturation level? 100 2nd order impulse integrator leapfrog integrator Test particle velocity
43 rel. energy error rel. energy error e-05 1e-06 1e-07 1e-08 1e e-05 1e-06 1e-07 1e-08 1e-09 N = 100,000 Impulse-approximation integrator time [crossing times] 2nd order impulse integrator leapfrog integratorr N = 10,000, can the impulse integrator be developed into a useful tool for N-body integrations? time [crossing times] 2nd order impulse integrator leapfrog integratorr with Hanno Rein ape velocity neareast star why does the error saturate? N = 10,000,000 why does it saturate at the same timestep for both methods? 1000 what determines the saturation level? 100 2nd order impulse integrator leapfrog integrator Test particle velocity
44 4. How do we know if an N-body integration is correct?
45 How do we know if an N-body integration is correct? small errors in an N-body integration grow exponentially (Miller 1964, 1971; Hayli 1970) e-folding or Liapunov time ~ crossing time tcr (Heggie 1988, Goodman et al. 1992). Not related to the relaxation time and independent of N therefore detailed final positions and velocities of particles in a simulation are not reproducible this has not deterred subsequent investigations from being carried out in an optimistic spirit (Aarseth & Lecar 1975) a pessimist might ask, why should we believe the results?
46 How do we know if an N-body integration is correct? small errors in an N-body integration grow exponentially (Miller 1964, 1971; Hayli 1970) e-folding or Liapunov time ~ crossing time tcr (Heggie 1988, Goodman et al. 1992). Not related to the relaxation time and independent of N therefore detailed final positions and velocities of particles in a simulation are not reproducible this has not deterred subsequent investigations from being carried out in an optimistic spirit (Aarseth & Lecar 1975) a pessimist might ask, why should we believe the results? forward error analysis: measures the divergence of the numerical solution from the exact solution with the same initial conditions backward error analysis: asks if the numerical solution is the same as the exact solution for a slightly different physical system
47 When should we believe the results of an N-body simulation? never because of the growth of computational errors N-body calculations in stellar dynamics are of limited usefulness (Miller 1964) - forward error analysis always, so long as the integrator is symplectic the argument is that then we are following a Hamiltonian system that is close to the real one - backward error analysis when the integration is time-reversible this is easy if it uses a time-reversible integrator and roundoff is negligible if it conserves energy and angular momentum this is very crude if the result looks OK implies we should use bigger and bigger timesteps until something goes grossly wrong if the result is reproduced in a statistical sample of simulations maybe they give consistent results only because they are all equally inaccurate (Heggie 1991) if the result looks the same when the integration is repeated with twice the particles, half the timestep, etc.
48 Shadowing (Quinlan & Tremaine 1992) forward error analysis: measures the divergence of the numerical solution from the exact solution with the same initial conditions backward error analysis: asks if the numerical solution is the same as the exact solution for a slightly different physical system shadowing: asks if the numerical solution is the same as an exact solution with slightly different initial conditions
49 Shadowing (Quinlan & Tremaine 1992) shadowing: asks if the numerical solution is the same as an exact solution with slightly different initial conditions a shadow orbit is an exact solution that stays within a small distance δ of the numerical solution over the integration time 0 < t < T existence of shadow orbits is not a rigorous proof that the N-body integration is realistic since the shadow orbits may have peculiar initial conditions finding shadow orbits is numerically challenging. Best approach is to use shooting: start the shadow in the contracting manifold of phase space at t=0 and integrate forward; start another one in the expanding manifold at t=t and integrate backward, then match in the middle works better in chaotic dynamical systems than integrable systems: it s easier to shadow ẍ = x than ẍ = const typically the existence of a shadow orbit terminates abruptly at a glitch ; these are real, not numerical artifacts
50 t=20.1 shadow exact t=3.25 numerical start, t=0 N=100 Plummer model, 99 particles fixed no softening numerical orbit followed with a standard Aarsethtype N-body integrator shadow distance in position δ = typical shadow times T ~ 10 and get longer if accuracy of numerical orbit is improved or softening is added Quinlan & Tremaine (1992
51 Hayes (2003): shadowed a wide variety of N-body integrations using leapfrog non-symplectic integrators are shdowed less well glitches appear to be a Poisson process the half-life for shadowing is crossing time softening timestep system size shadow duration is insensitive to N softening is needed to obtain long-lasting shadow orbits with any practical integration algorithm shadowing for many crossing times requires that particles travel a small fraction of the softening length in one timestep
52 Hayes (2003): shadowed a wide variety of N-body integrations using leapfrog non-symplectic integrators are shdowed less well glitches appear to be a Poisson process the half-life for shadowing is crossing time softening timestep system size how badly are we allowed to integrate? shadow duration is insensitive to N softening is needed to obtain long-lasting shadow orbits with any practical integration algorithm shadowing for many crossing times requires that particles travel a small fraction of the softening length in one timestep
High performance computing and numerical modeling
High performance computing and numerical modeling Volker Springel Plan for my lectures Lecture 1: Collisional and collisionless N-body dynamics Lecture 2: Gravitational force calculation Lecture 3: Basic
More informationDynamics of Galaxies: Basics. Frontiers in Numerical Gravitational Astrophysics July 3, 2008
Dynamics of Galaxies: Basics Frontiers in Numerical Gravitational Astrophysics July 3, 2008 c = G = 1 h = 0 Outline (a) Galaxies as physical systems. (b) The Collisionless Boltzmann Equation. (c) N-body
More informationDynamical modelling of galaxies. Eugene Vasiliev
Dynamical modelling of galaxies Eugene Vasiliev DAGAL summer school, Marseille, 2014 Plan Overview of dynamical modelling problem Theoretical and observational approaches The Schwarzschild orbit superposition
More informationThe Los Cabos Lectures
January 2009 The Los Cabos Lectures Dark Matter Halos: 3 Simon White Max Planck Institute for Astrophysics Shapes of halo equidensity surfaces Group Jing & Suto 2002 Galaxy δ 100 2500 6250 Shapes of halo
More informationGalactic Dynamics. James Binney and Scott Tremaine SECOND EDITION PRINCETON UNIVERSITY,PRESS PRINCETON AND OXFORD
Galactic Dynamics SECOND EDITION James Binney and Scott Tremaine PRINCETON UNIVERSITY,PRESS PRINCETON AND OXFORD Preface 1 Introduction - 1 1.1 An overview of the observations 5 Stars 5 The Galaxy 11 Other
More informationThe Worst-Case Scenario. The Final Parsec Problem. and. Milos Milosavljevic. California Institute of Technology. Collaborator: David Merritt
The Final Parsec Problem and The Worst-Case Scenario Milos Milosavljevic California Institute of Technology Collaborator: David Merritt NSF AST 00-71099 NASA NAG5-6037, NAG5-9046 Sherman Fairchild Foundation
More information3 The lives of galaxies
Discovering Astronomy : Galaxies and Cosmology 24 3 The lives of galaxies In this section, we look at how galaxies formed and evolved, and likewise how the large scale pattern of galaxies formed. But before
More informationOverview spherical accretion
Spherical accretion - AGN generates energy by accretion, i.e., capture of ambient matter in gravitational potential of black hole -Potential energy can be released as radiation, and (some of) this can
More informationEnsenada March 2008 Galaxy halos at (very) high resolution
Ensenada March 2008 Galaxy halos at (very) high resolution Simon White & Mark Vogelsberger Max Planck Institute for Astrophysics Milky Way halo seen in DM annihilation radiation Aquarius simulation: N200
More informationBlack Hole Mergers at Galactic. The Final Parsec: Supermassive. Centers. Milos Milosavljevic. California Institute of Technology
The Final Parsec: Supermassive Black Hole Mergers at Galactic Centers Milos Milosavljevic California Institute of Technology MBH Binaries Form in Galaxy Mergers Borne et al 2000 by transferring binding
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 informationThe fine-scale structure of dark matter halos
COSMO11, Porto, August 2011 The fine-scale structure of dark matter halos Simon White Max-Planck-Institute for Astrophysics COSMO11, Porto, August 2011 Mark Vogelsberger The fine-scale structure of dark
More informationTreecodes for Cosmology Thomas Quinn University of Washington N-Body Shop
Treecodes for Cosmology Thomas Quinn University of Washington N-Body Shop Outline Motivation Multipole Expansions Tree Algorithms Periodic Boundaries Time integration Gravitational Softening SPH Parallel
More informationAdvection Dominated Accretion Flows. A Toy Disk Model. Bohdan P a c z y ń s k i
ACTA ASTRONOMICA Vol. 48 (1998) pp. 667 676 Advection Dominated Accretion Flows. A Toy Disk Model by Bohdan P a c z y ń s k i Princeton University Observatory, Princeton, NJ 8544-11, USA e-mail: bp@astro.princeton.edu
More informationNumerical Investigations of a New N-Body Simulation Method
International Journal of Astronomy and Astrophysics, 2012, 2, 119-124 http://dx.doi.org/10.4236/ijaa.2012.23016 Published Online September 2012 (http://www.scirp.org/journal/ijaa) Numerical Investigations
More informationFrom An Apple To Black Holes Gravity in General Relativity
From An Apple To Black Holes Gravity in General Relativity Gravity as Geometry Central Idea of General Relativity Gravitational field vs magnetic field Uniqueness of trajectory in space and time Uniqueness
More informationGalaxy interaction and transformation
Galaxy interaction and transformation Houjun Mo April 13, 2004 A lot of mergers expected in hierarchical models. The main issues: The phenomena of galaxy interaction: tidal tails, mergers, starbursts When
More informationWhat can M2M do for Milky Way Models?
What can M2M do for Milky Way Models? Ortwin Gerhard Max-Planck-Institut für Extraterrestrische Physik, Garching gerhard@mpe.mpg.de Thanks to F. de Lorenzi, V. Debattista, P. Das, L. Morganti I. Made-to-Measure
More informationScott A. Hughes, MIT SSI, 28 July The basic concepts and properties of black holes in general relativity
The basic concepts and properties of black holes in general relativity For the duration of this talk ħ=0 Heuristic idea: object with gravity so strong that light cannot escape Key concepts from general
More informationAstronomy 102: Stars and Galaxies Final Exam Review Problems Revision 2
Astronomy 102: Stars and Galaxies Final Exam Review Problems Revision 2 Multiple Choice Questions: The first eight questions are multiple choice. Except where explicitly noted, only one answer is correct
More informationAstronomy 421. Lecture 24: Black Holes
Astronomy 421 Lecture 24: Black Holes 1 Outline General Relativity Equivalence Principle and its Consequences The Schwarzschild Metric The Kerr Metric for rotating black holes Black holes Black hole candidates
More informationOrbits, Integrals, and Chaos
Chapter 7 Orbits, Integrals, and Chaos In n space dimensions, some orbits can be formally decomposed into n independent periodic motions. These are the regular orbits; they may be represented as winding
More informationComponents of Galaxies: Dark Matter
Components of Galaxies: Dark Matter Dark Matter: Any Form of matter whose existence is inferred solely through its gravitational effects. -B&T, pg 590 Nature of Major Component of Universe Galaxy Formation
More informationConnecting observations to simulations arxiv: Joe Wolf (UC Irvine)
Connecting observations to simulations arxiv: 0908.2995 Joe Wolf (UC Irvine) University of Maryland December 8 th, 2009 Team Irvine: Greg Martinez James Bullock Manoj Kaplinghat Erik Tollerud Quinn Minor
More informationN-Body Methods. Chapter Numerical Stellar Dynamics Monte carlo methods Representing f (r;v;t)
Chapter 9 N-Body Methods The Collisionless Boltzmann and Poisson equations can be effectively solved with a Monte-Carlo technique which represents the distribution function f (r;v;t) as a set of N bodies.
More informationChapter 18 Lecture. The Cosmic Perspective Seventh Edition. The Bizarre Stellar Graveyard Pearson Education, Inc.
Chapter 18 Lecture The Cosmic Perspective Seventh Edition The Bizarre Stellar Graveyard The Bizarre Stellar Graveyard 18.1 White Dwarfs Our goals for learning: What is a white dwarf? What can happen to
More informationChapter 14: The Bizarre Stellar Graveyard
Lecture Outline Chapter 14: The Bizarre Stellar Graveyard 14.1 White Dwarfs Our goals for learning: What is a white dwarf? What can happen to a white dwarf in a close binary system? What is a white dwarf?
More informationAM 205 Final Project The N-Body Problem
AM 205 Final Project The N-Body Problem Leah Birch Elizabeth Finn Karen Yu December 14, 2012 Abstract The N-Body Problem can be solved using a variety of numeric integrators. Newton s Law of Universal
More informationWhite dwarfs are the remaining cores of dead stars. Electron degeneracy pressure supports them against the crush of gravity. The White Dwarf Limit
The Bizarre Stellar Graveyard Chapter 18 Lecture The Cosmic Perspective 18.1 White Dwarfs Our goals for learning: What is a white dwarf? What can happen to a white dwarf in a close binary system? Seventh
More informationOverview of Dynamical Modeling. Glenn van de Ven
Overview of Dynamical Modeling Glenn van de Ven glenn@mpia.de 1 Why dynamical modeling? -- mass total mass stellar systems key is to their evolution compare luminous mass: constrain DM and/or IMF DM radial
More informationAS1001:Extra-Galactic Astronomy
AS1001:Extra-Galactic Astronomy Lecture 5: Dark Matter Simon Driver Theatre B spd3@st-andrews.ac.uk http://www-star.st-and.ac.uk/~spd3 Stars and Gas in Galaxies Stars form from gas in galaxy In the high-density
More informationCHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,
More informationDwarf Galaxy Dispersion Profile Calculations Using a Simplified MOND External Field Effect
Dwarf Galaxy Dispersion Profile Calculations Using a Simplified MOND External Field Effect On the Shoulders of Giants Workshop Case Western Reserve University June 7, 2017 Stephen Alexander Physics Department
More informationPressure corrected SPH for fluid animation
Pressure corrected SPH for fluid animation Kai Bao, Hui Zhang, Lili Zheng and Enhua Wu Analyzed by Po-Ram Kim 2 March 2010 Abstract We present pressure scheme for the SPH for fluid animation In conventional
More information11/1/16. Important Stuff (Section 001: 9:45 am) Important Stuff (Section 002, 1:00 pm) 14.1 White Dwarfs. Chapter 14: The Bizarre Stellar Graveyard
Important Stuff (Section 001: 9:45 am) The Second Midterm is Thursday, November 10 The Second Midterm will be given in a different room: Willey 175 Bring 2 pencils and a photo-id. In accordance with the
More information10/25/2010. Stars, Galaxies & the Universe Announcements. Stars, Galaxies & the Universe Lecture Outline. Reading Quiz #9 Wednesday (10/27)
Stars, Galaxies & the Universe Announcements Reading Quiz #9 Wednesday (10/27) HW#8 in ICON due Friday (10/29) by 5 pm - available Wednesday 1 Stars, Galaxies & the Universe Lecture Outline 1. Black Holes
More informationNMAGIC Made-to-Measure Modeling of Elliptical Galaxies NMAGIC - 2 M2M
NMAGIC Made-to-Measure Modeling of Elliptical Galaxies Ortwin Gerhard, MPE, Garching gerhard@mpe.mpg.de I. Made-to-measure dynamical modeling with NMAGIC II Elliptical galaxy halo kinematics with planetary
More informationThe motions of stars in the Galaxy
The motions of stars in the Galaxy The stars in the Galaxy define various components, that do not only differ in their spatial distribution but also in their kinematics. The dominant motion of stars (and
More information11/1/17. Important Stuff (Section 001: 9:45 am) Important Stuff (Section 002, 1:00 pm) 14.1 White Dwarfs. Chapter 14: The Bizarre Stellar Graveyard
11/1/17 Important Stuff (Section 001: 9:45 am) The Second Midterm is Thursday, November 9 The Second Midterm will be given in a different room: Willey 175 Bring 2 pencils and a photo-id. In accordance
More informationChapter 18 The Bizarre Stellar Graveyard
Chapter 18 The Bizarre Stellar Graveyard 18.1 White Dwarfs Our goals for learning What is a white dwarf? What can happen to a white dwarf in a close binary system? What is a white dwarf? White Dwarfs White
More informationAstr 5465 Feb. 6, 2018 Today s Topics
Astr 5465 Feb. 6, 2018 Today s Topics Stars: Binary Stars Determination of Stellar Properties via Binary Stars Classification of Binary Stars Visual Binaries Both stars visible Only one star visible Spectroscopic
More informationChapter 18 The Bizarre Stellar Graveyard. White Dwarfs. What is a white dwarf? Size of a White Dwarf White Dwarfs
Chapter 18 The Bizarre Stellar Graveyard 18.1 White Dwarfs Our goals for learning What is a white dwarf? What can happen to a white dwarf in a close binary system? What is a white dwarf? White Dwarfs White
More informationSelected Topics in Mathematical Physics Prof. Balakrishnan Department of Physics Indian Institute of Technology, Madras
Selected Topics in Mathematical Physics Prof. Balakrishnan Department of Physics Indian Institute of Technology, Madras Module - 11 Lecture - 29 Green Function for (Del Squared plus K Squared): Nonrelativistic
More informationIntroduction: Special Relativity
Introduction: Special Relativity Observation: The speed c e.g., the speed of light is the same in all coordinate systems i.e. an object moving with c in S will be moving with c in S Therefore: If " r!
More informationAstronomy 8824: Numerical Methods Notes 2 Ordinary Differential Equations
Astronomy 8824: Numerical Methods Notes 2 Ordinary Differential Equations Reading: Numerical Recipes, chapter on Integration of Ordinary Differential Equations (which is ch. 15, 16, or 17 depending on
More informationStellar Dynamics and Structure of Galaxies
Stellar Dynamics and Structure of Galaxies Gerry Gilmore H47 email: gil@ast.cam.ac.uk Lectures: Monday 12:10-13:00 Wednesday 11:15-12:05 Friday 12:10-13:00 Books: Binney & Tremaine Galactic Dynamics Princeton
More informationModule 8: The Cosmos in Motion. UNC-TFA H.S. Astronomy Collaboration, Copyright 2011
Objectives/Key Points Module 8: The Cosmos in Motion UNC-TFA H.S. Astronomy Collaboration, Copyright 2011 1. Differentiate between classical motions across space (peculiar velocities) and cosmological
More informationMade-to-Measure (M2M) Dynamical modelling of external galaxies and the Milky Way
Made-to-Measure (M2M) Dynamical modelling of external galaxies and the Milky Way Shude Mao (NAOC/Manchester) Oct. 31, 2012@KIAS Collaborators: Richard Long, Yougang Wang, Juntai Shen Motivation Important
More informationPhysics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion
Physics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion Peter Young (Dated: April 14, 2009) I. INTRODUCTION One frequently obtains detailed dynamical
More informationClusters: Context and Background
Clusters: Context and Background We re about to embark on a subject rather different from what we ve treated before, so it is useful to step back and think again about what we want to accomplish in this
More informationBrad Gibson Centre for Astrophysics & Supercomputing Swinburne University
N-body Simulations: Tree and Mesh Approaches Brad Gibson Centre for Astrophysics & Supercomputing Swinburne University Demystifying the Jargon Trees Collisionless Boltzmann Equation Symplectic Time Integration
More informationA A + B. ra + A + 1. We now want to solve the Einstein equations in the following cases:
Lecture 29: Cosmology Cosmology Reading: Weinberg, Ch A metric tensor appropriate to infalling matter In general (see, eg, Weinberg, Ch ) we may write a spherically symmetric, time-dependent metric in
More informationComputer Simulation Using Particles
Computer Simulation Using Particles R W Hockney Emeritus Professor, University ofreading and J W Eastwood Culham Laboratory, Abingdon Institute of Physics Publishing Bristol and Philadelphia CONTENTS Foreword
More informationEinstein s Gravity. Understanding space-time and the gravitational effects of mass
Einstein s Gravity Understanding space-time and the gravitational effects of mass Albert Einstein (1879-1955) One of the iconic figures of the 20 th century, Einstein revolutionized our understanding of
More informationFor β = 0, we have from our previous equations: d ln ν d ln r + d ln v2 r
Isotropic models Now simplify the equation, and assume β = 0, or isotropic velocity dispersions - i.e., they are the same in all directions. Below we show how the observed light distribution and observed
More information2 Frederic A. Rasio. 2. Summary of Recent Results
2 Frederic A. Rasio dynamical evolution of a self-gravitating spherical cluster of N point masses whose orbits in the cluster are specified by an energy E and angular momentum J, with perturbations E and
More informationPotential/density pairs and Gauss s law
Potential/density pairs and Gauss s law We showed last time that the motion of a particle in a cluster will evolve gradually, on the relaxation time scale. This time, however, is much longer than the typical
More informationReview Questions for the new topics that will be on the Final Exam
Review Questions for the new topics that will be on the Final Exam Be sure to review the lecture-tutorials and the material we covered on the first three exams. How does speed differ from velocity? Give
More informationAP Physics QUIZ Gravitation
AP Physics QUIZ Gravitation Name: 1. If F1 is the magnitude of the force exerted by the Earth on a satellite in orbit about the Earth and F2 is the magnitude of the force exerted by the satellite on the
More information3 The Friedmann-Robertson-Walker metric
3 The Friedmann-Robertson-Walker metric 3.1 Three dimensions The most general isotropic and homogeneous metric in three dimensions is similar to the two dimensional result of eq. (43): ( ) dr ds 2 = a
More informationHW#6 is due; please pass it in. HW#7 is posted; it requires you to write out answers clearly and completely explaining your logic.
Oct 21, 2015 Basic properties of groups and clusters of galaxies Number density, structure X-ray emission from hot ICM Estimating cluster masses Cluster scaling relations Impact of environment on galaxies
More informationCollisionless Boltzmann Eq (Vlasov eq)" S+G sec 3.4! Collisionless Boltzmann Eq S&G 3.4!
Collisionless Boltzmann Eq (Vlasov eq)" S+G sec 3.4 When considering the structure of galaxies cannot follow each individual star (10 11 of them), Consider instead stellar density and velocity distributions.
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Numerics Michael Bader SCCS Technical University of Munich Summer 018 Recall: Molecular Dynamics System of ODEs resulting force acting on a molecule: F i = j
More informationA Guide to the Next Few Lectures!
Dynamics and how to use the orbits of stars to do interesting things chapter 3 of S+G- parts of Ch 11 of MWB (Mo, van den Bosch, White) READ S&G Ch 3 sec 3.1, 3.2, 3.4 we are skipping over epicycles 1
More informationAn Introduction to Galaxies and Cosmology. Jun 29, 2005 Chap.2.1~2.3
An Introduction to Galaxies and Cosmology Jun 29, 2005 Chap.2.1~2.3 2.1 Introduction external galaxies normal galaxies - majority active galaxies - 2% high luminosity (non-stellar origin) variability
More informationASTR 610 Theory of Galaxy Formation Lecture 18: Disk Galaxies
ASTR 610 Theory of Galaxy Formation Lecture 18: Disk Galaxies Frank van den Bosch Yale University, spring 2017 The Structure & Formation of Disk Galaxies In this lecture we discuss the structure and formation
More informationIf we see a blueshift on one side and a redshift on the other, this is a sign of rotation.
Galaxies : dynamics, masses, and formation Prof Andy Lawrence Astronomy 1G 2011-12 Overview Spiral galaxies rotate; this allows us to measure masses But there is also a problem : spiral arm winding Elliptical
More informationStars: basic observations
Stars: basic observations Basic properties of stars we would like to know in order to compare theory against observations: Stellar mass M Stellar radius R Surface temperature - effective temperature T
More informationDynamical friction, galaxy merging, and radial-orbit instability in MOND
Dynamical friction, galaxy merging, and radial-orbit instability in MOND Carlo Nipoti Dipartimento di Astronomia Università di Bologna Strasbourg, 29 June 2010 Outline MOND and Equivalent Newtonian Systems
More informationA = 6561 times greater. B. 81 times greater. C. equally strong. D. 1/81 as great. E. (1/81) 2 = 1/6561 as great Pearson Education, Inc.
Q13.1 The mass of the Moon is 1/81 of the mass of the Earth. Compared to the gravitational force that the Earth exerts on the Moon, the gravitational force that the Moon exerts on the Earth is A. 81 2
More informationMichela Mapelli. LECTURES on COLLISIONAL DYNAMICS: 1. RELEVANT TIMESCALES, FORMATION OF STAR CLUSTERS, EQUILIBRIUM MODELS
Michela Mapelli LECTURES on COLLISIONAL DYNAMICS: 1. RELEVANT TIMESCALES, FORMATION OF STAR CLUSTERS, EQUILIBRIUM MODELS COLLISIONAL/COLLISIONLESS? Collisional systems are systems where interactions between
More informationThe dynamical evolution of massive black hole binaries - II. Self-consistent N-body integrations.
The dynamical evolution of massive black hole binaries - II. Self-consistent N-body integrations. Gerald D. Quinlan 1 and Lars Hernquist 2 1 Dept. of Physics and Astronomy, Rutgers University, PO Box 849,
More informationOBSERVATIONAL EVIDENCE FOR DARK MATTER AND DARK ENERGY. Marco Roncadelli INFN Pavia (Italy)
OBSERVATIONAL EVIDENCE FOR DARK MATTER AND DARK ENERGY Marco Roncadelli INFN Pavia (Italy) ABSTRACT Assuming KNOWN physical laws, I first discuss OBSERVATIONAL evidence for dark matter in galaxies and
More informationToward Binary Black Hole Simulations in Numerical Relativity
Toward Binary Black Hole Simulations in Numerical Relativity Frans Pretorius California Institute of Technology BIRS Workshop on Numerical Relativity Banff, April 19 2005 Outline generalized harmonic coordinates
More information18.3 Black Holes: Gravity's Ultimate Victory
18.3 Black Holes: Gravity's Ultimate Victory Our goals for learning: What is a black hole? What would it be like to visit a black hole? Do black holes really exist? What is a black hole? Gravity, Newton,
More informationKozai-Lidov oscillations
Kozai-Lidov oscillations Kozai (1962 - asteroids); Lidov (1962 - artificial satellites) arise most simply in restricted three-body problem (two massive bodies on a Kepler orbit + a test particle) e.g.,
More informationChapter 16 Dark Matter, Dark Energy, & The Fate of the Universe
16.1 Unseen Influences Chapter 16 Dark Matter, Dark Energy, & The Fate of the Universe Dark Matter: An undetected form of mass that emits little or no light but whose existence we infer from its gravitational
More informationConservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body
Continuum Mechanics We ll stick with the Lagrangian viewpoint for now Let s look at a deformable object World space: points x in the object as we see it Object space (or rest pose): points p in some reference
More informationSummary So Far! M87van der Maerl! NGC4342! van den Bosch! rotation velocity!
Summary So Far Fundamental plane connects luminosity, scale length, surface brightness, stellar dynamics. age and chemical composition Elliptical galaxies are not randomly distributed within the 3D space
More informationPhysics 351, Spring 2017, Homework #4. Due at start of class, Friday, February 10, 2017
Physics 351, Spring 2017, Homework #4. Due at start of class, Friday, February 10, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page
More informationThe physical origin of stellar envelopes around globular clusters
The physical origin of stellar envelopes around globular clusters Phil Breen University of Edinburgh in collaboration with A. L. Varri, J. Peñarrubia and D. C. Heggie Current observational evidence Example:
More informationDynamical Models of the Globular Clusters M4 and NGC 6397
Dynamical Models of the Globular Clusters M4 and NGC 6397 Douglas Heggie University of Edinburgh d.c.heggie@ed.ac.uk in collaboration with Mirek Giersz CAMK, Warsaw, Poland Outline Introduction: M4 and
More informationLecture 14: Cosmological Principles
Lecture 14: Cosmological Principles The basic Cosmological Principles The geometry of the Universe The scale factor R and curvature constant k Comoving coordinates Einstein s initial solutions 3/28/11
More informationPhys 100 Astronomy (Dr. Ilias Fernini) Review Questions for Chapter 9
Phys 0 Astronomy (Dr. Ilias Fernini) Review Questions for Chapter 9 MULTIPLE CHOICE 1. We know that giant stars are larger in diameter than the sun because * a. they are more luminous but have about the
More informationInstability and different burning regimes
1 X-ray bursts Last time we talked about one of the major differences between NS and BH: NS have strong magnetic fields. That means that hot spots can be produced near the magnetic poles, leading to pulsations
More informationCOLOR MAGNITUDE DIAGRAMS
COLOR MAGNITUDE DIAGRAMS What will you learn in this Lab? This lab will introduce you to Color-Magnitude, or Hertzsprung-Russell, Diagrams: one of the most useful diagnostic tools developed in 20 th century
More informationGalaxies: Structure, Dynamics, and Evolution. Elliptical Galaxies (II)
Galaxies: Structure, Dynamics, and Evolution Elliptical Galaxies (II) Layout of the Course Feb 5: Review: Galaxies and Cosmology Feb 12: Review: Disk Galaxies and Galaxy Formation Basics Feb 19: Disk Galaxies
More informationLecture XIX: Particle motion exterior to a spherical star
Lecture XIX: Particle motion exterior to a spherical star Christopher M. Hirata Caltech M/C 350-7, Pasadena CA 95, USA Dated: January 8, 0 I. OVERVIEW Our next objective is to consider the motion of test
More informationN-body Simulations and Dark energy
N-Body Simulations and models of Dark Energy Elise Jennings Supported by a Marie Curie Early Stage Training Fellowship N-body Simulations and Dark energy elise jennings Introduction N-Body simulations
More informationGalaxies and Cosmology
4/28/17 The Discovery of Galaxies Up to the 1920 s, astronomers were not sure exactly how far away galaxies were, and thus didn t know how big they are! Spiral Nebulae could be assumed to be inside our
More informationModified Gravity (MOG) and Dark Matter: Can Dark Matter be Detected in the Present Universe?
Modified Gravity (MOG) and Dark Matter: Can Dark Matter be Detected in the Present Universe? John Moffat Perimeter Institute, Waterloo, Ontario, Canada Talk given at the Miami 2014 topical conference on
More informationHierarchical Bayesian Modeling
Hierarchical Bayesian Modeling Making scientific inferences about a population based on many individuals Angie Wolfgang NSF Postdoctoral Fellow, Penn State Astronomical Populations Once we discover an
More informationName Final Exam December 7, 2015
Name Final Exam December 7, 015 This test consists of five parts. Please note that in parts II through V, you can skip one question of those offered. Part I: Multiple Choice (mixed new and review questions)
More informationThe cosmic distance scale
The cosmic distance scale Distance information is often crucial to understand the physics of astrophysical objects. This requires knowing the basic properties of such an object, like its size, its environment,
More informationDark Matter and the DRIFT. Experiment. Why do we believe. How shall we identify it? What is Dark. Matter? in it? C. J. Martoff, Professor of Physics
Dark Matter and the DRIFT Experiment What is Dark Matter? Why do we believe in it? How shall we identify it? C J Martoff, Professor of Physics Physics C056, Astronomy Temple University Spring '03 What
More informationDynamics of Galaxies: Practice. Frontiers in Numerical Gravitational Astrophysics July 3, 2008
Dynamics of Galaxies: Practice Frontiers in Numerical Gravitational Astrophysics July 3, 2008 http://xkcd.com/323/ Outline 1. Initial Conditions a) Jeans Theorem b) Exact solutions for spheres c) Approximations
More informationAnnouncements. Homework. Set 8now open. due late at night Friday, Dec 10 (3AM Saturday Nov. 11) Set 7 answers on course web site.
Homework. Set 8now. due late at night Friday, Dec 10 (3AM Saturday Nov. 11) Set 7 answers on course web site. Review for Final. In class on Thursday. Course Evaluation. https://rateyourclass.msu.edu /
More informationDARK MATTER IN UNIVERSE. edited by. John Bahcall Institute for Advanced Study, Princeton, USA. Tsvi Piran The Hebrew University, Israel
DARK MATTER IN UNIVERSE n d edited by John Bahcall Institute for Advanced Study, Princeton, USA Tsvi Piran The Hebrew University, Israel Steven Weinberg University of Texas, Austin, USA TECHNiSCHE INFORMATIONSBIBLIOTHEK
More informationA5682: Introduction to Cosmology Course Notes. 2. General Relativity
2. General Relativity Reading: Chapter 3 (sections 3.1 and 3.2) Special Relativity Postulates of theory: 1. There is no state of absolute rest. 2. The speed of light in vacuum is constant, independent
More information