← Bibliothèque d'exercices← Tous les exercices du thème
Thème 3 — Postulats et applications

Transformée de Fourier d'une fonction d'onde exponentielle

Mots-clés: transformée de Fourier · représentation impulsion · distribution de Lorentz · relation d'incertitude · Plancherel

Exercice 1 : Transformée de Fourier d'une fonction d'onde exponentielle

Une particule sur une droite est décrite, à un instant donné, par la fonction d'onde

ψ(x)=1a e−∣x∣/a,a>0.\psi(x)=\frac{1}{\sqrt a}\,e^{-|x|/a},\qquad a>0 .

On rappelle que la fonction d'onde en représentation impulsion est

ψ~(p)=12πℏ∫Re−ipx/ℏψ(x) dx,\widetilde\psi(p)=\frac{1}{\sqrt{2\pi\hbar}}\int_{\R}e^{-ipx/\hbar}\psi(x)\,dx,

et que ∣ψ~(p)∣2|\widetilde\psi(p)|^2 est la densité de probabilité de l'impulsion. On pourra utiliser les intégrales

∫Rdu(1+u2)2=π2,∫Ru2 du(1+u2)2=π2,∫0∞xqe−x/bdx=q! bq+1.\int_{\R}\frac{du}{(1+u^2)^2}=\frac\pi2,\qquad\int_{\R}\frac{u^2\,du}{(1+u^2)^2}=\frac\pi2,\qquad\int_0^\infty x^qe^{-x/b}dx=q!\,b^{q+1}.
Question 1
Vérifier que ψ\psi est normalisée, et calculer ⟨X^⟩\langle\hat X\rangle et σX\sigma_X.

Solution
La fonction ∣ψ∣2=e−2∣x∣/a/a|\psi|^2=e^{-2|x|/a}/a est paire, et

∫R∣ψ∣2dx=2a∫0∞e−2x/adx=2a⋅a2=1.\int_{\R}|\psi|^2dx=\frac2a\int_0^\infty e^{-2x/a}dx=\frac2a\cdot\frac a2=1 .

Par parité, ⟨X^⟩=0\langle\hat X\rangle=0. La variance est donc

σX2=⟨X^2⟩=2a∫0∞x2e−2x/adx=2a×2!(a2)3=a22,σX=a2.\sigma_X^2=\langle\hat X^2\rangle=\frac2a\int_0^\infty x^2e^{-2x/a}dx=\frac2a\times2!\Bigl(\frac a2\Bigr)^3=\frac{a^2}{2}, \qquad\sigma_X=\frac{a}{\sqrt2}.

Question 2
Calculer ψ~(p)\widetilde\psi(p). Vérifier qu'elle est normalisée, conformément au théorème de Plancherel.

Solution
On découpe l'intégrale en deux, selon le signe de xx, en posant k=p/ℏk=p/\hbar :

∫Re−ikxe−∣x∣/adx=∫0∞e−(1/a+ik)xdx+∫−∞0e(1/a−ik)xdx=11/a+ik+11/a−ik=2/a1/a2+k2=2a1+k2a2.\begin{aligned} \int_{\R}e^{-ikx}e^{-|x|/a}dx&=\int_0^\infty e^{-(1/a+ik)x}dx+\int_{-\infty}^0e^{(1/a-ik)x}dx\\ &=\frac{1}{1/a+ik}+\frac{1}{1/a-ik}=\frac{2/a}{1/a^2+k^2}=\frac{2a}{1+k^2a^2}. \end{aligned}

Les deux intégrales convergent car les parties réelles des exposants sont négatives. On obtient

ψ~(p)=12πℏ⋅1a⋅2a1+p2a2/ℏ2=2aπℏ  11+p2a2/ℏ2.\widetilde\psi(p)=\frac{1}{\sqrt{2\pi\hbar}}\cdot\frac{1}{\sqrt a}\cdot\frac{2a}{1+p^2a^2/\hbar^2}=\sqrt{\frac{2a}{\pi\hbar}}\;\frac{1}{1+p^2a^2/\hbar^2}.

La densité d'impulsion est une fonction de Lorentz au carré, centrée en p=0p=0, de largeur caractéristique ℏ/a\hbar/a. Pour vérifier la normalisation, posons u=pa/ℏu=pa/\hbar, d'où dp=ℏ du/adp=\hbar\,du/a :

∫R∣ψ~(p)∣2dp=2aπℏ⋅ℏa∫Rdu(1+u2)2=2π⋅π2=1.\int_{\R}|\widetilde\psi(p)|^2dp=\frac{2a}{\pi\hbar}\cdot\frac\hbar a\int_{\R}\frac{du}{(1+u^2)^2}=\frac2\pi\cdot\frac\pi2=1 .

Question 3
Calculer ⟨P^2⟩\langle\hat P^2\rangle de deux façons : à partir de ψ~\widetilde\psi, puis en représentation position, à l'aide de ⟨P^2⟩=ℏ2∫∣ψ′(x)∣2dx\langle\hat P^2\rangle=\hbar^2\int|\psi'(x)|^2dx. En déduire σP\sigma_P.

Solution
En représentation impulsion, avec le même changement de variable,

⟨P^2⟩=∫Rp2∣ψ~(p)∣2dp=2aπℏ(ℏa)3∫Ru2du(1+u2)2=2ℏ2πa2⋅π2=ℏ2a2.\langle\hat P^2\rangle=\int_{\R}p^2|\widetilde\psi(p)|^2dp=\frac{2a}{\pi\hbar}\Bigl(\frac\hbar a\Bigr)^3\int_{\R}\frac{u^2du}{(1+u^2)^2}=\frac{2\hbar^2}{\pi a^2}\cdot\frac\pi2=\frac{\hbar^2}{a^2}.

En représentation position, la dérivée ψ′(x)=−sgn(x)aψ(x)\psi'(x)=-\frac{\mathrm{sgn}(x)}{a}\psi(x) est définie partout sauf en x=0x=0, où elle saute ; on a donc ∣ψ′∣2=∣ψ∣2/a2|\psi'|^2=|\psi|^2/a^2 presque partout, et

⟨P^2⟩=ℏ2∫R∣ψ′∣2dx=ℏ2a2∫R∣ψ∣2dx=ℏ2a2.\langle\hat P^2\rangle=\hbar^2\int_{\R}|\psi'|^2dx=\frac{\hbar^2}{a^2}\int_{\R}|\psi|^2dx=\frac{\hbar^2}{a^2}.

Les deux méthodes concordent. Comme ψ~\widetilde\psi est paire, ⟨P^⟩=0\langle\hat P\rangle=0, et σP=ℏ/a\sigma_P=\hbar/a. La formule ℏ2∫∣ψ′∣2\hbar^2\int|\psi'|^2 vient de l'intégration par parties ⟨ψ|P^2ψ⟩=⟨P^ψ|P^ψ⟩\braket{\psi}{\hat P^2\psi}=\braket{\hat P\psi}{\hat P\psi}, valable ici car ψ\psi est continue.

Question 4
Vérifier la relation d'incertitude de Heisenberg. Calculer la probabilité que l'impulsion soit comprise entre −ℏ/a-\hbar/a et ℏ/a\hbar/a.

Indication
Une primitive de 1/(1+u2)21/(1+u^2)^2 est u2(1+u2)+12arctan⁡u\frac{u}{2(1+u^2)}+\frac12\arctan u.

Solution
Le produit des écarts-types vaut σXσP=a2⋅ℏa=ℏ2≃0,71 ℏ\sigma_X\sigma_P=\frac{a}{\sqrt2}\cdot\frac\hbar a=\frac{\hbar}{\sqrt2}\simeq0{,}71\,\hbar, supérieur à ℏ/2\hbar/2 comme il se doit. Il ne dépend pas de aa : élargir la fonction d'onde en position réduit la dispersion en impulsion dans la même proportion. La probabilité demandée vaut

P(∣p∣<ℏa)=2π∫−11du(1+u2)2=4π[u2(1+u2)+12arctan⁡u]01=4π(14+π8)=12+1π≃0,82.P\Bigl(|p|<\frac\hbar a\Bigr)=\frac2\pi\int_{-1}^1\frac{du}{(1+u^2)^2}=\frac4\pi\Bigl[\frac{u}{2(1+u^2)}+\frac12\arctan u\Bigr]_0^1=\frac4\pi\Bigl(\frac14+\frac\pi8\Bigr)=\frac12+\frac1\pi\simeq0{,}82 .

Question 5
Montrer que ⟨P^4⟩\langle\hat P^4\rangle est infinie. Relier ce résultat à la régularité de ψ\psi en x=0x=0.

Solution
On a ⟨P^4⟩=2aπℏ(ℏa)5∫Ru4 du(1+u2)2\langle\hat P^4\rangle=\frac{2a}{\pi\hbar}\bigl(\frac\hbar a\bigr)^5\int_{\R}\frac{u^4\,du}{(1+u^2)^2}. Pour ∣u∣→∞|u|\to\infty, l'intégrande tend vers 11, et l'intégrale diverge. La distribution d'impulsion décroît trop lentement, comme 1/p41/p^4, pour que son moment d'ordre quatre existe. Ce comportement est lié au point anguleux de ψ\psi en x=0x=0 : la dérivée y présente un saut, donc la dérivée seconde contient un terme en delta de Dirac, ψ"=ψ/a2−2a3/2δ(x)\psi"=\psi/a^2-\frac{2}{a^{3/2}}\delta(x), qui n'est pas de carré intégrable. Ainsi P^2ψ\hat P^2\psi n'appartient pas à L2(R)L^2(\R), et ⟨P^4⟩=∥P^2ψ∥2\langle\hat P^4\rangle=\norm{\hat P^2\psi}^2 est infinie. Plus une fonction d'onde est irrégulière, plus sa transformée de Fourier décroît lentement, et plus elle contient de grandes impulsions.