Vertex Corrections Research Articles

Discontinuous Galerkin solutions for elliptic problems on polygonal grids using arbitrary-order Bernstein-Bézier functions

In this paper, we present a symmetric interior penalty discontinuous Galerkin finite element discretization for the numerical solution of second-order elliptic partial differential equations on general polygonal meshes using a reduced-space polygonal Bernstein-Bézier functions. They form a space of interpolatory functions (vice local polynomial functions), which satisfy the Lagrange property at their interpolation points. This can enable the use of efficient DGFEM physics-based diffusion preconditioners for the first-order linear transport equation on arbitrary polygons in the computationally-challenging asymptotic diffusion limit. On a polygonal element with n vertices and polynomial degree p, there are n vertex functions, (p−1) functions per edge, and (p−1)(p−2)/2 interior functions that span the {xayb}(a+b)≤p space of bivariate monomials. Numerical experiments highlighting the performance of the proposed method are presented through uniform and adaptive h, p, and hp refinement strategies. Appropriate error rates are observed including hp adaptivity yielding exponential convergence in the presence of singularities in the solution.

Are multi-quasiparticle interactions important in molecular ionization?

Photo-emission spectroscopy directly probes individual electronic states, ranging from single excitations to high-energy satellites, which simultaneously represent multiple quasiparticles (QPs) and encode information about electronic correlation. The first-principles description of the spectra requires an efficient and accurate treatment of all many-body effects. This is especially challenging for inner valence excitations where the single QP picture breaks down. Here, we provide the full valence spectra of small closed-shell molecules, exploring the independent and interacting quasiparticle regimes, computed with the fully correlated adaptive sampling configuration interaction method. We critically compare these results to calculations with the many-body perturbation theory, based on the GW and vertex corrected GWΓ approaches. The latter explicitly accounts for two-QP quantum interactions, which have often been neglected. We demonstrate that for molecular systems, the vertex correction universally improves the theoretical spectra, and it is crucial for the accurate prediction of QPs as well as capturing the rich satellite structures of high-energy excitations. GWΓ offers a unified description across all relevant energy scales. Our results suggest that the multi-QP regime corresponds to dynamical correlations, which can be described via perturbation theory.

Diagrammatic study of optical excitations in correlated systems

The optical conductivity contains relevant information on the properties of correlated electron systems. In infinite dimensions, where dynamical mean field theory becomes exact, vertex corrections can be neglected and the conductivity computed from particle-hole bubbles. An interesting question concerns the nature and effect of the most relevant vertex corrections in finite-dimensional systems. A recent numerical study showed that the dominant vertex correction near an ordering instability with wave vector {\pi} comes from a vertical ladder, analogous to the Maki-Thompson diagram. Since the RPA version of this ladder diagram, dubbed {\pi}-ton, can be easily evaluated, this suggests a simple procedure for incorporating antiferromagnetic or charge density wave fluctuations into dynamical mean field estimates of the optical conductivity and related susceptibilities. We implement this procedure for the half-filled Hubbard model, considering the {\pi}-ton and a double-ladder extension of the {\pi}-ton, and reveal the spectral signatures of these vertex corrections.

Anisotropic magnetotransport in tilted magnetic topological semimetals

We systematically investigate the anisotropic magneto-conductivity and planar Hall effect in tilted magnetic topological semimetals in the frame of Kubo formula by considering the vertex correction of velocity. The nonzero anisotropic magneto-conductivity is due to the intrinsic magnetization by magnetic doping rather than the external magnetic field previously studied in literatures. In the scenario, tilt, which is common in band structure, plays a key role in this anisotropic magneto-conductivity. This anisotropic magneto-conductivity completely originates from the band anisotropy. The vertex correction only amends the magnitude, but does not change the qualitative behavior of the longitudinal conductivity. However, the planar Hall conductivity is always zero for tilt along both x and y directions in vivid contrast to the case arising from an external magnetic field. Our work offers a route to understand the anisotropic magneto-conductivity effect in tilted Weyl semimetals with Zeeman field. It will be helpful to get a deeper understanding of the transport characteristic in tilted topological semimetals.

Towards exact solutions for the superconducting Tc induced by electron-phonon interaction

Electron-phonon interaction plays an important role in metals and can lead to superconductivity and other instabilities. Previous theoretical studies on superconductivity are largely based on the Migdal-Eliashberg theory, which neglects all the vertex corrections to electron-phonon coupling and breaks down in many unconventional superconductors. Here, we go beyond the Migdal-Eliashberg approximation and develop a nonperturbative Dyson-Schwinger equation approach to deal with the superconducting transition. Remarkably, we take into account all the vertex corrections by solving two coupled Ward-Takahashi identities derived from two global U(1) symmetries and rigorously prove that the fully renormalized electron propagator satisfies a self-closed integral equation that is directly amenable to numerical computations. Our approach works equally well in the weak and strong coupling regimes and provides an efficient method to determine superconducting $T_c$ and other quantities. As an application, our approach is used to investigate the high-$T_c$ superconductivity in one-unit-cell FeSe/SrTiO$_3$.

Probing new U(1) gauge symmetries via exotic Z \u2192 Z\u2032\u03b3 decays

New U(1) gauge theories involving Standard Model (SM) fermions typically require additional electroweak fermions for anomaly cancellation. We study the non-decoupling properties of these new fermions, called anomalons, in the Z − Z′ − γ vertex function, reviewing the connection between the full model and the effective Wess-Zumino operator. We calculate the exotic Z → Z′γ decay width in U(1)B−L and U(1)B models, where B and L denote the SM baryon and lepton number symmetries. For U(1)B−L gauge symmetry, each generation of SM fermions is anomaly free and the exotic Z → {Z}_{BL}^{prime}gamma decay width is entirely induced by intragenerational mass splittings. In contrast, for U(1)B gauge symmetry, the existence of two distinct sources of chiral symmetry breaking enables a heavy, anomaly-free set of fermions to have an irreducible contribution to the Z → {Z}_B^{prime}gamma decay width. We show that the current LEP limits on the exotic Z → {Z}_B^{prime}gamma decay are weaker than previously estimated, and low-mass {Z}_B^{prime } dijet resonance searches are currently more constraining. We present a summary of the current collider bounds on U(1)B and a projection for a TeraZ factory on the Z → {Z}_B^{prime}gamma exotic decay, and emphasize how the Z → Z′γ decay is emblematic of new anomalous U(1) gauge symmetries.

Tiling with triangles: parquet and GWγ methods unified

The parquet formalism and Hedin's $GW\gamma$ approach are unified into a single theory of vertex corrections, corresponding to an exact reformulation of the parquet equations in terms of boson exchange. The method has no drawbacks compared to previous parquet solvers but has the significant advantage that the vertex functions decay quickly with frequencies and with respect to distances in real space. These properties coincide with the respective separation of the length and energy scales of the two-particle correlations into long/short-ranged and high/low-energetic.

Weak Localization Correction to Linear Absorption in Conventional Superconductors

The weak localization effect on a linear absorption spectrum is investigated for disordered s-wave superconductors. The vertex correction is incorporated into the response function in a way that is consistent with the weak localization correction to a one-particle spectrum. The effect of interactions between electrons enhanced by their diffusive motion makes a major contribution to the correction term. Numerical calculations show that the weak localization effect suppresses the absorption by the excitation across the gap to a greater extent than that by thermal excitation, which is a similar tendency to the results of experiments. The ratio of the linear absorption with vertex corrections to that without the corrections shows a large variation at the frequency that is around twice the energy of the superconducting gap. This behavior represents a relationship between the weak localization effect on the linear absorption and that on the one-particle spectrum.

Dyson’s Equations for Quantum Gravity in the Hartree–Fock Approximation

Unlike scalar and gauge field theories in four dimensions, gravity is not perturbatively renormalizable and as a result perturbation theory is badly divergent. Often the method of choice for investigating nonperturbative effects has been the lattice formulation, and in the case of gravity the Regge–Wheeler lattice path integral lends itself well for that purpose. Nevertheless, lattice methods ultimately rely on extensive numerical calculations, leaving a desire for alternate methods that can be pursued analytically. In this work, we outline the Hartree–Fock approximation to quantum gravity, along lines which are analogous to what is done for scalar fields and gauge theories. The starting point is Dyson’s equations, a closed set of integral equations which relate various physical amplitudes involving graviton propagators, vertex functions, and proper self-energies. Such equations are in general difficult to solve, and as a result they are not very useful in practice, but nevertheless provide a basis for subsequent approximations. This is where the Hartree–Fock approximation comes in, whereby lowest order diagrams get partially dressed by the use of fully interacting Green’s function and self-energies, which then lead to a set of self-consistent integral equations. The resulting nonlinear equations for the graviton self-energy show some remarkable features that clearly distinguish it from the scalar and gauge theory cases. Specifically, for quantum gravity one finds a nontrivial ultraviolet fixed point in Newton’s constant G for spacetime dimensions greater than two, and nontrivial scaling dimensions between d=2 and d=4, above which one obtains Gaussian exponents. In addition, the Hartree–Fock approximation gives an explicit analytic expression for the renormalization group running of Newton’s constant, suggesting gravitational antiscreening with Newton’s constant slowly increasing on cosmological scales.

The full vertex functions constrained by the longitudinal and transverse Ward–Takahashi identities in massless QED3 theory

In this paper, we find apart from the Ward–Takahashi (WT) identity, the identity between gamma matrices can also constrain the vertex functions in low-dimensional gauge theories. In (1 + 1) dimensions, the identity between gamma matrices gives the identity between vector and axial-vector vertex functions while in (2 + 1) dimensions it leads to the identity between vector and tensor vertex functions. Then, we derive the expressions of the full scalar, vector and tensor vertex functions in (2 + 1) dimensions Quantum Electrodynamics (QED3) by using the longitudinal and transverse WT identities for vector and tensor currents. Furthermore, we find that in the chiral limit with zero fermion masses, the contribution of Wilson line in full vector vertex function is eliminated and the full vector vertex function is strictly expressed in terms of the fermion propagators when using the identity between vector and tensor vertex functions to further constraint the vertex functions.

Schwarzschild-Tangherlini metric from scattering amplitudes in various dimensions

We derive the static Schwarzschild-Tangherlini metric by extracting the classical contributions from the multiloop vertex functions of a graviton emitted from a massive scalar field. At each loop order the classical contribution is proportional to a unique master integral given by the massless sunset integral. By computing the scattering amplitudes up to three-loop order in general dimension, we explicitly derive the expansion of the metric up to the fourth post-Minkowskian order $O({G}_{N}^{4})$ in four, five and six dimensions. There are ultraviolet divergences that are cancelled with the introduction of higher-derivative nonminimal couplings. The standard Schwarzschild-Tangherlini is recovered by absorbing their effects by an appropriate coordinate transformation induced from the de Donder gauge condition.

QED effects on the spin asymmetry in elastic $$^{12}$$C($$e,e'$$) collisions

The influence of quantum electrodynamic (QED) effects on the Sherman function S for the scattering of polarized electrons from a $$^{12}$$ C nucleus is investigated. While vacuum polarization is treated to all orders via the Uehling potential, the vertex correction, self energy, soft bremsstrahlung and dispersion corrections are accounted for to lowest order within the plane-wave Born approximation. The latter are estimated with the help of a closure approximation. The additional emission of hard bremsstrahlung is calculated by means of the Sommerfeld–Maue theory. Photons up to 4 MeV are included, corresponding to the detector resolution in future experiments. In the collision energy regime considered, 20–180 MeV, the combined QED effects lead to considerable changes in S, amounting to 5–40% at the larger scattering angles.

Real time ultrasoft fermion self energy at next to-leading order in hot QED

Subsequent studies of the behavior of the gluon and quark damping rates in the imaginary-time formalism have indicated that there are difficulties in the infrared sector [1, 2, 3, 4, 5, 6, 7]. To look further into the infrared behavior, we propose to calculate the next-to-leading order dispersion relations for slow-moving Fermions at high-temperature quantum electrodynamics (QED) in real-time formalism. We determine a compact analytic expression for the complete next-to-leading contribution to the retarded fermion self-energy with ultrasoft momentum in the framework of hard-thermal-loop (HTL)-summed perturbation of massless QED at high temperature. The calculation is done using real-time formalism. The next-to-leading order fermion self-energy is written in terms of three and four HTL-dressed vertex functions. The real part and the opposite of the imaginary part of the retarded fermion self-energy are related to the next-to-leading order contributions of energy and damping rate respectively.

Fluctuation diagnostic of the nodal/antinodal dichotomy in the Hubbard model at weak coupling: A parquet dual fermion approach

We apply the boson exchange parquet solver for dual fermions to the half-filled Hubbard model on a square lattice at small interaction. Our results establish that, in this regime, nonlocal vertex corrections play an important role in the formation of the pseudogap. Namely, in comparison to the simpler ladder approximation, these additional vertex corrections included in the parquet equations enhance the coupling of spin fluctuations with the quasiparticles. The pseudogap thus opens already at a higher temperature, in quantitative agreement with the numerically exact diagrammatic Monte Carlo. The representation of the parquet diagrams in terms of boson exchange facilitates large lattice sizes and gives rise to an unbiased fluctuation diagnostic of the self-energy, which does not rely on the Fierz ambiguity. The fluctuation diagnostic implies that nodal and antinodal fermions are affected equally by spin fluctuations with the exact commensurate nesting vector $(\pi,\pi)$. However, the antinode couples more efficiently to incommensurate fluctuations than the node, leading to the nodal/antinodal dichotomy. We corroborate this finding in terms of a spin-fermion-like calculation.

Bulk spin conductivity of three-dimensional topological insulators

We study the spin conductivity of the bulk states of three-dimensional topological insulators within Kubo formalism. Spin Hall effect is the generation of the spin current that is perpendicular to the applied voltage. In the case of a three-dimensional topological insulator, applied voltage along x direction generates transverse spin currents along y and z directions with comparable values. We found that finite non-universal value of the spin conductivity exists in the gapped region due to the inversion of bands. Contribution to the spin conductivity from the vertex corrections enhances the spin conductivity from the filled states. These findings explain large spin conductivity that has been observed in topological insulators.

Boson-exchange parquet solver for dual fermions

We present and implement a parquet approximation within the dual-fermion formalism based on a partial bosonization of the dual vertex function which substantially reduces the computational cost of the calculation. The method relies on splitting the vertex exactly into single-boson exchange contributions and a residual four-fermion vertex, which physically embody respectively long-range and short-range spatial correlations. After recasting the parquet equations in terms of the residual vertex, these are solved using the truncated unity method of Eckhardt et al. [Phys. Rev. B 101, 155104 (2020)], which allows for a rapid convergence with the number of form factors in different regimes. While our numerical treatment of the parquet equations can be restricted to only a few Matsubara frequencies, reminiscent of Astretsov et al. [Phys. Rev. B 101, 075109 (2020)], the one- and two-particle spectral information is fully retained. In applications to the two-dimensional Hubbard model the method agrees quantitatively with a stochastic summation of diagrams over a wide range of parameters.

Anomalous Hall effect in single-band chiral superconductors from impurity superlattices

Unlike anomalous quantum Hall insulators, clean single-band chiral superconductors do not exhibit intrinsic Hall effect at the one-loop approximation. Finite ac Hall conductance was found to emerge beyond one-loop, such as with vertex corrections associated with extrinsic random impurity scatterings. In this paper, we investigate the effect of impurities embedded in single-band chiral superconductors in a superlattice pattern, instead of in random distributions. The impurity-induced Bogoliubov quasiparticle bound states hybridize to form subgap bands, constituting an emergent low-energy effective theory whose Hall effect can be studied with ease. We demonstrate that the occurrence of the Hall effect depends on the superlattice geometry and on the parity of the chiral pairing. In particular, due to the mixed particle-hole character of the subgap states, the Hall conductance may arise at the one-loop level of the current-current correlator in our effective model. Our theory provides a new insight into the impurity-induced anomalous Hall effect in chiral superconductors.

On degree sequence optimization

We consider the problem of finding a subgraph of a given graph maximizing a given function evaluated at its degree sequence. While it is intractable already for convex functions, we show it is polynomial time solvable for convex multi-criteria objectives. We also consider a colored extension of the problem with separable objectives, which includes the notorious exact matching problem as a special case, and show that it is polynomial time solvable on graphs of bounded tree-depth for any vertex functions.

One-loop vertex correction in a plane wave

We compute the general expression of the one-loop vertex correction in an arbitrary plane-wave background field for the case of two on-shell external electrons and an off-shell external photon. The properties of the vertex corrections under gauge transformations of the plane-wave background field and of the radiation field are studied. Concerning the divergences of the vertex correction, the infrared one is cured by assigning a finite mass to the photon, whereas the ultraviolet one is shown to be renormalized exactly as in vacuum. Finally, the corresponding expression of the vertex correction within the locally-constant crossed field is also derived and the high-field asymptotic is shown to scale according to the Ritus-Narozhny conjecture.

Kinetics and luminescence of the excitations of a nonequilibrium polariton condensate

We extend the description of the nonequilibrium Gross-Pitaevskii equation for microcavity polaritons by including not only the kinetics of the nonresonantly pumped exciton reservoir but also that of the collective excitations---called bogolons---above the condensate. Quasiequilibrium Hartree-Fock-Bogoliubov-Popov self-energies are used to describe the condensate and the bogolons. With nonequilibrium Green functions, the quantum Boltzmann equation of the bogolons in a frame moving with the condensate, together with the coupled Gross-Pitaevskii equation can be derived. The resulting scattering rates contain vertex corrections in form of symmetrized products of the $u$ and $v$ pairing functions of all involved scattering states. Corresponding results have been derived for ultracold atoms by A. Griffin, T. Nikuni, and E. Zaremba [Bose-Condensed Gases at Finite Temperatures, 1st ed. (Cambridge University Press, Cambridge, New York, 2009)] Furthermore, we show that the results of the quantum kinetics also determine the luminescence spectrum including the ghost branch. The kinetics is solved numerically, the resulting condensate, the bogolon distribution, and the luminescence are calculated temporally and spatially resolved after the switch-on of a $\mathit{cw}$ nonresonant pump beam with a Gaussian profile. It is studied in detail to which extent the bogolons approach a local thermal equilibrium.