Source author record

Markus Holzmann

Markus Holzmann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

35works
18topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2021arXiv

A note on the three dimensional Dirac operator with zigzag type boundary conditions

In this note the three dimensional Dirac operator $A_m$ with boundary conditions, which are the analogue of the two dimensional zigzag boundary conditions, is investigated. It is shown that $A_m$ is self-adjoint in $L^2(Ω;\mathbb{C}^4)$ for any open set $Ω\subset \mathbb{R}^3$ and its spectrum is described explicitly in terms of the spectrum of the Dirichlet Laplacian in $Ω$. In particular, whenever the spectrum of the Dirichlet Laplacian is purely discrete, then also the spectrum of $A_m$ consists of discrete eigenvalues that accumulate at $\pm \infty$ and one additional eigenvalue of infinite multiplicity.

preprint2021arXiv

Efficient computations of continuous action densities of states for lattice models

The Logarithmic Linear Relaxation (LLR) algorithm is an efficient method for computing densities of states for systems with a continuous spectrum. A key feature of this method is exponential error reduction, which allows us to evaluate the density of states of a system over hundreds of thousands of orders of magnitude with a fixed level of relative accuracy. As a consequence of exponential error reduction, the LLR method provides a robust alternative to traditional Monte Carlo calculations in cases in which states suppressed by the Boltzmann weight play nevertheless a relevant role, e.g., as transition regions between dominant configuration sets. After reviewing the algorithm, we will show an application in U(1) Lattice Gauge Theory that has enabled us to obtain the most accurate estimate of the critical coupling with modest computational resources, defeating exponential tunneling times between metastable vacua. As a further showcase, we will then present an application of the LLR method to the decorrelation of the topological charge in SU(3) Lattice Gauge Theory near the continuum limit. Finally, we will review in general applications of the LLR algorithm to systems affected by a strong sign problem and discuss the case of the Bose gas at finite chemical potential.

preprint2020arXiv

Electronic structure and optical properties of quantum crystals from first principles calculations in the Born-Oppenheimer approximation

We develop a formalism to accurately account for the renormalization of electronic structure due to quantum and thermal nuclear motions within the Born-Oppenheimer approximation. We focus on the fundamental energy gap obtained from electronic addition and removal energies from Quantum Monte Carlo calculations in either the canonical or grand canonical ensembles. The formalism applies as well to effective single electron theories such as those based on Density Functional Theory. We show that electronic (Bloch) crystal momentum can be restored by marginalizing the total electron-ion wave function with respect to the nuclear equilibrium distribution, and we describe an explicit procedure to establish the band structure of electronic excitations for quantum crystals within the Born-Oppenheimer approximation. Based on the Kubo-Greenwood equation, we discuss the effects of nuclear motion on optical conductivity. Our methodology applies to the low temperature regime where nuclear motion is quantized and in general differs from the semi-classical approximation. We apply our method to study the electronic structure of C2/c-24 crystalline hydrogen at 200K and 250 GPa and discuss the optical absorption profile of hydrogen crystal at 200K and carbon diamond at 297K.

preprint2020arXiv

Energy gap closure of crystalline molecular hydrogen with pressure

We study the gap closure with pressure of crystalline molecular hydrogen. The gaps are obtained from grand-canonical Quantum Monte Carlo methods properly extended to quantum and thermal crystals, simulated by Coupled Electron Ion Monte Carlo. Nuclear zero point effects cause a large reduction in the gap ($\sim 2eV$). \CP{Depending on the structure,} the fundamental indirect gap closes \CP{between 380GPa and} 530GPa for ideal crystals and 330-380GPa for quantum crystals. Beyond this pressure the system enters into a bad metal phase where the density of states at the Fermi level increases with pressure up to $\sim$450\CP{-500} GPa when the direct gap closes. Our work partially supports the interpretation of recent experiments in high pressure hydrogen.

preprint2020arXiv

Quasi-condensation of bilayer excitons in a periodic potential

We study two-dimensional excitons confined in a lattice potential, for high fillings of the lattice sites. We show that a quasi-condensate is possibly formed for small values of the lattice depth, but for larger ones the critical phase-space density for quasi-condensation rapidly exceeds our experimental reach, due to the increase of the excitons effective mass. On the other hand, in the regime of a deep lattice potential where excitons are strongly localised at the lattice sites, we show that an array of phase-independent quasi-condensates, different from a Mott insulating phase, is realised.

preprint2019arXiv

Electronic band gaps from Quantum Monte Carlo methods

We develop a method for calculating the fundamental electronic gap of semiconductors and insulators using grand canonical Quantum Monte Carlo simulations. We discuss the origin of the bias introduced by supercell calculations of finite size and show how to correct the leading and subleading finite size errors either based on observables accessible in the finite-sized simulations or from DFT calculations. Our procedure is applied to solid molecular hydrogen and compared to experiment for carbon and silicon crystals.

preprint2019arXiv

Free energy of the self-interacting relativistic lattice Bose gas at finite density

The density of state approach has recently been proposed as a potential route to circumvent the sign problem in systems at finite density. In this study, using the Linear Logarithmic Relaxation (LLR) algorithm, we extract the generalised density of states, which is defined in terms of the imaginary part of the action, for the self-interacting relativistic lattice Bose gas at finite density. After discussing the implementation and testing the reliability of our approach, we focus on the determination of the free energy difference between the full system and its phase-quenched counterpart. Using a set of lattices ranging from $4^4$ to $16^4$ , we show that in the low density phase, this overlap free energy can be reliably extrapolated to the thermodynamic limit. The numerical precision we obtain with the LLR method allows us to determine with sufficient accuracy the expectation value of the phase factor, which is used in the calculation of the overlap free energy, down to values of ${\cal O}(10^{-480})$. When phase factor measurements are extended to the dense phase, a change of behaviour of the overlap free energy is clearly visible as the chemical potential crosses a critical value. Using fits inspired by the approximate validity of mean-field theory, which is confirmed by our simulations, we extract the critical chemical potential as the non-analyticity point in the overlap free energy, obtaining a value that is in agreement with other determinations. Implications of our findings and potential improvements of our methodology are also discussed.

preprint2019arXiv

Quantum Monte Carlo Compton profiles of solid and liquid lithium

We computed the Compton profile of solid and liquid lithium using quantum Monte Carlo (QMC) and compared with recent experimental measurements obtaining good agreement. Importantly, we find it crucial to account for proper core-valence orthogonalization and to address density differences when comparing with experiment. To account for disorder effects, we sampled finite-temperature configurations using molecular dynamics (MD), then performed diffusion Monte Carlo (DMC) simulations on each configuration. We used Slater-Jastrow wavefunctions and grand-canonical twist-averaged boundary conditions. A QMC pseudopotential correction, derived from an all-electron DMC simulation of the perfect crystal was also used. Our calculations provide the first all-electron QMC benchmark for the Compton profile of lithium crystal and pseudopotential-corrected QMC Compton profiles for both the liquid and solid.

preprint2019arXiv

Self-adjoint Dirac operators on domains in $\mathbb{R}^3$

In this paper the spectral and scattering properties of a family of self-adjoint Dirac operators in $L^2(Ω; \mathbb{C}^4)$, where $Ω\subset \mathbb{R}^3$ is either a bounded or an unbounded domain with a compact $C^2$-smooth boundary, are studied in a systematic way. These operators can be viewed as the natural relativistic counterpart of Laplacians with Robin boundary conditions. Among the Dirac operators treated here is also the so-called MIT bag operator, which has been used by physicists and more recently was discussed in the mathematical literature. Our approach is based on abstract boundary triple techniques from extension theory of symmetric operators and a thorough study of certain classes of (boundary) integral operators, that appear in a Krein-type resolvent formula. The analysis of the perturbation term in this formula leads to a description of the spectrum and a Birman-Schwinger principle, a qualitative understanding of the scattering properties in the case that $Ω$ is unbounded, and corresponding trace formulas.

preprint2019arXiv

Two-dimensional Dirac operators with singular interactions supported on closed curves

We study the two-dimensional Dirac operator with an arbitrary combination of electrostatic and Lorentz scalar $δ$-interactions of constant strengths supported on a smooth closed curve. For any combination of the coupling constants a rigorous description of the self-adjoint realizations of the operators is given and the qualitative spectral properties are described. The analysis covers also all so-called critical combinations of coupling constants, for which there is a loss of regularity in the operator domain. In this case, if the mass is non-zero, the resulting operator has an additional point in the essential spectrum, and the position of this point inside the central gap can be made arbitrary by a suitable choice of the coupling constants. The analysis is based on a combination of the extension theory of symmetric operators with a detailed study of boundary integral operators viewed as periodic pseudodifferential operators.

preprint2017arXiv

Electron localization properties in high pressure hydrogen at the liquid-liquid phase transition by Coupled Electron-Ion Monte Carlo

We analyze in detail the electronic properties of high pressure hydrogen around the liquid-liquid phase transition based on Coupled Electron-Ion Monte Carlo calculations. Computing the off-diagonal single particle density matrix and the momentum distribution we discuss localization properties of the electrons. The abrupt changes of these distributions indicate a metal to insulator transition occurring together with the structural transition from the atomic to molecular fluid. We further discuss the electron-proton and electron-electron pair correlation functions, which also change abruptly at the transition.

preprint2016arXiv

On the Fermion Sign Problem in Imaginary-Time Projection Continuum Quantum Monte Carlo with Local Interaction

We use the Shadow Wave Function formalism as a convenient model to study the fermion sign problem affecting all projector Quantum Monte Carlo methods in continuum space. We demonstrate that the efficiency of imaginary time projection algorithms decays exponentially with increasing number of particles and/or imaginary-time propagation. Moreover, we derive an analytical expression that connects the localization of the system with the magnitude of the sign problem, illustrating this prediction through some numerical results. Finally, we discuss the fermion sign problem computational complexity and methods for alleviating its severity.

preprint2016arXiv

Quasi-long-range order in trapped two-dimensional Bose gases

We study the fate of algebraic decay of correlations in a harmonically trapped two-dimensional degenerate Bose gas. The analysis is inspired by recent experiments on ultracold atoms where power-law correlations have been observed despite the presence of the external potential. We generalize the spin wave description of phase fluctuations to the trapped case and obtain an analytical expression for the one-body density matrix within this approximation. We show that algebraic decay of the central correlation function persists to lengths of about 20% of the Thomas--Fermi radius. We establish that the trap-averaged correlation function decays algebraically with a strictly larger exponent weakly changing with trap size and find indications that the recently observed enhanced scaling exponents receive significant contributions from the normal component of the gas. We discuss radial and angular correlations and propose a local correlation approximation which captures the correlations very well. Our analysis goes beyond the usual local density approximation and the developed summation techniques constitute a powerful tool to investigate correlations in inhomogeneous systems.

preprint2016arXiv

Theory of Finite Size Effects for Electronic Quantum Monte Carlo Calculations of Liquids and Solids

Concentrating on zero temperature Quantum Monte Carlo calculations of electronic systems, we give a general description of the theory of finite size extrapolations of energies to the thermodynamic limit based on one and two-body correlation functions. We introduce new effective procedures, such as using the potential and wavefunction split-up into long and short range functions to simplify the method and we discuss how to treat backflow wavefunctions. Then we explicitly test the accuracy of our method to correct finite size errors on example hydrogen and helium many-body systems and show that the finite size bias can be drastically reduced for even small systems.

preprint2015arXiv

Benchmarking Hydrogen-Helium Mixtures with QMC: Energetics, Pressures, and Forces

An accurate understanding of the phase diagram of dense hydrogen and helium mixtures is a crucial component in the construction of accurate models of Jupiter, Saturn, and Jovian extrasolar planets. Though DFT based first principles methods have the potential to provide the accuracy and computational efficiency required for this task, recent benchmarking in hydrogen has shown that achieving this accuracy requires a judicious choice of functional, and a quantification of the errors introduced. In this work, we present a quantum Monte Carlo based benchmarking study of a wide range of density functionals for use in hydrogen-helium mixtures at thermodynamic conditions relevant for Jovian planets. Not only do we continue our program of benchmarking energetics and pressures, but we deploy QMC based force estimators and use them to gain insights into how well the local liquid structure is captured by different density functionals. We find that TPSS, BLYP and vdW-DF are the most accurate functionals by most metrics, and that the enthalpy, energy, and pressure errors are very well behaved as a function of helium concentration. Beyond this, we highlight and analyze the major error trends and relative differences exhibited by the major classes of functionals, and estimate the magnitudes of these effects when possible.

preprint2015arXiv

Iterative backflow renormalization procedure for many-body ground state wave functions of strongly interacting normal Fermi liquids

We show how a ground state trial wavefunction of a Fermi liquid can be systematically improved introducing a sequence of renormalized coordinates through an iterative backflow transformation. We apply this scheme to calculate the ground state energy of liquid $^3$He in two dimensions at freezing density using variational and fixed-node diffusion Monte Carlo. Comparing with exact transient estimate results for systems with small number of particles, we find that variance extrapolations provide accurate results for the true ground state together with stringent lower bounds. For larger systems these bounds can in turn be used to quantify the systematic bias of fixed-node calculations. These wave functions are size consistent and the scaling of their computational complexity with the number of particles is the same as for standard backflow wave functions.

preprint2014arXiv

Molecular-Atomic Transition in the Deuterium Hugoniot with Coupled Electron Ion Monte Carlo

We have performed accurate simulations of the Deuterium Hugoniot using Coupled Electron Ion Monte Carlo (CEIMC). Using highly accurate quantum Monte Carlo methods for the electrons, we study the region of maximum compression along the principal Hugoniot, where the system undergoes a continuous transition from a molecular fluid to a monatomic fluid. We include all relevant physical corrections so that a direct comparison to experiment can be made. Around 50 GPa we found a maximum compression of 4.85, roughly 10% larger than previous theoretical predictions and experimental data but still compatible with the latter because of their large uncertainty.

preprint2013arXiv

Exchange-correlation energy for the 3D homogeneous electron gas at arbitrary temperature

We fit finite-temperature path integral Monte Carlo calculations of the exchange-correlation energy of the 3D finite-temperature homogeneous electron gas in the warm-dense regime (r_{s} = (3/4πn)^{1/3} a_{B}^{-1} < 40 and Θ= T/T_{F} > 0.0625). In doing so, we construct a Padé approximant which collapses to Debye-Hückel theory in the high-temperature, low-density limit. Likewise, the zero-temperature limit matches the numerical results of ground-state quantum Monte Carlo, as well as analytical results in the high-density limit.

preprint2013arXiv

Hartree-Fock Ground State Phase Diagram of Jellium

We calculate the ground state phase diagram of the homogeneous electron gas in three dimensions within the Hartree-Fock approximation and show that broken symmetry states are energetically favored at any density against the homogeneous Fermi gas state with isotropic Fermi surface. At high density, we find metallic spin-unpolarized solutions where electronic charge and spin density form an incommensurate crystal having more crystal sites than electrons. For $r_s \to 0$, our solutions approach pure spin-density waves, whereas the commensurate Wigner crystal is favored at lower densities, $r_s > 3.4$. Decreasing the density, the system undergoes several structural phase transitions with different lattice symmetries. The polarization transition occurs around $r_s \approx 8.5$.

preprint2013arXiv

Universal Superfluid Transition and Transport Properties of Two-Dimensional Dirty Bosons

We study the phase diagram of two-dimensional, interacting bosons in the presence of a correlated disorder in continuous space, using large-scale finite temperature quantum Monte Carlo simulations. We show that the superfluid transition is strongly protected against disorder. It remains of the Berezinskii-Kosterlitz-Thouless type up to disorder strengths comparable to the chemical potential. Moreover, we study the transport properties in the strong disorder regime where a zero-temperature Bose-glass phase is expected. We show that the conductance exhibits a thermally activated behavior vanishing only at zero temperature. Our results point towards the existence of Bose bad-metal phase as a precursor of the Bose-glass phase.

preprint2012arXiv

Effect of disorder close to the superfluid transition in a two-dimensional Bose gas

We experimentally study the effect of disorder on trapped quasi two-dimensional (2D) 87Rb clouds in the vicinity of the Berezinskii-Kosterlitz-Thouless (BKT) phase transition. The disorder correlation length is of the order of the Bose gas characteristic length scales (thermal de Broglie wavelength, healing length) and disorder thus modifies the physics at a microscopic level. We analyze the coherence properties of the cloud through measurements of the momentum distributions, for two disorder strengths, as a function of its degeneracy. For moderate disorder, the emergence of coherence remains steep but is shifted to a lower entropy. In contrast, for strong disorder, the growth of coherence is hindered. Our study is an experimental realization of the dirty boson problem in a well controlled atomic system suitable for quantitative analysis.

preprint2012arXiv

Momentum distribution and Compton profile by the ab initio GW approximation

We present two possible approaches to calculate the momentum distribution n(p) and the Compton profile within the framework of the ab initio GW approximation on the self-energy. The approaches are based on integration of the Green's function along either the real or the imaginary axes. Examples will be presented on the jellium model and on real bulk sodium. Advantages and drawbacks of both methods are discussed in comparison with accurate quantum Monte Carlo (QMC) calculations and x-ray Compton scattering experiments. We illustrate the effect of many-body correlations and disentangle them from band-structure and anisotropy effects by a comparison with density functional theory in the local density approximation. Our results suggest the use of G0W0 momentum distributions as reference for future experiments and theory developments.

preprint2011arXiv

Condensed Fraction of an Atomic Bose Gas Induced by Critical Correlations

We study the condensed fraction of a harmonically-trapped atomic Bose gas at the critical point predicted by mean-field (MF) theory. The non-zero condensed fraction $f_0$ is induced by critical correlations which increase the transition temperature $T_c$ above $\T_c^{MF}$. Unlike the $T_c$ shift in a trapped gas, $f_0$ is sensitive only to the critical behaviour in the quasi-uniform part of the cloud near the trap centre. To leading order in the interaction parameter $a/λ_0$, where $a$ is the s-wave scattering length and $λ_0$ the thermal wavelength, we expect a universal scaling $f_0 \propto (a/λ_0)^4$. We experimentally verify this scaling using a Feshbach resonance to tune $a/λ_0$. Further, using the local density approximation, we compare our measurements with the universal result obtained from Monte-Carlo simulations for a uniform system, and find excellent quantitative agreement.

preprint2011arXiv

Finite-size analysis of the Fermi liquid properties of the homogeneous electron gas

We analyze the extrapolation to the thermodynamic limit of Fermi liquid properties of the homogeneous electron gas in two and three dimensions. Using field theory, we explicitly calculate finite-size effects of the total energy, the renormalization factor, and the effective mass at the Fermi surface within the random phase approximation (RPA) and discuss the validity for general metallic systems.

preprint2011arXiv

The momentum distribution of the homogeneous electron gas

We calculate the off-diagonal density matrix of the homogeneous electron gas at zero temperature using unbiased Reptation Monte Carlo for various densities and extrapolate the momentum distribution, and the kinetic and potential energies to the thermodynamic limit. Our results on the renormalization factor allows us to validate approximate G_0W_0 calculations concerning quasiparticle properties over a broad density region (1 <= r_s <= 10) and show that near the Fermi surface, vertex corrections and self-consistency aspects almost cancel each other out.

preprint2010arXiv

High-energy excitonic effects in graphite and graphene

We present ab initio many-body calculations of the optical absorption in bulk graphite, graphene and bilayer of graphene. Electron-hole interaction is included solving the Bethe-Salpeter equation on top of a GW quasiparticle electronic structure. For all three systems, we observe strong excitonic effects at high energy, well beyond the continuum of π-> π* transitions. In graphite, these affect the onset of σ-> σ* transitions. In graphene, we predict an excitonic resonance at 8.3 eV arising from a background continuum of dipole forbidden transitions. In the graphene bilayer, the resonance is shifted to 9.6 eV. Our results for graphite are in good agreement with experiments.

preprint2010arXiv

Momentum Distribution and Renormalization Factor in Sodium and the Electron Gas

We present experimental and theoretical results on the momentum distribution and the quasiparticle renormalization factor in sodium. From an x-ray Compton-profile measurement of the valence-electron momentum density, we derive its discontinuity at the Fermi wavevector. This yields an accurate measure of the renormalization factor that we compare with quantum Monte Carlo and GW calculations performed both on crystalline sodium and on the homogeneous electron gas. Our calculated results are in good agreement with the experiment.

preprint2010arXiv

The equilibrium state of a trapped two-dimensional Bose gas

We study experimentally and numerically the equilibrium density profiles of a trapped two-dimensional $^{87}$Rb Bose gas, and investigate the equation of state of the homogeneous system using the local density approximation. We find a clear discrepancy between in-situ measurements and Quantum Monte Carlo simulations, which we attribute to a non-linear variation of the optical density of the atomic cloud with its spatial density. However, good agreement between experiment and theory is recovered for the density profiles measured after time-of-flight, taking advantage of their self-similarity in a two-dimensional expansion.

preprint2009arXiv

Universal correlations and coherence in quasi-two-dimensional trapped Bose gases

We study the quasi-two-dimensional Bose gas in harmonic traps at temperatures above the Kosterlitz-Thouless transition, where the gas is in the normal phase. We show that mean-field theory takes into account the dominant interaction effects for experimentally relevant trap geometries. Comparing with Quantum Monte Carlo calculations, we quantify the onset of the fluctuation regime, where correlations beyond mean-field become important. Although the density profile depends on the microscopic parameters of the system, we show that the correlation density (the difference between the exact and the mean-field density) is accurately described by a universal expression, obtained from classical-field calculations of the homogeneous strictly two-dimensional gas. Deviations from universality, due to the finite value of the interaction or to the trap geometry, are shown to be small for current experiments. We further study coherence and pair correlations on a microscopic scale. Finite-size effects in the off-diagonal density matrix allows us to characterize the cross-over from Kosterlitz-Thouless to Bose-Einstein behavior for small particle numbers. Bose-Einstein condensation occurs below a characteristic number of particles which rapidly diverges with vanishing interactions.

preprint2008arXiv

Renormalization factor and effective mass of the two-dimensional electron gas

We calculate the momentum distribution of the Fermi liquid phase of the homogeneous, two-dimensional electron gas. We show that, close to the Fermi surface, the momentum distribution of a finite system with $N$ electrons approaches its thermodynamic limit slowly, with leading order corrections scaling as $N^{-1/4}$. These corrections dominate the extrapolation of the renormalization factor, $Z$, and the single particle effective mass, $m^*$, to the infinite system size. We show how convergence can be improved analytically. In the range $1 \le r_s \le 10$, we get a lower renormalization factor $Z$ and a higher effective mass, $m^*>m$, compared to the perturbative RPA values.

preprint2006arXiv

Many-body wavefunctions for normal liquid $^3$He

We present new trial wave-functions which include 3-body correlations into the backflow coordinates and a 4-body symmetric potential. We show that our wavefunctions lower the energy enough to stabilize the ground state energies of normal liquid $^3$He in the unpolarized state at all pressures in agreement with experiment; however, quantitative discrepancies remain. Further, we include strong spin coupling into the Fermi liquid by adapting pairing wave functions. We demonstrate a new, numerically stable method to evaluate pairing functions which is also useful for Path Integrals calculations at low, but non-zero temperatures.

preprint2005arXiv

The Kosterlitz-Thouless-Berezinskii transition of homogeneous and trapped Bose gases in two dimensions

We derive the scaling structure of the Kosterlitz-Thouless-Berezinskii (KTB) transition temperature of a homogeneous Bose gas in two dimensions within diagrammatic perturbation theory. Approaching the system from above the transition, we calculate the critical temperature, $T_{KT}$, and show how the superfluid mass density emerges from Josephson's relation as an interplay between the condensate density in a finite size system, and the infrared structure of the single particle Green's function. We then discuss the trapped two-dimensional Bose gas, where the interaction changes the transition qualitatively from Bose-Einstein in an ideal gas to a KTB transition in the thermodynamic limit. We show that the transition temperature lies below the ideal Bose-Einstein transition temperature, and calculate the first correction in terms of the interparticle interactions. The jump of the total superfluid mass at the transition is suppressed in a trapped system.

preprint2004arXiv

Optimized Periodic Coulomb Potential in Two Dimension

The 1/r Coulomb potential is calculated for a two dimensional system with periodic boundary conditions. Using polynomial splines in real space and a summation in reciprocal space we obtain numerically optimized potentials which allow us efficient calculations of any periodic (long-ranged) potential up to high precision. We discuss the parameter space of the optimized potential for the periodic Coulomb potential. Compared to the analytic Ewald potential, the optimized potentials can reach higher precisions by up to several orders of magnitude. We explicitly give simple expressions for fast calculations of the periodic Coulomb potential where the summation in reciprocal space is reduced to a few terms.

preprint1995arXiv

Relation between energy shifts and relaxation rates for a small system coupled to a reservoir

For a small system the coupling to a reservoir causes energy shifts as well as transitions between the system's energy levels. We show for a general stationary situation that the energy shifts can essentially be reduced to the relaxation rates. The effects of reservoir fluctuations and self reaction are treated separately. We apply the results to a two-level atom coupled to a reservoir which may be the vacuum of a radiation field.