Wichmann-Kroll vacuum polarization density in a finite Gaussian basis set

• URL: http://arxiv.org/abs/2512.16569v1
• Date: Thu, 18 Dec 2025 14:11:30 GMT
• Title: Wichmann-Kroll vacuum polarization density in a finite Gaussian basis set
• Authors: Ryan Benazzouk, Maen Salman, Trond Saue,
• Abstract summary: This work further develops the calculation of QED effects in a finite Gaussian basis.<n>We focus on the non-linear $(Z)nge 3$ contribution to the vacuum polarization density.<n>We then report a strategy for computing the energy shift with sufficient precision to enable a sensible extrapolation of the partial-wave expansion.
• Score: 0.0
• License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/
• Abstract: This work further develops the calculation of QED effects in a finite Gaussian basis. We focus on the non-linear $α(Zα)^{n\ge 3}$ contribution to the vacuum polarization density, computing the energy shift of 1s$_{1/2}$ states of hydrogen-like ions. Our goal is to improve the numerical computations to achieve a precision comparable to that of Green's function methods reported in the literature. To do so, an analytic expression for the linear contribution to the vacuum polarization density is derived using Riesz projectors. Alternative formulations of the vacuum polarization density and their relation is discussed. The convergence of the finite Gaussian basis scheme is investigated, and the numerical difficulties that arise are characterized. In particular, an error analysis is performed to assess the method's robustness to numerical noise. We then report a strategy for computing the energy shift with sufficient precision to enable a sensible extrapolation of the partial-wave expansion. A key feature of the procedure is the use of even-tempered basis sets, allowing for an extrapolation towards the complete basis set limit.

Related papers

• Determination of Particle-Size Distributions from Light-Scattering Measurement Using Constrained Gaussian Process Regression [7.5715377940231114]
  We propose a novel methodology for robustly estimating particle size distributions from optical scattering measurements.<n>The proposed constrained Gaussian process regression framework accurately reconstructs particle size distributions.
  arXiv  Detail & Related papers  (2025-07-04T17:56:16Z)
• Vacuum polarization and Wichmann-Kroll correction in the finite basis set approximation [0.0]
  We study the convergence of the finite basis set method with different types and sizes of basis sets. We consider several heavy hydrogen-like ions and evaluate the vacuum polarization correction for $S$ and $P$ electron orbitals.
  arXiv  Detail & Related papers  (2024-06-10T19:36:11Z)
• A score-based particle method for homogeneous Landau equation [7.600098227248821]
  We propose a novel score-based particle method for solving the Landau equation in plasmas.<n>Our primary innovation lies in recognizing that this nonlinearity is in the form of the score function.<n>We provide a theoretical estimate by demonstrating that the KL divergence between our approximation and the true solution can be effectively controlled by the score-matching loss.
  arXiv  Detail & Related papers  (2024-05-08T16:22:47Z)
• Noise-Free Sampling Algorithms via Regularized Wasserstein Proximals [3.4240632942024685]
  We consider the problem of sampling from a distribution governed by a potential function. This work proposes an explicit score based MCMC method that is deterministic, resulting in a deterministic evolution for particles.
  arXiv  Detail & Related papers  (2023-08-28T23:51:33Z)
• Sobolev Space Regularised Pre Density Models [51.558848491038916]
  We propose a new approach to non-parametric density estimation that is based on regularizing a Sobolev norm of the density. This method is statistically consistent, and makes the inductive validation model clear and consistent.
  arXiv  Detail & Related papers  (2023-07-25T18:47:53Z)
• A mathematical analysis of the adiabatic Dyson equation from time-dependent density functional theory [0.0]
  We analyze the Dyson equation for the density-density response function (DDRF) that plays a central role in linear response time-dependent density functional theory. We derive a representation formula for the solution of the Dyson equation in terms of an operator version of the Casida matrix. We show that for adiabatic approximations satisfying a suitable compactness condition, the maximal domains of meromorphic continuation of the initial density-density response function and the solution of the Dyson equation are the same.
  arXiv  Detail & Related papers  (2023-05-15T15:46:33Z)
• Calculating the many-potential vacuum polarization density of the Dirac equation in the finite-basis approximation [0.0]
  We propose an efficient and accurate computational method to evaluate the many-potential vacuum polarization density of hydrogen-like atoms. To prove the performance of our computational method, we choose to work with the one-electron $_,,,92238textU$ atom.
  arXiv  Detail & Related papers  (2023-04-18T14:23:06Z)
• Sampling with Mollified Interaction Energy Descent [57.00583139477843]
  We present a new optimization-based method for sampling called mollified interaction energy descent (MIED) MIED minimizes a new class of energies on probability measures called mollified interaction energies (MIEs) We show experimentally that for unconstrained sampling problems our algorithm performs on par with existing particle-based algorithms like SVGD.
  arXiv  Detail & Related papers  (2022-10-24T16:54:18Z)
• Mean-Square Analysis with An Application to Optimal Dimension Dependence of Langevin Monte Carlo [60.785586069299356]
  This work provides a general framework for the non-asymotic analysis of sampling error in 2-Wasserstein distance. Our theoretical analysis is further validated by numerical experiments.
  arXiv  Detail & Related papers  (2021-09-08T18:00:05Z)
• Large-Scale Wasserstein Gradient Flows [84.73670288608025]
  We introduce a scalable scheme to approximate Wasserstein gradient flows. Our approach relies on input neural networks (ICNNs) to discretize the JKO steps. As a result, we can sample from the measure at each step of the gradient diffusion and compute its density.
  arXiv  Detail & Related papers  (2021-06-01T19:21:48Z)
• Variational Transport: A Convergent Particle-BasedAlgorithm for Distributional Optimization [106.70006655990176]
  A distributional optimization problem arises widely in machine learning and statistics. We propose a novel particle-based algorithm, dubbed as variational transport, which approximately performs Wasserstein gradient descent. We prove that when the objective function satisfies a functional version of the Polyak-Lojasiewicz (PL) (Polyak, 1963) and smoothness conditions, variational transport converges linearly.
  arXiv  Detail & Related papers  (2020-12-21T18:33:13Z)
• Local optimization on pure Gaussian state manifolds [63.76263875368856]
  We exploit insights into the geometry of bosonic and fermionic Gaussian states to develop an efficient local optimization algorithm. The method is based on notions of descent gradient attuned to the local geometry. We use the presented methods to collect numerical and analytical evidence for the conjecture that Gaussian purifications are sufficient to compute the entanglement of purification of arbitrary mixed Gaussian states.
  arXiv  Detail & Related papers  (2020-09-24T18:00:36Z)

This list is automatically generated from the titles and abstracts of the papers in this site.

This site does not guarantee the quality of this site (including all information) and is not responsible for any consequences.