Path Integral Formula Research Articles

An Impulse Ghost Fluid Method for Simulating Two-Phase Flows

This paper introduces a two-phase interfacial fluid model based on the impulse variable to capture complex vorticity-interface interactions. Our key idea is to leverage bidirectional flow map theory to enhance the transport accuracy of both vorticity and interfaces simultaneously and address their coupling within a unified Eulerian framework. At the heart of our framework is an impulse ghost fluid method to solve the two-phase incompressible fluid characterized by its interfacial dynamics. To deal with the history-dependent jump of gauge variables across a dynamic interface, we develop a novel path integral formula empowered by spatiotemporal buffers to convert the history-dependent jump condition into a geometry-dependent jump condition when projecting impulse to velocity. We demonstrate the efficacy of our approach in simulating and visualizing several interface-vorticity interaction problems with cross-phase vortical evolution, including interfacial whirlpool, vortex ring reflection, and leapfrogging bubble rings.

Particle-Laden Fluid on Flow Maps

We propose a novel framework for simulating ink as a particle-laden flow using particle flow maps. Our method addresses the limitations of existing flow-map techniques, which struggle with dissipative forces like viscosity and drag, thereby extending the application scope from solving the Euler equations to solving the Navier-Stokes equations with accurate viscosity and laden-particle treatment. Our key contribution lies in a coupling mechanism for two particle systems, coupling physical sediment particles and virtual flow-map particles on a background grid by solving a Poisson system. We implemented a novel path integral formula to incorporate viscosity and drag forces into the particle flow map process. Our approach enables state-of-the-art simulation of various particle-laden flow phenomena, exemplified by the bulging and breakup of suspension drop tails, torus formation, torus disintegration, and the coalescence of sedimenting drops. In particular, our method delivered high-fidelity ink diffusion simulations by accurately capturing vortex bulbs, viscous tails, fractal branching, and hierarchical structures.

A path integral formula of quantum gravity emergent from entangled local structures

We couple to group field theory (GFT) a scalar field that encodes the entanglement between manifold sites. The scalar field provides a relational clock that enables the derivation of the Hamiltonian of the system from the GFT action. Inspecting the Hamiltonian, we show that a theory of emergent gravity arises, and that this can be recast according to the Ashtekar’s formulation of general relativity. The evolution of the GFT observables is regulated by the Shrödinger equation generated by the Hamiltonian. This is achieved by imposing a renormalization group (RG) flow that corresponds to a simplified Ricci flow. As a consequence of the quantization procedure, the Hamiltonian is recovered to be non-Hermitian, and can be related to the complex action formalism, in which the initial conditions and the related future evolution of the systems are dictated by the imaginary part of the action.

The path integral formula for the stochastic evolutionary game dynamics

Although the long-term behavior of stochastic evolutionary game dynamics in finite populations has been fully investigated, its evolutionary characteristics in a limited period of time is still unclear. In order to answer this question, we introduce the formulation of the path integral approach for evolutionary game theory. In this framework, the transition probability is the sum of all the evolutionary paths. The path integral formula of the transition probability is expected to be a new mathematical tool to explore the stochastic game evolutionary dynamics. As an example, we derive the transition probability for stochastic evolutionary game dynamics by the path integral in a limited period of time with the updating rule of the Wright-Fisher process.

Matrix-valued Schrödinger operators over finite adeles

Let [Formula: see text] be an algebraic number field. With [Formula: see text] we associate the ring of finite adeles [Formula: see text] In this paper we give a path integral formula for the propagator of a quantum mechanical system over the abelian group [Formula: see text] Specifically, we consider matrix-valued Hamiltonian operators [Formula: see text] where [Formula: see text] is the Vladimirov operator and [Formula: see text] is a non-negative definite potential. The free part of the Hamiltonian gives rise to a measure on the Skorokhod space of paths which allows us to prove the Feynman–Kac formula for the Schrödinger semigroup generated by [Formula: see text] This formula is given in terms of the ordered time exponentials.

From Dual Unitarity to Generic Quantum Operator Spreading.

Dual-unitary circuits are paradigmatic examples of exactly solvable yet chaotic quantum many-body systems, but solvability naturally goes along with a degree of nongeneric behavior. By investigating the effect of weakly broken dual unitarity on the spreading of local operators, we study whether, and how, small deviations from dual unitarity recover fully generic many-body dynamics. We present a discrete path-integral formula for the out-of-time-order correlator and recover a butterfly velocity smaller than the light-cone velocity, v_{B}<v_{LC}, and a diffusively broadening operator front, two generic features of ergodic quantum spin chains absent in dual-unitary circuit dynamics. The butterfly velocity and diffusion constant are determined by a small set of microscopic quantities, and the operator entanglement of the gates has a crucial role.

A new algorithm for computing path integrals and weak approximation of SDEs inspired by large deviations and Malliavin calculus

The paper gives a novel path integral formula inspired by large deviation theory and Malliavin calculus. The proposed finite-dimensional approximation of integrals on path space will be a new higher-order weak approximation of multidimensional stochastic differential equations where the dominant part of the local expansion is governed by Varadhan's geodesic distance and the correction terms are given as Malliavin weights. An optimal truncation of asymptotic expansion is used to reduce computational effort. Kusuoka's estimate is applied to justify the finite-dimensional approximation of path integrals. An efficient simulation method is provided with the algorithm. Numerical results are shown to verify the effectiveness.

K-space theory and convective gains of stimulated Raman side scattering.

We have developed a k-space theory for stimulated Raman scattering. The theory is used to calculate the convective gain of stimulated Raman side scattering (SRSS) in order to clarify the discrepancies found between previous gain formulas proposed in the literature. The gains are dramatically modified by the eigenvalue of SRSS, and the maximum gain occurs not at the perfect wave-number matching condition but instead at a wave number having a small deviation that is associated with the eigenvalue. Numerical solutions of the k-space theory equationsare compared with, and used to verify, analytically derived gains. We demonstrate connections to the existing path integral theories, and we derive a similar path integral formula in k space.

Medium-modifications of g → cc splitting

We study medium-modifications of the gluon splitting into a quark and anti-quark pair. Applying the Baier-Dokshitzer-Mueller-Peigné-Schiff and Zakharov (BDMPS-Z) formalism, we derive a path-integral formula for the in-medium g → qq¯ splitting function in the close-to-eikonal limit. Our analysis shows that there are two qualitatively different medium effects: transverse momentum broadening of qq¯ pairs and enhanced production of such pairs. We note that both effects are numerically sizeable if the average momentum transfer from the medium to the parton is at the quark mass scale. In ultra-relativistic heavy-ion collisions, this condition is realized by charm quarks, therefore we focus our numerical analysis on the medium-modifications of g → cc¯ splitting.

Exact Solutions of Fractional Order Oscillation Equation with Two Fractional Derivative Terms

The fractional oscillation equation with two fractional derivative terms in the sense of Caputo, where the orders alpha and beta satisfy 1<alpha le 2 and 0<beta le 1, is investigated, and the unit step response and two initial value responses are obtained in two forms by using different methods of inverse Laplace transform. The first method yields series solutions with the nonnegative powers of t, which converge fast for small t. The second method is our emphasis, where the complex path integral formula of the inverse Laplace transform is used. In order to determine singularities of integrand we first seek for the roots of the characteristic equation, which is a transcendental equation with four parameters, two coefficients and two noninteger power exponents. The existence conditions and properties of the roots on the principal Riemann surface are given. Based on the results on the characteristic equation, we derive these responses as a sum of a classical exponentially damped oscillation, which vanishes in an indicated case, and an infinite integral of the Laplace type, which converges fast for large t, with a steady component in the unit step response. Asymptotic behaviors of solutions for large t are derived as algebraic decays in negative power laws characterized by the orders alpha and beta. The fractional system exhibits a magical transition from oscillation to monotonic decay in negative power law.

A path integral formula for certain fourth order elliptic operators

A path integral formula for certain fourth order elliptic operators

Discriminant and Root Trajectories of Characteristic Equation of Fractional Vibration Equation and Their Effects on Solution Components

The impulsive response of the fractional vibration equation z′′(t)+bDtαz(t)+cz(t)=F(t), b&gt;0,c&gt;0,0≤α≤2, is investigated by using the complex path-integral formula of the inverse Laplace transform. Similar to the integer-order case, the roots of the characteristic equation s2+bsα+c=0 must be considered. It is proved that for any b&gt;0, c&gt;0 and α∈(0,1)∪(1,2), the characteristic equation always has a pair of conjugated simple complex roots with a negative real part on the principal Riemann surface. Particular attention is paid to the problem as to how the couple conjugated complex roots approach the two roots of the integer case α=1, especially to the two different real roots in the case of b2−4c&gt;0. On the upper-half complex plane, the root s(α) is investigated as a function of order α and with parameters b and c, and so are the argument θ(α), modulus r(α), real part λ(α) and imaginary part ω(α) of the root s(α). For the three cases of the discriminant b2−4c: &gt;0, =0 and &lt;0, variations of the argument and modulus of the roots according to α are clarified, and the trajectories of the roots are simulated. For the case of b2−4c&lt;0, the trajectories of the roots are further clarified according to the change rates of the argument, real part and imaginary part of root s(α) at α=1. The solution components, i.e., the residue contribution and the Hankel integral contribution to the impulsive response, are distinguished for the three cases of the discriminant.

Relativistic quantum field theory of stochastic dynamics in the Hilbert space

We develop an action formulation of stochastic dynamics in the Hilbert space. By generalizing the Wiener process into 1+3-dimensional spacetime, we define a Lorentz-invariant random field. By coupling the random to quantum fields, we obtain a random-number action which has the statistical spacetime translation and Lorentz symmetries. The canonical quantization of the theory results in a Lorentz-invariant equation of motion for the state vector or density matrix. We derive the path integral formula of $S$-matrix and the diagrammatic rules for both the stochastic free field theory and stochastic $\phi^4$-theory. The Lorentz invariance of the random $S$-matrix is strictly proved. We then develop a diagrammatic technique for calculating the density matrix. Without interaction, we obtain the exact $S$-matrix and density matrix. With interaction, we prove a simple relation between the density matrices of stochastic and conventional $\phi^4$-theory. Our formalism leads to an ultraviolet divergence which has the similar origin as that in QFT. The divergence is canceled by renormalizing the coupling strength to random field. We prove that the stochastic QFT is renormalizable even in the presence of interaction. In the models with a linear coupling between random and quantum fields, the random field excites particles out of the vacuum, driving the universe towards an infinite-temperature state. The number of excited particles follows the Poisson distribution. The collision between particles is not affected by the random field. But the signals of colliding particles are gradually covered by the background excitations caused by random field.

Unimodular gauge and ADM gravity path integral

This paper proposes a definition of gravitational observables and of their path integral formula within the framework of ADM foliation and the choice of unimodular gauge classes. The method enforces a BRST invariant gauge fixing of the lapse and shift fields. It yields the quantum level extension of the known classical property that the conformal classes of internal metrics of constant Lorentz time leafs define the gravitational physical degrees of freedom.

FKN Formula and Ground State Energy for the Spin Boson Model with External Magnetic Field

We consider the spin boson model with external magnetic field. We prove a path integral formula for the heat kernel, known as Feynman–Kac–Nelson (FKN) formula. We use this path integral representation to express the ground state energy as a stochastic integral. Based on this connection, we determine the expansion coefficients of the ground state energy with respect to the magnetic field strength and express them in terms of correlation functions of a continuous Ising model. From a recently proven correlation inequality, we can then deduce that the second order derivative is finite. As an application, we show existence of ground states in infrared-singular situations.

Inspiraling corrugation-induced quantum effects on neutron star binary plane

We use the path-integral formula and investigate some dynamical quantum effects induced by the inspiraling lateral corrugation of orbital plane in gravitationally bound neutron star (NS) binaries, with orbital separation of 109 m. Based on Dewitt's approach, we calculate the gravitational Casimir energy cost of the binary plane, which consists of statically gravitational effects and deformation-induced effects. It is found that the static effects include a term coming from the self-gravity of the orbital plane and the contribution of Newtonian gravitational potential of the binary system. While the deformation-induced effect also results from two parts, i.e. the instability of orbital binding energy, scaling as 1(R−r)2, and the dynamically Casimir energy cost of the orbital binding energy, decaying as 1(R−r)4. The dynamically gravitational Casimir phenomena and the corresponding energy cost modify the spiral-in orbital motion of the binary and thus the frequency of released gravitational waves (GWs). We consider the mechanical response of two NS components and qualitatively study the corrections to the orbital motion of the system and the GW frequencies. It is found that the dynamical Casimir effects exert a dissipative force on the binary plane, depending on the frequency of GWs. The resultant dissipation may enhance with the decaying separation and increasing GW frequencies, which subsequently accelerates the orbital decay of the binary. However, the dissipation rate just has an order of 10−70 eV/s. So the corrections to the dynamics of NS binaries are very marginal, by considering the wide separation, the cosmological coalescence time, and low-frequency GWs of the system.

Revisiting relativistic magnetohydrodynamics from quantum electrodynamics

We provide a statistical mechanical derivation of relativistic magnetohydrodynamics on the basis of (3 + 1)-dimensional quantum electrodynamics; the system endowed with a magnetic one-form symmetry. The conservation laws and constitutive relations are presented in a manifestly covariant way with respect to the general coordinate transformation. The method of the local Gibbs ensemble (or nonequilibrium statistical operator) combined with the path-integral formula for a thermodynamic functional enables us to obtain exact forms of constitutive relations. Applying the derivative expansion to exact formulas, we derive the first-order constitutive relations for nonlinear relativistic magnetohydrodynamics. Our results for the QED plasma preserving parity and charge-conjugation symmetries are equipped with two electrical resistivities and five (three bulk and two shear) viscosities. We also show that those transport coefficients satisfy the Onsager’s reciprocal relation and a set of inequalities, indicating semi-positivity of the entropy production rate consistent with the local second law of thermodynamics.

Feynman path integrals for magnetic Schrödinger operators on infinite weighted graphs

We prove a Feynman path integral formula for the unitary group exp(—itLυ,θ), t ≥ 0, associated with a discrete magnetic Schrodinger operator Lυ,θ on a large class of weighted infinite graphs. As a consequence, we get a new Kato-Simon estimate $$\left| {\exp \left({- it{L_{v,\theta}}} \right)\left({x,y} \right)} \right| \le \exp \left({- t{L_{- \deg ,0}}} \right)\left({x,y} \right),$$ which controls the unitary group uniformly in the potentials in terms of a Schrodinger semigroup, where the potential deg is the weighted degree function of the graph.

Mini-Workshop: Recent Progress in Path Integration on Graphs and Manifolds

Ever since Richard Feynman’s PhD thesis, path integrals have played a decisive role in mathematical physics. While it is well-known that such formulae can hold only formally, it was Mark Kac who realized that by replacing the unitary group by the heat semigroup, one obtains well-defined and rigorous formulae. Following this pioneering work, Feynman-Kac path integral formulae have been adapted to several situations and generalized into several directions providing the central focus of this workshop.

A Rigorous Mathematical Construction of Feynman Path Integrals for the Schrödinger Equation with Magnetic Field

A Feynman path integral formula for the Schrödinger equation with magnetic field is rigorously mathematically realized in terms of infinite dimensional oscillatory integrals. We show (by the example of a linear vector potential) that the requirement of the independence of the integral on the approximation procedure forces the introduction of a counterterm to be added to the classical action functional. This provides a natural explanation for the appearance of a Stratonovich integral in the path integral formula for both the Schrödinger and heat equation with magnetic field.