LOG#122. Basic Neutrinology(VII).

The observed mass and mixing both in the neutrino and quark cases could be evidence for some interfamily hierarchy hinting that the lepton and quark sectors were, indeed, a result of the existence of a new quantum number related to “family”. We could name this family symmetry as U(1)_F. It was speculated by people like Froggatt long ago. The actual intrafamily hierarchy, i.e., the fact that m_u>>m_d in the quark sector, seem to require one of these symmetries to be anomalous.

A simple model with one family dependent anomalous U(1) beyond the SM was first proposed long ago to produce the given Yukawa coupling and their hierarchies, and the anomalies could be canceled by the Green-Schwarz mechanism which as by-product is able to fix the Weinberg angle as well. Several developments include the models inspired by the E_6\times E_8 GUT or the E_8\times E_8 heterotic superstring theory. The gauge structure of the model is that of the SM but enlarged by 3 abelian U(1) symmetries and their respective fields, sometimes denoted by X,Y^{1,2}. The first one is anomalous and family independent. Two of these fields, the non-anomalous, have specific dependencies on the 3 chiral families designed to reproduce the Yukawa hierarchies. There are right-handed neutrinos which “trigger” neutrino masses by some special types of seesaw mechanisms.

The 3 symmetries and their fields X,Y^{1,2} are usually spontaneously broken at some high energy scale M_X by stringy effects. It is assumed that 3 fields, \theta_i, with i=1,2,3, develop a non-null vev. These \theta_i fields are singlets under the SM gauge group but not under the abelian symmetries carried by X, Y^{1,2}. Thus, the Yukawa couplings appear as some effective operators after the U(1)_F spontaneous symmetry breaking. In the case of neutrinos, we have the mass lagrangian (at effective level):

\mathcal{L}_m\sim h_{ij}L_iH_uN_j^c\lambda^{q_i+n_j}+M_N\xi_{ij}N_i^cN_j^c\lambda^{n_i+n_j}

and where h_ {ij},\xi_{ij}\sim \mathcal{O}(1). The parameters \lambda determine the mass and mixing hierarchy with the aid of some simple relationships:

\lambda=\dfrac{\langle \theta\rangle}{M_X}\sim\sin\theta_c

and where \theta_c is the Cabibblo angle. The q_i,n_i are the U(1)_F charges assigned to the left handed leptons L and the right handed neutrinos N. These couplings generate the following mass matrices for neutrinos:

m_\nu^D=\mbox{diag}(\lambda^{q_1},\lambda^{q_2},\lambda^{q_3})\hat{h}\mbox{diag}(\lambda^{n_1},\lambda^{n_2},\lambda^{n_3})\langle H_u\rangle


From these matrices, the associated seesaw mechanism gives the formula for light neutrinos:

m_\nu\approx \dfrac{\langle H_u\rangle^2}{M_X}\mbox{diag}(\lambda^{q_1},\lambda^{q_2},\lambda^{q_3})\hat{h}\hat{\xi}^{-1}\hat{h}^T\mbox{diag}(\lambda^{q_1},\lambda^{q_2},\lambda^{q_3})

The neutrino mass mixing matrix depends only on the charges we assign to the LH neutrinos due to cancelation of RH neutrino charges and the seesaw mechanism. There is freedom in the assignment of the charges q_i. If the charges of the second and the third generation of leptos are equal (i.e., if q_2=q_3), then one is lead to a mass matrix with the following structure (or “texture”):

m_\nu\sim \begin{pmatrix}\lambda^6 & \lambda^3 & \lambda^3\\ \lambda^3 & a & b\\ \lambda^3 & b & c\end{pmatrix}

and where a,b,c\sim \mathcal{O}(1). This matrix can be diagonalized in a straightforward fashion by a large \nu_2-\nu_3 rotation. It is consistent (more or less), with a large \mu-\tau mixing. In this theory or model, the explanation of the large neutrino mixing angles is reduced to a theory of prefactors in front of powers of the parameters \lambda, related with the vev after the family group spontaneous symmetry breaking!

LOG#106. Basic Cosmology (I).


The next thread is devoted to Cosmology. I will intend to be clear and simple about equations and principles of current Cosmology with a General Relativity background.

First of all…I will review the basic concepts of natural units I am going to use here. We will be using the following natural units:


We will take the Planck mass to be given by

M_P=\sqrt{8\pi G_N}\approx 1\mbox{.}2\cdot 10^{19}GeV

The solar mass is M_\odot=2\cdot 10^{30}kg and the parsec is given by the value

1pc=3\mbox{.}26lyr=3\mbox{.}1\cdot 10^{16}m

Well, current Cosmology is based on General Relativity. Even if I have not reviewed this theory with detail in this blog, the nice thing is that most of Cosmology can be learned with only a very little knowledge of this fenomenal theory. The most important ideas are: metric field, geodesics, Einstein equations and no much more…

In fact, newtonian gravity is a good approximation in some particular cases! And we do know that even in this pre-relativistic theory

\mbox{Gravitational force}=\mbox{Matter/Mass density}

via the Poisson’s equation

\nabla^2\phi =4\pi G_N\rho

This idea, due to the equivalence principle, is generalized a little bit in the general relativistic framework

\mbox{Spacetime geometry}=\mbox{Matter content/Energy-momentum}

The spacetime geometry is determined by the metric tensor g_{\mu\nu}(x). The matter content is given by the stress-energy-momentum tensor T_{\mu\nu}. As we know one of these two elements, we can know, via Eisntein’s field equations the another. That is, given a metric tensor, we can tell how energy-momentum “moves” in space-time. Given the energy-momentum tensor, we can know what is the metric tensor in spacetime and we can guess how the spacetime bends… This is the origin of the famous motto: “Spacetime says matter how to move, energy-momentum says spacetime how to curve”! Remember that we have “deduced” the Einstein’s field equations in the previous post. Without a cosmological constant term, we get

G_{\mu\nu}=\kappa^2T_{\mu\nu}=8\pi G_NT_{\mu\nu}

Given a spacetime metric g_{\mu\nu}, we can calculate the (affine/Levi-Civita) connection

\Gamma^\sigma_{\;\;\mu\nu}=\dfrac{1}{2}g^{\sigma\rho}\left(\partial_\mu g_{\nu\rho}+\partial_\nu g_{\rho\mu}-\partial_\rho g_{\mu\nu}\right)

The Riemann tensor that measures the spacetime curvature is provided by the equation

R^\rho_{\;\; \sigma\mu\nu}=\partial_\mu \Gamma^\rho_{\;\;\mu\sigma}-\partial_\mu \Gamma^\rho_{\;\; \mu \sigma}+\Gamma^\rho_{\;\;\mu\lambda}\Gamma^\lambda_{\;\;\nu\sigma}-\Gamma^\rho_{\;\;\nu\lambda}\Gamma^\lambda_{\;\;\mu\sigma}

The Ricci tensor is defined to be the following “trace” of the Riemann tensor

R_{\mu\nu}=R^\lambda_{\;\;\mu\lambda \nu}

The Einstein tensor is related to the above tensors in the well-known manner


The Einstein’s equations can be derived from the Einstein-Hilbert action we learned in the previous post, using the action principle and the integral

\boxed{S_{EH}=\int d^4x \sqrt{-g} \left(\kappa^{-2}R+\mathcal{L}_M\right)}

The geodesic equation is the path of a freely falling particle. It gives a “condensation” of the Einstein’s equivalence principle too and it is also a generalization of Newton’s law of “no force”. That is, the geodesic equation is the feynmanity

\boxed{\dfrac{d^2 x^\mu}{d\tau^2}+\Gamma^\mu _{\rho\sigma}\dfrac{dx^\rho}{d\tau}\dfrac{dx^\sigma}{d\tau}=0}

Finally, an important concept in General Relativity is that of isometry. The symmetry of the “spacetime manifold” is provided by a Killing vector that preserves transformations (isometries) of that manifold. Mathematically speaking, the Killing vector fields satisfy certain equation called the Killing equation

\boxed{\xi_{\mu ; \nu}+\xi_{\nu ; \mu}=0}

Maximally symmetric spaces have n(n+1)/2 Killing vectors in n-dimensional (nD) spacetime. There are 3 main classes or types of 2D maximally symmetric that can be generalized to higher dimensions:

1. The euclidean plane E^2.

2. The pseudo-sphere H^2. This is a certain “hyperbolic” space.

3. The spehre S^2. This is a certain “elliptic” space.

The Friedmann-Robertson-Walker Cosmology

Current cosmological models are based in General Relativity AND  a simplification of the possible metrics due to the so-called Copernican (or cosmological) principle: the Universe is pretty much the same “everywhere” you are in the whole Universe! Remarkbly, the old “perfect” Copernican (cosmological) principle that states that the Universe is the same “everywhere” and “every time” is wrong. Phenomenologically, we have found that the Universe has evolved and it evolves, so the Universe was “different” when it was “young”. Therefore, the perfect cosmological principle is flawed. In fact, this experimental fact allows us to neglect some old theories like the “stationary state” and many other “crazy theories”.

What are the observational facts to keep the Copernican principle? It seems that:

1st. The distribution of matter (mainly galaxies, clusters,…) and radiation (the cosmic microwave background/CMB) in the observable Universe is homogenous and isotropic.

2nd. The Universe is NOT static. From Hubble’s pioneer works/observations, we do know that galaxies are receeding from us!

Therefore, these observations imply that our “local” Hubble volume during the Hubble time is similar to some spacetime with homogenous and isotropic spatial sections, i.e., it is a spacetime manifold M=\mathbb{R}\times \Sigma. Here, \mathbb{R} denotes the time “slice” and \Sigma represents a 3D maximally symmetric space.

The geometry of a locally isotropic and homogeneous Universe is represented by the so-called Friedmann-Robertson-Walker metric


Here, a(t) is the called the scale factor.  The parameter k determines the geometry type (plane, hyperbolic or elliptical/spherical):

1) If k=0, then the Universe is “flat”. The manifold is E^3.

2) If k=-1, then the Universe is “open”/hyperbolic. The manifold would be H^3.

3) If k=+1, then the Universe is “closed”/spherical or elliptical. The manifold is then S^3.

Remark: The ansatz of local homogeneity and istoropy only implies that the spatial metric is locally one of the above three spaces, i.e., E^3,H^3,S^3. It could be possible that these 3 spaces had different global (likely topological) properties beyond these two properties.

Kinematical features of a FRW Universe

The first property we are interested in Cosmology/Astrophysics is “distance”. Measuring distance in a expanding Universe like a FRW metric is “tricky”. There are several notions of “useful” distances. They can be measured by different methods/approaches and they provide something called sometimes “the cosmologidal distance ladder”:

1st. Comoving distance. It is a measure in which the distance is “taken” by a fixed coordinate.

2nd. Physical distance. It is essentially the comoving distance times the scale factor.

3rd. Luminosity distance. It uses the light emitted by some object to calculate its distance (provided the speed of light is taken constant, i.e., special relativity holds and we have a constant speed of light)

4th. Angular diameter distance. Another measure of distance using the notion of parallax and the “extension” of the physical object we measure somehow.

There is an important (complementary) idea in FRW Cosmology: the particle horizon. Consider a light-like particle with ds^2=0. Then,




The total comoving distance that light have traveled since a time t=0 is equal to


It shows that NO information could have propagated further and thus, there is a “comoving horizon” with every light-like particle! Here, this time is generally used as a “conformal time” as a convenient tiem variable for the particle. The physical distance to the particle horizon can be calculated


There are some important kinematical equations to be known

A) For the geodesic equation, the free falling particle, we have




for the FRW metric and, moreover, the energy-momentum vector P^\mu=(E,\mathbf{p}) is defined by the usual invariant equation


This definition defines, in fact, the proper “time” \lambda implicitely, since


and the 0th component of the geodesic equation becomes


g_{\mu\nu}p^\mu p^\nu=-E^2+a^2\delta_{ij}p^ip^j=-m^2

EdE=a^2\vert \mathbf{p}\vert d\vert \mathbf{p}\vert

a^2 p\dfrac{dp}{dt}=-a\dot{a} p^2

\dfrac{1}{\vert \mathbf{p}\vert }\dfrac{d\vert \mathbf{p}\vert}{dt}+\dfrac{\dot{a}}{a}=0

Therefore we have deduced that \vert \mathbf{p}\vert \propto a^{-1}. This is, in fact, the socalled “redshift”.  The cosmological  redshift parameter is more generally defined through the equation


B) The Hubble’s law.

The luminosity distance measures the flux of light from a distant object of known luminosity (if it is not expanding). The flux and luminosity distance are bound into a single equation

\boxed{F=\dfrac{L}{4\pi d^2_L}}

If we use the comoving distance between a distant emitter and us, we get

\chi (a)=\int_t^{t_0}\dfrac{dt'}{a(t')}=\int_a^1\dfrac{da'}{a'^2 H(a')}

for a expanding Universe! That is, we have used the fact that luminosity itself goes through a comoving spherical shell of radius \chi (a). Moreover, it shows that

F=\dfrac{L (\chi)}{4 \pi \chi (a)^2 a_0^2}=\dfrac{L}{4\pi (\chi (a)/a)^2}

The luminosity distance in the expanding shell is

d_L=\dfrac{\chi (a)}{a}=\left(\dfrac{L}{4\pi F}\right)^{1/2}

and this is what we MEASURE in Astrophysics/Cosmology. Knowing a(t), we can express the luminosity distance in terms of the redshift. Taylor expansion provides something like this:


where higher order terms are sometimes referred as “statefinder parameters/variables”. In particular, we have




C) Angular diameter distance.

If we know that some object has a known length l, and it gives some angular “aperture” or separation \theta, the angular diameter distance is given by


The comoving size is defined as l/a, and the coming distance is again \chi (a). For “flat” space, we obtain that

\theta=\dfrac{l/a}{\chi (a)}

that is

d_A=a\chi (a)=\dfrac{\chi}{1+z}

In the case of “curved” spaces, we get

d_A=\dfrac{a}{H_0\sqrt{\vert \omega_k\vert}}\cdot\begin{cases}\sinh \left( \sqrt{\Omega_k}H_0\chi\right),\Omega_k>0\\ \sin \left( \sqrt{-\Omega_k}H_0\chi\right),\Omega_k<0\end{cases}

FRW dynamics

Gravity in General Relativity, a misnomer for the (locally) relativistic theory of gravitation, is described by a metric field, i.e., by a second range tensor (covariant tensor if we are purist with the nature of components). The metric field is related to the matter-energy-momentum content through the Einstein’s equations

G_{\mu\nu}=-\kappa^2 T_{\mu\ nu}

The left-handed side can be calculated for a FRW Universe as follows




The right-handed side is the energy-momentum of the Universe. In order to be fully consistent with the symmetries of the metric, the energy-momentum tensor MUST be diagonal and T_{11}=T_{22}=T_{33}=T. In fact, this type of tensor describes a perfect fluid with

T_{\mu\nu}=(\rho+p)U_\mu U_\nu+pg_{\mu\nu}

Here, \rho, p are functions of t (cosmological time) only. They are “state variables” somehow. Moreover, we have

U_\mu =(1,0,0,0)

for the fluid at rest in the comoving frame. The Friedmann equations are indeed the EFE for a FRW metric Universe

3\left(\dfrac{\dot{a}^2}{a^2}+\dfrac{k}{a^2}\right)=\kappa^2\rho for the 00th compoent as “constraint equation.

2\dfrac{\ddot{a}}{a}+\dfrac{\dot{a}^2}{a^2}+\dfrac{k}{a^2}=-\kappa^2p for the iith components.

Moreover, we also have


and this conservation law implies that


Therefore, we have got two independent equations for three unknowns (a, \rho, p). We need an additional equation. In fact, the equation of state for p=p(\rho) provides such an additional equation. It gives the “dynamics of matter”!

In summary, the basic equations for Cosmology in a FRW metric, via EFE, are the Friedmann’s equations (they are secretly the EFE for the FRW metric) supplemented with the energy-momentum conservations law and the equation of state for the pressure p=p(\rho):

1) \boxed{\dfrac{\dot{a}^2}{a^2}+\dfrac{k^2}{a^2}=\dfrac{\kappa^2}{3}\rho}

2) \boxed{\dot{\rho}+3\dfrac{\dot{a}}{a}(\rho+p)=0}

3) \boxed{p=p(\rho)}

There are many kinds of “matter-energy” content of our interest in Cosmology. Some of them can be described by a simple equation of state:

\boxed{p=\omega \rho}

Energy-momentum conservation implies that \rho\propto a^{-3(\omega +1)}. 3 special cases are used often:

1st. Radiation (relativistic “matter”). \omega=1/3 and thus, p=1/3\rho and \rho\propto a^{-4}

2nd. Dust (non-relativistic matter). \omega=0. Then, p=0 and \rho\propto a^{-3}

3rd. Vacuum energy (cosmological constant). \omega=-1. Then, p=-\rho and \rho=\mbox{constant}

Remark (I): Particle physics enters Cosmology here! Matter dynamics or matter fields ARE the matter content of the Universe.

Remark (II): Existence of a Big Bang (and a spacetime singularity). Using the Friedmann’s equation


if we have that (\rho+3p)>0, the so-called weak energy condition, then a=0 should have been reached at some finite time in the past! That is the “Big Bang” and EFE are “singular” there. There is no scape in the framework of GR. Thus, we need a quantum theory of gravity to solve this problem OR give up the FRW metric at the very early Universe by some other type of metric or structure.

Particles and the chemical equilibrium of the early Universe

Today, we have DIRECT evidence for the existence of a “thermal” equilibrium in the early Universe: the cosmic microwave background (CMB). The CMB is an isotropic, accurate and non-homogeneous (over certain scales) blackbody spectrum about T\approx 3K!

Then, we know that the early Universe was filled with a hot dieal gas in thermal equilibrium (a temperature T_e can be defined there) such as the energy density and pressure can be written in terms of this temperature. This temperature generates a distribution f(\mathbf{x},\mathbf{p}). The number of phase space elements in d^3xd^3p is


and where the RHS is due to the uncertainty principle. Using homogeneity, we get that, indeed, f(x,p)=f(p), and where we can write the volume d^3x=dV. The energy density and the pressure are given by (natural units are used)

\rho_i=g_i\int \dfrac{d^3p}{(2\pi)^3}f_i(p)E(p)

p_i=g_i\int \dfrac{d^3p}{(2\pi)^3}f_i (p)\dfrac{p^2}{3E(p)}

When we are in the thermal equilibrium at temperature T, we have the Bose-Einstein/Fermi-Dirac distribution

f(p)=\dfrac{1}{e^{(E-\mu)/T}\pm 1}

and where the + is for the Fermi-Dirac distribution (particles) and the - is for the Bose-Einstein distribution (particles). The number density, the energy density and the pressure are the following integrals

\boxed{\mbox{Number density}=n=\dfrac{N}{V}=\dfrac{g}{2\pi^2}\int_m^\infty \dfrac{(E^2-m^2)^{1/2}}{e^{(E-\mu)/T}\pm 1}dE}

\boxed{\mbox{Density energy}=\rho=\dfrac{E}{V}=\dfrac{g}{2\pi^2}\int_m^\infty \dfrac{(E^2-m^2)^{1/2}E^2}{e^{(E-\mu)/T}\pm 1}dE}

\boxed{\mbox{Pressure}=p=\dfrac{g}{6\pi^2}\int_m^\infty \dfrac{(E^2-m^2)^{3/2}}{e^{(E-\mu)/T}\pm 1}dE}

And now, we find some special cases of matter-energy for the above variables:

1st. Relativistic, non-degenerate matter (e.g. the known neutrino species). It means that T>>m and T>>\mu. Thus,

n=\left(\dfrac{3}{4}\right)\dfrac{\zeta (3)}{\pi^2}gT^3



2nd. Non-relativistic matter with m>>T only. Then,


\rho= mn+\dfrac{3}{2}p, and p=nT<<\rho

The total energy density is a very important quantity. In the thermal equilibrium, the energy density of non-relativistic species is exponentially smaller (suppressed) than that of the relativistic particles! In fact,

\rho_R=\dfrac{\pi^2}{30}g_\star T^4 for radiation with p_R=\dfrac{1}{3}\rho_R

and the effective degrees of freedom are


Remark: The factor 7/8 in the DOF and the variables above is due to the relation between the Bose-Einstein and the Fermi-Dirac integral in d=3 space dimensions. In general d, the factor would be


Entropy conservation and the early Universe

The entropy in a comoving volume IS a conserved quantity IN THE THERMAL EQUILIBRIUM. Therefore, we have that

\dfrac{\partial p_i}{\partial T}=g_i\int \dfrac{d^3p}{(2\pi)^3}\dfrac{df}{dT}\dfrac{p^2}{3E(p)}=g_i\int \dfrac{4\pi pE dE}{(2\pi)^3}\dfrac{df}{dE}\left(-\dfrac{E}{T}\right)\dfrac{p^2}{3E}

and then

\dfrac{\partial p_i}{\partial T}=\dfrac{g_i}{2\pi^2}\int \left(-\dfrac{d}{dE}\left(f\dfrac{p^3E}{3T}\right)+f\dfrac{d}{dE}\left(\dfrac{p^3E}{3T}\right)\right)dE


\dfrac{\partial p_i}{\partial T}=\dfrac{1}{T}(\rho_i+p_i)

Now, since

\dfrac{\partial \rho}{\partial t}+3\dfrac{\dot{a}}{a}(\rho+p)=0


\dfrac{\partial}{\partial t}\left(a^3(\rho+p)\right)-a^3\dfrac{\partial p}{\partial t}=0

\dfrac{1}{a^3}\dfrac{\partial (a^3(\rho +p))}{\partial t}-\dfrac{\partial \rho}{\partial t}=0

if we multiply by T and use the chain rule for \rho, we obtain

\dfrac{1}{a^3}\dfrac{\partial}{\partial t}\left(\dfrac{a^3(\rho+p)}{T}\right)=0

but it means that a^3s=\mbox{constant}, where s is the entropy density defined by

\boxed{s\equiv \dfrac{\rho+p}{T}}

Well, the fact is that we know that the entropy or more precisely the entropy density is the early Universe is dominated by relativistic particles ( this is “common knowledge” in the Stantard Cosmological Model, also called \Lambda CDM). Thus,

\boxed{s=\dfrac{2\pi^2}{45}g_\star T^3}

It implies the evolution of temperature with the redshift in the following way:

T\propto g_\star^{-1/3}a^{-1}

Indeed, since we have that n\propto a^{-3}, s\propto a^{-3}, the yield variable

Y_i\equiv \dfrac{n_i}{s}

is a convenient quantity that represents the “abundance” of decoupled particles.

See you in my next cosmological post!