Schwinger Equation Research Articles

Tubings, chord diagrams, and Dyson–Schwinger equations

AbstractWe give series solutions to single insertion place propagator‐type systems of Dyson–Schwinger equations using binary tubings of rooted trees. These solutions are combinatorially transparent in the sense that each tubing has a straightforward contribution. The Dyson–Schwinger equations solved here are more general than those previously solved by chord diagram techniques, including systems and noninteger values of the insertion parameter . We remark on interesting combinatorial connections and properties.

Lee–Yang edge singularities in QCD via the Dyson–Schwinger equations

We take the Dyson–Schwinger Equation approach of QCD for the quark propagator at complex chemical potential to study the QCD phase transition. The phase transition line of the (2+1)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$(2+1)$$\\end{document}-flavor QCD matter in the imaginary chemical potential region is computed via a simplified truncation scheme, whose curvature is found to be consistent with the one at real chemical potential. Moreover, the computation in the complex chemical potential plane allows us to determine the location of the Lee–Yang edge singularities. We show explicitly that the critical end point coincides with the Lee–Yang edge singularities on the real μB\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mu _{B} $$\\end{document} axis. We also investigate the scaling behavior of the singularities and discuss the possibility of extrapolating the CEP from a certain range of chemical potential.

On the conformal limit of a QED-inspired model

A conformal invariant QED-inspired model is solved for a general covariant linear gauge using the Dyson–Schwinger equations for the propagators assuming a pure vector like interaction. The leading corrections to the asymptotic solutions and the exponents, that characterize the corrections to each of the two fermion propagator functions, are computed as a function of the coupling and gauge fixing parameter ξ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\xi $$\\end{document}. For the scalar component of the fermion propagator our findings generalizes for linear covariant gauges previous results found in the literature and reproduce the outcome of the perturbative analysis of quenched QED in the Landau gauge. Our solution for the exponent associated with vector component of the fermion propagator is new and, in the weak coupling regime, agrees with the estimation based on the perturbative analysis of quenched QED. Of the two critical exponents describing the conformal limit of the vector interaction, one of them is, in QED, associated with the regime where chiral symmetry is broken dynamically, which demands one mass scale, namely the Miransky scaling. A second mass scale has to be introduced at larger coupling constants and is associated with a change on the nature of the fermion wave function. This provides one example, that it is possible to find two interwoven cycles in Quantum Field Theory, albeit in a truncated framework, as it is known in the quantum few-body problem in the limit of a zero-range interaction.

Confronting the Lippmann–Schwinger equation and the [formula omitted] method for coupled-wave separable potentials

We study a family of separable potentials with and without added contact interactions by solving the associated Lippmann–Schwinger equation with two coupled partial waves. The matching of the resulting amplitude matrix with the effective-range expansion is studied in detail. When a counterterm is included in the potential we also carefully discuss its renormalization. Next, we use the matrix N/D method and study whether the amplitude matrices from the potentials considered admit an N/D representation in matrix form. As a novel result we show that it is typically not possible to find such matrix representation for the coupled partial-wave case. However, a separate N/D representation for each coupled partial wave — a valid option known in the literature — is explicitly implemented and numerically solved in cases where the matrix N/D method is unavailable.

Flavor mixing and solution structures in Dyson–Schwinger equations for a two-flavor system

We solved the Dyson–Schwinger (DS) equations for a two-flavor system with symmetry to study its flavor mixing effects. Initially, we employed the point interaction model and bare vertex approximation to reveal the structure of the solutions. Using the point interaction model, the DS equations can be solved analytically, and we found that these solutions can be classified into three groups, each forming an ellipse. These solutions exhibit SO(2) symmetry, while the original SU(2) symmetry at the Lagrangian level is dynamically broken to SO(2), corresponding to the emergence of flavor mixing effects. However, this flavor mixing effect does not manifest in the final physical state. By utilizing the system’s SO(2) symmetry, we can diagonalize the propagators of the DS equations, eliminating the flavor mixing effect but causing the originally degenerate masses at the Lagrangian level to split. These mass eigenstates have identical quantum numbers but different masses. If we can correspond these to quark particles of different generations, we can explain why the three generations of quarks have different masses and obtain the corresponding quark mass spectrum. Finally, we provide the corresponding numerical results using a more realistic interaction model.

Experimental study of the solutions of the Lippmann–Schwinger equation for an elliptical billiard with an intense laser beam

This work presents an experimental study to emulate solutions of Lippmann-Schwinger equations for scattering by an elliptical billiard. The probability density patterns are recovered in the intensity of a laser beam scattered by a Spatial Light Modulator (SLM). This result opens an interesting perspective to explore scattering solutions, mainly when physical distributions with specific properties are required. Optical simulations performed by experiments with intense laser beams have been applied for many quantum systems and in this work, we expand this application.

FFT-based homogenisation for efficient concurrent multiscale modelling of thin plate structures

AbstractThis paper presents a novel approach to concurrent multiscale analysis, where structures are formulated at both microscopic and macro levels for simulation purposes. The proposed method employs a plate model to formulate structures at both scales, and homogenisation is performed using the FFT-based approach, offering higher efficiency compared to conventional methods. Additionally, the macroscopic tangent operator of the microscopic model is derived through an algorithmically consistent process within the FFT-based framework, incorporating the application of Lippman–Schwinger equations as outlined in this work. The effectiveness of the proposed method is demonstrated through case studies in real simulations, revealing comparable results to traditional multiscale schemes in addressing multiscale thin plate structures. Importantly, the method significantly reduces computing time and memory usage, attributed to the efficiency of plate modelling and the FFT-based homogenisation strategy.

Computation of Green’s Function in a Strongly Heterogeneous Medium Using the Lippmann–Schwinger Equation: A Generalized Successive Over-Relaxion plus Preconditioning Scheme

The computation of Green’s function is a basic and time-consuming task in realizing seismic imaging using integral operators because the function is the kernel of the integral operators and because every image point functions as the source point of Green’s function. If the perturbation theory is used, the problem of the computation of Green’s function can be transformed into one of solving the Lippmann–Schwinger (L–S) equation. However, if the velocity model under consideration has large scale and strong heterogeneity, solving the L–S equation may become difficult because only numerical or successive approximate (iterative) methods can be used in this case. In the literature, one of these methods is the generalized successive over-relaxation (GSOR) iterative method, which can effectively solve the L–S equation and obtain the desired convergent iterative series. However, the GSOR iterative method may encounter slow convergence when calculating the high-frequency Green’s function. In this paper, we propose a new scheme that utilizes the GSOR iterative with a precondition to solve the complex wavenumber L–S equation in a slightly attenuated medium. The complex wavenumber with imaginary components localizes the energy of the background Green’s function and reduces its singularity by enabling exponential decay. Introduction of the preconditioning operator can further improve the convergence speed of the GSOR iterative series. Then, we provide a preconditioned generalized successive over-relaxation (pre-GSOR) iterative format. Our numerical results show that if an appropriate damping factor and a proper preconditioning operator are selected, the method presented here outperforms the GSOR iterative for the real wavenumber L–S equation in terms of the convergence speed, accuracy, and adaptation to high frequencies.

An FFT-based adaptive polarization method for infinitely contrasted media with guaranteed convergence

We propose an FFT-based iterative algorithm for solving the Lippmann–Schwinger equation in the context of periodic homogenization of infinitely contrasted linear elastic composites. Our work initially reformulates the Moulinec–Suquet, Eyre–Milton and Monchiet–Bonnet schemes using a residual formulation. Subsequently, we introduce an enhanced scheme, termed Adaptive Eyre–Milton (AEM), as a natural extension of the EM scheme where we optimize a relaxation parameter to minimize the residual. We demonstrate the unconditional linear convergence of the AEM scheme, regardless of initialization and the chosen reference material. The paper further extends the AEM scheme to handle composites with both pores and infinitely rigid inclusions. Practical implementation aspects and illustrative applications in two- and three-dimensional settings are discussed, highlighting the efficiency of the proposed AEM scheme. We particularly emphasize the scheme’s robustness for materials with infinitely large contrasts in elastic properties.

Fast Fourier transform method in peridynamic micromechanics of composites

We consider a static linear bond–based peridynamic (proposed by Silling, see J. Mech. Phys. Solids 2000; 48:175–209) composite materials (CMs) of a periodic structure. In the framework of the second background of micromechanics (also called computational analytical micromechanics), one proved that local micromechanics (LM) and peridynamic micromechanics (PM) are formally analogous to each other for CM of both random and periodic structures. It allows a straightforward generalization of LM methods (including fast Fourier transform, FFT) to their PM counterparts. So, in the PM counterpart of the implicit periodic Lippmann–Schwinger (L-S) equation in LM, we have three convolution kernels corresponding to the properties of the matrix, inclusions, and interactive interface. Eshelby tensor in LM, depending on the inclusion shape, is replaced by PM counterparts depending on the shapes of inclusions, and the interaction interface (although the Eshelby concept of homogeneous eigenfields does not work in PM). The mentioned tensors are estimated once (as in LM) in a frequency domain (also by the FFT method). The possible incorrectness of FFT applications to PM is analyzed and corrected. The polarization schemes of the solution of the L-S equation in the Fourier space have one primary unknown variable (polarization), whereas the PM counterpart contains three primary ones estimated at each step, which are formally similar to the LM case. A description of the generalized basic scheme and the Krylov approach is presented. Computational complexities O(N log2 N) of the FFT methods are the same in both LM and PM.

Renormalisable Non-Local Quark–Gluon Interaction: Mass Gap, Chiral Symmetry Breaking and Scale Invariance

We derive a Nambu–Jona-Lasinio (NJL) model from a non-local gauge theory and show that it has confining properties at low energies. In particular, we present an extended approach to non-local QCD and a complete revision of the technique of Bender, Milton and Savage applied to non-local theories, providing a set of Dyson–Schwinger equations in differential form. In the local case, we obtain closed-form solutions in the simplest case of the scalar field and extend it to the Yang–Mills field. In general, for non-local theories, we use a perturbative technique and a Fourier series and show how higher-order harmonics are heavily damped due to the presence of the non-local factor. The spectrum of the theory is analysed for the non-local Yang–Mills sector and found to be in agreement with the local results on the lattice in the limit of the non-locality mass parameter running to infinity. In the non-local case, we confine ourselves to a non-locality mass that is sufficiently large compared to the mass scale arising from the integration of the Dyson–Schwinger equations. Such a choice results in good agreement, in the proper limit, with the spectrum of the local theory. We derive a gap equation for the fermions in the theory that gives some indication of quark confinement in the non-local NJL case as well. Confinement seems to be a rather ubiquitous effect that removes some degrees of freedom in the original action, favouring the appearance of new observable states, as seen, e.g., for quantum chromodynamics at lower energies.

Quantum scattering in a wavy circle

We study the non-relativistic quantum mechanical scattering problem of a plane-wave by a Chebyshev particle in two spatial dimensions. The Chebyshev particle is a wavy, or undulating, circle. We formulate the problem using the Lippmann–Schwinger (LS) equation and the Chebyshev particle as a boundary-wall potential. This potential is naturally symmetric as a function of the so-called undulation parameter, and we use these symmetries to solve the LS equation and calculate three important physical quantities: the differential cross-section, the total cross-section and the quantum refractive index. We analyze these quantities as functions of the wave number k, the coupling strength γ as well as the undulation and deformation parameters of the Chebyshev particle.

Analysis of Strong Coupling Constant with Machine Learning and Its Application

We investigate the nature of the strong coupling constant and related physics. Through the analysis of accumulated experimental data around the world, we employ the ability of machine learning to unravel its physical laws. The result of our efforts is a formula that captures the expansive panorama of the distribution of the strong coupling constant across the entire energy range. Importantly, this newly derived expression is very similar to the formula derived from the Dyson–Schwinger equations based on the framework of Yang–Mills theory. By introducing the Euler number e into the functional formula of the strong coupling constant at high energies, we successfully solve the puzzle of the infrared divergence, which allows for a seamless transition of the strong coupling constant from the perturbative to the non-perturbative energy regime. Moreover, the obtained ghost and gluon dressing function distribution results confirm that the obtained strong coupling constant formula can well describe the physical properties of the non-perturbed regime. In addition, we study the quantum-chromodynamics strong coupling constant result of the Bjorken sum rule Γ1p−n and the quark–quark static energy E 0(r), and find that the global energy scale can effectively interpret the experimental data. The present results shed light on the puzzling properties of quantum chromodynamics and the intricate interplay of strong coupling constants at both low and high energy scales.

Disassociation of a one-dimensional cold molecule via quantum scattering

Motivated by the recent experimental developments in ultracold molecules and atoms, we propose a simple theoretical model to address the disassociation, reflection, and transmission probability of a one-dimensional cold molecule via quantum scattering. First, we show the Born approximation results in the weak interaction regime. Then, by employing the Lippmann–Schwinger equation, we give the numerical solution and investigate the disassociation’s dependence on the injection momentum and the interaction strengths. We find that the maximum disassociation rate has a limit when increasing the interaction strengths and injection momentum. We expect that our model can be realized in experiments in the near future.

Fourier Transform of the Lippmann-Schwinger Equation: Solving Vectorial Electromagnetic Scattering by Arbitrary Shapes

In Electromagnetics, the field scattered by an ensemble of particles—of arbitrary size, shape, and material—can be obtained by solving the Lippmann–Schwinger equation. This singular vectorial integral equation is generally formulated in the direct space Rn (typically n=2 or n=3). In the article, we rigorously computed the Fourier transform of the vectorial Lippmann–Schwinger equation in the space of tempered distributions, S′(R3), splitting it in a singular and a regular contribution. One eventually obtains a simple equation for the scattered field in the Fourier space. This permits to draw an explicit link between the shape of the scatterer and the field through the Fourier Transform of the body indicator function. We compare our results with accurate calculations based on the T-matrix method and find a good agreement.

Binding properties of some alkali metal vapors: specifically lithium-7 (7Li), soduim-23 (23Na), and potassium-39 (39K)

In this work, the binding energy of three alkali metal molecules, namely, lithium-7 (7Li2), sodium-23 (23Na2), and potassium-39 (39K2), is calculated in both free space and the vapor phase. In free space, the Lippmann–Schwinger equation is solved for negative-definite eigenenergies using a highly effective matrix-inversion method with a symmetrized kernel. In the vapor, two medium effects are taken into account: the effective atomic mass, which is somewhat larger than the “bare” mass and therefore enhances the binding; and the effective binary interaction, with less overall attraction than the “bare” interaction. In free space, the binding energies of the respective molecules are, in units of (103 K), 5.304, 2.936, and 2.522. The corresponding experimental results are 6.150 × 103, 4.275 ×103, and 4.275 × 103 K. The results in the vapor phase are by and large somewhat lower than those in free space.

Neutrino Oscillations in Finite Time Path Out-of-Equilibrium Thermal Field Theory

We demonstrate that the Finite-Time-Path Field Theory is an adequate tool for calculating neutrino oscillations. We apply this theory using a mass-mixing Lagrangian which involves the correct Dirac spin and chirality structure and a Pontecorvo–Maki–Nakagawa–Sakata (PMNS)-like mixing matrix. The model is exactly solvable. The Dyson–Schwinger equations transform propagators of the input free (massless) flavor neutrinos into a linear combination of oscillating (massive) neutrinos. The results are consistent with the predictions of the PMNS matrix while allowing for extrapolation to early times.

New insights into the oscillations of the nucleon electromagnetic form factors

The electromagnetic form factors of the proton and the neutron in the timelike region are investigated. Electron–positron annihilation into antinucleon-nucleon (N¯N) pairs is treated in distorted wave Born approximation, including the final-state interaction in the N¯N system. The latter is obtained by a Lippmann–Schwinger equation for N¯N potentials derived within SU(3) chiral effective field theory. By fitting to the phase shifts and (differential) cross section data, a high quality description is achieved. With these amplitudes, the oscillations of the electromagnetic form factors of the proton and the neutron are studied. It is found that each of them can be described by two fractional oscillators. One is characterized as “overdamped” and dominates near the threshold, while the other is “underdamped” and plays an important role in the high-energy region. These two oscillators are essential to understand the distributions of polarized electric charges induced by hard photons for the nucleons.

Fermionic greybody factors and quasinormal modes of black holes in Kalb–Ramond gravity

In non-minimally coupled Kalb–Ramond (KR) fields, the vacuum expectation value (VEV) leads to spontaneous local Lorentz symmetry violation, resulting in static spherically symmetric solution. In this solution, fermionic greybody factors (GFs) and quasinormal modes (QNMs) are investigated. First, we study the wave dynamics of fermions using the Dirac and Rarita–Schwinger equations. We have derived the equations of effective potential for fermions with different spins. The fermionic GFs radiated by the black holes (BHs) are calculated using rigorous bounds as well as semi-analytic methods. The fermionic QNM of spin- 1/2 of the BHs in the KR gravity was analyzed using the sixth order WKB method. Furthermore, we discuss how Lorentz violating (LV) parameters impact GFs and QNMs. According to our results, the effective potential curve increases significantly when the LV parameter value γ is high. Thus, as the LV parameter γ increases, the bound of the GF and the QNM’s wave oscillation frequency will increase, while the decay rate will slow down.

An efficient iterative method for reconstructing the refractive index in complex domains from far field data

We present a computational method for reconstructing the refractive index of an unknown complex-shaped two-dimensional medium with embedded metal inclusions from its transverse magnetic electromagnetic scattering properties. We present a novel hybrid surface-volume integral equation that generalises the Lippmann–Schwinger equation to media containing embedded scatterers. Using this hybrid equation, we show that the Frechet derivative of the forward-mapping satisfies a particular inhomogeneous wave scattering problem. We solve the inhomogeneous wave scattering problem using a novel coupled FEM-BEM formulation. Our numerical scheme is based on a thin plate spline ansatz for the refractive index, which can be constructed using a relatively small number of control points, and can be efficiently constructed and evaluated even for the complex-shaped multiply-connected domains of interest. Numerical experiments demonstrate the effectiveness of our method by reconstructing several challenging media.