Neutrino 0.1.4

1.3.4. Halo Effect - Analytical Methods

1.3.4.1. Linear Stability

1.3.4.1.1. Single Beam Both Start with Electron Flavor

We assume that both the forward and backward beams are starting from electron flavor.

\[\begin{split}\rho_F &= \frac{1}{2}\begin{pmatrix} 1 & \epsilon_F \\ \epsilon_F^* & -1 \end{pmatrix} \\ \rho_B &= \frac{1}{2}\begin{pmatrix} 1 & \epsilon_B \\ \epsilon_B^* & -1 \end{pmatrix}.\end{split}\]

The linearized equation of motion is

\[\begin{split}i\partial_z \begin{pmatrix} \epsilon_F \\ \epsilon_B \end{pmatrix} = \begin{pmatrix} R \mu + \eta & - R \mu \\ \mu & -(\mu+ \eta) \end{pmatrix} \begin{pmatrix} \epsilon_F \\ \epsilon_B \end{pmatrix},\end{split}\]

where \(\eta\) is the hierarchy and all quantities are scaled by \(\omega_v\).

We can solve the eigenvalues of this problem. To obtain complex solutions, we require \(\Delta < 0\), which translates to

\[\mu_-< \mu <\mu_+,\]

where

\[\mu_\pm = \frac{ -2 (1+R) \eta \pm 8 \sqrt{ R } }{ (1-R)^2 }.\]
../../_images/table-of-mu-values-for-instability.jpg

Fig. 1.182 Table of {\(\mu_{-}\), \(\mu_{+}\)}.

We plot out growth rate as a function of \(R\) and \(\mu\), which shows that such so called instabilities only exist for normal hierarchy.

../../_images/growth-rate-mu-refl-nh.jpg

Fig. 1.183 For normal hierarchy. Vacuum mixing angle is set to 0.

../../_images/growth-rate-mu-refl-ih.jpg

Fig. 1.184 For inverted hierarchy. Vacuum mixing angle is set to 0.

1.3.4.1.2. Single Beam Undetermined Back Beam

The equation of motion is

\[\begin{split}i\partial_z \rho_F &= [ H_F, \rho_F ] \\ i\partial_z \rho_B &= [ H_B, \rho_B ] ,\end{split}\]

where

\[\begin{split}H_F &= R \mu \rho_B\\ H_B &= \mu \rho_F.\end{split}\]

Conventions

The angle contribution is \(1-\cos\theta=2\). We absorbe this factor into \(\mu\).

As for the initial states, we can assume some general forms

\[\begin{split}\rho_F &= \frac{1}{2} \begin{pmatrix} 1 & \epsilon_F \\ \epsilon_F^* & -1 \end{pmatrix}\\ \rho_B &= \frac{1}{2} \begin{pmatrix} a & \epsilon_B \\ \epsilon_B^* & -a \end{pmatrix}.\end{split}\]

The forward beam start from neutrino flavors. However, the state of the backward beam is to be determined. The exact state is not necessary since we are interested in linear stability.

We linearize the equation of motion by keeping only first order of the small off-diagonal elements.

\[\begin{split}i\partial_z \begin{pmatrix} 1 & \epsilon_F \\ \epsilon_F^* & -1 \end{pmatrix} &= \frac{ R \mu}{2} \begin{pmatrix} 0 & 2(a \epsilon_F -\epsilon_B) \\ c.c. & 0 \end{pmatrix} + \frac{R\mu}{2} \begin{pmatrix} 0 & 2(a\epsilon_F - b) \\ c.c. & 0 \end{pmatrix} \\ - i\partial_z \begin{pmatrix} 1 & \epsilon_B \\ \epsilon_B^* & -1 \end{pmatrix} &= -\frac{\mu}{2} \begin{pmatrix} 0 & 2(a \epsilon_F -\epsilon_B) \\ c.c. & 0 \end{pmatrix} - \frac{\mu}{2} \begin{pmatrix} 0 & 2 (a\epsilon_F - b) \\ c.c. & 0 \end{pmatrix}\end{split}\]

The system of equations can be easily solved and we found that the exponentials are \(e^{i(2a R -1) \mu z}\). No exponent growth is possible.

eqnF = I D[epsF[z], z] == alpha mu (a epsF[z] - epsB[z]) + alpha mu (a epsF[z] - b)
eqnB = I D[epsB[z], z] == mu (a epsF[z] - epsB[z]) + mu (a epsF[z] - b)
DSolve[{eqnF,eqnB},{epsF,epsB},z]

The Two Beams have Synchronized Evolution

Notice that the two equations can be combined to

\[\partial_z \epsilon_F = \partial_z ( R \epsilon_B ).\]

This result indicates that the initial development of the forward and backward beams are synchronized.

Dealing with Inhomogenous Differential Equations

  1. The solution to inhomogenous systems is the general solution plus the particular solution. At least we can know that if the general solution contains exponential growth terms, the final solution should also contain exponential growth terms, if the particular solution doesn’t cancel them.

  2. Floquet theory for periodic matrix systems might be useful.

1.3.4.1.2.1. Linear Regime Behavior

The solutions to such a problem shows that we always obtain decreasing modes. It’s not easy to comphrend. But the linear could be solved exactly.

In linear regime, we define the density matrices for forward and backward beams to be

\[\begin{split}\rho_F &= \frac{1}{2} \begin{pmatrix} 1 & \epsilon_F \\ \epsilon_B^* & -1 \end{pmatrix} \\ \rho_B &= \frac{1}{2} \begin{pmatrix} 1 & \epsilon_B \\ \epsilon_B^* & -1 \end{pmatrix}.\end{split}\]

The Hamiltonians are

\[\begin{split}H_F &= H_v + \mu \rho_F \\ H_B &= H_v + R \mu \rho_B.\end{split}\]

We will investigate the instability for zero mixing angle.

The linearized equation of motion can be simplified to

\[\begin{split}i\partial_z \begin{pmatrix} \epsilon_F \\ \epsilon_B \end{pmatrix} = \begin{pmatrix} -\omega_v + R \xi \mu & - R \xi \mu \\ \xi \mu & \omega_v - \xi \mu \end{pmatrix} \begin{pmatrix} \epsilon_F \\ \epsilon_B \end{pmatrix}.\end{split}\]

This equation can be easily solved. The eigenvalues are

\[\begin{split}\Omega_+ &= \frac{1}{2} ( (R-1)\xi\mu + \sqrt{\Delta} ) \\ \Omega_- &= \frac{1}{2} ( (R-1)\xi\mu - \sqrt{\Delta} ),\end{split}\]

where

\[\Delta = (1-R)^2 \mu^2 \xi^2 - 4\mu\xi \omega_v (1+R) + 4\omega_v^2.\]

The corresponding eigenvectors are

\[\begin{split}V_+ &=\begin{pmatrix} \frac{ -2\omega_v + \xi \mu (1+R) + \sqrt{\Delta} }{2\xi\mu} \\ 1 \end{pmatrix} \\ V_- &=\begin{pmatrix} \frac{ -2\omega_v + \xi \mu (1+R) - \sqrt{\Delta} }{2\xi\mu} \\ 1 \end{pmatrix}.\end{split}\]

The general solution to the equation is

\[\begin{split}\begin{pmatrix} \epsilon_F(z) \\ \epsilon_B(z) \end{pmatrix} = C_+ V_+ e^{-i \Omega_+ z} + C_- V_- e^{-i \Omega_- z}.\end{split}\]

The special property about this reflection prolem is that the density matrices for the forward and backward beams should be the same at the reflection point, say \(L\). With such a simple relation, we can find the relations between \(C_\pm\) by setting \(\epsilon_F(L)=\epsilon_B(L)\),

\[\frac{C_+}{C_-} = e^{-i(\Omega_- -\Omega_+)L} \frac{ \sqrt{\Delta} + 2\omega_v + \mu \xi (1-R) }{\sqrt{\Delta} - 2\omega_v - \mu \xi (1-R)}.\]

The solution to be problem can be simplified,

\[\begin{split}\begin{pmatrix} \epsilon_F(z) \\ \epsilon_B(z) \end{pmatrix} = C_- e^{-i\Omega_- L} \left( \frac{ \sqrt{\Delta} + 2\omega_v + \mu \xi (1-R) }{\sqrt{\Delta} - 2\omega_v - \mu \xi (1-R)} V_+ e^{-i \Omega_+ (z-L)} + V_- e^{-i \Omega_- (z-L)} \right).\end{split}\]

We are interested in the absolute values of each elements so that the overall factors can be neglected. The forward beam evolution is obtained by taking the absolute value of \(\epsilon_F\),

\[\left\vert \epsilon_F \right\vert \propto \lvert (2\omega_v +\xi\mu(1-R) +i \delta ) ( -2\omega_v + \xi\mu(1+R) + i \delta ) e^{delta(z-L)} + ( -2\omega_v - \xi\mu(1-R) +i \delta ) ( -2\omega_v + \xi\mu(1+R) - i \delta ) e^{-delta(z-L)} \rvert,\]

in which \(\sqrt{\Delta}\) is replaced by \(i \delta\).

Collecting terms, we could verify that it has the form

\[\left\vert \epsilon_F \right\vert \propto A + B \cosh( 2\delta(L-z) ).\]

The only \(z\) dependent term is cosh( 2delta(L-z) ), which is decreasing within \([0,L]\) and is increasing in \([L,2L]\). The slope at \(z=L\) is 0.

../../_images/cosh-function.png

Fig. 1.185 \(\cosh(2\delta(z-L))\) with \(\delta=1\).

1.3.4.2. Nonlocal Boundary Value Problem

Some References

  1. Springer Open has a curated list of NBVP. DENFC.

  2. To solve coupled systems with nonlocal boundary conditions: Goodrich, C. S. (2017). Pointwise conditions in discrete boundary value problems with nonlocal boundary conditions. Applied Mathematics Letters, 67, 7–15.

1.3.4.3. References and Notes

Back to top

© 2021, Lei Ma | Created with Sphinx and . | On GitHub | Physics Notebook Statistical Mechanics Notebook Index | Page Source