Quasiprobability distribution functions from fractional Fourier transforms

We show, in a formal way, how a class of complex quasiprobability distribution functions may be introduced by using the fractional Fourier transform. This leads to the Fresnel transform of a characteristic function instead of the usual Fourier transform. We end the manuscript by showing a way in which the distribution we are introducing may be reconstructed by using atom-field interactions.


Introduction
It has been already shown that quasiprobability distribution functions may be reconstructed by the measurement of atomic properties in ion-laser interactions [1] and two-level atoms interacting with quantized fields [2,3].Such measurements of the wave function are realized usually by measuring atomic observables, namely, the atomic inversion and polarization.
Ideal interactions, i.e., without taking into account an environment, have shown to lead to the reconstruction of the Wigner function [3] by taking advantage of its expression in terms of the parity operator.However, the interaction of a system with its environment [11] leads to s-prametrized quasiprobability distribution functions [12][13][14] where D(α) = exp(αâ † − α * â), with â and â † the annihilation and creation operators of the harmonic oscillator, respectively, is the Glauber displacement operator [15].The state D(α)|k = |α, k is a so-called displaced number state [16].Note that, in order to reconstruct a given quasiprobability function it is needed to do displace the system by an amplitude α and then measure the diagonal elements of the displaced density matrix.The parameter s defines different orderings and therefore different quasiprobability distribution functions (QDF).The Glauber-Sudarshan P -function [15,17] is given for s = 1, and is used to obtain averages of functions of normal ordered creation and annihilation operators; s = −1 gives the Husimi Q-function, used to obtain averages of functions of anti-normal ordered creation and annihilation operators, while s = 0 is used for the symmetric ordering and gives the Wigner function.
Equation (1) may be rewritten as that, by using the commutation properties under the symbol of trace, and the system in a pure state |ψ , may be casted into Recent studies have openned the possibility of measuring, instead of observables, non-Hermitian operators [18].It would be plausible then that such measurements could be related to complex quasiprobability distributions like the McCoy-Kirkwood-Rihaczek-Dirac distribution functions [5,6,8,19].
In this contribution we would like to introduce other kind of complex quasiprobabilities that, although they could be introduced simply by taking s as a complex number, we introduce them in a formal way by considering the fractional Fourier transform (FrFT) [20][21][22] of a signal.Then, by writing the Dirac-delta function in terms of its FrFT, we are able to write a general expression for complex quasiprobability distributions in terms of the Fresnel transform.Indeed, the representation of these complex quasiprobability distributions in terms of a Fresnel transform implies that they are solutions of a paraxial wave equation [3].Finally, by using an effective Hamiltonian for the atom-field interaction, we show how this quasiprobability distribution function may be reconstructed.

Fractional Fourier Transform
Up to a phase, the fractional Fourier Transform of a signal ψ(x) can be written by the following expression [20][21][22] F that may be expressed in terms of an integral transform as where Then, if we consider equation ( 6) as a propagator, Dirac's delta distribution function takes the form Now, if we apply the fractional Fourier transform to the Dirac delta function we obtain Then, applying the inverse fractional Fourier transform to equation ( 8) we obtain an alternative representation of the Dirac delta distribution function From the above equation it may be seen that there is a phase multiplying the usual integral representation of the Dirac delta function, that although could be omitted by using properties of the delta function, we keep in order to obtain a quasiprobability distribution function as a fractional Fourier (Fresnel) transform of the characteristic function.

Probability distribution in the phase space
We define J(q, p) a probability distribution in the phase space as then, using equation ( 9), the distribution J(q, p) may be rewritten as and because equation ( 11) takes the form that by using the equivalence may be casted into the expression The above quasiprobability distribution function is defined for a range of parameters α and β, however, for the sake of simplicity, we will consider the case cot α = − cot β = π.We may relate the quasiprobability distribution function J(q, p) to the Wigner function, by noting that, for cot α = − cot β = π, equation ( 15) has the form According to trace representation of Wigner function [14] we write the distribution J(q, p) as the Fresnel transform of the Wigner function It is easy to show that the quasiprobability distribution ( 18) can be normalized Therefore, for normalization reasons, the quasiprobability distribution is finally given in the form that, by applying the change of variables β = u/ √ 2 + iv/ √ 2 takes the form with α = q/ √ 2 + ip/ √ 2. From the above expression it is direct to show that the Wigner function and the function J(α) may be easily related by the differential relation The above quasiprobability function may be written as a trace by noting that that leads to the trace representation of J(q, p) Last equation allows us to show that J(q, p) is correctly normalized, for this we do the double integration where we have defined and with By replacing equation (29) into equation ( 28) we obtain Â = e −iθn , (30) that shows that equation ( 26) is correctly normalized 4. Kirkwood distribution and J(q, p) distribution Being the QDF J(q, p) and Kirkwood distributions complex functions we show now some differences between them.The Kirkwood distribution is defined as [8,19,23,24] du dv e iup−ivq e i uv 2 T r ρe iv q−iup , or an alternative way to write it as an expectation value [25] is

Number state
The Kirkwood K(q, p) and J(q, p) distributions for number state |n , are represented by the following equations and where, H n (x) and L n (x) are Hermite and Laguerre polynomials, respectively.

Superposition of two coherent states
Now, we consider a superposition of two coherent states as: where , such that the Kirkwood K(q, p) and the J(q, p) distributions for the superposition of two coherent states, |ψ ± , is given by and respectively.We plot both distribution in Figures 1 and 2. In both figures a more  and d) we see the distribution J(q, p), for the same number state, again, the real and imaginary parts, respectively.
uniform behaviour may be seen in the QDF J ± (q, p) than in the Kirkwood function.In fact, the real and imaginary parts of the distribution we have introduced here, look like Wigner function for number states (Fig. 1) and Scrhödinger cat states (Fig. 2).

Reconstruction of distribution J(α)
It is not difficult to show that the real part of QDF J(α) may be measured.This can be achieved by measuring the atomic polarization in the dispersive interaction between an atom and a quantized field [3], whose Hamiltonian reads Therefore, by measuring the polarizations σx and σy we are able to measure the QDF J(α).

Conclusions
We have introduced a set of parametrized (in terms of α and β) quasiprobability distribution functions, equation (15), by using the fractional Fourier transform.This has lead us to generalize QDF to Fresnel transforms of the characteristic function instead of their usual Fourier transforms.We have also shown how such QDF may be recosntructed in the dispersive atom-field interaction.We have also given a (differential) relation that allows the calculation of the newly introduced QDF from the Wigner function.

Figure 1 .
Figure 1.In figures a) and c) we can see the phase space distribution of the real and imaginary parts of the Kirkwood function for a number state |n = 3 .In figures b)and d) we see the distribution J(q, p), for the same number state, again, the real and imaginary parts, respectively.
with σz = |e e|−|g g|, the Pauli matrix corresponding to the atomic inversion operator, where |g and |e represent the ground and excited states of the two-level atom.The

Figure 2 . 1 √ 2 (π 2 − 4 2χ,
Figure 2. In figures a) andc) we can see the phase space distribution of the real and imaginary parts of the Kirkwood function for two superposition of coherent states |ψ + wiht q 1 = −q 2 = 4 and p 1 = p 2 = 0.In figures b) and d) we see the distribution J(q, p), again, the real and imaginary parts, respectively.