Solving Einstein s Equation Numerically III
|
|
- Scot Kenneth Marsh
- 6 years ago
- Views:
Transcription
1 Solving Einstein s Equation Numerically III Lee Lindblom Center for Astrophysics and Space Sciences University of California at San Diego Mathematical Sciences Center Lecture Series Tsinghua University 21 November 2014 Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 1 / 17
2 Gauge and Hyperbolicity in Electromagnetism The usual representation of the vacuum Maxwell equations split into evolution equations and constraints: t E = B, E = 0, t B = E, B = 0. These equations are often written in the more compact 4-dimensional form a F ab = 0 and [a F bc] = 0, where F ab has components E and B. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 2 / 17
3 Gauge and Hyperbolicity in Electromagnetism The usual representation of the vacuum Maxwell equations split into evolution equations and constraints: t E = B, E = 0, t B = E, B = 0. These equations are often written in the more compact 4-dimensional form a F ab = 0 and [a F bc] = 0, where F ab has components E and B. Maxwell s equations can be solved in part by introducing a vector potential F ab = a A b b A a. This reduces the system to the single equation: a a A b b a A a = 0. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 2 / 17
4 Gauge and Hyperbolicity in Electromagnetism The usual representation of the vacuum Maxwell equations split into evolution equations and constraints: t E = B, E = 0, t B = E, B = 0. These equations are often written in the more compact 4-dimensional form a F ab = 0 and [a F bc] = 0, where F ab has components E and B. Maxwell s equations can be solved in part by introducing a vector potential F ab = a A b b A a. This reduces the system to the single equation: a a A b b a A a = 0. This form of the equations can be made manifestly hyperbolic by choosing the gauge correctly, e.g., let a A a = H(x, t, A), giving: a a A b = ( 2 t + 2 x + 2 y + 2 z ) Ab = b H. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 2 / 17
5 Gauge and Hyperbolicity in General Relativity The spacetime Ricci curvature tensor can be written as: R ab = 1 2 ψcd c d ψ ab + (a Γ b) + Q ab (ψ, ψ), where ψ ab is the 4-metric, and Γ a = ψ ad ψ bc Γ d bc. Like Maxwell s equations, these equation can not be solved without specifying suitable gauge conditions. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 3 / 17
6 Gauge and Hyperbolicity in General Relativity The spacetime Ricci curvature tensor can be written as: R ab = 1 2 ψcd c d ψ ab + (a Γ b) + Q ab (ψ, ψ), where ψ ab is the 4-metric, and Γ a = ψ ad ψ bc Γ d bc. Like Maxwell s equations, these equation can not be solved without specifying suitable gauge conditions. The gauge freedom in general relativity theory is the freedom to represent the equations using any coordinates x a on spacetime. Solving the equations requires some specific choice of coordinates be made. Gauge conditions are used to impose the desired choice. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 3 / 17
7 Gauge and Hyperbolicity in General Relativity The spacetime Ricci curvature tensor can be written as: R ab = 1 2 ψcd c d ψ ab + (a Γ b) + Q ab (ψ, ψ), where ψ ab is the 4-metric, and Γ a = ψ ad ψ bc Γ d bc. Like Maxwell s equations, these equation can not be solved without specifying suitable gauge conditions. The gauge freedom in general relativity theory is the freedom to represent the equations using any coordinates x a on spacetime. Solving the equations requires some specific choice of coordinates be made. Gauge conditions are used to impose the desired choice. One way to impose the needed gauge conditions is to specify H a, the source term for a wave equation for each coordinate x a : H a = c c x a = ψ bc ( b c x a Γ e bc e x a ) = Γ a, where Γ a = ψ bc Γ a bc and ψ ab is the 4-metric. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 3 / 17
8 Gauge Conditions in General Relativity Specifying coordinates by the generalized harmonic (GH) method is accomplished by choosing a gauge-source function H a (x, ψ), e.g. H a = ψ ab H b (x), and requiring that H a (x, ψ) = Γ a = 1 2 ψad ψ bc ( b ψ dc + c ψ db d ψ bc ). Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 4 / 17
9 Gauge Conditions in General Relativity Specifying coordinates by the generalized harmonic (GH) method is accomplished by choosing a gauge-source function H a (x, ψ), e.g. H a = ψ ab H b (x), and requiring that H a (x, ψ) = Γ a = 1 2 ψad ψ bc ( b ψ dc + c ψ db d ψ bc ). Recall that the spacetime Ricci tensor is given by R ab = 1 2 ψcd c d ψ ab + (a Γ b) + Q ab (ψ, ψ). The Generalized Harmonic Einstein equation is obtained by replacing Γ a = ψ ab Γ b with H a (x, ψ) = ψ ab H b (x, ψ): R ab (a [ Γb) + H b) ] = 1 2 ψcd c d ψ ab (a H b) + Q ab (ψ, ψ). Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 4 / 17
10 Gauge Conditions in General Relativity Specifying coordinates by the generalized harmonic (GH) method is accomplished by choosing a gauge-source function H a (x, ψ), e.g. H a = ψ ab H b (x), and requiring that H a (x, ψ) = Γ a = 1 2 ψad ψ bc ( b ψ dc + c ψ db d ψ bc ). Recall that the spacetime Ricci tensor is given by R ab = 1 2 ψcd c d ψ ab + (a Γ b) + Q ab (ψ, ψ). The Generalized Harmonic Einstein equation is obtained by replacing Γ a = ψ ab Γ b with H a (x, ψ) = ψ ab H b (x, ψ): R ab (a [ Γb) + H b) ] = 1 2 ψcd c d ψ ab (a H b) + Q ab (ψ, ψ). The vacuum GH Einstein equation, R ab = 0 with Γ a + H a = 0, is therefore manifestly hyperbolic, having the same principal part as the scalar wave equation: 0 = a a ψ = ψ ab a b ψ + Q( ψ). Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 4 / 17
11 ADM 3+1 Approach to Fixing Coordinates Coordinates must be chosen to label points in spacetime before the Einstein equations can be solved. For some purposes it is convenient to split the spacetime coordinates x a into separate time and space components: x a = {t, x i }. (t + δt, x k ) Construct spacetime foliation by spacelike slices. Choose time function with t = const. on these slices. Choose spatial coordinates, x k, on each slice. t = τ (t, x k ) k t Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 5 / 17
12 ADM 3+1 Approach to Fixing Coordinates Coordinates must be chosen to label points in spacetime before the Einstein equations can be solved. For some purposes it is convenient to split the spacetime coordinates x a into separate time and space components: x a = {t, x i }. (t + δt, x k ) Construct spacetime foliation by spacelike slices. Choose time function with t = const. on these slices. t = τ Choose spatial coordinates, x k, on each slice. (t, x k ) Decompose the 4-metric ψ ab into its 3+1 parts: ds 2 = ψ ab dx a dx b = N 2 dt 2 + g ij (dx i + N i dt)(dx j + N j dt). k t Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 5 / 17
13 ADM 3+1 Approach to Fixing Coordinates Coordinates must be chosen to label points in spacetime before the Einstein equations can be solved. For some purposes it is convenient to split the spacetime coordinates x a into separate time and space components: x a = {t, x i }. (t + δt, x k ) Construct spacetime foliation by spacelike slices. Choose time function with t = const. on these slices. t = τ Choose spatial coordinates, x k, on each slice. (t, x k ) Decompose the 4-metric ψ ab into its 3+1 parts: ds 2 = ψ ab dx a dx b = N 2 dt 2 + g ij (dx i + N i dt)(dx j + N j dt). The unit vector t a normal to the t =constant slices depends only on the lapse N and shift N i : t = τ = x a τ a = 1 N t Nk N k. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 5 / 17 k t
14 ADM Approach to the Einstein Evolution System Decompose the Einstein equations R ab = 0 using the ADM 3+1 coordinate splitting. The resulting system includes evolution equations for the spatial metric g ij and extrinsic curvature K ij : t g ij N k k g ij = 2NK ij + g jk i N k + g ik j N k, t K ij N k k K ij = NR (3) ij + K jk i N k + K ik j N k i j N 2NK ik K k j + NK k kk ij. The resulting system also includes constraints: 0 = R (3) K ij K ij + (K k k) 2, 0 = k K ki i K k k. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 6 / 17
15 ADM Approach to the Einstein Evolution System Decompose the Einstein equations R ab = 0 using the ADM 3+1 coordinate splitting. The resulting system includes evolution equations for the spatial metric g ij and extrinsic curvature K ij : t g ij N k k g ij = 2NK ij + g jk i N k + g ik j N k, t K ij N k k K ij = NR (3) ij + K jk i N k + K ik j N k i j N 2NK ik K k j + NK k kk ij. The resulting system also includes constraints: 0 = R (3) K ij K ij + (K k k) 2, 0 = k K ki i K k k. System includes no evolution equations for lapse N or shift N i. These quanties can be specified freely to fix the gauge. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 6 / 17
16 ADM Approach to the Einstein Evolution System Decompose the Einstein equations R ab = 0 using the ADM 3+1 coordinate splitting. The resulting system includes evolution equations for the spatial metric g ij and extrinsic curvature K ij : t g ij N k k g ij = 2NK ij + g jk i N k + g ik j N k, t K ij N k k K ij = NR (3) ij + K jk i N k + K ik j N k i j N 2NK ik K k j + NK k kk ij. The resulting system also includes constraints: 0 = R (3) K ij K ij + (K k k) 2, 0 = k K ki i K k k. System includes no evolution equations for lapse N or shift N i. These quanties can be specified freely to fix the gauge. Resolving the issues of hyperbolicity (i.e. well posedness of the initial value problem) and constraint stability are much more complicated in this approach. The most successful version is the BSSN evolution system used by many (most) codes. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 6 / 17
17 Dynamical GH Gauge Conditions The spacetime coordinates x b are fixed in the generalized harmonic Einstein equations by specifying H b : a a x b H b. The generalized harmonic Einstein equations remain hyperbolic as long as the gauge source functions H b are taken to be functions of the coordinates x b and the spacetime metric ψ ab. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 7 / 17
18 Dynamical GH Gauge Conditions The spacetime coordinates x b are fixed in the generalized harmonic Einstein equations by specifying H b : a a x b H b. The generalized harmonic Einstein equations remain hyperbolic as long as the gauge source functions H b are taken to be functions of the coordinates x b and the spacetime metric ψ ab. The simplest choice H b = 0 (harmonic gauge) fails for very dynamical spacetimes, like binary black hole mergers. This failure seems to occur because the coordinates themselves become very dynamical solutions of the wave equation a a x b = 0 in these situations. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 7 / 17
19 Dynamical GH Gauge Conditions The spacetime coordinates x b are fixed in the generalized harmonic Einstein equations by specifying H b : a a x b H b. The generalized harmonic Einstein equations remain hyperbolic as long as the gauge source functions H b are taken to be functions of the coordinates x b and the spacetime metric ψ ab. The simplest choice H b = 0 (harmonic gauge) fails for very dynamical spacetimes, like binary black hole mergers. This failure seems to occur because the coordinates themselves become very dynamical solutions of the wave equation a a x b = 0 in these situations. Another simple choice keeping H b fixed in the co-moving frame of the black holes also works well during the long inspiral phase, but fails when the black holes begin to merge. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 7 / 17
20 Dynamical GH Gauge Conditions II Some of the extraneous gauge dynamics could be removed by adding a damping term to the harmonic gauge condition: a a x b = H b = µt a a x b = µt b = µnψ tb. This works well for the spatial coordinates x i, driving them toward solutions of the spatial Laplace equation on the timescale 1/µ. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 8 / 17
21 Dynamical GH Gauge Conditions II Some of the extraneous gauge dynamics could be removed by adding a damping term to the harmonic gauge condition: a a x b = H b = µt a a x b = µt b = µnψ tb. This works well for the spatial coordinates x i, driving them toward solutions of the spatial Laplace equation on the timescale 1/µ. For the time coordinate t, this damped wave condition drives t to a time independent constant, which is not a good coordinate. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 8 / 17
22 Dynamical GH Gauge Conditions II Some of the extraneous gauge dynamics could be removed by adding a damping term to the harmonic gauge condition: a a x b = H b = µt a a x b = µt b = µnψ tb. This works well for the spatial coordinates x i, driving them toward solutions of the spatial Laplace equation on the timescale 1/µ. For the time coordinate t, this damped wave condition drives t to a time independent constant, which is not a good coordinate. A better choice sets t a H a = µ log g/n 2. The gauge condition in this case becomes t a a log g/n 2 = µ log g/n 2 + N 1 k N k This coordinate condition keeps g/n 2 close to unity, even during binary black hole mergers (where it became of order 100 using simpler gauge conditions). Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 8 / 17
23 The Constraint Problem Fixing the gauge in an appropriate way makes the Einstein equations hyperbolic, so the initial value problem becomes well-posed mathematically. In a well-posed representation, the constraints, C = 0, remain satisfied for all time if they are satisfied initially. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 9 / 17
24 The Constraint Problem Fixing the gauge in an appropriate way makes the Einstein equations hyperbolic, so the initial value problem becomes well-posed mathematically. In a well-posed representation, the constraints, C = 0, remain satisfied for all time if they are satisfied initially. There is no guarantee, however, that constraints that are small initially will remain small. Constraint violating instabilities were one of the major problems that made progress on solving the binary black hole problem so slow. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 9 / 17
25 The Constraint Problem Fixing the gauge in an appropriate way makes the Einstein equations hyperbolic, so the initial value problem becomes well-posed mathematically. In a well-posed representation, the constraints, C = 0, remain satisfied for all time if they are satisfied initially. There is no guarantee, however, that constraints that are small initially will remain small. Constraint violating instabilities were one of the major problems that made progress on solving the binary black hole problem so slow. Special representations of the Einstein equations are needed that control the growth of any constraint violations. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 9 / 17
26 Constraint Damping in Electromagnetism Electromagnetism is described by the hyperbolic evolution equation a a A b = b H. Are there any constraints? Where have the usual E = B = 0 constraints gone? Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 10 / 17
27 Constraint Damping in Electromagnetism Electromagnetism is described by the hyperbolic evolution equation a a A b = b H. Are there any constraints? Where have the usual E = B = 0 constraints gone? Gauge condition becomes a constraint: 0 = C b A b H. Maxwell s equations imply that this constraint is preserved: a a ( b A b H) = a a C = 0. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 10 / 17
28 Constraint Damping in Electromagnetism Electromagnetism is described by the hyperbolic evolution equation a a A b = b H. Are there any constraints? Where have the usual E = B = 0 constraints gone? Gauge condition becomes a constraint: 0 = C b A b H. Maxwell s equations imply that this constraint is preserved: a a ( b A b H) = a a C = 0. Modify evolution equations by adding multiples of the constraints: a a A b = b H+γ 0 t b C = b H+γ 0 t b ( a A a H). These changes effect the constraint evolution equation, a a C γ 0 t b b C = 0, so constraint violations are damped when γ 0 > 0. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 10 / 17
29 Constraints in the GH Evolution System The GH evolution system has the form, 0 = R ab (a Γ b) (a H b), = R ab (a C b), where C a = H a + Γ a plays the role of a constraint. Without constraint damping, these equations are very unstable to constraint violating instabilities. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 11 / 17
30 Constraints in the GH Evolution System The GH evolution system has the form, 0 = R ab (a Γ b) (a H b), = R ab (a C b), where C a = H a + Γ a plays the role of a constraint. Without constraint damping, these equations are very unstable to constraint violating instabilities. Imposing coordinates using a GH gauge function profoundly changes the constraints. The GH constraint, C a = 0, where C a = H a + Γ a, depends only on first derivatives of the metric. The standard Hamiltonian and momentum constraints, M a = 0, are determined by derivatives of the gauge constraint C a : M a [ R ab 1 2 ψ abr ] t b = [ (a C b) 12 ψ ab c C c ] t b. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 11 / 17
31 Constraint Damping Generalized Harmonic System Pretorius (based on a suggestion from Gundlach, et al.) modified the GH system by adding terms proportional to the gauge constraints: 0 = R ab (a C b) + γ 0 [ t(a C b) 1 2 ψ ab t c C c ], where t a is a unit timelike vector field. Since C a = H a + Γ a depends only on first derivatives of the metric, these additional terms do not change the hyperbolic structure of the system. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 12 / 17
32 Constraint Damping Generalized Harmonic System Pretorius (based on a suggestion from Gundlach, et al.) modified the GH system by adding terms proportional to the gauge constraints: 0 = R ab (a C b) + γ 0 [ t(a C b) 1 2 ψ ab t c C c ], where t a is a unit timelike vector field. Since C a = H a + Γ a depends only on first derivatives of the metric, these additional terms do not change the hyperbolic structure of the system. Evolution of the constraints C a follow from the Bianchi identities: 0 = c c C a 2γ 0 c[ t (c C a) ] + C c (c C a) 1 2 γ 0 t a C c C c. This is a damped wave equation for C a, that drives all small short-wavelength constraint violations toward zero as the system evolves (for γ 0 > 0). Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 12 / 17
33 Summary of the GH Einstein System Choose coordinates by fixing a gauge-source function H a (x, ψ), e.g. H a = ψ ab H b (x), and requiring that H a (x, ψ) = c c x a = Γ a = 1 2 ψad ψ bc ( b ψ dc + c ψ db d ψ bc ). Gauge condition H a = Γ a is a constraint: C a = H a + Γ a = 0. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 13 / 17
34 Summary of the GH Einstein System Choose coordinates by fixing a gauge-source function H a (x, ψ), e.g. H a = ψ ab H b (x), and requiring that H a (x, ψ) = c c x a = Γ a = 1 2 ψad ψ bc ( b ψ dc + c ψ db d ψ bc ). Gauge condition H a = Γ a is a constraint: C a = H a + Γ a = 0. Principal part of evolution system becomes manifestly hyperbolic: R ab (a C b) = 1 2 ψcd c d ψ ab (a H b) + Q ab (ψ, ψ). Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 13 / 17
35 Summary of the GH Einstein System Choose coordinates by fixing a gauge-source function H a (x, ψ), e.g. H a = ψ ab H b (x), and requiring that H a (x, ψ) = c c x a = Γ a = 1 2 ψad ψ bc ( b ψ dc + c ψ db d ψ bc ). Gauge condition H a = Γ a is a constraint: C a = H a + Γ a = 0. Principal part of evolution system becomes manifestly hyperbolic: R ab (a C b) = 1 2 ψcd c d ψ ab (a H b) + Q ab (ψ, ψ). Add constraint damping terms for stability: 0 = R ab (a C b) + γ 0 [ t(a C b) 1 2 ψ ab t c C c ], where t a is a unit timelike vector field. Since C a = H a + Γ a depends only on first derivatives of the metric, these additional terms do not change the hyperbolic structure of the system. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 13 / 17
36 Numerical Tests of the GH Evolution System 3D numerical evolutions of static black-hole spacetimes illustrate the constraint damping properties of the GH evolution system. These evolutions are stable and convergent when γ 0 = 1. The boundary conditions used for this simple test problem freeze the incoming characteristic fields to their initial analytical values. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 14 / 17
37 What Do We Mean By Hyperbolic? We have argued that Einstein s equation is manifestly hyperbolic because its principal part is the same as the scalar wave equation. Exactly what does this mean? Does this make sense? Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 15 / 17
38 What Do We Mean By Hyperbolic? We have argued that Einstein s equation is manifestly hyperbolic because its principal part is the same as the scalar wave equation. Exactly what does this mean? Does this make sense? From a pragmatic physicist s point of view, hyperbolic means anything that acts like the wave equation, i.e. any system of equations having a well posed initial-boundary value problem. Symmetric hyperbolic systems are one class of equations for which suitable well-posedness theorems exist, and which are general enough to include Einstein s equations together with most of the other dynamical field equations used by physicists. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 15 / 17
39 What Do We Mean By Hyperbolic? We have argued that Einstein s equation is manifestly hyperbolic because its principal part is the same as the scalar wave equation. Exactly what does this mean? Does this make sense? From a pragmatic physicist s point of view, hyperbolic means anything that acts like the wave equation, i.e. any system of equations having a well posed initial-boundary value problem. Symmetric hyperbolic systems are one class of equations for which suitable well-posedness theorems exist, and which are general enough to include Einstein s equations together with most of the other dynamical field equations used by physicists. Evolution equations of the form, t u α + A k α β(u, x, t) k u β = F α (u, x, t), for a collection of dynamical fields u α, are called symmetric hyperbolic if there exists a positive definite S αβ having the property that S αγ A k γ β A k αβ = Ak βα. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 15 / 17
40 Example: Scalar Wave Equation Consider the scalar wave equation in flat space, expressed in terms of arbitrary spatial coordinates: 0 = 2 t ψ + k k ψ = 2 t ψ + gkl ( k l ψ Γ n kl nψ). Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 16 / 17
41 Example: Scalar Wave Equation Consider the scalar wave equation in flat space, expressed in terms of arbitrary spatial coordinates: 0 = 2 t ψ + k k ψ = 2 t ψ + gkl ( k l ψ Γ n kl nψ). Define the first-order dynamical fields, u α = {ψ, Π, Φ k }, which satisfy the following evolution equations: t ψ = Π, t Π + k Φ k = 0, t Φ k + k Π = 0. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 16 / 17
42 Example: Scalar Wave Equation Consider the scalar wave equation in flat space, expressed in terms of arbitrary spatial coordinates: 0 = 2 t ψ + k k ψ = 2 t ψ + gkl ( k l ψ Γ n kl nψ). Define the first-order dynamical fields, u α = {ψ, Π, Φ k }, which satisfy the following evolution equations: t ψ = Π, t Π + k Φ k = 0, t Φ k + k Π = 0. The principal part of this system, t u α + A k α β k u β 0, is given: t ψ Π Φ x Φ y Φ z g xx g xy g xz x ψ Π Φ x Φ y Φ z Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 16 / 17
43 Example: Scalar Wave Equation II The symmetrizer for the first-order scalar field system is: ds 2 = S αβ du α du β = Λ 2 dψ 2 + dπ 2 + g mn dφ m dφ n. Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 17 / 17
44 Example: Scalar Wave Equation II The symmetrizer for the first-order scalar field system is: ds 2 = S αβ du α du β = Λ 2 dψ 2 + dπ 2 + g mn dφ m dφ n. Check the symmetrization of the characteristic matrices: Λ S αγ A x γ g xx g xy g xz β = 0 0 g xx g xy g xz g yx g yy g yz g zx g zy g zz = g xx g xy g xz 0 g xx g yx g zx = Ax αβ Lee Lindblom (CASS UCSD) Solving Einstein s Equation Numerically 2014/11/21 MSC Tsinghua U 17 / 17
Solving Einstein s Equations: PDE Issues
Solving Einstein s Equations: PDE Issues Lee Lindblom Theoretical Astrophysics, Caltech Mathematical and Numerical General Relativity Seminar University of California at San Diego 22 September 2011 Lee
More informationGeneralized Harmonic Evolutions of Binary Black Hole Spacetimes
Generalized Harmonic Evolutions of Binary Black Hole Spacetimes Lee Lindblom California Institute of Technology AMS Meeting :: New Orleans :: 7 January 2007 Lee Lindblom (Caltech) Generalized Harmonic
More informationCalculating Accurate Waveforms for LIGO and LISA Data Analysis
Calculating Accurate Waveforms for LIGO and LISA Data Analysis Lee Lindblom Theoretical Astrophysics, Caltech HEPL-KIPAC Seminar, Stanford 17 November 2009 Results from the Caltech/Cornell Numerical Relativity
More informationSolving Einstein s Equations for Binary Black Hole Spacetimes
Solving Einstein s Equations for Binary Black Hole Spacetimes Lee Lindblom Theoretical Astrophysics, Caltech University of Wisconsin at Milwaukee Department of Physics Colloquium 14 October 2011 Lee Lindblom
More informationIntroduction to Numerical Relativity I. Erik Schnetter, Pohang, July 2007
Introduction to Numerical Relativity I Erik Schnetter, Pohang, July 2007 Lectures Overview I. The Einstein Equations (Formulations and Gauge Conditions) II. Analysis Methods (Horizons and Gravitational
More informationGeneralized Harmonic Gauge Drivers
Generalized Harmonic Gauge Drivers Lee Lindblom Caltech Collaborators: Keith Matthews, Oliver Rinne, and Mark Scheel GR18 Sydney, Australia 13 July 2007 Gauge conditions are specified in the GH Einstein
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 informationNumerical Relativity
Numerical Relativity Mark A. Scheel Walter Burke Institute for Theoretical Physics Caltech July 7, 2015 Mark A. Scheel Numerical Relativity July 7, 2015 1 / 34 Outline: Motivation 3+1 split and ADM equations
More informationGeneralized Harmonic Coordinates Using Abigel
Outline Introduction The Abigel Code Jennifer Seiler jennifer.seiler@aei.mpg.de Max-Planck-Institut für Gravitationsphysik (Albert-Einstein-Institut) Oberjoch, Germany 10th June 2006 Outline Introduction
More informationNumerical Simulations of Black Hole Spacetimes
Numerical Simulations of Black Hole Spacetimes Lee Lindblom Senior Research Associate Theoretical Astrophysics Physics Research Conference California Institute of Technology 24 May 2007 Lee Lindblom (Caltech)
More informationGeneral Relativity and Gravitational Waveforms
General Relativity and Gravitational Waveforms Deirdre Shoemaker Center for Relativistic Astrophysics School of Physics Georgia Institute of Technology Kavli Summer Program in Astrophysics 2017 Astrophysics
More informationGeneral Relativity ASTR 2110 Sarazin. Einstein s Equation
General Relativity ASTR 2110 Sarazin Einstein s Equation Curvature of Spacetime 1. Principle of Equvalence: gravity acceleration locally 2. Acceleration curved path in spacetime In gravitational field,
More informationTHE INITIAL VALUE FORMULATION OF GENERAL RELATIVITY
THE INITIAL VALUE FORMULATION OF GENERAL RELATIVITY SAM KAUFMAN Abstract. The (Cauchy) initial value formulation of General Relativity is developed, and the maximal vacuum Cauchy development theorem is
More informationImproved constrained scheme for the Einstein equations: An approach to the uniqueness issue
Improved constrained scheme for the Einstein equations: An approach to the uniqueness issue Jérôme Novak (Jerome.Novak@obspm.fr) Laboratoire Univers et Théories (LUTH) CNRS / Observatoire de Paris / Université
More informationKadath: a spectral solver and its application to black hole spacetimes
Kadath: a spectral solver and its application to black hole spacetimes Philippe Grandclément Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris F-92195 Meudon, France philippe.grandclement@obspm.fr
More informationPAPER 52 GENERAL RELATIVITY
MATHEMATICAL TRIPOS Part III Monday, 1 June, 2015 9:00 am to 12:00 pm PAPER 52 GENERAL RELATIVITY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight.
More informationOptimal Constraint Projection for Hyperbolic Evolution Systems
Optimal Constraint Projection for Hyperbolic Evolution Systems Michael Holst 1,2, Lee Lindblom 1, Robert Owen 1, Harald P. Pfeiffer 1, Mark A. Scheel 1, and Lawrence E. Kidder 3 1 Theoretical Astrophysics
More informationA non-strictly hyperbolic system for the Einstein equations with arbitrary lapse and shift
IFP-UNC-518 TAR-UNC-054 gr-qc/9607006 A non-strictly hyperbolic system for the Einstein equations with arbitrary lapse and shift Andrew Abrahams, Arlen Anderson, Yvonne Choquet-Bruhat[*] and James W. York,
More informationFirst order BSSN formulation of Einstein s field equations
David Brown 1 Peter Diener 2 3 Jan Hesthaven 4 Frank Herrmann 3 Abdul Mroué 5 Olivier Sarbach 6 Erik Schnetter 7 Manuel Tiglio 3 Michael Wagman 4 1 North Carolina State University 2 Louisiana State University
More informationInitial-Value Problems in General Relativity
Initial-Value Problems in General Relativity Michael Horbatsch March 30, 2006 1 Introduction In this paper the initial-value formulation of general relativity is reviewed. In section (2) domains of dependence,
More informationSolving PDEs Numerically on Manifolds with Arbitrary Spatial Topologies
Solving PDEs Numerically on Manifolds with Arbitrary Spatial Topologies Lee Lindblom Theoretical Astrophysics, Caltech Center for Astrophysics and Space Sciences, UC San Diego. Collaborators: Béla Szilágyi,
More informationCurved Spacetime III Einstein's field equations
Curved Spacetime III Einstein's field equations Dr. Naylor Note that in this lecture we will work in SI units: namely c 1 Last Week s class: Curved spacetime II Riemann curvature tensor: This is a tensor
More informationInitial Data for Black-Hole Binaries
Initial Data for Black-Hole Binaries Gregory B. Cook Wake Forest University June 11/1, 004 Abstract We will examine the current state of our efforts to generate astrophysically realistic initial data for
More informationGeneral Relativity (225A) Fall 2013 Assignment 8 Solutions
University of California at San Diego Department of Physics Prof. John McGreevy General Relativity (5A) Fall 013 Assignment 8 Solutions Posted November 13, 013 Due Monday, December, 013 In the first two
More informationStability and Instability of Black Holes
Stability and Instability of Black Holes Stefanos Aretakis September 24, 2013 General relativity is a successful theory of gravitation. Objects of study: (4-dimensional) Lorentzian manifolds (M, g) which
More informationWell-Posed Initial-Boundary Value Problem for Constrained Evolution
Well-Posed Initial-Boundary Value Problem for Constrained Evolution Alexander Alekseenko California State University, Northridge Design of parallel algorithms for Einstein s equations requires a better
More information8 Symmetries and the Hamiltonian
8 Symmetries and the Hamiltonian Throughout the discussion of black hole thermodynamics, we have always assumed energy = M. Now we will introduce the Hamiltonian formulation of GR and show how to define
More informationGetting The Spin Right In Black-Hole Binaries
Getting The Spin Right In Black-Hole Binaries Gregory B. Cook Wake Forest University July 25, 2005 Abstract We will take a detailed look at the issues involved in setting the spin of a black hole during
More informationUNIVERSITY OF DUBLIN
UNIVERSITY OF DUBLIN TRINITY COLLEGE JS & SS Mathematics SS Theoretical Physics SS TSM Mathematics Faculty of Engineering, Mathematics and Science school of mathematics Trinity Term 2015 Module MA3429
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 informationPAPER 311 BLACK HOLES
MATHEMATICAL TRIPOS Part III Friday, 8 June, 018 9:00 am to 1:00 pm PAPER 311 BLACK HOLES Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY
More informationInflation in Flatland
Inflation in Flatland Austin Joyce Center for Theoretical Physics Columbia University Kurt Hinterbichler, AJ, Justin Khoury, 1609.09497 Theoretical advances in particle cosmology, University of Chicago,
More informationLecture VIII: Linearized gravity
Lecture VIII: Linearized gravity Christopher M. Hirata Caltech M/C 350-17, Pasadena CA 91125, USA (Dated: November 5, 2012) I. OVERVIEW We are now ready to consider the solutions of GR for the case of
More informationCurved Spacetime I. Dr. Naylor
Curved Spacetime I Dr. Naylor Last Week Einstein's principle of equivalence We discussed how in the frame of reference of a freely falling object we can construct a locally inertial frame (LIF) Space tells
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 informationCausality, hyperbolicity, and shock formation in Lovelock theories
Causality, hyperbolicity, and shock formation in Lovelock theories Harvey Reall DAMTP, Cambridge University HSR, N. Tanahashi and B. Way, arxiv:1406.3379, 1409.3874 G. Papallo, HSR arxiv:1508.05303 Lovelock
More informationAn introduction to General Relativity and the positive mass theorem
An introduction to General Relativity and the positive mass theorem National Center for Theoretical Sciences, Mathematics Division March 2 nd, 2007 Wen-ling Huang Department of Mathematics University of
More informationGeneral Relativity I
General Relativity I presented by John T. Whelan The University of Texas at Brownsville whelan@phys.utb.edu LIGO Livingston SURF Lecture 2002 July 5 General Relativity Lectures I. Today (JTW): Special
More informationRigidity of Black Holes
Rigidity of Black Holes Sergiu Klainerman Princeton University February 24, 2011 Rigidity of Black Holes PREAMBLES I, II PREAMBLE I General setting Assume S B two different connected, open, domains and
More informationarxiv:gr-qc/ v1 3 Feb 1998
A New Method for Solving the Initial Value Problem with Application to Multiple Black Holes by John Baker and Raymond Puzio arxiv:gr-qc/980006v 3 Feb 998.Introduction The approaching date for broad band
More informationOptimal constraint projection for hyperbolic evolution systems
PHYSICAL REVIEW D, VOLUME 70, 084017 Optimal constraint projection for hyperbolic evolution systems Michael Holst, 1,2 Lee Lindblom, 1 Robert Owen, 1 Harald P. Pfeiffer, 1 Mark A. Scheel, 1 and Lawrence
More informationRIGIDITY OF STATIONARY BLACK HOLES WITH SMALL ANGULAR MOMENTUM ON THE HORIZON
RIGIDITY OF STATIONARY BLACK HOLES WITH SMALL ANGULAR MOMENTUM ON THE HORIZON S. ALEXAKIS, A. D. IONESCU, AND S. KLAINERMAN Abstract. We prove a black hole rigidity result for slowly rotating stationary
More informationSolving Einstein s equations with dual coordinate frames
PHYSICAL REVIEW D 74, 104006 (2006) Solving Einstein s equations with dual coordinate frames Mark A. Scheel, 1 Harald P. Pfeiffer, 1 Lee Lindblom, 1 Lawrence E. Kidder, 2 Oliver Rinne, 1 and Saul A. Teukolsky
More informationA Simple, Direct Finite Differencing of the Einstein Equations
A Simple, Direct Finite Differencing of the Einstein Equations Travis M. Garrett Louisiana State University (Dated: 2008.11.12) We investigate a simple variation of the Generalized Harmonic method for
More informationUmbilic cylinders in General Relativity or the very weird path of trapped photons
Umbilic cylinders in General Relativity or the very weird path of trapped photons Carla Cederbaum Universität Tübingen European Women in Mathematics @ Schloss Rauischholzhausen 2015 Carla Cederbaum (Tübingen)
More informationUniformity of the Universe
Outline Universe is homogenous and isotropic Spacetime metrics Friedmann-Walker-Robertson metric Number of numbers needed to specify a physical quantity. Energy-momentum tensor Energy-momentum tensor of
More informationTime-Periodic Solutions of the Einstein s Field Equations II: Geometric Singularities
Science in China Series A: Mathematics 2009 Vol. 52 No. 11 1 16 www.scichina.com www.springerlink.com Time-Periodic Solutions of the Einstein s Field Equations II: Geometric Singularities KONG DeXing 1,
More informationarxiv: v2 [gr-qc] 24 Apr 2012
Boundary Conditions for the Gravitational Field Jeffrey Winicour 1,2 1 Department of Physics and Astronomy University of Pittsburgh, Pittsburgh, PA 1526, USA 2 Max-Planck-Institut für Gravitationsphysik,
More informationNonlinear wave-wave interactions involving gravitational waves
Nonlinear wave-wave interactions involving gravitational waves ANDREAS KÄLLBERG Department of Physics, Umeå University, Umeå, Sweden Thessaloniki, 30/8-5/9 2004 p. 1/38 Outline Orthonormal frames. Thessaloniki,
More informationExact Solutions of the Einstein Equations
Notes from phz 6607, Special and General Relativity University of Florida, Fall 2004, Detweiler Exact Solutions of the Einstein Equations These notes are not a substitute in any manner for class lectures.
More informationSpecial Theory of Relativity
June 17, 2008 1 1 J.D.Jackson, Classical Electrodynamics, 3rd Edition, Chapter 11 Introduction Einstein s theory of special relativity is based on the assumption (which might be a deep-rooted superstition
More informationEinstein Toolkit Workshop. Joshua Faber Apr
Einstein Toolkit Workshop Joshua Faber Apr 05 2012 Outline Space, time, and special relativity The metric tensor and geometry Curvature Geodesics Einstein s equations The Stress-energy tensor 3+1 formalisms
More informationA Description of the Initial Value Formulation of. Mark Miller, Syracuse University. October 10, 1994
A Description of the Initial Value Formulation of Vacuum General Relativity for the Non-Specialist 1 Mark Miller, Syracuse University October 10, 1994 1 Work supported by NSF ASC 93 18152/PHY 93 18152
More informationNewman-Penrose formalism in higher dimensions
Newman-Penrose formalism in higher dimensions V. Pravda various parts in collaboration with: A. Coley, R. Milson, M. Ortaggio and A. Pravdová Introduction - algebraic classification in four dimensions
More informationLecture 9: RR-sector and D-branes
Lecture 9: RR-sector and D-branes José D. Edelstein University of Santiago de Compostela STRING THEORY Santiago de Compostela, March 6, 2013 José D. Edelstein (USC) Lecture 9: RR-sector and D-branes 6-mar-2013
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Preliminary Math Concept of Stress Stress Components Equilibrium
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 informationA DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any
Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»
More informationReview of General Relativity
Lecture 3 Review of General Relativity Jolien Creighton University of Wisconsin Milwaukee July 16, 2012 Whirlwind review of differential geometry Coordinates and distances Vectors and connections Lie derivative
More informationBoundary conditions for the Einstein evolution system
PHYSICAL REVIEW D 71, 064020 (2005) Boundary conditions for the Einstein evolution system Lawrence E. Kidder, 1 Lee Lindblom, 2 Mark A. Scheel, 2 Luisa T. Buchman, 3 and Harald P. Pfeiffer 2 1 Center for
More informationLevel sets of the lapse function in static GR
Level sets of the lapse function in static GR Carla Cederbaum Mathematisches Institut Universität Tübingen Auf der Morgenstelle 10 72076 Tübingen, Germany September 4, 2014 Abstract We present a novel
More informationIntroduction to Numerical Relativity
Introduction to Numerical Relativity Given a manifold M describing a spacetime with 4-metric we want to foliate it via space-like, threedimensional hypersurfaces: {Σ i }. We label such hypersurfaces with
More informationGG303 Lecture 6 8/27/09 1 SCALARS, VECTORS, AND TENSORS
GG303 Lecture 6 8/27/09 1 SCALARS, VECTORS, AND TENSORS I Main Topics A Why deal with tensors? B Order of scalars, vectors, and tensors C Linear transformation of scalars and vectors (and tensors) II Why
More informationBlack-Hole Binary Initial Data: Getting the Spin Right
Black-Hole Binary Initial Data: Getting the Spin Right Gregory B. Cook Wake Forest University November 4, 2005 Abstract Using the conformal thin-sandwich approach for constructing initial data together
More informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationEXCISION TECHNIQUE IN CONSTRAINED FORMULATIONS OF EINSTEIN EQUATIONS
EXCISION TECHNIQUE IN CONSTRAINED FORMULATIONS OF EINSTEIN EQUATIONS Journée Gravitation et Physique Fondamentale Meudon, 27 May 2014 Isabel Cordero-Carrión Laboratoire Univers et Théories (LUTh), Observatory
More informationGravity as Machian Shape Dynamics
Gravity as Machian Shape Dynamics Julian Barbour Alternative Title: Was Spacetime a Glorious Historical Accident? Time will not be fused with space but emerge from the timeless shape dynamics of space.
More informationM E 320 Professor John M. Cimbala Lecture 10
M E 320 Professor John M. Cimbala Lecture 10 Today, we will: Finish our example problem rates of motion and deformation of fluid particles Discuss the Reynolds Transport Theorem (RTT) Show how the RTT
More informationDynamics of star clusters containing stellar mass black holes: 1. Introduction to Gravitational Waves
Dynamics of star clusters containing stellar mass black holes: 1. Introduction to Gravitational Waves July 25, 2017 Bonn Seoul National University Outline What are the gravitational waves? Generation of
More informationLagrange Multipliers
Optimization with Constraints As long as algebra and geometry have been separated, their progress have been slow and their uses limited; but when these two sciences have been united, they have lent each
More informationKerr black hole and rotating wormhole
Kerr Fest (Christchurch, August 26-28, 2004) Kerr black hole and rotating wormhole Sung-Won Kim(Ewha Womans Univ.) August 27, 2004 INTRODUCTION STATIC WORMHOLE ROTATING WORMHOLE KERR METRIC SUMMARY AND
More informationBlack Hole Formation in Randall-Sundrum II Braneworlds
Black Hole Formation in Randall-Sundrum II Braneworlds Matt Choptuik CIFAR Cosmology & Gravity Program Dept of Physics & Astronomy, UBC Workshop on Numerical and Mathematical Relativity Oppurg, DE September
More informationarxiv:gr-qc/ v2 28 Nov 2001
Boundary conditions in linearized harmonic gravity Béla Szilágyi 1, Bernd Schmidt 2, Jeffrey Winicour 1,2 1 Department of Physics and Astronomy University of Pittsburgh, Pittsburgh, PA 15260, USA 2 Max-Planck-Institut
More informationTheory. V H Satheeshkumar. XXVII Texas Symposium, Dallas, TX December 8 13, 2013
Department of Physics Baylor University Waco, TX 76798-7316, based on my paper with J Greenwald, J Lenells and A Wang Phys. Rev. D 88 (2013) 024044 with XXVII Texas Symposium, Dallas, TX December 8 13,
More informationGroup, Rings, and Fields Rahul Pandharipande. I. Sets Let S be a set. The Cartesian product S S is the set of ordered pairs of elements of S,
Group, Rings, and Fields Rahul Pandharipande I. Sets Let S be a set. The Cartesian product S S is the set of ordered pairs of elements of S, A binary operation φ is a function, S S = {(x, y) x, y S}. φ
More informationIntroduction to Numerical Relativity I Intro & Geometry
Ë ¼ Å ¼ Â ¼ Introduction to Numerical Relativity I Intro & Geometry S. Husa University of the Balearic Islands, sascha.husa@uib.es August 5, 2010 Ë ¾ Å ¾ Â ¾ Ë ½ Å ½ Â S. Husa (UIB) Numerical Relativity
More informationON THE FOLIATION OF SPACE-TIME BY CONSTANT MEAN CURVATURE HYPERSURFACES
ON THE FOLIATION OF SPACE-TIME BY CONSTANT MEAN CURVATURE HYPERSURFACES CLAUS GERHARDT Abstract. We prove that the mean curvature τ of the slices given by a constant mean curvature foliation can be used
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 informationRotational & Rigid-Body Mechanics. Lectures 3+4
Rotational & Rigid-Body Mechanics Lectures 3+4 Rotational Motion So far: point objects moving through a trajectory. Next: moving actual dimensional objects and rotating them. 2 Circular Motion - Definitions
More informationBOUNDARY CONDITIONS FOR THE EINSTEIN-CHRISTOFFEL FORMULATION OF EINSTEIN S EQUATIONS
Sixth Mississippi State Conference on Differential Equations and Computational Simulations, Electronic Journal of Differential Equations, Conference 15 (2007), pp. 11 27. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
More informationConnection Variables in General Relativity
Connection Variables in General Relativity Mauricio Bustamante Londoño Instituto de Matemáticas UNAM Morelia 28/06/2008 Mauricio Bustamante Londoño (UNAM) Connection Variables in General Relativity 28/06/2008
More informationHyperbolic Geometric Flow
Hyperbolic Geometric Flow Kefeng Liu Zhejiang University UCLA Page 1 of 41 Outline Introduction Hyperbolic geometric flow Local existence and nonlinear stability Wave character of metrics and curvatures
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 informationProblem Set 2 Due Tuesday, September 27, ; p : 0. (b) Construct a representation using five d orbitals that sit on the origin as a basis: 1
Problem Set 2 Due Tuesday, September 27, 211 Problems from Carter: Chapter 2: 2a-d,g,h,j 2.6, 2.9; Chapter 3: 1a-d,f,g 3.3, 3.6, 3.7 Additional problems: (1) Consider the D 4 point group and use a coordinate
More informationLecture 13 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell
Lecture 13 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell 1. Covariant Geometry - We would like to develop a mathematical framework
More informationGeneral Relativity (225A) Fall 2013 Assignment 2 Solutions
University of California at San Diego Department of Physics Prof. John McGreevy General Relativity 5A) Fall 13 Assignment Solutions Posted October 3, 13 Due Monday, October 15, 13 1. Special relativity
More informationCausal RG equation for Quantum Einstein Gravity
Causal RG equation for Quantum Einstein Gravity Stefan Rechenberger Uni Mainz 14.03.2011 arxiv:1102.5012v1 [hep-th] with Elisa Manrique and Frank Saueressig Stefan Rechenberger (Uni Mainz) Causal RGE for
More informationAn Overview of Mathematical General Relativity
An Overview of Mathematical General Relativity José Natário (Instituto Superior Técnico) Geometria em Lisboa, 8 March 2005 Outline Lorentzian manifolds Einstein s equation The Schwarzschild solution Initial
More informationStatus of Hořava Gravity
Status of Institut d Astrophysique de Paris based on DV & T. P. Sotiriou, PRD 85, 064003 (2012) [arxiv:1112.3385 [hep-th]] DV & T. P. Sotiriou, JPCS 453, 012022 (2013) [arxiv:1212.4402 [hep-th]] DV, arxiv:1502.06607
More informationCosmology on Simplicial Complexes
Gravitation and Regge Calculus Astro Coee, Frankfurt, April 2015 Outline Gravitation and Regge Calculus 1 Gravitation and Regge Calculus Foundations of General Relativity Geometric Structure of Regge Calculus
More informationLecture XIV: Global structure, acceleration, and the initial singularity
Lecture XIV: Global structure, acceleration, and the initial singularity Christopher M. Hirata Caltech M/C 350-17, Pasadena CA 91125, USA (Dated: December 5, 2012) I. OVERVIEW In this lecture, we will
More informationAttempts at relativistic QM
Attempts at relativistic QM based on S-1 A proper description of particle physics should incorporate both quantum mechanics and special relativity. However historically combining quantum mechanics and
More informationClassification theorem for the static and asymptotically flat Einstein-Maxwell-dilaton spacetimes possessing a photon sphere
Classification theorem for the static and asymptotically flat Einstein-Maxwell-dilaton spacetimes possessing a photon sphere Boian Lazov and Stoytcho Yazadjiev Varna, 2017 Outline 1 Motivation 2 Preliminaries
More informationThe Klein-Gordon Equation Meets the Cauchy Horizon
Enrico Fermi Institute and Department of Physics University of Chicago University of Mississippi May 10, 2005 Relativistic Wave Equations At the present time, our best theory for describing nature is Quantum
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 informationHorava-Lifshitz. Based on: work with ( ), ( ) arxiv: , JCAP 0911:015 (2009) arxiv:
@ 2010 2 18 Horava-Lifshitz Based on: work with ( ), ( ) arxiv:0908.1005, JCAP 0911:015 (2009) arxiv:1002.3101 Motivation A quantum gravity candidate Recently Horava proposed a power-counting renormalizable
More informationTheoretical Concepts of Spin-Orbit Splitting
Chapter 9 Theoretical Concepts of Spin-Orbit Splitting 9.1 Free-electron model In order to understand the basic origin of spin-orbit coupling at the surface of a crystal, it is a natural starting point
More informationSPECIAL RELATIVITY AND ELECTROMAGNETISM
SPECIAL RELATIVITY AND ELECTROMAGNETISM MATH 460, SECTION 500 The following problems (composed by Professor P.B. Yasskin) will lead you through the construction of the theory of electromagnetism in special
More informationEinstein boundary conditions in relation to constraint propagation for the initial-boundary value problem of the Einstein equations
PHYSICAL REVIEW D 69, 24020 2004 Einstein boundary conditions in relation to constraint propagation for the initial-boundary value problem of the Einstein equations Simonetta Frittelli* Department of Physics,
More informationCourse 2BA1: Hilary Term 2007 Section 8: Quaternions and Rotations
Course BA1: Hilary Term 007 Section 8: Quaternions and Rotations David R. Wilkins Copyright c David R. Wilkins 005 Contents 8 Quaternions and Rotations 1 8.1 Quaternions............................ 1 8.
More information