Path Integral Description Research Articles

Spectral density of the Wigner path integral operator correlation function representation. Monte Carlo simulation of the fermion dynamic structure factor

An important mathematical function that describes the interparticle correlations and their time evolution is the dynamic structure factor (DSF). According to Van Hove the dynamic structure factor (DSF) is the space and time Fourier transform of the density-density correlation function of the unperturbed system. In the case of classical systems, computations of the DSF typically employ classical molecular dynamics techniques. However, the present study focuses on quantum systems, in particular on the quantum dynamic structure factor. The main idea of this article is to employ the Wigner formulation of quantum mechanics to derive a new path integral description of the quantum DSF. The next important idea is to consider paths in this representation as stationary random processes and apply the Wiener–Khinchin theorem to express the DSF through the expected values of the power spectrum of trajectories in this representation. To calculate the DSF, spin-resolved radial distribution functions (RDFs) and other thermodynamic parameters of a 3D system of soft sphere fermions we developed a Wigner Path Integral Monte Carlo (WPIMC) method. The observed peaks on the RDFs and DSFs were attributed to wave interference resulting from a multiple-scattering and perturbation effects due to particle exchange and interaction. This phenomenon, observed in a soft sphere scattering system, may serve as a precursor to Anderson localisation, a mechanism in which wave interference across multiple scattering channels prevents wave propagation.

The magnetic maze: a system with tunable scale invariance

Random magnetic field configurations are ubiquitous in nature. Such fields lead to a variety of dynamical phenomena, including localization and glassy physics in some condensed matter systems and novel transport processes in astrophysical systems. Here we consider the physics of a charged quantum particle moving in a “magnetic maze”: a high-dimensional space filled with a randomly chosen vector potential and a corresponding magnetic field. We derive a path integral description of the model by introducing appropriate collective variables and integrating out the random vector potential, and we solve for the dynamics in the limit of large dimensionality. We derive and analyze the equations of motion for Euclidean and real-time dynamics, and we calculate out-of-time-order correlators. We show that a special choice of vector potential correlations gives rise, in the low temperature limit, to a novel scale-invariant quantum theory with a tunable dynamical exponent. Moreover, we show that the theory is chaotic with a tunable chaos exponent which approaches the chaos bound at low temperature and strong coupling.

The open effective field theory of inflation

In our quest to understand the generation of cosmological perturbations, we face two serious obstacles: we do not have direct information about the environment experienced by primordial perturbations during inflation, and our observables are practically limited to correlators of massless fields, heavier fields and derivatives decaying exponentially in the number of e-foldings. The flexible and general framework of open systems has been developed precisely to face similar challenges. Building on previous work, we develop a Schwinger-Keldysh path integral description for an open effective field theory of inflation, describing the possibly dissipative and non-unitary evolution of the Goldstone boson of time translations interacting with an unspecified environment, under the key assumption of locality in space and time. Working in the decoupling limit, we study the linear and interacting theory in de Sitter and derive predictions for the power spectrum and bispectrum that depend on a finite number of effective couplings organised in a derivative expansion. The smoking gun of interactions with the environment is an enhanced but finite bispectrum close to the folded kinematical limit. We demonstrate the generality of our approach by matching our open effective theory to an explicit model. Our construction provides a standard model to simultaneously study phenomenological predictions as well as quantum information aspects of the inflationary dynamics.

Dynamical actions and q-representation theory for double-scaled SYK

We show that DSSYK amplitudes are reproduced by considering the quantum mechanics of a constrained particle on the quantum group SUq(1, 1). We construct its left-and right-regular representations, and show that the representation matrices reproduce two-sided wavefunctions and correlation functions of DSSYK. We then construct a dynamical action and path integral for a particle on SUq(1, 1), whose quantization reproduces the aforementioned representation theory. By imposing boundary conditions or constraining the system we find the q-analog of the Schwarzian and Liouville boundary path integral descriptions. This lays the technical groundwork for identifying the gravitational bulk description of DSSYK. We find evidence the theory in question is a sine dilaton gravity, which interestingly is capable of describing both AdS and dS quantum gravity.

Quantum information geometry of driven CFTs

Driven quantum systems exhibit a large variety of interesting and sometimes exotic phenomena. Of particular interest are driven conformal field theories (CFTs) which describe quantum many-body systems at criticality. In this paper, we develop both a spacetime and a quantum information geometry perspective on driven 2d CFTs. We show that for a large class of driving protocols the theories admit an alternative but equivalent formulation in terms of a CFT defined on a spacetime with a time-dependent metric. We prove this equivalence both in the operator formulation as well as in the path integral description of the theory. A complementary quantum information geometric perspective for driven 2d CFTs employs the so-called Bogoliubov-Kubo-Mori (BKM) metric, which is the counterpart of the Fisher metric of classical information theory, and which is obtained from a perturbative expansion of relative entropy. We compute the BKM metric for the universal sector of Virasoro excitations of a thermal state, which captures a large class of driving protocols, and find it to be a useful tool to classify and characterize different types of driving. For Möbius driving by the SL(2, ℝ) subgroup, the BKM metric becomes the hyperbolic metric on the disk. We show how the non-trivial dynamics of Floquet driven CFTs is encoded in the BKM geometry via Möbius transformations. This allows us to identify ergodic and non-ergodic regimes in the driving. We also explain how holographic driven CFTs are dual to driven BTZ black holes with evolving horizons. The deformation of the black hole horizon towards and away from the asymptotic boundary provides a holographic understanding of heating and cooling in Floquet CFTs.

Asymmetric transverse momentum broadening in an inhomogeneous medium

Gradient jet tomography in high-energy heavy-ion collisions utilizes the asymmetric transverse momentum broadening of a propagating parton in an inhomogeneous medium. Such broadening is studied within a path-integral description of the evolution of the Wigner distribution for a propagating parton in medium. Going beyond the eikonal approximation of multiple scattering, the evolution operator in the transverse direction can be expressed as the functional integration over all classical trajectories of a massive particle with the light-cone momentum $\ensuremath{\omega}$ as its mass. With a dipole approximation of the Wilson line correlation function, evolution with the light-cone time $t$ is determined by the jet transport coefficient $\stackrel{^}{q}$ that can vary with space and time. In a uniform medium with a constant ${\stackrel{^}{q}}_{0}$, the analytical solution to the Wigner distribution becomes a typical drifted Gaussian in both transverse momentum and coordinate with the diffusion width $\sqrt{{\stackrel{^}{q}}_{0}t}$ and $\sqrt{{\stackrel{^}{q}}_{0}{t}^{3}/3{\ensuremath{\omega}}^{2}}$, respectively. In the case of a simple Gaussian-like transverse inhomogeneity with a spatial width $\ensuremath{\sigma}$ on top of a uniform medium, the final asymmetrical momentum distribution can be calculated semianalytically. The transverse asymmetry defined for jet gradient tomography that characterizes the asymmetrical distribution is found to linearly correlate with the initial transverse position of the propagating parton within the domain of the inhomogeneity. It decreases with the parton energy $\ensuremath{\omega}$, increases with the propagation time initially and saturates when the diffusion distance is much larger than the size of the inhomogeneity or ${t}^{3}\ensuremath{\gg}3{\ensuremath{\omega}}^{2}{\ensuremath{\sigma}}^{2}/{\stackrel{^}{q}}_{0}$. The transverse momentum broadening due to the inhomogeneity also saturates at late times in contrast to the continued increase with time if the drifted diffusion in space is ignored.

Non-adiabatic ring polymer molecular dynamics in the phase space of the SU(N) Lie group.

We derive the non-adiabatic ring polymer molecular dynamics (RPMD) approach in the phase space of the SU(N) Lie Group. This method, which we refer to as the spin mapping non-adiabatic RPMD (SM-NRPMD), is based on the spin-mapping formalism for the electronic degrees of freedom (DOFs) and ring polymer path-integral description for the nuclear DOFs. Using the Stratonovich-Weyl transform for the electronic DOFs and the Wigner transform for the nuclear DOFs, we derived an exact expression of the Kubo-transformed time-correlation function (TCF). We further derive the spin mapping non-adiabatic Matsubara dynamics using the Matsubara approximation that removes the high frequency nuclear normal modes in the TCF and derive the SM-NRPMD approach from the non-adiabatic Matsubara dynamics by discarding the imaginary part of the Liouvillian. The SM-NRPMD method has numerical advantages compared to the original NRPMD method based on the Meyer-Miller-Stock-Thoss (MMST) mapping formalism due to a more natural mapping using the SU(N) Lie Group that preserves the symmetry of the original system. We numerically compute the Kubo-transformed position auto-correlation function and electronic population correlation function for three-state model systems. The numerical results demonstrate the accuracy of the SM-NRPMD method, which outperforms the original MMST-based NRPMD. We envision that the SM-NRPMD method will be a powerful approach to simulate electronic non-adiabatic dynamics and nuclear quantum effects accurately.

Path integral description and direct interaction approximation for elastic plate turbulence

In this work, we apply the Martin–Siggia–Rose path integral formalism to the equations of a thin elastic plate. Using a diagrammatic technique, we obtain the direct interaction approximation (DIA) equations to describe the evolutions of the correlation function and the response function of the fields. Consistent with previous results, we show that DIA equations for elastic plates can be derived from a non-markovian stochastic process and that in the weakly nonlinear limit, the DIA equations lead to the kinetic equation of wave turbulence theory. We expect that this approach will allow a better understanding of the statistical properties of wave turbulence and that DIA equations can open new avenues for understanding the breakdown of weakly nonlinear turbulence for elastic plates.

Clustering, collision, and relaxation dynamics in pure and doped helium nanoclusters: Density- vs particle-based approaches.

The clustering, collision, and relaxation dynamics of pristine and doped helium nanodroplets is theoretically investigated in cases of pickup and clustering of heliophilic argon, collision of heliophobic cesium atoms, and coalescence of two droplets brought into contact by their mutual long-range van der Waals interaction. Three approaches are used and compared with each other. The He time-dependent density functional theory method considers the droplet as a continuous medium and accounts for its superfluid character. The ring-polymer molecular dynamics method uses a path-integral description of nuclear motion and incorporates zero-point delocalization while bosonic exchange effects are ignored. Finally, the zero-point averaged dynamics approach is a mixed quantum-classical method in which quantum delocalization is described by attaching a frozen wavefunction to each He atom, equivalent to classical dynamics with effective interaction potentials. All three methods predict that the growth of argon clusters is significantly hindered by the helium host droplet due to the impeding shell structure around the dopants and kinematic effects freezing the growing cluster in metastable configurations. The effects of superfluidity are qualitatively manifested by different collision dynamics of the heliophilic atom at high velocities, as well as quadrupole oscillations that are not seen with particle-based methods, for droplets experiencing a collision with cesium atoms or merging with each other.

Real-time thermal self-energies: In the variational bases and spaces

In this work, we introduce a systematic study for studying the scalar propagator and tadpole self-energy by considering an arbitrary parameter σ that allows for a path integral description in real-time formalism (RTF). The closed time path formalism (CTP) and Thermofield Dynamics (TFD) are two popular choices for the parameter σ in the Feynman rules. We have constructed a scalar propagator and a tadpole self-energy in two different bases in the momentum space as well as the mixed space. The results show that the diagonal components of self-energy in both spaces for the 1/2 basis are the same in both approaches within RTF, whereas the other off-diagonal components of self-energy are different because they depend on the path parameter. On the other hand, the diagonal components of self-energy in both spaces for the new basis are not the same in both approaches within RTF, whereas the off-diagonal components of self-energy are vanishing. That means the new basis allows one to reduce the components for the quantities studied, like self-energy or other.

A study of non-linear Langevin dynamics under non-Gaussian noise with quartic cumulant

We consider a non-linear Langevin equation in presence of non-Gaussian noise originating from non-linear bath. We claim, the parameters in the Langevin equation are not physical. The physical parameters are obtained from a path-integral description of the system, where the Langevin parameters are related to the physical parameters by renormalisation flow equations. Then we compute both numerically and analytically the velocity two point function and show that it saturates to the bath temperature even in presence of non-linearity. We also find the velocity four point function numerically and show that it saturates to the analytically evaluated thermal velocity four point function when the non-linear fluctuation–dissipation relation (FDR) Chakrabarty and Chaudhuri (2019 SciPost Phys. 7 013) is satisfied. When the non-linear FDR, which is a manifestation of time reversibility of thermal bath, is violated then the system does not seem to thermalise. Rather, its velocity four point function settles to a steady state.

Influence Matrix Approach to Many-Body Floquet Dynamics

In this work, we introduce an approach to study quantum many-body dynamics, inspired by the Feynman-Vernon influence functional. Focusing on a family of interacting, Floquet spin chains, we consider a Keldysh path-integral description of the dynamics. The central object in our approach is the influence matrix (IM), which describes the effect of the system on the dynamics of a local subsystem. For translationally invariant models, we formulate a self-consistency equation for the influence matrix. For certain special values of the model parameters, we obtain an exact solution which represents a perfect dephaser (PD). Physically, a PD corresponds to a many-body system that acts as a perfectly Markovian bath on itself: at each period, it measures every spin. For the models considered here, we establish that PD points include dual-unitary circuits investigated in recent works. In the vicinity of PD points, the system is not perfectly Markovian, but rather acts as a bath with a short memory time. In this case, we demonstrate that the self-consistency equation can be solved using matrix-product states (MPS) methods, as the IM temporal entanglement is low. A combination of analytical insights and MPS computations allows us to characterize the structure of the influence matrix in terms of an effective "statistical-mechanics" description. We finally illustrate the predictive power of this description by analytically computing how quickly an embedded impurity spin thermalizes. The influence matrix approach formulated here provides an intuitive view of the quantum many-body dynamics problem, opening a path to constructing models of thermalizing dynamics that are solvable or can be efficiently treated by MPS-based methods, and to further characterizing quantum ergodicity or lack thereof.

Chemical Langevin equation: A path-integral view of Gillespie's derivation.

In 2000, Gillespie rehabilitated the chemical Langevin equation (CLE) by describing two conditions that must be satisfied for it to yield a valid approximation of the chemical master equation (CME). In this work, we construct an original path-integral description of the CME and show how applying Gillespie's two conditions to it directly leads to a path-integral equivalent to the CLE. We compare this approach to the path-integral equivalent of a large system size derivation and show that they are qualitatively different. In particular, both approaches involve converting many sums into many integrals, and the difference between the two methods is essentially the difference between using the Euler-Maclaurin formula and using Riemann sums. Our results shed light on how path integrals can be used to conceptualize coarse-graining biochemical systems and are readily generalizable.

Path-integral description of quantum nonlinear optics in arbitrary media

We present a method, based on Feynman path integrals, to describe the propagation and properties of the quantised electromagnetic field in an arbitrary, nonlinear medium. We provide a general theory, valid for any order of optical nonlinearity, and we then specialise the the case of second order nonlinear processes. In particular, we show, that second-order nonlinear processes in arbitrary media, under the undepleted pump approximation, can be described by an effective free electromagnetic field, propagating in a vacuum, dressed by the medium itself. Moreover, we show,that the probability of such processes to occur is related to the biphoton propagator, which contains informations about the structure of the medium, its nonlinear properties, and the structure of the pump beam.

Quantum incompleteness of inflation

Inflation is most often described using quantum field theory (QFT) on a fixed, curved spacetime background. Such a description is valid only if the spatial volume of the region considered is so large that its size and shape moduli behave classically. However, if we trace an inflating universe back to early times, the volume of any comoving region of interest -- for example the present Hubble volume -- becomes exponentially small. Hence, quantum fluctuations in the trajectory of the background cannot be neglected at early times. In this paper, we develop a path integral description of a flat, inflating patch (approximated as de Sitter spacetime), treating both the background scale factor and the gravitational wave perturbations quantum mechanically. We find this description fails at small values of the initial scale factor, because \emph{two} background saddle point solutions contribute to the path integral. This leads to a breakdown of QFT in curved spacetime, causing the fluctuations to be unstable and out of control. We show the problem may be alleviated by a careful choice of quantum initial conditions, for the background and the fluctuations, provided that the volume of the initial, inflating patch is larger than $\gg H^{-1}$ in Planck units with $H$ the Hubble constant at the start of inflation. The price of the remedy is high: not only the inflating background, but also the stable, Bunch-Davies fluctuations must be input by hand. Our discussion emphasizes that, even if the inflationary scale is far below the Planck mass, new physics is required to explain the initial quantum state of the universe.

Squeezed-field path-integral description of second sound in Bose-Einstein condensates

We propose a generalization of the Feynman path integral using squeezed coherent states. We apply this approach to the dynamics of Bose-Einstein condensates, which gives an effective low energy description that contains both a coherent field and a squeezing field. We derive the classical trajectory of this action, which constitutes a generalization of the Gross Pitaevskii equation, at linear order. We derive the low energy excitations, which provides a description of second sound in weakly interacting condensates as a squeezing oscillation of the order parameter. This interpretation is also supported by a comparison to a numerical c-field method.

State dependent ring polymer molecular dynamics for investigating excited nonadiabatic dynamics

A recently proposed nonadiabatic ring polymer molecular dynamics (NRPMD) approach has shown to provide accurate quantum dynamics by incorporating explicit state descriptions and nuclear quantizations. Here, we present a rigorous derivation of the NRPMD Hamiltonian and investigate its performance on simulating excited state nonadiabatic dynamics. Our derivation is based on the Meyer-Miller-Stock-Thoss mapping representation for electronic states and the ring-polymer path-integral description for nuclei, resulting in the same Hamiltonian proposed in the original NRPMD approach. In addition, we investigate the accuracy of using NRPMD to simulate the photoinduced nonadiabatic dynamics in simple model systems. These model calculations suggest that NRPMD can alleviate the zero-point energy leakage problem that is commonly encountered in the classical Wigner dynamics and provide accurate excited state nonadiabatic dynamics. This work provides a solid theoretical foundation of the promising NRPMD Hamiltonian and demonstrates the possibility of using the state-dependent RPMD approach to accurately simulate electronic nonadiabatic dynamics while explicitly quantizing nuclei.

Effective driven dynamics for one-dimensional conditioned Langevin processes in the weak-noise limit

In this work we focus on fluctuations of time-integrated observables for a particle diffusing in a one-dimensional periodic potential in the weak-noise asymptotics. Our interest goes to rare trajectories presenting an atypical value of the observable, that we study through a biased dynamics in a large-deviation framework. We determine explicitly the effective probability-conserving dynamics which makes rare trajectories of the original dynamics become typical trajectories of the effective one. Our approach makes use of a weak-noise path-integral description in which the action is minimised by the rare trajectories of interest. For ‘current-type’ additive observables, we find criteria for the emergence of a propagative trajectory minimising the action for large enough deviations, revealing the existence of a dynamical phase transition at a fluctuating level, whose singular behaviour is between first and second order. In addition, we provide a new method to determine the scaled cumulant generating function of the observable without having to optimise the action. It allows one to show that the weak-noise and the large-time limits commute in this problem. Finally, we show how the biased dynamics can be mapped in practice to an explicit effective driven dynamics, which takes the form of a driven Langevin dynamics in an effective potential. The non-trivial shape of this effective potential is key to understand the link between the dynamical phase transition in the large deviations of current and the standard depinning transition of a particle in a tilted potential.

Causality in the Classical Limit for Quantum Electrodynamics

We use the path integral form of quantum electrodynamics (QED) to show that a causal classical limit to QED can be derived by functionally integrating over the photon coordinates, starting from an initial photon vacuum and ending in a final coherent radiation state driven by the anticipated classical charged particle trajectories. The resulting charged particle transition amplitude depends only on particle coordinates. When the \( {\hbar} \, \to \,0 \) limit is taken, only those particle paths that are not constrained by the final radiation state are varied. These results demonstrate that the collapse from an infinity of charged particle paths, a path integral description, to causally interacting classical trajectories, a stationary-action description, is critically dependent on including final coherent state radiation and maintaining the distinction between particle paths that are free to vary and those trajectories that can be monitored by the final state radiation.

Simultaneous continuous measurement of noncommuting observables: Quantum state correlations

We consider the temporal correlations of the quantum state of a qubit subject to simultaneous continuous measurement of two non-commuting qubit observables. Such qubit state correlators are defined for an ensemble of qubit trajectories, which has the same fixed initial state and can also be optionally constrained by a fixed final state. We develop a stochastic path integral description for the continuous quantum measurement and use it to calculate the considered correlators. Exact analytic results are possible in the case of ideal measurements of equal strength and are also shown to agree with solutions obtained using the Fokker-Planck equation. For a more general case with decoherence effects and inefficiency, we use a diagrammatic approach to find the correlators perturbatively in the quantum efficiency. We also calculate the state correlators for the quantum trajectories which are extracted from readout signals measured in a transmon qubit experiment, by means of the quantum Bayesian state update. We find an excellent agreement between the correlators based on the experimental data and those obtained from our analytical and numerical results.