358 words
2 minutes
Correlation Function and the Classical Limit

In quantum mechanics, the correlation function C(t)C(t) provides deep insights into the evolution of a state and its departure from classical behavior. This problem examines a state prepared as a superposition of energy eigenstates with a Gaussian distribution of energies.

Problem Statement#

Suppose you have prepared a state α|\alpha\rangle as a superposition of NN energy eigenstates wk|w_k\rangle with a mean energy ω0\omega_0 and a width σω\sigma_\omega: cke(ωkω0)22σω2c_k \propto e^{-\frac{(\omega_k - \omega_0)^2}{2\sigma_\omega^2}} α=kckwk|\alpha\rangle = \sum_k c_k |w_k\rangle

Our goal is to compute the normalized state and the correlation amplitude C(t)=αU(t)αC(t) = \langle \alpha | U(t) | \alpha \rangle.

Normalization#

To normalize the state, we require αα=1\langle \alpha | \alpha \rangle = 1: αα=k,lckclwkwl=kck2=1\langle \alpha | \alpha \rangle = \sum_{k,l} c_k^* c_l \langle w_k | w_l \rangle = \sum_k |c_k|^2 = 1

Given the proportionality ck=Ae(ωkω0)22σω2c_k = A e^{-\frac{(\omega_k - \omega_0)^2}{2\sigma_\omega^2}}, we find the normalization constant AA: A=1ke(ωkω0)2σω2A = \frac{1}{\sqrt{\sum_k e^{-\frac{(\omega_k - \omega_0)^2}{\sigma_\omega^2}}}}

The Correlation Function#

The time evolution operator is U(t)=eiHt/U(t) = e^{-iHt/\hbar}. Acting on the state α|\alpha\rangle: U(t)α=kckeiωktwkU(t) |\alpha\rangle = \sum_k c_k e^{-i\omega_k t} |w_k\rangle

The correlation function C(t)C(t) is defined as the overlap between the initial state and the state at time tt: C(t)=αU(t)α=kck2eiωktC(t) = \langle \alpha | U(t) | \alpha \rangle = \sum_k |c_k|^2 e^{-i\omega_k t}

Using the normalization and the Gaussian coefficients: C(t)=ke(ωkω0)2σω2eiωktke(ωkω0)2σω2C(t) = \frac{\sum_k e^{-\frac{(\omega_k - \omega_0)^2}{\sigma_\omega^2}} e^{-i\omega_k t}}{\sum_k e^{-\frac{(\omega_k - \omega_0)^2}{\sigma_\omega^2}}}

Classical Limit#

In the limit of a large number of states (NN \to \infty) and small spacing, we can approximate the sum as an integral: C(t)e(ωω0)2σω2eiωtdωe(ωω0)2σω2dωC(t) \approx \frac{\int e^{-\frac{(\omega - \omega_0)^2}{\sigma_\omega^2}} e^{-i\omega t} d\omega}{\int e^{-\frac{(\omega - \omega_0)^2}{\sigma_\omega^2}} d\omega}

The integral in the numerator is the Fourier transform of a Gaussian: eaω2+bωdω=πaeb2/4a\int_{-\infty}^{\infty} e^{-a\omega^2 + b\omega} d\omega = \sqrt{\frac{\pi}{a}} e^{b^2/4a}

After completing the square and performing the integration, we find: C(t)eiω0teσω2t2/4C(t) \approx e^{-i\omega_0 t} e^{-\sigma_\omega^2 t^2 / 4}

Interpretation#

The absolute value of the correlation function, C(t)=eσω2t2/4|C(t)| = e^{-\sigma_\omega^2 t^2 / 4}, shows a Gaussian decay. The “correlation time” τ\tau can be defined as the time it takes for C(t)|C(t)| to drop by 1/e1/e: τ=2σω\tau = \frac{2}{\sigma_\omega}

This relationship, τσω=2\tau \sigma_\omega = 2, is a manifestation of the energy-time uncertainty principle. As the width of the energy distribution σω\sigma_\omega increases (more “classical” or less well-defined energy), the state de-correlates much faster.

Correlation Function and the Classical Limit
https://rohankulkarni.me/posts/notes/correlation-function-classical-limit/
Author
Rohan Kulkarni
Published at
2023-07-20
License
CC BY-NC-SA 4.0