emris kicks tails from black hole perturbation theory
play

EMRIs, Kicks & Tails from Black Hole Perturbation Theory (using - PowerPoint PPT Presentation

Jun 26, 2023 •324 likes •787 views

EMRIs, Kicks & Tails from Black Hole Perturbation Theory (using GPUs) Gaurav Khanna Associate Professor, UMass Dartmouth 7/29/15 AEI Golm 1 Collaborators & Support Alessandra Buonanno, AEI Golm Lior Burko, Georgia G. College

  1. EMRIs, Kicks & Tails from Black Hole Perturbation Theory (using GPUs) Gaurav Khanna Associate Professor, UMass Dartmouth 7/29/15 AEI Golm 1

  2. Collaborators & Support • Alessandra Buonanno, AEI Golm • Lior Burko, Georgia G. College • Scott Hughes, MIT • Richard Price, UTexas / UMassD • P. Sundararajan, MIT, Morgan-Stanley • Andrea Taracchini, AEI Golm • Anil Zenginoglu, UMaryland Funding: National Science Foundation 7/29/15 AEI Golm 2

  3. Talk Outline • EMRIs using BH perturbation theory • Teukolsky equation solver with particle-source in the time-domain • Waveforms and kicks & anti-kicks • Scalar-field “tails” in BH spacetimes • Kerr BH “tails” controversy • Computing with gaming devices • OpenCL / GPU computing • Summary and Future Outlook 7/29/15 AEI Golm 3

  4. Kerr (rotating) black hole perturbation theory • Write Einstein’s GR field equations to linear order expanding about a BH solution • Teukolsky equation -- a wave-equation like PDE that describes how generic fields (scalar, vector, tensor) in the space-time of a Kerr BH behave (propagate/evolve) • In the gravitational field context -- describes the behavior of GWs emitted from Kerr black holes • Relatively simple: linear, hyperbolic, (3+1)D PDE .. (can be reduced down to (2+1)D) 7/29/15 AEI Golm 4

  5. Teukolsky equation Here -- a: spin of the Kerr black hole M: mass of the Kerr black hole s: spin-weight of the field considered (s is -2 for outgoing GWs) 7/29/15 AEI Golm 5

  6. EMRI: Extreme-Mass Ratio (binary) Inspiral EMRI 7/29/15 AEI Golm 6

  7. Point-like compact object .. • Source of the GWs (perturbation) in EMRI is the inspiraling compact object • How to model a point-like compact object (technically a Dirac-delta function) on a numerical (finite-difference) grid? • Obvious approach would be to take a narrow Gaussian distribution; do several runs with successively narrower profiles; take some sort of a limit .. • Decent results, but very expensive 7/29/15 AEI Golm 7

  8. Discrete Dirac-delta • Alternative approach to representing Dirac- delta on numerical finite-difference grid • Particularly well suited for numerical computation (time-domain, finite-difference) • Work done in collaboration with Pranesh Sundararajan & Scott Hughes (MIT) • Excellent (highly accurate) results • Extremely efficient (by over an order-of- magnitude) 7/29/15 AEI Golm 8

  9. 7/29/15 AEI Golm 9

  10. Teukolsky Equation 7/29/15 AEI Golm 10

  11. 7/29/15 AEI Golm 11

  12. “source-term” T4 s4 = 1.D0/2.D0 s15 = -(r+cmplx(0.D0,1.D0)*a*ctheta)/(r**2+a**2*ctheta**2)*sqrt(2. s6 = 1/(r**2+a**2*ctheta**2) #D0)*(-a*stheta*nmu*(E*(a**2+rp**2)-a*lz)**2/dtdT/rp**6/pie**2/wr*e s9 = r**2+a**2 #xp(-(r-rp)**2/wr**2/2)/wt*exp(-ctheta**2/wt**2/2)*mm*dphidt*exp(cm s12 = (r+cmplx(0.D0,1.D0)*a*ctheta)/(r**2+a**2*ctheta**2)*sqrt(2.D #plx(0.D0,-1.D0)*mm*phip)/16+nmu*(E*(a**2+rp**2)-a*lz)**2/dtdT/rp** #0)*(a*stheta*nmu/rp**5*sqrt(2.D0)/dtdT*(E*(a**2+rp**2)-a*lz)*(a*E- #6/pie**2/wr*exp(-(r-rp)**2/wr**2/2)/wt**3*ctheta*stheta*exp(-cthet #lz)/pie**2/wr*exp(-(r-rp)**2/wr**2/2)/wt*exp(-ctheta**2/wt**2/2)*m #a**2/wt**2/2)*exp(cmplx(0.D0,-1.D0)*mm*phip)/16+1/stheta*nmu*(E*(a #m**2*dphidt**2*exp(cmplx(0.D0,-1.D0)*mm*phip)/16-nmu/rp**5*sqrt(2. #**2+rp**2)-a*lz)**2/dtdT/rp**6/pie**2/wr*exp(-(r-rp)**2/wr**2/2)/w #D0)/dtdT*(E*(a**2+rp**2)-a*lz)*(a*E-lz)/pie**2/wr*exp(-(r-rp)**2/w #t*exp(-ctheta**2/wt**2/2)*mm*exp(cmplx(0.D0,-1.D0)*mm*phip)/16)/2 #r**2/2)/wt**3*ctheta*stheta*exp(-ctheta**2/wt**2/2)*mm*dphidt*exp( s16 = -(cmplx(0.D0,-1.D0/2.D0)*a*(r+cmplx(0.D0,1.D0)*a*ctheta)/(r* #cmplx(0.D0,-1.D0)*mm*phip)/16-1/stheta*nmu/rp**5*sqrt(2.D0)/dtdT*( #*2+a**2*ctheta**2)**2*(r+cmplx(0.D0,-1.D0)*a*ctheta)*stheta*sqrt(2 #E*(a**2+rp**2)-a*lz)*(a*E-lz)/pie**2/wr*exp(-(r-rp)**2/wr**2/2)/wt #.D0)+cmplx(0.D0,1.D0)*a*(r+cmplx(0.D0,1.D0)*a*ctheta)**2/(r**2+a** #*exp(-ctheta**2/wt**2/2)*mm**2*dphidt*exp(cmplx(0.D0,-1.D0)*mm*phi #2*ctheta**2)**2*stheta*sqrt(2.D0))*nmu*(E*(a**2+rp**2)-a*lz)**2/dt #p)/16)/2 #dT/rp**6/pie**2/wr*exp(-(r-rp)**2/wr**2/2)/wt*exp(-ctheta**2/wt**2 s13 = cmplx(0.D0,1.D0/8.D0)*a*(r+cmplx(0.D0,1.D0)*a*ctheta)/(r**2+ #/2)*exp(cmplx(0.D0,-1.D0)*mm*phip)/16 #a**2*ctheta**2)**2*(r+cmplx(0.D0,-1.D0)*a*ctheta)*stheta*nmu/rp**5 s14 = s15+s16 #/dtdT*(E*(a**2+rp**2)-a*lz)*(a*E-lz)/pie**2/wr*exp(-(r-rp)**2/wr** s10 = s13+s14 #2/2)/wt*exp(-ctheta**2/wt**2/2)*mm*dphidt*exp(cmplx(0.D0,-1.D0)*mm s8 = s9*s10 #*phip)-(cmplx(0.D0,1.D0/2.D0)*a*(r+cmplx(0.D0,1.D0)*a*ctheta)**2/( s6 = s7*s8 #r**2+a**2*ctheta**2)**2*stheta*sqrt(2.D0)-(r+cmplx(0.D0,1.D0)*a*ct s8 = cmplx(0.D0,2.D0)*a*(r+cmplx(0.D0,1.D0)*a*ctheta)**2 #heta)/(r**2+a**2*ctheta**2)/tan(th)*sqrt(2.D0)/4)*nmu/rp**5*sqrt(2 s10 = 1/((r**2+a**2*ctheta**2)**2)*stheta #.D0)/dtdT*(E*(a**2+rp**2)-a*lz)*(a*E-lz)/pie**2/wr*exp(-(r-rp)**2/ s12 = sqrt(2.D0) #wr**2/2)/wt*exp(-ctheta**2/wt**2/2)*mm*dphidt*exp(cmplx(0.D0,-1.D0 s15 = 1/(r**2+a**2*ctheta**2)*(-(r**2+a**2)*nmu/rp**5*sqrt(2.D0)/d #)*mm*phip)/8 #tdT*(E*(a**2+rp**2)-a*lz)*(a*E-lz)/pie**2/wr*exp(-(r-rp)**2/wr**2/ s11 = s12+s13 #2)/wt*exp(-ctheta**2/wt**2/2)*mm*dphidt*exp(cmplx(0.D0,-1.D0)*mm*p s12 = s11-1/(r**2+a**2*ctheta**2)*((r**2+a**2)*nmu/rp**4/dtdT*(a*E #hip)/16+cmplx(0.D0,-1.D0/16.D0)*(r**2+a**2-2*M*r)*nmu/rp**5*sqrt(2 #-lz)**2/pie**2/wr*exp(-(r-rp)**2/wr**2/2)/wt*exp(-ctheta**2/wt**2/ #.D0)/dtdT*(E*(a**2+rp**2)-a*lz)*(a*E-lz)/pie**2/wr**3*(r-rp)*exp(- #2)*mm**2*dphidt**2*exp(cmplx(0.D0,-1.D0)*mm*phip)/8+cmplx(0.D0,1.D #(r-rp)**2/wr**2/2)/wt*exp(-ctheta**2/wt**2/2)*exp(cmplx(0.D0,-1.D0 #0/8.D0)*(r**2+a**2-2*M*r)*nmu/rp**4/dtdT*(a*E-lz)**2/pie**2/wr**3* #)*mm*phip)+a*nmu/rp**5*sqrt(2.D0)/dtdT*(E*(a**2+rp**2)-a*lz)*(a*E- #(r-rp)*exp(-(r-rp)**2/wr**2/2)/wt*exp(-ctheta**2/wt**2/2)*mm*dphid #lz)/pie**2/wr*exp(-(r-rp)**2/wr**2/2)/wt*exp(-ctheta**2/wt**2/2)*m #t*exp(cmplx(0.D0,-1.D0)*mm*phip)-a*nmu/rp**4/dtdT*(a*E-lz)**2/pie* #m*exp(cmplx(0.D0,-1.D0)*mm*phip)/16)/2 #*2/wr*exp(-(r-rp)**2/wr**2/2)/wt*exp(-ctheta**2/wt**2/2)*mm**2*dph s16 = cmplx(0.D0,-1.D0/8.D0)*(-(r+cmplx(0.D0,1.D0)*a*ctheta)**2/(r #idt*exp(cmplx(0.D0,-1.D0)*mm*phip)/8)/2 #**2+a**2*ctheta**2)**3*(r+cmplx(0.D0,-1.D0)*a*ctheta)*(r**2+a**2-2 s13 = s12+cmplx(0.D0,-1.D0/4.D0)*(-(r+cmplx(0.D0,1.D0)*a*ctheta)** #*M*r)/2+(r+cmplx(0.D0,1.D0)*a*ctheta)/(r**2+a**2*ctheta**2)**2*(r+ #2/(r**2+a**2*ctheta**2)**3*(r+cmplx(0.D0,-1.D0)*a*ctheta)*(r**2+a* #cmplx(0.D0,-1.D0)*a*ctheta)*(r-M)/2)*nmu/rp**5*sqrt(2.D0)/dtdT*(E* #*2-2*M*r)/2+(r+cmplx(0.D0,1.D0)*a*ctheta)/(r**2+a**2*ctheta**2)**2 #(a**2+rp**2)-a*lz)*(a*E-lz)/pie**2/wr*exp(-(r-rp)**2/wr**2/2)/wt*e #*(r+cmplx(0.D0,-1.D0)*a*ctheta)*(r-M)/2)*nmu/rp**4/dtdT*(a*E-lz)** #xp(-ctheta**2/wt**2/2)*exp(cmplx(0.D0,-1.D0)*mm*phip) #2/pie**2/wr*exp(-(r-rp)**2/wr**2/2)/wt*exp(-ctheta**2/wt**2/2)*mm* s14 = s15+s16 #dphidt*exp(cmplx(0.D0,-1.D0)*mm*phip) 7/29/15 AEI Golm 12

Download Presentation
Download Policy: The content available on the website is offered to you 'AS IS' for your personal information and use only. It cannot be commercialized, licensed, or distributed on other websites without prior consent from the author. To download a presentation, simply click this link. If you encounter any difficulties during the download process, it's possible that the publisher has removed the file from their server.

Recommend

Studying Black-hole kicks BH-LMXB The Stars for All Meeting observed pop. Supernova kicks

Studying Black-hole kicks BH-LMXB The Stars for All Meeting observed pop. Supernova kicks

BH kicks S. Repetto Issue Studying Black-hole kicks BH-LMXB The Stars for All Meeting observed pop. Supernova kicks Serena Repetto 1 , 2 , Melvyn B. Davies 1 and S. Sigurdsson 3 Galactic orbits BPS, How 1 Lund Observatory to 2 Dept. of

622 views • 20 slides
A Summary of the Black Hole Perturbation Theory Steven Hochman Introduction Many frameworks

A Summary of the Black Hole Perturbation Theory Steven Hochman Introduction Many frameworks

A Summary of the Black Hole Perturbation Theory Steven Hochman Introduction Many frameworks for doing perturbation theory The two most popular ones Direct examination of the Einstein equations -> Zerilli-Regge-Wheeler equations for

407 views • 30 slides
Lectures on black-hole perturbation theory Leor Barack University of Southampton Kavli-RISE

Lectures on black-hole perturbation theory Leor Barack University of Southampton Kavli-RISE

Lectures on black-hole perturbation theory Leor Barack University of Southampton Kavli-RISE Summer School on Gravitational Waves Kavli-RISE Summer School on GWs () Black-hole perturbation theory 1 / 73 Plan Overview Types of perturbative

1.22k views • 92 slides
( ) ( | z ) P ( z ) P Y P Y z 3 Bayes Rule Inference We will write

( ) ( | z ) P ( z ) P Y P Y z 3 Bayes Rule Inference We will write

Full Joint Probability Distributions Coin Card Candy P(Coin, Card, Candy) tails black 1 0.15 tails black 2 0.06 The probabilities CS 331: Artificial Intelligence in the last column tails black 3 0.09 sum to 1 tails red 1 0.02

71 views • 4 slides
( ) ( | z ) P ( z ) P Y P Y z 3 Inference Independence We will write the

( ) ( | z ) P ( z ) P Y P Y z 3 Inference Independence We will write the

Full Joint Probability Distributions Coin Card Candy P(Coin, Card, Candy) tails black 1 0.15 tails black 2 0.06 The probabilities CS 331: Artificial Intelligence in the last column tails black 3 0.09 sum to 1 tails red 1 0.02

153 views • 4 slides
CS 331: Artificial Intelligence in the last column tails black 3 0.09 sum to 1 tails red 1

CS 331: Artificial Intelligence in the last column tails black 3 0.09 sum to 1 tails red 1

Full Joint Probability Distributions Coin Card Candy P(Coin, Card, Candy) tails black 1 0.15 tails black 2 0.06 The probabilities CS 331: Artificial Intelligence in the last column tails black 3 0.09 sum to 1 tails red 1 0.02

459 views • 5 slides
Black Hole Perturbation Toolkit Barry Wardell University College Dublin Advances in Computational

Black Hole Perturbation Toolkit Barry Wardell University College Dublin Advances in Computational

Black Hole Perturbation Toolkit Barry Wardell University College Dublin Advances in Computational Relativity, bhptoolkit.org ICERM, 16 th September 2020 Workshop outline - Overview of the toolkit (~30 minutes) - Tutorial: Getting started with

734 views • 25 slides
Self-consistent orbital evolution of a particle around a Schwarzschild black hole Barry Wardell

Self-consistent orbital evolution of a particle around a Schwarzschild black hole Barry Wardell

Self-consistent orbital evolution of a particle around a Schwarzschild black hole Barry Wardell University College Dublin Collaborators: Peter Diener, Ian Vega, Steven Detweiler BritGrav Conference, Southampton, 4th April 2012 EMRIs and

457 views • 14 slides
Black hole Degeneracies from the effective action Guillaume Bossard CPHT, Ecole Polytechnique

Black hole Degeneracies from the effective action Guillaume Bossard CPHT, Ecole Polytechnique

Black hole Degeneracies from the effective action Guillaume Bossard CPHT, Ecole Polytechnique Texas A & M, April 2018 Outline Heterotic CHL strings and U-duality Protected couplings beyond perturbation theory Decompactification limit and

741 views • 35 slides
Constructing black holes and black hole microstates String theory and the fuzzball proposal Cl

Constructing black holes and black hole microstates String theory and the fuzzball proposal Cl

Introduction : black hole issues and entropy counting The fuzzball proposal Constructing three-charge supersymmetric solutions Non-BPS extremal black holes Conclusion and perspectives Constructing black holes and black hole microstates String

1.37k views • 95 slides
Electromagnetic-gravitational perturbations of Kerr-Newman black hole Elena Giorgi Gravity

Electromagnetic-gravitational perturbations of Kerr-Newman black hole Elena Giorgi Gravity

Kerr-Newman BH perturbation theory The Teukolsky equations The Regge-Wheeler equations Electromagnetic-gravitational perturbations of Kerr-Newman black hole Elena Giorgi Gravity Initiative, Princeton University October 5th, 2020 ICERM

823 views • 59 slides
Geometry of Higher-Dimensional Black Hole Thermodynamics (hep-th/0510139) Narit Pidokrajt (work

Geometry of Higher-Dimensional Black Hole Thermodynamics (hep-th/0510139) Narit Pidokrajt (work

Geometry of Higher-Dimensional Black Hole Thermodynamics (hep-th/0510139) Narit Pidokrajt (work with Jan E. Aman) Quantum and Field Theory Fysikum, Stockholms Universitet Abstract We apply the Ruppeiner theory to black hole

298 views • 14 slides
CNN Black Hole-based 'Black Holes' In Ocean Exist, Malaysian airplane theory' Scientists Say

CNN Black Hole-based 'Black Holes' In Ocean Exist, Malaysian airplane theory' Scientists Say

CNN Black Hole-based 'Black Holes' In Ocean Exist, Malaysian airplane theory' Scientists Say The many faces of BHs Astrophysics Gravitation Supermassive BHs BHs are elementary particles of gravity final stage of stellar collapse

872 views • 66 slides
Shaving the Black Hole Yogesh K. Srivastava Work with Dileep Jatkar and Ashoke Sen KEK

Shaving the Black Hole Yogesh K. Srivastava Work with Dileep Jatkar and Ashoke Sen KEK

Prologue String Theory Precision counting of microstates Black Hole Hair Analysis of the BMPV BH Entropy Analysis of the 4D BH Entropy Hair modes in Supergravity Regularity of Hair Modes Conclusion Shaving the Black Hole Yogesh K.

876 views • 54 slides
Black hole X-ray binaries V: Formation and evolution of black hole binaries Thomas J. Maccarone

Black hole X-ray binaries V: Formation and evolution of black hole binaries Thomas J. Maccarone

Black hole X-ray binaries V: Formation and evolution of black hole binaries Thomas J. Maccarone Jan 14, 2017, 2nd Fudan Winter School on Black Hole Astrophysics Overview Phenomena associated with binary evolution Formation of close binaries

651 views • 30 slides
A Conceptual Analysis of Black Hole Entropy in String Theory Sebastian De Haro University of

A Conceptual Analysis of Black Hole Entropy in String Theory Sebastian De Haro University of

A Bit of String Theory The Strominger-Vafa Calculation Further questions about the Strominger-Vafa Black Hole A Conceptual Analysis of Black Hole Entropy in String Theory Sebastian De Haro University of Amsterdam and University of Cambridge

592 views • 47 slides
Simulations at the nanoscale on the GRID using Quantum ESPRESSO P. Giannozzi Universit` a di

Simulations at the nanoscale on the GRID using Quantum ESPRESSO P. Giannozzi Universit` a di

Simulations at the nanoscale on the GRID using Quantum ESPRESSO P. Giannozzi Universit` a di Udine and Democritos CNR-IOM Trieste, Italy Hands on Training School on Molecular and Material Science GRID Applications, Trieste, 2010/03/31

494 views • 30 slides
Physical predictions from lattice QCD Christian Hoelbling Bergische Universitt Wuppertal

Physical predictions from lattice QCD Christian Hoelbling Bergische Universitt Wuppertal

Introduction Hadron spectrum Quark masses etc. Summary Physical predictions from lattice QCD Christian Hoelbling Bergische Universitt Wuppertal Budapest: S. Katz Marseille: L. Lellouch, A. Portelli, A. Sastre Wuppertal: Sz. Borsanyi, S.

1.03k views • 81 slides
Constraints on the equation of state of dense matter from experiments and observations

Constraints on the equation of state of dense matter from experiments and observations

Constraints on the equation of state of dense matter from experiments and observations Constana Providncia CFisUC, Universidade de Coimbra, Portugal NewCompStar 2017, 31 March, 2017 Outline Constraints from Heavy ion Collisions

848 views • 35 slides
Autodifferentiation CMSC 678 UMBC Recap from last time Maximum Entropy (Log-linear) Models

Autodifferentiation CMSC 678 UMBC Recap from last time Maximum Entropy (Log-linear) Models

Neural Networks and Autodifferentiation CMSC 678 UMBC Recap from last time Maximum Entropy (Log-linear) Models ) exp( , ) model the posterior probabilities of the K classes via linear functions

735 views • 69 slides
WELCOME TO THE 4% Code Liberation x Indiecade 2014 A trans-inclusive, women-only pro g rammin g

WELCOME TO THE 4% Code Liberation x Indiecade 2014 A trans-inclusive, women-only pro g rammin g

WELCOME TO THE 4% Code Liberation x Indiecade 2014 A trans-inclusive, women-only pro g rammin g workshop 1 WHO ARE YOU? 2 WHO AM I? 3 Jane Friedhoff Creative Coder & Game Designer janefriedhoff.com @jfriedhoff c/o

1.86k views • 153 slides
Lecture no: 7 Channel Coding Ove Edfors, Department of Electrical and Information Technology

Lecture no: 7 Channel Coding Ove Edfors, Department of Electrical and Information Technology

RADIO SYSTEMS ETIN15 Lecture no: 7 Channel Coding Ove Edfors, Department of Electrical and Information Technology Ove.Edfors@eit.lth.se 2012-04-23 Ove Edfors - ETIN15 1 Contents (CHANNEL CODING) Overview Block codes

580 views • 34 slides
Computing and Communications 1. Introduction Ying Cui Department of Electronic Engineering

Computing and Communications 1. Introduction Ying Cui Department of Electronic Engineering

1896 1920 1987 2006 Computing and Communications 1. Introduction Ying Cui Department of Electronic Engineering Shanghai Jiao Tong University, China 2017, Autumn 1 COURSE INFORMATION 2 Lecture Time: Monday 8:00-10:00am, Sep 11-Dec 25

1.18k views • 38 slides
Relative generalized Hamming weights of one-point algebraic geometric codes: an application to

Relative generalized Hamming weights of one-point algebraic geometric codes: an application to

Relative generalized Hamming weights of one-point algebraic geometric codes: an application to secret sharing INdAM meeting: International meeting on numerical semigroups Cortona 2014, September 10th. Diego Ruano

589 views • 57 slides

More recommend