Physlib

PhyslibAlpha.QuantumMechanics.HarmonicOscillator.Vacuum

A vacuum state for the harmonic oscillator

i. Overview

For a `d`-dimensional oscillator `Q`, `vacuumGaussian` is the Gaussian `exp(-∑ᵢ(xᵢ/ξᵢ)²/2)`, rescaled coordinatewise by the characteristic lengths `ξᵢ` through the diagonal continuous linear equivalence `diagEquiv`. Every annihilation operator kills this Schwartz function, so `hasVacuum_vacuumGaussian` supplies a `LadderSystem.HasVacuum` witness without an isotropy assumption.

The isotropic unit-length case (`Q.ξ i = 1` for all `i`) is represented by `stdGaussian (Space d) ℂ`, the function `x ↦ exp(-‖x‖²/2)` from `Physlib.Mathematics.InnerProductSpace.Gaussian`. `hasVacuum_stdGaussian_of_forall_xi_eq_one` derives the unit-length case from `vacuumGaussian`. The theorem `annihilationCLM_stdGaussian_of_xi_eq_one` assumes unit length only in the selected mode, so it is proved separately.

ii. Key results

- `diagEquiv` : the diagonal rescaling of `Space d` by the characteristic lengths. - `vacuumGaussian` : the anisotropic vacuum Gaussian, `exp(-∑ᵢ(xᵢ/ξᵢ)²/2)`. - `annihilationCLM_vacuumGaussian` : every `annihilationCLM i` kills the vacuum Gaussian. - `hasVacuum_vacuumGaussian` : `Q.vacuumGaussian` is a vacuum for `Q.toLadderSystem`. - `annihilationCLM_stdGaussian_of_xi_eq_one`, `hasVacuum_stdGaussian_of_forall_xi_eq_one` : the isotropic-unit-length special case.

iii. References

The standard Gaussian: derivative and momentum action

`Q`-free prerequisites about the plain, unscaled Gaussian, reused by both the general anisotropic case below and the isotropic-unit-length special case at the end of this file.

The general, anisotropic vacuum

No isotropy assumption needed: every oscillator `Q` (any masses/frequencies, hence any characteristic lengths `ξᵢ`) has an explicit vacuum, the Gaussian rescaled coordinatewise by `ξ`.

The isotropic-unit-length special case

18 declarations

theorem

iexp(12x2)=xiexp(12x2)\partial_i \exp(-\frac{1}{2} \|x\|^2) = -x_i \exp(-\frac{1}{2} \|x\|^2)

For any vector xSpace dx \in \text{Space } d and any coordinate index ii, the partial derivative of the real-valued Gaussian function f(x)=exp(12x2)f(x) = \exp\left(-\frac{1}{2} \|x\|^2\right) with respect to the ii-th coordinate is given by: iexp(12x2)=xiexp(12x2) \partial_i \exp\left(-\frac{1}{2} \|x\|^2\right) = -x_i \exp\left(-\frac{1}{2} \|x\|^2\right) where x\|x\| is the Euclidean norm and xix_i is the ii-th component of xx.

theorem

Pointwise value of the standard Gaussian G(x)=exp(12x2)\mathcal{G}(x) = \exp\left(-\frac{1}{2} \|x\|^2\right)

For any vector xx in the dd-dimensional real inner product space Space d\text{Space } d, the pointwise value of the standard Gaussian Schwartz function GS(Space d,C)\mathcal{G} \in \mathcal{S}(\text{Space } d, \mathbb{C}) is given by the real-valued Gaussian G(x)=exp(12x2)\mathcal{G}(x) = \exp\left(-\frac{1}{2} \|x\|^2\right) where x\|x\| is the Euclidean norm on Space d\text{Space } d.

theorem

Partial derivative of the standard Gaussian iΩ(x)=xiΩ(x)\partial_i \Omega(x) = -x_i \Omega(x)

Let Ω:Space dC\Omega: \text{Space } d \to \mathbb{C} be the standard Gaussian function defined by Ω(x)=exp(12x2)\Omega(x) = \exp\left(-\frac{1}{2} \|x\|^2\right). For any point xSpace dx \in \text{Space } d and any coordinate index i{0,,d1}i \in \{0, \dots, d-1\}, the partial derivative of Ω\Omega with respect to the ii-th coordinate satisfies: iΩ(x)=xiΩ(x)\partial_i \Omega(x) = -x_i \Omega(x) where xix_i denotes the ii-th component of the vector xx.

theorem

p^jΩ=ixjΩ\hat{p}_j \Omega = i \hbar x_j \Omega

Let ΩS(Rd,C)\Omega \in \mathcal{S}(\mathbb{R}^d, \mathbb{C}) be the standard Gaussian function defined by Ω(x)=exp(12x2)\Omega(x) = \exp\left(-\frac{1}{2} \|x\|^2\right). For any index j{0,,d1}j \in \{0, \dots, d-1\}, the jj-th component of the momentum operator p^j\hat{p}_j acting on Ω\Omega satisfies (p^jΩ)(x)=ixjΩ(x) (\hat{p}_j \Omega)(x) = i \hbar x_j \Omega(x) for all xRdx \in \mathbb{R}^d, where \hbar is the reduced Planck constant, ii is the imaginary unit, and xjx_j is the jj-th coordinate of the position vector xx.

definition

Diagonal rescaling map xiξixix_i \mapsto \xi_i x_i

For a dd-dimensional harmonic oscillator QQ, `diagEquiv` is the continuous linear equivalence from Space d\text{Space } d to itself that performs a diagonal rescaling of coordinates. For any point xSpace dx \in \text{Space } d and coordinate index i{0,,d1}i \in \{0, \dots, d-1\}, the ii-th component of the transformed vector is given by (Q.diagEquiv x)i=ξixi(Q.\text{diagEquiv } x)_i = \xi_i x_i, where ξi\xi_i is the characteristic length of the oscillator in the ii-th dimension.

theorem

(Q.diagEquiv x)i=ξixi(Q.\text{diagEquiv } x)_i = \xi_i x_i

Let QQ be a dd-dimensional harmonic oscillator. The diagonal rescaling operator Q.diagEquiv:Space dSpace dQ.\text{diagEquiv} : \text{Space } d \to \text{Space } d is defined such that for any vector xSpace dx \in \text{Space } d and any coordinate index i{0,,d1}i \in \{0, \dots, d-1\}, the ii-th component of the transformed vector is given by (Q.diagEquiv x)i=ξixi(Q.\text{diagEquiv } x)_i = \xi_i x_i where ξi\xi_i is the characteristic length of the oscillator in the ii-th dimension and xix_i is the ii-th component of the vector xx.

theorem

(Q.diagEquiv1x)i=ξi1xi(Q.\text{diagEquiv}^{-1} x)_i = \xi_i^{-1} x_i

Let QQ be a dd-dimensional harmonic oscillator with characteristic lengths ξi\xi_i for i{0,,d1}i \in \{0, \dots, d-1\}. Let Q.diagEquiv1Q.\text{diagEquiv}^{-1} be the inverse of the diagonal rescaling map. For any vector xSpace dx \in \text{Space } d, the ii-th component of the transformed vector is given by: (Q.diagEquiv1x)i=ξi1xi(Q.\text{diagEquiv}^{-1} x)_i = \xi_i^{-1} x_i where xix_i denotes the ii-th component of xx.

definition

Anisotropic vacuum Gaussian exp(12i(xi/ξi)2)\exp\left(-\frac{1}{2} \sum_i (x_i/\xi_i)^2\right) for a harmonic oscillator QQ

For a dd-dimensional harmonic oscillator QQ with characteristic lengths ξi\xi_i for i{0,,d1}i \in \{0, \dots, d-1\}, the vacuum Gaussian is the complex-valued Schwartz function ψS(Space d,C)\psi \in \mathcal{S}(\text{Space } d, \mathbb{C}) defined by: ψ(x)=exp(12i=0d1(xiξi)2)\psi(x) = \exp\left( -\frac{1}{2} \sum_{i=0}^{d-1} \left( \frac{x_i}{\xi_i} \right)^2 \right) for all xSpace dx \in \text{Space } d. This function is the anisotropic Gaussian obtained by applying the diagonal continuous linear equivalence Q.diagEquivQ.\text{diagEquiv} (which rescales coordinates as xiξixix_i \mapsto \xi_i x_i) to the standard Gaussian xexp(12x2)x \mapsto \exp(-\frac{1}{2}\|x\|^2).

theorem

Value of the vacuum Gaussian ψQ(x)=exp(12Q.diagEquiv1x2)\psi_Q(x) = \exp\left(-\frac{1}{2} \|Q.\text{diagEquiv}^{-1} x\|^2\right)

For a dd-dimensional harmonic oscillator QQ and any vector xSpace dx \in \text{Space } d, the value of the vacuum Gaussian function ψQ(x)\psi_Q(x) is given by ψQ(x)=exp(12Q.diagEquiv1(x)2) \psi_Q(x) = \exp\left( -\frac{1}{2} \|Q.\text{diagEquiv}^{-1}(x)\|^2 \right) where Q.diagEquiv1Q.\text{diagEquiv}^{-1} is the inverse of the diagonal rescaling map xiξixix_i \mapsto \xi_i x_i associated with the oscillator's characteristic lengths.

theorem

diagEquiv1(ei)=ξi1ei\text{diagEquiv}^{-1}(e_i) = \xi_i^{-1} e_i

For a dd-dimensional harmonic oscillator QQ, let eie_i be the ii-th vector of the standard orthonormal basis for Space d\text{Space } d, and let ξi\xi_i be the characteristic length of the oscillator in the ii-th dimension. The inverse of the diagonal rescaling map diagEquiv1\text{diagEquiv}^{-1} acts on the basis vector eie_i by scaling it by the reciprocal of ξi\xi_i: diagEquiv1(ei)=ξi1ei \text{diagEquiv}^{-1}(e_i) = \xi_i^{-1} e_i

theorem

Fréchet derivative of Q.diagEquiv1x2\|Q.\text{diagEquiv}^{-1} x\|^2 is 2Q.diagEquiv1x,Q.diagEquiv1()2 \langle Q.\text{diagEquiv}^{-1} x, Q.\text{diagEquiv}^{-1} (\cdot) \rangle

For a dd-dimensional harmonic oscillator QQ and any point xx in the dd-dimensional Euclidean space Space d\text{Space } d, the Fréchet derivative of the squared norm of the rescaled coordinate vector, f(x)=Q.diagEquiv1(x)2f(x) = \|Q.\text{diagEquiv}^{-1}(x)\|^2, is the continuous linear map given by: Df(x)=(2Q.diagEquiv1(x),)Q.diagEquiv1 Df(x) = (2 \langle Q.\text{diagEquiv}^{-1}(x), \cdot \rangle) \circ Q.\text{diagEquiv}^{-1} where Q.diagEquiv1Q.\text{diagEquiv}^{-1} is the inverse of the diagonal rescaling map xiξixix_i \mapsto \xi_i x_i, and ,\langle \cdot, \cdot \rangle denotes the standard real inner product on Space d\text{Space } d.

theorem

ig=(xi/ξi2)g\partial_i g = -(x_i/\xi_i^2)g for the Anisotropic Real Gaussian

For a dd-dimensional harmonic oscillator QQ with characteristic lengths ξj\xi_j, let g:Space dRg: \text{Space } d \to \mathbb{R} be the real-valued anisotropic Gaussian defined by g(x)=exp(12j=1dxj2ξj2). g(x) = \exp\left( -\frac{1}{2} \sum_{j=1}^d \frac{x_j^2}{\xi_j^2} \right). Then for any point xSpace dx \in \text{Space } d and coordinate index i{0,,d1}i \in \{0, \dots, d-1\}, the partial derivative of gg with respect to the ii-th coordinate is: ig(x)=xiξi2exp(12j=1dxj2ξj2). \partial_i g(x) = -\frac{x_i}{\xi_i^2} \exp\left( -\frac{1}{2} \sum_{j=1}^d \frac{x_j^2}{\xi_j^2} \right). Note that the expression Q.diagEquiv1(x)2\left\| Q.\text{diagEquiv}^{-1}(x) \right\|^2 in the formal statement corresponds to the sum j=1d(xj/ξj)2\sum_{j=1}^d (x_j/\xi_j)^2.

theorem

Coordinate Derivative of the Vacuum Gaussian: iΩ=xiξi2Ω\partial_i \Omega = -\frac{x_i}{\xi_i^2} \Omega

For a dd-dimensional harmonic oscillator QQ with characteristic lengths ξj\xi_j for j{0,,d1}j \in \{0, \dots, d-1\}, let Ω:Space dC\Omega: \text{Space } d \to \mathbb{C} be the vacuum Gaussian function defined by Ω(x)=exp(12j=0d1(xj/ξj)2)\Omega(x) = \exp\left( -\frac{1}{2} \sum_{j=0}^{d-1} (x_j/\xi_j)^2 \right). For any point xSpace dx \in \text{Space } d and coordinate index ii, the partial derivative of the vacuum Gaussian with respect to the ii-th coordinate is: iΩ(x)=xiξi2Ω(x)\partial_i \Omega(x) = -\frac{x_i}{\xi_i^2} \Omega(x) where xix_i is the ii-th coordinate of xx and ξi\xi_i is the characteristic length of the oscillator in that dimension.

theorem

Action of the Momentum Operator on the Vacuum Gaussian: p^iΩ=ixiξi2Ω\hat{p}_i \Omega = i \hbar \frac{x_i}{\xi_i^2} \Omega

For a dd-dimensional harmonic oscillator QQ, let Ω\Omega be the vacuum Gaussian state defined by Ω(x)=exp(12j=0d1(xj/ξj)2)\Omega(x) = \exp\left( -\frac{1}{2} \sum_{j=0}^{d-1} (x_j/\xi_j)^2 \right), where ξj\xi_j are the characteristic lengths. For any coordinate index i{0,,d1}i \in \{0, \dots, d-1\} and position xx, the action of the ii-th component of the momentum operator p^i\hat{p}_i on the vacuum Gaussian is given by: (p^iΩ)(x)=ixiξi2Ω(x) (\hat{p}_i \Omega)(x) = i \hbar \frac{x_i}{\xi_i^2} \Omega(x) where ii is the imaginary unit and \hbar is the reduced Planck constant.

theorem

The ii-th annihilation operator kills the vacuum Gaussian: a^iψQ=0\hat{a}_i \psi_Q = 0

For a dd-dimensional quantum harmonic oscillator QQ with characteristic lengths ξj\xi_j, let a^i\hat{a}_i be the annihilation operator for the ii-th mode and let ψQ\psi_Q be the vacuum Gaussian function defined by ψQ(x)=exp(12j=0d1(xjξj)2)\psi_Q(x) = \exp\left( -\frac{1}{2} \sum_{j=0}^{d-1} \left( \frac{x_j}{\xi_j} \right)^2 \right) for xRdx \in \mathbb{R}^d. Then for any i{0,,d1}i \in \{0, \dots, d-1\}, applying the annihilation operator to the vacuum Gaussian yields zero: a^iψQ=0.\hat{a}_i \psi_Q = 0.

theorem

The Vacuum Gaussian ψQ\psi_Q is a Vacuum State for the Ladder System of QQ

For a dd-dimensional harmonic oscillator QQ with characteristic lengths ξi\xi_i for i{0,,d1}i \in \{0, \dots, d-1\}, let ψQS(Space d,C)\psi_Q \in \mathcal{S}(\text{Space } d, \mathbb{C}) be the anisotropic vacuum Gaussian defined by ψQ(x)=exp(12i=0d1(xiξi)2).\psi_Q(x) = \exp\left( -\frac{1}{2} \sum_{i=0}^{d-1} \left( \frac{x_i}{\xi_i} \right)^2 \right). Then ψQ\psi_Q is a vacuum state for the ladder system of QQ, which means that for every mode ii, the annihilation operator a^i\hat{a}_i satisfies a^iψQ=0\hat{a}_i \psi_Q = 0. This result holds for any choice of characteristic lengths ξi\xi_i, without requiring isotropy.

theorem

a^i\hat{a}_i annihilates the standard Gaussian when ξi=1\xi_i = 1

Let QQ be a dd-dimensional quantum harmonic oscillator. For any mode index i{0,,d1}i \in \{0, \dots, d-1\}, let ξi\xi_i be the characteristic length of the ii-th mode and a^i\hat{a}_i be its corresponding annihilation operator. If ξi=1\xi_i = 1, then applying the annihilation operator a^i\hat{a}_i to the standard Gaussian function g(x)=exp(12x2)g(x) = \exp\left(-\frac{1}{2} \|x\|^2\right) yields zero: a^i(exp(12x2))=0\hat{a}_i \left( \exp\left(-\frac{1}{2} \|x\|^2\right) \right) = 0

theorem

Standard Gaussian is a vacuum if i,ξi=1\forall i, \xi_i = 1

Let QQ be a dd-dimensional quantum harmonic oscillator. If for every mode i{0,,d1}i \in \{0, \dots, d-1\} the characteristic length ξi\xi_i is equal to 11, then the standard Gaussian function g(x)=exp(12x2)g(x) = \exp\left(-\frac{1}{2} \|x\|^2\right) (viewed as a Schwartz map in S(Rd,C)\mathcal{S}(\mathbb{R}^d, \mathbb{C})) is a vacuum state for the ladder system associated with QQ. This implies that g(x)g(x) is annihilated by the annihilation operators a^i\hat{a}_i for all ii.