Approximate density matrices and Husimi functions using the maximum entropy formulation with constraints

This paper considers a Husimi representation of quantum mechanics in which the (stationary) state of a system or ensemble is described by a Husimi function and an observable is described by a phase space function or distribution such that the expectation value of the observable is given by an integral over phase space of the product of that function or distribution and the Husimi function. The density matrix, Wigner function, and Husimi function are considered to be alternative ways of describing the state of a system or ensemble, and methods of recovering the Wigner function or density matrix from the Husimi function are discussed. The classical limits of the Wigner and Husimi functions and of the relationship between them are considered. © 1993 John Wiley & Sons, Inc.

We propose localization measures in phase space of the ground state of bilayer quantum Hall systems at fractional filling factors , to characterize the three quantum phases (shortly denoted by spin, canted and ppin) for arbitrary -isospin λ. We use a coherent state (Bargmann) representation of quantum states, as holomorphic functions in the 8-dimensional Grassmannian phase-space (a higher-dimensional generalization of the Haldane’s 2-dimensional sphere ). We quantify the localization (inverse volume) of the ground state wave function in phase-space throughout the phase diagram (i.e. as a function of Zeeman, tunneling, layer distance, etc, control parameters) with the Husimi function second moment, a kind of inverse participation ratio that behaves as an order parameter. Then we visualize the different ground state structure in phase space of the three quantum phases, the canted phase displaying a much higher delocalization (a Schrödinger cat structure) than the spin and ppin phases, where the ground state is highly coherent. We find a good agreement between analytic (variational) and numeric diagonalization results.

For the representation of eigenstates on a Poincare section at the boundary of a billiard different variants have been proposed. We compare these Poincare Husimi functions, discuss their properties, and based on this select one particularly suited definition. For the mean behavior of these Poincare Husimi functions an asymptotic expression is derived, including a uniform approximation. We establish the relation between the Poincare Husimi functions and the Husimi function in phase space from which a direct physical interpretation follows. Using this, a quantum ergodicity theorem for the Poincare Husimi functions in the case of ergodic systems is shown.

We study the efficiency of quantum algorithms which aim at obtaining phase-space distribution functions of quantum systems. Wigner and Husimi functions are considered. Different quantum algorithms are envisioned to build these functions, and compared with the classical computation. Different procedures to extract more efficiently information from the final wave function of these algorithms are studied, including coarse-grained measurements, amplitude amplification, and measure of wavelet-transformed wave function. The algorithms are analyzed and numerically tested on a complex quantum system showing different behavior depending on parameters: namely, the kicked rotator. The results for the Wigner function show in particular that the use of the quantum wavelet transform gives a polynomial gain over classical computation. For the Husimi distribution, the gain is much larger than for the Wigner function and is larger with the help of amplitude amplification and wavelet transforms. We discuss the generalization of these results to the simulation of other quantum systems. We also apply the same set of techniques to the analysis of real images. The results show that the use of the quantum wavelet transform allows one to lower dramatically the number of measurements needed, but at the cost of a large loss of information.

The Wigner formalism is a convenient way of describing quantum mechanical effects through a framework of distribution functions in phase space. Currently, there are stochastic and deterministic approaches in use. In our deterministic method, the critical discretization of the diffusion term is done through the utilization of an integral formulation of the Wigner equation. This deterministic method is studied in the context of superposed quantum states as a precursor to simulations of entangled states.

A nonseparable wave-like integro-differential equation for the time evolution of the Wigner distribution function in phase space is educed from the corresponding separable kinetic equation. By employing the quantum hydrodynamical description, a non-local evolution wave equation is also derived by synthesizing the Hamilton‐Jacobi equation with that of continuity, which predicts the generation of nonlocal and quadrupole quantum phenomena in the propagation of the spatial probability density.

A plasma consists of charged particles moving relatively freely, but interacting through the electric and magnetic fields produced by each particle. It is well known that all except a very small (but nevertheless important) fraction of the detectable universe is a plasma. Most plasma physics can be expressed in terms of five fairly straightforward equations, the four Maxwell equations for electromagnetism and the equation of motion for charged particles in electromagnetic fields. Why then is plasma physics such a mathematically difficult subject and why is numerical simulation so important? The most basic reason is that any experimental plasma contains a very large number of particles, 1018–1020 in a typical fusion experiment. Moreover, particles in plasmas have a tendancy to act independently of each other. In conventional hydrodynamics, large numbers of particles can be described by a few simple numbers such as their temperature, density and mean velocity. Sometimes we can do this with a plasma but even then we have a few more numbers to include such as current and charge density, along with the electromagnetic fields. More usually, plasmas are not so simply described. Plasmas are not usually Maxwellian and we have to think in terms of distribution functions in phase space. Whenever we can, we model the non—Maxwellian nature of the plasma by transport models, e.g. for heat flow, but this is a specialised skill in itself. When one considers that conventional hydrodynamics is very complex and involves advanced numerical simulation, it is hardly surprising that simulation plays such an important role in plasma physics.

A new scheme for the calculation of the collision integral of the nonlinear Boltzmann kinetic equation is presented, based on a piecewise constant discrete approximation of the distribution function in phase space and corresponding precalculated collision rate coefficients. The computing effort required for evaluation of the fivefold collision integral is significantly reduced by the use of precalculated rate coefficients. The method was combined with the conservative splitting method to solve the Boltzmann equation for a number of test cases to examine the performance of this method and demonstrate its efficiency. 10 refs., 6 figs.

Analytical solution for self-consistent particle equilibrium distribution in RF field with transverse focusing is found. Solution is attained in approximation of high brightness beam. Distribution function in phase space is determined as a stationary function of the energy integral. Equipartitioning for beam distribution between degrees of freedom follows directly from the choice of stationary distribution function. Analytical expressions for r-z equilibrium beam profile and maximum beam current in RF field are obtained.

The Rouse model for dilute polymer solutions undergoing homogeneous flows has been generalized to include the inertia of the beads in the equations of motion. To obtain the correct ‘‘diffusion equation’’ for the probability density distribution function in phase space, we generalize the diffusion equation derived by Murphy and Aguirre [J. Chem. Phys. 57, 2098 (1972)] from Hamilton’s equations of motion for an arbitrary number of interacting Brownian particles at equilibrium. Material functions are found, and the noninertial case is seen to be obtained as the zero mass limit in all steps of the development. In particular, the steady-state shear results are unaffected by the inclusion of inertia. It is also shown how two assumptions, ‘‘equilibration in momentum space,’’ and the neglect of acceleration, made independently by Curtiss, Bird, and Hassager in their phase-space kinetic theory, are actually the result of assuming zero mass.

Self-accelerating beam is a kind of light beam capable of self-bending in free space without any external potential, of which a typical one is the well-known Airy beam. Such a beam has gained great attention for its extraordinary properties, including nondiffracting, self-accelerating and self-healing, which may have versatile applications in the delivery and guiding of energy, information and objects using light, such as particle manipulation, micro-machining, optical routing, super-resolution imaging, etc. However, since Airy beam can only propagate along parabolic trajectory, which reduces the flexibility in practical applications, thus how to design accelerating beams propagating along arbitrary trajectory is still a crucial problem in this area. One scheme is to keep on finding other analytical solutions of the wave equation besides Airy beam, such as semi-Bessel accelerating beams, Mathius beams, and Weber beams, moving along circular, elliptical, or parabolic trajectories, but it becomes increasingly difficult to find out any more solutions. A more effective solution to this problem is based on the caustic method, which associates the predesigned trajectory with an optical caustics and then obtains the necessary initial field distribution by performing a light-ray analysis of the caustics. This method has been implemented in real space and Fourier space based on Fresnel diffraction integral and angular-spectrum integral, respectively. It has been found recently that they can be unified by constructing Wigner distribution function in phase space. Based on the caustic method, accelerating beams were constructed to propagate along arbitrary convex trajectories in two-dimensional space at first. With continuous development of this method, the types of accelerating beams available have been extending from convex trajectories to nonconvex trajectories, from two-dimensional trajectories to three-dimensional trajectories, and from one main lobe to multiple main lobes, which opens up more possibilities for emerging applications based on accelerating beams. In future, previous researches and applications based on Airy beams will certainly be generalized to all these new types of accelerating beams, and owing to the great flexibility in designing accelerating beams, more application scenarios may emerge in this process with huge development potential. Thus in this paper, we review the principle and progress of the caustic method in designing accelerating beams.

A semiclassical independent electron model for ion-atom collisions

Landau theory of Fermi liquids and nuclear fluid dynamics

New approach for solving the Wigner function from the Husimi function in the two-mode entangled state representation

The squeezed light field is a kind of important continuous variable quantum resource.It has wide applications in precision measurement and quantum information processing.Quantum storage is the foundations of quantum repeater and long distance quantum communication,and alkali metal atoms are an ideal quantum storage medium due to long ground state coherent time. With the rapid development of quantum storage technology in atomic medium,the preparation of the squeezed light which resonates with alkali metal atoms has become one of the research hotspots in the field of quantum information.In this paper,we report the generation of squeezed vacuum at 795 nm (resonant on the rubidium D1 transition line) by using an optical parametric oscillation based on a periodically poled KTiOPO4 crystal. The generated squeezed light field is detected by a balanced homodyne detector,and the squeezing of-3 dB and anti-squeezing of 5.8 dB are observed at a pump power of 45 mW.By using a maximum likelihood estimation,the density matrix of the squeezed light field is reconstructed.The time-domain signals from the balanced homodyne detector are collected to acquire the noise distribution of the squeezed light under different phase angles.The likelihood function is established for the measured quadrature components.An identity matrix is chosen as an initial density matrix,and the density matrix of the squeezed field is obtained through an iterative algorithm.The diagonal elements of the density matrix denote the photon number distribution,which includes not only even photon number states but also odd photon number states.The occurrence of odd photon number states mainly comes from the system losses and the imperfect quantum efficiency of detector.The Wigner function in phase space is calculated through the density matrix,and the maximum value of the Wigner function is 0.309.The standard deviation of the squeezed component is 64.4% of that of the vacuum state,corresponding to the squeezing degree of-3.8 dB.The standard deviation of the anti-squeezing component is 1.64 times that of the vacuum state,corresponding to the anti-squeezing degree of 4.3 dB.We theoretically calculate the photon number distribution and the Wigner function of the vacuum squeezed field,and compare the results obtained by theoretical calculation with those obtained by maximum likelihood reconstruction.The probability of vacuum state|0 obtained by maximum likelihood reconstruction is greater,and the probability of photon number state|n(n=1,2,) is smaller than the corresponding theoretical calculation results.From the theoretical calculation,the maximum value of Wigner function is 0.231,and the short axis and long axis of noise range deduced from the contours of the Wigner function are larger than the results from the maximum likelihood reconstruction.The possible reasons for the discrepancy are as follows. 1) The phase scanning is nonuniform during the measurement of the quadrature components.2) The low-frequency electronic noise is not completely filtered out in the datum acquisition process.3) The datum points of measured quadrature components are not enough.In conclusion,we produce a vacuum squeezed field of 795 nm,and obtain the photon number distribution and the Wigner function in phase space through maximum likelihood estimation and theoretical calculation,respectively.This work will provide an experimental basis for generating the Schrodinger cat state.