black hole universe -bh in an expanding box- yoo, chulmoon ( yitp) hiroyuki abe (osaka city univ.)...
DESCRIPTION
Chulmoon Yoo 3 Lattice Universe “Dynamics of a Lattice Universe by the Schwarzschild-Cell Method” [Lindquist and Wheeler(1957)] The maximum radius asymptotically agrees with the dust universe case Putting N equal mass Sch. BHs on a 3-sphere, requiring a matching condition, we get a dynamics of the lattice universe maximum radius of lattice universe number of BHsTRANSCRIPT
![Page 1: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/1.jpg)
Black Hole Universe-BH in an expanding box-
Yoo, Chulmoon ( YITP)
Hiroyuki Abe (Osaka City Univ.)Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka City Univ.)
![Page 2: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/2.jpg)
Chulmoon Yoo
2Cluster of Many BHs ~ Dust Fluids?
Naively thinking, we can treat the cluster of a number of BHs as a dust fluid on average
In this work, as a simplest case, we try to construct “the BH universe” which would be approximated by the EdS universe on average
But, it is very difficult to show it from the first principle. Because we need to solve the N-body dynamics with the Einstein equations.
dust fluid~~
![Page 3: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/3.jpg)
Chulmoon Yoo
3
Lattice Universe“Dynamics of a Lattice Universe by the Schwarzschild-Cell Method” [Lindquist and Wheeler(1957)]
The maximum radius asymptotically agrees with the dust universe case
Putting N equal mass Sch. BHs on a 3-sphere, requiring a matching condition, we get a dynamics of the lattice universe
maximum radius of lattice universe
number of BHs
![Page 4: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/4.jpg)
Chulmoon Yoo
4
Swiss-cheese Universe
Expand
Homogeneous dust universe
Cutting spherical regions, put Schwarzschild BHs with the same mass
Swiss-cheese universe
We want to make it without cheese(“Swiss universe” ?)
![Page 5: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/5.jpg)
Chulmoon Yoo
5
Some Aspects of This Work
If perturbations of metric components are small enough, we don’t need to treat full GR but perturbation theory is applicable. Perhaps, even if the density perturbation is nonlinear in small scales, we could handle the inhomogeneities without full numerical relativity.
1. “Cosmological Numerical Relativity (CNR)”In which situation, CNR may be significant?
(In this sense, for late time cosmology, CNR might not be significant.)CNR may play a role in an extreme situation where the metric perturbation is full nonlinear on cosmological scales (e.g. primordial BH formation)
2. BH simulation without asymptotic flatness-In higher-dimensional theory, compactified directions often exist, and they are not asymptotically flat. -BH physics might be applied to other fields (e.g. AdS/CFT,QCD,CMP) without asymptotic flatness
Their dynamical simulations might have common feature?
![Page 6: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/6.jpg)
Chulmoon Yoo
6
Contents
◎Part 1 : “A recipe for the BH universe”How to construct the initial data for the BH universe
◎Part 2 “ Structure of the BH universe”- Horizons- Effective Hubble equation with an averaging
![Page 7: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/7.jpg)
Chulmoon Yoo
7
Part 1A recipe for the BH
universe
![Page 8: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/8.jpg)
Chulmoon Yoo
8
…
What We Want to Do
◎Vacuum solution for the Einstein eqs.
First, we construct the puncture initial data
◎Expansion of the universe is crucial to avoid the potential divergence
Periodic boundary
ExpandingBH…
…
…
![Page 9: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/9.jpg)
Chulmoon Yoo
9
PunctureBoundary
Infinity of the other world
![Page 10: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/10.jpg)
Chulmoon Yoo
10
Constraint Eqs.
We assume
Setting trK by hand, we solve these eqs.How should we choose trK?
We construct the initial data.
where
![Page 11: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/11.jpg)
Chulmoon Yoo
11Expansion of the universe
tr K must be a finite value around the boundary
Expand
finite Hubble parameter HH =-tr K / 3
→Swiss-cheese case
![Page 12: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/12.jpg)
Chulmoon Yoo
12
CMC (constant mean curvature) Slicetr K = const. ⇔ ∇ana=const.
induced metric
isotropic coordinate
CMC slice?
![Page 13: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/13.jpg)
Chulmoon Yoo
13
r=∞ r=∞
R=Rc
For K≠0, we have a finite R at r=∞We need to take care of the inner boundaryTo avoid this, we choose K=0 near the infinity(maximal slice)
r=∞R=0
Difficulty to use CMC slice
![Page 14: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/14.jpg)
Chulmoon Yoo
14
trK
CMC sliceMaximal slice
trK
/Kc
R
![Page 15: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/15.jpg)
Chulmoon Yoo
15
Constraint Eqs.
Extraction of 1/R divergence
Near the center R=0 (trK=0)
ψ is regular at R=0
Periodic boundary condition for ψ and Xi
1
* f=0 at the boundary
r=∞
R=0
![Page 16: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/16.jpg)
Chulmoon Yoo
16
Equations
xy
z
L
R:=(x2+y2+z2)1/2
Source terms must vanish by integrating in a box Poisson equation with periodic boundary condition
![Page 17: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/17.jpg)
Chulmoon Yoo
17
Integration of source terms
vanishes by integrating in the box because ∂x Z and ∂x K are odd function of x
Vanishes by integrating in the box because K=const. at the boundary
Integration of this part also must vanish
![Page 18: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/18.jpg)
Chulmoon Yoo
18Effective Hubble Equation
Integrating in a box, we have
Hubble parameter H
effective mass density
![Page 19: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/19.jpg)
Chulmoon Yoo
19
Parameters•BH mass•Box size (isotropic coord.)
•Hubble radius
We set Kc so that the following equation is satisfied
This is just the integration of the constraint equation. We update the value of Kc at each step of the numerical iteration.
Free parameter is only
other than and
![Page 20: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/20.jpg)
Chulmoon Yoo
20
Part 2Structure of the BH
universe
![Page 21: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/21.jpg)
Chulmoon Yoo
21
trK
trK
/Kc
R
0.1
L-0.1
![Page 22: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/22.jpg)
Chulmoon Yoo
22
Numerical Solutions(1)
xy
z
L
ψ(x,y,L) for L=2M
ψ(x,y,0) for L=2M
![Page 23: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/23.jpg)
Chulmoon Yoo
23
Numerical Solutions(2)Z(x,y,L) for L=2M
Z(x,y,0) for L=2M
xy
z
L
![Page 24: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/24.jpg)
Chulmoon Yoo
24
Numerical Solutions(3)Xx(x,y,L) for L=2M
Xx(x,y,0) for L=2M
xy
z
L
![Page 25: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/25.jpg)
Chulmoon Yoo
25
Convergence Test
◎Beautiful quadratic convergence!◎We cannot find the solution for L<1.4M
![Page 26: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/26.jpg)
Chulmoon Yoo
26
Horizons ◎To see Horizons, we calculate outgoing(+) and ingoing(-) null expansions of spheres
◎We plot the value of χ for three independent directions (χ is not spherically symmetric in general)
: unit normal vector to sphere
◎Horizons (approximate position): Black hole horizon: White hole horizon
![Page 27: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/27.jpg)
Chulmoon Yoo
27
Expansion ◎parameter : L=1.4M
χ+
χ-
◎Horizons are almost spherically symmetric◎BH horizon exists outside WH horizon or they are almost identical
R
expa
nsio
n
WH BH
![Page 28: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/28.jpg)
Chulmoon Yoo
28
Time slice
BH horizon
“WH horizon”
Bifurcation point
![Page 29: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/29.jpg)
Chulmoon Yoo
29
Inhomogeneity
(x,y,L) for L=2M
(x,y,0) for L=2M
xy
z
L
◎Square of the traceless part of 3-dim Ricci curvature
homogeneous ⇒
homogeneous and empty⇒Milne universe (ΩK=1)
![Page 30: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/30.jpg)
Chulmoon Yoo
30
Inhomogeneity(x,y,L) for L=2M
(x,y,L) for L=4M(x,y,L) for L=5M0.6
0.6
0.7
Not homogeneous around the center of a boundary face
![Page 31: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/31.jpg)
Chulmoon Yoo
31
An Averaging◎Effective density
xy
z
L
Area: Effective volume of a box ( )
Effective density
◎Hubble parameter (defined by the boundary value of trK)
◎We may expect (?)This relation is
nontrivial!No dust, No matter, No symmetry, but additional gravitational energy other than “the point mass”
![Page 32: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/32.jpg)
Chulmoon Yoo
32
Effective Hubble◎Effective Hubble parameter
◎It asymptotically agrees with the expected value!
![Page 33: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/33.jpg)
Chulmoon Yoo
33
Conclusion◎We constructed initial data for the BH universe
◎BH horizon exists outside WH horizon or they are almost identical
◎When the box size is sufficiently larger than the Schwarzschild radius of the mass M, an effective density and an effective Hubble parameter satisfy Hubble equation of the EdS universe, that is, the BH universe is the EdS universe on Average!
◎Around vertices, it is well described by the Milne universe
![Page 34: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/34.jpg)
Chulmoon Yoo
34
Including Λ<0◎We may have momentarily static initial data
integrate in a box
◎Probably, It will collapse when we consider the time evolution because essentially it is dust + negative Λ universe
◎Negative Λ can compensate the mass term
◎It seems very difficult to get stable solution without exotic matter other than negative Λ
![Page 35: Black Hole Universe -BH in an expanding box- Yoo, Chulmoon ( YITP) Hiroyuki Abe (Osaka City Univ.) Ken-ichi Nakao (Osaka City Univ.) Yohsuke Takamori (Osaka](https://reader036.vdocument.in/reader036/viewer/2022062503/5a4d1aeb7f8b9ab05997aae7/html5/thumbnails/35.jpg)
Chulmoon Yoo
35
Thank you very much!