How to understand Dispersion Relation

We should come back to the wave picture. DR tells us how easily a wave will disperse as it travels. The things that can be research on are the disperse itself. Stability is one important topic but it also tells us how the wave disperse.

Box Spectrum Defined Spectra Code

```
spectWC1 = {{{-1, -0.3}, 1}, {{-0.3, 1}, 1}};
spectWC2 = {{{-1, -0.3}, 0.5}, {{-0.3, 1}, 1}};
spectWC3 = {{{-1, -0.3}, 0.1}, {{-0.3, 1}, 1}};
spectWC4 = {{{-1, -0.3}, 0}, {{-0.3, 1}, 1}};
spectC1 = {{{-1, -0.3}, -0.1}, {{-0.3, 1}, 1}};
spectC2 = {{{-1, -0.3}, -0.5}, {{-0.3, 1}, 1}};
spectC3 = {{{-1, -0.3}, -1}, {{-0.3, 1}, 1}};
spectC4 = {{{-1, -0.3}, -1.5}, {{-0.3, 1}, 1}};
spectC5 = {{{-1, -0.3}, -15}, {{-0.3, 1}, 1}};
spectC6 = {{{-1, -0.3}, -0.01}, {{-0.3, 1}, 1}};
spectC7 = {{{-1, -0.3}, 0.01}, {{-0.3, 1}, 1}};
spectWC5 = {{{0.1, 0.3}, 1}, {{0.3, 1}, 1}};
```

Numerical Issues with Box Spectrum in LSA

Mathematica seems to fail the mission when finding the complex k’s given real \(\omega\).

I tested bot NIntegrate of the integrals and Integrate first. It seems that Integrate first works better and faster.

Since I could find the analytical expressions for the integrals, I simple plug in the expression then calculate numbers.

I could plot the values of \(\omega - (I_0-I_4)/4\) (MAA solution) for real \(\omega\) to check if they have interceptions with the 0 plane.

I also ploted out the MZA solutions. I found that interceptions exsits for only some values of \(\omega\). Examples of such plots are shown below.

Similarly for MZA- solution we have the following.

After some test, my thought is that the integral is nicely behaved for some range of real \(\omega\), which makes sense since the line of instability can not extend to arbitary large and small \(\omega\).

Thus the solution to the warnings and errors in FindRoot is to plugin the suitable range of \(\omega\).

Meanwhile for real \(k\), we have only solutions at \(\omega =0\).

MAA solutions also has similar shape.

In any case, we managed to solve instabilities for MAA and MZA. In Mathematica package dispersion-relation.wl, ConAxialSymOmegaNMZApEqnLHSComplex[-omegareal, k, spectrum] returns the value of \(\omega - \left( -\frac{1}{4}\left( I_0 - I_2 + \sqrt{ (I_0+I_2-2I_1)(I_0+I_2+2I_1) } \right) \right)\). Simply using FindRoot from Mathematica should return the solutions.

Since box spectrum is easier and faster to calculate, we’ll explore the phenomena within box spectrum. For example, the spectrum notation

```
{{{-1, -0.3}, -0.2}, {{-0.3, 1}, 0.7}}
```

indicates tht the spectrum has a constant value -0.2 within \(\cos\theta \in [-1,-0.3]\) and value 0.7 within \(\cos\theta \in [-0.3,1]\).

Observations and Questions

Any 0 in spectrum make it possible go cross the singularity line.

MZA disappears at some point of parameters. What are the constraints.

MZA disappears

The analytical expression for it to dispear shows that it involves all the elements of the spectrum. Numerically, I can show that for a specific one \(\{\{\{-1, -0.3\}, g_1\}, \{\{-0.3, 1\}, 1\}\}\), the region of \(g_1\) that leads to no MZA solution for a given \(n\) is shown in Fig. 1.141.

This figure shows that the MZA solution will come back if we choose really small \(g_1\). It is verified using another spectrum which is labeled as C5 spectrum.

However examples of numerical calculations won’t be enough. Alternatively, I could derive an analytical expression. MZA solution is

\[\omega = -\frac{1}{4}\left(I_0-I_2\pm \sqrt{ (I_0-2I_1+I_2)(I_0+2I_1+I_2) }\right).\]

We examine the term inside square root and set it to be negative which will show us the region where real omega disppears.

Comment: the upper limit is not a horizontal line.

For the MZA solutions to disappear, we need to find the smallest value of the upper limit (-0.73160 in this case) and the largest value of lower limit (-7.16327 in this case).

Taylor Expansion of MZA solutions

What about Taylor expansion around small values of \(g_1\)? Taylor expanding the most general abstract form doesn’t give us any useful results since it is super complicated.

I am not sure how to use these results.

Crossing has an effect on \(d\omega/dn\)

By observing the \(\omega\sim n\) plot, we noticed that crossing is changes the behavior of it, especially at the point \(n=-1\). It seems that crossing changes whether the \(\omega\) goes to \(\infty\) or \(-\infty\) at \(n=-1\).

More specifically, we plot out the \(d\omega/dn\) for three different \(g_1\).

MZA solutions Disapears?

The DR plots show that MZA solutions disappears for `spectC3`

and `spectC4`

. Is this true?

I tested `Exclusions->None`

in Mathematica to make sure this is not a plotting issue. The results do not change with this option.

However, I found that the \(\lvert f(\omega=1,k)\rvert\) is approaching 0 for some none zero real k.

This makes the DR plots questionable. However, no bug has been found. In fact I have been using the same functions to make this density plot.

The I realized that using only real \(n\)’s could be a problem since we might have complex \(n\)’s that give us real \(\omega\), in principle.

Indeed this is true for `spectC3`

.

Meanwhile the real parts of these solutions do not go to 0.

But there can never be real k for these real \(\omega\) lines because real \(k=\omega n\) where \(n\) is complex.

For `spectC3`

, we have the following.

So things must be special for `spectWC1`

etc. Here we plot the real and imaginary part of `spectWC3`

.

But there can never be real k for these real \(\omega\) lines because real \(k=\omega n\) where \(n\) is complex.

So, where did the MZA solutions go?

Crossing seems to have little effects on instabilities

Instabilities on DR plot seems to be NOT affected by crossing. Probably because of the lines in the forbidden region (using abs for log argument in the results of integral for I’s) doesn’t seem to change a lot.

What does backward emission mean? Is it somewhat equivalent to crossing spectrum without backward emission?

Analytically, crossing and backward emission are not simply related.

The MZA solution has some interesting behaviors. It seems that a crossing in spectrum can bend the DR and make it cross 0.

We can prove that the crossing indeed needed but not guarantee to make DR cross 0.

The MZA solution is

\[(-4\omega -(I_0 - I_2) )^2 = (I_0 -2 I_1 + I_2)(I_0 + 2 I_1 + I_2).\]

For \(\omega\to 0\), we have

\[I_0I_2 = I_1^2.\]

We restore the integral form, which is

\[\int du \frac{ G(u) u^2 }{ 1- n u} = \left( \int du \frac{ G(u) u }{ 1- n u} \right)^2.\]

We actually have \(1-n u>0\). Thus

\[\tilde G (u) = \frac{G(u)}{1-n u},\]

always preserves the sign of the spectrum. If the spectrum is always positive, then we are calculating some average of \(u\) and \(u^2\), the equation becomes

\[\langle u^2 \rangle = \left(\langle u \rangle \right)^2,\]

which can not be true. However, a crossing brings some negative regions and makes it possible but not guaranteed to have solutions.

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