Path Integrals: From Quantum Mechanics to Photonics Charles W. Robson,1 Yaraslau Tamashevich,1 Tapio T. Rantala,2 and Marco Ornigotti1, a) 1)Photonics Laboratory, Tampere University, P.O. Box 692, FI-33014, Tampere, Finland 2)Department of Physics, Tampere University, P.O. Box 692, FI-33014, Tampere, Finland (Dated: 25 June 2021) The path integral formulation of quantum mechanics, i.e., the idea that the evolution of a quantum system is determined as a sum over all the possible trajectories that would take the system from the initial to its final state of its dynamical evolution, is perhaps the most elegant and universal framework developed in theoretical physics, second only to the Standard Model of particle physics. In this tutorial, we retrace the steps that led to the creation of such a remarkable framework, discuss its foundations, and present some of the classical examples of problems that can be solved using the path integral formalism, as a way to introduce the readers to the topic, and help them get familiar with the formalism. Then, we focus our attention on the use of path integrals in optics and photonics, and discuss in detail how they have been used in the past to approach several problems, ranging from the propagation of light in inhomogeneous media, to parametric amplification, and quantum nonlinear optics in arbitrary media. To complement this, we also briefly present the Path Integral Monte Carlo (PIMC) method, as a valuable computational resource for condensed matter physics, and discuss its potential applications and advantages if used in photonics. I. INTRODUCTION proposing the Lagrangian as a more natural basis for a theory of quantum mechanics, as opposed to a Hamiltonian-based The path integral formulation of quantum mechanics, de- method which he argued was less fundamental due to its non- 45 veloped in the mid-twentieth century, is not only a remark- relativistic form . This led him to propose the key idea that able synthesis of several of the core ideas of theoretical a quantum mechanical transition amplitude for a particle is physics, but it is also a powerful computational technique for given by a phase factor controlled by the action along that 46 the analysis of a huge variety of physical systems in very particle’s path . This also led him to state, that the classical different contexts, such as quantum mechanics1–6, quantum path of a quantum system can be interpreted as resulting from field theory7–10, gauge field theory11,12, black hole physics6, the constructive interference of all such paths. In other terms, quantum gravity13,14, string theory15, topology6,16condensed the action S of a given system counts, de-facto, the number of matter physics17–23, statistical mechanics2,6,24,25, polymer waves of the path in units of Planck’s constant h, and, there- 47 physics5,6, financial markets5, optical communications26, fore, for each path, the phase of the related wave is given atomic physics27,28, spectroscopy29, light propagation in by exp(i2p S=h) = exp(iS=h¯). This allows the evaluation of turbid media30 and classical31–36 and quantum37–39 optics, the interference pattern of the particle dynamics, and essen- amongst others. tially represents what is nowadays commonly understood as 2,5 The underlying concepts of the path integral approach are path integral . sometimes considered difficult to grasp. Indeed, from a philo- The pivotal step in the development of the theory occurred sophical standpoint, its underpinnings are extraordinary – in in the 1940s, when Feynman formulated his version of quan- describing a system, all possible paths between its initial and tum mechanics, a ‘third way’ following the earlier, well- 4 final state must be taken into account mathematically. This known Schrödinger and Heisenberg alternatives . Based on idea of summing over all paths has previously been charac- integration over all paths between initial and final physical terised ontologically as everything that can happen does hap- states, with each path contributing an action-dependent phase pen40. as Dirac had proposed, Feynman’s formulation culminated As befits the subject, the history of path integrals is rather in his important 1948 paper entitled Space-time approach to 1 circuitous and does not follow a straight line from Feyn- non-relativistic quantum mechanics . man’s seminal work on the subject, published in 19481, to From a conceptual point of view, one of the most interest- arXiv:2105.00948v2 [quant-ph] 24 Jun 2021 the present day. Already in the 1920s, the mathematician ing consequences of path integrals is that they can provide a Norbert Wiener had developed a method to treat Brownian deep understanding of the relation between quantum and clas- motion and diffusion using a technique of integrating over sical mechanics, as the limit h ! 0 emerges naturally from the 3 paths, with formal similarities to Feynman’s eventual con- formalism as the classical limit of the theory . struction but in a purely classical context41–43. During the Although the concepts underlying the path integral formal- early part of the following decade, the general ideas of de ism may at first appear quite alien, they are deeply profound. Broglie and Schrödinger, that waves can be associated with Even if its utility were limited, the beauty of the idea would particle dynamics44, motivated Dirac to publish a paper in the still merit a wide audience. The fact is, however, that the path Physikalische Zeitschrift der Sowjetunion (Physical Journal of integral has immense value as a practical computational tech- the Soviet Union) setting the stage for future developments by nique in the physical sciences, and it can be applied to solve problems in many diverse areas, even beyond physics. Two examples showing the remarkable breadth of its applicability are its use in quantum gravity, where the “sum over histories” a)Electronic mail: marco.ornigotti@tuni.fi is interpreted as a sum over all different spacetime configu- 2 rations interpolating between an initial and final state of the density matrix, which in space basis means finite closed loop universe13, and in financial market modelling, where the for- paths. Thus, in terms of path integrals each of the particles malism has proven useful, since the time dependence of asset propagate from a position in space in imaginary time back prices can be represented by fluctuating paths5. to the same position, the time period being inversely propor- The universality of the technique has allowed scientists tional to the temperature. Then, with Metropolis Monte Carlo to tackle many problems, and gain tremendous physical in- it is possible to sample particle paths with correct weight in sight into them. In statistical physics, for example, path in- predefined temperature and collect data enough for conver- tegrals conveyed the basic framework for the first formula- gence of expectation values of relevant operators. tion of the renormalisation group transformation, and they are This approach based on imaginary time propagation is largely employed to study systems with random distribution called Path Integral Monte Carlo method (PIMC). David of impurities6. In particle physics, they allow one to under- Ceperley and coworkers have carried out seminal develop- stand and properly account for the presence of instantons48. ment work and a sizable number of PIMC simulations of In quantum field theory they provide the natural framework to various many-particle systems ranging from superfluid He17 quantise gauge fields11,12,49. In chemical, atomic, and nuclear and neutron matter59 to electrons and hydrogen in extreme physics, on the other hand, they have been applied to vari- conditions60, including both bosonic and fermonic particles. ous semiclassical schemes for scattering theory6. Through In recent years, one of the authors of this tutorial (TTR) path integrals, topological and geometrical features of clas- has significantly contributed to taking PIMC simulations to sical and quantum fields can be readily investigated and be new though simple quantum particle systems, such as small used to create novel forms of perturbative and nonperturba- atoms61,62, molecules63–66, a chemical reaction67 and quan- tive analysis of fundamental processes of Nature15,50,51. In tum dots68–70, with the ultimate goal of providing a more addition to that, path integrals allow one to re-interpret es- accurate description of their electronic structure and related tablished results, such as, for example, the BCS theory of properties including many-body effects, and how they change superconductivity7,52, from a novel, more insightful, perspec- with temperature62,71,72. In this context, PIMC has proven to tive. be a very reliable and excellent method to calculate the elec- Explicit, analytical solutions to problems formulated in tric polarisabilities of atomic and molecular systems, leading terms of path integrals, however, are scarce, and only available therefore to an accurate estimation of the optical properties for very simple systems, such as a free particle, or the ubiq- of both individual small quantum systems and collections of uitous harmonic oscillator2. The complexity of the path inte- them, in the form of dilute gases66,73,74. gral formalism, in fact, increases very rapidly to overwhelm- A different approach, based on a real, rather than imagi- ing levels of difficulties for many simple problems. As an nary, time path integral (RTPI) has been recently proposed as example of that, the simplest quantum system, i.e., a single a way to describe the full quantum dynamics of a quantum hydrogen atom, required nearly 40 years to be fully solved in system at zero-Kelvin, and to also characterise the evolution terms of path integrals53. of its eigenstates69,73,75–78. A combination of PIMC and RTPI On the other hand, it is amidst complex and computa- therefore gives the possibility to have a comprehensive tool tionally challenging problems that path integrals show their to study the properties of complex systems and their classical true potential, providing a simple, insightful and intuitive and quantum evolution.