User:Phenomenam/sandbox

Introduction
In order to capture the interference phenomena that occur during a scattering process that the classical theory fails to address properly, yet to avoid a fully quantum treatment that is complicated, a semiclassical approach may be employed (see Quantum Scattering Theory for a fully quantum treatment and Classical Scattering Theory for a classical treatment of a scattering process). Semiclassical scattering theory makes use of trajectories, impact parameters, and the like that are classical and the superposition of scattering waves that are quantum to approximate or model a scattering phenomenon fairly accurately under the condition where each classical trajectories carry unique phases.

Definition
Semiclassical physics is "a vast discipline of theoretical physics that attempts to compute spectra and wave functions of quantum (wave) systems on the basis of their classical trajectories (rays)" with the addition of quantum superposition principle and Planck's constant $$h$$. (see also Semiclassical physics)

Usefulness
In general, semiclassical approach is said to be useful since very few analytically solvable quantum systems are known to date. Usually, semiclassical methods are incorporated as they allow us to compute spectra, i.e. energy levels, and wave functions of quantum systems based on classical trajectories even in the absence of the exact quantum approach. Note that the method still involves quantum phases and the quantum superposition principle.

Background
This section largely follows section 5.2 of the reference.

In order to obtain approximate solutions to a quantum system with potential $$V(x)$$, we assume stationary wave function $$\psi(x)$$



\psi(x) = A(x) e^{\frac{i}{\hbar} B(x)}, $$

where $$A(x)$$ and $$B(x)$$ are chosen such that the wave function satisfies Schrödinger equation. Inserting the form into Schrödinger equation and solving for the equation generally, we obtain a "classical" (absence of $$\hbar$$) equation



\frac{1}{2m} (B')^2 + V - E = 0 $$

and a "quantum" equation

\hbar A + 2iA'B' + iAB = 0. $$ Solving for $$B(x)$$ using the classical equation,



B(x) = \pm \int_{a}^{x} p(x) {\mathrm dx} = \pm S(x;E), $$

where $$a$$ is an arbitrary location in the classically allowed region, $$p(x)$$ is the momentum $$\sqrt{2m[V(x)-E]} $$, and $$ S(x;E) $$ is the classical action. The importance of the result lies in the fact that the quantum phase angle shift $$B(x)$$ of the wave function $$\psi(x)$$ is now directly related to the classical action $$S(x;E)$$, which serves as the backbone of semiclassical treatments.

Semiclassical Approximation
Among many semiclassical approximations (see Semiclassical physics), we focus on Bohr-Sommerfeld Quantization and JWKB Approximation that are most useful in dealing with scattering problems.

Bohr-Sommerfeld Quantization
Dating from the early 20th century, Bohr-Sommerfeld Quantization procedure forms the basis of the "old quantum mechanics''. The approximation taken in this procedure is to assume the amplitude $$ A(x) $$ to be constant, namely,



\psi(x) \approx N e^{\frac{i}{\hbar}B(x)}. $$

Computing the accumulated phase angle for the wave that takes off from the starting position, say $$a$$, then returns to its initial position, we obtain

\Delta B(a) = \oint p(x) {\mathrm dx} = S(E). $$ Since the wave function $$\psi(x)$$ is stationary and unique,

S(E_n) = 2\pi n \hbar, $$ where $$n$$ is a positive integer.

JWKB Approximation
In the limit $$\hbar$$ is small, which is justified by the fact that we are in semiclassical realm, the "quantum" equation (see Background section above) may be approximated to be

2A'B' + A'B'' = 0, $$ which has a solution $$A(x)$$ of the form

A(x) = \frac{C}{\sqrt{|p(x)|}}, $$ where $$C$$ is a complex number. With $$B(x) = \pm S(x;E)$$ found in Background section, the wave function $$\psi(x)$$ in the classically allowed region may be written as

\psi(x) = \frac{Ce^{iS/\hbar} + De^{-iS/\hbar}}{\sqrt{|p(x)|}}, $$ and the wave function $$\psi(x)$$ in the classically forbidden region may be written as

\psi(x) = \frac{Fe^{S/\hbar} + Ge^{-S/\hbar}}{\sqrt{|p(x)|}}, $$ where $$D$$, $$F$$, and $$G$$ are complex numbers. At this point, we notice that the obtained wave functions are not valid at the classical turning point due to their non-zero momenta. Using linear approximation of the shape of $$V(x)$$ at the turning point, with the help of Airy functions and the boundary conditions on $$\psi(x)$$ and $$\psi'(x)$$ (see WKB approximation,, or for details), we obtain the wave function in the classically allowed region to be

\psi(x) \sim \sin\left(\frac{\pi}{4} + \frac{1}{\hbar} \int_{x_0}^{x} |p(x)| {\mathrm dx} \right), $$ where $$x_0$$ denotes the location of the classical turning point. For the sake of completeness, we note in passing that in the case where there are two classical turning points, say $$x_1$$ and $$x_2$$, which would not be of interest for scattering purpose, one may show

\Delta B = S(E_n) = 2\pi\hbar \left(n-\frac{1}{2}\right), $$ as an analogous expression to the classical action found in Bohr-Sommerfeld Quantization. In the case where there is only one classical turning point $$x_0$$, which is of interest for investigating a scattering system, a more involved approach is needed (see below).

Applications to scattering theory
Applying the semiclassical approximation techniques to a scattering system is relatively straightforward. The theoretical treatment of a scattering system states (i) mass $$m$$ $$\rightarrow$$ reduced mass $$\mu$$, (ii) potential $$V$$ $$\rightarrow$$ effective potential $$V_$$, and (iii) starting position $$a$$ $$\rightarrow$$ far limit $$\infty$$. Applying these changes to Bohr-Sommerfeld Quantization, we have a situation where the phase of the wave function oscilates on a nearly constant background provided by the potential, with the fictitious particle with mass $$\mu$$ coming in from $$\infty$$, rebounding from the inner wall of the effective potential $$V_$$, and returning to $$\infty$$. During this process, the particle or the wave function, then, picks up an extra phase angle $$\Delta B$$ whose 1-D analytical computation is shown above. When we apply the same changes to JWKB approximation, in analogy to the 1-D case shown above, we obtain the 3-D wave function at far limit to be

\psi^{(>)}_l (r) \approx \sin\left(\frac{\pi}{4} + \int_{r_0}^{r} |p_l(r')| {\mathrm dr'} \right), $$ where

\frac{p^2_l(r')}{2\mu} = E - V_(r') = E - V(r') - \frac{\left(l+\frac{1}{2}\right)^2}{r'^2}. $$ Rewriting the far limit wave function as

\psi^{(>)}_l (r) = \sin\left(\frac{\pi}{4} + \frac{1}{\hbar} \left\{\int_{r_0}^{r} \left[|p_l(r')| - \sqrt{2\mu E} \right] {\mathrm dr'} + \sqrt{2\mu E}(r-r_0)\right\} \right), $$ and comparing it with the form of the far limit of the exact quantum solution

\psi^{(>),({\rm exact})}_l(r) \approx \frac{\sqrt{2\mu E}}{\hbar}r \times j_l\left(\frac{\sqrt{2\mu E}}{\hbar} r\right), $$ where $$j_l$$ is a spherical Bessel function of the first kind, together with the large $$r$$ approximation of $$j_l$$

j_l^{(>)}(kr) \approx \frac{\sin\left(kr - \frac{l\pi}{2}\right)}{kr}, $$ we see that the accumulated or extra phase picked up via scattering in the limit $$r \rightarrow \infty$$ is given by

\delta_l^ = \left(l+\frac{1}{2}\right)\frac{\pi}{2} -kr_0 + \int_{r_0}^{\infty} [k(r')-k] {\mathrm dr'}, $$ where we used

\frac{|p(r')|}{\hbar} = k(r') $$ and

\frac{\sqrt{2\mu E}}{\hbar} = k. $$ The usefulness of the semiclassical $$\delta_l$$ expression above is apparent given the wide use of $$\delta_l$$ in finding the following scattering related quantities. The list includes but is not limited to the total cross section

\sigma = \frac{4\pi}{k^2} \sum_{l=0}^{\infty} (2l+1) \sin^2 \delta_l, $$ the scattering length

a_l(k) = \frac{\tan \delta_l(k)}{k}, $$ etc.

Connection to Classical and Quantum Scattering
Previously, we have shown that the classical action in semiclassical physics corresponds to the integral of the momentum of a given system over displacement, which, in 3-D, is written as

S(\vec{r};E,b) = \int_{\vec{r}_}^{\vec{r}} \vec{p}(\vec{r}') \cdot {\mathrm d\vec{r}'}, $$ where $$b$$ is the impact parameter. Dividing the coordinates into the radial and the angular parts, namely,

\vec{p} \cdot {\mathrm d\vec{r}} = p_{r} {\mathrm dr} + p_{\varphi} r{\mathrm d\varphi}, $$ and using the fact that the difference between the classical action in the presence of the potential and the classical action in the absence of the potential is the scattering phase shift $$S(E,b)$$ when integrated along a closed path, we write in general, with $$\chi'(E,b)$$ the deflection angle and $$L$$ the angular momentum $$p_{\varphi} r$$ of the system,

S(E,b) = \Delta (E,b) - L\chi'(E,b), $$ where

\Delta (E,b) = \oint p_r(r') {\mathrm dr'} - \oint p^{V=0}_r(r') {\mathrm dr'} = \frac{\delta^_l}{2\hbar}, $$ which has the name "classical phase", and

L\chi' = \int_{0}^{\pi-\chi'} p_{\varphi} r {\mathrm d\varphi} - \int_{0}^{\pi} p_{\varphi} r {\mathrm d\varphi}. $$ We note that, with

L = \mu v_{\infty} b, $$ where $$v_{\infty}$$ is the asymptotic relative velocity of the two bodies participating in a scattering process, and

p_r(r') = 2\mu v_{\infty} \sqrt{1-\frac{b^2}{r'^2}-\frac{V(r')}{E}}, $$ the first $$L$$ derivative of $$\Delta$$ becomes

\frac{\partial \Delta}{\partial L}\bigg|_E = \chi(E,b), $$ so that

\frac{\partial S}{\partial L}\bigg|_E = \chi(E,b) - \chi'. $$ Applying the principle of least action, the action derivative, then, must read zero, meaning $$\chi(E,b) = \chi'$$, hence it is the classically allowed path. Said differently, the waves that are on the paths near the classically allowed one constructively interfere with one another whereas those that are not interfere destructively.

Comparing the two momenta expressions $$p_r(r')$$ above and $$p_l(r')$$ in the previous section, we obtain

b = \frac{\lambda}{2\pi} \left(l+\frac{1}{2}\right), $$ where

\lambda \equiv \frac{2\pi}{k} $$ with $$k$$ defined in the previous section. The expression bears an importance as it connects the specific impact parameters with corresponding angular momentum quantum numbers, which may be used in the following ways.

Recalling the quantum mechanical expression for the total cross section above, we may now rewrite $$\sigma$$ as

\sigma = 8\pi \int_0^{\infty} b \sin^2\left(\frac{\Delta(E,b)}{2\hbar}\right) {\mathrm db}. $$ In order to see the difference that semiclassical methods make, suppose we have a hard sphere (radius $$R$$) scattering, for instance. The evaluation of the integral results in $$\sigma = 2\pi R^2$$, which is twice the classical value. The existence of additional contribution to $$\sigma$$, then, demonstrates the ability of semiclassical method capturing the effect of diffraction of the incident plane wave around the target.

This time, we take a look at the JWKB phase shift $$\delta_l^$$. Instead of comparing the semiclassical wave function with the form of the exact quantum solution as in previous section, we write directly;

\delta_l^ = 2\mu v_{\infty} \left[ \int_{r_0}^{\infty} \left(1-\frac{b^2}{r^2}-\frac{V(r)}{E}\right)^{\frac{1}{2}} {\mathrm dr} -\int_{b}^{\infty} \left(1-\frac{b^2}{r^2}\right)^{\frac{1}{2}} {\mathrm dr} \right]. $$ In the limit of large angular momentum, since the contribution from the potential $$V(r)$$ would be small, the turning point $$r_0$$ would approximately be $$b$$. Combining the two integrals with the common factor $$1-b^2/r^2$$ and using Taylor expansion up to first order in $$\frac{V(r)/E}{(1-b^2/r^2)^{1/2}}$$, one can show

\delta_l^ \approx -\frac{1}{2} \int_{b}^{\infty} \frac{V(r)/E}{\left(1-\frac{b^2}{r^2}\right)^{\frac{1}{2}}} {\mathrm dr}, $$ which was obtained by Massey and Mohr in 1934 . The expression may then be used to compute the total cross section given the long-range potential of the form $$\sim \frac{1}{r^n}$$, which is known as the Massey-Mohr formula.

Finally, consider the partial wave expansion of the (axially-symmetric) scattering amplitude, the amplitude of the scattered wave function whose radial dependence is $$\sim \frac{e^{ikr}}{r}$$,

f(\theta) = \frac{1}{2ik} \sum_{l=0}^{\infty} (2l+1)(e^{2i\delta_l}-1) P_l(\cos\theta), $$ where $$\delta_l$$ is the extra phase picked up by the presence of potential in the scattering process associated with the angular quantum number $$l$$ and $$P_l$$ is the Legendre polynomials. Assuming those terms in the sum with $$l\theta \leq 1$$ are small, replacing $$P_l(\cos\theta)$$ in the sum with an approximate formula that holds for $$l\theta \geq 1$$

P_l(\cos\theta) \approx \frac{\sin\left[\left(l+\frac{1}{2}\right)\theta+\frac{\pi}{4}\right]} {\sqrt{\frac{\pi}{2}\left(l+\frac{1}{2}\right)\sin\theta}} $$ results in

f(\theta) = \frac{-1}{\sqrt{\lambda\sin\theta}} \int_{0}^{\infty} \sqrt{b}(e^{i\varphi_{+}}-e^{i\varphi_{-}}){\mathrm db}, $$ where we first replaced $$l$$ and $$k$$ with the correspondent $$b$$ and $$\lambda$$ expressions respectively (see above) then used integral approximation on $$b$$ and

\varphi_{\pm} = 2\delta(b) \pm \left(\frac{2\pi b\theta}{\lambda}+\frac{\pi}{4}\right). $$ Now relaxing the positivity condition of $$b$$, with $$\theta = |\chi'|$$, we obtain

f(\theta) = \frac{e^{-i\frac{\pi}{4}}}{\sqrt{\lambda\sin\theta}} \int_{-\infty}^{+\infty} \sqrt{b} e^{i\frac{S(\theta,b)}{\hbar}} {\mathrm db}, $$ where we have $$\theta = \chi'>0$$ for the repulsive trajectory with $$b>0$$ and $$-\theta = \chi' < 0$$ for the attractive trajectory with $$b<0$$. We note that we have replaced the argument $$E$$ of $$S$$ with $$\theta$$ since $$\theta(E,b) = |\chi'(E,b)|$$ is a function of $$E$$ and $$b$$. Incorporating the stationary phase or the steepest descent method, justified by the fact that the action must not change rapidly in $$b$$ if the integral were not to be 0 in consistency with the least action principle, $$f(\theta)$$ now evaluates to

f(\theta) = \sum_{b_i(\theta)} \sqrt{\frac{2\pi b_i(\theta)}{\lambda\sin\theta\Delta''(b_i)/\hbar}} e^{iS(\theta,b_i)/\hbar}, $$ where we used

\frac{\partial^2S}{\partial b^2} = \frac{\partial^2\Delta}{\partial b^2} $$ and $$b_i$$ denotes the stationary point. Obviously, if there are more than one $$b_i$$, the different paths will cause interference phenomena. Solving for the spacing of the maxima given two different $$b_i$$'s, namely $$b_1(\theta)$$ and $$b_2(\theta)$$, we obtain

\Delta\theta = \frac{\lambda}{|b_1(\theta) - b_2(\theta)|}, $$ which is nothing but the spacing of interference maxima in Young's double-slit experiment, for which $$\Delta\theta = \lambda/2\pi d$$; multiple amplitudes will lead to multiple sets of interference minima and maxima. In analogy to the double-slit experiment, we conclude that the greater the spacing between the $$b_i$$'s are, the closer the spacing of the fringes, a.k.a. diffraction oscillations.

Limitations
Semiclassical approximations taken here do have their own limitations. Among many other possible limitations, the notable one includes the breakdown of $$f(\theta)$$ expression at $$\theta=\theta_R$$, the rainbow angle. This is so, because the stationary phase approximation does not hold at the rainbow angle as the second derivative of action disappears. This failure of the primitive semiclassical approximation may be saved by

S(\theta,b) \approx S(\theta_R,b_R) + \frac{1}{6} \frac{\partial^3S}{\partial b^3}\bigg|_{b=b_R} (b-b_R)^3, $$ where $$b_R$$ is the rainbow impact parameter. Because now the integral that needs to be done includes an exponential in a cubic power of $$b$$, the final form of the integral is an Airy function.