Post-Newtonian N-body Codes. Sverre Aarseth. Institute of Astronomy, Cambridge
|
|
- Matthew Paul
- 6 years ago
- Views:
Transcription
1 Post-Newtonian N-body Codes Sverre Aarseth Institute of Astronomy, Cambridge Introduction N-body tools Three-body formulation PN implementations Numerical examples Discussion
2 Hermite Integration Taylor series for F and F (1) F = F 0 + F (1) 0 t F(2) 0 t F(3) 0 t3 F (1) = F (1) 0 + F (2) 0 t F(3) 0 t2 Prediction r j = (( 1 6 F(1) 0 δt j F 0)δt j + v 0)δt j +r 0 v j = (( 1 2 F(1) 0 δt j + F 0)δt j + v 0; δt j = t t 0 New forces F, F (1) Higher derivatives F (3) 0 = (2(F 0 F) + (F (1) 0 +F (1) )t) 6 t 3 F (2) 0 = ( 3(F 0 F) (2F (1) 0 +F (1) )t) 2 t 2 Corrector for i r i = 1 24 F(2) 0 t F(3) 0 t5 v i = 1 6 F(2) 0 t F(3) 0 t4
3 AC Neighbour Scheme Total force F(t) = n j=1 F j +F d (t) Prediction F(t) = F n +Ḟ d (t t 0 )+F d (t 0 ) Ḟ = Ḟ n +Ḟ d Time-scales t n << t d, n << N Neighbour sphere R new s = R old s ( ) np 1/3 n Neighbour selection r i r j < R s Derivative corrections F (2) ij, F(3) ij
4 Basic KS Relations New coordinates R = u 2 1+u 2 2+u 2 3+u 2 4 Time transformation Coordinate transformation Generalized Levi Civita matrix Inverse relations L(u) = dt = R dτ R = L(u) u u 1 u 2 u 3 u 4 u 2 u 1 u 4 u 3 u 3 u 4 u 1 u 2 u 4 u 3 u 2 u 1 R 1 = u 2 1 u 2 2 u 2 3+u 2 4 R 2 = 2(u 1 u 2 u 3 u 4 ) R 3 = 2(u 1 u 3 +u 2 u 4 ) R 4 = 0 Velocity Ṙ = 2Lu R, R = t = 2u u Equations of motion u = 1 2 hu+ 1 2 RLT P h = 2u L T P t = u u
5 KS Decision-Making Close encounter Time-step criterion Neighbour list search R cl = 4r h N C 1/3, t cl = β t k < t cl list all r 2 kj, t j < 2 t cl cl R 3 1/2 m Two-body selection R kl < R cl, Ṙ kl < 0 Dominant motion m k + m l R 2 kl > m k + m j R 2 kj KS initialization Initialization of c.m. F U, F U, τ & t (n) t r cm = m kr k + m l r l m k + m l Perturber search r p < 2m p m b γ min 1/3 a(1+e) Slow-down adjustment Termination test Delayed termination Final iteration Polynomial initialization γ < γ 0, τ κ t R > R 0, γ > γ T block t > t i τ from τ, τ,... and δt F j, Ḟ j, t j, j = k, l
6 Hermite KS with block-steps Predict u, u, h to highest order u pred = u+ τu + τ2 2 u + τ3 6 u + τ4 24 u 4 + τ 120 u pred = u + τu + τ2 2 u + τ3 6 u 4 + τ4 24 u 5, h pred = h+ τh + τ 2 h + τ 6 h + τ 24 h. Force evaluation P = L T uf, P = L T u F +rl T u F. (2) Derivatives for corrector u = 1 2 (hu+rp), u = 1 2 (h u+hu +r P +rp ), h = 2(u P), h = 2(u P +u P ). (3)
7 Fourth-order corrector (Hut, Makino & McMillan 1995) u corr = u old + τ 2 (u corr +u τ2 old ) 12 (u new u old ), u corr = u old + τ 2 (u new +u τ2 old ) 12 (u new u old ), h corr = h old + τ 2 (h new +h τ2 old ) 12 (h new h old ). (4) Velocity derivatives for prediction u 4 mid = u new u τ u 5 mid = 12 τ 2 old, ( u new +u old 2 u new u ) old. (5) τ Shifted from midpoint u 4 new = u 4 mid + τ 2 u 5 mid, u 5 new = u 5 mid, h new = h mid + τ h new = h 2 h mid, mid. (6)
8 Taylor series for physical time t = u u, t = 2(u u ), t = 2(u u +u u ), t 4 = 2(u u )+6(u u ), t 5 = 2(u u 4 +8(u u )+6(u u ), t 6 = 2(u u 5 +10(u u 4 +20(u u ). (7) Next look-up time t next = t+ 6 n=1 t n τ n. (8) n! Regularized time-step by inversion t+ 6 n=1 t n τ n n! = 0. (9) Newton-Raphson iteration to solve for τ x i+1 = x i f(x i) f (x i ).
9 Initial guess τ 1 = t bs t τ orig, with quantized time-step t bs (= 2 n t orig ). Commensurability condition mod (t sys, t bs ) = 0. Solve equation (9) to obtain τ bs.
10 Three-Body Regularization Initialconditions r i, p i, p 3 = (p 1 + p 2 ) Basic Hamiltonian µ k3 = m km 3 m k +m 3 H = 2 k=1 1 p 2 2µ k+ 1 p 1 p 2 m 1m 3 m 2m 3 m 1m 2 k3 m 3 R 1 R 2 R KS coordinate transformation Q 2 k = R k, (k = 1,2) Time transformation dt = R 1 R 2 dτ Regularized Hamiltonian Γ = R 1 R 2 (H E 0 ) Γ = 2 k=1 1 R 3 k P 2 k + 1 P T 8µ k3 16m 1A 1 A T 2P 2 3 m 1 m 3 R 2 m 2 m 3 R 1 m 1m 2 R 1 R 2 R 1 R 2 E 0R 1 R 2 Equations of motion dq k dτ = Γ P k, dp k dτ = Γ Q k Regular solutions: R 1 0 or R 2 0 Singularterm<regularterms: R 1 R 2 > max(r 1, R 2 )
11 Three-Body Transformations Coordinates & momenta q k = q k q 3, p k = p k Regularized coordinates (q 1 0) Q 1 = [ 1 2 ( q 1 +q 1 )] 1/2 Q 2 = 1 2 q 2/Q 1 Q 3 = 1 2 q 3/Q 1 Q 4 = 0 Regularized momenta Basic matrix A 1 = 2 P k = A k p k Q 1 Q 2 Q 3 Q 4 Q 2 Q 1 Q 4 Q 3 Q 3 Q 4 Q 1 Q 2 Q 4 Q 3 Q 2 Q 1 KS transformations Physical momenta q k = 1 2 A k T Q k p k = 1 4 A k T P k /R k Coordinates & momenta q 3 = 2 m k q k /M k=1 q k = q 3 +q k p k = p k p 3 = (p 1 +p 2 ) (k = 1,2)
12 PN Codes Two-body Hermite KS toy code alternative PN formulations Three-body AZ regularization dominant binary, sequential terms homework: Einstein shift Chain Algorithmic regularization compact sub-system sequential PN, unperturbed KS
13 6x10-8 Two-body KS 5x10-8 4x10-8 SEMI 3x10-8 2x10-8 1x x10-9 1x x10-8 2x x10-8 3x x10-8 4x x10-8 5x10-8 Time
14 V/C = V/C = ECC Time
15 0.006 V/C = V/C = SEMI Time
16 0.006 IPN = 1 IPN = SEMI Time
17
18 Initialization Two-body R < R cl, Ṙ < 0, KS Hierarchy B + S H, Stability test Triple chain a 1 (1 e 1 ) < 3a 0, ṙ < 0, T Four-body chain B +B, ṙ < 0, B-B Higher orders T +B, ṙ < 0, Q
19 Chain Initialization Unperturbed exit ω > Activate chain N ch = 2 or > 2 Alternative path perturbed KS search Selection chain members & c.m. Initialization c.m. force & chain vectors Synchronization block time-step Interface perturber list
20 Termination Hyperbolic KS R > R 0 or γ > Binary KS γ > 0.1, γ = (F 1 F 2 )R 2 m 1 +m 2 Collision KS S Triple d 3 > 20a, B +S Hierarchy a 1 (1 e 1 ) < Sa 0, γ > 0.01 ARchain R 1 << R 2 +R 3 B +2S
21 Relativistic Effects Einstein shift ω = 6πM ac 2 (1 e 2 ) GR radiation time-scale t GR = 5 c 5 g(e)a 4 64X(1+X)m 3 1 X = m 2, g(e) (1 e2 ) 7/2 m Speed of light & Schwarzschild radius c = V, R Sch = 2M c 2 Velocities ( M V = α R pc ) 1/2, v max = ( M a ) 1/2 ( ) 1/2 1+e 1 e PN conditions ω > , t GR < t, r < n R Sch
22 Post-Newtonian Terms Equation of motion d 2 r dt = M [ ( 1+A) r ] 2 r 2 r +Bv First-order precession M = m 1 +m 2, η = m 1m 2 M 2 A 1 = 2(2+η) M r (1+3η)v ηṙ2 B 1 = 2(2 η)ṙ Higher-order precession A 2 =..., B 2 =..., A 3 =..., B 3 =... Gravitational radiation A 5/2 = 8 ( ) 17M 5 ηm r ṙ 3r +3v2 (3 Mr ) +v2 B 5/2 = 8 5 ηm r Total GR perturbation P GR = M c 2 r 2 [ (A 1 + A 2 c + A 5/2 2 c )r 3 r + (B 1 + B 2 c + B ] 5/2 2 c )v 3 GR spin dominant BH ARC energy change E = m 1m 2 M (P GR +P dyn ) vdt
23 Decision-Making Increasing GR effect t GR < 500, 50, 1 IPN = 1, 2, 3 GR time-scale T z = 150a4 M 3 (1 e2 ) 7/2 Myr Graduated PN w > (1,10,100) 10 4 IPN = 1, 2, 3 Coalescence R < 4R Sch = 8M c 2 Alternative merging IPN = 3, N = 2, N p = 0 IPN 2, a(1 e) < R Sch IPN = 3, N = 3, a 1 (1 e 1 ) > 100a Unperturbed KS Tidal disruption t GR < 500 & w > , γ < KS or CHAIN
24 PN Elements Energy ǫ b = ǫ 0 + ǫ 1 c 2 + ǫ 2 c 4 + ǫ 3 c 6, a = M 2ǫ b ǫ 0 = 1 2 V 2 M R, η = m 1m 2 M 2 ǫ 1 = 1 2 M R +3 8 (1 3η)V ( (3+η)V 2 +ηṙ 2) M R Angular momentum J = J 0 (1+f 1 /c 2 +f 2 /c 4 ) Eccentricity e 2 = (1 J2 Ma ) Lenz vector e = V R V/M R/R Periapse advance ω = 6πM c 2 a(1 e 2 ) Kozai-Lidov cycles T Kozai = T2 out T in ( ) 1+qout q out (1 e 2 out )3/2 g(e in,ω in,ψ)
25 Unperturbed GR Orbit Compact objects max (K 1,K 2) > 12 Derivatives ȧ = 64M 1M 2 (M 1 +M 2 ) 5c 5 a 3 (1 e 2 ) 7/2 ė = 304eM 1M 2 (M 1 +M 2 ) g(e) 15c 5 a 4 (1 e 2 ) 5/2 Einstein shift ω = 6π(M 1+M 2 ) ac 2 (1 e 2 ) KS rotation New elements θ = ω t T K, rotate by θ/2 a 1 = a+ȧ t, e 1 = max (e+ė t,0) Energy updating E = µ(h old H new ) Modify KS orbit Shrink, e = const; decrease, H = const Validity Unperturbed KS < ω < γ <
26 Membership Reduction escape R 1 or R n 1 > R, ṙ > 0 Double KS n = 4, two well separated pairs Injection γ > & r p < R, ṙ p < 0 Retention n = 2, IPN > 0, ω > Perturbers massive binary, N p 3 10 Perturbation γ 2m p M bh max ( ) 3 Rk d p
27 PN1 and PN Energy Time
28 e-06 Unperturbed PN 7.144e e e e e-06 A e e e e e e Time
29 Energy Corrections Unperturbed KS < ω < Perturbed KS P 2.5 & Ṗ 2.5 added to F & Ḟ Energy derivatives Ė = P v, Ë = P v+ṗ v Hermite method h = 1 2 (Ė 0 +Ė) t (Ë 0 Ë) t 2 Energy loss µ h Conservation E ch = E sub E E tot = E N +E ch
30 1 ECC
31 Tidal Disruption ( mbh ) 1/3 Classical condition r t = r m Case A Two-body regularization Corrections Case B Mass transfer Ejection Updating Initialization Pericentre test: a(1 e) < r t Transform to apocentre θ = π/2 Mass transfer m = 0.1m Escape velocity kick KS elements and energy Chain regularization m = 0.1m Energy corrections Chain membership and tables Chain c.m. polynomials Requirements Effects Special cases Astrophysics Observations Large eccentricity Mass segregation & Kozai-Lidov WD & AGB cores Instantaneous accretion Light flickering & emission lines
32 0.378 Model ECC Time
33 3.625e-08 Model e e e-08 SEMI 3.605e e e e Time
34 0.001 Model BH e-05 A 1e-06 1e-07 1e Time
35 Model BH8 1e-05 1e-06 a 1e-07 1e-08 1e-09 1e Time
36 1000 Model BH TGR e-05 1e-06 1e Time
37 Discussion IMF Mass outcome upper limit, Z BH from BSE PN simulations speed of light c = V Dynamics binary shrinkage & Kozai-Lidov Slingshot Retention Depletion Observations Astropysics fast escapers velocity kicks ejections of BH & BBH BH + S binaries, hyper-velocities tidal disruption & BH accretion
Parallel KS Block-Step Method. Sverre Aarseth. Institute of Astronomy, Cambridge
Parallel KS Block-Step Method Sverre Aarseth Institute of Astronomy, Cambridge Code Overview Hermite KS Prediction & Correction Iteration Time-Steps Decision-Making Binary Project Code Overview Directories
More informationAveraging the average: Morphology transitions in spin precession of black-hole binaries
Averaging the average: Morphology transitions in spin precession of black-hole binaries U. Sperhake DAMTP, University of Cambridge M. Kesden, D. Gerosa, R. O Shaughnessy, E. Berti VII Black Holes Workshop
More informationDynamics of Stars and Black Holes in Dense Stellar Systems:
Michela Mapelli INAF - Padova Dynamics of Stars and Black Holes in Dense Stellar Systems: Lecture VI: DYNAMICS AROUND SUPER-MASSIVE BHs 0. nuclear star clusters (NSCs) 1. dynamics around super-massive
More informationMassive star clusters
Massive star clusters as a host of compact binaries Michiko Fujii ( 藤井通子 ) The University of Tokyo Outline Massive star clusters and compact binaries Dynamical evolution of star clusters Distribution of
More informationStellar Dynamics and Structure of Galaxies
Stellar Dynamics and Structure of Galaxies in a given potential Vasily Belokurov vasily@ast.cam.ac.uk Institute of Astronomy Lent Term 2016 1 / 59 1 Collisions Model requirements 2 in spherical 3 4 Orbital
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 informationPost-Keplerian effects in binary systems
Post-Keplerian effects in binary systems Laboratoire Univers et Théories Observatoire de Paris / CNRS The problem of binary pulsar timing (Credit: N. Wex) Some classical tests of General Relativity Gravitational
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 informationStellar mass black holes in young massive and open clusters and their role in gravitational-wave generation
Stellar mass black holes in young massive and open clusters and their role in gravitational-wave generation Sambaran Banerjee Argelander-Institut für Astronomie (AIfA) and Helmholtz-Instituts für Strahlen-
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 informationBlue Straggler Stars Formation Channels
Blue Straggler Stars Formation Channels Melvyn B. Davies Department of Astronomy and Theoretical Physics Lund University www.astro.lu.se Ways to make blue stragglers Produced in collisions/mergers between
More informationStellar Black Hole Binary Mergers in Open ClusterS
Stellar Black Hole Binary Mergers in Open ClusterS SARA RASTELLO Collaborators: P. Amaro-Seoane, M. Arca-Sedda, R. Capuzzo Dolcetta, G. Fragione, I.Tosta e Melo Modest 18, Santorini, greece 27/06/2018
More information2.5.1 Static tides Tidal dissipation Dynamical tides Bibliographical notes Exercises 118
ii Contents Preface xiii 1 Foundations of Newtonian gravity 1 1.1 Newtonian gravity 2 1.2 Equations of Newtonian gravity 3 1.3 Newtonian field equation 7 1.4 Equations of hydrodynamics 9 1.4.1 Motion of
More informationOpen problems in compact object dynamics
Open problems in compact object dynamics Melvyn B. Davies Department of Astronomy and Theoretical Physics Lund University www.astro.lu.se Key ideas and open questions in compact object dynamics Melvyn
More informationImplementing Few-Body Algorithmic Regularization with Post-Newtonian terms
Rochester Institute of Technology RIT Scholar Works Articles 5-14-2008 Implementing Few-Body Algorithmic Regularization with Post-Newtonian terms Seppo Mikkola University of Turku David Merritt Rochester
More informationPost-Newtonian N-body Dynamics
Minor Master Project Post-Newtonian N-body Dynamics Author: Emiel Por por@strw.leidenuniv.nl Supervisors: Simon F. Portegies Zwart Adrian S. Hamers October 2014 1 Post-Newtonian N-body Dynamics Emiel Por
More informationForming Intermediate-Mass Black Holes in Dense Clusters Through Collisional Run-away
Forming Intermediate-Mass Black Holes in Dense Clusters Through Collisional Run-away Marc Freitag ARI, Heidelberg freitag@ari.uni-heidelberg.de http://obswww.unige.ch/~freitag Collaboration with Atakan
More informationBinary sources of gravitational waves
Binary sources of gravitational waves For our final two lectures we will explore binary systems in general and the Advanced LIGO detections in particular. Binaries obviously have a large and varying quadrupole
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 informationEvolution of binary supermassive black holes via chain regularization
Rochester Institute of Technology RIT Scholar Works Presentations and other scholarship 2005 Evolution of binary supermassive black holes via chain regularization Andras Szell David Merritt Seppo Mikkola
More informationarxiv:astro-ph/ v1 15 Feb 2002
DRAFT VERSION OCTOBER 6, 018 Preprint typeset using L A TEX style emulateapj FOUR-BODY EFFECTS IN GLOBULAR CLUSTER BLACK HOLE COALESCENCE M. COLEMAN MILLER AND DOUGLAS P. HAMILTON Department of Astronomy,
More informationDynamic Exoplanets. Alexander James Mustill
Dynamic Exoplanets Alexander James Mustill Exoplanets: not (all) like the Solar System Exoplanets: not (all) like the Solar System Solar System Lissauer et al 14 Key questions to bear in mind What is role
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 informationThe Role of Binary-Binary Interactions in Inducing Eccentric Black Hole Mergers
The Role of Binary-Binary Interactions in Inducing Eccentric Black Hole Mergers - CIERA / Northwestern University with Johan Samsing (Princeton), Carl Rodriguez (MIT), Carl-Johan Haster (CITA), & Enrico
More informationStellar collisions and their products
Stellar collisions and their products Melvyn B. Davies Department of Astronomy and Theoretical Physics Lund University www.astro.lu.se KEY IDEA #1 Collision rate depends on V. Stellar encounter timescales
More informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2: The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
More informationarxiv: v1 [math.ds] 6 Dec 2015
THE BINARY RETURNS! CONNOR JACKMAN arxiv:1512.01852v1 [math.ds] 6 Dec 2015 Abstract. Consider the spatial Newtonian three body problem at fixed negative energy and fixed angular momentum. The moment of
More informationAdvisor : Prof. Hagai Perets. Israel Institute of Technology, Haifa
Secular Dynamics in Hierarchical Three-Body Systems with Mass Loss and MassErez Michaely Advisor : Prof. Hagai Perets Israel Institute of Technology, Haifa Context Systematically exploring triple systems
More informationThe effect of f - modes on the gravitational waves during a binary inspiral
The effect of f - modes on the gravitational waves during a binary inspiral Tanja Hinderer (AEI Potsdam) PRL 116, 181101 (2016), arxiv:1602.00599 and arxiv:1608.01907? A. Taracchini F. Foucart K. Hotokezaka
More informationA Numerical Study of Boson Star Binaries
A Numerical Study of Boson Star Binaries Bruno C. Mundim with Matthew W. Choptuik (UBC) 12th Eastern Gravity Meeting Center for Computational Relativity and Gravitation Rochester Institute of Technology
More informationComputational Problem: Keplerian Orbits
Computational Problem: Keplerian Orbits April 10, 2006 1 Part 1 1.1 Problem For the case of an infinite central mass and an orbiting test mass, integrate a circular orbit and an eccentric orbit. Carry
More informationFORMATION AND EVOLUTION OF COMPACT BINARY SYSTEMS
FORMATION AND EVOLUTION OF COMPACT BINARY SYSTEMS Main Categories of Compact Systems Formation of Compact Objects Mass and Angular Momentum Loss Evolutionary Links to Classes of Binary Systems Future Work
More informationHow black holes get their kicks! Gravitational radiation recoil from binary inspiral and plunge into a rapidly-rotating black hole.
How black holes get their kicks! Gravitational radiation recoil from binary inspiral and plunge into a rapidly-rotating black hole. Marc Favata (Cornell) Daniel Holz (U. Chicago) Scott Hughes (MIT) The
More informationSPIN PRECESSION IN A 2 BODY SYSTEM: A NEW TEST OF GENERAL RELATIVITY R. F. O CONNELL DEPT. OF PHYSICS & ASTRONOMY LOUISIANA STATE UNIVERSITY
SPIN PRECESSION IN A 2 BODY SYSTEM: A NEW TEST OF GENERAL RELATIVITY R. F. O CONNELL DEPT. OF PHYSICS & ASTRONOMY LOUISIANA STATE UNIVERSITY 1 1. Newtonian Theory (p. 2) 2. General Relativistic Corrections
More informationTesting GR with Compact Object Binary Mergers
Testing GR with Compact Object Binary Mergers Frans Pretorius Princeton University The Seventh Harvard-Smithsonian Conference on Theoretical Astrophysics : Testing GR with Astrophysical Systems May 16,
More informationGravitational waves from binary black holes
Gravitational waves from binary black holes Hiroyuki Nakano YITP, Kyoto University DECIGO workshop, October 27, 2013 Hiroyuki Nakano Gravitational waves from binary black holes Binary black holes (BBHs)
More information6 th lecture of Compact Object and Accretion, Master Programme at Leiden Observatory
6 th lecture of Compact Object and Accretion, Master Programme at Leiden Observatory Accretion 1st class study material: Chapter 1 & 4, accretion power in astrophysics these slides at http://home.strw.leidenuniv.nl/~emr/coa/
More informationGravitational waves from compact objects inspiralling into massive black holes
Gravitational waves from compact objects inspiralling into massive black holes Éanna Flanagan, Cornell University American Physical Society Meeting Tampa, Florida, 16 April 2005 Outline Extreme mass-ratio
More informationClassical and Quantum Dynamics in a Black Hole Background. Chris Doran
Classical and Quantum Dynamics in a Black Hole Background Chris Doran Thanks etc. Work in collaboration with Anthony Lasenby Steve Gull Jonathan Pritchard Alejandro Caceres Anthony Challinor Ian Hinder
More informationChapter 11. Angular Momentum
Chapter 11 Angular Momentum Angular Momentum Angular momentum plays a key role in rotational dynamics. There is a principle of conservation of angular momentum. In analogy to the principle of conservation
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 informationCelestial Mechanics of Asteroid Systems
Celestial Mechanics of Asteroid Systems D.J. Scheeres Department of Aerospace Engineering Sciences The University of Colorado scheeres@colorado.edu 1 Granular Mechanics and Asteroids Asteroid systems are
More informationPhys 7221, Fall 2006: Homework # 5
Phys 722, Fal006: Homewor # 5 Gabriela González October, 2006 Prob 3-: Collapse of an orbital system Consier two particles falling into each other ue to gravitational forces, starting from rest at a istance
More informationCoalescing Binary Black Holes Originating from Globular Clusters
Coalescing Binary Black Holes Originating from Globular Clusters Abbas Askar (askar@camk.edu.pl) Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, Warsaw, Poland Abbas Askar Globular
More informationEvolution of SMBH-SMBH and SMBH-IMBH Binaries: Effect of Large Mass Ratio
Evolution of SMBH-SMBH and SMBH-IMBH Binaries: Effect of Large Mass Ratio Jun Makino Center for Computational Astrophysics and Division Theoretical Astronomy National Astronomical Observatory of Japan
More informationALGORITHMIC REGULARIZATION CHAIN (state of the art 2010) Seppo Mikkola
ALGORITHMIC REGULARIZATION CHAIN (state of the art 2010) Seppo Mikkola University of Turku, Finland Department of Physics and Astronomy Tuorla Observatory e-mail: Seppo.Mikkola@utu.fi H = T U = k p 2 k
More informationA map approximation for the restricted three-body problem
A map approximation for the restricted three-body problem Shane Ross Engineering Science and Mechanics, Virginia Tech www.esm.vt.edu/ sdross Collaborators: P. Grover (Virginia Tech) & D. J. Scheeres (U
More informationThe Gravitational Radiation Rocket Effect. Marc Favata Cornell University GR17, Dublin, July 2004
The Gravitational Radiation Rocket Effect recoil Marc Favata Cornell University GR17, Dublin, July 004 Favata, Hughes, & Holz, ApJL 607, L5, astro-ph/040056 Merritt, Milosavljevic, Favata, Hughes, & Holz,
More informationHOMEWORK 10. Applications: special relativity, Newtonian limit, gravitational waves, gravitational lensing, cosmology, 1 black holes
General Relativity 8.96 (Petters, spring 003) HOMEWORK 10. Applications: special relativity, Newtonian limit, gravitational waves, gravitational lensing, cosmology, 1 black holes 1. Special Relativity
More informationPost-Newtonian Approximation
Post-Newtonian Approximation Piotr Jaranowski Faculty of Physcis, University of Bia lystok, Poland 01.07.2013 1 Post-Newtonian gravity and gravitational-wave astronomy 2 3 4 EOB-improved 3PN-accurate Hamiltonian
More information3.4 Numerical orbit integration
196 Chapter 3: The Orbits of Stars about the center of the galaxy. For small-amplitude librations, the libration frequency is p, consistent with our assumption that the oscillation frequency is of order
More informationarxiv:astro-ph/ v1 27 Mar 2001
ApJ, in press On the origin of the distribution of binary-star periods Pavel Kroupa 1 and Andreas Burkert 2 arxiv:astro-ph/0103429v1 27 Mar 2001 1 Institut für Theoretische Physik und Astrophysik Universität
More informationStellar-mass black holes in a globular cluster
Stellar-mass black holes in a globular cluster Douglas Heggie University of Edinburgh, UK d.c.heggie@ed.ac.uk Motivation Discovery of two stellar-mass BH in the Galactic globular cluster M22 (Strader+
More informationCurved spacetime tells matter how to move
Curved spacetime tells matter how to move Continuous matter, stress energy tensor Perfect fluid: T 1st law of Thermodynamics Relativistic Euler equation Compare with Newton =( c 2 + + p)u u /c 2 + pg j
More informationSpin and quadrupole moment effects in the post-newtonian dynamics of compact binaries. László Á. Gergely University of Szeged, Hungary
Spin and quadrupole moment effects in the post-newtonian dynamics of compact binaries László Á. Gergely University of Szeged, Hungary University of Cardiff, UK Spinning Coalescing Binaries Workshop - September
More informationPulsars ASTR2110 Sarazin. Crab Pulsar in X-rays
Pulsars ASTR2110 Sarazin Crab Pulsar in X-rays Test #2 Monday, November 13, 11-11:50 am Ruffner G006 (classroom) Bring pencils, paper, calculator You may not consult the text, your notes, or any other
More informationκ = f (r 0 ) k µ µ k ν = κk ν (5)
1. Horizon regularity and surface gravity Consider a static, spherically symmetric metric of the form where f(r) vanishes at r = r 0 linearly, and g(r 0 ) 0. Show that near r = r 0 the metric is approximately
More informationGravitational Waves from Coalescing Binaries and the post-newtonian Theory
Gravitational Waves from Coalescing Binaries and the post-newtonian Theory Riccardo Sturani Instituto de Física Teórica UNESP/ICTP-SAIFR São Paulo (Brazil) Ubu - Anchieta, April 16 th 2015 Riccardo Sturani
More information2 Ivanova, Fregeau, & Rasio
2 Ivanova, Fregeau, & Rasio Figure 1. Orbital period vs companion mass for all observed binary MSPs in the galactic field (open circles), in 47 Tuc (solid squares) and in other globular clusters (solid
More informationFrom space-time to gravitation waves. Bubu 2008 Oct. 24
From space-time to gravitation waves Bubu 008 Oct. 4 Do you know what the hardest thing in nature is? and that s not diamond. Space-time! Because it s almost impossible for you to change its structure.
More informationarxiv:astro-ph/ v2 15 Jan 2007
Mon. Not. R. Astron. Soc. 000, 000 000 (0000) Printed 16 January 2007 (MN LATEX style file v2.2) A hypervelocity star from the Large Magellanic Cloud Alessia Gualandris 1 and Simon Portegies Zwart 1 1
More informationarxiv: v2 [astro-ph.ga] 21 Dec 2017
Mon. Not. R. Astron. Soc. 000, 1 14 (2017) Printed 25 December 2017 (MN LATEX style file v2.2) arxiv:1710.10671v2 [astro-ph.ga] 21 Dec 2017 The fate of close encounters between binary stars and binary
More informationThe two-body Kepler problem
The two-body Kepler problem set center of mass at the origin (X = 0) ignore all multipole moments (spherical bodies or point masses) define r := r 1 r 2,r:= r,m:= m 1 + m 2,µ:= m 1 m 2 /m reduces to effective
More informationSavvas Nesseris. IFT/UAM-CSIC, Madrid, Spain
Savvas Nesseris IFT/UAM-CSIC, Madrid, Spain What are the GWs (history, description) Formalism in GR (linearization, gauges, emission) Detection techniques (interferometry, LIGO) Recent observations (BH-BH,
More informationGravitational Waves. Masaru Shibata U. Tokyo
Gravitational Waves Masaru Shibata U. Tokyo 1. Gravitational wave theory briefly 2. Sources of gravitational waves 2A: High frequency (f > 10 Hz) 2B: Low frequency (f < 10 Hz) (talk 2B only in the case
More informationMonte Carlo Models of Dense Star Clusters
Monte Carlo Models of Dense Star Clusters Sourav Chatterjee MODEST 18, Firá, Santorini June 26, 2018 Hénon-type Monte Carlo N-body model spread all masses in spherical shells (E, J) 2-body relaxation approximate
More informationRelativistic loss-cone dynamics: Implications for the Galactic Center
Relativistic loss-cone dynamics: Implications for the Galactic Center Relativistic loss-cone dynamics: Implications for the Galactic Center Tal Alexander Weizmann Institute of cience 3 RR < NR Plunge log
More informationBlack Hole-Neutron Star Binaries in General Relativity. Thomas Baumgarte Bowdoin College
Black Hole-Neutron Star Binaries in General Relativity Thomas Baumgarte Bowdoin College 1 Why do we care? Compact binaries (containing neutron stars and/or black holes) are promising sources of gravitational
More informationGravitational radiation
Lecture 28: Gravitational radiation Gravitational radiation Reading: Ohanian and Ruffini, Gravitation and Spacetime, 2nd ed., Ch. 5. Gravitational equations in empty space The linearized field equations
More informationGravitational Waves Theory - Sources - Detection
Gravitational Waves Theory - Sources - Detection Kostas Glampedakis Contents Part I: Theory of gravitational waves. Properties. Wave generation/the quadrupole formula. Basic estimates. Part II: Gravitational
More informationGravitational Waves from Boson Stars
Gravitational Waves from Boson Stars Ruxandra Bondarescu (ICG, Portsmouth) Gregory Daues (NCSA) Jayashree Balakrishna (Harris Stowe State U) Edward Seidel (NCSA) CQG 23, 2631 (2006), gr-qc/0602078, Phys.
More informationPinhole Cam Visualisations of Accretion Disks around Kerr BH
Pinhole Camera Visualisations of Accretion Disks around Kerr Black Holes March 22nd, 2016 Contents 1 General relativity Einstein equations and equations of motion 2 Tetrads Defining the pinhole camera
More informationThe first order formalism and the transition to the
Problem 1. Hamiltonian The first order formalism and the transition to the This problem uses the notion of Lagrange multipliers and Legendre transforms to understand the action in the Hamiltonian formalism.
More information3. Dynamics and test particle motion.
. Dynamics and test particle motion... X-ray Binaries ystems with stars orbiting around their common center of mass are relatively abundant. If one of the components is a BH we have a chance to determine
More informationBinary black holes in nuclei of extragalactic radio sources 3C 345. VLBI component = radio blob superluminal motion = apparent motion
Binary black holes in nuclei of extragalactic radio sources 3C 345 VLBI component = radio blob superluminal motion = apparent motion VLBI observations of compact radio sources show that the ejection of
More informationASTRONOMY AND ASTROPHYSICS. Escape with the formation of a binary in two-dimensional three-body problem. I. Navin Chandra and K.B.
Astron. Astrophys. 46 65 66 1999 ASTRONOMY AND ASTROPHYSICS Escape with the formation of a binary in two-dimensional three-body problem. I Navin Chandra and K.B. Bhatnagar Centre for Fundamental Research
More informationWhy is the Universe Expanding?
Why is the Universe Expanding? In general relativity, mass warps space. Warped space makes matter move, which changes the structure of space. Thus the universe should be dynamic! Gravity tries to collapse
More informationEinstein s Equations. July 1, 2008
July 1, 2008 Newtonian Gravity I Poisson equation 2 U( x) = 4πGρ( x) U( x) = G d 3 x ρ( x) x x For a spherically symmetric mass distribution of radius R U(r) = 1 r U(r) = 1 r R 0 r 0 r 2 ρ(r )dr for r
More informationCentral force motion/kepler problem. 1 Reducing 2-body motion to effective 1-body, that too with 2 d.o.f and 1st order differential equations
Central force motion/kepler problem This short note summarizes our discussion in the lectures of various aspects of the motion under central force, in particular, the Kepler problem of inverse square-law
More informationQuantum Gravity with Linear Action and Gravitational Singularities. George Savvidy. Testing Fundamental Physical Principles
Quantum Gravity with Linear Action and Gravitational Singularities George Savvidy Demokritos National Research Centre Greece, Athens Testing Fundamental Physical Principles Corfu, September 22-28 (T he
More informationAnalytic methods in the age of numerical relativity
Analytic methods in the age of numerical relativity vs. Marc Favata Kavli Institute for Theoretical Physics University of California, Santa Barbara Motivation: Modeling the emission of gravitational waves
More informationAST1100 Lecture Notes
AST00 Lecture Notes Part E General Relativity: Gravitational lensing Questions to ponder before the lecture. Newton s law of gravitation shows the dependence of the gravitational force on the mass. In
More informationA choreographic solution to the general relativistic three-body problem. Abstract
A choreographic solution to the general relativistic three-body problem Tatsunori Imai, Takamasa Chiba, and Hideki Asada Faculty of Science and Technology, Hirosaki University, Hirosaki 036-8561, Japan
More informationarxiv: v1 [astro-ph.ep] 11 Aug 2016
Noname manuscript No. (will be inserted by the editor) Secular and tidal evolution of circumbinary systems. Alexandre C. M. Correia Gwenaël Boué Jacques Laskar arxiv:168.3484v1 [astro-ph.ep] 11 Aug 216
More informationA hybrid N-body code incorporating algorithmic regularization and post-newtonian forces
Mon. Not. R. Astron. Soc. 389, 1 (008) doi:10.1111/j.1365-966.008.13557.x A hybrid N-body code incorporating algorithmic regularization and post-newtonian forces S. Harfst, 1 A. Gualandris, D. Merritt
More informationDetectability of extrasolar debris. Mark Wyatt Institute of Astronomy, University of Cambridge
Detectability of extrasolar debris Mark Wyatt Institute of Astronomy, University of Cambridge Why image extrasolar debris? Emission spectrum shows dust thermal emission, used to infer radius of parent
More informationBlack Hole Physics via Gravitational Waves
Black Hole Physics via Gravitational Waves Image: Steve Drasco, California Polytechnic State University and MIT How to use gravitational wave observations to probe astrophysical black holes In my entire
More informationBlack Hole-Neutron Star Binaries in General Relativity. Thomas Baumgarte Bowdoin College
Black Hole-Neutron Star Binaries in General Relativity Thomas Baumgarte Bowdoin College Keisuke Taniguchi, Joshua Faber, Stu Shapiro University of Illinois Numerical Relativity Solve Einstein s equations
More informationTHE ORIGIN AND EVOLUTION OF FREE-FLOATING PLANETS IN STAR CLUSTERS
THE ORIGIN AND EVOLUTION OF FREE-FLOATING PLANETS IN STAR CLUSTERS M.B.N. (Thijs) Kouwenhoven Kavli Institute for Astronomy and Astrophysics, Peking University Hao Wei (MPIA), Li Yun (KIAA), Wang Long
More informationarxiv: v2 [astro-ph.he] 6 Aug 2018
Draft version August 7, 018 Preprint typeset using L A TEX style AASTeX6 v. 1.0 A TRIPLE ORIGIN FOR THE HEAVY AND LOW-SPIN BINARY BLACK HOLES DETECTED BY LIGO/VIRGO Carl L. Rodriguez 1 and Fabio Antonini
More informationPulsars as probes for the existence of IMBHs
Universidad de Valencia 15 November 2010 Pulsars as probes for the existence of IMBHs ANDREA POSSENTI Layout Known Black Hole classes Formation scenarios for the IMBHs IMBH candidates IMBH candidates (?)
More informationExoplanets: a dynamic field
Exoplanets: a dynamic field Alexander James Mustill Amy Bonsor, Melvyn B. Davies, Boris Gänsicke, Anders Johansen, Dimitri Veras, Eva Villaver The (transiting) exoplanet population Solar System Hot Jupiters:
More informationDynamics of spinning particles in Schwarzschild spacetime
Dynamics of spinning particles in Schwarzschild spacetime, Volker Perlick, Claus Lämmerzahl Center of Space Technology and Microgravity University of Bremen, Germany 08.05.2014 RTG Workshop, Bielefeld
More informationLECTURES on COLLISIONAL DYNAMICS: 4. HOT TOPICS on COLLISIONAL DYNAMICS
Michela Mapelli LECTURES on COLLISIONAL DYNAMICS: 4. HOT TOPICS on COLLISIONAL DYNAMICS Part 2 1) IMBHs: runaway collapse, repeated mergers,... 2) BHs eject each other? 3) Effects of 3-body on X-ray binaries
More informationBlack Holes ASTR 2110 Sarazin. Calculation of Curved Spacetime near Merging Black Holes
Black Holes ASTR 2110 Sarazin Calculation of Curved Spacetime near Merging Black Holes Test #2 Monday, November 13, 11-11:50 am Ruffner G006 (classroom) Bring pencils, paper, calculator You may not consult
More informationHigh Energy Astrophysics
High Energy Astrophysics Accretion Giampaolo Pisano Jodrell Bank Centre for Astrophysics - University of Manchester giampaolo.pisano@manchester.ac.uk April 01 Accretion - Accretion efficiency - Eddington
More informationBlack Hole Subsystems in Galactic Globular Clusters: Unravelling BH Populations in GCs using MOCCA Star Cluster Simulations
Black Hole Subsystems in Galactic Globular Clusters: Unravelling BH Populations in GCs using MOCCA Star Cluster Simulations Artist impression of a stellar mass BH detected in NGC 3201 (Giesers et al. 2018)
More informationLecture 13: Binary evolution
Lecture 13: Binary evolution Senior Astrophysics 2017-04-12 Senior Astrophysics Lecture 13: Binary evolution 2017-04-12 1 / 37 Outline 1 Conservative mass transfer 2 Non-conservative mass transfer 3 Cataclysmic
More informationDYNAMICS OF MIXED BINARIES
DYNAMICS OF MIXED BINARIES Luciano Rezzolla Albert Einstein Institute, Golm, Germany In collaboration with Frank Löffler & Marcus Ansorg [Phys. Rev. D 74 104018 (2006)] SISSA (Trieste, Italy), AEI (Golm,
More informationStructure of black holes in theories beyond general relativity
Structure of black holes in theories beyond general relativity Weiming Wayne Zhao LIGO SURF Project Caltech TAPIR August 18, 2016 Wayne Zhao (LIGO SURF) Structure of BHs beyond GR August 18, 2016 1 / 16
More information