summaryrefslogtreecommitdiff
path: root/source/know/concept/modulational-instability
diff options
context:
space:
mode:
authorPrefetch2022-10-20 18:25:31 +0200
committerPrefetch2022-10-20 18:25:31 +0200
commit16555851b6514a736c5c9d8e73de7da7fc9b6288 (patch)
tree76b8bfd30f8941d0d85365990bcdbc5d0643cabc /source/know/concept/modulational-instability
parente5b9bce79b68a68ddd2e51daa16d2fea73b84fdb (diff)
Migrate from 'jekyll-katex' to 'kramdown-math-sskatex'
Diffstat (limited to 'source/know/concept/modulational-instability')
-rw-r--r--source/know/concept/modulational-instability/index.md52
1 files changed, 26 insertions, 26 deletions
diff --git a/source/know/concept/modulational-instability/index.md b/source/know/concept/modulational-instability/index.md
index 628cbb4..a01293c 100644
--- a/source/know/concept/modulational-instability/index.md
+++ b/source/know/concept/modulational-instability/index.md
@@ -15,17 +15,17 @@ In fiber optics, **modulational instability** (MI)
is a nonlinear effect that leads to the exponential amplification
of background noise in certain frequency regions.
It only occurs in the [anomalous dispersion regime](/know/concept/dispersive-broadening/)
-($\beta_2 < 0$), which we will prove shortly.
+($$\beta_2 < 0$$), which we will prove shortly.
Consider the following simple solution to the nonlinear Schrödinger equation:
-a time-invariant constant power $P_0$ at the carrier frequency $\omega_0$,
+a time-invariant constant power $$P_0$$ at the carrier frequency $$\omega_0$$,
which is experiencing [self-phase modulation](/know/concept/self-phase-modulation/):
$$\begin{aligned}
A(z,t) = \sqrt{P_0} \exp( i \gamma P_0 z)
\end{aligned}$$
-We add a small perturbation $\varepsilon(z,t)$ to this signal,
+We add a small perturbation $$\varepsilon(z,t)$$ to this signal,
representing background noise:
$$\begin{aligned}
@@ -33,8 +33,8 @@ $$\begin{aligned}
\end{aligned}$$
We insert this into the nonlinear Schrödinger equation to get a perturbation equation,
-which we linearize by assuming that $|\varepsilon|^2$ is negligible compared to $P_0$,
-such that all higher-order terms of $\varepsilon$ can be dropped, yielding:
+which we linearize by assuming that $$|\varepsilon|^2$$ is negligible compared to $$P_0$$,
+such that all higher-order terms of $$\varepsilon$$ can be dropped, yielding:
$$\begin{aligned}
0
@@ -51,9 +51,9 @@ $$\begin{aligned}
\end{aligned}$$
We split the perturbation into real and imaginary parts
-$\varepsilon(z,t) = \varepsilon_r(z,t) + i \varepsilon_i(z,t)$,
+$$\varepsilon(z,t) = \varepsilon_r(z,t) + i \varepsilon_i(z,t)$$,
which we fill in in this equation.
-The point is that $\varepsilon_r$ and $\varepsilon_i$ are real functions:
+The point is that $$\varepsilon_r$$ and $$\varepsilon_i$$ are real functions:
$$\begin{aligned}
0
@@ -63,7 +63,7 @@ $$\begin{aligned}
\end{aligned}$$
Splitting this into its real and imaginary parts gives two PDEs
-relating $\varepsilon_r$ and $\varepsilon_i$:
+relating $$\varepsilon_r$$ and $$\varepsilon_i$$:
$$\begin{aligned}
\pdv{\varepsilon_r}{z} = \frac{\beta_2}{2} \pdvn{2}{\varepsilon_i}{t}
@@ -72,8 +72,8 @@ $$\begin{aligned}
\end{aligned}$$
We [Fourier transform](/know/concept/fourier-transform/)
-these in $t$ to turn them into ODEs relating
-$\tilde{\varepsilon}_r(z,\omega)$ and $\tilde{\varepsilon}_i(z,\omega)$:
+these in $$t$$ to turn them into ODEs relating
+$$\tilde{\varepsilon}_r(z,\omega)$$ and $$\tilde{\varepsilon}_i(z,\omega)$$:
$$\begin{aligned}
\pdv{\tilde{\varepsilon}_r}{z} = - \frac{\beta_2}{2} \omega^2 \tilde{\varepsilon}_i
@@ -82,7 +82,7 @@ $$\begin{aligned}
\end{aligned}$$
We are interested in exponential growth, so let us make the following ansatz,
-where $k$ may be a function of $\omega$, as long as it is $z$-invariant:
+where $$k$$ may be a function of $$\omega$$, as long as it is $$z$$-invariant:
$$\begin{aligned}
\tilde{\varepsilon}_r(z, \omega) = \tilde{\varepsilon}_r(0, \omega) \exp(k z)
@@ -91,7 +91,7 @@ $$\begin{aligned}
\end{aligned}$$
With this, we can write the system of ODEs for
-$\tilde{\varepsilon}_r(z,\omega)$ and $\tilde{\varepsilon}_i(z,\omega)$
+$$\tilde{\varepsilon}_r(z,\omega)$$ and $$\tilde{\varepsilon}_i(z,\omega)$$
in matrix form:
$$\begin{aligned}
@@ -112,7 +112,7 @@ $$\begin{aligned}
k = \pm \sqrt{ - \frac{\beta_2}{2} \omega^2 \Big( \frac{\beta_2}{2} \omega^2 + 2 \gamma P_0 \Big) }
\end{aligned}$$
-To get exponential growth, it is essential that $\mathrm{Re}\{k\} > 0$,
+To get exponential growth, it is essential that $$\mathrm{Re}\{k\} > 0$$,
so we discard the negative sign,
and get the following condition for MI:
@@ -124,22 +124,22 @@ $$\begin{aligned}
}
\end{aligned}$$
-Since $\omega^2$ is positive, $\beta_2$ must be negative,
+Since $$\omega^2$$ is positive, $$\beta_2$$ must be negative,
so MI can only occur in the ADR.
-It is worth noting that $\beta_2 = \beta_2(\omega_0)$,
+It is worth noting that $$\beta_2 = \beta_2(\omega_0)$$,
meaning there can only be exponential
noise growth when the "parent pulse" is in the anomalous dispersion regime,
but that growth may appear in areas of normal dispersion,
as long as the above condition is satisfied by the parent.
This result has been derived using perturbation,
-so only holds as long as $|\varepsilon|^2 \ll P_0$.
+so only holds as long as $$|\varepsilon|^2 \ll P_0$$.
Over time, the noise gets amplified so greatly
that this approximation breaks down.
-Next, we define the **gain** $g(\omega)$,
+Next, we define the **gain** $$g(\omega)$$,
which expresses how quickly the
-perturbation grows as a function of the frequency offset $\omega$:
+perturbation grows as a function of the frequency offset $$\omega$$:
$$\begin{aligned}
\boxed{
@@ -149,7 +149,7 @@ $$\begin{aligned}
}
\end{aligned}$$
-The frequencies with maximum gain are then found as extrema of $g(\omega)$,
+The frequencies with maximum gain are then found as extrema of $$g(\omega)$$,
which satisfy:
$$\begin{aligned}
@@ -162,10 +162,10 @@ $$\begin{aligned}
A simulation of MI is illustrated below.
The pulse considered was a solition of the following form
-with settings $T_0 = 10\:\mathrm{ps}$, $P_0 = 10\:\mathrm{kW}$,
-$\beta = -10\:\mathrm{ps}^2/\mathrm{m}$ and $\gamma = 0.1/\mathrm{W}/\mathrm{m}$,
+with settings $$T_0 = 10\:\mathrm{ps}$$, $$P_0 = 10\:\mathrm{kW}$$,
+$$\beta = -10\:\mathrm{ps}^2/\mathrm{m}$$ and $$\gamma = 0.1/\mathrm{W}/\mathrm{m}$$,
whose peak is approximately flat, so our derivation is valid there,
-hence it "wrinkles" in the $t$-domain:
+hence it "wrinkles" in the $$t$$-domain:
$$\begin{aligned}
A(0, t)
@@ -176,19 +176,19 @@ $$\begin{aligned}
<img src="pheno-mi-small.jpg" style="width:100%">
</a>
-Where $L_\mathrm{NL} = 1/(\gamma P_0)$ is the characteristic length of nonlinear effects.
+Where $$L_\mathrm{NL} = 1/(\gamma P_0)$$ is the characteristic length of nonlinear effects.
Note that no noise was added to the simulation;
what you are seeing are pure numerical errors getting amplified.
-If one of the gain peaks accumulates a lot of energy quickly ($L_\mathrm{NL}$ is small),
+If one of the gain peaks accumulates a lot of energy quickly ($$L_\mathrm{NL}$$ is small),
and that peak is in the anomalous dispersion regime,
then it can in turn also cause MI in its own surroundings,
leading to a cascade of secondary and tertiary gain areas.
-This is seen above for $z > 30 L_\mathrm{NL}$.
+This is seen above for $$z > 30 L_\mathrm{NL}$$.
What we described is "pure" MI, but there also exists
a different type caused by Raman scattering.
-In that case, amplification occurs at the strongest peak of the Raman gain $\tilde{g}_R(\omega)$,
+In that case, amplification occurs at the strongest peak of the Raman gain $$\tilde{g}_R(\omega)$$,
even when the parent pulse is in the NDR.
This is an example of stimulated Raman scattering (SRS).