Correspondence between quantum-optical transform and classical-optical transform explored by developing Dirac’s symbolic method

Hong-yi Fan , Li-yun Hu

Front. Phys. ›› 2012, Vol. 7 ›› Issue (3) : 261 -310.

PDF (696KB)
Front. Phys. ›› 2012, Vol. 7 ›› Issue (3) :261 -310. DOI: 10.1007/s11467-011-0206-z
REVIEW ARTICLE
Correspondence between quantum-optical transform and classical-optical transform explored by developing Dirac’s symbolic method
Author information +
History +
PDF (696KB)

Abstract

By virtue of the new technique of performing integration over Dirac’s ket–bra operators, we explore quantum optical version of classical optical transformations such as optical Fresnel transform, Hankel transform, fractional Fourier transform, Wigner transform, wavelet transform and Fresnel–Hadmard combinatorial transform etc. In this way one may gain benefit for developing classical optics theory from the research in quantum optics, or vice-versa. We cannot only find some new quantum mechanical unitary operators which correspond to the known optical transformations, deriving a new theorem for calculating quantum tomogram of density operators, but also can reveal some new classical optical transformations. For examples, we find the generalized Fresnel operator (GFO) to correspond to the generalized Fresnel transform (GFT) in classical optics. We derive GFO’s normal product form and its canonical coherent state representation and find that GFO is the loyal representation of symplectic group multiplication rule. We show that GFT is just the transformation matrix element of GFO in the coordinate representation such that two successive GFTs is still a GFT. The ABCD rule of the Gaussian beam propagation is directly demonstrated in the context of quantum optics. Especially, the introduction of quantum mechanical entangled state representations opens up a new area in finding new classical optical transformations. The complex wavelet transform and the condition of mother wavelet are studied in the context of quantum optics too. Throughout our discussions, the coherent state, the entangled state representation of the two-mode squeezing operators and the technique of integration within an ordered product (IWOP) of operators are fully used. All these have confirmed Dirac’s assertion: “…for a quantum dynamic system that has a classical analogue, unitary transformation in the quantum theory is the analogue of contact transformation in the classical theory”.

Graphical abstract

Keywords

Dirac’s symbolic method / IWOP technique / entangled state of continuum variables / entangled Fresnel transform / Collins formula / Generalized Fresnel operator / complex wavelet transform / complex Wigner transform / complex fractional Fourier transform / symplectic wavelet transform / entangled symplectic wavelet transform / Symplectic-dilation mixed wavelet transform / fractional Radon transform / new eigenmodes of fractional Fourier transform

Cite this article

Download citation ▾
Hong-yi Fan, Li-yun Hu. Correspondence between quantum-optical transform and classical-optical transform explored by developing Dirac’s symbolic method. Front. Phys., 2012, 7 (3) : 261-310 DOI:10.1007/s11467-011-0206-z

登录浏览全文

4963

注册一个新账户 忘记密码

1 Introduction

The history of quantum mechanics records that from the very beginning the founders of the quantum theory realized that there might exist formal connection between classical optics and quantum mechanics. For example, Schrödinger considered that classical dynamics of a point particle should be the “geometrical optics” approximation of a linear wave equation, in the same way as ray optics is a limiting approximation of wave optics; Schrödinger also searched for some quantum mechanical state which behaves like a classical ‘particle’, and this state was later recognized as the coherent state [1–6], which plays an essential role in quantum optics theory and laser physics; As Dirac wrote in his famous book “Principles of Quantum Mechanics” [7]: “. . . for a quantum dynamic system that has a classical analogue, unitary transformation in the quantum theory is the analogue of contact transformation in the classical theory”. According to Dirac, there should exist a formal correspondence between quantum optics unitary-transform operators and classical optics transformations. Indeed, in the last century physicists also found some rigorous mathematical analogies between classical optics and quantum mechanics, i.e., the similarity between the optical Helmholtz equation and the time-independent Schrödinger equation; Since 1960s, the advent of a laser and the appearance of coherent state theory of radiation field [1–5, 8], quantum optics has experienced rapid development and achieved great success in revealing and explaining the quantum mechanical features of optical field and non-classical behavior (for instance, Hanbury–Brown–Twiss effect, photon antibunching, squeezing, sub-Poissonian photon statistics) of photons in various photon–atom interactions [9–11]. The relationship between classical and quantum coherence has been discussed in the book of Mandel and Wolf [9]; The Hermite–Gauss or the Laguerre–Gauss modes of a laser beam are described using the bosonic operator algebra by Nienhuis and Allen [12]. In addition, displaced light beams refracted by lenses according to the law of geometrical optics were found to be the paraxial optics analog of a coherent state. Besides, phase space correspondence between classical optics and quantum mechanics, say for example, the Wigner function theory, is inspected in the literature [13].

On the other hand, classical optics, which tackles vast majority of physical-optics experiments and is based on Maxwell’s equations, has never ceased its own evolving steps, physicists have endeavored to develop various optical transforms in light propagation through lens systems and various continuous media. The two research fields, quantum optics and classical optics, have their own physical objects and conceptions. From the point of view of mathematics, classical optics is framed in the group transform and associated representations on appropriate function space, while quantum optics deals with operators and state vectors, and their overlap seems little at first glance. It seems to us that if one wants to further relate them to each other, one needs some new theoretical method to “bridge” them. For example, what is the quantum mechanical unitary operator corresponding to the Fresnel transform in Fourier optics? Is there any so-called Fresnel operator as the image of classical generalized Fresnel transform? Since generalized Fresnel transforms are very popularly used in optical instrument design and optical propagation through lenses and various media, it is worthy to study these transforms in the context of quantum optics theory, especially based on coherent state, squeezed state [14–17] and the newly invented entangled state theory [18–27].

Fortunately, the recently invented technique of integration within an ordered product (IWOP) of operators [28, 29], which present a novel method for performing integration over Dirac’s ket-bra operators, is of great aid to studying quantum optical version of classical optical transformations. Using the IWOP technique one may gain benefit for classical optics from quantum optics’ research, or vice-versa.

Our present review is arranged as follows: in Section 2 we briefly recall the classical diffraction theory [30, 31]. This prepares for later sections in which we shall show that most frequently employed classical optical transforms have their counterparts in quantum optics theory. In Section 3 we introduce the IWOP technique and demonstrate that the completeness relation of fundamental quantum mechanical representations can be recast into normally ordered Gaussian operator form. Using the IWOP technique we can directly perform the asymmetric ket–bra integration μ−1/2∫−∞∞dq|q/μ⟩⟨q| in the coordinate representation, which leads to the normally ordered single-mode squeezing operator, which seems to be a direct way of understanding the squeezing mechanism as a mapping from the classical scaling q→q/μ. In Section 4 with the help of IWOP technique and based on the concept of quantum entanglement of Einstein–Podolsky–Rosen [32] we construct two mutually conjugate entangled states of continuum variables, |η⟩ versus |ξ⟩, and their deduced entangled states (or named correlated-amplitude — number-difference entangled states), which are all qualified to make up quantum mechanical representations. It is remarkable that using the IWOP technique in performing the asymmetric ket–bra integration μ−1∫d2η|η/μ⟩⟨η| leads to the normally ordered two-mode squeezing operator, which implies that the two-mode squeezed state is simultaneously an entangled state. We point out that the entangled state |η⟩ also embodies entanglement in the aspect of correlative amplitude and the phase. We felt encouraged from that the overlap between two mutually conjugate deduced entangled states is just the Bessel function — the optical Hankel transform kernel [33], which again shows that the new representations in the context of physics theory match beautiful mathematical formalism exactly. We then employ the deduced entangled states to derive quantum optical version of classical circular harmonic correlation. Section 5 is devoted to finding a quantum operator which corresponds to the optical Fresnel transform, with use of the coherent state representation and by projecting the classical sympletic transform z→sz−rz*(|s|2−|r|2=1) in phase space onto the quantum mechanical Hilbert space. We are able to recognize which operator is the single-mode Fresnel operator (FO). It turns out that the one-dimensional (1D) optical Fresnel transform is just the matrix element of the Fresnel operator F in the coordinate eigenstates. Besides, the coherent state projection operator representation of FO constitutes a loyal realization of symplectic group, which coincides with the fact that two successive optical Fresnel transforms make up a new Fresnel transform. Then in Section 6, based on the coherent state projection representation of FO, we prove ABCD rule for optical propagation in the context of quantum optics. In Section 7 the quadratic operator form of FO is also presented and the four fundamental optical operators are derived by decomposing the FO. In Section 8 we discuss how to apply the Fresnel operator to quantum tomography theory. By introducing the Fresnel quadrature phase FXF†=XF, we point out that Wigner operator’s Radon transformation is just the pure state projection operator |x⟩s,rs,r⟨x|, where |x⟩s,r=F|x⟩ and |x⟩ is the position eigenstate, so the probability distribution for the Fresnel quadrature phase is the Radon transform of the Wigner function. Moreover, the tomogram of quantum state |ψ⟩ is just the squared modulus of the wave function ⟨x|ψ⟩s,r. This new relation between quantum tomography and optical Fresnel transform may provide experimentalists to figure out new approach for testing tomography. In addition, we propose another new theorem for calculating tomogram, i.e., the tomogram of a density operator ρ is equal to the marginal integration of the classical Weyl correspondence function of F†ρF. In Section 9, by virtue of the coherent state and IWOP method we propose two-mode generalized Fresnel operator (GFO). In this case we employ the entangled state representation to relate the 2-mode GFO to classical transforms, since the 2-mode GFO is not simply the direct product of two 1-mode GFOs. The corresponding quantum optics ABCD rule for two-mode case is also proved. The 2-mode GFO can also be expressed in quadratic operators form in entangled way. The relation between optical Fresnel transform and quantum tomography in the two-mode case is also revealed. In Section 10 we propose a kind of integration transformation, ∫∫−∞∞dpdqπe2i(p−x)(q−y)h(p,q)≡f(x,y), which is invertible and obeys Parseval theorem. Remarkably, it can convert chirplet function to the kernel of fractional Fourier transform (FrFT). This transformation can also serve for solving some operator ordering problems. In Section 11 we employ the entangled state representation to introduce the complex FrFT (CFrFT), which is not the direct product of two independent 1D FrFT transform. The eigenmodes on CFrFT is derived. New eigen-modes for light propagation in graded-index medium and the fractional Hankel transform are presented. The Wigner transform theory is extended to the complex form and its relation to CFrFT is shown; The integration transformation in Section 10 is also extended to the entangled case. In Section 12 we shall treat the adaption problem of Collins diffraction formula to the CFrFT with the use of 2-mode (3 parameters) squeezing operator and in the entangled state representation of continuous variables, in so doing the quantum mechanical version of associated theory of classical diffraction and classical CFrFT is obtained, which connects classical optics and quantum optics in this aspect. In Section 13 we discuss quantum optical version of classical wavelet transforms (WTs), including how to recast the condition of mother wavelet into the context of quantum optics. Some properties, such as Parseval theorem, Inversion formula, and orthogonal property, the relation between WT and Wigner–Husimi distribution function are also discussed. In Section 14, we generalize the usual wavelet transform to symplectic wavelet transformation (SWT) by using the coherent state representation and making transformation z→s(z−κ)−r(z*−κ*) (|s|2−|r|2=1) in phase space. The relation between SWT and optical Fresnel transformation is revealed. Then the SWT is extended to the entangled case by mapping the classical mixed transformation (z,z′)→(sz+rz′⁣*,sz′+rz*) in 2-mode coherent state |z,z′⟩ representation. At the end of this section, we introduce a new symplectic-dilation mixed WT by employing a new entangled-coherent state representation |α,x⟩. The corresponding classical optical transform is also presented. In the last section, we introduce the Fresnel–Hadamard combinatorial operator by virtue of the IWOP technique and |α,x⟩. This unitary operator plays the role of both Fresnel transformation for mode a1−a22 and Hadamard transformation for mode a1+a22, respectively, and the two transformations are combinatorial. All these sections are used to prove the existence of a one-to-one correspondence between quantum optical operators that transform state vectors in Hilbert space and the classical optical transforms that change the distribution of optical field.

2 Some typical classical optical transformations

Here we briefly review some typical optical transforms based on light diffraction theory. These transformations, as one can see in later sections, are just the correspondence of some representation transformations between certain quantum mechanical states of which some are newly constructed.

It was Huygens who gave a first illustrative explanation to wave theory by proposing every point in the propagating space as a sub-excitation source of a new sub-wave. An intuitive theory mathematically supporting Huygens’ principle is the scalar diffraction approximation, so named because optical fields (electromagnetic fields) actually are vector fields, whereby the theory is valid approximately. This theory is based on the superposition of the combined radiation field of multiple re-emission sources initiated by Huygens. Light diffraction phenomena has played an important role in the development of the wave theory of light, and now underlies the Fourier optics and information optics. The formulation of a diffraction problem essentially considers an incident free-space wave whose propagation is interrupted by an obstacle or mask which changes the phase and/or amplitude of the wave locally by a well determined factor [34]. A more rigorous, but still in the scheme of scalar wave, derivation has been given by Kirchhoff who reformulated the diffraction problem as a boundary-value problem, which essentially justifies the use of Huygens principle. The Fresnel–Kirchhoff (or Rayleigh–Sommerfeld) diffraction formula is practically reduced to the Fresnel integral formula in paraxial and far-field approximation [30, 31] that reads

(1)U2(x2,y2)=exp⁡(ikz)iλz∫∫−∞∞U1(x1,y1)×exp⁡{ik2z[(x2−x1)2+(y2−y1)2]}×dx1dy1

where U1(x1,y1) is the optical distribution of a 2D light source and U2(x2,y2) is its image on the observation plane, λ is the optical wavelength, k=2πλ is the wave number in the vacuum and z is the propagation distance. When z2≫k2(x12+y12)max is satisfied, Eq. (1) is reduced to

(2)U2(x2,y2)=exp⁡(ikz)exp⁡[ik2z(x12+y12)]iλz×∫∫−∞∞U1(x1,y1)×exp⁡[−i2πλz(x1x2+y1y2)]dx1dy1

which is named the Fraunhofer diffraction formula.

The Fresnel integral is closely related to the fractional Fourier transform (FrFT). Actually, it has been proved that the Fresnel transform can be interpreted as a scaled FrFT with a residual phase curvature [35]. The FrFT is a very useful tool in Fourier optics and information optics. This concept was firstly introduced in 1980 by Namias [36], but had not attracted enough attention until FrFT was defined physically, based on propagation in quadratic graded-index media (GRIN media). Mendlovic and Ozaktas [37–39] defined the αth FrFT as follows: Let the original function be input from one side of quadratic GRIN medium, at z=0. Then, the light distribution observed at the plane z=z0 corresponds to the α equal to the (z0/L)-th FrFT of the input fraction, where L≡(π/2)(n1/n2)1/2 is a characteristic distance. The FrFT can also be implemented by lenses. Another approach for introducing FrFT was made by Lohmann who pointed out the algorithmic isomorphism among image rotation, rotation of the Wigner distribution function [42], and fractional Fourier transforming [40]. Lohmann proposed the FrFT as the transform performed on a function that leads to a rotation with an angle of the associated Wigner distribution function. In this sense, the FrFT bridges the gap between the classical optics and the optical Wigner distribution theory. Recently, the FrFT has been paid more and more attention within different contexts of both mathematics and physics [36–41]. The FrFT is defined as

(3)Fα[U1](x2,y2)=ei(1−α)π22sin⁡(π2α)exp⁡[−i(x22+y22)2tan⁡(π2α)]×∫∫−∞∞dx1dy1πexp[−i(x12+y12)2tan⁡(π2α)]×exp⁡[i(x2x1+y2y1)sin⁡(π2α)]U1(x1,y1)

We can see that F0 is the identity operator and Fπ/2 is just the Fourier transform. The most important property of FrFT is that Fα obeys the semigroup property, i.e., two successive FrFTs of order α and β makes up the FrFT of order α+β. A more general form describing the light propagation in an optical system characterized by the [A,B;C,D] ray transfer matrix is the Collins diffraction integral formula [43],

(4)U2(x2,y2)=kexp⁡(ikz)2πiB∫∫−∞∞dx1dy1U1(x1,y1)×exp{ik2B[A(x12+y12)−2(x1x2+y1y2)+D(x22+y22)]}dx1dy1

where AD−BC=1 if the system is lossless. One can easily find the similarity between Collins formula and the FrFT by some scaling transform and relating the [A,B,C,D] matrix to α in the FrFT [44]. Note that M=(ABCD) is a ray transfer matrix describing optical systems belonging to the unimodular symplectic group. When treating the light propagation in optical elements in near-axis approximation, matrices M representing linear transformations are a convenient mathematical tool for calculating the fundamental properties of optical systems, which is the origin of the name of matrix optics. In cylindrical coordinates the Collins formula is expressed as [43, 44]

(5)U2(r2,φ)=iλB∫0∞∫02πexp{−iπλB[Ar12+Dr22−2r1r2cos(θ−φ)]}U1(r1,θ)r1dr1dθ

where x1=r1cos⁡θ,y1=r1sin⁡θ,x2=r2cos⁡φ and y2=r2sin⁡φ. When U1(r1,θ) has rotational symmetry

(6)U1(r1,θ)=u1(r1)exp⁡(imθ)U2(r2,φ)=u2(r2)exp⁡(imφ)

then Eq. (5) becomes

(7)u2(r2)=2πλBexp⁡[i(1+m)π2]×∫0∞exp[−iπλB(Ar12+Dr22)]×Jm(2πr1r2λB)u1(r1)r1dr1

where we have used the m-order Bessel function

(8)Jm(x)=12π∫02πexp⁡[ixcos⁡θ+im(θ−π2)]dθ

When A=0, Eq. (7) reduces to the standard Hankel transform (up to a phase factor)

(9)u2(r2)→2πλB∫0∞Jm(2πr1r2λB)u1(r1)r1dr1

The compact form of 1D Collins formula is

(10)g(x2)=∫−∞∞KM(x2,x1)f(x1)dx1

where the transform kernel is

(11)KM(x2,x1)=12πiBexp⁡[i2B(Ax12−2x2x1+Dx22)]

M is the parameter matrix [A,B,C,D]. Eq. (11) is called generalized Fresnel transform [46–50]. In the following sections we will show how we find the quantum optical counterpart for those transformations of classical optics. For this purpose in the next chapter we introduce the IWOP technique to demonstrate how Dirac’s symbolic method can be developed and be applied to quantum optics theory. Also, we briefly review some properties of the entangled state [18–27] and reveal the connection between the mutual transform generated by these entangled states and the Hankel transform in classical optics.

3 The IWOP technique and two mutually conjugate entangled states

3.1 The IWOP technique

The history of mathematics tells us that whenever there appears a new important mathematical symbol, there coexists certain operational rules for it, the quantum mechanical operators in ket-bra projective form (the core of Dirac’s symbolic method) also need their own operational rules. The terminology “symbolic method” was first shown in the preface of Dirac’s book “The Principle of Quantum Mechanics”: “The symbolic method, which deals directly in an abstract way with the quantities of fundamental importance⋯, however, seems to go more deeply into the nature of things. It enables one to express the physical law in a neat and concise way, and will probably be increasingly used in the future as it becomes better understood and its own special mathematics gets developed.” [7] Then two questions naturally arise: How to better understand the symbolic method? How to develop Dirac’s symbolic method, especially its mathematics? We noticed that Newton–Leibniz integration rule only applies to commuting functions of continuum variables, while operators made of Dirac’s symbols (ket versus bra, e.g., |q/μ⟩⟨q| of continuous parameter q) in quantum mechanics are usually not commutative. Therefore integrations over the operators of type |⟩⟨| (where ket- and bra-state vectors need not to be Hermitian-conjugate to each other) can not be directly performed by the Newton–Leibniz rule. Thus we invented an innovative technique of integration within an ordered product (IWOP) of operators that made the integration of non-commutative operators possible. The core of IWOP technique is to arrange non-commutable quantum operators within an ordered product (say, normal ordering) in a way that they become commutable. In this sense the gap between q-numbers and c-numbers is “narrowed”. However, the nature of operators which are within : : is not changed. They are still q-numbers, not c-numbers. After the integration over c-numbers within ordered product is performed, we can get rid of the normal ordering symbol after putting the integration result in normal ordering [51]. The IWOP technique thus bridges this mathematical gap between classical mechanics and quantum mechanics, and further reveals the beauty and elegance of Dirac’s symbolic method and transformation theory. This technique develops symbolic method significantly, i.e., makes Dirac’s representation theory and the transformation theory more plentiful, and consequently to be better understood. The beauty and elegance of Dirac’s symbolic method are further revealed. Various applications of the IWOP technique, including constructing the entangled state, developing the nonlinear coherent state theory, Wigner function theory, etc., are found; many new unitary operators and operator-identities as well as new quantum mechanical representations can be derived too, which are partly summarized in the review articles [22–24].

We begin with listing some properties of normal product of operators which means all the bosonic creation operators a† are standing on the left of annihilation operators a in a monomial of a† and a.

(i) The order of Bose operators a and a† within a normally ordered product can be permuted. That is to say, even though [a,a†]=1, we can have: aa†:=:a†a:=a†a, where : : denotes normal ordering.

(ii) c-numbers can be taken out of the symbol : : as one wishes.

(iii) The symbol : : which is within another symbol : : can be deleted.

(iv) The vacuum projection operator |0⟩⟨0| has the normal product form

(12)|0⟩⟨0|=:e−a†a:

(v) A normally ordered product can be integrated or differentiated with respect to a c-number provided the integration is convergent.

3.2 The IWOP technique for deriving normally ordered Gaussian form of the completeness relations of fundamental quantum mechanical representations

As an application of IWOP, (in the following, unless particularly mentioned, we take ℏ=ω=m=1 for convenience.), using the Fock representation of the coordinate eigenvector Q|q⟩=q|q⟩,(Q=(a+a†)/2)

(13)|q⟩=π−1/4e−q22+2qa†−a†22|0⟩

we perform the integration below:

(14)S1≡∫−∞∞dqμ|qμ⟩⟨q|=∫−∞∞dqπμe−q22μ2+2qμa†−a†22|0⟩⟨0|e−q22+2qa−a22

Substituting Eq. (12) into Eq. (14), we see

(15)S1=∫−∞∞dqπμe−q22μ2+2qμa†−a†22:e−a†a:e−q22+2qa−a22

Note that on the left of : e−a†a: are all creation operators, while on its right are all annihilation operators, so the whole integral is in normal ordering, thus using property i) we have

(16)S1=∫−∞∞dqπμ:e−q22(1+1μ2)+2q(a†μ+a)−12(a+a†)2:

As a commutes with a† within : :, a† and a can be considered as if they were parameters while the integration is being performed. Therefore, by setting μ=eλ, sech λ=2μ1+μ2, tanh λ=μ2+1μ2−1, we are able to perform the integration and obtain

(17)S1=(sechλ)1/2e−a†22tanh⁡λ:e(sechλ−1)a†a:ea22tanh⁡λ

which is just the single-mode squeezing operator in normal ordering appearing in many references. It is worth mentioning that we have not used the SU (1,1) Lie algebra method in the derivation. The integration automatically arranges the squeezing operator in normal ordering. With

(18)eλa†a=:exp[(eλ−1)a†a]:

Eq. (17) becomes

(19)∫−∞∞dqμ|qμ⟩⟨q|=e−a†22tanh⁡λe(a†a+12)ln⁡sechλea22tanh⁡λ

This shows that the classical dilation q→qμ is mapped into the normally ordered squeezing operator manifestly. It also exhibits that the fundamental representation theory can be formulated in a less abstract way, as we can now directly perform the integral over ket–bra projection operators. Moreover, the IWOP technique can be employed to perform many complicated integrations for ket–bra projection operators.

There is a deep ditch between quantum mechanical operators (q-numbers) theory and classical numbers (c-numbers) theory. The IWOP technique arranges non-commutable operators within an ordered product symbol in a way that they become commutable. In this sense the “ditch” between q-numbers and c-numbers is “shoaled”. However, the nature of operators are not changed, they are still q-numbers, not c-numbers. After the integration over c-numbers within ordered product is performed, we can finally get rid of the normal ordering symbol by using Eq. (18).

When μ=1, Eq. (19) becomes

(20)∫−∞∞dq|q⟩⟨q|=∫−∞∞dqπ:e−(q−Q)2:=1

a real simple Gaussian integration! This immediately leads us to put the completeness relation of the momentum representation into the normally ordered Gaussian form

(21)∫−∞∞dp|p⟩⟨p|=∫−∞∞dpπ:e−(p−P)2:=1

where P=(a−a†)/(i2),and |p⟩ is the momentum eigenvector P|p⟩=p|p⟩,

(22)|p⟩=π−14exp⁡(−12p2+i2pa†+12a†2)|0⟩

In addition, we should notice that |q⟩ and |p⟩ are related by the Fourier transform (FT), i.e., ⟨p|q⟩=12πexp⁡(−iqp). The integral kernel of the Fraunhofer diffraction formula in 1D is such an FT, so FT in classical optics has its correspondence in quantum mechanical representations’ transform. This enlightens us that in order to find more general analogy between unitary operators in quantum optics and transformations in classical optics we should construct new representations for quantum optics theory, and these are the bi-partite entangled state and many-particle entangled state. These ideal states can be implemented by optical devices and optical network [52]. In the following we will focus on the bi-partite entangled state.

3.3 Single-mode Wigner operator

When we combine Eqs. (20) and (21) we can obtain

(23)π−1:e−(q−Q)2−(p−P)2:≡Δ(q,p)

which is just the normally ordered Wigner operator since its marginal integration gives |q⟩⟨q| and |p⟩⟨p| respectively, i.e.,

(24)∫−∞∞dqΔ(q,p)=1π:e−(p−P)2:=|p⟩⟨p|

(25)∫−∞∞dpΔ(q,p)=1π:e−(q−Q)2:=|q⟩⟨q|

Thus the Wigner function of quantum state ρ can be calculated as W(q,p)=Tr[ρΔ(q,p)]. On the other hand, the Wigner operator (23) can be recast into the coherent state representation,

(26)Δ(q,p)→Δ(α,α*)=∫d2zπ2|α+z⟩⟨α−z|eαz*−α*z

where |z⟩ is a coherent state. In fact, using the IWOP technique we can obtain

(27)Δ(α,α*)=1π:exp⁡{−2(a−α)(a†−α*)}:

which is the same as Eq. (23).

3.4 Entangled state |η⟩and its Fourier transform in complex form

The concept of quantum entanglement was first employed by Einstein, Rosen and Poldosky (EPR) to challenge that quantum mechanics is incomplete when they observed that two particles’ relative position Q1−Q2 and the total momentum P1+P2 are commutable. Hinted by EPR, the bipartite entangled state |η⟩ is introduced as [53–55]

(28)|η⟩=exp⁡(−12|η|2+ηa1†−η*a2†+a1†a2†)|00⟩

|η=η1+iη2⟩ is the common eigenstate of relative coordinate Q1−Q2 and the total momentum P1+P2,

(29)(Q1−Q2)|η⟩=2η1|η⟩,(P1+P2)|η⟩=2η2|η⟩

where Qi=(aj+aj†)/2,Pj=(aj−aj†)/(i2),j=1,2. Using the IWOP technique, we can immediately prove that |η⟩ possesses the completeness relation

(30)∫d2ηπ|η⟩⟨η|=1,d2η=dη1dη2

and orthonormal relation

(31)⟨η|η′⟩=πδ(η1−η1′)δ(η2−η2′)

The Schmidt decomposition of |η⟩ is

(32)|η⟩=e−iη2η1∫−∞∞dq|q⟩1⊗|q−2η1⟩2ei2η2q

The |η⟩ state can also be Schmidt-decomposed in momentum eigenvector space as

(33)|η⟩=eiη1η2∫−∞∞dp|p⟩1⊗|2η2−p⟩2e−i2η1p

The |η⟩ is physically appealing in quantum optics theory, because the two-mode squeezing operator has its natural representation on ⟨η| basis [55]

(34)∫d2ηπμ|η/μ⟩⟨η|=ea1†a2†tanh⁡λe(a1†a1+a2†a2+1)ln⁡sechλ×e−a1a2tanh⁡λ≡S2,   μ=eλ

The proof of Eq. (34) is proceeded by virtue of the IWOP technique. So the necessity of introducing |η⟩ into quantum optics is clear. S2 squeezes |η⟩ in the manifest way

(35)S2|η⟩=1μ|η/μ⟩,μ=eλ

and the two-mode squeezed state itself is an entangled state which entangles the idle mode and signal mode as an outcome of a parametric-down conversion process [56–58].

We can also introduce the conjugate state of |η⟩ [59],

(36)|ξ⟩=exp⁡(−12|ξ|2+ξa1†+ξ*a2†−a1†a2†)|00⟩

which obeys the eigen-equations

(37)(Q1+Q2)|ξ⟩=2ξ1|ξ⟩,(P1−P2)|ξ⟩=2ξ2|ξ⟩

where ξ=ξ1+iξ2. Because [(Q1−Q2),(P1−P2)]=2i, we name the conjugacy between |ξ⟩ and |η⟩. The completeness and orthonormal relations of |ξ⟩ are

(38)∫d2ξπ|ξ⟩⟨ξ|=1,d2ξ=dξ1dξ2

(39)⟨ξ|ξ′⟩=πδ(ξ1−ξ1′)δ(ξ2−ξ2′)

respectively. |η⟩ and |ξ⟩ can be related to each other by

(40)⟨η|ξ⟩=12exp⁡(ξη*−ξ*η2)

Since ξ*η−ξη* is a pure imaginary number, and Eq. (40) is the Fourier transform kernel in complex form (or named entangled Fourier transform. This concept should also be extended to multipartite entangled states.). It will be shown in later sections that departing from entangled states |η⟩ and |ξ⟩, the generalized Fresnel operator would lead to a new entangled Fresnel transform in classical optics.

3.5 Two-mode Wigner operator in the |η⟩ representation

Combining Eqs. (30) and (38) we can construct the following operator:

(41)1π2:e−[σ*−(a1†−a2)][σ−(a1−a2†)]−[γ*−(a1†+a2)][γ−(a1+a2†)]:=Δ(α,α*)⊗Δ(β,β*)≡Δ(σ,γ)

where

(42)σ=α−β*,γ=α+β*

Eq. (41) is just equal to the direct product of two single-mode Wigner operators. It is convenient to express the Wigner operator in the |η⟩ representation as [60]

(43)Δ(σ,γ)=∫d2ηπ3|σ−η⟩⟨σ+η|eηγ*−η*γ

For two-mode correlated system, it prefers to using Δ(σ,γ) to calculate quantum states’ Wigner function. For example, with ⟨η|00⟩=exp⁡(−|η|2/2), the two-mode squeezed states’ Wigner function is

(44)⟨00|S2†(μ)Δ(σ,γ)S2(μ)|00⟩=π−2exp⁡(−μ2|σ|2−|γ|2/μ2)

4 Two deduced entangled state representations and Hankel transform

4.1 Deduced entangled states

Starting from the entangled state |η=reiθ⟩ and introducing an integer m, we can deduce new states [61],

(45)|m,r⟩=12π∫02πdθ|η=reiθ⟩e−imθ

which is worthy of attention because when we operate the number-difference operator, D≡a1†a1−a2†a2 on |η⟩, using Eq. (28) we see

(46)D|η⟩=(ηa1†+η*a2†)|η⟩=−i∂∂θ|η⟩,η=|η|eiθ

so the number-difference operator corresponds to a differential operation i∂∂θ in the ⟨η| representation, which is a remarkable property of ⟨η|. It then follows that

(47)D|m,r⟩=∫02πdθ2πe−imθ(−i∂∂θ|η=reiθ⟩)=m|m,r⟩

On the other hand, by defining K≡(a1−a2†)(a1†−a2), we see [D,K]=0, and |m,r⟩ is its eigenstate,

(48)K|m,r⟩=r2|m,r⟩

where K is named correlated-amplitude operator since K|η⟩=r2|η⟩. Thus we name |m,r⟩ correlated-amplitude— number-difference entangled states. It is not difficult to prove the completeness and orthonormal property of |m,r⟩,

(49)∑m=−∞∞∫0∞d(r2)|m,r⟩⟨m,r|=1

(50)⟨m,r|m′,r′′⟩=δm,m′12rδ(r−r′)

On the other hand, from |ξ⟩ we can derive another state

(51)|s,r′⟩=12π∫02πdφ|ξ=r′eiφ⟩e−isφ

which satisfies

(52)D|ξ⟩=(a1†ξ−a2†ξ*)|ξ⟩=−i∂∂φ|ξ=r′eiφ⟩

So D in ⟨ξ=r′eiφ| representation is equal to i∂∂φ. Consequently,

(53)D|s,r′⟩=∫02πdθ2πe−isθ(−i∂∂θ|ξ=r′eiθ⟩)=s|s,r′⟩

Note [D,(a1†+a2)(a1+a2†)]=0 and

(54)(a1†+a2)(a1+a2†)|s,r′⟩=r′⁣2|s,r′⟩

|s,r′⟩ is qualified to be a new representation since

(55)∑s=−∞∞∫0∞d(r′⁣2)|s,r′⟩⟨s,r′|=1⟨s,r′|s′′,r′′⟩=δs,s′12r′δ(r′−r′′)

4.2 Hankel transform between two deduced entangled state representations

Since |ξ⟩ and |η⟩ are mutual conjugate, |s,r′⟩ is the conjugate state of |m,r⟩. From the definition of |m,r⟩ and |s,r′⟩ and Eq. (40) we calculate the overlap [33]

(56)⟨s,r′|q,r⟩=12∑l=−∞∞δl,mδl,sJl(rr′)=12δs,mJs(rr′)

where we have identified the generating function of the s-order Bessel function Jl,eixsin⁡t=∑l=−∞∞Jl(x)eilt, and

(57)Jl(x)=∑k=0∞(−1)lk!(l+k)!(x2)l+2k

Eq. (56) is remarkable, because Js(rr′′) is just the integral kernel of Hankel transform. In fact, if we define

(58)⟨m,r|g⟩≡g(m,r),⟨s,r′|g⟩≡G(s,r′)

and use Eq. (49) as well as Eq. (56), we obtain

(59)G(s,r′)=12∫0∞d(r2)Js(rr′′)g(s,r)≡H[g(s,r)]

which is just the Hankel transform of g(m,r) [or it can be regarded as a simplified form of the Collins formula in cylindrical coordinate, see Eq. (9)]. The inverse transform of Eq. (59) is

(60)g(m,r)=12∫0∞d(r′⁣2)Jq(rr′)G(m,r′)≡H−1[G(m,r′)]

Now we know that the quantum optical image of classical Hankel transform just corresponds to the representation transformation between two mutually conjugate entangled states ⟨s,r′| and |m,r⟩, which is like the case that the Fourier transform kernel is just the matrix element between the coordinate state and the momentum state, a wonderful result unnoticed before. Therefore the bipartite entangled state representations’ transforms, which can lead us to the Hankel transform that was proposed first in classical optics, can find their way back in quantum optics.

4.3 Quantum optical version of classical circular harmonic correlation

From Eq. (45) we can see that its reciprocal relation is the circular harmonic expansion,

(61)|η=reiθ⟩=∑m=−∞∞|m,r⟩eimθ

or the correlated-amplitude — number-difference entangled state |m,r⟩ can be considered as circular harmonic decomposition of |η=reiθ⟩. Let g(r,θ), a general 2D function expressed in polar coordinates, be periodic in the variable θ, it can be looked as the wavefunction of the state vector |g⟩ in the ⟨η=reiθ| representation

(62)g(r,θ)=⟨η=reiθ|g⟩

Using Eq. (61) we have

(63)g(r,θ)=∑m=−∞∞gm(r)e−imθ,gm(r)=⟨m,r|g⟩

where gm(r) is the wavefunction of |g⟩ in ⟨m,r| representation. We see from Eq. (46)

(64)e−iα(a1†a1−a2†a2)|η=reiθ⟩=e−α∂∂θ|η=reiθ⟩=|η=rei(θ−α)⟩

so e−iα(a1†a1−a2†a2) behaves a rotation operator in |η⟩ representation, the expectation value of e−iα(a1†a1−a2†a2) in |g⟩ is

(65)π⟨g|e−iα(a1†a1−a2†a2)|g⟩=π⟨g|∫d2ηπe−iα(a1†a1−a2†a2)|η⟩⟨η|g⟩=∫0∞rdr∫02πg*(r,θ−α)g(r,θ)dθ≡Rα

which is just the cross-correlation between g(r,θ) and an angularly rotated version of the same function, g*(r,θ−α). On the other hand, using Eq. (49) we have

(66)|g⟩=∑m=−∞∞∫0∞d(r2)|m,r⟩gm(r)

Substituting Eq. (66) into Eq. (65) and using the eigenvector equation (47) as well as Eq. (50), we obtain

(67)Rα=2π∑m=−∞∞e−imα∫0∞r|gm(r)|2dr

from which we see that each of the circular harmonic components of the crosscorrelation undergoes a different phase shift −mα, so Rα is not rotation invariant. However, when we consider only one harmonic component

(68)Rα,M=2πe−iMα∫0∞r|gM(r)|2dr

is extracted digitally, then from the phase associated with this component it is possible to determine the angular shift that one version of the object has undergone. When an optical filter that is matched to Rα,M of a particular object is constructed, if that some object is entered as an input to the system with any angular rotation, a correlation peak of strength proportional to ∫0∞r|gM(r)|2dr will be produced, independent of rotation. Hence an optical correlator can be constructed that will recognize that object independent of rotation [31].

So far we have studied the circular harmonic correlation in the context of quantum optics, we have endowed the crosscorrelation Rα with a definite quantum mechanical meaning, i.e., the overlap between ⟨g| and the rotated state eiα(a1†a1−a2†a2)|g⟩, in the entangled state representation.

Note that Fourier-based correlators is also very sensitive to magnification. However, the magnitude of Mellin transform is independent of scale-size changes in the input [31]. Now we examine when |g⟩ is both rotated and squeezed (by a two-mode squeezing operator S2(λ)=exp[λ(a1†a2†−a1a2)]), then from Eqs. (35) and (62) we have

(69)S2(λ)|g⟩=∫d2ηπμ|η/μ⟩⟨η|g⟩=∫d2ηπμ|η/μ⟩g(r,θ)

it follows the overlap between ⟨g| and the state e−iα(a1†a1−a2†a2)S(λ)|g⟩,

(70)Wα,λ≡π⟨g|e−iα(a1†a1−a2†a2)S(λ)|g⟩=∫0∞rdrμ∫02πg*(r/μ,θ−α)g(r,θ)dθ

which corresponds to the cross-correlation arising from combination of squeezing and rotation (joint transform correlator). On the other hand, from Eqs. (45) and (35) we see

(71)S2(λ)|m,r⟩=12πμ∫02πdθ|η=rμeiθ⟩e−imθ=1μ|m,rμ⟩

and therefore

(72)Wα,λ=2πμ∑m=−∞∞e−imα∫0∞rgm(r)gm′*(re−λ)dr

from which one can see that to achieve simultaneous scale and rotation invariance, a two-dimensional object g(r,θ) should be entered into the optical system in a distorted polar coordinate system, the distortation arising from the fact that the radial coordinate is stretched by a logarithmic transformation (λ=−ln⁡μ), which coincides with Ref. [62]. The quantum optical version is thus established which is a new tie connecting Fourier optics and quantum optics [63].

At the end of this section, using the two-variable Hermite polynomials’ definition [64]

(73)Hm,n(ξ,ξ*)=∑l=0min⁡(m,n)m!n!(−1)lξm−lξ*n−ll!(m−l)!(n−l)!

which is quite different from the product of two single-variable Hermite polynomials, and its generating function formula is

(74)∑m,n=0∞tmt′⁣nm!n!Hm,n(ξ,ξ*)=exp⁡(−tt′+tξ+t′ξ*)

and noting Hm,n(ξ,ξ*)=ei(m−n)φHm,n(r′,r′), we can directly perform the integral in Eq. (51) and derive the explicit form of |s,r′⟩,

(75)|s,r′⟩=e−r′⁣2/2∑n=0∞1(n+s)!n!Hn+s,n(r′,r′)|n+s,n⟩

which is really an entangled state in two-mode Fock space. Eqs. (74) and (75) will be used often in the following discussions.

In the following we concentrate on finding the generalized Fresnel operators in both one- and two-mode cases with use of the IWOP technique.

5 Single-mode Fresnel operator as the image of the classical optical Fresnel transform

In this section we shall mainly introduce so-called generalized Fresnel operators (GFO) (in one-mode and two-mode cases both) [65] and some appropriate quantum optical representations (e.g. coherent state representation and entangled state representation) to manifestly link the formalisms in quantum optics to those in classical optics. In so doing, we find that the various transforms in classical optics are just the result of generalized Fresnel operators inducing transforms on appropriate quantum state vectors, i.e., classical optical Fresnel transforms have their counterpart in quantum optics. Besides, we can study the important ABCD rule obeyed by Gaussian beam propagation (also the ray propagation in matrix optics) [66] in the domain of quantum optics.

5.1 Single-mode GFO gained via coherent state method

For the coherent state |z⟩ in quantum optics [1–5]

(76)|z⟩=exp⁡(za†−z*a)|0⟩≡|(zz*)⟩

which is the eigenstate of annihilation operator a, a|z⟩=z|z⟩, using the IWOP and (12), we can put the over-completeness relation of |z⟩ into normal ordering ∫d2zπ|z⟩⟨z|=∫d2zπ:e−(z*−a†)(z−a):=1. The canonical form of coherent state |z⟩ is expressed as

(77)|z⟩=|p,q⟩=exp⁡[i(pQ−qP)]|0⟩≡|(qp)⟩

where z=(q+ip)/2. It follows that ⟨p,q|Q|p,q⟩=q,⟨p,q|P|p,q⟩=p, which indicates that the states |p,q⟩ generate a canonical phase-space representation for a state |Ψ⟩,Ψ(p,q)=⟨p,q|Ψ⟩. Thus the coherent state is a good candidate for providing with classical phase-space description of quantum systems. Considering that the Fresnel transform’s parameters (A,B,C,D) are elements of a ray transfer matrix M describing optical systems, M belongs to the unimodular symplectic group; and the coherent state |p,q⟩ is a good candidate for providing with classical phase-space description of quantum systems, we naturally think of that the symplectic transformation (ABCD)(qp) in classical phase space may be mapped onto a generalized Fresnel operator in Hilbert space through the coherent state basis. Thus we construct the following ket–bra projection operator

(78)∫∫−∞∞dxdp|(ABCD)(qp)⟩⟨(qp)|

as the GFO. In fact, using notation of |z⟩ (coherent state), and introducing complex numbers s and r,

(79)s=12[A+D−i(B−C)]r=−12[A−D+i(B+C)],|s|2−|r|2=1

from Eq. (77) we know

(80)|(ABCD)(qp)⟩=|(s−r−r*s*)(zz*)⟩≡|sz−rz*⟩=exp⁡[−12|sz−rz*|2+(sz−rz*)a†]|0⟩

(s−r−r*s*) is still a symplectic group element, so Eq. (78) becomes [67–70]

(81)F1(s,r)=s∫d2zπ|sz−rz*⟩⟨z|

where the factor s is attached to anticipate the unitarity of the operator F1. Eq. (78) tells us that c-number transform (q,p)→(Aq+Bp,Cq+Dp) in coherent state basis is mapped into F1(s,r). Now we prove F1(s,r) is really the FO we want. Using the IWOP technique and Eqs. (80) and (12) we can perform the integral

(82)F1(s,r)=exp⁡(−r2s*a†2)exp⁡[(a†a+12)ln⁡1s*]×exp⁡(r*2s*a2)

where we have used the mathematical formula [71, 72]

(83)∫d2zπexp⁡(ζ|z|2+ξz+ηz*+fz2+gz*2)=1ζ2−4fgexp⁡(−ζξη+ξ2g+η2fζ2−4fg)

with the convergent condition Re(ζ±f±g)<0, Re(ζ2−4fgζ±f±g)<0. It then follows

(84)⟨z|F1(s,r)|z′⟩=1s*exp⁡(−|z|2+|z′|22−rz*22s*+r*z′⁣22s*+z*z′s*)

Then using

(85)⟨xi|z⟩=π−1/4exp⁡(−xi22+2xiz−z22−|z|22)

and the completeness relation of coherent state as well as Eq. (79), we obtain the matrix element of F1(s,r))≡F1(A,B,C)) in coordinate representation ⟨xi|,

(86)⟨x2|F1(s,r)|x1⟩=12πiBexp⁡[i2B(Ax12−2x2x1+Dx22)]≡K(x2,x1)

which is just the kernel of generalized Fresnel transform (GFT) K(x2,x1) in Eq. (11). The above discussions demonstrate how to transit classical Fresnel transform to GFO through the coherent state and the IWOP technique.

Now if we define g(x2)=⟨x2|g⟩, f(x1)=⟨x1|f⟩ and using Eq. (16), we can rewrite Eq. (10) as

(87)⟨x2|g⟩=∫−∞∞dx1⟨x2|F1(A,B,C)|x1⟩⟨x1|f⟩=⟨x2|F1(A,B,C)|f⟩

which is just the quantum mechanical version of GFO. Therefore, the 1D GFT in classical optics corresponds to the 1-mode GFO F1(A,B,C) operating on state vector |f⟩ in Hilbert space, i.e., |g⟩=F1(A,B,C)|f⟩. One merit of GFO is: using coordinate-momentum representation transform we can immediately obtain GFT in “frequency” domain, i.e.,

(88)⟨p2|F|p1⟩=12πi(−C)exp⁡[i2(−C)(Dp12−2p2p1+Ap22)]

Obviously, F1(A,B,C) induces the following transform

(89)F1−1(A,B,C)(QP)F1(A,B,C)=(ABCD)(QP)

5.2 Group multiplication rule for single-mode GFO

Because two successive optical Fresnel transforms is still a Fresnel transform, we wonder if the product of two GFO is still a GFO. On the other hand, we have known that the GFO is the image of the symplectic transform z→sz−rz*, we expect that the product of two symplectic transforms are mapped into the GFO which is just the product of two GFOs. If so, then the correspondence between GFT and GFO is perfect. Using Eq. (81), ⟨z|z′⟩=exp⁡[−12(|z|2+|z′|2)+z*z′] and the IWOP technique we can directly perform the following integrals

(90)F1(s,r)F1(s′,r′)=ss′∫d2zd2z′π2|sz−rz*⟩⟨z|s′z′−r′z′⁣*⟩⟨z′|=s′′∫d2zπ|s′′⁣*z−r′′z*⟩⟨z|=F1(s′′,r′′)

where we have set

(91)s′′=ss′+rr′⁣*,r′′=r′s+rs′⁣*

or

(92)M′′≡(s−r−r*s*)(s′−r′−r′⁣*s′⁣*)=MM′|s′′|2−|r′′|2=1

From which we see that it is just the mapping of the above (A,B,C,D) matrices multiplication. Hence F1(s,r)F1(s′,r′) is the loyal representation of the product of two symplectic group elements shown in Eq. (92). The above discussion actually reveals an important property of coherent states, though two coherent state vectors are not orthogonal, but the equation

(93)ss′∫d2zd2z′π2|sz−rz*⟩⟨z|s′z′−r′z′⁣*⟩⟨z′|=s′′∫d2zπ|s′′⁣*z−r′′z*⟩⟨z|

seems as if their overlap ⟨z|s′z′−r′z′⁣*⟩ was a δ-function. The coherent state representation for GFOs’ product may be visualized very easily, but it achieves striking importance, because it does not change its form when treating symplectic transform according to z→sz−rz*.

As a result of this group multiplication rule of GFO, we immediately obtain

(94)KM′′(x2,x1)=⟨x2|F1(s′′,r′′)|x1⟩=∫−∞∞dx3KM(x2,x3)KM′(x3,x1)

provided that the parameter matrices (s′′,r′′) satisfy Eq. (91). Thus by virtue of the group multiplication property of GFO we immediately find the successive transform property of GFTs.

6 Quantum optical ABCD law for optical propagation — single-mode case

In classical optics, ray-transfer matrices, M=(ABCD),AD−BC=1 have been used to describe the geometrical formation of images by a centered lens system. For an optical ray (a centered spherical wavefront) passing through optical instruments there is a famous law, named ABCD law, governing the relation between input ray (r1,α1) and output ray (r2,α2), i.e.,

(95)(r2α2)=N(r1α1)

where r1 is the ray height from the optical axis, and α1 is named the optical direction-cosine; r1/α1≡R1 specifies the ray’s wavefront shape. Eq. (95) implies

(96)R2≡r2α2=AR1+BCR1+D

This law is the core of matrix optics, since it tells us how the curvature of a centered spherical wavefront changes from one reference plane to the next. Besides, the multiplication rule of matrix optics implies that if the ray-transfer matrices of the n optical components are M1,M2,M3,⋯,Mn, respectively, then the whole system is determined by a matrix M=M1M2M3⋯Mn.

One of the remarkable things of modern optics is the case in which geometrical ray-transfer methods, constituting the matrix optics, can be adapted to describe the generation and propagation of Laser beams. In 1965, Kogelnik [73] pointed out that propagation of Gaussian beam also obeys ABCD law via optical diffraction integration, i.e., the input light field f(x1) and output light field g(x2) are related to each other by so-called Fresnel integration [31] g(x2)=∫−∞∞K(A,B,C;x2,x1)f(x1)dx1, where

(96)K(A,B,C;x2,x1)=12πiBexp⁡[i2B(Ax12−2x2x1+Dx22)]

The ABCD law for Gaussian beam passing through an optical system is [74, 75]

(97)q2=Aq1+BCq1+D

where q1(q2) represents the complex curvature of the input (output) Gaussian beam, Eq. (97) has the similar form as Eq. (96). An interesting and important question naturally arises [76]: Does ABCD law also exhibit in quantum optics? Since classical Fresnel transform should have its quantum optical counterpart?

To see the ABCD law more explicitly, using Eq. (79) we can re-express Eq. (82) as

(98)F1(A,B,C)=2A+D+i(B−C):exp{A−D+i(B+C)2[A+D+i(B−C)]a†2+[2A+D+i(B−C)−1]a†a−A−D−i(B+C)2[A+D+i(B−C)]a2}:

and the multiplication rule for F1 is F(A′,B′,C′,D′)F(A,B,C,D)=F(A′′,B′′,C′′,D′′), where

(99)(A′′B′′C′′D′′)=(A′B′C′D′)(ABCD)

Next, we directly use the GFO to derive ABCD law in quantum optics. From Eq. (99) we see that the GFO generates

(100)F1(A,B,C)|0⟩=2A+iB−i(C+iD)×exp⁡{A−D+i(B+C)2[A+D+i(B−C)]a†2}|0⟩

If we identify

(101)A−D+i(B+C)A+D+i(B−C)=q1−iq1+i

then

(102)F1(A,B,C)|0⟩=−2/(C+iD)q1+iexp⁡[q1−i2(q1+i)a†2]|0⟩

The solution of Eq. (101) is

(103)q1≡−A+iBC+iD

Let F1(A,B,C)|0⟩ expressed by Eq. (102) be an input state for an optical system which is characteristic by parameters A′,B′,C′,D′, then the quantum optical ABCD law states that the output state is

(104)F1(A′,B′,C′)F1(A,B,C)|0⟩=−2/(C′′+iD′′)q2+iexp⁡[q2−i2(q2+i)a†2]|0⟩

which has the similar form as Eq. (102), where (C′′,D′′) is determined by Eq. (99), and

(105)q¯2=A′q¯1+B′C′q¯1+D′,q¯j≡−qj,j=1,2

which resembles Eq. (97).

Proof: According to the multiplication rule of two GFOs and Eqs. (98) and (99), we have

(106)F1(A′,B′,C′)F1(A,B,C)|0⟩=−2/(C+iD)A′q1−B′−i(C′q1−D′)×exp⁡{A′q1−B′+i(C′q1−D′)2[A′q1−B′−i(C′q1−D′)]a†2}|0⟩

Using Eq. (103) we see 2/(C+iD)C′q1−D′=−2/(C′′+iD′′), along with Eq. (105) we can reach Eq. (104), thus the law is proved. Using Eq. (102) we can re-express Eq. (105) as

(107)q2=−A′(A+iB)+B′(C+iD)C′(A+iB)+D′(C+iD)=−A′′+iB′′C′′+iD′′

which is in consistent to Eq. (103). Eqs. (102)–(107) are therefore self-consistent.

As an application of quantum optical ABCD law, we apply it to tackle the time-evolution of a time-dependent harmonic oscillator whose Hamiltonian is

(108)H=12e−2γtP2+12ω02e2γtQ2,ℏ=1

where we have set the initial mass m0=1, γ denotes damping. Using u(t)=eiγ2Q2e−iγt2(QP+PQ) to perform the transformation

(109)u(t)Qu−1(t)=e−γtQu(t)Pu−1(t)=eγtP−γeγtQ

then i∂|ψ(t)⟩∂t=H|ψ(t)⟩ leads to i∂|ϕ⟩∂t=H|ϕ⟩,|ϕ⟩=u(t)|ψ(t)⟩,

(110)H→H=u(t)Hu−1(t)−iu(t)∂u−1(t)∂t=12P2+12ω2Q2

where ω2=ω02−γ2. H does not contain t explicitly. The dynamic evolution of a mass-varying harmonic oscillator from the Fock state |0⟩ at initial time to a squeezed state at time t is

(111)|ψ(t)⟩0=u−1(t)|0⟩=eiγt2(QP+PQ)e−iγ2Q2|0⟩

if we let A=D=1,B=0,C=−γ; and A′=e−γt,D′=eγt,B′=C′=0, then q1=1γ−i,q2=e−2γtγ−i, according to Eq. (104) we directly obtain

(112)u−1(t)|0⟩=2e−γte−2γt+iγ+1×exp⁡[e−2γt−1−iγ2(e−2γt+1+iγ)a†2]|0⟩

so the time evolution of the damping oscillator embodies the quantum optical ABCD law.

7 Optical operator method studied via GFO’s decomposition

Fresnel diffraction is the core of Fourier optics [31, 43, 48–50]. Fresnel transform is frequently used in optical imaging, optical propagation and optical instrument design. The GFT represents a class of optical transforms which are of great importance for their applications to describe various optical systems. It is easily seen that when we let the transform kernel K(x2,x1)=exp⁡(ix2x1), the GFT changes into the well-known Fourier transform, which is adapted to express mathematically the Fraunhofer diffraction. And if K(x2,x1)=exp⁡[i(x2−x1)2], the GFT then describes a Fresnel diffraction. In studying various optical transformations someone also proposed so-called optical operator method [77–79] which used quantum mechanical operators’ ordered product to express the mechanism of optical systems, such that the ray transfer through optical instruments and the diffraction can be discussed by virtue of the commutative relations of operators and the matrix algebra. Two important questions thus naturally arises: how to directly map the classical optical transformations to the optical operator method? How to combine the usual optical transformation operators, such as the square phase operators, scaling operator, Fourier transform operator and the propagation operator in free space, into a concise and unified form? In this section we shall solve these two problems and develop the optical operator method onto a new stage.

7.1 Four fundamental optical operators derived by decomposing GFO

The GFO F1(A,B,C) can also be expressed in the form of quadratic combination of canonical operators Q and P [80], i.e.,

(113)F1(A,B,C)=exp⁡(iC2AQ2)exp⁡[−i2(QP+PQ)ln⁡A]×exp⁡(−iB2AP2)

where we have set ℏ=1,A≠0. To confirm this, we first calculate matrix element

(114)⟨x|F1(A,B,C)|p⟩=12πAexp⁡(iC2Ax2−iB2Ap2+ipxA)

where we have used the squeezing property

(115)exp⁡[−i2(QP+PQ)ln⁡A]|p⟩=1A|p/A⟩

It then follows from Eq. (114) and AD−BC=1, we have

(116)⟨x2|F1(A,B,C)|x1⟩=∫−∞∞dp⟨x2|F1(A,B,C)|p⟩⟨p|x1⟩=KM(x2,x1)

Thus F1(A,B,C) in Eq. (113) is really the expected GFO. Next we directly use Eq. (78) and the canonical operator (Q,P) representation (113) to develop the optical operator method.

Noticing the matrix decompositions [44]

(117)(ABCD)=(10C/A1)(A00A−1)(1B/A01)

and comparing Eqs. (78) and (113) as well as using (99) we know

(118)F1(A,B,C)=F1(1,0,C/A)F1(A,0,0)F1(1,B/A,0)

where

(119)F1(1,0,C/A)=exp⁡(iC2AQ2)

(120)F1(1,B/A,0)=exp⁡(−iB2AP2)

(121)F1(A,0,0)=exp⁡[−i2(QP+PQ)ln⁡A]

which are named quadrature phase operator, Fresnel propagator in free space, and scaling operator (squeezed operator [14–17]), respectively. When A=D=0,B=1,C=−1, from Eq. (98) we see

(122)F1(0,1,−1)=exp⁡[−iπ2(a†a+12)]

which is named the Fourier operator, since it quantum mechanically transforms [71, 72]

(123)F1†(0,1,−1)QF1(0,1,−1)=PF1†(0,1,−1)PF1(0,1,−1)=−Q

7.2 Alternate decompositions of GFO

Note that when A=0, the decomposition (113) is not available, instead, from

(124)(ABCD)−1=(D−B−CA)

and Eqs. (98), (99) and (113) we have

(125)F1−1(A,B,C)=exp⁡(−iC2DQ2)exp⁡[−i2(QP+PQ)ln⁡D]exp⁡(iB2DP2)

It then follows that

(126)F1(A,B,C)=exp⁡(−iB2DP2)exp⁡[i2(QP+PQ)ln⁡D]exp⁡(iC2DQ2),D≠0

Besides, when we notice

(127)(ABCD)=(10D/B1)(B001/B)(01−10)(10A/B1)

and

(128)(ABCD)=(1A/C01)(−1/C00−C)(01−10)(1D/C01)

we have another decomposition for B≠0,

(129)F1⁢(A,B,C)=exp⁡(i⁢D2⁢B⁢Q2)⁢exp⁡[−i2⁢(Q⁢P+P⁢Q)⁢ln⁡B]×exp⁡[−i⁢π2⁢(a†⁢a+12)]⁢exp⁡(i⁢A2⁢B⁢Q2)

and for C≠0,

(130)F1⁢(A,B,C)=exp⁡(−i⁢A2⁢C⁢P2)×exp⁡[−i2⁢(Q⁢P+P⁢Q)⁢ln⁡(−1C)]×exp⁡[−i⁢π2⁢(a†⁢a+12)]⁢exp⁡(−i⁢D2⁢C⁢P2)

7.3 Some optical operator identities

For a special optical systems with the parameter A=0,C=−B−1,

(131)(0B−B−1D)=(10D/B1)⁢(B00B−1)⁢(01−10)

we have

(132)exp⁡(−i⁢B2⁢D⁢P2)⁢exp⁡[i2⁢(Q⁢P+P⁢Q)⁢ln⁡D]×exp⁡(−i2⁢D⁢B⁢Q2)=exp⁡(−i⁢D2⁢B⁢Q2)⁢exp⁡[i2⁢(Q⁢P+P⁢Q)⁢ln⁡B]×exp⁡[−i⁢π2⁢(a†⁢a+12)]

In particular, when A=D=0,C=−B−1, from

(133)(0B−1B0)=(B001B)⁢(01−10)

we have

(134)exp⁡[−B2−12⁢(B2+1)⁢a†2]⁢exp⁡[(a†⁢a+12)⁢ln⁡(−2⁢B⁢iB2+1)]×exp⁡[−B2−12⁢(B2+1)⁢a2]=exp⁡[−i2⁢(Q⁢P+P⁢Q)⁢ln⁡B]⁢exp⁡[−i⁢π2⁢(a†⁢a+12)]

Using the following relations

(135)(ABCD)=(1(A−1)/C01)⁢(10C1)×(1(D−1)/C01)

it then follows that

(136)F1⁢(A,B,C)=exp⁡[−i⁢(A−1)2⁢C⁢P2]⁢exp⁡(i⁢C2⁢Q2)×exp⁡[−i⁢(D−1)2⁢C⁢P2]

While from

(137)(ABCD)=(10(D−1)/B1)⁢(1B01)×(10(A−1)/B1)

we obtain

(138)F1⁢(A,B,C)=exp⁡[i⁢(D−1)2⁢B⁢Q2]⁢exp⁡(−i⁢B2⁢P2)×exp⁡[i⁢(A−1)2⁢B⁢Q2]

so we have

(139)exp⁡[i⁢(D−1)2⁢B⁢Q2]⁢exp⁡(−i⁢B2⁢P2)⁢exp⁡[i⁢(A−1)2⁢B⁢Q2]=exp⁡[−i⁢(A−1)2⁢C⁢P2]⁢exp⁡(i⁢C2⁢Q2)⁢exp⁡[−i⁢(D−1)2⁢C⁢P2]

In this section, based on a one-to-one correspondence between classical Fresnel transform in phase space and quantum unitary transform in state-vector space and the IWOP technique as well as the coherent state representation, we have found a way to directly map the classical optical transformations to the optical operator method. We have combined the usual optical transformation operators, such as the square phase operators, scaling operator, Fourier transform operator and the propagation operator in free space, into a concise and unified form. The various decompositions of Fresnel operator into the exponential canonical operators are also obtained.

8 Quantum tomography and probability distribution for the Fresnel quadrature phase

In quantum optics theory all possible linear combinations of quadratures Q and P of the oscillator field mode a and a† can be measured with the homodyne measurement just by varying the phase of the local oscillator. The average of the random outcomes of the measurement, at a given local oscillator phase, is connected with the marginal distribution of Wigner function (WF). Thus the homodyne measurement of light field permits the reconstruction of the WF of a quantum system by varying the phase shift between two oscillators. In Ref. [81] Vogel and Risken pointed out that the probability distribution for the rotated quadrature phase Qθ≡[a†⁢exp⁡(i⁢θ)+a⁢exp⁡(−i⁢θ)]/2,[a,a†]=1, which depends on only one θ angle, can be expressed in terms of WF, and that the reverse is also true (named as Vogel–Risken relation). That is, one can obtain the Wigner distribution by tomographic inversion of a set of measured probability distributions, Pθ⁢(qθ), of the quadrature amplitude. Once the distribution Pθ⁢(qθ) are obtained, one can use the inverse Radon transformation familiar in tomographic imaging to obtain the Wigner distribution and density matrix. The Radon transform of the WF is closely related to the expectation values or densities formed with the eigenstates to the rotated canonical observables. The field of problems of the reconstruction of the density operator from such data is called quantum tomography. (Optical tomographic imaging techniques derive 2D data from a three-dimensional object to obtain a slice image of the internal structure and thus have the ability to peer inside the object noninvasively. The slice image is equated with tomogram.) The theoretical development in quantum tomography in the last decade has progressed in the direction of determining more physical relevant parameters of the density from tomographic data [81–86].

8.1 Relation between Fresnel transform and Radon transform of WF

In Refs. [87, 88], the Radon transform of WF which depends on two continuous parameters is introduced. This has the advantage in conveniently associating quantum tomography theory with squeezed coherent state theory. In this subsection we want to derive relations between the Fresnel transform and the Radon transform of WF in quantum optics in tomography theory.

By extending the rotated quadrature phase Qθ to the Fresnel quadrature phase

(140)QF≡(s*⁢a+r⁢a†+s⁢a†+r*⁢a)/2=F1⁢Q⁢F1†

where s and r are related to ABCD through Eq. (79), we shall prove that the (D,B) related Radon transform of Wigner operator Δ⁢(q,p) is just the pure state density operator |q⟩s,r⁢s,r⁢⟨q| (named as the tomographic density operator) formed with the eigenstates belonging to the quadrature QF (|q⟩s,r=F1⁢|q⟩; Q is the coordinate operator),

(141)F1⁢|q⟩⁢⟨q|⁢F1†=|q⟩s,r⁢s,r⁢⟨q|=∫−∞∞dq′⁢dp′⁢δ⁢[q−(D⁢q′−B⁢p′)]⁢Δ⁢(q′,p′)

Since F corresponds to classical Fresnel transform in optical diffraction theory, so Eq. (141) indicates that the probability distribution for the Fresnel quadrature phase is the Radon transform of WF [89].

Proof: First, from Eq. (82) we see

(142)F1⁢(s,r)⁢a⁢F1†⁢(s,r)=s*⁢a+r⁢a†

so from Q=a+a†2,P=i⁢a†−a2, we have indeed

(143)F1⁢Q⁢F1†=F1⁢a+a†2⁢F1†=(s*⁢a+r⁢a†+s⁢a†+r*⁢a)/2=QF

Second, we can derive the explicit form of |q⟩s,r. Starting from s*+r*=D+i⁢B, s*−r*=A−i⁢C, we set up the eigenvector equation

(144)QF⁢|q⟩s,r=(D⁢Q−B⁢P)⁢|q⟩s,r=q⁢|q⟩s,r

It follows that

(145)|q⟩s,r=F1⁢(s,r)⁢|q⟩

In the coordinate and momentum representations we have

(146)⟨q′|QF|q⟩s,r=(D⁢q′+i⁢B⁢dd⁢q′)⁢⟨q′|q⟩s,r=q⁢⟨q′|q⟩s,r

(147)⟨p|QF|q⟩s,r=(i⁢D⁢dd⁢p−B⁢p)⁢⟨p|q⟩s,r=q⁢⟨p|q⟩s,r

The normalizable solutions to Eqs. (146) and (147) are

(148)⟨q′|q⟩s,r=c⁢(q)⁢exp⁡[i⁢q′⁢(D⁢q′−2⁢q)2⁢B]

(149)⟨p|q⟩s,r=d⁢(q)⁢exp⁡[i⁢p⁢(−B⁢p−2⁢q)2⁢D]

Using the Fock representation of |q⟩ and |p⟩ in Eqs. (13) and (22), we obtain

(150)|q⟩s,r=π−1/4c(q)2⁢B⁢πB−i⁢Dexp[−q22⁢B⁢(B−i⁢D)+2⁢a†⁢qD+i⁢B−D−i⁢BD+i⁢Ba†22]|0⟩

and

(151)|q⟩s,r=d(q)π−1/42⁢π⁢DD+i⁢Bexp[−q22⁢D⁢(D+i⁢B)+2⁢a†⁢qD+i⁢B−D−i⁢BD+i⁢Ba†22]|0⟩

Comparing Eq. (150) with Eq. (151) we see

(152)c⁢(q)d⁢(q)=Di⁢B⁢exp⁡(i⁢A2⁢B⁢q2−i⁢C⁢q22⁢D)

On the other hand, according to the orthogonalization of |q⟩s,r, ⟨q′|q′′⟩s,rs,r=δ(q′−q′′), we have

(153)|c⁢(q)|2=12⁢π⁢B,|d⁢(q)|2=12⁢π⁢D

Thus, combining Eqs. (152) and (153) we deduce

(154)c⁢(q)=12⁢π⁢i⁢B⁢exp⁡(i⁢A2⁢B⁢q2)d⁢(q)=12⁢π⁢D⁢exp⁡(i⁢C⁢q22⁢D)

and

(155)|q⟩s,r=π−1/4D+i⁢Bexp(−A−i⁢CD+i⁢Bq22+2⁢qD+i⁢Ba†−D−i⁢BD+i⁢Ba†22)|0⟩

or

(156)|q⟩s,r≡π−1/4s*+r*exp(−s*−r*s*+r*q22+2⁢qs*+r*a†−s+rs*+r*a†22)|0⟩

It is easily seen that that |q⟩s,r make up a complete set ∫−∞∞dq⁢|q⟩ss,r,r⁢⟨q|=1 (so |q⟩s,r can be named as the tomography representation). Then according to the Weyl quantization scheme [51]

(157)H⁢(Q,P)=∫−∞∞dp⁢dq⁢Δ⁢(q,p)⁢h⁢(q,p)

where h⁢(q,p) is the Weyl correspondence of H⁢(Q,P),

(158)h⁢(q,p)=2⁢π⁢Tr⁢[H⁢(Q,P)⁢Δ⁢(q,p)]

Δ⁢(q,p) is the Wigner operator [42, 90–94],

(159)Δ⁢(q,p)=12⁢π⁢∫−∞∞du⁢ei⁢p⁢u⁢|q+u2⟩⁢⟨q−u2|

and using Eqs. (158), (159) and (150) we know that the classical Weyl correspondence (Weyl image) of the projection operator |q⟩s,r⁢s,r⁢⟨q| is

(160)2⁢π⁢Tr⁢[Δ⁢(q′,p′)⁢|q⟩s,r⁢s,r⁢⟨q|]=δ⁢[q−(D⁢q′−B⁢p′)]

which means

(161)|q⟩s,r⁢s,r⁢⟨q|=∫−∞∞dq′⁢dp′⁢δ⁢[q−(D⁢q′−B⁢p′)]⁢Δ⁢(q′,p′)

Combining Eqs. (142)–(145) together we complete the proof. Therefore, the probability distribution for the Fresnel quadrature phase is the Radon transform of WF

(162)|⟨q|F1†|ψ⟩|2=|s,r⟨q|ψ⟩|2=∫−∞∞dq′⁢dp′⁢δ⁢[q−(D⁢q′−B⁢p′)]⁢⟨ψ|Δ⁢(q′,p′)|ψ⟩

So we name |q⟩s,r⁢s,r⁢⟨q|the tomographic density. Moreover, the tomogram of quantum state |ψ⟩ is just the squared modulus of the wave function ⟨q|ψ⟩s,r. This new relation between quantum tomography and optical Fresnel transform may provide experimentalists to figure out new approach for generating tomography.

The introduction of |q⟩s,r also brings convenience to obtaining the inverse of Radon transformation. Using the complete relation of |q⟩s,r and Eq. (144) we have

(163)e−i⁢g⁢QF=∫−∞∞dq⁢|q⟩ss,r,r⁢⟨q|⁢e−i⁢g⁢q=∫−∞∞dq⁢dp⁢Δ⁢(q,p)⁢e−i⁢g⁢(D⁢q−B⁢p)

Considering its right hand-side as a Fourier transformation, its reciprocal transform is

(164)Δ(q,p)=14⁢π2∫−∞∞dq′∫−∞∞dg′|g′|∫0πdφ|q′⟩s,rs,r⟨q′|×exp⁡[−i⁢g′⁢(q′D2+B2−q⁢cos⁡φ−p⁢sin⁡φ)]

where g′=g⁢D2+B2,cos⁡φ=DD2+B2,sin⁡φ=−BD2+B2. So once the distribution |s,r⟨q|ψ⟩|2 are obtained, one can use the inverse Radon transformation familiar in tomographic imaging to obtain the Wigner distribution.

By analogy, we can conclude that the (A,C) related Radon transform of Δ⁢(q,p) is just the pure state density operator |p⟩s,r⁢s,r⁢⟨p| formed with the eigenstates belonging to the conjugate quadrature of QF,

(165)F1⁢|p⟩⁢⟨p|⁢F1†=|p⟩s,r⁢s,r⁢⟨p|=∫−∞∞dq′⁢dp′⁢δ⁢[p−(A⁢q′−C⁢p′)]⁢Δ⁢(q′,p′)

Similarly, we find that for the momentum density,

(166)F1⁢|p⟩⁢⟨p|⁢F1†=|p⟩s,r⁢s,r⁢⟨p|=∫−∞∞dq′⁢dp′⁢δ⁢[p−(A⁢p−C⁢q′)]⁢Δ⁢(q′,p′)

where

(167)F1|p⟩=|p⟩s,r=π−1/4A−i⁢Cexp(−D+i⁢BA−i⁢Cp22+2⁢i⁢pA−i⁢Ca†+A+i⁢CA−i⁢Ca†22)|0⟩

As an application of the relation (141), recalling that the F1⁢(r,s) makes up a faithful representation of the symplectic group, it then follows from Eq. (141) that

(168)F1′⁢(r′,s′)⁢F1⁢(r,s)⁢|q⟩⁢⟨q|⁢F1†⁢(r,s)⁢F1′⁣†⁢(r′,s′)=|q⟩s′′,r′′⁢s′′,r′′⁢⟨q|=∫∫−∞∞dq′dp′δ[q−((B′C+DD′)q′−(AB′+BD′)p′)]Δ(q′,p′)

In this way a complicated Radon transform of tomography can be viewed as the sequential operation of two Fresnel transforms. This confirms that the continuous Radon transformation corresponds to the symplectic group transformation [87, 88]. This is an advantage of introducing the Fresnel operator. The group property of Fresnel operators help us to analyze complicated Radon transforms in terms of some sequential Fresnel transformations. The new relation may provide experimentalists to figure out new approach for realizing tomography.

8.2 A theorem for the calculation of the tomogram

In this subsection, we introduce a new theorem, which states that the tomogram of a density operator ρ is equal to the marginal integration of the classical Weyl correspondence function of F†⁢ρ⁢F1, where F1 is the Fresnel operator.

Multiply the both sides of Eq. (161) by a density matrix ρ. Then perform the trace. Note the Wigner function W⁢(p,q)=Tr⁢[ρ⁢Δ⁢(p,q)]. One can see

(169)Tr⁢[∫∫−∞∞⁢d⁢q′⁢d⁢p′⁢δ⁢[q−(D⁢q′−B⁢p′)]⁢Δ⁢(q′,p′)⁢ρ]=Tr⁢(|q⟩s,r⁢s,r⁢⟨q|⁢ρ)=s,r⟨q|ρ|q⟩s,r=⟨q|F†⁢ρ⁢F1|q⟩=∫∫−∞∞⁢d⁢q′⁢d⁢p′⁢δ⁢[q−(D⁢q′−B⁢p)]⁢W⁢(p,q)

The right hand side of Eq. (169) is commonly defined as the tomogram of quantum states in (B,D) direction, so in our view the calculation of tomogram in (B,D) direction is ascribed to the calculation of

(170)⟨q|F†⁢ρ⁢F1|q⟩≡Ξ

This is a concise and neat formula. Similarly, the tomogram in (A,C) direction is ascribed to ⟨p|F†⁢ρ⁢F1|p⟩.

According to the Weyl correspondence rule

(171)H⁢(X,P)=∫∫−∞∞⁢d⁢p⁢d⁢x⁢h⁢(p,x)⁢Δ⁢(p,x)

and the Weyl ordering form of Δ⁢(p,q)

(172)Δ⁢(p,q)=::⁢δ⁢(q−Q)⁢δ⁢(p−P)⁢::

where the symbol ::⁢:: denotes Weyl ordering, the classical correspondence of a Weyl ordered operator ::⁢h⁢(Q,P)⁢:: is obtained just by replacing Q→q,P→p in h, i.e.,

(173)::⁢h⁢(Q,P)⁢::=∫∫−∞∞⁢d⁢p⁢d⁢q⁢h⁢(p,q)⁢Δ⁢(p,q)

Let the classical Weyl correspondence of F†⁢ρ⁢F1 be h⁢(p,q)

(173)F†⁢ρ⁢F=∫∫−∞∞⁢d⁢p⁢d⁢q⁢h⁢(p,q)⁢Δ⁢(p,q)

then using Eqs. (170) and (159) we have

(174)Ξ=⟨q|∫∫⁢d⁢p⁢d⁢q′⁢h⁢(p,q′)⁢Δ⁢(p,q′)|q⟩=∫−∞∞d⁢p2⁢π⁢h⁢(p,q)

Thus we reach a theorem as below.

The tomogram of a density operator ρ is equal to the marginal integration of the classical Weyl correspondence h⁢(p,q) of F†⁢ρ⁢F1, where F is the Fresnel operator, expressed by

(175)Tr⁢[ρ⁢|q⟩s,r⁢s,r⁢⟨q|]=∫−∞∞d⁢p2⁢π⁢h⁢(p,q)

or

(176)Tr⁢[ρ⁢|p⟩s,r⁢s,r⁢⟨p|]=∫−∞∞d⁢x2⁢π⁢h⁢(p,q)

In this way the relationship between tomogram of a density operator ρ and the Fresnel transformed ρ′s classical Weyl function is established.

9 Two-mode GFO and its application

For two-dimensional optical Fresnel transforms [see Eq. (4)] in the x–y plane one may naturally think that the 2-mode GFO is just the direct product of two independent 1-mode GFOs, but with the same (A,B,C,D) matrix. However, here we present another 2-mode Fresnel operator which cannot only lead to the usual 2-dimensional optical Fresnel transforms in some appropriate quantum mechanical representations, but also provide us with some new classical transformations (we name them entangled Fresnel transformations.).

9.1 Two-mode GFO gained via coherent state representation

Similar to the single-mode case, we introduce the two-mode GFO F2⁢(r,s) through the following 2-mode coherent state representation [67–70]

(177)F2⁢(r,s)=s⁢∫d2⁢z1⁢d2⁢z2π2⁢|s⁢z1+r⁢z2*,r⁢z1*+s⁢z2⟩⁢⟨z1,z2|

which indicates that F2⁢(r,s) is a mapping of classical sympletic transform (z1,z2)→(s⁢z1+r⁢z2*,r⁢z1*+s⁢z2) in phase space. Concretely, the ket in Eq. (177) is

(178)|s⁢z1+r⁢z2*,r⁢z1*+s⁢z2⟩≡|s⁢z1+r⁢z2*⟩1⊗|r⁢z1*+s⁢z2⟩2

s and r are complex and satisfy the unimodularity condition. Using the IWOP technique we perform the integral in Eq. (177) and obtain

(179)F2(r,s)=exp(rs*a1†a2†)exp[(a1†a1+a2†a2+1)ln(s*)−1]×exp⁡(−r*s*⁢a1⁢a2)

Thus F2⁢(r,s) induces the transform

(180)F2⁢(r,s)⁢a1⁢F2−1⁢(r,s)=s*⁢a1−r⁢a2†F2⁢(r,s)⁢a2⁢F2−1⁢(r,s)=s*⁢a2−r⁢a1†

and F2 is actually a general 2-mode squeezing operator. Recall that Eq. (35) implies the intrinsic relation between the EPR entangled state and the two-mode squeezed state, which has physical implementation. That is, in the output of a parametric down-conversion the idler-mode and the signal-mode constitute a two-mode squeezed state. Meanwhile, they are entangled with each other in frequency domain. We naturally select the entangled state representation to relate F2⁢(r,s) to two-dimensional GFT. Letting |g⟩=F2⁢(r,s)⁢|f⟩, and then project |ψ⟩ onto the entangled state ⟨η′| defined by Eq. (28), and using the completeness relation (30) of |η⟩, we obtain

(181)g⁢(η′)≡⟨η′|g⟩=⟨η′|F2⁢(r,s)|f⟩=∫d2⁢ηπ⁢⟨η′|F2⁢(r,s)|η⟩⁢⟨η|f⟩≡∫d2⁢η⁢K2(r,s)⁢(η′,η)⁢f⁢(η)

Then using the overcompleteness relation of the coherent state, we can calculate the integral kernel

(182)K2(r,s)⁢(η′,η)=12⁢i⁢B⁢π×exp⁡[i2⁢B⁢(A⁢|η|2−(η⁢η′⁣*+η*⁢η′)+D⁢|η′|2)]≡K2M⁢(η′,η)

where the relation between s,r and (A,B,C,D) in Eq. (79) is used and the superscript M only means that the parameters of K2M are [A,B;C,D], and the subscript 2 means the 2D kernel. Eq. (182) has the similar form as Eq. (11) except for its complex form. Taking η1=x1,η2=x2 and η1′=x1′, η2′=x2′, we have

(183)K2M⁢(η′,η)=K2M⁢(x1′,x2′;x1,x2)=K1M⁢(x1,x1′)⊗K1M⁢(x2,x2′)

This shows that F2⁢(r,s) is really the counterpart of the 2D GFT.

If taking the matrix element of F2⁢(r,s) in the |ξ⟩ representation which is conjugate to |η⟩, we obtain the 2-dimensional GFT in its “frequency domain”, i.e.,

(184)⟨ξ′|F2(r,s)|ξ⟩=12⁢i⁢(−C)⁢πexp[i2⁢(−C)(D|ξ|2+A|ξ′|2−ξ′⁣*ξ−ξ′ξ*)]≡K2N(ξ′,ξ)

where the superscript N means that this transform kernel corresponds to the parameter matrix M=[D,−C,−B,A].

The two-mode GFO also abides by the group multiplication rule. Using the IWOP technique and Eq. (177) we can also obtain

(185)F2⁢(r,s)⁢F2⁢(r′,s′)=F2⁢(r′′,s′′)

where (r′′,s′′) are given by Eq. (91) or Eq. (92). Therefore, Eq. (185) is a loyal representation of the multiplication rule for ray transfer matrices in the sense of Matrix Optics.

9.2 Quantum optical ABCD law for two-mode GFO

Next we extended quantum optical ABCD law to the two-mode case. Operating with F2⁢(r,s) on two-mode number state |m,n⟩ and using the overlap between coherent state and number state, i.e., ⟨z1,z2|m,n⟩=z1*m⁢z2*n/m!⁢n!⁢exp⁡[−12⁢(|z1|2+|z2|2)], and the integral formula [71, 72]

(186)Hm,n⁢(ξ,η)=(−1)n⁢eξ⁢η⁢∫d2⁢zπ⁢zn⁢z*m⁢e−|z|2+ξ⁢z−η⁢z*

we can calculate

(187)F2⁢(r,s)⁢|m,n⟩=r*ns*n+1⁢m!⁢n!⁢Hm,n⁢[a1†s*,a2†r*]×exp⁡(r⁢a1†⁢a2†s*)⁢|00⟩

where Hm,n⁢(ϵ,ε)is the two variables Hermite polynomial [87, 88], shown in Eqs. (73) and (74).

Using Eqs. (79) and (103), we recast Eq. (187) into

(188)F2⁢(r,s)⁢|m,n⟩=−2/(C+i⁢D)(q1+i)⁢m!⁢n!⁢(−q1*+iq1+i⁢C−i⁢DC+i⁢D)n×Hm,n⁢[−2⁢a1†/(C+i⁢D)q1+i,2⁢a2†/(C−i⁢D)q1*+i]×exp⁡(−q1−iq1+i⁢a1†⁢a2†)⁢|00⟩

Taking notice of the multiplication rule of F2⁢(r,s) in Eq. (185), which is equivalent to

(189)F2⁢(A′,B′,C′)⁢F2⁢(A,B,C)=F2⁢(A′′,B′′,C′′)

where (A′,B′,C′),(A,B,C) and (A′′,B′′,C′′) are related to each other by Eq. (99).

Next we directly use the GFO to derive ABCD rule in quantum optics for Gaussian beam in the two-mode case. According to Eq. (187) and Eq. (189) we obtain

(190)F2⁢(A′,B′,C′)⁢F2⁢(A,B,C)⁢|m,n⟩=r′′⁣*ns′′⁣*n+1⁢m!⁢n!⁢Hm,n⁢[a1†s′′⁣*,a2†r′′⁣*]⁢exp⁡(r′′⁢a1†⁢a2†s′′⁣*)⁢|00⟩

Similar to the way of deriving Eq. (188), we can simplify Eq. (190) as

(191)F2⁢(A′,B′,C′)⁢F2⁢(A,B,C)⁢|00⟩=−2/(C′′+i⁢D′′)(q2+i)⁢m!⁢n!⁢(−q2*+iq2+i⁢C′′−i⁢D′′C′′+i⁢D′′)n×Hm,n⁢[−2⁢a1†/(C′′+i⁢D′′)q2+i,2⁢a2†/(C′′−i⁢D′′)q2*+i]×exp⁡[−q2−iq2+i⁢a1†⁢a2†]⁢|00⟩

where the relation between q2 and q1 are determined by Eq. (105) which resembles Eq. (97). This is just the new ABCD law for two-mode case in quantum optics.

9.3 Optical operators derived by decomposing GFO

9.3.1 GFO as quadratic combinations of canonical operators

In order to obtain the quadratic combinations of canonical operators, first derive an operator identity. Note Qi=(ai+ai†)/2,Pi=(ai−ai†)/(2⁢i), and Eqs. (29), (30), and we can prove the operator identity

(192)eλ2⁢[(Q1−Q2)2+(P1+P2)2]=11−λ:exp⁡(2⁢λ1−λ⁢K+):

where we have set

(193)K+≡14⁢[(Q1−Q2)2+(P1+P2)2]

When B=0,A=1,C→C/A,D=1, and using Eq. (179) we see that

(194)F2⁢(1,0,C/A)=exp⁡(i⁢CA⁢K+)

which is corresponding to the square phase operator in single-mode case. In a similar way, using Eqs. (37) and (39) we can derive another operator identity

(195)eλ2⁢[(Q1+Q2)2+(P1−P2)2]=11−λ:exp⁡[2⁢λ1−λ⁢K−]:

where

(196)K−=14⁢[(Q1+Q2)2+(P1−P2)2]

It then follows from Eqs. (179) and (195) that

(197)F2⁢(1,B/A,0)=exp⁡(−i⁢BA⁢K−)

which is corresponding to Fresnel propagator in free space (single-mode case). In particular, when B=C=0, and D=A−1, Eq. (179) becomes

(198)F2(A,0,0)=sechλ:exp[−a1†a2†tanhλ+(sechλ−1)×(a1†a1+a2†a2)+a1a2tanhλ]:

where A−A−1A+A−1=tanh⁡λ,A=eλ. Eq. (198) is just the two-mode squeezing operator,

(199)F2⁢(A,0,0)=exp⁡(−2⁢K0⁢ln⁡A)

(200)K0≡−i2⁢(Q1⁢P2+Q2⁢P1)

which actually squeezes the entangled state |ξ⟩ (its conjugate state is |η⟩),

(201)F2⁢(A,0,0)⁢|ξ⟩=∫d2⁢ξ′π⁢A⁢|ξ′/A⟩⁢⟨ξ′|ξ⟩=1A⁢|ξ/A⟩

Using the decomposition (117) of the matrix and combining equations (194), (197) and (199) together, we see that

(202)F2⁢(A,B,C)=exp⁡(i⁢CA⁢K+)⁢exp⁡(−2⁢K0⁢ln⁡A)×exp⁡(−i⁢BA⁢K−)

This is the two-mode quadratic canonical operator representation of F2⁢(A,B,C).

To prove Eq. (202), using Eqs. (201) and (40) we have

(203)⟨η|F2⁢(A,B,C)|ξ⟩=12⁢A⁢exp⁡(i⁢C2⁢A⁢|η|2−i⁢B2⁢A⁢|ξ|2)×exp⁡[12⁢A⁢(η*⁢ξ−η⁢ξ*)]

It then follows that

(204)⟨η′|F2⁢(A,B,C)|η⟩=12⁢i⁢B⁢exp⁡[i2⁢B⁢(A⁢|η|2−i⁢(η⁢η′⁣*+η*⁢η′)+D⁢|η′|2)]≡K2M⁢(η′,η)

which is just the transform kernel of a 2-dimensional GFT and the definition given in Eq. (202) is true.

Note that the quadratic combinations in Eqs. (193), (196) and (200) of the four canonical operators (Q1,Q2;P1,P2) obey the commutative relations [K+,K−]=2⁢K0,[K0,K±]=±K±, so F2⁢(A,B,C) involves a SU (2) Lie algebra structure. (This structure is also compiled by Q2/2,P2/2 and −i⁢(Q⁢P+P⁢Q)/2 that have been used in constructing F1⁢(A,B,C)).

9.3.2 Alternate decompositions of GFO and new optical operator identities

When A=D=0,B=1,C=−1, from Eq. (179) we see

(205)F2⁢(0,1,−1)=exp⁡[−i⁢π2⁢(a1†⁢a1+a2†⁢a2+1)]≡F

which can also be named the Fourier operator, since it induces the quantum mechanically transforms [71, 72]

(206)F†⁢Qi⁢F=Pi,F†⁢Pi⁢F†=−Qi

It then follows that

(207)F†⁢K+⁢F=K−

On the other hand, in order to obtain the decomposition of F2⁢(A,B,C) for A=0, similar to deriving Eq. (126), we have

(208)F2⁢(A,B,C)=exp⁡(−i⁢BD⁢K−)⁢exp⁡(2⁢K0⁢ln⁡D)×exp(i⁢CDK+),D≠0

While for B≠0 or C≠0, using Eqs. (127) and (128) we have another decomposition of F2⁢(A,B,C), i.e.,

(209)F2⁢(A,B,C)=exp⁡(i⁢DB⁢K+)⁢exp⁡(−2⁢K0⁢ln⁡B)⁢F×exp(i⁢ABK+),B≠0

and

(210)F2⁢(A,B,C)=exp⁡(−i⁢AC⁢K−)⁢exp⁡(−2⁢K0⁢ln⁡−1C)⁢F×exp(−i⁢DCK−),C≠0

In addition, noticing Eqs. (137) and (135), we can rewrite Eqs. (209) and (210) as follows:

(211)F2⁢(A,B,C)=exp⁡[iB⁢(D−1)⁢K+]⁢exp⁡(−i⁢B⁢K−)×exp[iB(A−1)K+],B≠0

and

(212)F2⁢(A,B,C)=exp⁡[−iC⁢(A−1)⁢K−]⁢exp⁡(i⁢C⁢K+)×exp[−iC(D−1)K−],C≠0

respectively.

Next, according to some optical systems used frequently in physical optics, we derive some new entangled optical operator identities. For a special optical system Eq. (131) with the parameter A=0,C=−B−1, which corresponds to the Fourier transform system, we have

(213)exp⁡(−i⁢BD⁢K−)⁢exp⁡(2⁢K0⁢ln⁡D)⁢exp⁡(−iB⁢D⁢K+)=exp⁡(i⁢DB⁢K+)⁢exp⁡(−2⁢K0⁢ln⁡B)⁢F

In particular, when A=D=0,C=−B−1, and with Eq. (133) corresponding to the ideal spectrum analyzer, we have

(214)exp⁡(−i⁢B⁢K−)⁢exp⁡(−iB⁢K+)⁢exp⁡(−i⁢B⁢K−)=exp⁡(−2⁢K0⁢ln⁡B)⁢F

When B=0,D=A−1,

(214)(A0CA−1)=(10C/A1)⁢(A00A−1)

which corresponds to the form of the image system, another operator identity is given by

(215)exp⁡(−2⁢K0⁢ln⁡A)⁢exp⁡(i⁢A⁢C⁢K+)=exp⁡(i⁢CA⁢K+)⁢exp⁡(−2⁢K0⁢ln⁡A)

When C=0,A=D−1,(D−1B0D)=(1B/D01)⁢(D−100D), which corresponds to the far foci system,

(216)exp⁡(i⁢DB⁢K+)⁢exp⁡(−2⁢K0⁢ln⁡B)⁢F⁢exp⁡(iB⁢D⁢K+)=exp⁡(−i⁢BD⁢K−)⁢exp⁡(2⁢K0⁢ln⁡D)

When D=0,C=−B−1, corresponding to the Fresnel transform system,

(216)(AB−B−10)=(B00B−1)⁢(01−10)⁢(10A/B1)

we have

(217)exp⁡(−iA⁢B⁢K+)⁢exp⁡(−2⁢K0⁢ln⁡A)⁢exp⁡(−i⁢BA⁢K−)=exp⁡(−2⁢K0⁢ln⁡B)⁢F⁢exp⁡(i⁢AB⁢K+)

The GFO can unify those optical operators in the two-mode case. Various decompositions of the GFO into the exponential canonical operators, corresponding to the decomposition of ray transfer matrix [A,B,C,D], are also derived. In our derivation, the entangled state representation has been fully used.

9.4 Quantum tomography and probability distribution for the Fresnel quadrature phase — two-mode entangled case

In Section 8 we have found that under the Fresnel transformation the pure position density |q⟩⁢⟨q| becomes the tomographic density |q⟩s,r⁢s,r⁢⟨q|, which is just the Radon transform of the Wigner operator Δ⁢(q,p). In this section we want to generalize the above conclusion to the two-mode entangled case. Here we shall prove

(218)F2|η⟩⟨η|F2†=|η⟩s,r⁢s,r⟨η|=π∫d2γd2σδ(η2−Dσ2+Bγ1)δ(η1−Dσ1−Bγ2)Δ(σ,γ)

That is, we show that |η⟩s,r⁢s,r⁢⟨η| is just the Radon transform of the entangled Wigner operator Δ⁢(σ,γ).

Similar to the single-mode case, operating F2⁢(r,s) on entangled state representation |η⟩, we obtain (F2⁢(r,s)⁢|η⟩≡|η⟩s,r)

(219)|η⟩s,r=1s*+r*exp[−s*−r*2⁢(s*+r*)|η|2+η⁢a1†s*+r*−η*⁢a2†s*+r*+s+rs*+r*a1†a2†]|00⟩

or

(220)|η⟩s,r=1D+i⁢Bexp[−A−i⁢C2⁢(D+i⁢B)|η|2+η⁢a1†D+i⁢B−η*⁢a2†D+i⁢B+D−i⁢BD+i⁢Ba1†a2†]|00⟩

where we have used the integration formula

(221)∫d2⁢zπ⁢exp⁡(ζ⁢|z|2+ξ⁢z+η⁢z*)=−1ζ⁢e−ξ⁢ηζ,Re⁢(ζ)<0

Noticing the completeness relation and the orthogonality of |η⟩, we immediately derive

(222)∫d2⁢ηπ|η⟩s,r⁢s,r⟨η|=1,s,r⟨η|η′⟩s,r=πδ(η−η′)δ(η*−η′⁣*)

a generalized entangled state representation |η⟩s,r with the completeness relation (222). From Eq. (220) we can see that

(223)a1⁢|η⟩s,r=(ηD+i⁢B+D−i⁢BD+i⁢B⁢a2†)⁢|η⟩s,r

(224)a2⁢|η⟩s,r=(−η*D+i⁢B+D−i⁢BD+i⁢B⁢a1†)⁢|η⟩s,r

so we have the eigen-equations for |η⟩s,r as follows:

(225)[D⁢(Q1−Q2)−B⁢(P1−P2)]⁢|η⟩s,r=2⁢η1⁢|η⟩s,r

(226)[B⁢(Q1+Q2)+D⁢(P1+P2)]⁢|η⟩s,r=2⁢η2⁢|η⟩s,r

We can also check Eqs. (223)–(226) by another way.

9.4.1 |η⟩s,r⁢s,r⁢⟨η| as Radon transform of the entangled Wigner operator

For a two-mode correlated system, we have introduced the Wigner operator in Eq. (43). According to the Wely correspondence rule [51]

(227)H⁢(a1†,a2†;a1,a2)=∫d2⁢γ⁢d2⁢σ⁢h⁢(σ,γ)⁢Δ⁢(σ,γ)

where h⁢(σ,γ) is the Weyl correspondence of H(a1†,a2†;a1,a2), and

(228)h⁢(σ,γ)=4⁢π2⁢Tr⁢[H⁢(a1†,a2†;a1,a2)⁢Δ⁢(σ,γ)]

The classical Weyl correspondence of the projection operator |η⟩r,s⁢r,s⁢⟨η| can be calculated as

(229)4⁢π2⁢Tr⁢[|η⟩r,s⁢r,s⁢⟨η|⁢Δ⁢(σ,γ)]=4⁢π2⁢∫d2⁢η′π3⁢⟨η|F2†|σ−η′⟩⁢⟨σ+η′|F2|η⟩×exp⁡(η′⁢γ*−η′⁣*⁢γ)

Then using Eq. (182) we have

(230)4⁢π2⁢Tr⁢[|η⟩s,r⁢s,r⁢⟨η|⁢Δ⁢(σ,γ)]=π⁢δ⁢(η2−D⁢σ2+B⁢γ1)⁢δ⁢(η1−D⁢σ1−B⁢γ2)

which means the following Weyl correspondence

(231)|η⟩s,r⁢s,r⁢⟨η|=π⁢∫d2⁢γ⁢d2⁢σ⁢δ⁢(η2−D⁢σ2+B⁢γ1)×δ⁢(η1−D⁢σ1−B⁢γ2)⁢Δ⁢(σ,γ)

so the projector operator |η⟩s,r⁢s,r⁢⟨η| is just the Radon transformation of Δ⁢(σ,γ), and D and B are the Radon transformation parameter. Combining Eqs. (219)–(231) together we complete the proof (218). Therefore, the quantum tomography in the two-mode entangled case is expressed as

(232)|s,r⟨η|ψ⟩|2=|⟨η|F†|ψ⟩|2=π∫d2γd2σδ(η2−Dσ2+Bγ1)δ(η1−Dσ1−Bγ2)⟨ψ|Δ(σ,γ)|ψ⟩

where ⟨ψ|Δ⁢(σ,γ)|ψ⟩ is the Wigner function. So the probability distribution for the Fresnel quadrature phase is the tomography (Radon transform of the two-mode Wigner function). This new relation between quantum tomography and optical Fresnel transform may provide experimentalists to figure out new approach for generating tomography.

Next we turn to the “frequency” domain, that is to say, we shall prove that the (A,C) related Radon transform of entangled Wigner operator Δ⁢(σ,γ) is just the pure state density operator |ξ⟩s,r⁢s,r⁢⟨ξ|, i.e.,

(233)F2⁢|ξ⟩⁢⟨ξ|⁢F2†=|ξ⟩s,r⁢s,r⁢⟨ξ|=π⁢∫δ⁢(ξ1−A⁢σ1−C⁢γ2)×δ⁢(ξ2−A⁢σ2+C⁢γ1)⁢Δ⁢(σ,γ)⁢d2⁢σ⁢d2⁢γ

where |ξ⟩ is the conjugated entangled state to |η⟩.

By analogy with the above procedure, we obtain the 2-dimensional Fresnel transformation in its “frequency domain”, i.e.,

(234)K2N⁢(ξ′,ξ)≡1π⁢⟨ξ′|F2⁢(r,s)|ξ⟩=12⁢i⁢(−C)⁢π×exp⁡[i2⁢(−C)⁢(D⁢|ξ|2+A⁢|ξ′|2−ξ′⁣*⁢ξ−ξ′⁢ξ*)]

where the superscript N means that this transform kernel corresponds to the parameter matrix N=[D,−C,−B,A]. Thus the 2D Fresnel transformation in its “frequency domain” is given by

(235)Ψ⁢(ξ′)=∫K2N⁢(ξ′,ξ)⁢Φ⁢(ξ)⁢d2⁢ξ

Operating F2⁢(r,s) on |ξ⟩ we also have

(236)|ξ⟩s,r=1A−i⁢Cexp[−D+i⁢B2⁢(A−i⁢C)|η|2+ξ⁢a1†A−i⁢C+ξ*⁢a2†A−i⁢C−A+i⁢CA−i⁢Ca1†a2†]|00⟩

or

(237)|ξ⟩s,r=1s*−r*exp[−s*+r*2⁢(s*−r*)|ξ|2+ξ⁢a1†s*−r*+ξ*⁢a2†s*−r*−s−rs*−r*a1†a2†]|00⟩

Noticing that the entangled Wigner operator in ⟨ξ| representation is expressed as

(238)Δ⁢(σ,γ)=∫d2⁢ξπ3⁢|γ+ξ⟩⁢⟨γ−ξ|⁢exp⁡(ξ*⁢σ−σ*⁢ξ)

and using the classical correspondence of |ξ⟩s,r⁢s,r⁢⟨ξ| which is calculated by

(239)h⁢(σ,γ)=4⁢π2⁢Tr⁢[|ξ⟩s,r⁢s,r⁢⟨ξ|⁢Δ⁢(σ,γ)]=π⁢δ⁢(ξ1−A⁢σ1−C⁢γ2)⁢δ⁢(ξ2−A⁢σ2+C⁢γ1)

we obtain

(240)|ξ⟩s,r⁢s,r⁢⟨ξ|=π⁢∫δ⁢(ξ1−A⁢σ1−C⁢γ2)×δ⁢(ξ2−A⁢σ2+C⁢γ1)⁢Δ⁢(σ,γ)⁢d2⁢σ⁢d2⁢γ

so the projector operator |ξ⟩s,r⁢s,r⁢⟨ξ| is another Radon transformation of the two-mode Wigner operator, with A and C being the Radon transformation parameter (“frequency domain”). Therefore, the quantum tomography in ⟨ξ|s,r representation is expressed as the Radon transformation of the Wigner function

(241)|⟨ξ|F2†|ψ⟩|2=|s,r⟨ξ|ψ⟩|2=π∫d2γd2σδ(ξ1−Aσ1−Cγ2)×δ(ξ2−Aσ2+Cγ1)⟨ψ|Δ(σ,γ)|ψ⟩

and ⟨ξ|=⟨ξ|F2†.s,r

9.4.2 Inverse Radon transformation

Now we consider the inverse Radon transformation. For instance, using Eq. (231) we see that the Fourier transformation of |η⟩s,r⁢s,r⁢⟨η| is

(242)∫d2⁢η⁢|η⟩s,r⁢s,r⁢⟨η|⁢exp⁡(−i⁢ζ1⁢η1−i⁢ζ2⁢η2)=π⁢∫d2⁢γ⁢d2⁢σ⁢Δ⁢(σ,γ)×exp⁡[−i⁢ζ1⁢(D⁢σ1+B⁢γ2)−i⁢ζ2⁢(D⁢σ2−B⁢γ1)]

The right-hand side of Eq. (242) can be regarded as a special Fourier transformation of Δ⁢(σ,γ), so by making its inverse Fourier transformation, we get

(243)Δ⁢(σ,γ)=1(2⁢π)4⁢∫−∞∞dr1⁢|r1|⁢∫−∞∞dr2⁢|r2|⁢∫0πdθ1⁢dθ2×∫−∞∞d2⁢ηπ|η⟩s,r⁢s,r⟨η|K(r1,r2,θ1,θ2)

where cos⁡θ1=cos⁡θ2=DB2+D2,r1=ζ1⁢B2+D2,r2=ζ2⁢B2+D2 and

(244)K⁢(r1,r2,θ1,θ2)≡exp⁡[−i⁢r1⁢(η1B2+D2−σ1⁢cos⁡θ1−γ2⁢sin⁡θ1)]×exp⁡[−i⁢r2⁢(η2B2+D2−σ2⁢cos⁡θ2+γ1⁢sin⁡θ2)]

Eq. (243) is just the inverse Radon transformation of entangled Wigner operator in the entangled state representation. This is different from the two independent Radon transformations’ direct product of the two independent single-mode Wigner operators, because in Eq. (220) the |η⟩s,r is an entangled state. Therefore the Wigner function of quantum state |ψ⟩ can be reconstructed from the tomographic inversion of a set of measured probability distributions |s,r⟨η|ψ⟩|2, i.e.,

(245)Wψ=1(2⁢π)4⁢∫−∞∞dr1⁢|r1|⁢∫−∞∞dr2⁢|r2|⁢∫0πdθ1⁢dθ2×∫−∞∞d2⁢ηπ|s,r⟨η|ψ⟩|2K(r1,r2,θ1,θ2)

Thus, based on the previous section, we have further extended the relation connecting optical Fresnel transformation with quantum tomography to the entangled case. The tomography representation ⟨η|=⟨η|F2†s,r is set up, based on which the tomogram of quantum state |ψ⟩ is just the squared modulus of the wave function ⟨η|ψ⟩s,r. In other words, the probability distribution for the Fresnel quadrature phase is the tomogram (Radon transform of the Wigner function).

10 Fractional Fourier transformation (FrFT) for a 1D case

The fractional Fourier transform (FrFT) has been shown to be a very useful tool in Fourier optics and information optics. The concept of FrFT was first introduced mathematically in 1980 by Namias [36] as a mathematical tool to solve theoretical physical problems [95], but did not bring enough attention until Mendlovic and Ozaktas [37–39] defined the α-th FrFT physically, based on propagation in quadratic graded-index media (GRIN media with medium parameters n⁢(r)=n1−n2⁢r2/2). Since then a lot of work have been done on its properties, optical implementations and applications [96–100].

10.1 Quantum version of FrFT

The FrFT of θ-order is defined in such a manner as

(246)Fθ⁢[f⁢(x)]=ei⁢(π2−θ)2⁢π⁢sin⁡θ⁢∫−∞∞×exp⁡[−i⁢x2+y22⁢tan⁡θ+i⁢x⁢x′sin⁡θ]⁢f⁢(x)⁢d⁢x

where the exponential function is an integral kernel. In order to find the quantum correspondence of FrFT, multiply the function exp⁡(−i⁢x2+y22⁢tan⁡θ+i⁢x⁢x′sin⁡θ) by the ket ∫dy⁢|y⟩ and bra ∫dx⁢⟨x| from left and right, respectively, where |y⟩ and |x⟩ are coordinate eigenvectors, X⁢|x⟩=x⁢|x⟩, and then use Eq. (13) and the IWOP technique to perform the integration. We obtain

(247)∫−∞∞dx⁢dy⁢|y⟩⁢exp⁡(−i⁢x2+y22⁢tan⁡θ+i⁢x⁢ysin⁡θ)⁢⟨x|=−2⁢π⁢i⁢sin⁡θ⁢ei⁢θ⁢exp⁡(i⁢θ⁢a†⁢a)

where we have used the operator identity exp(fa†a)=:exp[(ef−1)a†a]: in the last step of Eq. (247).

From the orthogonal relation ⟨x′|x⟩=δ⁢(x−x′), we know that Eq. (247) indicates

(248)ei⁢(π2−θ)2⁢π⁢sin⁡θ⁢exp⁡(−i⁢x2+y22⁢tan⁡θ+i⁢x⁢ysin⁡θ)=⟨y|ei⁢θ⁢a†⁢a|x⟩

which implies that the integral kernel in Eq. (246) is just the matrix element of operator exp⁡(i⁢θ⁢a†⁢a) in coordinate state exp⁡(i⁢θ⁢a†⁢a) called as Fractional Fourier Operator [70]. Therefore, if we consider f⁢(x) as ⟨x|f⟩, the wave function of quantum state |f⟩ in the coordinate representation, from Eqs. (246) and (248) it then follows that

(249)Fθ⁢[f⁢(x)]=∫−∞∞dx⁢⟨y|ei⁢θ⁢a†⁢a|x⟩⁢f⁢(x)=⟨y|ei⁢θ⁢a†⁢a|f⟩≡g⁢(y)

which suggests

(250)|g⟩=ei⁢θ⁢a†⁢a⁢|f⟩

From Eqs. (249) and (246) one can see that the FrFT in Eq. (246) corresponds actually to the rotating operator (ei⁢θ⁢a†⁢a) transform in Eq. (249) between two quantum states, which is just the quantum version of FrFT.

In fact, using quantum version of FrFT, one can directly derive various properties of the FrFTs. An important feature of the FrFT is that they are composed according to Fθ′⁢Fθ=Fθ′+θ (the additivity property). Without losing generality, we examine

(251)Fθ+θ′⁢[f⁢(x)]≡∫−∞∞d2⁢ηπ⁢⟨y|ei⁢(θ+θ′)⁢a†⁢a|x⟩⁢f⁢(x)

According to the completeness relation of coordinate eigenvector, ∫−∞∞dx′⁢|x′⟩⁢⟨x′|=1, Eq. (251) yields

(252)Fθ+θ′⁢[f⁢(x)]=∫−∞∞dx′⁢⟨y|ei⁢θ⁢a†⁢a|x′⟩⁢∫−∞∞dx⁢⟨x′|ei⁢θ′⁢a†⁢a|x⟩⁢f⁢(x)=∫−∞∞dx′⁢⟨y|ei⁢θ⁢a†⁢a|x′⟩⁢Fθ⁢[f⁢(x)]=Fθ⁢Fθ′⁢[f⁢(x)]

which is just the additivity of FrFT.

In particular, when |f⟩ is the number state, |f⟩=|n⟩=a†nn!⁢|0⟩, its wavefunction in coordinate representation is

(253)f⁢(x)=⟨x|n⟩=12n⁢n!⁢π⁢e−x2/2⁢Hn⁢(x)

The FrFT of ⟨x|n⟩ is Fθ⁢[⟨x|n⟩]=⟨y|ei⁢θ⁢a†⁢a|n⟩=ei⁢n⁢θ⁢⟨y|n⟩, or

(254)Fθ⁢[e−x2/2⁢Hn⁢(x)]=ei⁢n⁢θ⁢e−y2/2⁢Hn⁢(y)

which indicates that the eigenfunction is Hermite-Gaussian function with the corresponding eigenvalue being ei⁢n⁢θ.

10.2 On the scaled FrFT operator

In studying various optical transformations the optical operator method is proposed [77–79] as mapping of ray-transfer ABCD matrix, such that the ray transfer through optical instruments and the diffraction can be discussed by virtue of the commutative relations of operators and the matrix algebra. The square phase operators, scaling operator, Fourier transform operator and the propagation operator in free space have been proposed in the literature. Two important questions thus naturally arise: (i) What is the scaled FrFT (SFrFT) operator which corresponds to the SFrFT’s integration kernel [36]

(255)12⁢π⁢i⁢fe⁢sin⁡ϕ⁢exp⁡[i⁢(x2+x′⁣2)2⁢fe⁢tan⁡ϕ−i⁢x′⁢xfe⁢sin⁡ϕ]

where fe is standard focal length (or a scaled parameter)? (ii) If this operator is found, can it be further decomposed into simpler operators and what are their physical meaning? Since SFrFT has wide application in optical information detection and can be implemented even by using a thick lens [44], so our questions are worth paying attention to Ref. [102].

Let us start with a thick lens (shown in Fig. 1) which represents a transfer matrix [44]

(256)(ABCD)=(1−(1−1/n)⁢lR1ln−[(n−1)⁢R1+R2R1⁢R2−l⁢(n−1)2n⁢R1⁢R2]1−(1−1/n)⁢lR2)

where n is the reflective index; l is the thickness of thick lens; R1 and R2 denotes the curvature radius of the two surfaces of the lens, respectively. When we choose R1=R2=R, Eq. (256) is reduced to

(257)(ABCD)=(1−(1−1/n)⁢lRln−[(n−1)⁢2R−l⁢(n−1)2n⁢R2]1−(1−1/n)⁢lR)

By defining 1−(1−1/n)⁢lR=cos⁡ϕ,ln=fe⁢sin⁡ϕ, and lR=n⁢(1−cos⁡ϕ)n−1,l=n⁢fe⁢sin⁡ϕ, we can recast (257) into the simple form

(258)(ABCD)=(cos⁡ϕfe⁢sin⁡ϕ−sin⁡ϕ/fecos⁡ϕ)det(ABCD)=1

According to Eq. (98) we immediately know that the operator of SFrFT is

(259)F1⁢(A,B,C)=exp⁡[i⁢(fe−1/fe)⁢tan⁡ϕ2⁢V⁢a†2]×exp⁡[(a†⁢a+12)⁢ln⁡2⁢ϕV]×exp⁡[i⁢(fe−1/fe)⁢tan⁡ϕ2⁢V⁢a2]V=[2+i⁢(fe+1/fe)⁢tan⁡ϕ]

Noting that the matrix (ABCD) can be decomposed into

(260)(ABCD)=(10−1fe⁢tan⁡ϕ1)⁢(cos⁡ϕ00cos⁡ϕ)×(1fe⁢tan⁡ϕ01)

according to the previous section, we have

(261)F1⁢(A,B,C)=exp⁡(tan⁡ϕ2⁢i⁢fe⁢Q2)×exp⁡[−i2⁢(Q⁢P+P⁢Q)⁢ln⁡cos⁡ϕ]×exp⁡(fe⁢tan⁡ϕ2⁢i⁢P2)

where Q=(a+a†)/2,P=(a−a†)/(2⁢i) and exp⁡(−i⁢tan⁡ϕ2⁢fe⁢Q2),exp⁡[−i2⁢(Q⁢P+P⁢Q)⁢ln⁡cos⁡ϕ] and exp⁡(−i⁢fe⁢tan⁡ϕ2⁢P2) are the quadrature phase operator, the squeezing operator and the free propagation operator, respectively. On the other hand, from (ABCD)−1=(D−B−CA), and Eq. (261) we see

(262)F1−1⁢(A,B,C)=exp⁡(−i⁢C2⁢D⁢Q2)×exp⁡[−i2⁢(Q⁢P+P⁢Q)⁢ln⁡D]×exp⁡(i⁢B2⁢D⁢P2)

It then follows that

(263)F1⁢(A,B,C)=exp⁡(fe⁢tan⁡ϕ2⁢i⁢P2)×exp⁡[i2⁢(Q⁢P+P⁢Q)⁢ln⁡cos⁡ϕ]×exp⁡(tan⁡ϕ2⁢i⁢fe⁢Q2)

Using the canonical operator form (261) or (263) of F1⁢(A,B,C) we can deduce its matrix element in the coordinate states |x⟩ (its conjugate state is |p⟩)

(264)⟨x′|F1⁢(A,B,C)|x⟩=12⁢π⁢i⁢fe⁢sin⁡ϕ⁢exp⁡[i⁢(x2+x′⁣2)2⁢fe⁢tan⁡ϕ−i⁢x′⁢xfe⁢sin⁡ϕ]

which is just the kernel of SFrFT. Thus we name F1⁢(A,B,C) SFrFT operator. Noticing that the Q2/2,P2/2 and i4⁢(Q⁢P+P⁢Q) construct a close SU (2) Lie algebra, we can put Eq. (263) into a more compact form, i.e.,

(265)F1⁢(A,B,C)=exp⁡[−i⁢ϕ⁢fe2⁢(P2+Q2fe2)]

Eqs. (261), (263) and (265) are different forms of the same operator of SFrFT; especially when fe=1,F1⁢(A,B,C)→exp⁡(−i⁢ϕ⁢a†⁢a), which is the usual FrFT operator.

Using Eq. (264) the SFrFT of f⁢(x)=⟨x|f⟩, denoted as Ffeϕ⁢[f⁢(x)], can be expressed as an matrix element in quantum optics context,

(266)Ffeϕ⁢[f⁢(x)]=∫dx⁢⟨x′|F1⁢(A,B,C)|x⟩⁢⟨x|f⟩=⟨x′|F1⁢(A,B,C)|f⟩

The above discussions are useful since any unimodular matrix (ABCD) can be decomposed into [103]

(267)(ABCD)=(10−P1)⁢(m00m−1)×(cos⁡ϕfe⁢sin⁡ϕ−sin⁡ϕ/fecos⁡ϕ)

where the parameters m,P,ϕ are all real,

(268)m2=A2+B2fe2,tan⁡ϕ=BA⁢fe,P=−A⁢C+D⁢B/fe2A2+B2fe2

Correspondingly, the operator of FrFT is given by

(269)F1⁢(A,B,C)=exp⁡(−i2⁢P⁢Q2)×exp⁡[−i2⁢(Q⁢P+P⁢Q)⁢ln⁡m]×exp⁡[ϕ⁢fe2⁢i⁢(P2+Q2fe2)]

where F1⁢(1,0,−P)=exp⁡(−i2⁢P⁢Q2) is the quadratic phase operator. Thus the general Fresnel transform can always be expressed by SFrFT as follows:

(270)g⁢(x′)=∫dy⁢⟨x′|F1⁢(1,0,−P)⁢F1⁢(m,0,0)|y⟩×∫dx⟨y|F1(A,B,C)|x⟩f(x)=1m⁢exp⁡(−i2⁢P⁢x′⁣2)⁢Ffeϕ⁢(f)⁢(x′m)

i.e., the output g⁢(x′) is the SFrFT of the input f⁢(x) plus a quadratic phase term exp⁡(−i2⁢P⁢x′⁣2).

10.3 An integration transformation from Chirplet to FrFT kernel

In the history of the development of optics we have known that each optical setup corresponds to an optical transformation. For example, thick lens as a fractional Fourier transformer. In turn, once a new integration transform is found, its experimental implementation is expected. In this subsection we report a new integration transformation which can convert chirplet function to FrFT kernel [104]. As this new transformation is invertible and obeys Parseval theorem, we expect it be realized by experimentalists.

The new transform we propose here is

(271)∬−∞∞⁢d⁢p⁢d⁢qπ⁢e2⁢i⁢(p−x)⁢(q−y)⁢h⁢(p,q)≡f⁢(x,y)

which differs from the usual two-fold Fourier transformation ∬−∞∞⁢d⁢x⁢d⁢y4⁢π2⁢ei⁢p⁢x+i⁢q⁢y⁢f⁢(x,y). In particular, when h⁢(p,q)=1, Eq. (271) is reduced to

(272)∬−∞∞⁢d⁢p⁢d⁢qπ⁢e2⁢i⁢(p−x)⁢(q−y)=1

so e2⁢i⁢(p−x)⁢(q−y) can be considered a basis funtion in p−q phase space, or Eq. (271) can be looked as an expansion of f⁢(x,y) with the expansion coefficient being h⁢(p,q). We can prove that the reciprocal transformation of Eq. (271) is

(273)∬−∞∞⁢d⁢x⁢d⁢yπ⁢e−2⁢i⁢(p−x)⁢(q−y)⁢f⁢(x,y)=h⁢(p,q)

In fact, substituting Eq. (271) into the left-hand side of Eq. (273), we get

(274)∬−∞∞⁢d⁢p′⁢d⁢q′π⁢h⁢(p′,q′)⁢∬⁢d⁢x⁢d⁢yπ×e2⁢i⁢[(p′−x)⁢(q′−y)−(p−x)⁢(q−y)]=∬−∞∞⁢d⁢p′⁢d⁢q′⁢h⁢(p′,q′)⁢e2⁢i⁢(p′⁢q′−p⁢q)×δ(p−p′)δ(q−q′)=h(p,q)

This transformation’s Parseval-like theorem can be proved

(275)∬−∞∞⁢d⁢p⁢d⁢qπ⁢|h⁢(p,q)|2=∬⁢d⁢x⁢d⁢yπ⁢|f⁢(x,y)|2⁢∬⁢d⁢x′⁢d⁢y′π⁢e2⁢i⁢(x′⁢y′−x⁢y)×∬−∞∞⁢d⁢p⁢d⁢qπ⁢e2⁢i⁢[(−y′⁢p−x′⁢q)+(p⁢y+x⁢q)]=∬⁢d⁢x⁢d⁢yπ⁢|f⁢(x,y)|2

Now we apply Eq. (271) to phase space transformation in quantum optics. Recall that a signal ψ⁢(q)’s Wigner transform [86, 90–94, 105] is

(276)ψ⁢(q)→∫d⁢u2⁢π⁢ei⁢p⁢u⁢ψ∗⁢(q+u2)⁢ψ⁢(q−u2)

Use Dirac’s symbol [7] to write ψ⁢(q)=⟨q|ψ⟩. |q⟩ is the eigenvector of coordinate Q, and the Wigner operator emerges from Eq. (276),

(277)12⁢π⁢∫−∞∞du⁢e−i⁢p⁢u⁢|q−u2⟩⁢⟨q+u2|=Δ⁢(p,q),ℏ=1

If h⁢(q,p) is quantized as the operator H^⁢(P,Q) through the Weyl–Wigner correspondence [51]

(278)H⁢(P,Q)=∬−∞∞⁢d⁢p⁢d⁢q⁢Δ⁢(p,q)⁢h⁢(q,p)

then

(279)h⁢(q,p)=∫−∞∞du⁢e−i⁢p⁢u⁢⟨q+u2|H^⁢(Q,P)|q−u2⟩

which in the literature is named the Weyl transform. h⁢(q,p) is the Weyl classical correspondence of the operator H^⁢(Q,P). Substituting Eq. (279) into Eq. (271), we have

(280)∬−∞∞⁢d⁢p⁢d⁢qπ⁢e2⁢i⁢(p−x)⁢(q−y)⁢h⁢(p,q)=∫−∞∞du⁢e−i⁢x⁢u⁢⟨y+u|H^⁢(Q,P)|y⟩

Using ⟨y+u|=⟨u|⁢ei⁢P⁢y and (2⁢π)−1⁢e−i⁢x⁢u=⟨p=x|u⟩, where ⟨p| is the momentum eigenvector, and

(280)∫−∞∞du⁢e−i⁢x⁢u⁢⟨y+u|=2⁢π⁢⟨p=x|⁢ei⁢x⁢y

then Eq. (280) becomes

(281)∬−∞∞⁢d⁢p⁢d⁢qπ⁢e2⁢i⁢(p−x)⁢(q−y)⁢h⁢(p,q)=2⁢π⁢⟨p=x|H^⁢(Q,P)|y⟩⁢ei⁢x⁢y

Thus through the new integration transformation a new relationship between a phase space function h⁢(p,q) and its Weyl–Wigner correspondence operator H^⁢(Q,P) has been revealed. The inverse of Eq. (281), according to Eq. (273), is

(282)∬−∞∞⁢d⁢x⁢d⁢yπ/2⁢e−2⁢i⁢(p−x)⁢(q−y)⁢⟨p=x|H^⁢(Q,P)|y⟩⁢ei⁢x⁢y=h⁢(p,q)

For example, when H^⁢(Q,P)=ef⁢(P2+Q2−1)/2, its classical correspondence is

(283)ef⁢(P2+Q2−1)/2→h⁢(p,q)=2ef+1⁢exp⁡[2⁢ef−1ef+1⁢(p2+q2)]

Substituting Eq. (283) into Eq. (281), we have

(284)2ef+1⁢∬−∞∞⁢d⁢p⁢d⁢qπ⁢e2⁢i⁢(p−x)⁢(q−y)⁢exp⁡[2⁢ef−1ef+1⁢(p2+q2)]=2⁢π⁢⟨p=x|ef⁢(P2+Q2−1)/2|y⟩⁢ei⁢x⁢y

Using the Gaussian integration formula

(285)∬−∞∞⁢d⁢p⁢d⁢qπ⁢e2⁢i⁢(p−x)⁢(q−y)⁢e−λ⁢(p2+q2)=1λ2+1⁢exp⁡[−λ⁢(x2+y2)λ2+1+2⁢i⁢λ2λ2+1⁢x⁢y]

in particular, when λ=−i⁢tan⁡(π4−α2), with −λλ2+1=i2⁢tan⁡α,2⁢λ2λ2+1=1−1sin⁡α, Eq. (285) becomes

(286)2ie−i⁢α+1⁢∬−∞∞⁢d⁢p⁢d⁢qπ⁢e2⁢i⁢(p−x)⁢(q−y)×exp⁡[i⁢(p2+q2)⁢tan⁡(π4−α2)]=1ie−i⁢α⁢sin⁡α⁢exp⁡[i⁢(x2+y2)2⁢tan⁡α−i⁢x⁢ysin⁡α]⁢ei⁢x⁢y

where exp⁡[i⁢tan⁡(π4−α2)⁢(p2+q2)] represents an infinite long chirplet function. Comparing Eq. (286) with Eq. (284) we see ie−i⁢α=ef,f=i⁢(π2−α), it then follows that

(287)⟨p=x|ei⁢(π2−α)⁢(P2+Q2−1)/2|y⟩=12⁢π⁢ie−i⁢α⁢sin⁡α⁢exp⁡[i⁢(x2+y2)2⁢tan⁡α−i⁢x⁢ysin⁡α]

where the right-hand side of Eq. (287) is just the FrFT kernel. Therefore the new integration transformation (271) can convert spherical wave to FrFT kernel. We expect that this transformation could be implemented by experimentalists.

Moreover, this transformation can also be used to solve some operator ordering problems. We notice

(288)1π⁢exp⁡[2⁢i⁢(p−x)⁢(q−y)]=∫−∞∞d⁢v2⁢π⁢δ⁢(q−y−v2)⁢exp⁡[i⁢(p−x)⁢v]

so the transformation (271) is equivalent to

(289)h⁢(p,q)→∬−∞∞⁢d⁢p⁢d⁢qπ⁢e2⁢i⁢(p−x)⁢(q−y)⁢h⁢(p,q)=∬−∞∞⁢d⁢p⁢d⁢q2⁢π⁢h⁢(p+x,y+q2)⁢ei⁢p⁢q

For example, using Eqs. (277) and (288), we have

(290)Δ⁢(p,q)→∬−∞∞⁢d⁢p⁢d⁢q2⁢π⁢Δ⁢(p+x,y+q2)⁢ei⁢p⁢q=δ⁢(y−Q)⁢δ⁢(x−P)

so

(291)1π⁢∬⁢d⁢p′⁢d⁢q′⁢Δ⁢(q′,p′)⁢e2⁢i⁢(p−p′)⁢(q−q′)=δ⁢(q−Q)⁢δ⁢(p−P)

Thus this new transformation can convert the Wigner operator to δ⁢(q−Q)⁢δ⁢(p−P). Similarly, we have

(291)1π⁢∬⁢d⁢p′⁢d⁢q′⁢Δ⁢(q′,p′)⁢e−2⁢i⁢(p−p′)⁢(q−q′)=δ⁢(p−P)⁢δ⁢(q−Q)

Then for the Wigner function of a density operator ρ,Wψ⁢(p,q)≡Tr⁢[ρ⁢Δ⁢(p,q)], we have

(292)∬−∞∞⁢d⁢p′⁢d⁢q′π⁢Tr⁢[ρ⁢Δ⁢(p′,q′)]⁢e2⁢i⁢(p−p′)⁢(q−q′)=Tr⁢[ρ⁢δ⁢(q−Q)⁢δ⁢(p−P)]=∫d⁢u⁢d⁢v4⁢π2⁢Tr⁢[ρ⁢ei⁢(q−Q)⁢u⁢ei⁢(p−P)⁢v]

We may define Tr⁢[ρ⁢ei⁢(q−Q)⁢u⁢ei⁢(p−P)⁢v] as the Q–Pcharacteristic function. Similarly,

(293)∬−∞∞⁢d⁢p′⁢d⁢q′π⁢Tr⁢[ρ⁢Δ⁢(p′,q′)]⁢e−2⁢i⁢(p−p′)⁢(q−q′)=Tr⁢[ρ⁢δ⁢(p−P)⁢δ⁢(q−Q)]=∫d⁢u⁢d⁢v4⁢π2⁢Tr⁢[ρ⁢ei⁢(p−P)⁢v⁢ei⁢(q−Q)⁢u]

We name Tr⁢[ρ⁢ei⁢(p−P)⁢v⁢ei⁢(q−Q)⁢u] as the P–Q characteristic function.

11 Complex fractional Fourier transformation

In this section, we extend 1D FrFT to the complex fractional Fourier transformation (CFrFT).

11.1 Quantum version of CFrFT

According to Ref. [106], based on the entangled state |η⟩ in two-mode Fock space and its orthonormal property, we can take the matrix element of exp⁡[−i⁢α⁢(a1†⁢a1+a2†⁢a2)] in the entangled state |η⟩,

(294)KF⁢(η′,η)=⟨η′|exp⁡[−i⁢α⁢(a1†⁢a1+a2†⁢a2)]|η⟩

as the integral transform kernel of CFrFT,

(295)Fα⁢[f]⁢(η)=∫d2⁢ηπ⁢KF⁢(η′,η)⁢f⁢(η)

Using the normally ordered expansion of eλ⁢a1†⁢a1=:exp[(eλ−1)a1†a1]: and the completeness relation of the coherent state representation, |zi⟩=exp⁡(−12⁢|zi|2+zi⁢ai†)⁢|0⟩i, we calculate that KF⁢(η′,η) is

(296)KF⁢(η′,η)=ei⁢(α−π2)2⁢sin⁡α⁢exp⁡[i⁢(|η′|2+|η|2)2⁢tan⁡α−i⁢(η′⁣∗⁢η+η∗⁢η′)2⁢sin⁡α]

which is just the integral kernel of the CFrFT in Ref. [107]. Thus we see that the matrix element of exp⁡[−i⁢α⁢(a1†⁢a1+a2†⁢a2)] between two entangled state representations |η⟩ and |η′⟩ corresponds to CFrFT. This is a new route from quantum optical transform to classical CFrFT transform. Let η=x2+i⁢y2,η′=x1+i⁢y1, Eq. (295) becomes

(297)Fα⁢[f]⁢(x2,y2)=ei⁢(α−π2)2⁢sin⁡α⁢exp⁡[i⁢(x22+y22)2⁢tan⁡α]×∫d⁢x1⁢d⁢y1πexp[i⁢(x12+y12)2⁢tan⁡α−i(x1⁢x2+y1⁢y2)sin⁡α]×f⁢(x1,y1)

In fact, let f⁢(η)=⟨η|f⟩. Then use Eqs. (30) and (296), and we have

(298)⟨η′|exp⁡[−i⁢α⁢(a1†⁢a1+a2†⁢a2)]|f⟩=∫d2⁢ηπ⁢KF⁢(η′,η)⁢f⁢(η)=ei⁢(α−π2)2⁢sin⁡α⁢∫d2⁢ηπ×exp⁡[i⁢(|η′|2+|η|2)2⁢tan⁡α−i⁢(η′⁣∗⁢η+η∗⁢η′)2⁢sin⁡α]⁢f⁢(η)

Thus the quantum mechanical version of CFrFT is given by

(299)Fα⁢[f]⁢(η′)≡⟨η′|exp⁡[−i⁢α⁢(a1†⁢a1+a2†⁢a2)]|f⟩

The standard complex Fourier transform is Fπ/2. F0 is the identity operator.

11.2 Additivity property and eigenmodes of CFrFT

We will show later that CFrFT can help us reveal some new property which has been overlooked in the formulation of the direct product of two real FrFTs [108]. The definition (309) is of course required to satisfy the basic postulate that Fα⁢Fβ⁢[f]⁢(η′)=Fα+β⁢[f⁢(η)] (the additivity property). For this purpose, using Eq. (309) and Eq. (30), we see

(300)Fα+β⁢[f⁢(η)]≡⟨η′|e−i⁢(α+β)⁢(a1†⁢a1+a2†⁢a2)|f⟩=∫−∞∞d2⁢η′′π⁢⟨η′|e−i⁢α⁢(a1†⁢a1+a2†⁢a2)|η′′⟩⁢Fβ⁢[f⁢(η)]=Fα⁢Fβ⁢[f⁢(η)]

This derivation is clear and concise by employing the |η⟩ representation and quantum mechanical version of CFrFT.

On the other hand, the formula (299) can help us to derive CFrFT of some wave functions easily. For example, when |f⟩ is a two-mode number state |m,n⟩=a1†m⁢a2†n/m!⁢n!⁢|00⟩, then the CFrFT of the wave function ⟨η|m,n⟩ is

(301)Fα⁢[⟨η|m,n⟩]=⟨η′|ei⁢(α+β)⁢(a1†⁢a1+a2†⁢a2)|m,n⟩=ei⁢(α+β)⁢(m+n)⁢⟨η′|m,n⟩

To calculate ⟨η′|m,n⟩, let us recall the definition of two-variable Hermite polynomial Hm,n⁢(ξ,ξ∗), and we can expand ⟨η′| as

(302)⟨η′|=⟨00|⁢∑m,n=0∞im+n⁢a1m⁢a2nm!⁢n!⁢Hm,n⁢(−i⁢η′⁣∗,i⁢η′)⁢e−|η′|2/2

thus

(303)⟨η|′m,n⟩=im+nm!⁢n!Hm,n(−iη′⁣∗,iη)′e−|η|2′/2

As a result of Eq. (303), we see that Eq. (301) becomes

(304)Fα⁢[Hm,n⁢(−i⁢η∗,i⁢η)⁢e−|η|2/2]=ei⁢(α+β)⁢(m+n)Hm,n(−iη′⁣∗,iη)′e−|η|2′/2

If we consider the operation Fα as an operator, one can say that the eigenfunction of Fα (the eigenmodes of CFrFT) is the two-variable Hermite polynomials Hm,n with the eigenvalue being ei⁢(α+β)⁢(m+n). This is a new property of CFrFT. Since the function space spanned by Hm,n⁢(η,η∗) is complete,

(305)∫d2⁢ηπ⁢e−|η|2⁢Hm,n⁢(η,η∗)⁢[Hm,n⁢(η,η∗)]∗=m!n!m!′n!′⁢δm,m′⁢δn,n′

and

(306)∑m,n=0∞1m!⁢n!Hm,n(η,η∗)[Hm,n(η,′η′⁣∗)]∗e−|η|2=πδ(η−η)′δ(η∗−η′⁣∗)

one can confirm that the eigenmodes of CFrFT form an orthogonal and complete basis set [109]. Note that the two variable Hermite polynomial Hm,n⁢(η,η∗) is not the direct product of two independent ordinary Hermite polynomials, so CFrFT differs from the direct product of two FrFTs.

11.3 From chirplet to CFrFT kernel

In this subsection, by developing Eq. (271) to a more general case which can be further related to the transformation between two mutually conjugate entangled state representations |ξ⟩ and |η⟩, we shall propose a new integration transformation in ξ–η phase space [see Eq. (307) below] and its inverse transformation. We find that Eq. (307) also possesses some well-behaved transformation properties and can be used to obtain the CFrFT kernel from a chirplet [110].

11.3.1 New complex integration transformation

Corresponding to the structure of phase space spanned by |ξ⟩ and |η⟩ and enlightened by Eq. (271), we propose a new complex integration transformation in ξ–η phase space

(307)∫d2⁢ξ⁢d2⁢ηπ2⁢e(ξ−μ)⁢(η∗−ν∗)−(η−ν)⁢(ξ∗−μ∗)⁢F⁢(η,ξ)≡D⁢(ν,μ)

When F⁢(ξ,η)=1, Eq. (307) becomes

(308)∫d2⁢ξ⁢d2⁢ηπ2⁢e(ξ−μ)⁢(η∗−ν∗)−(η−ν)⁢(ξ∗−μ∗)=∫d2⁢ξ⁢δ⁢(ξ−μ)⁢δ⁢(ξ∗−μ∗)⁢eν⁢(ξ∗−μ∗)−ν∗⁢(ξ−μ)=1

so e(ξ−μ)⁢(η∗−ν∗)−(η−ν)⁢(ξ∗−μ∗)can be considered a basis function in ξ−η phase space, or Eq. (307) can be looked as an expansion of D⁢(ν,μ) in terms of e(ξ−μ)⁢(η∗−ν∗)−(η−ν)⁢(ξ∗−μ∗), with the expansion coefficient being F⁢(η,ξ).

We can prove that the inverse transform of Eq. (307) is

(309)∫d2⁢μ⁢d2⁢νπ2⁢e(ξ∗−μ∗)⁢(η−ν)−(η∗−ν∗)⁢(ξ−μ)⁢D⁢(ν,μ)≡F⁢(η,ξ)

In fact, substituting Eq. (307) into the left-hand side of Eq. (309), and we obtain

(310)∫d2ξd′2η′π2F(η,′ξ)′∫d2μd2νπ2×e(ξ−′μ)(η′∗−ν∗)−(η−′ν)(ξ′∗−μ∗)+(ξ∗−μ∗)(η−ν)−(η∗−ν∗)(ξ−μ)=∫d2ξd′2ηF′(η,′ξ)′e(ξ′η′∗−ηξ′′∗+ξ∗η−η∗ξ)×δ(2)(η−′η)δ(2)δ(ξ−ξ)′=F(η,ξ)

The Parseval-like theorem for this transformation can also be demonstrated as

(311)∫d2⁢ξ⁢d2⁢ηπ2⁢|F⁢(η,ξ)|2=∫d2⁢μ⁢d2⁢νπ2⁢|D⁢(ν,μ)|2

11.3.2 Complex integration transformation and complex Weyl transformation

In Ref. [53, 54] for correlated two-body systems, we have successfully established the so-called entangled Wigner operator, expressed in the entangled state ⟨η| representation as Eq. (43),

(312)Δ⁢(σ,γ)→Δ⁢(η,ξ)=∫d2⁢σπ3⁢|η−σ⟩⁢⟨η+σ|⁢eσ⁢ξ∗−σ∗⁢ξ

The advantage of introducing Δ⁢(η,ξ) can be seen in Ref. [80]. The corresponding Wigner function for a density matrix ρ is

(313)Wρ⁢(η,ξ)=∫d2⁢ηπ3⁢⟨η+σ|ρ|η−σ⟩⁢eσ⁢ξ∗−σ∗⁢ξ

If F⁢(η,ξ) is quantized as the operator F⁢(Q1,Q2,P1,P2) through the Weyl–Wigner correspondence

(314)F⁢(Q1,Q2,P1,P2)=∫d2⁢η⁢d2⁢ξ⁢F⁢(η,ξ)⁢Δ⁢(η,ξ)

then using Eq. (313) we see

(315)F⁢(η,ξ)=4⁢π2⁢Tr⁢[F⁢(Q1,Q2,P1,P2)⁢Δ⁢(η,ξ)]=4⁢∫d2⁢σπ⁢eσ⁢ξ∗−σ∗⁢ξ⁢⟨η+σ|F⁢(Q1,Q2,P1,P2)|η−σ⟩

which is named as the complex Weyl transform, and F⁢(η,ξ) is the Weyl classical correspondence of F⁢(Q1,Q2,P1,P2). Substituting Eq. (315) into Eq. (307), we get

(316)∬⁢d2⁢ξ⁢d2⁢ηπ2⁢e(ξ−μ)⁢(η∗−ν∗)−(η−ν)⁢(ξ∗−μ∗)⁢F⁢(η,ξ)=4⁢∫d2⁢σπ⁢eμ∗⁢σ−μ⁢σ∗⁢⟨ν+2⁢σ|F⁢(Q1,Q2,P1,P2)|ν⟩

Using Eq. (29), we have

(317)⟨ν+2⁢σ|=⟨2⁢σ|⁢exp⁡{i2⁢[ν1⁢(P1−P2)−ν2⁢(Q1+Q2)]}ν=ν1+i⁢ν2

As a result of Eq. (317) and 12⁢eμ∗⁢σ−μ⁢σ∗=⟨ξ=μ|2⁢σ⟩, we see

(318)4⁢∫d2⁢σ⁢eμ∗⁢σ−μ⁢σ∗⁢⟨ν+2⁢σ|=2⁢π⁢⟨ξ=μ|⁢ei⁢(μ2⁢ν1−μ1⁢ν2)

Using Eq. (318), we convert Eq. (316) as

(319)∬⁢d2⁢ξ⁢d2⁢ηπ2⁢e(ξ−μ)⁢(η∗−ν∗)−(η−ν)⁢(ξ∗−μ∗)⁢F⁢(η,ξ)=2⁢π⁢⟨ξ=μ|F⁢(Q1,Q2,P1,P2)|ν⟩⁢ei⁢(ν1⁢μ2−ν2⁢μ1)

The inverse of Eq. (319), according to Eq. (309), is

(320)F⁢(η,ξ)=∬⁢2⁢d2⁢μ⁢d2⁢νπ⁢e(ξ∗−μ∗)⁢(η−ν)−(η∗−ν∗)⁢(ξ−μ)×⟨ξ=μ|F⁢(Q1,Q2,P1,P2)|ν⟩⁢ei⁢(ν1⁢μ2−ν2⁢μ1)

Thus through the new integration transformation, a new relationship between a phase space function F⁢(η,ξ) and its Weyl–Wigner correspondence operator F⁢(Q1,Q2,P1,P2) is revealed.

For example, from the following Weyl–Wigner correspondence

(321)4(ef+1)2⁢exp⁡[ef−1ef+1⁢(|η|2+|ξ|2)]→exp⁡{f⁢[K++K−−1]}

[K+ and K− are defined in Eqs. (193) and (196)] and Eq. (320) we have

(322)4(ef+1)2⁢∬⁢d2⁢ξ⁢d2⁢ηπ2⁢e(ξ−μ)⁢(η∗−ν∗)−(η−ν)⁢(ξ∗−μ∗)×exp⁡[ef−1ef+1⁢(|η|2+|ξ|2)]=2⁢π⁢⟨ξ=μ|F⁢(Q1,Q2,P1,P2)|ν⟩⁢ei⁢(ν1⁢μ2−ν2⁢μ1)

Using the Gaussian integration formula

(323)∬⁢d2⁢ξ⁢d2⁢ηπ2⁢e(ξ−μ)⁢(η∗−ν∗)−(η−ν)⁢(ξ∗−μ∗)⁢e−λ⁢(|η|2+|ξ|2)=11+λ2⁢exp⁡[−λ⁢(|μ|2+|ν|2)1+λ2+λ2⁢(μ⁢ν∗−μ∗⁢ν)1+λ2]

in particular, when λ=−i⁢tan⁡(π4−α2), with −λλ2+1=i2⁢tan⁡α,λ2λ2+1=12−12⁢sin⁡α, Eq. (323) becomes

(324)∬⁢d2⁢ξ⁢d2⁢ηπ2⁢e(ξ−μ)⁢(η∗−ν∗)−(η−ν)⁢(ξ∗−μ∗)×exp⁡[i⁢tan⁡(π4−α2)⁢(|η|2+|ξ|2)]=ei⁢αi⁢sin⁡αexp[i⁢(|μ|2+|ν|2)2⁢tan⁡α−μ⁢ν∗−μ∗⁢ν2⁢sin⁡α+iμ2ν1−iμ1ν2]

where exp⁡[i⁢tan⁡(π4−α2)⁢(|η|2+|ξ|2)] represents an infinite long chirplet function. By taking f=i⁢(π2−α) in Eq. (322), such that ie−i⁢α=ef, and comparing with Eq. (324) we obtain

(325)⟨ξ=μ|F⁢(Q1,Q2,P1,P2)|ν⟩=−iei⁢α2⁢π⁢sin⁡α⁢exp⁡[i⁢(|μ|2+|ν|2)2⁢tan⁡α−μ⁢ν∗−μ∗⁢ν2⁢sin⁡α]

where the right-hand side of Eq. (325) is just the CFrFT kernel whose properties can be seen in Ref. [111]. (One may compare the forms (296) and (325) to see their slight difference. For the relation between them we refer to Refs. [110, 111]). Dragoman has shown that the kernel of the CFrFT can be classically produced with rotated astigmatic optical systems that mimic the quantum entanglement. Therefore the new integration transformation (307) can convert spherical wave to CFrFT kernel. We expect this transformation could be implemented by experimentalists.

11.4 Squeezing for the generalized scaled FrFT

In some practical applications it is necessary to introduce input and output scale parameters [44, 103] into FrFT, i.e., scaled FrFT. The reason lies in two facts: (i) the scaled FrFT may be more useful and convenient for optical information processing due to the scale parameters (free parameters) introduced into FrFT; (ii) it can be reduced to the conventional FrFT under a given condition. In this subsection, by establishing the relation between the optical scaled FrFT and quantum mechanical squeezing-rotating operator transform in a one-mode case, we employ the IWOP technique and the bipartite entangled state representation of two-mode squeezing operator to extend the scaled FrFT to more general cases, such as scaled complex FrFT and entangled scaled FrFT. The properties of scaled FrFTs can be seen more clearly from the viewpoint of representation transform in quantum mechanics.

11.4.1 Quantum correspondence of the scaled FrFT

The scaled FrFT [44] of α-order is defined in such a manner that the usual FrFT is its special case, i.e.,

(326)Fα[f(x)]=ei⁢(π2−α)2⁢π⁢μ⁢ν⁢sin⁡α∫−∞∞exp[−ix2/μ2+y2/ν22⁢tan⁡α+ixx′μ⁢ν⁢sin⁡α]f(x)dx

where the exponential function is an integral kernel. In a similar way to deriving the quantum correspondence of FrFT in Eq. (248), and using the natural repression of single-mode squeezing operator S1 in coordinate representation [112–114],

(327)S1⁢(μ)=1μ⁢∫−∞∞dx⁢|xμ⟩⁢⟨x|

we have

(328)exp⁡(−i⁢x2/μ2+y2/ν22⁢tan⁡α+ixx′μ⁢ν⁢sin⁡α)=−2⁢π⁢i⁢μ⁢ν⁢ei⁢α⁢sin⁡α⁢⟨y|S1†⁢(ν)⁢exp⁡(i⁢α⁢a†⁢a)⁢S1⁢(μ)|x⟩

which implies that the integral kernel in Eq. (326) is just the matrix element of operator S1†⁢(ν)⁢exp⁡(i⁢θ⁢a†⁢a)⁢S1⁢(μ) in coordinate states. From Eq. (328) it then follows that

(329)Fα⁢[f⁢(x)]=∫−∞∞dx⁢⟨y|S1†⁢(ν)⁢ei⁢α⁢a†⁢a⁢S1⁢(μ)|x⟩⁢f⁢(x)=⟨y|S1†⁢(ν)⁢ei⁢α⁢a†⁢a⁢S1⁢(μ)|f⟩≡g⁢(y)

which suggests

(330)|g⟩=S1†⁢(ν)⁢ei⁢α⁢a†⁢a⁢S1⁢(μ)⁢|f⟩

From Eqs. (329) and (326) one can see that the scaled FrFT in Eq. (329) corresponds actually to the squeezing-rotating operator (S1†⁢(ν)⁢ei⁢θ⁢a†⁢a⁢S1⁢(μ)) transform in Eq. (329) between two quantum states.

11.4.2 The Scaled CFrFT

On the basis of quantum mechanical version of one-mode scaled FrFT, we generalize it a two-mode case, i.e., we can introduce the integral

(331)FαC[f(η)]≡⟨η|′S2†(ν)ei⁢α⁢(a1†⁢a1+a2†⁢a2)S2(μ)|f⟩=∫−∞∞d2⁢ηπ⟨η|′S2†(ν)ei⁢α⁢(a1†⁢a1+a2†⁢a2)S2(μ)|η⟩f(η)

where f⁢(η)=⟨η|f⟩. Using the natural expression of the two-mode squeezing operator S2 [Eq. (34)], and noticing that ⟨η|′ei⁢θ⁢(a1†⁢a1+a2†⁢a2)|η⟩ is just the integral kernel of CFrFT (296), we can reform (331) as

(332)FαC⁢[f⁢(η)]=ei⁢(π2−α)2⁢μ⁢ν⁢sin⁡α⁢∫d2⁢ηπ⁢f⁢(η)×exp⁡[−i(|η|2′/ν2+|η|2/μ2)2⁢tan⁡α+i(ηη′∗+η∗η)′2⁢μ⁢ν⁢sin⁡α]

It is obvious that Eq. (332) is just a generalized CFrFT with squeezing parameters. We name it the scaled CFrFT. Thus we link a two-mode squeezing-rotating operator transform to the scaled CFrFT of complex functions.

11.4.3 Entangled scaled FrFT

On the other hand, recall that the entangled state |η⟩ can be Schmidt-decomposed as [33]

(333)|η⟩=e−i⁢η1⁢η2⁢∫−∞∞dx⁢|x⟩1⊗|x−2⁢η1⟩2⁢ei⁢2⁢x⁢η2

We see that

(334)⟨x1',x2'|η⟩′=e−iη1'η2'δ(2η1'+x2'−x1')ei2x1'η2'⟨η|x1,x2⟩=ei⁢η1⁢η2⁢δ⁢(2⁢η1+x2−x1)⁢e−i⁢2⁢x1⁢η2

Using Eq. (30), we have

(335)K(x1',x2',x1,x2)≡⟨x1',x2'|S2†(ν)ei⁢α⁢(a1†⁢a1+a2†⁢a2)S2(μ)|x1,x2⟩=∫d2ηd2η′π2⟨x1′,x2′|η⟩′⟨η|′S2†(ν)ei⁢α⁢(a1†⁢a1+a2†⁢a2)×S2⁢(μ)⁢|η⟩⁢⟨η|x1,x2⟩

where |x1,x2⟩=|x1⟩⊗|x2⟩. Substituting Eqs. (333) and (334) into Eq. (335), we can derive

(336)K(x1',x2',x1,x2)=ei⁢(π2−α)2⁢π⁢sin⁡α⁢exp⁡(−i⁢λν2′+λμ22⁢tan⁡α+iλμλν′sin⁡α)×ei⁢(π2−α)2⁢π⁢sin⁡α⁢exp⁡(−i⁢κν2′+κμ22⁢tan⁡α+iκμκν′sin⁡α)

where λμ=x1−x22⁢μ,λν=′x1'−x2'2⁢ν;κμ=μ⁢(x2+x1)2,κν=′ν(x1'+x2')2.

From Eq. (336) one can see that a new 2D scaled FrFT can be composed of one 1D scaled FFT in its space domain and the other in its “frequency” domain, with the transform variables being the combination of two coordinates. Hence Eq. (336) is quite different from the direct product of the two 1D scaled FrFTs that are both in “space domain” indicated in Eq. (329). Note that the new 2D scaled FFT is still characterized by only 3-parameter. Therefore, for any function f⁢(x1,x2)=⟨x1,x2|f⟩ we can define an entangled scaled FrFT, i.e.,

(337)FαE[f(x1,x2)]=∫−∞∞K(x1',x2',x1,x2)f(x1,x2)dx1dx2=⟨x1′,x2'|S2†(ν)ei⁢α⁢(a1†⁢a1+a2†⁢a2)S2(μ)|f⟩

Next we examine the properties of these scaled FrFTs in the quantum optics context. Without losing generality, for the additivity property, we consider the scaled CFrFT,

(338)Fα+βC[f(η)]≡∫d2⁢ηπ⟨η|′S2†(ν)ei⁢(α+β)⁢(a1†⁢a1+a2†⁢a2)×S2⁢(μ)⁢|η⟩⁢f⁢(η)

Inserting the completeness relation of |η⟩ into Eq. (338), we get

(339)Fα+βC[f(η)]=∫d2η′′π⟨η|′S2†(ν)ei⁢α⁢(a1†⁢a1+a2†⁢a2)×S2(τ)|η⟩′′=τFβ[f(η)]=FαFβ[f(η)]

which is just the additivity property. It should be pointed out that the condition of additive operator for the scaled FrFTs is that the parameter ν′ of the prior cascade should be equal to the parameter μ′ of the next one, i.e., μ=′ν′. For other scaled FrFTs, the properties can also be discussed in the similar way (according to their quantum versions).

To this end, we should emphasize that different scaled FrFTs correspond to different quantum mechanical squeezing operators or representations. That is to say, it is possible that some other scaled FrFT can be presented by using different quantum mechanical squeezing operators or representations.

12 Adaption of Collins diffraction formula and CFrFT

The connection between the Fresnel diffraction in free space and the FrFT had been bridged by Pellat–Finet [115] who found that FrFTs are adapted to the mathematical expression of Fresnel diffraction, just as the standard Fourier transform is adapted to Fraunhofer diffraction. In previous sections, a new formulation of the CFrFT and the Collins diffraction formula are respectively derived in the context of representation transform of quantum optics. In this section we inquire if the adaption problem of Collins diffraction formula to the CFrFT can also be tackled in the context of quantum optics. We shall treat this topic with the use of two-mode (3 parameters) squeezing operator and in the entangled state representation of continuous variables. In so doing the quantum mechanical version of the associated theory of classical diffraction and classical CFrFT is obtained, which connects classical optics and quantum optics in this aspect.

For Gaussian beam, the A⁢B⁢C⁢D rule is equally derived via optical diffraction integral theory — the Collins integral formula. As shown in Fig. 2, if f⁢(η) represents the input field amplitude at point η on S1, and g(η)′ denotes the diffraction field amplitude at point η′ on S2. Then Collins formula in complex form takes the form (181) and (182). Next we shall examine the adaption of the Collins formula to the CFrFT by virtue of the entangled state representation in quantum optics [108].

12.1 Adaption of the Collins formula to CFrFT

Using the completeness relation of |η⟩, we can further put Eq. (181) into

(340)g(η)′=⟨η|′U2(r,s)|f⟩=⟨η′|U2(r,s)μ12∫d2⁢ηπ|μ1η⟩⟨μ1η|f⟩

Take η′=BDσK,μ1=BA/L and write out

(341)g(η)′→⟨BDσK|g⟩≡G(σ),f(μ1η)≡F(η)

where K and L are two constants to be determined later. Then according to Eqs. (181) and (340) we have

(342)G(σ)=μ12∫d2⁢ηπ⟨η|′U2(r,s)|μ1η⟩F(η)=12⁢i⁢A⁢L2⁢exp⁡[i⁢|σ|22⁢K2]⁢∫d2⁢ηπ×exp⁡[i⁢|η|22⁢L2−i⁢(σ∗⁢η+σ⁢η∗)2⁢L⁢K⁢A⁢D]⁢F⁢(η)

Comparing Eq. (342) with Eq. (182) leads us to choose

(343)L2=tan⁡α,K=sin⁡2⁢α/(2⁢A⁢D)

Then Eq. (342) becomes

(344)G⁢(σ)=cos⁡αi2⁢A⁢sin⁡α⁢exp⁡(i⁢A⁢D−cos2⁡αsin⁡2⁢α⁢|σ|2)×∫d2⁢ηπexp[i⁢(|η|2+|σ|2)2⁢tan⁡α−i⁢(σ∗⁢η+σ⁢η∗)2⁢sin⁡α]F(η)=cos⁡αA⁢e−i⁢α⁢exp⁡(i⁢A⁢D−cos2⁡αsin⁡2⁢α⁢|σ|2)⁢Fα⁢[F]⁢(σ)

so Eq. (344) is a standard CFrFT up to a quadratic phase term exp⁡(i⁢A⁢D−cos2⁡αsin⁡2⁢α⁢|σ|2), according to Eq. (343) and BDσK=η′, it can also be written as

(345)exp(iA⁢D−cos2⁡αsin⁡2⁢α|σ|2)=exp(iR|η|2′)

which represents a quadratic approximation to a sphere wave diverging from a luminous point at the distance

(346)R=2⁢A⁢BA⁢D−cos2⁡α

from S2. Let S be the sphere tangent to S2 with radius R (see Fig. 2). A point on S is located by its projection on S2, which means that the coordinates on S2 can also be used as those on S. Therefore, the quadratic phase term can be compensated if the output field is observed on S instead of S2. Then, after considering the phase compensation, the field transforms from S1 to S is

(347)GS⁢(σ)=cos⁡αA⁢e−i⁢α⁢Fα⁢[F]⁢(σ)

In this way, the field amplitude on S is the perfect α-th CFrFT of the field amplitude on S1.

12.2 Adaption of the additivity property of CFrFT to the Collins formula for two successive Fresnel diffractions

The most important property of FrFT is that Fα obeys the additivity rule, i.e., two successive FrFT of order α and β makes up the FrFT of order α+β. For the CFrFT, its additivity property is proven in Eq. (300). For Collins diffraction from S1 to S′ (see Fig. 1), the additivity means that the diffraction pattern observed on S′(η¯)  (the sphere tangent to S3 with radius R′) and associated with Fα+β should be the result of a first diffraction phenomenon (associated with Fα) on S (with η′), followed by a second diffraction phenomenon (associated with Fβ) from S to S′. This is a necessary consequence of the Huygens principle. Next we prove that such is indeed the case.

First, let us consider the field transform from S1 (with η) to S′ (see Fig. 2) described by the ray transfer matrix [A,′B,′C,′D′]. Similar to deriving Eq. (346), after the squeezing transform and the phase compensation,

(348)R′=2AB′′AD′−′cos2α′→exp(iR′|η¯|2)

Thus we can obtain the expression of CFrFT for Collins diffraction from S1 to S′ (not S3),

(349)GS′(σ)′=cosα′A′e−iα′Fα′[f(μ1'η)](σ)′

where η¯=B′D′⁢σ′K′ and μ1'=B′A′/L,′

(350)L′2=tan α′,K′=sin 2α′/(2A′D′)

Eq. (349) is the same in form as Eq. (347) but with primed variables. Using Eqs. (348) and (181) one can prove that the transform from S1 to S′ is [see Eqs. (181), (344)–(346)]

(351)gS′(η¯)=exp(−iR′|η¯|2)⟨η¯|U2(r,′s)′|f⟩≡GS′(σ)′

In Eq. (351) we have taken the phase compensation term (348) into account.

Second, let us consider the second diffraction from S to S3 determined by the ray transfer matrix [A,′′B,′′C,′′D′′]. For this purpose, using the group multiplication rule of F2⁢(r,s), we can decompose the diffraction from S1 to S′ into two parts: one is described as the matrix [A,B,C,D] from the plane S1 (with η) to S2⁢(S) (with η′), the other is [A,′′B,′′C,′′D]′′ from the plane S2 to S3(S)′ (with η¯), then the total matrix from S1 to S3 is

(352)(A′B′C′D′)=(A′′B′′C′′D′′)⁢(ABCD)

Using Eq. (30) and the group multiplication rule of F2(r,′s)′, we can further put Eq. (351) into another form

(353)GS′(σ)′=exp(−iR′|η¯|2)⟨η¯|U2(r,′′s)′′U2(r,s)|f⟩=BD⁢K2⁢∫d2⁢σπ⁢exp⁡(−i⁢|η¯|2R′)×⟨η¯|U2(r,′′s)′′|μ2'σ⟩G(σ)

where μ2'=BD1K and we have made a reasonable assumption that μ1'=μ1 [hence f⁢(μ1'η)=F⁢(η)], which means that there are same scaled variants for the input field amplitudes on S1 of the diffractions from S1 to S and from S1 to S′.

In order to examine the second diffraction domain from S to S′ (not S3), we need to translate the output field amplitude G⁢(σ) observed on the plane S2 to the field amplitude observed on sphere plane S, i.e., putting G⁢(σ) into GS⁢(σ) [see Eq. (347)] by taking the phase compensation [see Eq. (342)] into account. Thus the field transform from S to S′ is [GS(σ)→GS′(σ)′]

(354)GS′(σ)′=BD⁢K2∫d2⁢σπexp(iRB⁢|σ|2D⁢K2−i⁢|η¯|2R′)×⟨η¯|U2(r,′′s)′′|μ2'σ⟩GS(σ)=BD⁢K212iB′′∫d2⁢σπexp[iB′K2′D′(D′′2B′′−1R′)|σ|2′+i⁢B⁢|σ|2D⁢K2(A′′2B′′+1R)−i⁢(σ⁢σ′⁣∗+σ⁢σ∗′)2⁢B⁢K′′⁢K′⁢DD′BB′]GS(σ)

Comparing Eq. (354) with Eq. (182), we have

(355)sin⁡β=B⁢K′′⁢K′⁢DD′BB′=B′′2⁢sin⁡2⁢α⁢sin′⁡2⁢αAA′BB′

and noticing that μ1'=μ1 yields B′Atan′α′=BA⁢tan⁡α, we have

(356)A′=B′′Bcos⁡α⁢sin′⁡αsin⁡β,A=B′′B′sin⁡α′⁢cos⁡αsin⁡β

combine Eqs. (340) and (356), and it then follows that (let α′=α+β)

(357)BD⁢K2⁢12iB′′=A⁢Bsin⁡2⁢α⁢1iB′′=Ai⁢A′⁢cos⁡α⁢cosα′2⁢sin⁡β

and

(358)B′K2′D′⁢(D′′2B′′−1R′)=12⁢cot⁡βBD⁢K2⁢(A′′2B′′+1R)=12⁢cot⁡β

where we have used Eqs. (346), (348) and (352).

Substitution of Eqs. (347), (355), (357) and (358) into Eq. (354) yields

(359)GS′(σ)′=Ai⁢A′⁢cos⁡αcosα′2⁢sin⁡β∫d2⁢σπexp[i(|σ|2′+|σ|2)2⁢tan⁡β−i(σσ∗′+σσ∗′)2⁢sin⁡β]GS(σ)=cos⁡α⁢e−i⁢α′iA′⁢12⁢sin⁡β⁢∫d2⁢σπ⁢exp⁡[i(|σ|2′+|σ|2)2⁢tan⁡β−i⁢(σ⁢σ′⁣∗+σ⁢σ∗′)2⁢sin⁡β]⁢Fα⁢[F]⁢(σ)=cosα′A′e−i⁢(α+β)FβFα[F](σ)′

The first equation of Eq. (359) indicates that it is just a CFrFT of GS⁢(σ) from S to S′. Comparing Eq. (359) with Eq. (349), we see

(360)FβFα[F](σ)′=Fα+β[F](σ)′

So far we have completed the study of adaption of CFrFT to the mathematical representation of Collins diffraction formula in quantum optics context.

13 Wavelet transformation and the IWOP technique

In recent years wavelet transforms [116, 117] have been developed which can overcome some shortcomings of the classical Fourier analysis and therefore has been widely used in Fourier optics and information science since 1980s. Here we present a quantum optical version of classical wavelet transform (WT) by virtue of the IWOP technique.

13.1 Quantum optical version of classical WTs

A wavelet has its energy concentrated in time to give a tool for the analysis of transient, nonstationary, or time-varying phenomena. (It is a wavelet because it is localized and it resembles a wave because it oscillates.) Mathematically, wavelets are defined by starting with a function ψ of the real variable x, named a mother wavelet which is required to decrease rapidly to zero as |x| tends to infinity,

(361)∫−∞∞ψ⁢(x)⁢dx=0

A more general requirement for a mother wavelet is for demanded ψ⁢(x) to have vanishing moments ∫−∞∞xl⁢ψ⁢(x)⁢dx=0,l=0,1,2⁢…,L. (A greater degree of smoothness than continuity also leads to vanishing moments for the mother wavelet). The theory of wavelets is concerned with the representation of a function in terms of a two-parameter family of dilates and translates of a fixed function. The mother wavelet ψ generates the other wavelets of the family ψ(μ,s) (μ is scaling parameter, s is a translation parameter, s∈R), and the dilated-translated function is defined as

(362)ψ(μ,s)⁢(x)=1|μ|⁢ψ⁢(x−sμ)

while the wavelet integral transform of a signal function f⁢(x)∈L2⁢(R) by ψ is defined by

(363)Wψ⁢f⁢(μ,s)=1|μ|⁢∫−∞∞f⁢(x)⁢ψ∗⁢(x−sμ)⁢dx

We can express Eq. (363) as

(364)Wψ⁢f⁢(μ,s)=⟨ψ|U⁢(μ,s)|f⟩

where ⟨ψ| is the state vector corresponding to the given mother wavelet, |f⟩ is the state to be transformed, and

(365)U⁢(μ,s)≡1|μ|⁢∫−∞∞|x−sμ⟩⁢⟨x|⁢dx

is the squeezing-translating operator [116, 118, 119], ⟨x| is the eigenvector of coordinate operator. In order to combine the wavelet transform with quantum states transform more tightly, using the IWOP technique, we can directly perform the integral in Eq. (365) (Q=(a+a†)/2,μ>0)

(366)U⁢(μ,s)=exp⁡[−s22⁢(1+μ2)−a†22⁢tanh⁡λ−a†⁢s2⁢sech⁢λ]×exp⁡[(a†⁢a+12)⁢ln⁡sech⁢λ]×exp⁡(a22⁢tanh⁡λ+s⁢a2⁢sech⁢λ)

where μ=eλ, sech⁢λ=2⁢μ1+μ2, tanh⁡λ=μ2−1μ2+1, and the operator identity eg⁢a†⁢a=:exp[(eg−1)a†a]: is used. In particular, when s=0, it is reduced to the well-known squeezing operator,

(367)U⁢(μ,0)=1μ⁢∫−∞∞|xμ⟩⁢⟨x|⁢dx=exp⁡[λ2⁢(a2−a†2)]

For a review of the squeezed state theory we refer to Ref. [14].

13.2 The condition of mother wavelet in the context of quantum optics

Now we analyze the condition (361) for mother wavelet from the point of view of quantum optics. Due to

(368)∫−∞∞|x⟩⁢dx=|p=0⟩

where |p⟩ is the momentum eigenstate, we can recast the condition into quantum mechanics as

(369)∫−∞∞ψ⁢(x)⁢dx=0→⟨p=0|ψ⟩=0

which indicates that the probability of a measurement of |ψ⟩ by the projection operator |p⟩⁢⟨p| with value p=0 being zero. Withouta loss of generality, we suppose

(370)|ψ⟩M=G⁢(a†)⁢|0⟩=∑n=0∞gn⁢a†n⁢|0⟩

where gn are chosen as to let |ψ⟩ obey the condition (361). Using the coherent states’ overcompleteness relation, we have

(371)⟨p=0|ψ⟩=⟨p=0|∫d2⁢zπ|z⟩⁢⟨z|∑ngn⁢a†n|0⟩=∑ng2⁢n=0

Eq. (371) provides a general formalism in finding the qualified wavelets. For example, assuming g2⁢n=0 for n>3, the coefficients of the survived terms should satisfy

(372)g0+g2+3⁢g4+15⁢g6=0

and |ψ⟩becomes

(373)|ψ⟩=(g0+g2⁢a†2+g4⁢a†4+g6⁢a†6)⁢|0⟩

Projecting it onto the coordinate representation, we get the qualified wavelets

(374)ψ(x)=π−1/4e−x2/2[g0+g2(2x2−1)+g4⁢(4⁢x4−12⁢x2+3)+g6(8x6−60x4+90x2−15)]

where we have used ⟨x|n⟩=(2n⁢n!⁢π)−1/2⁢Hn⁢(x)⁢e−x2/2, and Hn⁢(x) is the Hermite polynomials. Now we take some examples.

Case 1: In Eq. (371), by taking g0=12, g2=−12, g2⁢n=0 (otherwise), we have

(375)|ψ⟩M=12⁢(1−a†2)⁢|0⟩

It then follows that

(376)ψM⁢(x)≡12⁢⟨x|(1−a†2)|0⟩=π−1/4⁢e−x2/2⁢(1−x2)

which is just the Maxican hat wavelet, satisfying the condition ∫−∞∞e−x2/2⁢(1−x2)⁢dx=0. Hence 12⁢(1−a†2)⁢|0⟩ is the state vector corresponding to the Maxican hat mother wavelet (see Fig. 3). Once the state vector ⟨ψ| corresponding to mother wavelet is known, for any state |f⟩ the matrix element ⟨ψ|U(μ,s)|f⟩M is just the wavelet transform of f⁢(x) with respect to ⟨ψ|.

Case 2: When g0=−2, g2=−1, g4=1 and g6=0, from Eq. (374) we obtain (see Fig. 4)

(377)ψ2⁢(x)=2⁢π−1/4⁢e−x2/2⁢(2⁢x4−7⁢x2+1)

which also satisfies ∫−∞∞ψ2⁢(x)⁢dx=0. Note that when g0=−1, g2=−2, g4=1 and g6=0, we obtain a slightly different wavelet (see Fig. 5). Therefore, as long as the parameters g2⁢n conforms to condition (372), we can adjust their values to control the shape of the wavelets.

Case 3: When g0=1, g2=2, g4=4 and g6=−1, we get (see Fig. 6)

(378)ψ3⁢(x)=π−1/4⁢e−x2/2⁢(−8⁢x6+76⁢x4−134⁢x2+26)

and ∫−∞∞ψ3⁢(x)⁢dx=0. From these figures we observe that the number of the nodes of the curves at the x-axis is equal to the highest power of the wavelet functions.

To further reveal the properties of the newly found wavelets, we compare the wavelet transform computed with the well-known Mexican hat wavelet ψ1⁢(x) and that with our new wavelet ψ2⁢(x). Concretely, we map a simple cosine signal cos⁡(π⁢x) by performing the wavelet transforms with ψi⁢[T⁢(x−X)],i=1,2, into a 2D space (X,T), where X denotes the location of a wavelet and a its size. The resulting wavelet transforms by ψ1⁢(x) (=ψM⁢(x)) and ψ2⁢(x) are

(379)Ω1⁢(X,T)=23⁢∫−∞∞dx⁢ψ1⁢[T⁢(x−X)]⁢cos⁡(π⁢x)

(380)Ω2⁢(X,T)=130⁢∫−∞∞dx⁢ψ2⁢[T⁢(x−X)]⁢cos⁡(π⁢x)

where 2/3 and 1/30 are the normalization factors for ψ1 and ψ2 respectively. The wavelet integral Ωi⁢(X,T) are also called wavelet coefficient which measures the variation of cosπ⁢x in a neighborhood of X, whose size is proportional to 1/T. The contour line representation of Ω1⁢(X,T) and Ω2⁢(X,T) are depicted in Fig. 7 and Fig. 8, respectively, where the transverse axis is X-axis (time axis), while the longitudinal axis (T-axis) is the frequency axis.

It is remarkable that although two overall shapes of the two contour lines look similar, there exist two notable differences between these two figures: (i) Along T-axis Ω1⁢(X,T) has one maximum, while Ω2⁢(X,T) has one main maximum and one subsidiary maximum (“two islands”), so when ψ2 scales its size people have one more chance to identify the frequency information of the cosine wave than using ψM. Interestingly enough, the “two islands” of Ω2⁢(X,T) in Fig. 8 can be imagined as if they were produced while the figure of Ω1⁢(X,T) deforms into two sub-structures along a-axis. (ii) Near the maximum of Ω2⁢(X,T) the density of the contour lines along a-axis is higher than that of Ω1⁢(X,T), which indicates that the new wavelet ψ2 is more sensitive in detecting frequency information of the signal at this point. Therefore, ψ2⁢(x) may be superior to ψM⁢(x) in analyzing some signals. Finally, we mention that there exist some remarkable qualitative similarities between the mother wavelets presented in Figs. 3 through 6 and some of the amplitude envelopes of higher order laser spatial modes and spatial supermodes of phase locked diode laser arrays [120–122], which are due to spatial coherence.

13.3 Quantum mechanical version of Parseval theorem for WT

In this subsection, we shall prove that the Parseval theorem of 1D WT [118, 119, 123]:

(381)∫−∞∞d⁢μμ2⁢∫−∞∞ds⁢Wψ⁢f1⁢(μ,s)⁢Wψ∗⁢f2⁢(μ,s)=2⁢Cψ⁢∫f1⁢(x)⁢f2∗⁢(x)⁢dx

where ψ⁢(x) is a mother wavelet whose Fourier transform is ψ⁢(p),Cψ=2⁢π⁢∫0∞|ψ⁢(p)|2p⁢dp<∞. In the context of quantum mechanics, according to Eq. (364), we see that the quantum mechanical version of Parseval theorem should be

(382)∫−∞∞d⁢μμ2⁢∫−∞∞ds⁢⟨ψ|U⁢(μ,s)|f1⟩⁢⟨f2|U†⁢(μ,s)|ψ⟩=2⁢Cψ⁢⟨f2|f1⟩

and since ψ⁢(x)=⟨x|ψ⟩,ψ⁢(p) involved in Cψ is ⟨p|ψ⟩, and ⟨p| is the momentum eigenvector

(383)ψ⁢(p)=⟨p|ψ⟩=12⁢π⁢∫−∞∞dx⁢ψ⁢(x)⁢e−i⁢p⁢x

Eq. (382) indicates that once the state vector ⟨ψ| corresponding to mother wavelet is known, for any two states |f1⟩ and |f2⟩, their overlap up to the factor Cψ [determined by Eq. (396)] is just their corresponding overlap of WTs in the (μ,s) parametric space.

Proof of Equation (382): In order to show Eq. (382), we calculate

(384)U†⁢(μ,s)⁢|p⟩=1|μ|⁢e−i⁢p⁢sμ⁢|pμ⟩

which leads to

(385)∫−∞∞dsU†(μ,s)|p⟩′⟨p|U(μ,s)=2πδ(p−p)′|p′μ⟩⟨pμ|

where we have used the formula

(386)∫−∞∞d⁢x2⁢πeix(p−p)′=δ(p−p)′

Inserting the completeness relation ∫−∞∞dp⁢|p⟩⁢⟨p|=1 into the left side of Eq. (381) and then using Eq. (385), we have

(387)L.H.S⁢of⁢Eq.(381)=2⁢π⁢∫−∞∞d⁢μμ2⁢∫−∞∞dp⁢ψ∗⁢(p)⁢ψ⁢(p)⁢⟨f2|pμ⟩⁢⟨pμ|f1⟩≡I1+I2

where

(388)I1=2⁢π⁢∫−∞∞dp⁢[∫0∞|ψ⁢(μ⁢p)|2⁢d⁢μμ]⁢⟨f2|p⟩⁢⟨p|f1⟩

and

(389)I2=2⁢π⁢∫−∞∞dp⁢[∫0∞|ψ⁢(−μ⁢p)|2⁢d⁢μμ]⁢⟨f2|p⟩⁢⟨p|f1⟩

Further, we can put Eqs. (388) and (389) into the following forms:

(390)I1=Cψ⁢∫0∞dp⁢⟨f2|p⟩⁢⟨p|f1⟩+Cψ′⁢∫−∞0dp⁢⟨f2|p⟩⁢⟨p|f1⟩

and

(391)I2=Cψ′⁢∫0∞dp⁢⟨f2|p⟩⁢⟨p|f1⟩+Cψ⁢∫−∞0dp⁢⟨f2|p⟩⁢⟨p|f1⟩

where

(392)Cψ=2⁢π⁢∫0∞|ψ⁢(μ⁢p)|2⁢d⁢μμ=2⁢π⁢∫0∞|ψ⁢(p)|2⁢d⁢ppCψ′=2π∫0−∞|ψ(p)′|2dp′p′=2π∫0∞|ψ(−p)|2d⁢pp

Thus, when the definite integration satisfies the admissible condition, i.e.,

(393)∫0∞|ψ⁢(p)|2⁢d⁢pp=∫0∞|ψ⁢(−p)|2⁢d⁢pp

which leads to

(394)2⁢π⁢∫−∞∞|ψ⁢(p)|2⁢d⁢p|p|=2⁢Cψ

Eq. (387) can be transformed to

(395)L.H.S.of⁢Eq.(381)=2⁢Cψ⁢∫−∞∞dp⁢⟨f2|p⟩⁢⟨p|f1⟩=R.H.S⁢of⁢Eq.(381)

where

(396)Cψ≡2⁢π⁢∫0∞|ψ⁢(p)|2p⁢dp<∞

Thus, the theorem is proved. Especially, when |f1⟩=|f2⟩, Eq. (382) becomes

(397)∫−∞∞d⁢μμ2⁢∫−∞∞ds⁢|⟨ψ|U⁢(μ,s)|f1⟩|2=2⁢Cψ⁢⟨f1|f1⟩

which is named isometry of energy.

13.4 Inversion formula of WT

Now we can directly derive the inversion formula of WT. In other words, we take ⟨f2|=⟨x| in Eq. (382), then using Eq. (364), we see that Eq. (382) is reduced to

(398)∫−∞∞d⁢μμ2⁢∫−∞∞ds⁢Wψ⁢f1⁢(μ,s)⁢⟨x|U†⁢(μ,s)|ψ⟩=2⁢Cψ⁢⟨x|f1⟩

Due to Eq. (365) we have

(399)〈x|U†(μ,s)=1|μ|〈x|∫−∞∞dx′|x′〉〈x′−sμ|=1|μ|〈x−sμ|

It then follows that

(400)∫−∞∞d⁢μμ2⁢∫−∞∞ds⁢Wψ⁢f1⁢(μ,s)⁢1|μ|⁢⟨x−sμ|ψ⟩=2⁢Cψ⁢⟨x|f1⟩

which means

(401)f1⁢(x)=12⁢Cψ⁢∫−∞∞d⁢μμ2⁢|μ|⁢∫−∞∞ds⁢ψ⁢(x−sμ)⁢Wψ⁢f1⁢(μ,s)

This is the inversion formula of WT.

13.5 New orthogonal property of mother wavelet in parameter space

Form the Parserval theorem (381) of WT in quantum mechanics we can derive some new properties of the mother wavelet [124]. Taking |f1⟩=|x⟩, |f2⟩=|x⟩′ in Eq. (381) one can see that

(402)∫−∞∞d⁢μμ2⁢|μ|⁢∫−∞∞ds⁢ψ⁢(x−sμ)⁢ψ∗⁢(x−′sμ)=2Cψδ(x−x)′

which is a new orthogonal property of mother wavelet in parameter space spanned by (μ,s). In a similar way, we take |f1⟩=|f2⟩=|n⟩, a number state, since ⟨n|n⟩=1, then we have

(403)∫−∞∞d⁢μμ2⁢∫−∞∞ds⁢|⟨ψ|U⁢(μ,s)|n⟩|2=2⁢Cψ

Or take |f1⟩=|f2⟩=|z⟩, and |z⟩=exp⁡(−|z|2/2+z⁢a†)⁢|0⟩ is the coherent state, then

(404)∫−∞∞d⁢μμ2⁢∫−∞∞ds⁢|⟨ψ|U⁢(μ,s)|z⟩|2=2⁢Cψ

This indicates that Cψ is |f1⟩-independent, which coincides with the expression in Eq. (396). Next, we consider a special example. When the mother wavelet is the Mexican hat (376), we have

(405)ψM⁢(p)≡⟨p|ψM⟩=12⁢(⟨p|0⟩−2⁢⟨p|2⟩)=π−1/4⁢p2⁢e−12⁢p2

where

(406)⟨p|n⟩=(−i)n2n⁢n!⁢π⁢e−p2/2⁢Hn⁢(p)

Here Hn⁢(p) is the single-variable Hermit polynomial. Substituting Eq. (405) into Eq. (396), we have

(407)Cψ≡2⁢π⁢∫0∞|ψM⁢(p)|2p⁢dp=π

Thus, for the Mexican hat wavelet (376), we see

(408)∫−∞∞d⁢μμ2⁢|μ|⁢∫−∞∞ds⁢ψM⁢(x−sμ)⁢ψM∗⁢(x−′sμ)=2πδ(x−x)′

Eq. (408) can be checked as follows. Using Eq. (376) and noticing that ψM⁢(x)=ψM⁢(−x), we can put the left-hand side of Eq. (408) into

(409)L.H.S.of⁢Eq.(408)=2∫0∞du∫−∞∞dsψM(ux−s)ψM∗(ux−′s)={0,x≠x′32⁢∫0∞du→∞,x=x′=R.H.S.of⁢Eq.(408)

where we have used the integration formulas

(410)∫−∞∞(1−s2)2⁢exp⁡(−s2)⁢ds=34⁢π

and

(411)∫−∞∞(1−s2)⁢[1−(s−b)2]⁢e−s2/2−(b−s)2/2⁢ds=π16⁢e−b24⁢[12+b2⁢(b2−12)]

Next, we examine if the Morlet wavelet obey the formalism (402). The Morlet wavelet is defined as [125–127]

(412)ψm⁢o⁢r⁢(x)=π−1/4⁢(ei⁢f⁢x−e−f2/2)⁢e−x2/2

Substituting Eq. (412) into the left hand side of (402), we get

(413)I≡∫−∞∞d⁢μμ2⁢|μ|⁢∫−∞∞ds⁢ψmor⁢(x−sμ)⁢ψmor∗⁢(x−′sμ)={0,x≠x′2⁢(1+e−f2−2⁢e−3⁢f2/4)⁢∫0∞d⁢μμ2→∞,x=x′

Thus the Morlet wavelet satisfies Eq. (402).

13.6 WT and Wigner–Husimi distribution function

Phase space technique has been proved very useful in various branches of physics. Distribution functions in phase space have been a major topic in studying quantum mechanics and quantum statistics. Among various phase space distributions the Wigner function Fw⁢(q,p) [90–94] is the most popular, since its two marginal distributions lead to measuring probability density in coordinate space and momentum space, respectively. But the Wigner distribution function itself is not a probability distribution because it is both positive and negative. In spite of its attractive formal properties, it needs improving. To overcome this inconvenience, the Husimi distribution function Fh(q,′p)′ is introduced [128], which is defined in a manner that guarantees its nonnegative ness. Its definition is smoothing out the Wigner function by averaging over a “coarse graining” function,

(414)Fh(q,p,κ)=∫∫−∞∞dqd′pFw′(q,′p)′×exp[−κ(q−′q)2−(p−′p)2κ]

κ>0 is the Gaussian spatial width parameter, which is free to be chosen and determines the relative resolution in p-space versus q-space.

In the following, we shall employ the optical wavelet transformation to study the Husimi distribution function. That is, we shall show that the Husimi distribution function of a quantum state |ψ⟩ can be obtained by making a WT of the Gaussian function e−x2/2, i.e.,

(415)⟨ψ|Δh⁢(q,p,κ)|ψ⟩=e−p2κπ⁢κ⁢|∫−∞∞dx⁢ψ∗⁢(x−sμ)⁢e−x2/2|2

where

(416)s=−1κ⁢(κ⁢q+i⁢p),μ=κ

and ⟨ψ|Δh⁢(q,p)|ψ⟩ is the Husimi distribution function as well as Δh⁢(q,p,κ) is the Husimi operator,

(417)Δh⁢(q,p,κ)=2⁢κ1+κ:exp⁡[−κ⁢(q−Q)21+κ−(p−P)21+κ]:

Here :: denotes normal ordering; Q and P are the coordinate and the momentum operator.

Proof of Eq. (415): According to Eqs. (363) and (364), when |f⟩ is the vacuum state |0⟩,e−x2/2=π1/4⁢⟨x|0⟩, we see that

(418)π−1/4⁢∫−∞∞d⁢xμ⁢ψ∗⁢(x−sμ)⁢e−x2/2⁢dx=⟨ψ|U⁢(μ,s)|0⟩

From Eq. (366) it then follows that

(419)U(μ,s)|0⟩=sech1/2λexp[−s22⁢(1+μ2)−a†⁢s2sechλ−a†22tanhλ]|0⟩

Substituting Eq. (416) and tanh⁡λ=κ−1κ+1, cosh⁡λ=1+κ2⁢κ into Eq. (419), we get

(420)e−p22⁢κ+i⁢p⁢qκ+1⁢U⁢(μ=κ,s=−κ⁢q−i⁢p/κ)⁢|0⟩=(2⁢κ1+κ)1/2exp[−κ⁢q22⁢(1+κ)−p22⁢(1+κ)+2⁢a†1+κ(κq+ip)+1−κ2⁢(1+κ)a†2]|0⟩≡|p,q⟩κ

Then the WT of Eq. (418) can be further expressed as

(421)e−p22⁢κ+i⁢p⁢qκ+1⁢∫−∞∞d⁢x(κ⁢π)1/4⁢ψ∗⁢(x−sμ)⁢e−x2/2=⟨ψ|p,q⟩κ

Using the normally ordered form of the vacuum state projector |0〉〈0|=:e−a†⁢a:, and the IWOP method as well as Eq. (420), we have

(422)|p,q⟩κ⁢κ⁢⟨p,q|=2⁢κ1+κ:exp⁡[−κ⁢(q−Q)21+κ−(p−P)21+κ]:=Δh⁢(q,p,κ)

Now we explain why Δh⁢(q,p,κ) is the Husimi operator. Use the formula for converting an operator A into its Weyl ordering form [129]

(423)A=2∫d2βπ〈−β|A|β〉::exp[2(β∗a−a†β+a†a)]::d2β=dβ1dβ2,β=β1+iβ2

where the symbol ::⁢:: denotes the Weyl ordering, and |β⟩ is the usual coherent state. Substitute Eq. (422) into Eq. (423) and perform the integration by virtue of the technique of integration within a Weyl ordered product of operators. Then we obtain

(424)|p,q⟩κ⁢κ⁢⟨p,q|=2⁢::⁢exp⁡[−κ⁢(q−Q)2−(p−P)2κ]⁢::

This is the Weyl ordering form of |p,q⟩κ⁢κ⁢⟨p,q|. Then according to Weyl quantization scheme [35] we know the classical corresponding function of a Weyl ordered operator is obtained by just replacing Q→q,′P→p,′

(425)::⁢exp⁡[−κ⁢(q−Q)2−(p−P)2κ]⁢::→exp[−κ(q−q)′2−(p−p)′2κ]

and in this case the Weyl rule is expressed as

(426)|p,q⟩κ⁢κ⟨p,q|=2∫dqd′pΔw′(q,′p)′×exp[−κ(q−′q)2−(p−′p)2κ]

where at the last step we used the Weyl ordering form of the Wigner operator Δw⁢(q,p) [130]

(427)Δw⁢(q,p)=::⁢δ⁢(q−Q)⁢δ⁢(p−P)⁢::

In reference to Eq. (414) in which the relation between the Husimi function and the WF is shown, we know that the right-hand side of Eq. (426) should be just the Husimi operator, i.e.,

(428)|p,q⟩κ⁢κ⁢⟨p,q|=Δh⁢(q,p,κ)

Thus Eq. (415) is proved by combining Eqs. (428) and (421).

Thus the optical WT can be used to study the Husimi distribution function in the quantum optics phase space theory [131].

14 Symplectic wavelet transformation (SWT)

In this section we shall generalize the usual wavelet transform to symplectic wavelet transformation (SWT) by using the coherent state representation [132].

14.1 Single-mode SWT

First we are motivated to generalize the usual wavelet transform, which concerns dilation, to optical Fresnel transform (we will explain this in detail in section below). In other words, we shall use the symplectic-transformed — translated versions of the mother wavelet

(429)ψr,s;κ⁢(z)=s∗⁢ψ⁢[s⁢(z−κ)−r⁢(z∗−κ∗)]

as a weighting function to synthesize the original complex signal f⁢(z),

(430)Wψf(r,s;κ)=∫d2zπf(z)ψr,s;κ∗(z)d2z=dxdy,z=x+iy

This is named the symplectic-transformed — translated wavelet transform. One can see that the mother wavelet ψ generates the other wavelets of the family ψ∗⁢[s⁢(z−κ)−r⁢(z∗−κ∗)] through a translating transform followed by a symplectic transform. (r,s are the symplectic transform parameter, |s|2−|r|2=1,κ is a translation parameter, s, r and κ∈C). This will be seen more clearly if we write the second transform in a matrix form

(431)(z−κz∗−κ∗)→M⁢(z−κz∗−κ∗)M≡(s−r−r∗s∗)

where M is a symplectic matrix satisfying MT⁢J⁢M=J, J=(0I−I0).

Symplectic matrices in Hamiltonian dynamics correspond to canonical transformations and keep the Poisson bracket invariant, while in matrix optics they represent ray transfer matrices of optical instruments, such as lenses and fibers.

14.1.1 Properties of symplectic-transformed — translated WT

It is straightforward to evaluate this transform and its reciprocal transform when f⁢(z) is the complex Fourier exponentials, f⁢(z)=exp⁡(z⁢β∗−z∗⁢β). (Note that z⁢β∗−z∗⁢β is pure imaginary)

(432)Wψ⁢f=s⁢∫d2⁢zπ⁢exp⁡[(z+κ)⁢β∗−(z∗+κ∗)⁢β]⁢ψ∗⁢(s⁢z−r⁢z∗)=exp⁡(κ⁢β∗−κ∗⁢β)⁢s⁢∫d2z′π×exp[z(s∗β∗−r∗β)′−z′⁣∗(sβ−rβ∗)]ψ∗(z)′=s⁢exp⁡(κ⁢β∗−κ∗⁢β)⁢Φ⁢(s∗⁢β∗−r∗⁢β)

where Φ is the complex Fourier transform of ψ∗. Then we form the adjoint operation

(433)Wψ∗⁢(Wψ⁢f)⁢(z)=s∗⁢∫d2⁢κπ⁢(Wψ⁢f)⁢(r,s;κ)×ψ⁢[s⁢(z−κ)−r⁢(z∗−κ∗)]=|s|Φ(s∗β∗−r∗β)∫d2⁢κπexp[(κ+z)β∗−(κ∗+z∗)β]ψ[−sκ+rκ∗]=|s|⁢exp⁡(z⁢β∗−z∗⁢β)⁢|Φ⁢(s∗⁢β∗−r∗⁢β)|2

from which we have

(434)∫Wψ∗⁢(Wψ⁢f)⁢(z)⁢d2⁢s|s|2=exp⁡(z⁢β∗−z∗⁢β)⁢∫d2⁢s⁢|Φ⁢(s∗⁢β∗−r∗⁢β)|2|s|

so we get the inversion formula

(435)f⁢(z)=ez⁢β∗−z∗⁢β=∫d2⁢s⁢Wψ∗⁢(Wψ⁢f)⁢(z)/|s|2∫d2⁢s⁢|Φ⁢(s∗⁢β∗−r∗⁢β)|2/|s|

Eq. (435) leads us to impose the normalization

(436)∫d2⁢s⁢|Φ⁢(s∗⁢β∗−r∗⁢β)|2/|s|=1

in order to get the wavelet representation

(437)f⁢(z)=∫d2⁢s⁢Wψ∗⁢(Wψ⁢f)⁢(z)/|s|2

Then we can have a form of Parseval’s theorem for this new wavelet transform:

Preposition: For any f and f′ we have

(438)∬Wψf(r,s;κ)Wψf′∗(r,s;κ)d2κd2s|s|2=∫d2z2πf(z)f′∗(z)

Proof: Let us assume F⁢(β) and F(β)′ are the complex Fourier transform of f⁢(z) and f(z)′ respectively,

(439)F⁢(β)=∫d2⁢z2⁢π⁢f⁢(z)⁢ez⁢β∗−z∗⁢β

Recall the convolution theorem defined on complex Fourier transform,

(440)∫d2zf(α−z,α∗−z∗)f(z)′=∫d2⁢z⁢∫d2⁢β2⁢π⁢F⁢(β)⁢e(α∗−z∗)⁢β−(α−z)⁢β∗×∫d2β′2⁢πF(β)′′ez∗β−′zβ′⁣∗=∫d2βF(β)F(β)′eα∗⁢β−α⁢β∗

Hence from Eqs. (440) and (429), as well as Eq. (430) we see that Wψ⁢f⁢(r,s;κ)=∫d2⁢zπ⁢f⁢(z)⁢ψr,s;κ∗⁢(z) can be considered as a convolution in the form

(441)∫d2zf(z)ψ∗[s(z−κ)−r(z∗−κ∗)]=∫d2βF(β)Φ∗(sβ−rβ∗)exp(κβ∗−κ∗β)

It then follows from Eq. (440) that

(442)∫Wψ⁢f⁢(r,s;κ)⁢Wψ⁢f′⁣∗⁢(r,s;κ)⁢d2⁢κ=|s|∫d2⁢β2⁢πF(β)F′∗(β)|Φ(sβ−rβ∗)|2

Therefore, using Eq. (436), we see that the further integration yields

(443)∫d2⁢s|s|2⁢∫Wψ⁢f⁢(r,s;κ)⁢Wψ⁢f′⁣∗⁢(r,s;κ)⁢d2⁢κ=∫d2⁢β2⁢π⁢F⁢(β)⁢F′⁣∗⁢(β)=∫d2⁢z2⁢π⁢f⁢(z)⁢f′⁣∗⁢(z)

which completes the proof.

Theorem: From the Proposition (438) we have

(444)∫∫Wψ⁢f⁢(r,s;κ)⁢ψr,s;κ⁢(z)⁢d2⁢κ⁢d2⁢s|s|2=f⁢(z)

that is, there exists an inversion formula for arbitrary function f⁢(z). In fact, in Eq. (430) when we take f(z)=δ(z−z)′, then

(445)Wψf(r,s;κ)=∫d2⁢zπf(z)ψr,s;κ∗(z)=ψr,s;κ∗(z)′

Substituting Eq. (445) into Eq. (443) we obtain Eq. (444).

14.1.2 Relation between Wψ⁢f⁢(r,s;κ) and optical Fresnel transform

Now we explain why the idea of Wψ⁢f⁢(r,s;κ) is originated from the optical Fresnel transform. We can visualize the symplectic-transformed — translated wavelet transform in the context of quantum mechanics. Let f⁢(z)≡⟨z|f⟩, ⟨z| is the coherent state, |z⟩=exp⁡[z⁢a†−z∗⁢a]≡|(zz∗)⟩,|0⟩ is the vacuum state in Fock space, then Eq. (429) can be expressed as

(446)Wψ⁢f⁢(r,s;κ)=s⁢∫d2⁢zπ⁢⟨ψ|(s−r−r∗s∗)⁢(z−κz∗−κ∗)⟩⁢⟨z|f⟩=⟨ψ|F1⁢(r,s,κ)|f⟩

where F(r,s,κ) is defined as

(447)F1⁢(r,s,κ)=s⁢∫d2⁢zπ⁢|s⁢z−r⁢z∗⟩⁢⟨z+κ|

and |s⁢z−r⁢z∗⟩≡|(s−r−r∗s∗)⁢(zz∗)⟩. Employing the normal ordering of the vacuum projector |0⟩⟨0|=:exp(−a†a): and the IWOP technique to perform the integration in Eq. (447), one can get the explicit form of F1⁢(r,s,κ). The transformation matrix element of F1⁢(r,s,κ=0) in the coordinate representation |x⟩ is just the kernel of optical diffraction integration (86) (Fresnel transform), which explains our motivation to introduce Wψ⁢f⁢(r,s;κ). In particular, when r∗=r≡sinh⁡λ,s=s∗≡cosh⁡λ, F1⁢(r,s,κ=0) reduces to the well-known single-mode squeezing operator exp⁡[λ2⁢(a2−a†2)], which corresponds to dilation in the usual WT.

14.2 Entangled SWT

In the above subsection, the mother wavelet is gained through a translating transform followed by a symplectic transform. This motivation arises from the consideration that symplectic transforms are more general than the dilated transform, and are useful in Fresnel transform of Fourier optics, e.g., ray transfer matrices of optical instruments, such as lenses and fibers in matrix optics, while in quantum optics symplectic transforms correspond to single-mode Fresnel operator (or generalized SU (1,1) squeezing operator).

Recalling that in Section 9 we have introduced the 2-mode entangled Fresnel operator which is a mapping of classical mixed transformation (z,z)′→(sz+rz′⁣∗,sz+′rz∗) in a two-mode coherent state |z,z⟩′ representation onto quantum operator F2⁢(r,s), we are naturally led to develop the SWT in Eq. (430) to the so-called entangled SWT (ESWT) [133] for signals g(z,z)′ defined in two complex planes,

(448)Wϕg(r,s;k,k)′=∬d2zd2z′π2g(z,z)′ϕr,s;k,k′∗(z,z)′

where

(449)ϕr,s;k.k′(z,z)′≡s∗ϕ[s(z−k)+r(z′⁣∗−k′⁣∗)×s(z−′k)′+r(z∗−k∗)]

is used as a weighting function to synthesize the signal g(z,z)′ regarding two complex planes. One can see that the mother wavelet φ generates the family ϕ∗[s(z−k)+r(z′⁣∗−k′⁣∗),s(z−′k)′+r(z∗−k∗)] through a translating transform followed by an entangled symplectic transform. We emphasize that this transform mixes the two complex planes, which is different from the tensor product of two independent transforms (z,z)′→[s(z−k)−r(z∗−k∗),s(z−′k)′−r(z′⁣∗−k′⁣∗)] given by Eq. (446). The new symplectic transform can be seen more clearly if we write it in a matrix form:

(450)(z−kz∗−k∗z−′k′z′⁣∗−k′⁣∗)⟶M⁢(z−kz∗−k∗z−′k′z′⁣∗−k′⁣∗)M=(s00r0s∗r∗00rs0r∗00s∗)

where M is symplectic satisfying MTJM=J,J=(0I−I0),I is the 2×2 unit matrix.

As Eq. (448) has been qualified as a new wavelet transform we must prove that it possesses fundamental properties of the usual wavelet transforms, such as the admissibility condition, the Parseval’s theorem and the inversion formula. When g(z,z)′ is the complex Fourier exponential,

(451)g1(z,z)′=exp(zβ∗−z∗β+zγ∗′−z′⁣∗γ)

According to Eqs. (448) and (449) we evaluate its ESWT

(452)Wϕg1=s∬d2zd2z′π2ϕ∗[sz+rz∗′,sz+′rz∗]×exp[(z+k)β∗−(z∗+k∗)β+(z+′k)′γ∗−(z′⁣∗+k′⁣∗)γ]

Making the integration variables transform s⁢z+r⁢z′⁣∗→w, sz+′rz∗→w,′ Eq. (452) becomes

(453)Wϕ⁢g1=s⁢exp⁡(k⁢β∗−k∗⁢β+k⁢γ∗′−k′⁣∗⁢γ)×∬d2wd2w′π2ϕ∗(w,w)′×exp[w(s∗β∗+r∗γ)−w∗(sβ+rγ∗)+w(s∗γ∗+r∗β)′−w′⁣∗(sγ+rβ∗)]

The last integration is just the complex Fourier transform (CFT) of ϕ∗. Denoting it as Φ∗, we have

(454)Wϕ⁢g1=s⁢exp⁡(k⁢β∗−k∗⁢β+k′⁢γ∗−k′⁣∗⁢γ)×Φ∗⁢(s∗⁢β∗+r∗⁢γ,s∗⁢γ∗+r∗⁢β)

Then we form the adjoint operation of Eq. (454),

(455)Wϕ∗(Wϕg1)(z,z)′=s∗∬d2kd2k′π2[(Wϕg)(r,s;k,k)′]ϕ[s(z−k)+r(z′⁣∗−k′⁣∗),s(z−′k)′+r(z∗−k∗)]=|s|2⁢Φ∗⁢(s∗⁢β∗+r∗⁢γ,s∗⁢γ∗+r∗⁢β)×exp(zβ∗−z∗β+zγ∗′−z∗′γ)∬d2vd2v′π2ϕ(v,v)′×exp[−v(s∗β∗+r∗γ)+v∗(sβ+rγ∗)−v(s∗γ∗+r∗β)′+v′⁣∗(sγ+rβ∗)]

where the integration in the last line is just the CFT of ϕ [comparing with Eq. (453)], thus Eq. (455) leads to

(456)Wϕ∗(Wϕg1)(z,z)′=|s|2exp(zβ∗−z∗β+zγ∗′−z′⁣∗γ)⋅|Φ⁢(s∗⁢β∗+r∗⁢γ,s∗⁢γ∗+r∗⁢β)|2

From Eq. (456) we have

(457)∫d2sWϕ∗(Wϕg1)(z,z)′/|s|4=exp⁡(z⁢β∗−z∗⁢β+z⁢γ∗′−z′⁣∗⁢γ)×∫d2s|Φ(s∗β∗+r∗γ,s∗γ∗+r∗β)|2/|s|2

which together with Eq. (451) lead to

(458)g1(z,z)′=∫d2sWϕ∗(Wϕg1)(z,z)′/|s|4∫d2⁢s⁢|Φ⁢(s∗⁢β∗+r∗⁢γ,s∗⁢γ∗+r∗⁢β)|2/|s|2

Eq. (458) implies that we should impose the normalization

(459)∫d2⁢s⁢|Φ⁢(s∗⁢β∗+r∗⁢γ,s∗⁢γ∗+r∗⁢β)|2/|s|2=1

such that the reproducing process exists

(460)g1(z,z)′=∫d2sWϕ∗(Wϕg1)(z,z)′/|s|4

Eq. (459) may be named the generalized admissibility condition. Now we have the corresponding Parseval theorem: For any g and g′, we have

(461)∭Wϕg(r,s;k,k)′Wϕ∗g′(r,s;k,k′)d2⁢k⁢d2⁢k⁢d2′⁢s|s|4=∬d2zd2z′g(z,z′)g′⁣∗(z,z)′

Proof: Assuming F⁢(β,γ) and F(β,γ)′ to be CFT of g(z,z)′ and g(z,z)′′, respectively,

(462)F(β,γ)=∬d2zd2z′π2g(z,z)′×exp⁡(z⁢β∗−z∗⁢β+z⁢γ∗′−z′⁣∗⁢γ)

Recall the corresponding convolution theorem

(463)∬d2zd2zg′(α−z,α∗−z∗;α−′z,′α∗′−z∗′)g(z,z)′′=∬d2βd2γF(β,γ)F(β,γ)′×exp⁡(α∗⁢β−α⁢β∗+α′⁣∗⁢γ−α⁢γ∗′)

and from Eqs. (448), (462) and (463), we see that Wϕg(r,s;k,k)′=∬d2⁢z⁢d2⁢z′π2g(z,z)′ϕr,s;k,k′∗(z,z)′ can be considered as a convolution in the form [Noting that the CFT of ϕ∗ is Φ∗, see Eqs. (453) and (454)]:

(464)∬d2zd2z′g(z,z′)ϕ∗[s(z−k)+r(z∗′−k′⁣∗),s(z−′k)′+r(z∗−k∗)]=∬⁢d2⁢β⁢d2⁢γ⁢F⁢(β,γ)⁢Φ∗⁢(s∗⁢β∗+r∗⁢γ,s∗⁢γ∗+r∗⁢β)×exp⁡(k⁢β∗−k∗⁢β+k′⁢γ∗−k′⁣∗⁢γ)

Using Eq. (464) we calculate

(465)∬Wϕg(r,s;k,k)′Wϕ∗g(r,s;k,k)′′d2kd2k′=|s|2⁢⨌⁢d2⁢β⁢d2⁢γ⁢d2⁢β⁢d2′⁢γ⁢F′⁢(β,γ)×Φ∗⁢(s∗⁢β∗+r∗⁢γ,s∗⁢γ∗+r∗⁢β)×F(β,′γ)′′⁣∗Φ(s∗β∗+r∗γ,s∗γ∗+r∗β)×δ(β−β)′δ(β∗−β′⁣∗)δ(γ−γ)′δ(γ∗−γ∗′)=|s|2⁢∬⁢d2⁢β⁢d2⁢γ⁢F⁢(β,γ)⁢F′⁣∗⁢(β,γ)×|Φ⁢(s∗⁢β∗+r∗⁢γ,s∗⁢γ∗+r∗⁢β)|2

As a consequence of Eqs. (459) and (465) the further integration yields

(466)∫d2⁢s|s|4∬Wϕg(r,s;k,k)′Wϕ∗g′(r,s;k,k)′d2kd2k′=∬⁢d2⁢β⁢d2⁢γ⁢F⁢(β,γ)⁢F′⁣∗⁢(β,γ)×∫d2s|Φ(s∗β∗+r∗γ,s∗γ∗+r∗β)|2/|s|2=∬⁢d2⁢β⁢d2⁢γ⁢F⁢(β,γ)⁢F′⁣∗⁢(β,γ)=∬d2zd2z′g(z,z)′g′⁣∗(z,z)′

which completes the proof.

Inversion formula: From Eq. (461) we have

(467)g(z,z)′=∭Wϕg(r,s;k,k)′ϕr,s;k,k′(z,z)′d2⁢k⁢d2⁢k⁢d2′⁢sπ2⁢|s|4

That is, there exists an inversion formula for g(z,z)′ which represents the original signal g(z,z)′ as a superposition of wavelet functions ϕr,s;k,k′, with the value of entangled wavelet transform Wϕg(r,s;k,k)′ serving as coefficients. In fact, in Eq. (448) when we take

(468)g(z,z)′=δ(z−u)δ(z∗−u∗)δ(z−′u)′δ(z∗′−u′⁣∗)

then

(469)Wϕg(r,s;k,k)′=1π2ϕr,s;k,k′∗(u,u)′

Substituting Eqs. (468) and (469) into Eq. (466), we obtain Eq. (467). We can visualize the ESWT in the context of quantum mechanics. Let g(z,z)′=⟨z,z|′g⟩ and use Eqs. (449) and (450), Eq. (448) is expressed as

(470)Wϕg(r,s;k,k)′=s∬d2zd2z′π2ϕ∗[s(z−k)+r(z∗′−k′⁣∗),s(z−′k)′+r(z∗−k∗)]g(z,z)′=s⁢∬⁢d2zd2z′π2⁢⟨ϕ|M⁢(z−kz∗−k∗z−′k′z′⁣∗−k′⁣∗)⟩×⟨z,z|′g⟩=⟨ϕ|F2(r,s;k,k)′|g⟩

where F2(r,s;k,k)′ is defined as

(471)F2(r,s;k,k)′=s∬d2zd2z′π2|sz+rz′⁣∗,sz+′rz∗⟩×⟨z+k,z+′k|′|sz+rz∗′,sz+′rz∗⟩=|sz+rz∗′⟩1⊗|sz+′rz∗⟩2

When k=0 and k′=0,F2(r,s;k=k′=0) is just the 2-mode Fresnel operator.

Thus, we have extended the SWT of signals in one complex plane to ESWT of signals defined in two complex planes. The latter is not the tensor product of two independent SWTs. This generalization is inevitable, since it resembles the extending from the single-mode squeezing transform (or Fresnel operator) to the two-mode squeezing transform (or entangled Fresnel operator) in quantum optics.

14.3 Symplectic-dilation mixed WT

Next we shall introduce a new kind of WT, i.e., symplectic-dilation mixed WT [134]. Recalling Ref. [135] we have constructed a new entangled-coherent state representation (ECSR) |α,x⟩,

(472)|α,x⟩=exp[−12x2−14|α|2+(x+α2)a1†+(x−α2)a2†−14(a1†+a2†)2]|00⟩

which is the common eigenvector of the operator (X1+X2)/2 and a1−a2, i.e., (a1−a2)⁢|α,x⟩=α⁢|α,x⟩ and 12⁢(X1+X2)⁢|α,x⟩=12⁢x⁢|α,x⟩, where Xi=12⁢(ai+ai†) is the coordinate operator, (i=1,2), and |α,x⟩ is complete,

(473)∫−∞∞d⁢xπ⁢∫d2⁢α2⁢π⁢|α,x⟩⁢⟨α,x|=1

and exhibits partly non-orthogonal property (for α) and orthonormal property (for x),

(474)⟨α,′x|′α,x⟩=πexp[−14(|α|2+|α|2′)+12αα′⁣∗]×δ(x−′x)

Therefore, |α,x⟩ possess the behaviors of both the coherent state and the entangled state. An interesting question is: Can we introduce a new kind of continuous WT which the |α,x⟩ representation underlies? The answer is affirmative. Our motivation of this issue comes from the mixed lens-Fresnel transform in classical optics [136].

By synthesizing Eqs. (363) and (430) and in reference to Eq. (473), we propose the mixed WT for g⁢(α,x) (α=α1+i⁢α2):

(475)Wψ⁢g⁢(s,r,κ;a,b)≡∫−∞∞d⁢xπ⁢∫d2⁢α2⁢π⁢g⁢(α,x)⁢ψs,r,κ;a,b∗⁢(α,x)

where d2⁢α=d⁢α1⁢d⁢α2. The family of mother wavelet ψ involves both the the symplectic transform of α and the dilation-transform of x,

(476)ψs,r,κ;a,b⁢(α,x)=s∗|a|⁢ψ⁢[s⁢(α−κ)−r⁢(α∗−κ∗),x−ba]

Let g⁢(α,x)≡⟨α,x|g⟩, then Eq. (475) can be expressed as the quantum mechanical version

(477)Wψ⁢g⁢(s,r,κ;a,b)=⟨ψ|U⁢(s,r,κ;a,b)|g⟩

where U⁢(s,r,κ;a,b) is defined as

(478)U⁢(s,r,κ;a,b)=s|a|⁢∫−∞∞d⁢xπ⁢∫d2⁢α2⁢π×|s⁢α−r⁢α∗,x−ba⟩⁢⟨α+κ,x|

U(s,r,κ=0;a,b=0) is just the generalized squeezing operator, which causes a lens-Fresnel mixed transform.

As Eq. (475) has been qualified as a new WT we must prove that it possesses some fundamental properties of the usual WTs, such as the admissibility condition, the Parseval theorem and the inversion formula. It is straightforward to evaluate the transform (475) and its reciprocal transform when g⁢(α,x) is the exponential g1⁢(α,x)=exp⁡(α∗⁢β−α⁢β∗−i⁢p⁢x),

(479)Wψ⁢g1=s|a|⁢eκ∗⁢β−κ⁢β∗−i⁢p⁢b⁢∫−∞∞d⁢xπ⁢∫d2⁢α2⁢π×ψ∗⁢(s⁢α−r⁢α∗,xa)⁢eα∗⁢β−α⁢β∗−i⁢p⁢x

Making the integration variables transform sα−rα∗→w,xa→x,′ which leads to d2⁢α→d2⁢w and ∫−∞∞dx→|a|∫−∞∞dx′, Eq. (479) becomes

(480)Wψ⁢g1=s⁢|a|⁢Φ∗⁢(s∗⁢β∗−r∗⁢β,a⁢p)⁢eκ∗⁢β−κ⁢β∗−i⁢p⁢b

where Φ∗ is just the Fourier transform of ψ∗,

(481)Φ∗(s∗β∗−r∗β,ap)=∫−∞∞dx′π∫d2⁢w2⁢πψ∗(w,x)′×ew∗(sβ−rβ∗)−w(s∗β∗−r∗β)−iapx′

Then we perform the adjoint WT of Eq. (475), using Eqs. (476) and (480) we see

(482)Wψ∗⁢(Wψ⁢g1)⁢(α,x)=s∗|a|⁢∫−∞∞d⁢bπ⁢∫d2⁢κ2⁢π⁢Wψ⁢g1×ψ⁢[s⁢(α−κ)−r⁢(α∗−κ∗),x−ba]=|s|⁢|a|⁢g1⁢(α,x)⁢Φ∗⁢(s∗⁢β∗−r∗⁢β,a⁢p)⁢∫−∞∞db′π×∫d2κ′2⁢πeκβ∗′−κ′⁣∗β+iapb′ψ(sκ−′rκ′⁣∗,b)′=|s|⁢|a|⁢g1⁢(α,x)⁢|Φ⁢(s∗⁢β∗−r∗⁢β,a⁢p)|2

From Eq. (482) we obtain

(483)∫−∞∞d⁢aa2⁢∫d2⁢s|s|2⁢Wψ∗⁢(Wψ⁢g1)⁢(α,x)=g1⁢(α,x)⁢∫−∞∞d⁢a|a|⁢∫d2⁢s|s|⁢|Φ⁢(s∗⁢β∗−r∗⁢β,a⁢p)|2

which leads to

(484)g1⁢(α,x)=∫−∞∞d⁢aa2⁢∫d2⁢s|s|2⁢Wψ∗⁢(Wψ⁢g1)⁢(α,x)∫−∞∞d⁢a|a|⁢∫d2⁢s|s|⁢|Φ⁢(s∗⁢β∗−r∗⁢β,a⁢p)|2

Eq. (484) implies that we should impose the normalization

(485)∫−∞∞d⁢a|a|⁢∫d2⁢s|s|⁢|Φ⁢(s∗⁢β∗−r∗⁢β,a⁢p)|2=1

such that the reproducing processexists

(486)g1⁢(α,x)=∫−∞∞d⁢aa2⁢∫d2⁢s|s|2⁢Wψ∗⁢(Wψ⁢g1)⁢(α,x)

Eq. (485) may be named the generalized admissibility condition. Now we have the corresponding Parseval theorem: For any g and g′, we have

(487)∫−∞∞d⁢a⁢d⁢ba2∫d2⁢κ⁢d2⁢s|s|2Wψg(s,r,κ;a,b)Wψ∗g(s,r,κ;a,b)′=∫−∞∞dx⁢∫d2⁢α⁢g⁢(α,x)⁢g′⁣∗⁢(α,x)

Proof: Assuming F⁢(β,p) and F(β,p)′ to be the Fourier transforms of g⁢(α,x) and g(α,x)′, respectively,

(488)F⁢(β,p)=∫−∞∞d⁢x2⁢π⁢∫d2⁢απ⁢g⁢(α,x)⁢eα⁢β∗−α∗⁢β+i⁢p⁢x

In order to prove Eq. (487), we first calculate Wψ⁢g⁢(s,r,κ;a,b). In a similar way of deriving Eq. (175), using Eqs. (170), (171) and the inversion formula of Eq. (488), we have

(489)Wψ⁢g⁢(s,r,κ;a,b)=s|a|⁢∫−∞∞d⁢p2⁢π⁢∫d2⁢βπ⁢F⁢(β,p)⁢∫−∞∞d⁢xπ⁢∫d2⁢α2⁢π×eα∗⁢β−α⁢β∗−i⁢p⁢x⁢ψ∗⁢[s⁢(α−κ)−r⁢(α∗−κ∗),x−ba]=s⁢|a|⁢∫−∞∞d⁢p2⁢π⁢∫d2⁢βπ⁢F⁢(β,p)×Φ∗⁢(s∗⁢β∗−r∗⁢β,a⁢p)⁢eκ∗⁢β−κ⁢β∗−i⁢p⁢b

It then follows that

(490)∫−∞∞db∫d2κWψg(s,r,κ;a,b)Wψ∗g(s,r,κ;a,b)′=|as|∫−∞∞dpdp′∫d2βd2βF′(β,p)F′∗(β,′p)′×Φ∗(s∗β∗−r∗β,ap)Φ(s∗β∗′−r∗β,′ap)′×∫−∞∞d⁢b2⁢π∫d2⁢κπ2eκ∗(β−β)′−κ(β∗−β′⁣∗)+i(p−′p)b=|a⁢s|⁢∫−∞∞dp⁢∫d2⁢β⁢F⁢(β,p)⁢F′⁣∗⁢(β,p)×|Φ⁢(s∗⁢β∗−r∗⁢β,a⁢p)|2

Substituting Eq. (490) into the left-hand side (L.H.S) of Eq. (487) and using Eq. (485), we see

(491)L.H.SofEq.(487)=∫−∞∞dp∫d2βF(β,p)F′∗(β,p)×∫−∞∞d⁢a|a|∫d2⁢s|s||Φ∗(s∗β∗−r∗β,ap)|2=∫−∞∞dp⁢∫d2⁢β⁢F⁢(β,p)⁢F′∗⁢(β,p)

Thus we complete the proof of Eq. (487).

Inversion Formula: From Eq. (487) we have

(492)g⁢(α,x)=∫−∞∞d⁢a⁢d⁢bπ⁢a2⁢∫d2⁢κ⁢d2⁢s2⁢π⁢|s|2⁢Wψ⁢g⁢(s,r,κ;a,b)×ψs,r,κ;a,b⁢(α,x)

which is the inversion formula for the original signal g⁢(α,x) expressed by a superposition of wavelet functions ψs,r,κ;a,b⁢(α,x), with the value of continuous WT Wψ⁢g⁢(s,r,κ;a,b) serving as coefficients. In fact, in Eq. (475), when we take g(α,x)=δ(α−α)′δ(α∗−α′⁣∗)δ(x−x)′, then

(493)Wψg(s,r,κ;a,b)=12⁢π⁢πψs,r,κ;a,b∗(α,′x)′

Substituting Eq. (493) into Eq. (491), we get (492).

We can visualize the new WT Wψ⁢g⁢(s,r,κ;a,b) in the context of quantum optics. Noticing that the generalized squeezing operator U(s,r,κ=0;a,b=0) in Eq. (478) is an image of the combined mapping of the classical real dilation transform x→x/a (a>0) and the classical complex symplectic transform (α,α∗)→(s⁢α−r⁢α∗,s∗⁢α∗−r∗⁢α) in |α,x⟩ representation, one can use the technique of integration within the normal product of operators to perform the integration in Eq. (478) to derive its explicit form [see Eq. (15) in Ref. [136]]. The transform matrix element of U(s,r,κ=0;a,b=0) in the entangled state representation |η⟩ is

(494)⟨η|U(s,r,κ=0;a,b=0)|η⟩′=sa∫−∞∞d⁢xπ∫d2⁢α2⁢π⟨η|sα−rα∗,xa⟩⟨α,x|η⟩′

In Fock space |η=η1+i⁢η2⟩ is a two-mode EPR entangled state in Eq. (28).

Then using Eqs. (472) and (28), we obtain

(495)⟨η|α,x⟩=12⁢exp⁡[−α2+|α|24−12⁢η12+η1⁢α−i⁢η2⁢x]

Substituting Eq. (495) into Eq. (494) and using Eq. (79), we obtain

(496)⟨η|U(s,r,κ=0;a,b=0)|η⟩′=πaδ(η2′η2/a)12⁢i⁢π⁢B×exp[i2⁢B(Aη12′−2η1η1′Dη12)]

which is just the kernel of a mixed lens – Fresnel transform, i.e., the variable η1 of the object experiences a generalized Fresnel transform, while η2 undergoes a lens transformation. Thus, based on |α,x⟩, we have introduced SDWT which involves both the real variable dilation-transform and complex variable symplectic transform, corresponding to the lens-Fresnel mixed transform in classical optics.

15 Fresnel–Hadamard combinatorial transformation

In the theoretical study of quantum computer, the Hadamard transform is of great importance. This operation is n Hadamard gates acting in parallel on n qubits. The Hadamard transform produces an equal superposition of all computational basis states. From the point of view of Deutsch – Jozsa quantum algorithm, the Hadamard transform is an example of the N=2n quantum Fourier transform, which can be expressed as [137, 138]

(497) |j⟩=12n⁢∑k=02n−1e2⁢π⁢i⁢j⁢k/2n⁢|k⟩

Now the continuous Hadamard transform, which used to go from the coordinate basis |x⟩ to the momentum basis, is defined as [139]

(498)F⁢|x⟩=1π⁢σ⁢∫−∞∞dy⁢exp⁡(2⁢i⁢x⁢y/σ2)⁢|y⟩

where σ is the scale length, and F is named Hadamard operator. Using the completeness of ∫−∞∞dx⁢|x⟩⁢⟨x|=1, we have

(499)F=1π⁢σ⁢∬−∞∞⁢d⁢x⁢d⁢y⁢exp⁡(2⁢i⁢x⁢y/σ2)⁢|y⟩⁢⟨x|

The above two transforms (Fresnel transform and Hadmard transform) are independent of each other. An interesting question thus naturally arises: can we combine the two transforms together? To put it in another way, can we construct a combinatorial operator which play the role of both Fresnel transform and Hadmard transform for two independent optical modes? The answer is affirmative. In this section we try to construct the so-called Fresnel–Hadmard combinatorial transform.

15.1 The Hadamard–Fresnel combinatorial operator

Based on the coherent-entangled representation |α,x⟩, and enlightened by Eqs. (81) and (498), we now construct the following ket–bra integration [140]

(500)U=sπ⁢σ⁢∫d2⁢απ⁢∬⁢dx⁢dy×exp⁡(2⁢i⁢x⁢y/σ2)⁢|s⁢α−r⁢α∗,y⟩⁢⟨α,x|

We name U the Hadamard–Fresnel combinatorial operator.

Substituting Eq. (498) into Eq. (500), and using the two-mode vacuum projector’s normally ordered form |00⟩⟨00|=:exp(−a1+a1−a2+a2): as well as the IWOP technique, we get

(501)U=sπ⁢σ:∫d2⁢zπ⁢A⁢(z,z∗)⁢∬⁢dx⁢dy⁢B⁢(x,y)⁢eC:

where

(501)C≡−(a1†+a2†)2+(a1+a2)24−a1†⁢a1−a2†⁢a2B(x,y)≡exp[−y2+x22+y(a1†+a2†)+x(a1+a2)+2⁢i⁢x⁢yσ2]

and

(501)A(z,z∗)≡exp[−|s⁢z−r⁢z∗|2+|z|24+s⁢z−r⁢z∗2(a1†−a2†)+z∗⁢(a1−a2)2]

They are all within the normal ordering symbol :: . Now perform the integration over d⁢x⁢d⁢y within ::, remembering that all creation operators commute with all annihilation operators (the essence of the IWOP technique) so that they can be considered c-number during the integration. We can finally obtain the normally ordered form of Hadamard–Fresnel combinatorial operator [140]

(502)U=1s∗4⁢π⁢σσ4+4:exp{−12rs∗(a1†−a2†2)2+σ4−42⁢(σ4+4)(a1†+a2†2)2+(1s∗−1)a1†−a2†2a1−a22+(4⁢i⁢σ2σ4+4−1)a1†+a2†2a1+a22+12r∗s∗(a1−a22)2+σ4−42⁢(σ4+4)(a1+a22)2}:

15.2 The properties of Hadamard–Fresnel operator

Note [a1−a2,a1†+a2†]=0 and [(a1−a2)/2,(a1†−a2†)/2]=1,[(a1+a2)/2,(a1†+a2†)/2]=1,(a1−a2)/2 can be considered a mode independent of another mode (a1†+a2†)/2, thus we have the operator identity

(503)exp⁡[f⁢(a1†±a2†)⁢(a1±a2)]=:exp[12(e2⁢f−1)(a1†±a2†)(a1±a2)]:

Using Eq. (503), we can rewrite Eq. (502) as

(504)U=U2⁢U1=U1⁢U2

where

(505)U1=4⁢π⁢σσ4+4⁢exp⁡[σ4−42⁢(σ4+4)⁢(a1†+a2†2)2]×exp⁡[a1†+a2†2⁢a1+a22⁢ln⁡4⁢i⁢σ2(σ4+4)]×exp⁡[σ4−42⁢(σ4+4)⁢(a1+a22)2]

and

(506)U2=exp⁡[−r2⁢s∗⁢(a1†−a2†2)2]×exp⁡[(a1†−a2†2⁢a1−a22+12)⁢ln⁡1s∗]×exp⁡[r∗2⁢s∗⁢(a1−a22)2]

while U2 is the Fresnel operator for mode (a1−a2)/2, and U1 is named the Hadamard operator for mode (a1+a2)/2.

It then follows that

(507)U⁢a1−a22⁢U−1=s∗⁢a1−a22+r⁢a1†−a2†2U⁢a1†−a2†2⁢U−1=r∗⁢a1−a22+s⁢a1†−a2†2

from which we see that the Hadamard–Fresnel combinatorial operator can play the role of Fresnel transformation for a1−a22. Physically, (a1−a2)/2 and (a1+a2)/2 can be two output fields of a beamsplitter.

In a similar way, we have

(508)Ua1+a22U−1=14⁢i⁢σ2[(σ4+4)a1+a22−(σ4−4)a1†+a2†2]Ua1†+a2†2U−1=14⁢i⁢σ2[−(σ4+4)a1†+a2†2+(σ4−4)a1+a22]

which for the quadrature Xi=(ai+ai†)/2,  Pi=(ai−ai†)/(2⁢i),i=1,2, leads to

(509)U⁢X1+X22⁢U−1=σ24⁢(P1+P2),U⁢(P1+P2)⁢U−1=−4σ2⁢X1+X22

from which we see that the Hadamard–Fresnel combinatorial operator also plays the role of exchanging the total momentum — average position followed by a squeezing transform, with the squeezing parameter being σ24.

The mutual transform in Eq. (509) can be realized by

(510)ei⁢π2⁢(a1†⁢a1+a2†⁢a2)⁢(X1+X2)⁢e−i⁢π2⁢(a1†⁢a1+a2†⁢a2)=P1+P2

(511)ei⁢π2⁢(a1†⁢a1+a2†⁢a2)⁢(P1+P2)⁢e−i⁢π2⁢(a1†⁢a1+a2†⁢a2)=−(X1+X2)

while the two-mode squeezing operator is S2=exp⁡[ln⁡2σ2⁢(a1†⁢a2†−a1⁢a2)], therefore

(512)U1=S2−1⁢ei⁢π2⁢(a1†⁢a1+a2†⁢a2)

From Eq. (504) and Eq. (512), we see that the Hadamard–Fresnel combinatorial operator can be decomposed as

(513)U=U2⁢S2−1⁢ei⁢π2⁢(a1†⁢a1+a2†⁢a2)=S2−1⁢ei⁢π2⁢(a1†⁢a1+a2†⁢a2)⁢U2

It can also be seen that U is unitary, and U+⁢U=U⁢U+=1.

In this section, we have introduced the Fresnel-Hadamard combinatorial operator by virtue of the IWOP technique. This unitary operator plays the role of both Fresnel transformation for mode a1−a22 and Hadamard transformation for mode a1+a22, respectively, and the two transformations are combinatorial. We have shown that the two transformations are concisely expressed in the coherent-entangled state representation as a projective operator in an integrative form. We also found that the Fresnel-Hadamard operator can be decomposed as U2⁢S2−1⁢ei⁢π2⁢(a1†⁢a1+a2†⁢a2), a Fresnel operator U2, a two-mode squeezing operator S2−1 and the total momentum-average position exchanging operator. Physically, a1−a22 and a1+a22 can be two output fields of a beamsplitter. If an optical device can be designed for Fresnel–Hadamard combinatorial transform, then it can be directly applied to these two output fields of the beamsplitter.

In summary, although quantum optics and classical optics are quite different, no matter in the mathematical tools they employed or in a conceptual view-quantum optics concerns the wave–particle duality of optical field with an emphasis on its nonclassical properties, whereas classical optics works on the distribution aqnd propagation of the light waves, it is a new exploration to link them systematically. In this review, via the route of developing Dirac’s symbolic method we have revealed some links between them by mapping classical symplectic transformation in the coherent state representation onto quantum unitary operators (GFO). Throughout our discussion the IWOP technique is indispensable for the derivation. We have resorted to the quantum optical interpretation of various classical optical transformations by adopting quantum optics concepts such as the coherent states, squeezed states, and entangled states, etc. Remarkably, we have endowed complex fractional Fourier transform, Hankel transform with quantum optical representation-transform interpretation. Our formalism, starting from quantum optics theory, not only provides quantum mechanical account of various classical optical transformations, but also have found their way back to some new classical transformations, e.g. entangled Fresnel transform, Fresnel-wavelet transform, etc, which may have realistic optical interpretation in the future. As Dirac predicted, functions that have been applied in classical optical problems may be translated in an operator language in quantum mechanics, and vice-versa. We expect that the content of this work may take its part in quantum states engineering, i.e., optical field states’ preparation and design.

Now the correspondence in this respect between the two distinct fields is established, the power of Dirac’s symbolic method can be fully displayed to solve some new problems in classical optics, e.g., to find new eigen-modes of some optical transforms; to extend the research region of classical optics theoretically by introducing new transforms (for example, the entangled Fresnel transforms), which may bring attention of experimentalists, who may get new ideas to implement these new classical optical transformations.

References

[1]

R. J. Glauber, Phys. Rev., 1963, 130: 2529

[2]

R. J. Glauber, Phys. Rev., 1963, 131: 2766

[3]

J. R. Klauder, Ann. Phys., 1960, 11: 123

[4]

J. R. Klauder and B.-S. Skagerstam, Coherent States, Singapore: World Scientific, 1985

[5]

J. R. Klauder and E. C. G. Sudarshan, Fundamentals of Quantum Optics, New York: W. A. Benjamin Inc., 1968

[6]

E. Schrödinger, Naturwiss, 1926, 14: 664

[7]

P. A. M. Dirac, The Principles of Quantum Mechanics, 3rd, Oxford: Clarendon Press, 1930

[8]

R. J. Glauber, in: Optique et Electronique Quantiques–Quantum Optics and Electronics, C. DeWitt, A. Blandin, and C. Cohen-Tannoudji (Eds.), New York: Gordon and Breach, 1965

[9]

L. Mandel and E. Wolf, Optical Coherence and Quantum Optics, Cambridge: Cambridge University Press, 1996

[10]

M. O. Scully and M. S. Zubairy, Quantum Optics, Cambridge: Cambridge University Press, 1997

[11]

D. F. Walls and G. J. Milburn, Quantum Optics, Berlin: Springer-Verlag, 1994

[12]

G. Nienhuis and L. Allen, Phys. Rev. A, 1993, 48: 656

[13]

D. Dragoman, Progress in Optics, 2002, 42: 424, and references therein

[14]

H. P. Yuen, Phys. Rev. A, 1976, 13: 2226

[15]

R. Loudon and P. L. Knight, J. Mod. Opt., 1987, 34: 709

[16]

D. F. Walls, Nature, 1983, 306: 141

[17]

M. C. Teich and B. E. A. Saleh, Quantum Opt., 1989, 1: 153

[18]

H. Y. Fan, H. R. Zaidi, and J. R. Klauder, Phys. Rev. D, 1987, 35: 1831

[19]

H. Y. Fan and H. R. Zaidi, Phys. Rev. A, 1988, 37: 2985

[20]

H. Y. Fan and J. R. Klauder, J. Phys. A, 1988, 21: L715

[21]

H. Y. Fan and H. Zou, Phys. Lett. A, 1999, 252: 281

[22]

H. Y. Fan, J. Opt. B, 2003, 5: R147

[23]

H. Y. Fan, Int. J. Mod. Phys. B, 2004, 18: 1387

[24]

H. Y. Fan, Int. J. Mod. Phys. B, 2004, 18: 2771

[25]

H. Y. Fan and J. R. Klauder, Phys. Rev. A, 1994, 49: 704

[26]

H. Y. Fan and X. Ye, Phys. Rev. A, 1995, 51: 3343

[27]

H. Y. Fan and A. Wünsche, J. Opt. B, 2000, 2: 464

[28]

H. Y. Fan, H. L. Lu, and Y. Fan, Ann. Phys., 2006, 321: 480

[29]

A. Wünsche, J. Opt. B, 1999, 1: R11

[30]

M. Born and E. Wolf, Principles of Optics, 7th edition, Singapore: World Scientific, 1999

[31]

J. W. Goodman, Introduction to Fourier Optics, New York: McGraw-Hill, 1972

[32]

A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev., 1935, 47: 777

[33]

H. Y. Fan, Phys. Lett. A, 2003, 313: 343

[34]

S. G. Lipson, H. Lipson, and D. S. Tannhauser, Optical Physics, 3rd Ed., Cambridge: Cambridge University Press, 1998

[35]

A. Torre, FrFT and Its Applications in Optics, Progress in Optics, Vol. 43, edited by E. Wolf, 2003, and references therein

[36]

V. Namias, J. Inst. Math. Appl., 1980, 25: 241

[37]

D. Mendlovic and H. M. Ozakatas, J. Opt. Soc. Am. A, 1993, 10: 1875

[38]

D. Mendlovic, H. M. Ozaktas, and A. W. Lohmann, Appl. Opt., 1994, 33: 6188

[39]

H. M. Ozakatas and D. Mendlovic, J. Opt. Soc. Am. A, 1993, 10: 2522

[40]

A. W. Lohmann, J. Opt. Soc. Am. A, 1993, 10: 2181

[41]

L. M. Bernardo and O. D. Soares, Opt. Commun., 1994, 110: 517

[42]

E. Wigner, Phys. Rev., 1932, 40: 749

[43]

S. A. Collins, J. Opt. Soc. Am., 1970, 60: 1168

[44]

S. Wang and D. M. Zhao, Metrix Optics, Berlin Heidelberg: Springer-Verlag, and Beijing: Higher Education Press, 2000

[45]

J. A. Arnaud, J. Opt. Soc. Am., 1971, 61: 751

[46]

C. Gomez-Reino, Int. J. Optoelectron., 1992, 7: 607

[47]

O. Seger, Ph. D. Dissertation 301, Linkoping University, 1993

[48]

T. Alieva and F. Agullo-Lopez, Opt. Commun., 1995, 114: 161

[49]

T. Alieva and F. Agullo-Lopez, Opt. Commun., 1996, 115: 267

[50]

D. F. V. James and G. S. Agarwal, Opt. Commun., 1996, 126: 207

[51]

H. Weyl, Z. Phys., 1927, 46: 1

[52]

H. Y. Fan and T. Q. Song, J. Phys. A, 2003, 36: 7803

[53]

H. Y. Fan, Phys. Rev. A, 2002, 65: 064102

[54]

H. Y. Fan, Phys. Lett. A, 2002, 294: 253

[55]

H. Y. Fan and Y. Fan, Phys. Rev. A, 1996, 54: 958

[56]

R. A. Campos, B. E. A. Saleh, and M. C. Teich, Phys. Rev. A, 1989, 50: 5274

[57]

C. Silberhorn, P. K. Lam, O. Weiss, F. König, N. Korolkova, and G. Leuchs, Phys. Rev. Lett., 2001, 86: 4267

[58]

E. C. G. Sudarshan, Phys. Rev. Lett., 1963, 10: 277

[59]

H. Y. Fan and B. Z. Chen, Phys. Rev. A, 1996, 53: 1948

[60]

H. J. Wu and H. Y. Fan, Mod. Phys. Lett. B, 1997, 11: 544

[61]

H. Y. Fan, H. Zou, and Y. Fan, Phys. Lett. A, 1999, 254: 137

[62]

D. Casasent and D. Psaltis, Proc. IEEE, 1977, 65: 77

[63]

H. Y. Fan, Opt. Commun., 2008, 281: 2023

[64]

A. Erdelyi, Higher Transcendental Functions, The Bateman Manuscript Project, New York: McGraw-Hill, 1953

[65]

H. Y. Fan and H. L. Lu, Phys. Lett. A, 2005, 334: 132

[66]

A. Gerrard and N. M. Burch, Introduction to Matrix Methods in Optics, London: John Wiley & Sons, 1975

[67]

H. Y. Fan, Commun. Theor. Phys., 2002, 38: 147

[68]

H. Y. Fan and J. H. Chen, Commun. Theor. Phys., 2003, 40: 589

[69]

H. Y. Fan and Z.-H. Xu, Phys. Rev. A, 1994, 50: 2921

[70]

H. Y. Fan and J. VanderLinde, Phys. Rev. A, 1989, 39: 2987

[71]

H. Y. Fan, Representation and Transformation Theory in Quantum Mechanics, Shanghai: Shanghai Scientific and Technical Publishers, 1997 (in Chinese)

[72]

H. Y. Fan, Entangled State Representations in Quantum Mechanics and Their Applications, Shanghai: Shanghai Jiao Tong University Press, 2001

[73]

H. Kogelnik and T. Li, Appl. Opt., 1966, 5: 1550

[74]

P. A. Bélanger, Opt. Lett., 1991, 16: 196

[75]

V. Magni, G. Cerullo, and S. D. Silvestri, Opt. Commun., 1993, 96: 348

[76]

H. Y. Fan and L. Y. Hu, Opt. Commun., 2008, 281: 1629

[77]

M. Nazarathy and J. Shamir, J. Opt. Soc. Am., 1982, 72: 398

[78]

M. Nazarathy and J. Shamir, J. Opt. Soc. Am., 1982, 72: 356

[79]

M. Nazarathy, A. Hardy, and J. Shamir, J. Opt. Soc. Am. A, 1986, 3: 360

[80]

H. Y. Fan and A. Wünsche, Commun. Theor. Phys., 2003, 39: 717

[81]

K. Vogel and H. Risken, Phys. Rev. A, 1989, 40: 2847

[82]

W. Vogel and W. P. Schleich, Phys. Rev. A, 1991, 44: 7642

[83]

U. Leonhardt, Measuring the Quantum State of Light, Cambridge: Cambridge University Press, 1997, and references therein

[84]

M. Aspelmeyer, Nat. Phys., 2009, 5: 11

[85]

I. L. Chuang and M. A. Nielson, Mod. Opt., 1997, 44: 2455

[86]

W. P. Schleich, Quantumm Optics in Phase Space, Berlin: Wiley-VCH, 2001

[87]

A. Wünsche, J. Mod. Opt., 1997, 44: 2293

[88]

A. Wünsche, Phys. Rev. A, 1996, 54: 5291

[89]

H. Y. Fan and L. Y. Hu, Opt. Commun., 2009, 282: 3734

[90]

R. F. O’Connell and E. P. Wigner, Phys. Lett. A, 1981, 83: 145

[91]

G. S. Agawal and E. Wolf, Phys. Rev. D, 1972, 2: 2161

[92]

G. S. Agawal and E. Wolf, Phys. Rev. D, 1972, 2: 2187

[93]

G. S. Agawal and E. Wolf, Phys. Rev. D, 1972, 2: 2206

[94]

M. Hillery, R. F. O’Connell, M. O. Scully, and E. P. Wigner, Phys. Rep., 1984, 106: 121

[95]

A. C. McBride, and F. H. Kerr, IMA J. Appl. Math., 1987, 39: 159

[96]

Y. G. Cai and S. Y. Zhu, Opt. Lett., 2005, 30: 388

[97]

S. Liu, J. Xu, Y. Zhang, L. Chen, and C. Li, Opt. Lett., 1995, 20: 1053

[98]

C. C. Shih, Opt. Lett., 1995, 20: 1178

[99]

L. M. Bernardo and O. D. Soares, Appl. Opt., 1996, 35: 3163

[100]

L. M. Bernardo and O. D. Soares, J. Opt. Soc. Am. A, 1994, 11: 2622

[101]

S. Chountasis, A. Vourdas, and C. Bendjaballah, Phys. Rev. A, 1999, 60: 3467

[102]

H. Y. Fan and L. Y. Hu, Chin. Phys. Lett., 2008, 25: 513

[103]

H. M. Ozaktas and M. F. Erden, Opt. Commun., 1997, 143: 75

[104]

H. Y. Fan and L. Y. Hu, J. Mod. Opt., 2009, 56: 1819

[105]

H. Lee, Phys. Rep., 1995, 259: 147

[106]

H. Y. Fan, Opt. Lett., 2003, 28: 2177

[107]

H. Y. Fan and H. L. Lu, Opt. Lett., 2006, 31: 2622

[108]

L. Y. Hu and H. Y. Fan, J. Mod. Opt., 2008, 55: 2429

[109]

H. Y. Fan, L. Y. Hu, and J. S. Wang, J. Opt. Soc. Am. A, 2009, 25: 974

[110]

H. Y. Fan and C. H. Lv, J. Opt. Soc. Am. A, 2009, 26: 2306

[111]

D. Dragoman, J. Opt. Soc. Am. A, 2009, 26: 274

[112]

H. Y. Fan, Commun. Theor. Phys., 1999, 31: 285

[113]

H. Y. Fan and Y. Fan, Commun. Theor. Phys., 2000, 33: 701

[114]

H. Y. Fan and L. S. Li, Commun. Theor. Phys., 1998, 29: 477

[115]

P. Pellat-Finet, Opt. Lett., 1994, 19: 1388

[116]

See: e.g., S. Jaffard, Y. Meyer, and R. D. Ryan, Wavelets, Tools for Science & Technology, SIAM, Philadelphia, 2001

[117]

M. A. Pinsky, Introduction to Fourier Analysis and Wavelets, Book/Cole, 2002

[118]

C. S. Burrus, R. A. Gopinath, and H. T. Guo, Introduction to Wavelets and Wavelet Transforms: A Primer, New Jersey: Prentice Hall, 1998

[119]

C. K. Chiu, Introduction to Wavelets, San Diego: Academic Press, 1992

[120]

D. R. Scifres, R. D. Burnham, and W. Streifer, Appl. Phys. Lett., 1978, 12: 33

[121]

E. Kapon, J. Matz, and A. Yariv, Opt. Lett., 1984, 10: 125

[122]

M. Oka, H. Masuda, Y. Kaneda, and W. Streifer, IEEE J. Quantum Electron., 1992, 28: 1142

[123]

I. Daubechies, Ten Lectures on Wavelets, CBMS-NSF Series in Applied Mathematics, Society for Industrial and Applied Mathematics, 1992

[124]

L. Y. Hu and H. Y. Fan, J. Mod. Opt., 2008, 55: 1835

[125]

P. S. Addison, The Illustrated Wavelet Transform Handbook: Introductory Theory and Applications in Science, Engineering, Medicine and Finance, London: Institute of Physics Publishing, 2002

[126]

P. Antoine, B. Piraux, and A. Maquet, Phys. Rev. A, 1995, 51: R1750

[127]

S. De Luca and E. Fiordilino, J. Phys. B, 1996, 29: 3277

[128]

K. Husimi, Proc. Phys.-Math. Soc. Jpn., 1940, 22: 264

[129]

C. L. Mehta, Phys. Rev. Lett., 1967, 18: 752

[130]

H. Y. Fan, Ann. Phys., 2007,

[131]

L. Y. Hu and H. Y. Fan, Int. J. Theor. Phys., 2009, 48: 1539

[132]

H. Y. Fan and H. L. Lu, Opt. Lett., 2006, 31: 3432

[133]

H. Y. Fan and S. G. Liu, Opt. Lett., 2007, 32: 1507

[134]

H. Y. Fan, S. G. Liu, and L. Y. Hu, Opt. Lett., 2009, 34: 551

[135]

H. Y. Fan and H. L. Lu, J. Phys. A, 2004, 37: 10993

[136]

H. Y. Fan, X. B. Tan, and H. L. Lu, Phys. Lett. A, 2006, 357: 163

[137]

M. A. Nielsen and I. L. Chuang, The Quantum Computation and Quantum Information, Cambridge: Cambridge University Press, 2000

[138]

See: e.g., J. Preskill, Quantum Information and Computation, California Institute of Technology, 1998

[139]

S. Parker, S. Bose, and M. Plenio, Phys. Rev. A, 2000, 61: 032305

[140]

C. M. Xie and H. Y. Fan, J. Mod. Opt., 2010, 57: 582

Rights & permissions

© Higher Education Press and Springer-Verlag Berlin Heidelberg 2012

PDF (696KB)

2308

Accesses

0

Citation

Detail

Sections
Recommended

/

〈 〉