Submitted:

26 October 2023

Posted:

30 October 2023

You are already at the latest version

A peer-reviewed article of this preprint also exists.

Abstract
By utilizing the Fourier transform, we present a practical method for evaluating a function of a derivative applied to any other arbitrary function of a single complex variable. As an illustration of this approach, we compute the actions of the displacement and squeeze operators on an arbitrary function, as well as the propagation of paraxial fields, Airy and number states, without the need for algebraic quantum operator techniques or the Fresnel integral.
Keywords: 
;  

Introduction

Fourier optics is a theoretical and computational framework used to simulate and understand the propagation of light in free space and optical elements [1]. Moreover, quantum optics and its mathematical methods can explain the connections between geometrical and wave optics, classical mechanics and quantum mechanics. In addition, it offers new connections and operating procedures for solving problems in classical optics. These ideas were based on coherent and squeezed states [2], systematically including group theory associated with displacement and squeeze operators [3,4].
The main purpose of this work is to address the following question: Is it possible to calculate the actions of quantum operators on arbitrary functions using Fourier transforms? As illustrative examples, we will determine the actions of the displacement and squeeze operators on arbitrary functions and compute the value of a freely propagating paraxial field using the quantum propagator for a free particle acting on Airy functions and number states. These propagations will be carried out through the use of Fourier transforms.

Fourier transform differentiator operator

The standard 1D Fourier transform and inverse Fourier of an integrable function f ( x ) are defined as follows:
F f ( x ) = F ( k ) = d x e i k x f ( x ) ,
F 1 F ( k ) = f ( x ) = 1 2 π d k e i k x F ( k ) ,
respectively. Indeed, if the Fourier transform of f ( x ) exists, it is not difficult to establish the differentiation property
F d n d x n f ( x ) = ( i k ) n F ( k ) .
If we assume that g ( x ) is an analytic function, which admits a power-series expansion, and we write g ( x ) = n = 0 A n x n , then the extension to any function of the derivative is given by
F g d d x f ( x ) = g ( i k ) F ( k ) .
Furthermore, assuming that inverse Fourier transform exists, we can obtain
g d d x f ( x ) = F 1 g ( i k ) F ( k ) .
We therefore call this function a Fourier transform differentiator operator.

Displacement and squeeze operators

Coherent and squeezed states are fundamental components of the theoretical framework of modern optics. One way to define coherent states is by using the displacement operator D ^ ( α ) = exp ( α a ^ α * a ^ ) , where α = ( x 0 + i p 0 ) / 2 represents a point in the optical phase space, and a ^ = ( x ^ i p ^ ) / 2 and a ^ = ( x ^ + i p ^ ) / 2 are the creation and annihilation operators. In the position representation, we have
x | D ^ ( α ) | f = e i p 0 x exp x 0 d d x f ( x ) .
By comparing equation (5) with (6), we obtain that g ( x ) = exp ( x 0 x ) . Using the relation (5), we can express the action of the displacement operator on an arbitrary function f ( x ) as:
exp x 0 d d x f ( x ) = 1 2 π d k e i k ( x x 0 ) F ( k ) = f ( x x 0 ) .
Another class of minimum-uncertainty states are the squeezed states, which can be obtained using the squeeze operator S ^ ( r ) = exp 1 2 r ( a ^ 2 a ^ 2 ) , where r is an arbitrary real number. In the position representation we have,
x | S ^ ( r ) | f = e r 2 exp r x d d x f ( x ) .
By changing x = exp ( u ) , we obtain
x | S ^ ( r ) | f = e r 2 exp r d d u f ( e u ) = e r 2 f ( e r x ) .
This expression was stated without proof in equation (4.106) of the reference [5].

Freely propagating paraxial field

Airy function

The Airy function is defined by the inverse Fourier transform as
Ai ( x ) = 1 2 π d k exp i 3 k 3 + i k x .
If we consider the Airy function as an initial state for the paraxial equation, the free coordinate evolution reads as follows:
Ai ( x , z ) = exp i κ z 2 p ^ 2 Ai ( x , 0 ) , = exp i κ z 2 d 2 d x 2 Ai ( x , 0 ) ,
where κ = λ / 2 π , and z represents the propagation coordinate, (for time evolution, we set κ = and z = t ). By comparing equation (5) with (11), we obtain that g ( x ) = exp [ i ( κ z / 2 ) x 2 ] . Coordinate free evolution of Airy function, reads
Ai ( x , z ) = 1 2 π d k exp i 3 k 3 i κ z 2 k 2 + i k x .
After rearranging terms and using the definition of the Airy function, we have
Ai ( x , z ) = Ai x κ z 2 2 exp i κ z 2 x κ 2 z 2 6 ,
which is the same simple analytical expression reported in reference [6].

Number states

The free evolution of the field E ( x , z ) can be expressed in terms of the eigenfunctions of the quantum harmonic oscillator ψ n ( x ) as follows:
E ( x , z ) = n = 0 E n exp i κ z 2 d 2 d x 2 ψ n ( x ) ,
where
ψ n ( x ) = N n exp 1 2 x 2 H n ( x ) .
with N n = 1 / 2 n n ! π . The amplitudes E m are determined by the initial field, such that
E m = d x ψ m * ( x ) E ( x , 0 ) .
By comparing equation (14) with (5), we can deduce that
exp i κ z 2 d 2 d x 2 ψ n ( x ) = F 1 g ( i k ) F ( k ) = P n ( x , z ) ,
where g ( x ) = exp [ i ( κ z / 2 ) x 2 ] and f ( x ) = ψ n ( x ) . We have evaluated the direct and inverse Fourier transforms using Cauchy’s integral formula, considering a suitable contour containing a real axis, (see Gradshteyn and Ryzhik, Eq. 7.375.8) [7], we write
F ( k ) = ( i ) n 2 π N n exp 1 2 k 2 H n ( k ) ,
and therefore,
P n ( x , z ) = N n 1 + i κ z 1 i κ z 1 + i κ z n 2 exp 1 1 + i κ z x 2 2 H n x 1 + κ 2 z 2 ,
and we have
E ( x , z ) = n = 0 E n P n ( x , z ) ,
which is the same simple analytical expression for free propagation of any initial field, based on the eigenstates of the quantum harmonic oscillator, as that obtained by operators algebra in reference [8].

Conclusions

We have presented a method for evaluating a one-derivative function applied to any other arbitrary function of a single complex variable. This approach is based only on Fourier transforms and does not revolve around an algebra of creation and annihilation operators.

References

  1. J. Goodman, Introduction to Fourier Optics, McGraw-Hill physical and quantum electronics series, W. H. Freeman, 2005.
  2. C. C. Gerry and P. L. Knight, Introductory Quantum Optics, Cambridge University Press, 2004.
  3. G. S. Agarwal, Quantum optics, Cambridge University Press, New York, 2013.
  4. H. Moya-Cessa, Physics Reports 432 (1) (2006) 1–41.
  5. S. Dong, Factorization Method in Quantum Mechanics, Fundamental Theories of Physics, Springer Netherlands, 2007.
  6. M. V. Berry, N. L. M. V. Berry, N. L. Balazs, Nonspreading wave packets, American Journal of Physics 47 (3) (1979) 264–267.
  7. D. Zwillinger, A. D. Zwillinger, A. Jeffrey, Table of Integrals, Series, and Products, Elsevier Science, 2007.
  8. J. A. Anaya-Contreras, A. J. A. Anaya-Contreras, A. Zúñiga-Segundo, H. M. Moya-Cessa, J. Opt. Soc. Am. A 38 (5) (2021) 711–718.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Downloads

86

Views

37

Comments

0

Subscription

Notify me about updates to this article or when a peer-reviewed version is published.

Email

Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2025 MDPI (Basel, Switzerland) unless otherwise stated