Arxiv:2002.00255V3 [Quant-Ph] 13 Feb 2021 Lem (BVP)
Total Page:16
File Type:pdf, Size:1020Kb
A Path Integral approach to Quantum Fluid Dynamics Sagnik Ghosh Indian Institute of Science Education and Research, Pune-411008, India Swapan K Ghosh UM-DAE Centre for Excellence in Basic Sciences, University of Mumbai, Kalina, Santacruz, Mumbai-400098, India ∗ (Dated: February 16, 2021) In this work we develop an alternative approach for solution of Quantum Trajectories using the Path Integral method. The state-of-the-art technique in the field is to solve a set of non-linear, coupled partial differential equations (PDEs) simultaneously. We opt for a fundamentally different route. We first derive a general closed form expression for the Path Integral propagator valid for any general potential as a functional of the corresponding classical path. The method is exact and is applicable in many dimensions as well as multi-particle cases. This, then, is used to compute the Quantum Potential (QP), which, in turn, can generate the Quantum Trajectories. For cases, where closed form solution is not possible, the problem is formally boiled down to solving the classical path as a boundary value problem. The work formally bridges the Path Integral approach with Quantum Fluid Dynamics. As a model application to illustrate the method, we work out a toy model viz. the double-well potential, where the boundary value problem for the classical path has been computed perturbatively, but the Quantum part is left exact. Using this we delve into seeking insight in one of the long standing debates with regard to Quantum Tunneling. Keywords: Path Integral, Quantum Fluid Dynamics, Analytical Solution, Quantum Potential, Quantum Tunneling, Quantum Trajectories Submitted to: J. Phys. A: Math. Theor. I. INTRODUCTION Our derivation sets off from the Feynman theory of Path Integrals.[8, 9].[10] Firstly, we develop a general solution for Quantum Fluid Dynamics (QFD) [1–5] is a tool for the propagator. We derive each term in the series as a closed the study of a quantum system as motion of the proba- form functional of the classical path, which then can be used bility fluid. It is widely employed in theoretical chem- to obtain the Quantum Trajectories. The method is similar to istry, condensed matter physics and material sciences. The Ursell decomposition of statistical mechanical partition func- main advantage of this approach starts with recasting of tions. the Schrodinger’s equation into two coupled differential The paper is structured as follows. We first discuss in short, equations involving real quantities, formally known as the the results of Path Integrals and QFD, relevant to our work. Madelung transformation[2]. This in particular often make Then we present our derivation, firstly for 1-D single particle it computationally more useful over the usual formulation in case, and then generalize it to higher dimensions and many- the study of dynamics of various molecular or quantum optics body systems. Then we present a toy example of an 1-D dou- systems. ble well potential, solved analytically using the ideas we de- The usual route for computation in QFD is in some or the velop for the purpose of illustration of the method. The clas- other way related to solving a set of coupled differential equa- sical boundary value problem, in this case, has been obtained tions simultaneously[6]. In this study, we seek to develop a perturbatively using Linsted-Poincare method. We show that fundamentally different route to achieve the quantum trajec- the solution obtained by our method exhibits tunneling even in tories by invoking the Path Integral method, which reduces the very first order of perturbation of the classical path, when the problem of finding the quantum trajectories to solving the the quantum part of the calculation is left exact. We conclude corresponding classical trajectory as a boundary value prob- indicating some applications and further directions of study. arXiv:2002.00255v3 [quant-ph] 13 Feb 2021 lem (BVP). We show that, the QP can be expressed in terms of a series, each term of which can be computed in a closed cl form as a functional of the classical path, x (t). These ana- II. FORMULTION: A GENERAL SOLUTION FOR PATH lytical functional forms can be provided directly as an input INTEGRAL PROPAGATOR to the machine simplifying the algorithm for computation of Quantum Trajectories to computation of a BVP. This method, In this section we present an exact closed form series solu- hence, is expected to be computationally superior to the Quan- tion for the Kernel of the Path Integral (the two-point corre- tum Trajectory Method (QTM)[7] due to Wyatt et. al, which lator) for a quantum particle for any general applied potential sets the current standard for time complexity. V(x). In Path Integral formalism, the Kernel can be solved exactly for some simple systems viz the free particle, Simple Harmonic Oscillator and approximately for a few more com- plicated systems. To the best of our knowledge, the solution ∗ [email protected] presented here is the first exact solution without assuming any 2 special property of the applied potential other than continuity In this work we look for footing of the QP in the Feynman and differentiabilty. We express the Kernel as a functional of Path Integrals, with the aim of developing formal functional the classical path. expression of the former in terms of the classical path and the The correlator gives the amplitude at a point xt , after a time Cauchy data. We first treat the one dimensional system and of flight t, given that it has started from an initial point x0. It then extend it to higher dimensions. is also the Green’s function to the Schrodinger equation [9], Let’s denote the classical path of the system, from a initial cl rendering it the formal equivalence with usual formulations. position x0, to a final point xt after a time of flight t, by x (t). Equipped with this, we pose the solution to quantum trajecto- Being solution to a second order Boundary Value Problem, it ries, in an integral equation framework. parametrically depends on, x0 and xt . According to Path Integral approach, the dynamics of a A general path of the system, denoted by x(t), is any con- quantum particle from one point to another, is contributed by tinuous function of t that is constrained by, x(0) = x0 and all the possible paths that connects them. The contribution x(t) = xt . Any such general path can be broken down as, from each path is weighed by exp(iS[x(t)]=h¯), where S[x(t)] is the corresponding Action or Hamilton’s principal function. x(t) = xcl(t) + y(t): (5) The net propagator K(xt ;t;x0;0) is obtained by adding the contributions from all paths. The integral varying y(t) to span the space of continuous functions on [t0;t] subject to the boundary conditions, y(t0) = y(t) = 0 tan- Z ¥ tamounts to span the whole space of allowed paths of Path y(xt ;t) = K(xt ;t;x0;0)y0(x0)dx0 (1) Integral formulation, over which we need to sum over to ob- −¥ tain the propagator.[13] One way to compute the Path Integral is to first discretise gives the time evolution of the wave-function. the time axis into n moments, separated by e units. Next, we exploit the Madelung transformation[2], to ex- press y(t) in polar form. [11]. The Madelung Transformation provides a way of looking at the quantum systems as proba- t j = je; bility fluids, by decoupling the Schrodinger equation into the t = 0 ; t = t following coupled equations. [12] 0 n For a particle of mass M in a potential V(x) we have, The Path Integral can then be recovered in the limit n ! 0;e ! ¥, such that the product remains constant, ne = t, after the n-dimensional integral is carried out. Equation (5) under ¶r ∇(rv) + = 0 (2) the discretisation becomes, ¶t cl x(t j) = x j = x j + y j: (6) ¶S 1 − = Mv2 +V(x) + Q(x;t) (3) ¶t 2 at each of these moments. By definition y0 = yn = 0. The propagator is defined as, Here r(x;t) = R2(x;t), v(x;t) = ∇ S(x;t)=M and the quan- tum potential, Q(x,t) is given by, Z Z t i 0 0 0 0 K(xt ;t;x0;0) = exp L (x(t );x˙(t );t ) dt Dx (7) h¯ 0 2 2 h¯ ∇ R(x;t) 0 Q(x;t) = − (4) where L (x;x˙;t ) is the classical Lagrangian of the system. 2M R(x;t) In the discretised version it is rephrased as, Z ¥ Z ¥ n 1 ie h M 2 i lim ··· exp (x j − x j−1) −V(x j) dx1 ···dxn (8) n!¥ ∑ 2 A −¥ −¥ h¯ j=1 2e where A is the normalization constant which, in general, depends on n. Throughout the paper this limit is understood ¥ ym m to be taken as n ! ¥ and e ! 0 simultaneously, with ne = t cl cl j ¶ V(x j) = V(x j + y j) = V(x j ) + V(x)j cl (9) ∑ m x=x j remaining constant. m=1 m! ¶x The potential over the general path can be Taylor expanded the around the classical position at each moment, as [14] Substituting (6),(9) into (8) we obtain the propagator as, 3 Z ¥ Z ¥ n 1 ie h M cl cl 2 M cl cl M 2 lim ··· exp (x j − x j− ) + (x j − x j− )(y j − y j−1) + (y j − y j−1) n!¥ ∑ 2 1 2 1 2 A −¥ −¥ h¯ j=1 2e 2e 2e ¥ m m cl y j ¶ i −V(x j ) − V(x)j cl dy1 ···dyn (10) ∑ m x=x j m=1 m! ¶x The terms in (10), that are independent of y j, can be taken corresponding case of the Harmonic Oscillator, that forms cl out of the integrals.