Non-Equilibrium Aspects of Relic Neutrinos: From Freeze-out to the Present Day
|
|
- Rudolph Heath
- 5 years ago
- Views:
Transcription
1 Non-Equilibrium Aspects of Relic Neutrinos: From Freeze-out to the Present Day Program in Applied Mathematics The Dissertation Defense Presentation Advisor J. Rafelski, Department of Physics, March 20, 2014
2 Outline Part 1: Eras of Cosmology Relic Neutrinos Overview of Cosmology Part 2: Study of Neutrino Distribution using Conservation Laws Chemical vs Kinetic Equilibrium Neutrino Freeze-out Neutrinos Today Part 3: General Relativistic Boltzmann Equation Boltzmann Equation and Scattering Operators Neutrino Freeze-out Revisited
3 Part 1: Relic Neutrino Background At a temperature of 5 MeV the Universe consisted of e ±, photons, and neutrinos. At around 1 MeV neutrinos stop interacting or freeze-out and free stream through the universe. Today they comprise the relic neutrino background. Photons freeze-out around 0.25 ev and today they make up the Cosmic Microwave Background (CMB), currently at T γ = mev. Direct measurement: Relic neutrinos have not been directly measured. Our results give some input into such efforts. Indirect measurement: Impact on speed of Universe expansion constrains neutrino mass and number of massless particles. Theory: This motivates interest in modeling energy flow from e ± annihilation to neutrinos during reheating.
4 Overview of Standard Cosmology Planck satellite measures CMB which froze out at 0.25 ev. We need to model the Universe evolution between these two eras. The Friedmann Lemaitre Robertson Walker (FRW) Universe assumes the properties Homogeneous Isotropic [ ] dr ds 2 = dt 2 a 2 2 (t) 1 kr 2 + r2 (dθ 2 + sin 2 (θ)dφ 2 ). (1) Standard ΛCDM model is spatially flat (k = 0) [1]. a(t) determines the distance between objects at rest in the Universe frame (comoving). [1] Planck Collaboration, arxiv:
5 Einstein s Equations Einstein s equations for a perfect fluid source ( ) R G µν = R µν 2 + Λ g µν = 8πG N T µν, T ν µ = diag(ρ, P, P, P), ρ ρ + Λ/8πG N, P P Λ/8πG N (2) In a spatially flat FRW universe these reduce to Hubble equation H 2 = ρ 3M 2 p, H ȧ a, M p 1/ 8πG N = GeV Conservation of stress energy tensor (3) ρ = 3(ρ + P)H. (4) Together with an equation of state, P(ρ), one can solve for the large scale dynamics of the Universe.
6 Reheating In the standard model of neutrino freeze-out, the relic neutrino background and CMB temperatures differ by a reheating factor R ν T ν /T γ = ( ) 4 1/3. (5) 11 This is the result of energy and entropy from e ± annihilation going into photons. In general, the reheating ratio connects the temperature redshift to the geometric redshift R 1 + z T/T now, 1 + z a now a(t). (6)
7 Reheating and Particle Disappearance History Conservation of entropy relates R at different times ( R1 R 2 ) 3 = i g i(t 1,i /T 1,γ ) 3 j g j(t 2,j /T 2,γ ) 3 (7) Reheating periods occurred regularly throughout the history of the Universe, whenever a particle species disappeared as the temperature dropped below the mass.
8 Reheating and Particle Disappearance History 3 (1 + z)tnow/t QGP phase transition, disappearance of resonances and muons Disappearance of bottom, tau, and charm Disappearance of top, Higgs, Z, and W 1.5 e ± annihilation, photon reheating, and neutrino freeze-out Time T [MeV] Figure: Disappearance of degrees of freedom. The Universe volume inflated by a factor 27 above the thermal red shift scale as massive particles disappeared successively from the inventory. [1] J. Rafelski and J. Birrell, Submitted to: proceedings of Strangeness in Quark Matter Conference July 2013, Birmingham, UK
9 Deceleration Parameter Definition: q äa ȧ 2. (8) Using Einsteins equations we can derive an expression in terms of matter content q = 1 ( P ). (9) 2 ρ
10 Deceleration Parameter In particular: Radiation dominated universe: P = ρ/3 = q = 1. Matter dominated universe: P ρ = q = 1/2. Dark energy (Λ) dominated universe: P = ρ = q = 1.
11 T [ev] q 10 3 Radiation dominated Matter dominated T γ Dark Energy dominated T ν q t recomb t reion t [yr] Figure: From the present day until near BBN.
12 T [MeV] q T γ T ν 0.96 q t k t i BBN t f BBN Figure: From the end of baryon antimatter annihilation through BBN. t [s]
13 Outline Part 1: Background Relic Neutrinos Overview of Cosmology Part 2: Study of Neutrino Distribution using Conservation Laws Chemical vs Kinetic Equilibrium Neutrino Freeze-out Neutrinos today Part 3: General Relativistic Boltzmann Equation Boltzmann Equation and Scattering Operators Neutrino Freeze-out Revisited
14 Part 2: Why Study Neutrino Freeze-out in Detail? Energy in neutrinos impacts speed of Universe expansion. N ν is defined by comparing the total neutrino energy density to the energy density of a massless fermion with two degrees of freedom and standard reheating ratio R ν N ν = ρ ν π2 (R ν T γ ) 4, R ν = ( ) 4 1/3. (10) 11 Planck satellite data gives N ν = 3.36 ± 0.34 (CMB no priors) and N ν = 3.62 ± 0.25 (CMB + H 0 ) [1]. The currently used theoretical value is N ν = [2]. [1] Planck Collaboration, arxiv: [2] G. Mangano et. al., Nucl. Phys. B 729, 221 (2005)
15 Why Study Neutrino Freeze-out in Detail? The tension between the Planck result and reheating study motivates our work. This tension has inspired exotic theories introducing new particles [1]. Most conservative explanation is that some extra energy went into neutrinos during e ± annihilation i.e. neutrino reheating. This is the direction we explore. [1] Steven Weinberg Phys. Rev. Lett. 110, (2013)
16 Boltzmann equation Particle dynamics given by Boltzmann equation p α x αf Γ j µνp µ p ν p jf = C[f ]. (11) Here Γ α µν is the Levi-Civita connection (Christoffel symbols), the distribution function f is a function on the mass shell g αβ p α p β = m 2, (12) and C[f ] is the collision operator and encodes all information about point interactions between particles. We will return to solution of the Boltzmann equation later, but for now we use a simpler matter model that is independent of the details of C[f ].
17 Kinetic Equilibrium Model First we perform a model independent study of neutrino freeze-out, assuming instantaneous chemical/kinetic equilibrium and sharp freeze-out transitions. Dynamics will be derived from conservation laws involving the moments T µν = g ν 8π 3 n ν = g ν 8π 3 s µ = g ν 8π 3 f pµ p ν p 0 gd 3 p, (13) f pν gd 3 p, (14) p 0 h pµ gd 3 p, (15) p 0 h = f ln(f ) + (1 f ) ln(1 f ).
18 Combined Chemical and Kinetic Equilibrium Both number changing processes (chemical) and momentum exchange (kinetic) allowed. Distribution function obtained by maximizing entropy subject to fixed energy. f ch = g ν 8π 3 1 e E/T + 1, T ch < T. (16) T ch is the chemical freeze-out temperature, below which number changing processes cease.
19 Kinetic (but not chemical) Equilibrium Particles exchange momentum but particle (and antiparticle) numbers fixed. Distribution function obtained by maximizing entropy, subject to fixed particle (and antiparticle) number. f k = g ν 8π 3 1 Υ 1 e E/T + 1 T k < T < T ch. (17) The fugacity, Υ, controls the phase space occupancy. This was first introduced in [1] in the context of QGP. T k is the kinetic freeze-out temperature. Below this, there are no (significant) interactions and the particles free stream. [1]J. Rafelski and B. Muller, Phys. Rev. Lett. 48,1066 (1982)
20 Relic Neutrino Fugacity is NOT a Chemical Potential Υ is the same for particles and antiparticles (we ignore the small matter-antimatter asymmetry which would yield Υ Υe ±µ/t ) Υ(t) is time dependent, not just an initial condition (initially Υ = 1). We find that Υ < 1 emerges dynamically from chemical equilibrium. It is not put in by hand.
21 Dynamics with Kinetic Equilibrium Ansatz Divergence freedom of the total stress-energy tensor implies d dt (a3 s) + i σ i d dt (a3 n i ) = 0, Υ i e σ i. (18) where s is entropy density and n i the particle number densities. Assuming instantaneous chemical/kinetic equilibrium and sharp freeze-out transitions, this, along with the Einstein equations, fixes the dynamics of a(t), T(t), Υ i (t). H 2 = ρ 3M 2 p, d dt (a3 s) = 0, d dt (a3 n i ) = 0 ( when T < T ch ). (19)
22 The Dynamical Emergence of Fugacity Distinct chemical and thermal freeze-out temperatures lead to Υ ν < 1 for T < T ch. This fact has been overlooked in the literature. The mechanism is as follows: If all particles are (effectively) massless there is no effect T 1/a, Υ = constant. (20) When a mass scale, such as m e, becomes important and e + e annihilation begins, conservation of entropy gives d/dt(at) > 0 (reheating). For T < T ch, preserving the comoving number of neutrinos when at is increasing implies Υ < 0. a 3 n ν a 3 T 3 1 Υ 1 ν e u + 1 u2 du (21)
23 Neutrino Freeze-out 1. Conservation of comoving neutrino number: T1 3V 1 Tk 3V k = 2 3ζ(3) 0 u 2 du Υ 1 ν (T k )e u + 1. (22) 2. Conservation of e ±, γ, neutrino entropy before neutrino freeze-out: ( 7 8 g ν + 7 ) 2π 2 8 g e ± + g γ 45 T3 1 V 1 = (23) ( ) 2π 2 45 g γtk 3 + S e ±(T k) + S ν (T k ) V k. 3. Conservation of e ±, γ entropy after neutrino freeze-out: 2π 2 ( ) 2π 45 g γtγ,2 3 2 V 2 = 45 g γtk 3 + S e ±(T k) V k. (24)
24 Deceleration Parameter and Neutrino Fugacity Υ T k = 0.3m e T k = m e q Υ q T/m e [1]J. Birrell, C.T. Yang, P. Chen, and J. Rafelski, Mod. Phys. Lett. A 28, (2013)
25 Effective Number of Neutrinos and Reheating Ratio Nν Nν = Time T k /m e [1]J. Birrell, C.T. Yang, P. Chen, and J. Rafelski, Phys. Rev. D 89, (2014) Nν = 3.36 Nν = N ν T γ /T ν Tγ/Tν
26 Reheating Ratio T γ = aυ b ( 1 + cσ 2 + O(σ 3 ) ), Υ = e σ T ν ( a = ) g 1/3 ( ) e 11 1/3 = , 8 g γ 4 b = π g ν+g e 8 g γ 3645 ζ(3) 2 g ν 8π 6 g γ 27ζ(3) g e , c g γ (25a) (25b) (25c) (25d) (25e)
27 Reheating Ratio T γ T ν = aυ b ( 1 + cσ 2 + O(σ 3 ) ) (26) The first order approximation has a maximum relative error of 2% relative error) in the region of interest.4 Υ ν 1. The second order approximation brings the relative error down to less than Power law behavior first observed numerically in [1] [1]J. Birrell, C.T. Yang, P. Chen, and J. Rafelski, Mod. Phys. Lett. A 28, (2013)
28 Effective Number of Neutrinos N ν = 360 7π 4 e 4bσ (1 + cσ 2 ) 4 0 u 3 e u σ + 1 du ( 1 + O(σ 3 ) ). The relative error of this approximation is less than over the range 0.9 σ 0. (27)
29 Equation of State First physically consistent derivation of equation of state of free-streaming neutrinos, including dependence on N ν and neutrino mass (β = m ν /T γ ). Fugacity is a crucial aspect. ρ EV /ρ 0 = N ν βi δN ν i i δNν 2 βi βi 4, i i P EV /P 0 = N ν βi δN ν i i δNν 2 βi βi 4. i i β 2 i β 2 i
30 Neutrinos Today number density [1/cm 3 ] neutrinos (exact) neutrinos (linearized) photons N ν Figure: Conservation of comoving neutrinos gives present day neutrino number density.
31 m ν = 0.05 ev m ν = 0.25 ev m ν = 0.5 ev m ν = 1 ev m ν = 2 ev fv v [km/s] Figure: Normalized neutrino velocity distribution in the Earth frame for N ν = (solid lines) and N ν = 3.62 (dashed lines). [1]J. Birrell and J. Rafelski arxiv: [hep-ph]
32 Part 2 Summary The effective number of neutrinos is probed by BBN element abundances and the CMB. N ν larger than predicted by two body scattering? We performed a model independent parametric study of the dependence of Υ, T γ /T ν, and N ν on T k < T ch. Measurement of N ν constitutes a measurement of T k. N ν > indicates that there are significant factor(s) other than two body scattering in neutrino freeze-out. Characterization of the effects of reheating on neutrino spectrum today is a small step towards detector design. Drag force on a millimeter sized coherent detector is encouraging, considering precision force measurements on the order of yocto-newtons (10 24 N).
33 Outline Part 1: Background Relic Neutrinos Overview of Cosmology Part 2: Study of Neutrino Distribution using Conservation Laws Chemical vs Kinetic Equilibrium Neutrino Freeze-out Neutrinos today Part 3: General Relativistic Boltzmann Equation Boltzmann Equation and Scattering Operators Neutrino Freeze-out Revisited
34 Part 3: Boltzmann Equation in an FRW Spacetime Particle dynamics given by Boltzmann equation p α x αf Γ j µνp µ p ν p jf = C[f ]. (29) No longer assuming chemical or kinetic equilibrium. For a general distribution in an FRW spacetime the Boltzmann equation is t f ph p f = 1 ȧ C[f ], H = E a For 2 2 reactions the collision operator is C[f 1 ] = 1 2 F(p 1, p 2, p 3, p 4 )S M 2 (2π) 4 δ( p) 4 i=2 (30) d 3 p i 2(2π) 3 E i, (31) F =f 3 (p 3 )f 4 (p 4 )f 1 (p 1 )f 2 (p 2 ) f 1 (p 1 )f 2 (p 2 )f 3 (p 3 )f 4 (p 4 ), (32) f i =1 f i (Fermi blocking factors).
35 Neutrino Scattering Matrix Elements Process S M 2 ν e + ν e ν e + ν e 128G 2 F (p 1 p 4 )(p 2 p 3 ) ν e + ν e ν e + ν e 64G 2 F (p 1 p 2 )(p 3 p 4 ) ν e + ν e ν j + ν j 32G 2 F (p 1 p 4 )(p 2 p 3 ) ν e + ν j ν e + ν j 32G 2 F (p 1 p 4 )(p 2 p 3 ) ν e + ν j ν e + ν j 32G 2 F (p 1 p 2 )(p 3 p 4 ) ν e + ν e e + + e 128G 2 F [g2 L (p 1 p 4 )(p 2 p 3 ) + g 2 R (p 1 p 3 )(p 2 p 4 ) + g L g R m 2 e(p 1 p 2 )] ν e + e ν e + e 128G 2 F [g2 L (p 1 p 2 )(p 3 p 4 ) + g 2 R (p 1 p 4 )(p 2 p 3 ) g L g R m 2 e(p 1 p 3 )] ν e + e + ν e + e + 128G 2 F [g2 R (p 1 p 2 )(p 3 p 4 ) + g 2 L (p 1 p 4 )(p 2 p 3 ) g L g R m 2 e (p 1 p 3 )] Table: Matrix elements for electronic neutrino processes where j = µ, τ, g L = sin2 θ W, g R = sin 2 θ W, sin 2 (θ W ) = 0.23 is the Weinberg angle, and G 2 F = GeV 2 is Fermi s constant. [1] A. Dolgov, S.H. Hansen, and D.V. Semikoz, Nucl.Phys. B543 (1999) 269.
36 Neutrino Scattering Matrix Elements Process S M 2 ν i + ν i ν i + ν i 128G 2 F (p 1 p 4 )(p 2 p 3 ) ν i + ν i ν i + ν i 64G 2 F (p 1 p 2 )(p 3 p 4 ) ν i + ν i ν j + ν j 32G 2 F (p 1 p 4 )(p 2 p 3 ) ν i + ν j ν i + ν j 32G 2 F (p 1 p 4 )(p 2 p 3 ) ν i + ν j ν i + ν j 32G 2 F (p 1 p 2 )(p 3 p 4 ) ν i + ν i e + + e 128G 2 F [ g2 L (p 1 p 4 )(p 2 p 3 ) + g 2 R (p 1 p 3 )(p 2 p 4 ) + g L g R m 2 e (p 1 p 2 )] ν i + e ν i + e 128G 2 F [ g2 L (p 1 p 2 )(p 3 p 4 ) + g 2 R (p 1 p 4 )(p 2 p 3 ) g L g R m 2 e (p 1 p 3 )] ν i + e + ν i + e + 128G 2 F [g2 R (p 1 p 2 )(p 3 p 4 ) + g 2 L (p 1 p 4 )(p 2 p 3 ) g L g R m 2 e(p 1 p 3 )] Table: Matrix elements for µ and τ neutrino processes where i = µ, τ, j = e, µ, τ, j i, g L = g L 1 = sin2 θ W, g R = sin 2 θ W, sin 2 (θ W ) = 0.23 is the Weinberg angle, and G 2 F = GeV 2 is Fermi s constant. [1] A. Dolgov, S.H. Hansen, and D.V. Semikoz, Nucl.Phys. B543 (1999) 269.
37 Chemical Equilibrium Spectral Method Has been solved via a spectral method based on the orthonormal basis of polynomials generated by the weight [1], [2] f ch = 1 e y, y = a(t)p. (33) + 1 Similar to method of moments of Grad in hydrodynamics, except with ultra-relativistic fermions. Implicit assumptions Near chemical equilibrium (Υ 1). Reheating is insignificant (T scales as 1/a). [1] S. Esposito et. al., Nucl. Phys. B590 (2000) 539 [2] G. Mangano et. al., Phys. Lett. B534 (2002) 8
38 Chemical Non-Equilibrium Spectral Method Change variables z = p/t(t), T is effective temperature. ( t f z H + Ṫ ) z f = 1 C[f ]. (34) T E To model a distribution near kinetic equilibrium decompose f (t, z) = f Υ (z, Υ(t))ψ(t, z), f Υ (z, Υ) = 1 Υ 1 e z + 1. (35) We will solve the Boltzmann equation using a dynamical basis of orthogonal polynomials that is adapted to systems where ψ 1.
39 Chemical Non-Equilibrium Spectral Method Expand ψ(t, z) in the orthonormal basis of polynomials, ˆψ(z, Υ), generated by the Υ-dependent weight via the recurrence relation [1] w Υ = z 2 f Υ (z, Υ). (36) ψ 0 = 1, ψ 1 = ψ 0 (x x ˆψ 0, ˆψ 0 ) ˆψ 0, (37) [( ψ n+1 = ψ n x x ˆψ n, ˆψ ) n ˆψn x ˆψ n, ˆψ n 1 ˆψ ] n 1. (38) [1] See e.g. F. Olver, Asymptotics and Special Functions. New York: Academic Press, 1974.
40 Chemical Non-Equilibrium Spectral Method Expand ψ(t, z) in the orthonormal basis of polynomials, ˆψ(z, Υ) ψ(t, z) = i b i (t) ˆψ i (z, Υ(t)) (39) ḃ k = i i ( b i H + Ṫ ) ( z ) ˆψi z f Υ, T f ˆψ k + z z ˆψi, ˆψ k Υ b i Υ ( 1 f Υ f Υ Υ ˆψ i, ˆψ k + ˆψ i Υ, ˆψ k ) + 1 f Υ E C[f ], ˆψ k [1] J. Birrell and J. Rafelski, arxiv:
41 Chemical Non-Equilibrium Spectral Method Up to this point, T and Υ were arbitrary functions. Fix the dynamics of T(t) and Υ(t) by imposing In other words b 0 (t) ˆψ 0 (z, Υ(t)) = 1, b 1 (t) = 0. (40) f (t, z) = f Υ (t, z)(1 + φ(t, z)), φ = N b i ˆψi. (41) i=2 The non-thermal distortions are contained in φ.
42 Chemical Non-Equilibrium Spectral Method These condition lead to evolution equations for T(t), Υ(t) (Ab) Υ/Υ 1 1 f = C[f ], ˆψ ΥE 0 (Ab) 0 1 f C[f ], ˆψ ΥE 1 [Υ Υ 1, 1 /(2 ψ 0 ) + (Bb) 0 ](Ab) 1 (Ab) 0 (Bb) 1, (42) Ṫ/T = (Bb)1 1 f C[f ], ˆψ ΥE 0 1 f C[f ], ˆψ ΥE 1 [Υ Υ 1, 1 /(2 ψ 0 ) + (Bb) 0 ] [Υ Υ 1, 1 /(2 ψ 0 ) + (Bb) 0 ](Ab) 1 (Ab) 0 (Bb) 1 H = 1 ( ) (Ab) 1 (Bb) 1 1 Υ/Υ f Υ E C[f ], ˆψ 1 H. (43) where (Ab) n = N j=0 An j bj and similarly for B.
43 Chemical Non-Equilibrium Spectral Method Physical motivation for the weight function: n = g νt 3 2π 2 1, 1 = g νt 3 2π 2 ρ = g ν 2π 2 T4 1, z = g νt 4 2π f Υ z 2 dz, f Υ z 3 dz. (44a) (44b) The particle number and energy densities are captured by f Υ and the non-thermal distortions from kinetic equilibrium are entirely contained in φ. Only two degrees of freedom are necessary. [1] J. Birrell and J. Rafelski, arxiv:
44 Chemical Non-Equilibrium Spectral Method n = g νt 3 2π 2 1, 1 = g νt 3 2π 2 ρ = g ν 2π 2 T4 1, z = g νt 4 2π f Υ z 2 dz, f Υ z 3 dz. (45a) (45b) Only the zeroth order term f Υ is used here, not the full distribution i.e. computed as if in kinetic equilibrium. This is only possible with the correct evolution of T(t), Υ(t) shown earlier. It is in this sense that T, Υ are the averaged temperature and fugacity of the non-equilibrium distribution f.
45 Key Differences and Improvements Adapted to kinetic equilibrium (i.e. non-thermal distortions φ small) but allows for arbitrary Υ. Generalized averaged temperature (i.e. T not required to scale as 1/a). Introduction of z 2 factor means only 2 modes are required to capture particle number and energy flow, compared to 4 for the chemical equilibrium method. Numerical evaluation of scattering integrals dominates the numerical cost, so this last point results in a large reduction in computer time.
46 Generalizations In general, one can introduce additional parameters in the weight to capture additional properties/impose conditions. For example pressure for massive particles local flow direction for non-isotropic distributions (e.g. hydrodynamics for fermions) We have begun work on the generalization to massive fermions/bosons.
47 Neutrino Scattering Integrals The collision operator inner products must be analytically reduced to fewer dimensions for numerical computation M T3 2π 2 1 f Υ E C[f ], ˆψ k = ˆψ k (p 1 )C[f 1 ](p 1 ) = G k (p 1, p 2, p 3, p 4 )S M 2 (2π) 4 δ( p) d 3 p 1 (2π) 3 E 1 4 i=1 δ 0 (p 2 i m 2 i ) d4 p i (2π) 3, (46) G k = ˆψ k (z 1 ) [ f 3 (p 3 )f 4 (p 4 )f 1 (p 1 )f 2 (p 2 ) f 1 (p 1 )f 2 (p 2 )f 3 (p 3 )f 4 (p 4 ) ], (47) f i = 1 f i. (48)
48 Neutrino Scattering Integrals Integral is initially 16 dimensional. Mass shell delta functions reduce it to 12 dimensions. 4-momentum conservation reduces it to 8 dimensions. Spend symmetries to reduce the dimensionality further. Using isotropy, we reduced this to three iterated integrals and improved on previous methods by obtaining a smooth integrand over a simple integration region to numerically integrate. Numerically, the matrix element inner products computed by our method agree with the method used in [1]. [1] A. Dolgov, S.H. Hansen, and D.V. Semikoz, Nucl.Phys. B543 (1999) 269.
49 K(s, p) = 8πrr s 1 M = 1 [ (2π) 7 1 ( 2π 0 s 0 0 K(s, p) p2 dpds, (49) p0 S M 2 (s, t(cos(ψ) ) 1 y 2 1 z2 + yz))dψ ] (50) G 34((q ) 0 α r δy)dy G 12(q 0 α rδz)dz. (51) For neutrino reactions, K splits into a linear combination of products of 1-d integrals. p 0 = p 2 + s, α = p0, δ = p, q 0 = m2 1 m2 2, (q ) 0 = m2 3 m2 4, (52) s s s s r = 1 s (s (m 1 + m 2 ) 2 )(s (m 1 m 2 ) 2 ), r = 1 s (s (m 3 + m 4 ) 2 )(s (m 3 m 4 ) 2 ), (53) t(x) = 1 4 ((q0 (q ) 0 ) 2 r 2 (r ) 2 + 2rr x), s 0 = max{(m 1 + m 2 ) 2, (m 3 + m 4 ) 2 }. (54) F(p, q, q ) =F 12 ((p + q) U/2, (p q) U/2)F 34 ((p + q ) U/2, (p q ) U/2) (55) G 12 (p U, q U)G 34 (p U, q U).
50 Solving the Boltzmann Equation We tested the chemical non-equilibrium method against the chemical equilibrium method to show their equivalence. Define the following z fin = T γ /a, ρ ν0 = π2 a 4, δ ρ ν = ρ ν ρ ν0 1. (56) Method Modes z fin δ ρ νe δ ρ νµ,τ N ν Chemical Eq Chemical Non-Eq Chemical Non-Eq Chemical non-equilibrium method with minimum number of modes 20 faster.
51 Applications: Dependence on Weinberg Angle Parameter that impacts strength of weak force processes. Measured value in vacuum θ W 30. cos(θ W ) = m W m Z. (57) Value not determined by Standard Model. Neutrino freeze-out is a probe of possible time or temperature dependence.
52 Dependence on Weinberg Angle Nν sin 2 (θ W ) Figure: Effective number of neutrinos as a function of Weinberg angle. Vertical line is sin 2 (θ W ) =.23.
53 Applications: Dependence on Interaction Strength The Einstein Boltzmann system depends on a single dimensionless parameter,which we call the interaction strength, η M p m 3 eg 2 F. (58) Neutrino freeze-out is a probe of possible time or temperature dependence of this combination of natural constants. We let η be the presently measured value.
54 Dependence on Interaction Strength 10 0 Nν η/η 0 Figure: Effective number of neutrinos as a function of interaction strength ratio.
55 Part 3 Summary We developed an improved method for analytically simplifying the scattering integrals that substantially reduced the numerical integration cost. We developed a spectral method for solving the Boltzmann equation that is adapted to chemical non-equilibrium. The method incorporates three key improvements: a modified weight function together with a dynamical effective temperature and fugacity. We investigated impact of varying sin 2 (θ W ) on neutrino freeze-out. Maximum value of N ν = 3.13 for sin 2 (θ W ) = 1. We investigated impact of varying the dimensionless parameter η = M p m 3 eg 2 F on neutrino freeze-out. N ν = 3.30 for η/η 0 = 36.3 and N ν = 3.62 for η/η 0 =
56 Looking to the Future Investigate systems of massive particles. Extend to spatially varying systems (generalized hydrodynamics for bosons/fermions). We focused on algorithmic optimization. Can also look at code optimization (will be more critical for spatially dependent systems).
57 Acknowledgments Thanks to my advisor Dr. Rafelski for his invaluable physical insight, providing the physical compass for our work. Thanks to Cheng-Tao Yang, a student from National Taiwan University who visited U of A and worked on evaluating matrix elements for neutrino processes, supervisor Dr. Pisin Chen. Special thanks to my parents for their unwavering emotional support. This work has been supported by the Department of Defense (DoD) through the National Defense Science & Engineering Graduate Fellowship (NDSEG) Program. Supported in part by a grant from the U.S. Department of Energy, DE-FG02-04ER41318 (PI Dr. Johann Rafelski).
58 10 0 Energy Density Fraction Dark Energy Dark Matter Baryons e ± γ ν µ ± π ±, T [ev] Figure: Current era: 69% dark energy, 26% dark matter, 5% baryons, < 1% photons and neutrinos, 1 massless and 2.1 ev neutrinos.
59 Effective Number of Neutrinos Nν = Nν = 3.62 Nν 3.5 Nν Nν = Nν = Nν = T γ /T ν Nν = σ [1]J. Birrell, C.T. Yang, P. Chen, and J. Rafelski, Phys. Rev. D 89, (2014)
60 Neutrinos Today Neutrino distribution in Earth rest frame: f (p µ 1 ) =. (59) Υ e 1 (p µ U µ) 2 m 2 ν /Tν + 1 The 4-vector characterizing the rest frame of the neutrino distribution: U µ = (γ, 0, 0, V γ), γ = 1/ 1 V 2, V = 300 ± 30 km/s c (60) Note that spatial inhomogeneity of gravity and of the neutrino distribution are not included in our considerations.
61 m ν = 0.05 ev m ν = 0.25 ev m ν = 0.5 ev m ν = 1 ev m ν = 2 ev fv v [km/s] Figure: Normalized neutrino velocity distribution in the Earth frame for N ν = (solid lines) and N ν = 3.62 (dashed lines). [1]J. Birrell and J. Rafelski arxiv: [hep-ph]
62 m ν = 0.05 ev m ν = 0.1 ev m ν = 0.25 ev m ν = 0.5 ev m ν = 1 ev m ν = 2 ev fe T 1 2 m νv 2 [µev] Figure: Neutrino energy distribution in the Earth frame for N ν = (solid lines) and N ν = 3.62 (dashed lines). [1]J. Birrell and J. Rafelski arxiv: [hep-ph]
63 10 0 m ν = 0.05 ev m ν = 0.25 ev m ν = 0.5 ev m ν = 1 ev m ν = 2 ev fλ λ [mm] Figure: Neutrino de Broglie wavelength distribution in the Earth frame for N ν = (solid lines) and N ν = 3.62 (dashed lines). [1]J. Birrell and J. Rafelski arxiv: [hep-ph]
64 10 2 m ν = 0.05 ev m ν = 0.25 ev m ν = 0.5 ev m ν = 1 ev m ν = 2 ev f λ 3.62 /f λ λ [mm] Figure: Ratio of neutrino de Broglie wavelength distribution with N ν = 3.62 to the distribution with N ν = [1]J. Birrell and J. Rafelski arxiv: [hep-ph]
65 Neutrino Drag Drag force on spherical detector due to hard sphere scattering 1 dp A dt = 2 π f (p)p 2 r d 3 p cos φ sin φdφẑ, m ν 0 p r<0 1 f (p) =, (61) Υ e 1 (E V p ẑ) 2 γ 2 m 2 ν/t ν + 1 p ẑ = p r cos φ p φ sin φ.
66 Least squares fit to drag force in hard sphere limit, valid to approximately 1% relative error in the region 270 km/s V 330 km/s; 3 N ν 3.7; 0.05 ev m ν 2 ev. Define δn ν = N ν 3, δv = V /300 km/s 1 and recalling the CMB temperature T CMB = mev, x = m ν /10 4 T CMB. F drag /A zn/mm 2 = K 1 + K 2 δv + K 3 δn ν + K 4 δv δn ν, K 1 (x) = x x x x 8, K 2 (x) = x x x 6, K 3 (x) = x x x 6, K 4 (x) = x x 4. [1]J. Birrell and J. Rafelski arxiv: [hep-ph]
67 Fdrag/A [zn/mm 2 ] m ν = 0.05 ev m ν = 0.25 ev m ν = 0.5 ev m ν = 1 ev m ν = 2 ev V [km/s] Figure: Drag per unit area on spherical detector due to hard sphere neutrino scattering for N ν = (solid) and N ν = 3.62 (dashed). [1]J. Birrell and J. Rafelski arxiv: [hep-ph]
68 Scattering Length We want to define a characteristic length between scatterings. We have 1 d a 3 dt (a3 n) + = g ν 2π 2 T2 C[f ]zdz, (63) where the plus indicates that we want the one way flow and hence the one way scattering operator (denoted by arrow) in order to count the number of reactions and not just the net flow. Define the scattering length a 3 n L v d dt (a3 n) = 1 0 Υ 1 e z +1 z2 dz 0 C[f ]z 2 /Edz. (64)
69 Scattering Length νe Scattering Length [m] 1e5 1 1e 05 1e 10 1e 15 Annihilation Length ν e ± Scattering Length ν ν Scattering Length c/h νµ,τ Scattering Length [m] 1e5 1 1e 05 1e 10 1e 15 Annihilation Length ν e ± Scattering Length ν ν Scattering Length c/h T γ [MeV] T γ [MeV] Figure: Comparison of Hubble parameter to neutrino scattering length for various types of processes.
70 Dependence on Weinberg Angle δρν Υ ν e ν µ,τ ν e ν µ,τ sin 2 (θ W ) sin 2 (θ W ) Figure: Fractional increase in neutrino energy and neutrino fugacities, as functions of Weinberg angle. Vertical line is sin 2 (θ W ) =.23.
71 Dependence on Weinberg Angle atγ Tγ/Tν ν e ν µ,τ sin 2 (θ W ) sin 2 (θ W ) Figure: Photon reheating and neutrino to photon temperature ratios, as functions of Weinberg angle. Vertical line is sin 2 (θ W ) =.23.
72 Dependence on Weinberg Angle νe Freeze-out Temperatures [MeV] ν e ± Chemical Freeze-out ν e ± Kinetic Freeze-out ν ν Kinetic Freeze-out sin 2 (θ W ) νµ,τ Freeze-out Temperatures [MeV] ν e ± Chemical Freeze-out ν e ± Kinetic Freeze-out ν ν Kinetic Freeze-out sin 2 (θ W ) Figure: Freeze-out temperatures for various types of processes, as functions of Weinberg angle. Vertical line is sin 2 (θ W ) =.23.
73 Dependence on Interaction Strength ν e ν µ,τ δρν 0.1 Υ ν e ν µ,τ η/η η/η 0 Figure: Fractional increase in neutrino energy and neutrino fugacities, as functions of interaction strength ratio.
74 Dependence on Interaction Strength ν e ν µ,τ 1.39 atγ Tγ/Tν η/η 0 η/η 0 Figure: Photon reheating and neutrino to photon temperature ratios, as functions of interaction strength ratio.
75 Dependence on Interaction Strength νe Freeze-out Temperatures [MeV] ν e ± Chemical Freeze-out ν e ± Kinetic Freeze-out ν ν Kinetic Freeze-out η/η 0 νµ,τ Freeze-out Temperatures [MeV] ν e ± Chemical Freeze-out ν e ± Kinetic Freeze-out ν ν Kinetic Freeze-out η/η 0 Figure: Freeze-out temperatures for various types of processes, as functions of interaction strength ratio.
Time Evolution of the Hot Hagedorn Universe
Time Evolution of the Results obtained in collaboration with Jeremiah Birrell The University of Arizona 1965: Penzias and Wilson discover Alpher-Gamov CMB 1966-1968: Hot Big-Bang becoming conventional
More informationIntroduction: Cosmic Neutrinos Dark Radiation and the QGP Era Darkness Production: Universe Laboratory
Outline of the talk 1. Introduction: CMB fluctuation analysis observes N ν, number of invisible Cosmic Neutrinos pushing the Universe apart Anything else out there adding to the pressure? 2. Impact of
More informationPAPER 71 COSMOLOGY. Attempt THREE questions There are seven questions in total The questions carry equal weight
MATHEMATICAL TRIPOS Part III Friday 31 May 00 9 to 1 PAPER 71 COSMOLOGY Attempt THREE questions There are seven questions in total The questions carry equal weight You may make free use of the information
More informationBoiling Quarks, Melting Hadrons
Department of Physics, The University of Arizona Tucson, AZ June 2, 2015 1964/65: Two new fundamental ideas Quarks Standard Model of Particle Physics Hagedorn Temperature New State of Elementary Matter
More informationCosmology Neutrinos Connect to QGP Dark Radiation and the QGP Era Darkness Production: Universe Laboratory
Cosmology Neutrinos Connect to QGP Dark Radiation and the QGP Era Darkness Production: Universe Laboratory What has CMB to do with QGP? Figure: Photons freeze-out around 0.25 ev and today they make up
More informationQGP in the Universe and in the Laboratory
QGP in the Universe and in the Laboratory Results obtained in collaboration with Jeremiah Birrell, Michael Fromerth, Inga Kuznetsowa, Michal Petran Former Students at The University of Arizona Johann Rafelski,
More informationWe can check experimentally that physical constants such as α have been sensibly constant for the past ~12 billion years
² ² ² The universe observed ² Relativistic world models ² Reconstructing the thermal history ² Big bang nucleosynthesis ² Dark matter: astrophysical observations ² Dark matter: relic particles ² Dark matter:
More informationThe early and late time acceleration of the Universe
The early and late time acceleration of the Universe Tomo Takahashi (Saga University) March 7, 2016 New Generation Quantum Theory -Particle Physics, Cosmology, and Chemistry- @Kyoto University The early
More informationAY127 Problem Set 3, Winter 2018
California Institute of Technology AY27 Problem Set 3, Winter 28 Instructor: Sterl Phinney & Chuck Steidel TA: Guochao (Jason) Sun February 7, 28 Problem Detailed Solution : (a) We work in natural units
More informationCosmology: An Introduction. Eung Jin Chun
Cosmology: An Introduction Eung Jin Chun Cosmology Hot Big Bang + Inflation. Theory of the evolution of the Universe described by General relativity (spacetime) Thermodynamics, Particle/nuclear physics
More informationCosmology ASTR 2120 Sarazin. Hubble Ultra-Deep Field
Cosmology ASTR 2120 Sarazin Hubble Ultra-Deep Field Cosmology - Da Facts! 1) Big Universe of Galaxies 2) Sky is Dark at Night 3) Isotropy of Universe Cosmological Principle = Universe Homogeneous 4) Hubble
More informationQGP in the Universe and in the Laboratory
QGP in the Universe and in the Laboratory Results obtained in collaboration with Jeremiah Birrell, Michael Fromerth, Inga Kuznetsowa, Michal Petran Graduate Students at The University of Arizona Johann
More informationPROBLEM SET 10 (The Last!)
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.286: The Early Universe December 5, 2013 Prof. Alan Guth PROBLEM SET 10 (The Last!) DUE DATE: Tuesday, December 10, 2013, at 5:00 pm.
More informationMATHEMATICAL TRIPOS Part III PAPER 53 COSMOLOGY
MATHEMATICAL TRIPOS Part III Wednesday, 8 June, 2011 9:00 am to 12:00 pm PAPER 53 COSMOLOGY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY
More information7 Relic particles from the early universe
7 Relic particles from the early universe 7.1 Neutrino density today (14 December 2009) We have now collected the ingredients required to calculate the density of relic particles surviving from the early
More informationCosmological Issues. Consider the stress tensor of a fluid in the local orthonormal frame where the metric is η ab
Cosmological Issues 1 Radiation dominated Universe Consider the stress tensor of a fluid in the local orthonormal frame where the metric is η ab ρ 0 0 0 T ab = 0 p 0 0 0 0 p 0 (1) 0 0 0 p We do not often
More informationAstronomy, Astrophysics, and Cosmology
Astronomy, Astrophysics, and Cosmology Luis A. Anchordoqui Department of Physics and Astronomy Lehman College, City University of New York Lesson VI March 15, 2016 arxiv:0706.1988 L. A. Anchordoqui (CUNY)
More informationThe Influence of DE on the Expansion Rate of the Universe and its Effects on DM Relic Abundance
The Influence of DE on the Expansion Rate of the Universe and its Effects on DM Relic Abundance Esteban Jimenez Texas A&M University XI International Conference on Interconnections Between Particle Physics
More informationLearning from WIMPs. Manuel Drees. Bonn University. Learning from WIMPs p. 1/29
Learning from WIMPs Manuel Drees Bonn University Learning from WIMPs p. 1/29 Contents 1 Introduction Learning from WIMPs p. 2/29 Contents 1 Introduction 2 Learning about the early Universe Learning from
More informationCHAPTER 3 THE INFLATIONARY PARADIGM. 3.1 The hot Big Bang paradise Homogeneity and isotropy
CHAPTER 3 THE INFLATIONARY PARADIGM Ubi materia, ibi geometria. Johannes Kepler 3.1 The hot Big Bang paradise In General Relativity, the Universe as a whole becomes a dynamical entity that can be modeled
More informationObservational evidence and cosmological constant. Kazuya Koyama University of Portsmouth
Observational evidence and cosmological constant Kazuya Koyama University of Portsmouth Basic assumptions (1) Isotropy and homogeneity Isotropy CMB fluctuation ESA Planck T 5 10 T Homogeneity galaxy distribution
More informationTraveling Through the Universe: Back in Time to the Quark-Gluon Plasma Era
Journal of Physics: Conference Series OPEN ACCESS Traveling Through the Universe: Back in Time to the Quark-Gluon Plasma Era To cite this article: Johann Rafelski and Jeremiah Birrell 2014 J. Phys.: Conf.
More informationIs inflation really necessary in a closed Universe? Branislav Vlahovic, Maxim Eingorn. Please see also arxiv:
Is inflation really necessary in a closed Universe? Branislav Vlahovic, Maxim Eingorn North Carolina Central University NASA University Research Centers, Durham NC Please see also arxiv:1303.3203 Chicago
More informationEquilibration and decoupling of a relativistic gas in a Friedmann-Robertson-Walker spacetime
Equilibration and decoupling of a relativistic gas in a Friedmann-Robertson-Walker spacetime Juan M. Torres-Rincon (Frankfurt Institute for Advanced Studies) in collaboration with J. Tindall, J.-B. Rosé,
More informationModern Cosmology / Scott Dodelson Contents
Modern Cosmology / Scott Dodelson Contents The Standard Model and Beyond p. 1 The Expanding Universe p. 1 The Hubble Diagram p. 7 Big Bang Nucleosynthesis p. 9 The Cosmic Microwave Background p. 13 Beyond
More informationWIMP Dark Matter and the QCD Equation of State
Doktoratskolleg Graz Karl-Franzens-Universität, May 2006 WIMP Dark Matter and the QCD Equation of State Owe Philipsen Universität Münster Motivation + overview Cosmology I: expansion and contents of the
More informationQUIZ 3. Quiz Date: December 7, 2016
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.286: The Early Universe December 21, 2016 Prof. Alan Guth QUIZ 3 Quiz Date: December 7, 2016 Reformatted to Remove Blank Pages Please
More information14 Lecture 14: Early Universe
PHYS 652: Astrophysics 70 14 Lecture 14: Early Universe True science teaches us to doubt and, in ignorance, to refrain. Claude Bernard The Big Picture: Today we introduce the Boltzmann equation for annihilation
More informationHot Big Bang model: early Universe and history of matter
Hot Big Bang model: early Universe and history of matter nitial soup with elementary particles and radiation in thermal equilibrium. adiation dominated era (recall energy density grows faster than matter
More informationCosmology and particle physics
Cosmology and particle physics Lecture notes Timm Wrase Lecture 5 The thermal universe - part I In the last lecture we have shown that our very early universe was in a very hot and dense state. During
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.286: The Early Universe October 27, 2013 Prof. Alan Guth PROBLEM SET 6
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.86: The Early Universe October 7, 013 Prof. Alan Guth PROBLEM SET 6 DUE DATE: Monday, November 4, 013 READING ASSIGNMENT: Steven Weinberg,
More informationLecture 2: Cosmological Background
Lecture 2: Cosmological Background Houjun Mo January 27, 2004 Goal: To establish the space-time frame within which cosmic events are to be described. The development of spacetime concept Absolute flat
More informationVU lecture Introduction to Particle Physics. Thomas Gajdosik, FI & VU. Big Bang (model)
Big Bang (model) What can be seen / measured? basically only light _ (and a few particles: e ±, p, p, ν x ) in different wave lengths: microwave to γ-rays in different intensities (measured in magnitudes)
More informationOddities of the Universe
Oddities of the Universe Koushik Dutta Theory Division, Saha Institute Physics Department, IISER, Kolkata 4th November, 2016 1 Outline - Basics of General Relativity - Expanding FRW Universe - Problems
More informationCosmology (Cont.) Lecture 19
Cosmology (Cont.) Lecture 19 1 General relativity General relativity is the classical theory of gravitation, and as the gravitational interaction is due to the structure of space-time, the mathematical
More informationGeneral Relativity Lecture 20
General Relativity Lecture 20 1 General relativity General relativity is the classical (not quantum mechanical) theory of gravitation. As the gravitational interaction is a result of the structure of space-time,
More informationIntroduction to Inflation
Introduction to Inflation Miguel Campos MPI für Kernphysik & Heidelberg Universität September 23, 2014 Index (Brief) historic background The Cosmological Principle Big-bang puzzles Flatness Horizons Monopoles
More informationChapter 22 Lecture. The Cosmic Perspective. Seventh Edition. The Birth of the Universe Pearson Education, Inc.
Chapter 22 Lecture The Cosmic Perspective Seventh Edition The Birth of the Universe The Birth of the Universe 22.1 The Big Bang Theory Our goals for learning: What were conditions like in the early universe?
More informationPRECISE NUCLEOSYNTHESIS LIMITS ON NEUTRINO MASSES. Kimmo Kainulainen. CERN, CH-1211, Geneva 23, Switzerland
PRECISE NUCLEOSYNTHESIS LIMITS ON NEUTRINO MASSES Kimmo Kainulainen CERN, CH-1211, Geneva 23, Switzerland A computation of nucleosynthesis bounds on the masses of long-lived Dirac and Majorana neutrinos
More informationThermalization of axion dark matter
Thermalization of axion dark matter Ken ichi Saikawa ICRR, The University of Tokyo Collaborate with M. Yamaguchi (Tokyo Institute of Technology) Reference: KS and M. Yamaguchi, arxiv:1210.7080 [hep-ph]
More informationTHERMAL HISTORY OF THE UNIVERSE
M. Pettini: Introduction to Cosmology Lecture 7 THERMAL HISTORY OF THE UNIVERSE The Universe today is bathed in an all-pervasive radiation field, the Cosmic Microwave Background (CMB) which we introduced
More information2.1 Basics of the Relativistic Cosmology: Global Geometry and the Dynamics of the Universe Part I
1 2.1 Basics of the Relativistic Cosmology: Global Geometry and the Dynamics of the Universe Part I 2 Special Relativity (1905) A fundamental change in viewing the physical space and time, now unified
More informationKinetic Theory of Dark Energy within General Relativity
Kinetic Theory of Dark Energy within General Relativity Author: Nikola Perkovic* percestyler@gmail.com University of Novi Sad, Faculty of Sciences, Institute of Physics and Mathematics Abstract: This paper
More informationDecaying Dark Matter, Bulk Viscosity, and Dark Energy
Decaying Dark Matter, Bulk Viscosity, and Dark Energy Dallas, SMU; April 5, 2010 Outline Outline Standard Views Dark Matter Standard Views of Dark Energy Alternative Views of Dark Energy/Dark Matter Dark
More informationPAPER 310 COSMOLOGY. Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight.
MATHEMATICAL TRIPOS Part III Wednesday, 1 June, 2016 9:00 am to 12:00 pm PAPER 310 COSMOLOGY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY
More informationFrom inflation to the CMB to today s universe. I - How it all begins
From inflation to the CMB to today s universe I - How it all begins Raul Abramo Physics Institute - University of São Paulo abramo@fma.if.usp.br redshift Very brief cosmic history 10 9 200 s BBN 1 MeV
More informationPHY 475/375. Lecture 5. (April 9, 2012)
PHY 475/375 Lecture 5 (April 9, 2012) Describing Curvature (contd.) So far, we have studied homogenous and isotropic surfaces in 2-dimensions. The results can be extended easily to three dimensions. As
More informationThe Four Basic Ways of Creating Dark Matter Through a Portal
The Four Basic Ways of Creating Dark Matter Through a Portal DISCRETE 2012: Third Symposium on Prospects in the Physics of Discrete Symmetries December 4th 2012, Lisboa Based on arxiv:1112.0493, with Thomas
More informationIoP. An Introduction to the Science of Cosmology. Derek Raine. Ted Thomas. Series in Astronomy and Astrophysics
Series in Astronomy and Astrophysics An Introduction to the Science of Cosmology Derek Raine Department of Physics and Astronomy University of Leicester, UK Ted Thomas Department of Physics and Astronomy
More informationASTROPHYSICAL PROPERTIES OF MIRROR DARK MATTER
16 December 2011 ASTROPHYSICAL PROPERTIES OF MIRROR DARK MATTER Paolo Ciarcelluti Motivation of this research We are now in the ERA OF PRECISION COSMOLOGY and... Motivation of this research We are now
More informationCosmology and particle physics
Cosmology and particle physics Lecture notes Timm Wrase Lecture 9 Inflation - part I Having discussed the thermal history of our universe and in particular its evolution at times larger than 10 14 seconds
More informationSupplement: Statistical Physics
Supplement: Statistical Physics Fitting in a Box. Counting momentum states with momentum q and de Broglie wavelength λ = h q = 2π h q In a discrete volume L 3 there is a discrete set of states that satisfy
More informationLecture 2: The First Second origin of neutrons and protons
Lecture 2: The First Second origin of neutrons and protons Hot Big Bang Expanding and cooling Soup of free particles + anti-particles Symmetry breaking Soup of free quarks Quarks confined into neutrons
More informationGeneral Relativity and Cosmology Mock exam
Physikalisches Institut Mock Exam Universität Bonn 29. June 2011 Theoretische Physik SS 2011 General Relativity and Cosmology Mock exam Priv. Doz. Dr. S. Förste Exercise 1: Overview Give short answers
More informationUNIVERSITY OF OSLO Faculty of Mathematics and Natural Sciences
UNIVERSITY OF OSLO Faculty of Mathematics and Natural Sciences Exam for AST5220 Cosmology II Date: Tuesday, June 4th, 2013 Time: 09.00 13.00 The exam set consists of 13 pages. Appendix: Equation summary
More informationAstr 2320 Tues. May 2, 2017 Today s Topics Chapter 23: Cosmology: The Big Bang and Beyond Introduction Newtonian Cosmology Solutions to Einstein s
Astr 0 Tues. May, 07 Today s Topics Chapter : Cosmology: The Big Bang and Beyond Introduction Newtonian Cosmology Solutions to Einstein s Field Equations The Primeval Fireball Standard Big Bang Model Chapter
More informationThe Unifying Dark Fluid Model
The Model Centre de Recherche Astrophysique de Lyon Invisible Universe Paris July 2nd, 2009 s Dark Matter Problem Dark Matter Dark Energy Dark Fluids? Different scales involved Galactic scale Galaxy Rotation
More informationCosmological Issues. Consider the stress tensor of a fluid in the local orthonormal frame where the metric is η ab
Cosmological Issues Radiation dominated Universe Consider the stress tensor of a fluid in the local orthonormal frame where the metric is η ab ρ 0 0 0 T ab = 0 p 0 0 0 0 p 0 () 0 0 0 p We do not often
More informationTHE DARK SIDE OF THE COSMOLOGICAL CONSTANT
THE DARK SIDE OF THE COSMOLOGICAL CONSTANT CAMILO POSADA AGUIRRE University of South Carolina Department of Physics and Astronomy 09/23/11 Outline 1 Einstein s Greatest Blunder 2 The FLRW Universe 3 A
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.286: The Early Universe November 12, 2013 Prof. Alan Guth PROBLEM SET 7
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.286: The Early Universe November 12, 2013 Prof. Alan Guth PROBLEM SET 7 DUE DATE: Friday, November 15, 2013 READING ASSIGNMENT: Steven
More informationPOST-INFLATIONARY HIGGS RELAXATION AND THE ORIGIN OF MATTER- ANTIMATTER ASYMMETRY
POST-INFLATIONARY HIGGS RELAXATION AND THE ORIGIN OF MATTER- ANTIMATTER ASYMMETRY LOUIS YANG ( 楊智軒 ) UNIVERSITY OF CALIFORNIA, LOS ANGELES (UCLA) DEC 27, 2016 NATIONAL TSING HUA UNIVERSITY OUTLINE Big
More informationPar$cle Astrophysics
Par$cle Astrophysics Produc$on (Early Universe) Signatures (Large Scale Structure & CMB) Accelerator Detector Neutrinos and Dark MaCer were produced in the early universe Star$ng Point: Cosmic Photons
More informationHeating up QGP: towards charm quark chemical equilibration
Heating up QGP: towards charm quark chemical equilibration Mikko Laine (University of Bern, Switzerland) 1 What is it? 2 Melting / recombination: Leptonic annihilation: q q l + l Chemical equilibration:
More informationCosmological Signatures of a Mirror Twin Higgs
Cosmological Signatures of a Mirror Twin Higgs Zackaria Chacko University of Maryland, College Park Curtin, Geller & Tsai Introduction The Twin Higgs framework is a promising approach to the naturalness
More informationTuesday: Special epochs of the universe (recombination, nucleosynthesis, inflation) Wednesday: Structure formation
Introduction to Cosmology Professor Barbara Ryden Department of Astronomy The Ohio State University ICTP Summer School on Cosmology 2016 June 6 Today: Observational evidence for the standard model of cosmology
More informationGravitation: Cosmology
An Introduction to General Relativity Center for Relativistic Astrophysics School of Physics Georgia Institute of Technology Notes based on textbook: Spacetime and Geometry by S.M. Carroll Spring 2013
More informationComputational Physics and Astrophysics
Cosmological Inflation Kostas Kokkotas University of Tübingen, Germany and Pablo Laguna Georgia Institute of Technology, USA Spring 2012 Our Universe Cosmic Expansion Co-moving coordinates expand at exactly
More information4 Evolution of density perturbations
Spring term 2014: Dark Matter lecture 3/9 Torsten Bringmann (torsten.bringmann@fys.uio.no) reading: Weinberg, chapters 5-8 4 Evolution of density perturbations 4.1 Statistical description The cosmological
More informationComputational Applications in Nuclear Astrophysics using JAVA
Computational Applications in Nuclear Astrophysics using JAVA Lecture: Friday 10:15-11:45 Room NB 7/67 Jim Ritman and Elisabetta Prencipe j.ritman@fz-juelich.de e.prencipe@fz-juelich.de Computer Lab: Friday
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 3 June, 2004 9 to 12 PAPER 67 PHYSICAL COSMOLOGY Attempt THREE questions. There are four questions in total. The questions carry equal weight. You may not start to
More informationAstronomy, Astrophysics, and Cosmology
Astronomy, Astrophysics, and Cosmology Luis A. Anchordoqui Department of Physics and Astronomy Lehman College, City University of New York Lesson IX April 12, 2016 arxiv:0706.1988 L. A. Anchordoqui (CUNY)
More informationIntroduction to Cosmology
Introduction to Cosmology João G. Rosa joao.rosa@ua.pt http://gravitation.web.ua.pt/cosmo LECTURE 2 - Newtonian cosmology I As a first approach to the Hot Big Bang model, in this lecture we will consider
More informationLeptogenesis via Higgs Condensate Relaxation
The Motivation Quantum Fluctuations Higgs Relaxation Leptogenesis Summary Leptogenesis via Higgs Condensate Relaxation Louis Yang Department of Physics and Astronomy University of California, Los Angeles
More informationTHE UNIVERSITY OF SYDNEY FACULTY OF SCIENCE INTERMEDIATE PHYSICS PHYS 2913 ASTROPHYSICS AND RELATIVITY (ADVANCED) ALL QUESTIONS HAVE THE VALUE SHOWN
CC0937 THE UNIVERSITY OF SYDNEY FACULTY OF SCIENCE INTERMEDIATE PHYSICS PHYS 2913 ASTROPHYSICS AND RELATIVITY (ADVANCED) SEMESTER 2, 2014 TIME ALLOWED: 2 HOURS ALL QUESTIONS HAVE THE VALUE SHOWN INSTRUCTIONS:
More informationGraceful exit from inflation for minimally coupled Bianchi A scalar field models
Graceful exit from inflation for minimally coupled Bianchi A scalar field models Florian Beyer Reference: F.B. and Leon Escobar (2013), CQG, 30(19), p.195020. University of Otago, Dunedin, New Zealand
More informationwith Matter and Radiation By: Michael Solway
Interactions of Dark Energy with Matter and Radiation By: Michael Solway Advisor: Professor Mike Berger What is Dark Energy? Dark energy is the energy needed to explain the observed accelerated expansion
More information3 The Friedmann-Robertson-Walker metric
3 The Friedmann-Robertson-Walker metric 3.1 Three dimensions The most general isotropic and homogeneous metric in three dimensions is similar to the two dimensional result of eq. (43): ( ) dr ds 2 = a
More informationLecture 05. Cosmology. Part I
Cosmology Part I What is Cosmology Cosmology is the study of the universe as a whole It asks the biggest questions in nature What is the content of the universe: Today? Long ago? In the far future? How
More informationConserved Quantities in Lemaître-Tolman-Bondi Cosmology
1/15 Section 1 Section 2 Section 3 Conserved Quantities in Lemaître-Tolman-Bondi Cosmology Alex Leithes - Blackboard Talk Outline ζ SMTP Evolution Equation: ζ SMTP = H X + 2H Y 3 ρ Valid on all scales.
More informationThe first one second of the early universe and physics beyond the Standard Model
The first one second of the early universe and physics beyond the Standard Model Koichi Hamaguchi (University of Tokyo) @ Colloquium at Yonsei University, November 9th, 2016. Credit: X-ray: NASA/CXC/CfA/M.Markevitch
More informationTachyonic dark matter
Tachyonic dark matter P.C.W. Davies Australian Centre for Astrobiology Macquarie University, New South Wales, Australia 2109 pdavies@els.mq.edu.au Abstract Recent attempts to explain the dark matter and
More informationScale symmetry a link from quantum gravity to cosmology
Scale symmetry a link from quantum gravity to cosmology scale symmetry fluctuations induce running couplings violation of scale symmetry well known in QCD or standard model Fixed Points Quantum scale symmetry
More informationChallenging the Cosmological Constant
Challenging the Cosmological Constant Based on: work to appear next week Overview Dark thoughts Where fields hide Environmental mass effects and chameleonic behavior Changeling A chameleon that actually
More informationCalculation of Momentum Distribution Function of a Non-Thermal Fermionic Dark Matter
Calculation of Momentum Distribution Function of a Non-Thermal Fermionic Dark Matter, March 8, 2017. arxiv:1612.02793, with Anirban Biswas. Aritra Gupta Why Non-Thermal? 1 / 31 The most widely studied
More informationUniversité Paris Diderot. Classification of inflationary models and constraints on fundamental physics. 28 September Ph.D.
i Université Paris Diderot Classification of inflationary models and constraints on fundamental physics. 28 September 2016 Ph.D. Thesis Candidate Mauro Pieroni mauropieroni@gmail.com Supervisors: Prof.
More informationNeutrino Mass Limits from Cosmology
Neutrino Physics and Beyond 2012 Shenzhen, September 24th, 2012 This review contains limits obtained in collaboration with: Emilio Ciuffoli, Hong Li and Xinmin Zhang Goal of the talk Cosmology provides
More informationCosmology. Assumptions in cosmology Olber s paradox Cosmology à la Newton Cosmology à la Einstein Cosmological constant Evolution of the Universe
Cosmology Assumptions in cosmology Olber s paradox Cosmology à la Newton Cosmology à la Einstein Cosmological constant Evolution of the Universe Assumptions in Cosmology Copernican principle: We do not
More informationA A + B. ra + A + 1. We now want to solve the Einstein equations in the following cases:
Lecture 29: Cosmology Cosmology Reading: Weinberg, Ch A metric tensor appropriate to infalling matter In general (see, eg, Weinberg, Ch ) we may write a spherically symmetric, time-dependent metric in
More informationA Curvature Primer. With Applications to Cosmology. Physics , General Relativity
With Applications to Cosmology Michael Dine Department of Physics University of California, Santa Cruz November/December, 2009 We have barely three lectures to cover about five chapters in your text. To
More informationThird Year: General Relativity and Cosmology. 1 Problem Sheet 1 - Newtonian Gravity and the Equivalence Principle
Third Year: General Relativity and Cosmology 2011/2012 Problem Sheets (Version 2) Prof. Pedro Ferreira: p.ferreira1@physics.ox.ac.uk 1 Problem Sheet 1 - Newtonian Gravity and the Equivalence Principle
More informationWeak interactions. Chapter 7
Chapter 7 Weak interactions As already discussed, weak interactions are responsible for many processes which involve the transformation of particles from one type to another. Weak interactions cause nuclear
More informationModel Universe Including Pressure
Model Universe Including Pressure The conservation of mass within the expanding shell was described by R 3 ( t ) ρ ( t ) = ρ 0 We now assume an Universe filled with a fluid (dust) of uniform density ρ,
More informationThermodynamics in Cosmology Nucleosynthesis
Thermodynamics in Cosmology Nucleosynthesis Thermodynamics Expansion Evolution of temperature Freeze out Nucleosynthesis Production of the light elements Potential barrier Primordial synthesis calculations
More informationFRW cosmology: an application of Einstein s equations to universe. 1. The metric of a FRW cosmology is given by (without proof)
FRW cosmology: an application of Einstein s equations to universe 1. The metric of a FRW cosmology is given by (without proof) [ ] dr = d(ct) R(t) 1 kr + r (dθ + sin θdφ ),. For generalized coordinates
More informationA glimpse on Cosmology: Mathematics meets the Data
Naples 09 Seminar A glimpse on Cosmology: Mathematics meets the Data by 10 November 2009 Monica Capone 1 Toward a unified epistemology of Sciences...As we know, There are known knowns. There are things
More informationParticle Cosmology. V.A. Rubakov. Institute for Nuclear Research of the Russian Academy of Sciences, Moscow and Moscow State University
Particle Cosmology V.A. Rubakov Institute for Nuclear Research of the Russian Academy of Sciences, Moscow and Moscow State University Topics Basics of Hot Big Bang cosmology Dark matter: WIMPs Axions Warm
More informationParticles in the Early Universe
Particles in the Early Universe David Morrissey Saturday Morning Physics, October 16, 2010 Using Little Stuff to Explain Big Stuff David Morrissey Saturday Morning Physics, October 16, 2010 Can we explain
More informationMaking Dark Matter. Manuel Drees. Bonn University & Bethe Center for Theoretical Physics. Making Dark Matter p. 1/35
Making Dark Matter Manuel Drees Bonn University & Bethe Center for Theoretical Physics Making Dark Matter p. 1/35 Contents 1 Introduction: The need for DM Making Dark Matter p. 2/35 Contents 1 Introduction:
More informationLecture 12. Inflation. What causes inflation. Horizon problem Flatness problem Monopole problem. Physical Cosmology 2011/2012
Lecture 1 Inflation Horizon problem Flatness problem Monopole problem What causes inflation Physical Cosmology 11/1 Inflation What is inflation good for? Inflation solves 1. horizon problem. flatness problem
More informationRelativity, Gravitation, and Cosmology
Relativity, Gravitation, and Cosmology A basic introduction TA-PEI CHENG University of Missouri St. Louis OXFORD UNIVERSITY PRESS Contents Parti RELATIVITY Metric Description of Spacetime 1 Introduction
More information