Quasi-probability distributions

From Wikipedia, the free encyclopedia

In the most general form, the dynamics of a quantum mechanical system are determined by a master equation - an equation of motion for the density operator (usually written \rho\,) of the system. Although it is possible to directly integrate this equation for very small systems (ie systems with few particles or degrees of freedom), this quickly becomes impossible for larger systems. For this reason, it is frequently useful to represent the density operator as a distribution over some (over-)complete operator basis. The evolution of the system is then completely determined by the evolution of a quasi-probability distribution function. This general technique has a long history, especially in the context of quantum optics. The most common examples of quasi-probability representations are the Wigner, P- and Q-functions. More recently, the positive P function and a wider class of generalized P functions have been used to solve complex problems in both quantum optics and the newer field of quantum atom optics.

Contents

Analogous to probability theory, quantum quasi-probability distributions can be written in terms of characteristic functions, from which all operator expectation values can be derived. The characteristic functions for the Wigner, Glauber P and Q distributions of an N mode system are as follows:

  • \chi_W(\mathbf{z},\mathbf{z}^*)=tr(\rho e^{i\mathbf{z}\cdot\hat{\mathbf{a}}+i\mathbf{z}^*\cdot\hat{\mathbf{a}}^{\dagger}})
  • \chi_P(\mathbf{z},\mathbf{z}^*)=tr(\rho e^{i\mathbf{z}\cdot\hat{\mathbf{a}}}e^{i\mathbf{z}^*\cdot\hat{\mathbf{a}}^{\dagger}})
  • \chi_Q(\mathbf{z},\mathbf{z}^*)=tr(\rho e^{i\mathbf{z}^*\cdot\hat{\mathbf{a}}^{\dagger}}e^{i\mathbf{z}\cdot\hat{\mathbf{a}}})

Here \hat{\mathbf{a}} and \hat{\mathbf{a}}^{\dagger} are vectors containing the annihilation and creation operators for each mode of the system. These characteristic functions can be used to directly evaluate expectation values of operator moments. The ordering of the annihilation and creation operators in these moments is specific to the particular characteristic function. For instance, normally ordered (annihilation operators preceding creation operators) moments can be evaluated in the following way from \chi_P\,:

\langle\hat{a}_j^{\dagger m}\hat{a}_k^n\rangle = \frac{\partial^{m+n}}{\partial(iz_j^*)^m\partial(iz_k^*)^n}\chi_P(\mathbf{z},\mathbf{z}^*)\Big|_{\mathbf{z}=\mathbf{z}^*=0}

In the same way, expectation values of anti-normally ordered and symmetrically ordered combinations of annihilation and creation operators can be evaluated from the characteristic functions for the Q and Wigner distributions, respectively.

The quasi-probability functions themselves are defined as Fourier transforms of the above characteristic functions. That is,

\{W|P|Q\}(\mathbf{\alpha},\mathbf{\alpha}^*)=\frac{1}{\pi^{2N}}\int d^{2N}\mathbf{z}\chi_{\{W|P|Q\}}(\mathbf{z},\mathbf{z}^*)e^{-i\mathbf{z}^*\cdot\mathbf{\alpha}^*}e^{-i\mathbf{z}\cdot\mathbf{\alpha}}

Here \alpha_j\, and \alpha^*_k may be identified as coherent state amplitudes in the case of the Glauber P and Q distributions, but simply c-numbers for the Wigner function. Since differentiation in normal space becomes multiplication in fourier space, moments can be calculated from these functions in the following way:

  • \langle\hat{\mathbf{a}}_j^{\dagger m}\hat{\mathbf{a}}_k^n\rangle=\int d^{2N}\mathbf{\alpha}P(\mathbf{\alpha},\mathbf{\alpha}^*)\alpha_j^n\alpha_k^{*m}
  • \langle\hat{\mathbf{a}}_j^m\hat{\mathbf{a}}_k^{\dagger n}\rangle=\int d^{2N}\mathbf{\alpha}Q(\mathbf{\alpha},\mathbf{\alpha}^*)\alpha_j^m\alpha_k^{*n}
  • \langle(\hat{\mathbf{a}}_j^{\dagger m}\hat{\mathbf{a}}_k^n)_S\rangle=\int d^{2N}\mathbf{\alpha}W(\mathbf{\alpha},\mathbf{\alpha}^*)\alpha_j^m\alpha_k^{*n}

Here (\ldots)_S denotes symmetric ordering.

These relationships motivate comparisons between the distribution functions and classical probability densities. The analogy - though strong - is not perfect, as the above functions are not necessarily positive for all \mathbf{\alpha}. Hence the term quasi-probability function.

Since each of the above transformations from \rho\, through to the distribution function is linear, the equation of motion for each distribution can be obtained by performing the same transformations to \dot{\rho}. Furthermore, as any master equation which can be expressed in Lindblad form is completely described by the action of combinations of annihilation and creation operators on the density operator, it is useful to consider the effect such operations have on each of the quasi-probability functions.

For instance, consider the annihilation operator \hat{a}_j\, acting on \rho\,. For the characteristic function of the P distribution we have

tr(\hat{a}_j\rho e^{i\mathbf{z}\cdot\hat{\mathbf{a}}}e^{i\mathbf{z}^*\cdot\hat{\mathbf{a}}^{\dagger}}) = \frac{\partial}{\partial(iz_j)}\chi_P(\mathbf{z},\mathbf{z}^*)

Taking the Fourier transform with respect to \mathbf{z}\, to find the action corresponding action on the Glauber P function, we find

\hat{a}_j\rho \rightarrow \alpha_j P(\mathbf{\alpha},\mathbf{\alpha}^*).

By following this procedure for each of the above distributions, the following operator correspondences can be identified:

  • \hat{a}_j\rho \rightarrow \left(\alpha_j + s\frac{\partial}{\partial\alpha_j^*}\right)P_s(\mathbf{\alpha},\mathbf{\alpha}^*)
  • \rho\hat{a}^{\dagger}_j \rightarrow \left(\alpha_j^* + s\frac{\partial}{\partial\alpha_j}\right)P_s(\mathbf{\alpha},\mathbf{\alpha}^*)
  • \hat{a}^{\dagger}_j\rho \rightarrow \left(\alpha_j^* - (1-s)\frac{\partial}{\partial\alpha_j}\right)P_s(\mathbf{\alpha},\mathbf{\alpha}^*)
  • \hat{a}_j\rho \rightarrow \left(\alpha_j - (1-s)\frac{\partial}{\partial\alpha_j^*}\right)P_s(\mathbf{\alpha},\mathbf{\alpha}^*)

Here s = 0, 1/2 or 1 for P, Wigner and Q distributions, respectively.

In this way, master equations can be expressed as an equations of motion of quasi-probability functions.

Consider a single-mode system evolving under the following Hamiltonian operator:

\hat{H} = \frac{\hbar}{2}\hat{a}^{\dagger 2}\hat{a}^2

Image:Bvn-small.png Probability distributionsview  talk  edit ]
Univariate Multivariate
Discrete: BenfordBernoullibinomialBoltzmanncategoricalcompound PoissondegenerateGauss-Kuzmingeometrichypergeometriclogarithmicnegative binomialparabolic fractalPoissonRademacherSkellamuniformYule-SimonzetaZipfZipf-Mandelbrot Ewensmultinomialmultivariate Polya
Continuous: BetaBeta primeCauchychi-squareDirac delta functionErlangexponentialexponential powerFfadingFisher's zFisher-TippettGammageneralized extreme valuegeneralized hyperbolicgeneralized inverse GaussianHalf-LogisticHotelling's T-squarehyperbolic secanthyper-exponentialhypoexponentialinverse chi-square (scaled inverse chi-square)• inverse Gaussianinverse gamma (scaled inverse gamma) • KumaraswamyLandauLaplaceLévyLévy skew alpha-stablelogisticlog-normalMaxwell-BoltzmannMaxwell speednormal (Gaussian)normal inverse GaussianParetoPearsonpolarraised cosineRayleighrelativistic Breit-WignerRiceshifted GompertzStudent's ttriangulartype-1 Gumbeltype-2 GumbeluniformVariance-GammaVoigtvon MisesWeibullWigner semicircleWilks' lambda Dirichletinverse-WishartKentmatrix normalmultivariate normalmultivariate Studentvon Mises-FisherWigner quasiWishart
Miscellaneous: Cantorconditionalexponential familyinfinitely divisiblelocation-scale familymarginalmaximum entropyphase-typeposteriorpriorquasisamplingsingular

  • H. J. Carmichael, Statistical Methods in Quantum Optics I: Master Equations and Fokker-Planck Equations, Springer-Verlag (2002).
  • C. W. Gardiner, Quantum Noise, Springer-Verlag (1991).
Advanced Search
Included Web Search Engines


Safe Search

close

Top Matching Results

Occasionally Search.com will highlight specialized results that are based on the context of your query. Examples of specialized results include specific links to news, images, or video.

Top Matching Results may highlight information from other Search.com pages, content from the CNET Network of sites, or third party content. The listings are based purely on relevance. Search.com does not receive payment for listings in this section but our partners that provide this data may get paid for listing these products.

Sponsored Links

This section contains paid listings which have been purchased by companies that want to have their sites appear for specific search terms and related content. These listings are administered, sorted and maintained by a third party and are not endorsed by Search.com.

Search Results

Search.com sends your search query to several search engines at one time and integrates the results into one list which has been sorted by relevance using Search.com's proprietary algorithm. You can customize the list of search engines included in your metasearch from the preferences.

The search engines that are used in your metasearch may allow companies to pay to have their Web sites included within the results. To view the Paid Inclusion policy for a specific search engine, please visit their Web site. Search.com does not accept payment or share revenue with any search engine partner for listings in this section.