Analytic Continuation Research Articles

Analytic Continuation of Multipoint Correlation Functions

AbstractConceptually, the Matsubara formalism (MF), using imaginary frequencies, and the Keldysh formalism (KF), formulated in real frequencies, give equivalent results for systems in thermal equilibrium. The MF has less complexity and is thus more convenient than the KF. However, computing dynamical observables in the MF requires the analytic continuation from imaginary to real frequencies. The analytic continuation is well‐known for two‐point correlation functions (having one frequency argument), but, for multipoint correlators, a straightforward recipe for deducing all Keldysh components from the MF correlator had not been formulated yet. Recently, a representation of MF and KF correlators in terms of formalism‐independent partial spectral functions and formalism‐specific kernels was introduced by Kugler, Lee, and von Delft [Phys. Rev. X 11, 041006 (2021)]. This representation is used to formally elucidate the connection between both formalisms. How a multipoint MF correlator can be analytically continued to recover all partial spectral functions and yield all Keldysh components of its KF counterpart is shown. The procedure is illustrated for various correlators of the Hubbard atom.

A new method for solving the linearized 1D Vlasov–Poisson system yielding a new class of solutions

We describe a new method for solving the linearized 1D Vlasov–Poisson system by using properties of Cauchy-type integrals. Our method remedies critical flaws of the two standard methods, reveals a previously unrecognized Gaussian-in-time-like decay, and can also account for an externally applied electric field. The Landau approximation involves deforming the Bromwich contour around the poles closest to the real axis due to the analytically continued dielectric function, finding the long-time behavior for a stable system: Landau damping. Jackson's generalization encircles all poles while sending the contour to infinity, assuming its contribution vanishes, which is not true in general. This gives incorrect solutions for physically reasonable configurations and can exhibit pathological behavior, of which we show examples. The van Kampen method expresses the solution for a stable equilibrium as a continuous superposition of waves, resulting in an opaque integral. Case's generalization includes unstable systems and predicts a decaying discrete mode for each growing discrete mode, an apparent contradiction to both the Jackson solution and ours. We show, without imposing additional constraints, that the decaying modes are never present in the time evolution due to an exact cancellation with part of the continuum. Our solution is free of integral expressions, is obtained using algebra and Laurent series expansions, does not rely on analytic continuations, and results in a correct asymptotically convergent form in the case of infinite sums. The analysis used can be readily applied in higher-dimensional, electromagnetic systems and also provides a new technique for evaluating certain inverse Laplace transforms.

Modular Invariance of (Logarithmic) Intertwining Operators

Let V be a C2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$C_2$$\\end{document}-cofinite vertex operator algebra without nonzero elements of negative weights. We prove the conjecture that the spaces spanned by analytic extensions of pseudo-q-traces (q=e2πiτ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$q=e^{2\\pi i\ au }$$\\end{document}) shifted by -c24\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$-\\frac{c}{24}$$\\end{document} of products of geometrically-modified (logarithmic) intertwining operators among grading-restricted generalized V-modules are invariant under modular transformations. The convergence and analytic extension result needed to formulate this conjecture and some consequences on such shifted pseudo-q-traces were proved by Fiordalisi (Logarithmic intertwining operator and genus-one correlation functions, 2015) and Fiordalisi (Commun Contemp Math 18:1650026, 2016) using the method developed in Huang (Commun Contemp Math 7:649–706, 2005). The method that we use to prove this conjecture is based on the theory of the associative algebras AN(V)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$A^{N}(V)$$\\end{document} for N∈N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N\\in \\mathbb {N}$$\\end{document}, their graded modules and their bimodules introduced and studied by the author in Huang (Associative algebras and the representation theory of grading-restricted vertex algebras, 2020) and Huang (Commun Math Phys 396:1–44, 2022). This modular invariance result gives a construction of C2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$C_2$$\\end{document}-cofinite genus-one logarithmic conformal field theories from the corresponding genus-zero logarithmic conformal field theories.

Increasing stability in the n‐dimensional inverse source problem with multifrequencies

In this paper, we show for the first time the increasing stability of the inverse source problem for the n‐dimensional Helmholtz equation at multiple wave numbers, which is different from the two‐ or three‐dimensional Helmholtz equation. In addition, we develop a new, unified approach to study increasing stability in any dimension. The method is based on the Fourier transform and explicit bounds for analytic continuation.

Field theory in Rindler frame and more on the correspondence with thermal field theory formalisms

Considering two accelerated observers with same acceleration in two timelike wedges of Rindler frame we calculate the Feynman-like propagators for a real scalar field in a thermal bath with respect to the Minkowski vacuum. Only the same wedge correlators are symmetric under the exchange of the real thermal bath and Unruh thermal bath, while the cross-wedge ones are not. Interestingly, they contain a cross term which is a collective effects of acceleration and thermal nature of field. Particularly the zero temperature description along with no analytic continuation between coordinates in right and left Rindler wedges, as expected, corresponds to usual thermofield-double formalism. However, unlike in later formulation, the two fields are now parts of the original system. Moreover it bears the features of a spacial case of closed-time formalism (CTP) where the Keldysh contour is along the increasing Rindler time in the respective Rindler wedges. Interestingly, we observe a new feature that the analytic continuation between the wedges provides the two more spacial cases of CTP. Hence Rindler-frame-field theory seems to be a viable candidate to deal thermal theory of fields and may illuminate the search for a bridge between the usual existing formalisms.

Exact WKB analysis for PT -symmetric quantum mechanics: Study of the Ai-Bender-Sarkar conjecture

We consider exact Wentzel-Kramers-Brillouin analysis to a PT symmetric quantum mechanics defined by the potential, V(x)=ω2x2+gx2(ix)ϵ=2 with ω∈R≥0, g∈R>0. We in particular aim to verify a conjecture proposed by Ai-Bender-Sarkar (ABS), that pertains to a relation between D-dimensional PT-symmetric theories and analytic continuation (AC) of Hermitian theories concerning the energy spectrum or Euclidean partition function. For the purpose, we construct energy quantization conditions by exact Wentzel-Kramers-Brillouin analysis and write down their transseries solution by solving the conditions. By performing alien calculus to the energy solutions, we verify validity of the ABS conjecture and seek a possibility of its alternative form by Borel resummation theory if it is violated. Our results claim that the validity of the ABS conjecture drastically changes depending on whether ω>0 or ω=0: If ω>0, then the ABS conjecture is violated when exceeding the semiclassical level of the first nonperturbative order, but its alternative form is constructable by Borel resummation theory. The PT and the AC energies are related to each other by a one-parameter Stokes automorphism, and a median resummed form, which corresponds to a formal exact solution, of the AC energy (resp. PT energy) is directly obtained by acting Borel resummation to a transseries solution of the PT energy (resp. AC energy). If ω=0, then, with respect to the inverse energy level-expansion, not only perturbative/nonperturbative structures of the PT and the AC energies but also their perturbative parts do not match with each other. These energies are independent solutions, and no alternative form of the ABS conjecture can be reformulated by Borel resummation theory. Published by the American Physical Society 2024

Black hole perturbation theory meets CFT2 : Kerr-Compton amplitudes from Nekrasov-Shatashvili functions

We present a novel study of Kerr Compton amplitudes in a partial wave basis in terms of the Nekrasov-Shatashvili (NS) function of the confluent Heun equation (CHE). Remarkably, NS-functions enjoy analytic properties and symmetries that are naturally inherited by the Compton amplitudes. Based on this, we characterize the analytic dependence of the Compton phase shift in the Kerr spin parameter and provide a direct comparison to the standard post-Minkowskian (PM) perturbative approach within general relativity (GR). We also analyze the universal large frequency behavior of the relevant characteristic exponent of the CHE—also known as the renormalized angular momentum—and find agreement with numerical computations. Moreover, we discuss the analytic continuation in the harmonics quantum number ℓ of the partial wave, and show that the limit to the physical integer values commutes with the PM expansion of the observables. Finally, we obtain the contributions to the tree-level, point-particle, gravitational Compton amplitude in a covariant basis through O(aBH8), without the need to take the superextremal limit for Kerr spin. Published by the American Physical Society 2024

Cut structures and an observable singularity in the three-body threshold dynamics: The Tcc+ case

The three-body threshold effect, the distinctive and intriguing nonperturbative dynamics in the low-energy hadron-hadron scattering, has acquired compelling significance in the wake of the recent observation of the double-charm tetraquark Tcc+. This dynamics is characterized by the emergence of singular points and branch cuts within the interaction potential, occurring when the on-shell condition of the mediated particle is satisfied. The presence of these potential singularities indicates that the system is no longer Hermitian and also poses intractable challenges in obtaining exact solutions for dynamical scattering amplitudes. In this work, we develop a complex scaled Lippmann-Schwinger equation as an operation of analytical continuation of the T matrix to resolve this problem. Through a practical application to the DD*→DD* process, we reveal complicated cut structures of the three-body threshold dynamics in the complex plane, primarily stemming from the one-pion exchange. Notably, our methodology succeeds in reproducing the Tcc+ structure, in alignment with the quasibound pole derived from the complex scaling method within the Schrödinger equation framework. More remarkably, after solving the on-shell T matrix on the positive real axis of momentum plane, we find an extra new structure in the DD* mass spectrum, which arises from a right-hand cut at a physical pion mass and should be observable in lattice QCD simulations and future high-energy experiments. Published by the American Physical Society 2024

A Precision Benchmark Suite for Nuclear Reactor Point Kinetics Equations via Converged Accelerated Taylor Series (CATS)

Extreme benchmarks of 10 or more places for the point kinetics equations for time-dependent nuclear reactor power transients are rare. Therefore, to establish an extreme benchmark, we employ a Taylor series (TS) with continuous analytical continuation to solve the ordinary differential equations of point kinetics including feedback. Nonlinear Wynn-epsilon convergence acceleration confirms the highly precise solutions for neutron and precursor densities. Through adaptive partitioning of time intervals, the proposed Converged Accelerated Taylor Series, or CATS algorithm in double precision, automatically performs successive mesh refinement to obtain high-precision initial conditions for each subinterval, with the intent to reduce propagation error. Confirmation of 10 to 12 places comes from comparison to the BEFD (Backward Euler Finite Difference) algorithm in quadruple precision also developed by the author. We report benchmark results for common cases found in the literature including step, ramp, zigzag, and sinusoidal prescribed reactivity insertions and insertions with nonlinear adiabatic Doppler feedback. We also establish a suite of new prescribed reactivity insertions and insertions with feedback, based on reactivities with Taylor series representations as suggested by the CATS algorithm.

Quantum soft waveguides with resonances induced by broken symmetry

We consider two-dimensional, non-relativistic quantum system with asymptotically straight soft waveguide. We show that the local deformation of the symmetric waveguide can lead to the emerging of the embedded eigenvalues living in the continuous spectrum. The main problem of this paper is devoted to the analysis of weak perturbation of the symmetric system. We show that the original embedded eigenvalues turn to the second sheet of the resolvent analytic continuation and constitute resonances. We describe the asymptotics of the real and imaginary components of the complex resonant pole depending on deformation. Finally, we generalize the problem to three dimensional system equipped with a soft layer.

Cosmology in Lorentzian Regge calculus: causality violations, massless scalar field and discrete dynamics

We develop a model of spatially flat, homogeneous and isotropic cosmology in Lorentzian Regge calculus, employing four-dimensional Lorentzian frusta as building blocks. By examining the causal structure of the discrete spacetimes obtained by gluing such four-frusta in spatial and temporal direction, we find causality violations if the sub-cells connecting spatial slices are spacelike. A Wick rotation to the Euclidean theory can be defined globally by a complexification of the variables and an analytic continuation of the action. Introducing a discrete free massless scalar field, we study its equations of motion and show that it evolves monotonically. Furthermore, in a continuum limit, we obtain the equations of a homogeneous scalar field on a spatially flat Friedmann background. Vacuum solutions to the causally regular Regge equations are static and flat and show a restoration of time reparametrisation invariance. In the presence of a scalar field, the height of a frustum is a dynamical variable that has a solution if causality violations are absent and if an inequality relating geometric and matter boundary data is satisfied. Edge lengths of cubes evolve monotonically, yielding a contracting or an expanding branch of the Universe. In a small deficit angle expansion, the system can be deparametrised via the scalar field and a continuum limit of the discrete theory can be defined which we show to yield the relational Friedmann equation. These properties are obstructed if higher orders of the deficit angle are taken into account. Our results suggest that the inclusion of timelike sub-cells is necessary for a causally regular classical evolution in this symmetry restricted setting. Ultimately, this works serves as a basis for forthcoming investigations on the cosmological path integral within the framework of effective spin foams.

Double-Step Shape Invariance of Radial Jacobi-Reference Potential and Breakdown of Conventional Rules of Supersymmetric Quantum Mechanics

The paper reveals some remarkable form-invariance features of the ‘Jacobi-reference’ canonical Sturm–Liouville equation (CSLE) in the particular case of the density function with the simple pole at the origin. It is proven that the CSLE under consideration preserves its form under the two second-order Darboux–Crum transformations (DCTs) with the seed functions represented by specially chosen pairs of ‘basic’ quasi-rational solutions (q-RSs), i.e., such that their analytical continuations do not have zeros in the complex plane. It is proven that both transformations generally either increase or decrease by 2 the exponent difference (ExpDiff) for the mentioned pole while keeping two other parameters unchanged. The change is more complicated in the latter case if the ExpDiff for the pole of the original CSLE at the origin is smaller than 2. It was observed that the DCTs in question do not preserve bound energy levels according to the conventional supersymmetry (SUSY) rules. To understand this anomaly, we split the DCT in question into the two sequential Darboux deformations of the Liouville potentials associated with the CSLEs of our interest. We found that the first Darboux transformation turns the initial CSLE into the Heun equation written in the canonical form while the second transformation brings us back to the canonical form of the hypergeometric equation. It is shown that the first of these transformations necessarily places the mentioned ExpDiff into the limit-circle (LC) range and then the second transformation keeps the pole within the LC region, violating the conventional prescriptions of SUSY quantum mechanics.

Large N Schur index of N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 SYM from semiclassical D3 brane

We consider the refined Schur superconformal index of 4d N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 U(N) SYM and the first term of its giant-graviton expansion, first predicted in arXiv:2001.11667 using indirect superconformal algebra considerations and analytic continuation of fugacities. This correction is the leading non-perturbative correction to the index at large N and we reproduce it from the semiclassical partition function of quantum D3 brane wrapped on S1 × S3 in a twisted modification of the AdS5 × S5 string background, depending on the index R-symmetry fugacity. Our calculation does not exploit directly supersymmetry. It is based on the determination of the partition function of the various bosonic and fermionic fluctuations on the wrapped brane whose action is conformal with specific constant holonomies along thermal cycle. We show how those partition functions may be obtained by adapting the operator counting method of Cardy to the twisted background.

A mathematical analysis of the adiabatic Dyson equation from time-dependent density functional theory

In this article, we analyse the Dyson equation for the density–density response function (DDRF) that plays a central role in linear response time-dependent density functional theory (LR-TDDFT). First, we present a functional analytic setting that allows for a unified treatment of the Dyson equation with general adiabatic approximations for discrete (finite and infinite) and continuum systems. In this setting, we derive a representation formula for the solution of the Dyson equation in terms of an operator version of the Casida matrix. While the Casida matrix is well-known in the physics literature, its general formulation as an (unbounded) operator in the N-body wavefunction space appears to be new. Moreover, we derive several consequences of the solution formula obtained here; in particular, we discuss the stability of the solution and characterise the maximal meromorphic extension of its Fourier transform. We then show that for adiabatic approximations satisfying a suitable compactness condition, the maximal domains of meromorphic continuation of the initial DDRF and the solution of the Dyson equation are the same. The results derived here apply to widely used adiabatic approximations such as (but not limited to) the random phase approximation and the adiabatic local density approximation. In particular, these results show that neither of these approximations can shift the ionisation threshold of the Kohn–Sham system.

Crossing beyond scattering amplitudes

We find that different asymptotic measurements in quantum field theory can be related to one another through new versions of crossing symmetry. Assuming analyticity, we conjecture generalized crossing relations for multi-particle processes and the corresponding paths of analytic continuation. We prove them to all multiplicity at tree-level in quantum field theory and string theory. We illustrate how to practically perform analytic continuations on loop-level examples using different methods, including unitarity cuts and differential equations. We study the extent to which anomalous thresholds away from the usual physical region can cause an analytic obstruction to crossing when massless particles are involved. In an appendix, we review and streamline historical proofs of four-particle crossing symmetry in gapped theories.

K + N Elastic Scatterings for Estimation of the In-Medium Quark Condensate with Strange Quarks

Abstract We revisit the low-energy K+N elastic scatterings in the context of the in-medium quark condensate with strange quarks. The chiral ward identity connects the in-medium quark condensate to the soft limit value of the pseudoscalar correlation function evaluated in nuclear matter. The in-medium correlation function of the pseudoscalar fields with strangeness describes in-medium kaon propagation and is obtained by kaon–nucleon scattering amplitudes in the low-density approximation. We construct the kaon–nucleon scattering amplitudes in chiral perturbation theory up to the next-to-leading order and add some terms of the next-to-next-to-leading order with the strange quark mass to improve expansion of the strange quark sector. We also consider the effect of a possible broad resonance state around Plab = 600 MeV/c for I = 0 reported in the previous study. The low-energy constants are determined by existing K+N scattering data. We obtain good reproduction of the K+p scattering amplitude by chiral perturbation theory, while the description of the KN amplitude with I = 0 is not so satisfactory due to the lack of low-energy data. Performing analytic continuation of the scattering amplitudes obtained by chiral perturbation theory to the soft limit, we estimate the in-medium strange quark condensate.

Further remarks on de Sitter space, extremal surfaces, and time entanglement

We develop further the investigations by Narayan [de Sitter space, extremal surfaces, and time entanglement, ] on de Sitter space, extremal surfaces, and time entanglement. We discuss the no-boundary de Sitter extremal surface areas as certain analytic continuations from AdS while also amounting to space-time rotations. The structure of the extremal surfaces suggests a geometric picture of the time-entanglement or pseudo-entanglement wedge. We also study some entropy relations for multiple subregions. The analytic continuation suggests a heuristic Lewkowycz-Maldacena formulation of the extremal surface areas. In the bulk, this is now a replica formulation on the Wavefunction, which suggests interpretation as pseudo-entropy. Finally, we also discuss aspects of future-past entangled states and time evolution. Published by the American Physical Society 2024

Meromorphic Continuation of Global Zeta Function for Number Fields

Meromorphic Continuation of Global Zeta Function for Number Fields

Non-Markovian dynamics of time-fractional open quantum systems

The Time-Fractional Schrödinger Equation (TFSE) is suitable to describe the non-Markovian dynamics of a quantum system exposed to an environment, which is instructive for understanding and characterizing the time behavior of actual physical systems. Three popular TFSEs, namely Naber’s TFSE I, Naber’s TFSE II, and XGF’s TFSE, have been introduced in the form of ∂β∂tβ with β∈0,1. However, they suffer from the drawbacks. In their respective descriptions, the total probability for finding a particle in the entire space equals one only when β=1, implying that time-fractional quantum mechanics violates quantum mechanical probability conservation. By applying the three TFSEs to a basic single-qubit open system model, we discover that in describing the non-Markovian evolution of the system, the three TFSEs are limited by their own applicable ranges of β. This indicates that the three TFSEs cannot effectively describe the non-Markovian quantum dynamics. To address these issues, we introduce a well-performed TFSE by constructing a new analytic continuation of time without iβ combined with the conformable fractional derivative. In its description, for one thing, the total probability for finding the system in any single-qubit state equals one when β∈0,1. For another, the system evolves correctly in the non-Markovian manner at all values of β. Furthermore, we study the performances of the four TFSEs applying to a two-qubit open system model and show that our TFSE still possesses the above two advantages compared with the other three TFSEs.

Time operators of harmonic oscillators and their representations

A time operator T̂ϵ of the one-dimensional harmonic oscillator ĥϵ=12(p2+ϵq2) is rigorously constructed. It is formally expressed as T̂ϵ=121ϵ(arctan(ϵt̂0)+arctan(ϵt̂1)) with t̂0=p−1q and t̂1=qp−1. It is shown that the canonical commutation relation [hϵ,T̂ϵ]=−i1 holds true on a dense domain in the sense of sesqui-linear forms, and the limit of T̂ϵ as ϵ → 0 is shown. Finally a matrix representation of T̂ϵ and its analytic continuation are given.