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
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
. 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
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
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
we point out that Wigner operator’s Radon transformation is just the pure state projection operator
, where
and
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
. 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
. 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,
, 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
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
in 2-mode coherent state
representation. At the end of this section, we introduce a new symplectic-dilation mixed WT by employing a new entangled-coherent state representation
. 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
. This unitary operator plays the role of both Fresnel transformation for mode
and Hadamard transformation for mode
, 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
where is the optical distribution of a 2D light source and is its image on the observation plane, is the optical wavelength, is the wave number in the vacuum and is the propagation distance. When is satisfied, Eq. (1) is reduced to
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
. Then, the light distribution observed at the plane
corresponds to the
α equal to the
-th FrFT of the input fraction, where
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
We can see that
F0 is the identity operator and
is just the Fourier transform. The most important property of FrFT is that
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
ray transfer matrix is the Collins diffraction integral formula [
43],
where
if the system is lossless. One can easily find the similarity between Collins formula and the FrFT by some scaling transform and relating the
matrix to
in the FrFT [
44]. Note that
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]
where and . When has rotational symmetry
then Eq. (5) becomes
where we have used the m-order Bessel function
When , Eq. (7) reduces to the standard Hankel transform (up to a phase factor)
The compact form of 1D Collins formula is
where the transform kernel is
M is the parameter matrix
. 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.,
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 are standing on the left of annihilation operators a in a monomial of and a.
(i) The order of Bose operators a and within a normally ordered product can be permuted. That is to say, even though , we can have: 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 has the normal product form
(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 for convenience.), using the Fock representation of the coordinate eigenvector
we perform the integration below:
Substituting Eq. (12) into Eq. (14), we see
Note that on the left of : 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
As a commutes with within : :, and a can be considered as if they were parameters while the integration is being performed. Therefore, by setting , sech tanh we are able to perform the integration and obtain
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 (1,1) Lie algebra method in the derivation. The integration automatically arranges the squeezing operator in normal ordering. With
Eq. (17) becomes
This shows that the classical dilation 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 Eq. (19) becomes
a real simple Gaussian integration! This immediately leads us to put the completeness relation of the momentum representation into the normally ordered Gaussian form
where and is the momentum eigenvector ,
In addition, we should notice that
and
are related by the Fourier transform (FT), i.e.,
. 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
which is just the normally ordered Wigner operator since its marginal integration gives and respectively, i.e.,
Thus the Wigner function of quantum state ρ can be calculated as . On the other hand, the Wigner operator (23) can be recast into the coherent state representation,
where is a coherent state. In fact, using the IWOP technique we can obtain
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
and the total momentum
are commutable. Hinted by EPR, the bipartite entangled state
is introduced as [
53–
55]
is the common eigenstate of relative coordinate and the total momentum ,
where . Using the IWOP technique, we can immediately prove that possesses the completeness relation
and orthonormal relation
The Schmidt decomposition of is
The state can also be Schmidt-decomposed in momentum eigenvector space as
The
is physically appealing in quantum optics theory, because the two-mode squeezing operator has its natural representation on
basis [
55]
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
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],
which obeys the eigen-equations
where . Because we name the conjugacy between and . The completeness and orthonormal relations of are
respectively. and can be related to each other by
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:
where
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]
For two-mode correlated system, it prefers to using to calculate quantum states’ Wigner function. For example, with the two-mode squeezed states’ Wigner function is
4 Two deduced entangled state representations and Hankel transform
4.1 Deduced entangled states
Starting from the entangled state
and introducing an integer
m, we can deduce new states [
61],
which is worthy of attention because when we operate the number-difference operator, on , using Eq. (28) we see
so the number-difference operator corresponds to a differential operation in the representation, which is a remarkable property of . It then follows that
On the other hand, by defining we see and is its eigenstate,
where K is named correlated-amplitude operator since . Thus we name correlated-amplitude— number-difference entangled states. It is not difficult to prove the completeness and orthonormal property of ,
On the other hand, from we can derive another state
which satisfies
So D in representation is equal to . Consequently,
Note and
is qualified to be a new representation since
4.2 Hankel transform between two deduced entangled state representations
Since
and
are mutual conjugate,
is the conjugate state of
. From the definition of
and
and Eq. (40) we calculate the overlap [
33]
where we have identified the generating function of the s-order Bessel function and
Eq. (56) is remarkable, because is just the integral kernel of Hankel transform. In fact, if we define
and use Eq. (49) as well as Eq. (56), we obtain
which is just the Hankel transform of [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
Now we know that the quantum optical image of classical Hankel transform just corresponds to the representation transformation between two mutually conjugate entangled states and 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,
or the correlated-amplitude — number-difference entangled state can be considered as circular harmonic decomposition of . Let , a general 2D function expressed in polar coordinates, be periodic in the variable it can be looked as the wavefunction of the state vector in the representation
Using Eq. (61) we have
where is the wavefunction of in representation. We see from Eq. (46)
so behaves a rotation operator in representation, the expectation value of in is
which is just the cross-correlation between and an angularly rotated version of the same function, . On the other hand, using Eq. (49) we have
Substituting Eq. (66) into Eq. (65) and using the eigenvector equation (47) as well as Eq. (50), we obtain
from which we see that each of the circular harmonic components of the crosscorrelation undergoes a different phase shift so is not rotation invariant. However, when we consider only one harmonic component
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
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
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 with a definite quantum mechanical meaning, i.e., the overlap between and the rotated state 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
is both rotated and squeezed (by a two-mode squeezing operator
, then from Eqs. (35) and (62) we have
it follows the overlap between and the state
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
and therefore
from which one can see that to achieve simultaneous scale and rotation invariance, a two-dimensional object
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
, 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]
which is quite different from the product of two single-variable Hermite polynomials, and its generating function formula is
and noting we can directly perform the integral in Eq. (51) and derive the explicit form of
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
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
in quantum optics [
1–
5]
which is the eigenstate of annihilation operator , using the IWOP and (12), we can put the over-completeness relation of into normal ordering . The canonical form of coherent state is expressed as
where . It follows that which indicates that the states generate a canonical phase-space representation for a state . 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 are elements of a ray transfer matrix describing optical systems, belongs to the unimodular symplectic group; and the coherent state is a good candidate for providing with classical phase-space description of quantum systems, we naturally think of that the symplectic transformation 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
as the GFO. In fact, using notation of (coherent state), and introducing complex numbers s and
from Eq. (77) we know
is still a symplectic group element, so Eq. (78) becomes [
67–
70]
where the factor is attached to anticipate the unitarity of the operator F1. Eq. (78) tells us that c-number transform in coherent state basis is mapped into . Now we prove is really the FO we want. Using the IWOP technique and Eqs. (80) and (12) we can perform the integral
where we have used the mathematical formula [
71,
72]
with the convergent condition , . It then follows
Then using
and the completeness relation of coherent state as well as Eq. (79), we obtain the matrix element of in coordinate representation ,
which is just the kernel of generalized Fresnel transform (GFT) 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 , and using Eq. (16), we can rewrite Eq. (10) as
which is just the quantum mechanical version of GFO. Therefore, the 1D GFT in classical optics corresponds to the 1-mode GFO operating on state vector in Hilbert space, i.e., . One merit of GFO is: using coordinate-momentum representation transform we can immediately obtain GFT in “frequency” domain, i.e.,
Obviously, induces the following transform
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 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), and the IWOP technique we can directly perform the following integrals
where we have set
or
From which we see that it is just the mapping of the above matrices multiplication. Hence 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
seems as if their overlap 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 .
As a result of this group multiplication rule of GFO, we immediately obtain
provided that the parameter matrices 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 law for optical propagation — single-mode case
In classical optics, ray-transfer matrices, 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 and output ray i.e.,
where r1 is the ray height from the optical axis, and is named the optical direction-cosine; specifies the ray’s wavefront shape. Eq. (95) implies
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 , respectively, then the whole system is determined by a matrix .
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
law via optical diffraction integration, i.e., the input light field
and output light field
are related to each other by so-called Fresnel integration [
31]
where
The ABCD law for Gaussian beam passing through an optical system is [
74,
75]
where
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
law also exhibit in quantum optics? Since classical Fresnel transform should have its quantum optical counterpart?
To see the law more explicitly, using Eq. (79) we can re-express Eq. (82) as
and the multiplication rule for F1 is where
Next, we directly use the GFO to derive law in quantum optics. From Eq. (99) we see that the GFO generates
If we identify
then
The solution of Eq. (101) is
Let expressed by Eq. (102) be an input state for an optical system which is characteristic by parameters then the quantum optical ABCD law states that the output state is
which has the similar form as Eq. (102), where is determined by Eq. (99), and
which resembles Eq. (97).
Proof: According to the multiplication rule of two GFOs and Eqs. (98) and (99), we have
Using Eq. (103) we see along with Eq. (105) we can reach Eq. (104), thus the law is proved. Using Eq. (102) we can re-express Eq. (105) as
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
where we have set the initial mass γ denotes damping. Using to perform the transformation
then leads to
where . does not contain t explicitly. The dynamic evolution of a mass-varying harmonic oscillator from the Fock state at initial time to a squeezed state at time t is
if we let and then , according to Eq. (104) we directly obtain
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
, the GFT changes into the well-known Fourier transform, which is adapted to express mathematically the Fraunhofer diffraction. And if
, 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
can also be expressed in the form of quadratic combination of canonical operators
Q and
P [
80], i.e.,
where we have set . To confirm this, we first calculate matrix element
where we have used the squeezing property
It then follows from Eq. (114) and , we have
Thus in Eq. (113) is really the expected GFO. Next we directly use Eq. (78) and the canonical operator representation (113) to develop the optical operator method.
Noticing the matrix decompositions [
44]
and comparing Eqs. (78) and (113) as well as using (99) we know
where
which are named quadrature phase operator, Fresnel propagator in free space, and scaling operator (squeezed operator [
14–
17]), respectively. When
from Eq. (98) we see
which is named the Fourier operator, since it quantum mechanically transforms [
71,
72]
7.2 Alternate decompositions of GFO
Note that when , the decomposition (113) is not available, instead, from
and Eqs. (98), (99) and (113) we have
It then follows that
Besides, when we notice
and
we have another decomposition for
and for ,
7.3 Some optical operator identities
For a special optical systems with the parameter
we have
In particular, when from
we have
Using the following relations
it then follows that
While from
we obtain
so we have
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
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
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,
of the quadrature amplitude. Once the distribution
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 to the Fresnel quadrature phase
where s and r are related to ABCD through Eq. (79), we shall prove that the related Radon transform of Wigner operator is just the pure state density operator (named as the tomographic density operator) formed with the eigenstates belonging to the quadrature QF ( Q is the coordinate operator),
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
so from we have indeed
Second, we can derive the explicit form of . Starting from we set up the eigenvector equation
It follows that
In the coordinate and momentum representations we have
The normalizable solutions to Eqs. (146) and (147) are
Using the Fock representation of and in Eqs. (13) and (22), we obtain
and
Comparing Eq. (150) with Eq. (151) we see
On the other hand, according to the orthogonalization of , we have
Thus, combining Eqs. (152) and (153) we deduce
and
or
It is easily seen that that
make up a complete set
(so
can be named as the tomography representation). Then according to the Weyl quantization scheme [
51]
where is the Weyl correspondence of
is the Wigner operator [
42,
90–
94],
and using Eqs. (158), (159) and (150) we know that the classical Weyl correspondence (Weyl image) of the projection operator is
which means
Combining Eqs. (142)–(145) together we complete the proof. Therefore, the probability distribution for the Fresnel quadrature phase is the Radon transform of WF
So we name the tomographic density. Moreover, the tomogram of quantum state is just the squared modulus of the wave function . This new relation between quantum tomography and optical Fresnel transform may provide experimentalists to figure out new approach for generating tomography.
The introduction of also brings convenience to obtaining the inverse of Radon transformation. Using the complete relation of and Eq. (144) we have
Considering its right hand-side as a Fourier transformation, its reciprocal transform is
where . So once the distribution 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 related Radon transform of is just the pure state density operator formed with the eigenstates belonging to the conjugate quadrature of
Similarly, we find that for the momentum density,
where
As an application of the relation (141), recalling that the makes up a faithful representation of the symplectic group, it then follows from Eq. (141) that
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 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 . One can see
The right hand side of Eq. (169) is commonly defined as the tomogram of quantum states in direction, so in our view the calculation of tomogram in direction is ascribed to the calculation of
This is a concise and neat formula. Similarly, the tomogram in direction is ascribed to .
According to the Weyl correspondence rule
and the Weyl ordering form of
where the symbol denotes Weyl ordering, the classical correspondence of a Weyl ordered operator is obtained just by replacing in i.e.,
Let the classical Weyl correspondence of be
then using Eqs. (170) and (159) we have
Thus we reach a theorem as below.
The tomogram of a density operator ρ is equal to the marginal integration of the classical Weyl correspondence of where F is the Fresnel operator, expressed by
or
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 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
through the following 2-mode coherent state representation [
67–
70]
which indicates that is a mapping of classical sympletic transform in phase space. Concretely, the ket in Eq. (177) is
s and r are complex and satisfy the unimodularity condition. Using the IWOP technique we perform the integral in Eq. (177) and obtain
Thus induces the transform
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 to two-dimensional GFT. Letting and then project onto the entangled state defined by Eq. (28), and using the completeness relation (30) of , we obtain
Then using the overcompleteness relation of the coherent state, we can calculate the integral kernel
where the relation between and in Eq. (79) is used and the superscript only means that the parameters of are , and the subscript means the 2D kernel. Eq. (182) has the similar form as Eq. (11) except for its complex form. Taking and , , we have
This shows that is really the counterpart of the 2D GFT.
If taking the matrix element of in the representation which is conjugate to , we obtain the 2-dimensional GFT in its “frequency domain”, i.e.,
where the superscript N means that this transform kernel corresponds to the parameter matrix .
The two-mode GFO also abides by the group multiplication rule. Using the IWOP technique and Eq. (177) we can also obtain
where 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
on two-mode number state
and using the overlap between coherent state and number state, i.e.,
and the integral formula [
71,
72]
we can calculate
where
is the two variables Hermite polynomial [
87,
88], shown in Eqs. (73) and (74).
Using Eqs. (79) and (103), we recast Eq. (187) into
Taking notice of the multiplication rule of in Eq. (185), which is equivalent to
where and 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
Similar to the way of deriving Eq. (188), we can simplify Eq. (190) as
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 and Eqs. (29), (30), and we can prove the operator identity
where we have set
When and using Eq. (179) we see that
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
where
It then follows from Eqs. (179) and (195) that
which is corresponding to Fresnel propagator in free space (single-mode case). In particular, when , and , Eq. (179) becomes
where . Eq. (198) is just the two-mode squeezing operator,
which actually squeezes the entangled state (its conjugate state is ),
Using the decomposition (117) of the matrix and combining equations (194), (197) and (199) together, we see that
This is the two-mode quadratic canonical operator representation of .
To prove Eq. (202), using Eqs. (201) and (40) we have
It then follows that
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 obey the commutative relations so involves a SU (2) Lie algebra structure. (This structure is also compiled by and that have been used in constructing ).
9.3.2 Alternate decompositions of GFO and new optical operator identities
When from Eq. (179) we see
which can also be named the Fourier operator, since it induces the quantum mechanically transforms [
71,
72]
It then follows that
On the other hand, in order to obtain the decomposition of for , similar to deriving Eq. (126), we have
While for or , using Eqs. (127) and (128) we have another decomposition of , i.e.,
and
In addition, noticing Eqs. (137) and (135), we can rewrite Eqs. (209) and (210) as follows:
and
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 , which corresponds to the Fourier transform system, we have
In particular, when , and with Eq. (133) corresponding to the ideal spectrum analyzer, we have
When
which corresponds to the form of the image system, another operator identity is given by
When , which corresponds to the far foci system,
When corresponding to the Fresnel transform system,
we have
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 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 becomes the tomographic density , which is just the Radon transform of the Wigner operator . In this section we want to generalize the above conclusion to the two-mode entangled case. Here we shall prove
That is, we show that is just the Radon transform of the entangled Wigner operator
Similar to the single-mode case, operating on entangled state representation we obtain ()
or
where we have used the integration formula
Noticing the completeness relation and the orthogonality of , we immediately derive
a generalized entangled state representation with the completeness relation (222). From Eq. (220) we can see that
so we have the eigen-equations for as follows:
We can also check Eqs. (223)–(226) by another way.
9.4.1 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]
where is the Weyl correspondence of , and
The classical Weyl correspondence of the projection operator can be calculated as
Then using Eq. (182) we have
which means the following Weyl correspondence
so the projector operator is just the Radon transformation of , and and 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
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 related Radon transform of entangled Wigner operator is just the pure state density operator i.e.,
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.,
where the superscript N means that this transform kernel corresponds to the parameter matrix . Thus the 2D Fresnel transformation in its “frequency domain” is given by
Operating on we also have
or
Noticing that the entangled Wigner operator in representation is expressed as
and using the classical correspondence of which is calculated by
we obtain
so the projector operator is another Radon transformation of the two-mode Wigner operator, with and being the Radon transformation parameter (“frequency domain”). Therefore, the quantum tomography in representation is expressed as the Radon transformation of the Wigner function
and
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 is
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
where and
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 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 , i.e.,
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 is set up, based on which the tomogram of quantum state is just the squared modulus of the wave function . 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
). 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
where the exponential function is an integral kernel. In order to find the quantum correspondence of FrFT, multiply the function by the ket and bra from left and right, respectively, where and are coordinate eigenvectors, , and then use Eq. (13) and the IWOP technique to perform the integration. We obtain
where we have used the operator identity in the last step of Eq. (247).
From the orthogonal relation , we know that Eq. (247) indicates
which implies that the integral kernel in Eq. (246) is just the matrix element of operator
in coordinate state
called as Fractional Fourier Operator [
70]. Therefore, if we consider
as
, the wave function of quantum state
in the coordinate representation, from Eqs. (246) and (248) it then follows that
which suggests
From Eqs. (249) and (246) one can see that the FrFT in Eq. (246) corresponds actually to the rotating operator 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 (the additivity property). Without losing generality, we examine
According to the completeness relation of coordinate eigenvector, Eq. (251) yields
which is just the additivity of FrFT.
In particular, when is the number state, , its wavefunction in coordinate representation is
The FrFT of is , or
which indicates that the eigenfunction is Hermite-Gaussian function with the corresponding eigenvalue being .
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]
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]
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 Eq. (256) is reduced to
By defining and we can recast (257) into the simple form
According to Eq. (98) we immediately know that the operator of SFrFT is
Noting that the matrix can be decomposed into
according to the previous section, we have
where and and are the quadrature phase operator, the squeezing operator and the free propagation operator, respectively. On the other hand, from , and Eq. (261) we see
It then follows that
Using the canonical operator form (261) or (263) of we can deduce its matrix element in the coordinate states (its conjugate state is )
which is just the kernel of SFrFT. Thus we name SFrFT operator. Noticing that the and construct a close SU (2) Lie algebra, we can put Eq. (263) into a more compact form, i.e.,
Eqs. (261), (263) and (265) are different forms of the same operator of SFrFT; especially when , which is the usual FrFT operator.
Using Eq. (264) the SFrFT of denoted as , can be expressed as an matrix element in quantum optics context,
The above discussions are useful since any unimodular matrix
can be decomposed into [
103]
where the parameters are all real,
Correspondingly, the operator of FrFT is given by
where is the quadratic phase operator. Thus the general Fresnel transform can always be expressed by SFrFT as follows:
i.e., the output is the SFrFT of the input plus a quadratic phase term .
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
which differs from the usual two-fold Fourier transformation . In particular, when Eq. (271) is reduced to
so can be considered a basis funtion in phase space, or Eq. (271) can be looked as an expansion of with the expansion coefficient being . We can prove that the reciprocal transformation of Eq. (271) is
In fact, substituting Eq. (271) into the left-hand side of Eq. (273), we get
This transformation’s Parseval-like theorem can be proved
Now we apply Eq. (271) to phase space transformation in quantum optics. Recall that a signal
’s Wigner transform [
86,
90–
94,
105] is
Use Dirac’s symbol [
7] to write
.
is the eigenvector of coordinate
, and the Wigner operator emerges from Eq. (276),
If
is quantized as the operator
through the Weyl–Wigner correspondence [
51]
then
which in the literature is named the Weyl transform. is the Weyl classical correspondence of the operator . Substituting Eq. (279) into Eq. (271), we have
Using and where is the momentum eigenvector, and
then Eq. (280) becomes
Thus through the new integration transformation a new relationship between a phase space function and its Weyl–Wigner correspondence operator has been revealed. The inverse of Eq. (281), according to Eq. (273), is
For example, when its classical correspondence is
Substituting Eq. (283) into Eq. (281), we have
Using the Gaussian integration formula
in particular, when , with , Eq. (285) becomes
where represents an infinite long chirplet function. Comparing Eq. (286) with Eq. (284) we see , it then follows that
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
so the transformation (271) is equivalent to
For example, using Eqs. (277) and (288), we have
so
Thus this new transformation can convert the Wigner operator to . Similarly, we have
Then for the Wigner function of a density operator , we have
We may define as the –characteristic function. Similarly,
We name as the – 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
in the entangled state
as the integral transform kernel of CFrFT,
Using the normally ordered expansion of and the completeness relation of the coherent state representation, we calculate that is
which is just the integral kernel of the CFrFT in Ref. [
107]. Thus we see that the matrix element of
between two entangled state representations
and
corresponds to CFrFT. This is a new route from quantum optical transform to classical CFrFT transform. Let
, Eq. (295) becomes
In fact, let . Then use Eqs. (30) and (296), and we have
Thus the quantum mechanical version of CFrFT is given by
The standard complex Fourier transform is . 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
(the additivity property). For this purpose, using Eq. (309) and Eq. (30), we see
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 is a two-mode number state , then the CFrFT of the wave function is
To calculate , let us recall the definition of two-variable Hermite polynomial , and we can expand as
thus
As a result of Eq. (303), we see that Eq. (301) becomes
If we consider the operation as an operator, one can say that the eigenfunction of (the eigenmodes of CFrFT) is the two-variable Hermite polynomials with the eigenvalue being . This is a new property of CFrFT. Since the function space spanned by is complete,
and
one can confirm that the eigenmodes of CFrFT form an orthogonal and complete basis set [
109]. Note that the two variable Hermite polynomial
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
When Eq. (307) becomes
so can be considered a basis function in phase space, or Eq. (307) can be looked as an expansion of in terms of with the expansion coefficient being
We can prove that the inverse transform of Eq. (307) is
In fact, substituting Eq. (307) into the left-hand side of Eq. (309), and we obtain
The Parseval-like theorem for this transformation can also be demonstrated as
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),
The advantage of introducing
can be seen in Ref. [
80]. The corresponding Wigner function for a density matrix
ρ is
If is quantized as the operator through the Weyl–Wigner correspondence
then using Eq. (313) we see
which is named as the complex Weyl transform, and is the Weyl classical correspondence of . Substituting Eq. (315) into Eq. (307), we get
Using Eq. (29), we have
As a result of Eq. (317) and we see
Using Eq. (318), we convert Eq. (316) as
The inverse of Eq. (319), according to Eq. (309), is
Thus through the new integration transformation, a new relationship between a phase space function and its Weyl–Wigner correspondence operator is revealed.
For example, from the following Weyl–Wigner correspondence
[ and are defined in Eqs. (193) and (196)] and Eq. (320) we have
Using the Gaussian integration formula
in particular, when with , Eq. (323) becomes
where represents an infinite long chirplet function. By taking in Eq. (322), such that and comparing with Eq. (324) we obtain
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.,
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],
we have
which implies that the integral kernel in Eq. (326) is just the matrix element of operator in coordinate states. From Eq. (328) it then follows that
which suggests
From Eqs. (329) and (326) one can see that the scaled FrFT in Eq. (329) corresponds actually to the squeezing-rotating operator 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
where . Using the natural expression of the two-mode squeezing operator [Eq. (34)], and noticing that is just the integral kernel of CFrFT (296), we can reform (331) as
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]
We see that
Using Eq. (30), we have
where . Substituting Eqs. (333) and (334) into Eq. (335), we can derive
where
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 we can define an entangled scaled FrFT, i.e.,
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,
Inserting the completeness relation of into Eq. (338), we get
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
rule is equally derived via optical diffraction integral theory — the Collins integral formula. As shown in Fig. 2, if
represents the input field amplitude at point
η on
S1, and
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
Take and write out
where K and L are two constants to be determined later. Then according to Eqs. (181) and (340) we have
Comparing Eq. (342) with Eq. (182) leads us to choose
Then Eq. (342) becomes
so Eq. (344) is a standard CFrFT up to a quadratic phase term , according to Eq. (343) and , it can also be written as
which represents a quadratic approximation to a sphere wave diverging from a luminous point at the distance
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
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 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 (see Fig. 1), the additivity means that the diffraction pattern observed on () (the sphere tangent to S3 with radius ) and associated with should be the result of a first diffraction phenomenon (associated with on S (with ), followed by a second diffraction phenomenon (associated with from S to . 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 (see Fig. 2) described by the ray transfer matrix []. Similar to deriving Eq. (346), after the squeezing transform and the phase compensation,
Thus we can obtain the expression of CFrFT for Collins diffraction from S1 to (not S3),
where and
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 is [see Eqs. (181), (344)–(346)]
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 []. For this purpose, using the group multiplication rule of , we can decompose the diffraction from S1 to into two parts: one is described as the matrix from the plane S1 (with η) to (with ), the other is from the plane S2 to (with ), then the total matrix from S1 to S3 is
Using Eq. (30) and the group multiplication rule of , we can further put Eq. (351) into another form
where and we have made a reasonable assumption that [hence ], 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 .
In order to examine the second diffraction domain from to (not ), we need to translate the output field amplitude observed on the plane to the field amplitude observed on sphere plane i.e., putting into [see Eq. (347)] by taking the phase compensation [see Eq. (342)] into account. Thus the field transform from S to is []
Comparing Eq. (354) with Eq. (182), we have
and noticing that yields we have
combine Eqs. (340) and (356), and it then follows that (let )
and
where we have used Eqs. (346), (348) and (352).
Substitution of Eqs. (347), (355), (357) and (358) into Eq. (354) yields
The first equation of Eq. (359) indicates that it is just a CFrFT of from S to . Comparing Eq. (359) with Eq. (349), we see
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 tends to infinity,
A more general requirement for a mother wavelet is for demanded to have vanishing moments (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 (μ is scaling parameter, s is a translation parameter, ), and the dilated-translated function is defined as
while the wavelet integral transform of a signal function by ψ is defined by
We can express Eq. (363) as
where is the state vector corresponding to the given mother wavelet, is the state to be transformed, and
is the squeezing-translating operator [
116,
118,
119],
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) (
)
where , and the operator identity is used. In particular, when , it is reduced to the well-known squeezing operator,
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
where is the momentum eigenstate, we can recast the condition into quantum mechanics as
which indicates that the probability of a measurement of by the projection operator with value being zero. Withouta loss of generality, we suppose
where are chosen as to let obey the condition (361). Using the coherent states’ overcompleteness relation, we have
Eq. (371) provides a general formalism in finding the qualified wavelets. For example, assuming for , the coefficients of the survived terms should satisfy
and becomes
Projecting it onto the coordinate representation, we get the qualified wavelets
where we have used , and is the Hermite polynomials. Now we take some examples.
Case 1: In Eq. (371), by taking (otherwise), we have
It then follows that
which is just the Maxican hat wavelet, satisfying the condition Hence 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 the matrix element is just the wavelet transform of with respect to
Case 2: When , , and , from Eq. (374) we obtain (see Fig. 4)
which also satisfies . Note that when , , and , we obtain a slightly different wavelet (see Fig. 5). Therefore, as long as the parameters conforms to condition (372), we can adjust their values to control the shape of the wavelets.
Case 3: When , , and , we get (see Fig. 6)
and . 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 and that with our new wavelet . Concretely, we map a simple cosine signal by performing the wavelet transforms with , into a 2D space , where X denotes the location of a wavelet and a its size. The resulting wavelet transforms by (=) and are
where and are the normalization factors for and respectively. The wavelet integral are also called wavelet coefficient which measures the variation of cos in a neighborhood of whose size is proportional to . The contour line representation of and 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
has one maximum, while
has one main maximum and one subsidiary maximum (“two islands”), so when
scales its size people have one more chance to identify the frequency information of the cosine wave than using
. Interestingly enough, the “two islands” of
in Fig. 8 can be imagined as if they were produced while the figure of
deforms into two sub-structures along
a-axis. (ii) Near the maximum of
the density of the contour lines along
a-axis is higher than that of
, which indicates that the new wavelet
is more sensitive in detecting frequency information of the signal at this point. Therefore,
may be superior to
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]:
where is a mother wavelet whose Fourier transform is . In the context of quantum mechanics, according to Eq. (364), we see that the quantum mechanical version of Parseval theorem should be
and since involved in is , and is the momentum eigenvector
Eq. (382) indicates that once the state vector corresponding to mother wavelet is known, for any two states and , their overlap up to the factor [determined by Eq. (396)] is just their corresponding overlap of WTs in the ( parametric space.
Proof of Equation (382): In order to show Eq. (382), we calculate
which leads to
where we have used the formula
Inserting the completeness relation into the left side of Eq. (381) and then using Eq. (385), we have
where
and
Further, we can put Eqs. (388) and (389) into the following forms:
and
where
Thus, when the definite integration satisfies the admissible condition, i.e.,
which leads to
Eq. (387) can be transformed to
where
Thus, the theorem is proved. Especially, when Eq. (382) becomes
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 in Eq. (382), then using Eq. (364), we see that Eq. (382) is reduced to
Due to Eq. (365) we have
It then follows that
which means
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
,
in Eq. (381) one can see that
which is a new orthogonal property of mother wavelet in parameter space spanned by . In a similar way, we take a number state, since then we have
Or take and is the coherent state, then
This indicates that is -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
where
Here is the single-variable Hermit polynomial. Substituting Eq. (405) into Eq. (396), we have
Thus, for the Mexican hat wavelet (376), we see
Eq. (408) can be checked as follows. Using Eq. (376) and noticing that , we can put the left-hand side of Eq. (408) into
where we have used the integration formulas
and
Next, we examine if the Morlet wavelet obey the formalism (402). The Morlet wavelet is defined as [
125–
127]
Substituting Eq. (412) into the left hand side of (402), we get
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
[
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
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,
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 i.e.,
where
and is the Husimi distribution function as well as is the Husimi operator,
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 is the vacuum state , we see that
From Eq. (366) it then follows that
Substituting Eq. (416) and into Eq. (419), we get
Then the WT of Eq. (418) can be further expressed as
Using the normally ordered form of the vacuum state projector and the IWOP method as well as Eq. (420), we have
Now we explain why
is the Husimi operator. Use the formula for converting an operator
into its Weyl ordering form [
129]
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
This is the Weyl ordering form of
Then according to Weyl quantization scheme [
35] we know the classical corresponding function of a Weyl ordered operator is obtained by just replacing
and in this case the Weyl rule is expressed as
where at the last step we used the Weyl ordering form of the Wigner operator
[
130]
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.,
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
as a weighting function to synthesize the original complex signal ,
This is named the symplectic-transformed — translated wavelet transform. One can see that the mother wavelet ψ generates the other wavelets of the family through a translating transform followed by a symplectic transform. ( are the symplectic transform parameter, is a translation parameter, s, r and ). This will be seen more clearly if we write the second transform in a matrix form
where M is a symplectic matrix satisfying , .
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 is the complex Fourier exponentials, (Note that is pure imaginary)
where is the complex Fourier transform of Then we form the adjoint operation
from which we have
so we get the inversion formula
Eq. (435) leads us to impose the normalization
in order to get the wavelet representation
Then we can have a form of Parseval’s theorem for this new wavelet transform:
Preposition: For any f and we have
Proof: Let us assume and are the complex Fourier transform of and respectively,
Recall the convolution theorem defined on complex Fourier transform,
Hence from Eqs. (440) and (429), as well as Eq. (430) we see that can be considered as a convolution in the form
It then follows from Eq. (440) that
Therefore, using Eq. (436), we see that the further integration yields
which completes the proof.
Theorem: From the Proposition (438) we have
that is, there exists an inversion formula for arbitrary function . In fact, in Eq. (430) when we take then
Substituting Eq. (445) into Eq. (443) we obtain Eq. (444).
14.1.2 Relation between and optical Fresnel transform
Now we explain why the idea of is originated from the optical Fresnel transform. We can visualize the symplectic-transformed — translated wavelet transform in the context of quantum mechanics. Let , is the coherent state, is the vacuum state in Fock space, then Eq. (429) can be expressed as
where is defined as
and Employing the normal ordering of the vacuum projector and the IWOP technique to perform the integration in Eq. (447), one can get the explicit form of . The transformation matrix element of in the coordinate representation is just the kernel of optical diffraction integration (86) (Fresnel transform), which explains our motivation to introduce . In particular, when reduces to the well-known single-mode squeezing operator , 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
in a two-mode coherent state
representation onto quantum operator
, we are naturally led to develop the SWT in Eq. (430) to the so-called entangled SWT (ESWT) [
133] for signals
defined in two complex planes,
where
is used as a weighting function to synthesize the signal regarding two complex planes. One can see that the mother wavelet φ generates the family 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 given by Eq. (446). The new symplectic transform can be seen more clearly if we write it in a matrix form:
where is symplectic satisfying is the 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 is the complex Fourier exponential,
According to Eqs. (448) and (449) we evaluate its ESWT
Making the integration variables transform Eq. (452) becomes
The last integration is just the complex Fourier transform (CFT) of . Denoting it as we have
Then we form the adjoint operation of Eq. (454),
where the integration in the last line is just the CFT of ϕ [comparing with Eq. (453)], thus Eq. (455) leads to
From Eq. (456) we have
which together with Eq. (451) lead to
Eq. (458) implies that we should impose the normalization
such that the reproducing process exists
Eq. (459) may be named the generalized admissibility condition. Now we have the corresponding Parseval theorem: For any g and , we have
Proof: Assuming and to be CFT of and , respectively,
Recall the corresponding convolution theorem
and from Eqs. (448), (462) and (463), we see that can be considered as a convolution in the form [Noting that the CFT of is see Eqs. (453) and (454)]:
Using Eq. (464) we calculate
As a consequence of Eqs. (459) and (465) the further integration yields
which completes the proof.
Inversion formula: From Eq. (461) we have
That is, there exists an inversion formula for which represents the original signal as a superposition of wavelet functions with the value of entangled wavelet transform serving as coefficients. In fact, in Eq. (448) when we take
then
Substituting Eqs. (468) and (469) into Eq. (466), we obtain Eq. (467). We can visualize the ESWT in the context of quantum mechanics. Let and use Eqs. (449) and (450), Eq. (448) is expressed as
where is defined as
When and 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)
,
which is the common eigenvector of the operator and i.e., and where is the coordinate operator, (, and is complete,
and exhibits partly non-orthogonal property (for and orthonormal property (for x),
Therefore,
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
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 ():
where The family of mother wavelet ψ involves both the the symplectic transform of α and the dilation-transform of x,
Let then Eq. (475) can be expressed as the quantum mechanical version
where is defined as
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 is the exponential
Making the integration variables transform which leads to and , Eq. (479) becomes
where is just the Fourier transform of
Then we perform the adjoint WT of Eq. (475), using Eqs. (476) and (480) we see
From Eq. (482) we obtain
which leads to
Eq. (484) implies that we should impose the normalization
such that the reproducing processexists
Eq. (485) may be named the generalized admissibility condition. Now we have the corresponding Parseval theorem: For any g and , we have
Proof: Assuming and to be the Fourier transforms of and , respectively,
In order to prove Eq. (487), we first calculate . In a similar way of deriving Eq. (175), using Eqs. (170), (171) and the inversion formula of Eq. (488), we have
It then follows that
Substituting Eq. (490) into the left-hand side (L.H.S) of Eq. (487) and using Eq. (485), we see
Thus we complete the proof of Eq. (487).
Inversion Formula: From Eq. (487) we have
which is the inversion formula for the original signal expressed by a superposition of wavelet functions with the value of continuous WT serving as coefficients. In fact, in Eq. (475), when we take then
Substituting Eq. (493) into Eq. (491), we get (492).
We can visualize the new WT
in the context of quantum optics. Noticing that the generalized squeezing operator
in Eq. (478) is an image of the combined mapping of the classical real dilation transform
(
) and the classical complex symplectic transform
in
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
in the entangled state representation
is
In Fock space is a two-mode EPR entangled state in Eq. (28).
Then using Eqs. (472) and (28), we obtain
Substituting Eq. (495) into Eq. (494) and using Eq. (79), we obtain
which is just the kernel of a mixed lens – Fresnel transform, i.e., the variable of the object experiences a generalized Fresnel transform, while undergoes a lens transformation. Thus, based on , 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
quantum Fourier transform, which can be expressed as [
137,
138]
Now the continuous Hadamard transform, which used to go from the coordinate basis
to the momentum basis, is defined as [
139]
where σ is the scale length, and is named Hadamard operator. Using the completeness of we have
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
, and enlightened by Eqs. (81) and (498), we now construct the following ket–bra integration [
140]
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 as well as the IWOP technique, we get
where
and
They are all within the normal ordering symbol
. Now perform the integration over
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]
15.2 The properties of Hadamard–Fresnel operator
Note and can be considered a mode independent of another mode thus we have the operator identity
Using Eq. (503), we can rewrite Eq. (502) as
where
and
while U2 is the Fresnel operator for mode and U1 is named the Hadamard operator for mode
It then follows that
from which we see that the Hadamard–Fresnel combinatorial operator can play the role of Fresnel transformation for Physically, and can be two output fields of a beamsplitter.
In a similar way, we have
which for the quadrature leads to
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
The mutual transform in Eq. (509) can be realized by
while the two-mode squeezing operator is therefore
From Eq. (504) and Eq. (512), we see that the Hadamard–Fresnel combinatorial operator can be decomposed as
It can also be seen that U is unitary, and .
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 and Hadamard transformation for mode 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 a Fresnel operator U2, a two-mode squeezing operator and the total momentum-average position exchanging operator. Physically, and 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.
© Higher Education Press and Springer-Verlag Berlin Heidelberg 2012