Density functional study of the adsorption of K on the Ag(111) surface

Full-potential gradient corrected density functional calculations of the adsorption of potassium on the Ag(111) surface have been performed. The considered structures are Ag(111) (root 3 x root 3) R30degree-K and Ag(111) (2 x 2)-K. For the lower coverage, fcc, hcp and bridge site; and for the higher coverage all considered sites are practically degenerate. Substrate rumpling is most important for the top adsorption site. The bond length is found to be nearly identical for the two coverages, in agreement with recent experiments. Results from Mulliken populations, bond lengths, core level shifts and work functions consistently indicate a small charge transfer from the potassium atom to the substrate, which is slightly larger for the lower coverage.Comment: to appear in Phys Rev

Ground state properties of heavy alkali halides

We extend previous work on alkali halides by calculations for the heavy-atom species RbF, RbCl, LiBr, NaBr, KBr, RbBr, LiI, NaI, KI, and RbI. Relativistic effects are included by means of energy-consistent pseudopotentials, correlations are treated at the coupled-cluster level. A striking deficiency of the Hartree-Fock approach are lattice constants deviating by up to 7.5 % from experimental values which is reduced to a maximum error of 2.4 % by taking into account electron correlation. Besides, we provide ab-initio data for in-crystal polarizabilities and van der Waals coefficients.Comment: accepted by Phys. Rev.

Moments of spectral functions: Monte Carlo evaluation and verification

The subject of the present study is the Monte Carlo path-integral evaluation of the moments of spectral functions. Such moments can be computed by formal differentiation of certain estimating functionals that are infinitely-differentiable against time whenever the potential function is arbitrarily smooth. Here, I demonstrate that the numerical differentiation of the estimating functionals can be more successfully implemented by means of pseudospectral methods (e.g., exact differentiation of a Chebyshev polynomial interpolant), which utilize information from the entire interval $(-\beta \hbar / 2, \beta \hbar/2)$. The algorithmic detail that leads to robust numerical approximations is the fact that the path integral action and not the actual estimating functional are interpolated. Although the resulting approximation to the estimating functional is non-linear, the derivatives can be computed from it in a fast and stable way by contour integration in the complex plane, with the help of the Cauchy integral formula (e.g., by Lyness' method). An interesting aspect of the present development is that Hamburger's conditions for a finite sequence of numbers to be a moment sequence provide the necessary and sufficient criteria for the computed data to be compatible with the existence of an inversion algorithm. Finally, the issue of appearance of the sign problem in the computation of moments, albeit in a milder form than for other quantities, is addressed.Comment: 13 pages, 2 figure

An Infinite Swapping Approach to the Rare-Event Sampling Problem

We describe a new approach to the rare-event Monte Carlo sampling problem. This technique utilizes a symmetrization strategy to create probability distributions that are more highly connected and thus more easily sampled than their original, potentially sparse counterparts. After discussing the formal outline of the approach and devising techniques for its practical implementation, we illustrate the utility of the technique with a series of numerical applications to Lennard-Jones clusters of varying complexity and rare-event character.Comment: 24 pages, 16 figure

Cohesive properties of alkali halides

We calculate cohesive properties of LiF, NaF, KF, LiCl, NaCl, and KCl with ab-initio quantum chemical methods. The coupled-cluster approach is used to correct the Hartree-Fock crystal results for correlations and to systematically improve cohesive energies, lattice constants and bulk moduli. After inclusion of correlations, we recover 95-98 % of the total cohesive energies. The lattice constants deviate from experiment by at most 1.1 %, bulk moduli by at most 8 %. We also find good agreement for spectroscopic properties of the corresponding diatomic molecules.Comment: LaTeX, 10 pages, 1 figure, accepted by Phys. Rev.

Precise Measurement of Magnetic Field Gradients from Free Spin Precession Signals of 3^{3}He and 129^{129}Xe Magnetometers

We report on precise measurements of magnetic field gradients extracted from transverse relaxation rates of precessing spin samples. The experimental approach is based on the free precession of gaseous, nuclear spin polarized $^3$He and $^{129}$Xe atoms in a spherical cell inside a magnetic guiding field of about 400 nT using LT$_C$ SQUIDs as low-noise magnetic flux detectors. The transverse relaxation rates of both spin species are simultaneously monitored as magnetic field gradients are varied. For transverse relaxation times reaching 100 h, the residual longitudinal field gradient across the spin sample could be deduced to be$|\vec{\nabla}B_z|=(5.6 \pm 0.4)$ pT/cm. The method takes advantage of the high signal-to-noise ratio with which the decaying spin precession signal can be monitored that finally leads to the exceptional accuracy to determine magnetic field gradients at the sub pT/cm scale

Upon the existence of short-time approximations of any polynomial order for the computation of density matrices by path integral methods

In this article, I provide significant mathematical evidence in support of the existence of short-time approximations of any polynomial order for the computation of density matrices of physical systems described by arbitrarily smooth and bounded from below potentials. While for Theorem 2, which is ``experimental'', I only provide a ``physicist's'' proof, I believe the present development is mathematically sound. As a verification, I explicitly construct two short-time approximations to the density matrix having convergence orders 3 and 4, respectively. Furthermore, in the Appendix, I derive the convergence constant for the trapezoidal Trotter path integral technique. The convergence orders and constants are then verified by numerical simulations. While the two short-time approximations constructed are of sure interest to physicists and chemists involved in Monte Carlo path integral simulations, the present article is also aimed at the mathematical community, who might find the results interesting and worth exploring. I conclude the paper by discussing the implications of the present findings with respect to the solvability of the dynamical sign problem appearing in real-time Feynman path integral simulations.Comment: 19 pages, 4 figures; the discrete short-time approximations are now treated as independent from their continuous version; new examples of discrete short-time approximations of order three and four are given; a new appendix containing a short review on Brownian motion has been added; also, some additional explanations are provided here and there; this is the last version; to appear in Phys. Rev.