First-principles energetics of water: a many-body analysis

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for solid, liquid and cluster forms of water. We use a many-body separation of the total energy into its 1-body, 2-body (2B) and beyond-2-body (B2B) components to analyze the deficiencies of two popular DFT approximations. We show how machine-learning methods make this analysis possible for ice structures as well as for water clusters. We find that the crucial energy balance between compact and extended geometries can be distorted by 2B and B2B errors, and that both types of first-principles error are important.

First-principles energetics of water clusters and ice: a many-body analysis.

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for cluster, solid, and liquid forms of water. Recent work has stressed the importance of DFT errors in describing dispersion, but we note that errors in other parts of the energy may also contribute. We obtain information about the nature of DFT errors by using a many-body separation of the total energy into its 1-body, 2-body, and beyond-2-body components to analyze the deficiencies of the popular PBE and BLYP approximations for the energetics of water clusters and ice structures. The errors of these approximations are computed by using accurate benchmark energies from the coupled-cluster technique of molecular quantum chemistry and from quantum Monte Carlo calculations. The systems studied are isomers of the water hexamer cluster, the crystal structures Ih, II, XV, and VIII of ice, and two clusters extracted from ice VIII. For the binding energies of these systems, we use the machine-learning technique of Gaussian Approximation Potentials to correct successively for 1-body and 2-body errors of the DFT approximations. We find that even after correction for these errors, substantial beyond-2-body errors remain. The characteristics of the 2-body and beyond-2-body errors of PBE are completely different from those of BLYP, but the errors of both approximations disfavor the close approach of non-hydrogen-bonded monomers. We note the possible relevance of our findings to the understanding of liquid water.

Reactive Many-Body Expansion for a Protonated Water Cluster.

We generalize the standard many-body expansion technique that is used to approximate the total energy of a molecular system to enable the treatment of chemical reactions by quantum chemical techniques. By considering all possible assignments of atoms to monomer units of the many-body expansion and associating suitable weights with each, we construct a potential energy surface that is a smooth function of the nuclear positions. We derive expressions for this reactive many-body expansion energy and describe an algorithm for its evaluation, which scales polynomially with system size, and therefore will make the method feasible for future condensed phase simulations. We demonstrate the accuracy and smoothness of the resulting potential energy surface on a molecular dynamics trajectory of the protonated water hexamer, using the Hartree-Fock method for the many-body term and Møller-Plesset theory for the low order terms of the many-body expansion.

Microscopic three-body forces and kaon condensation in cold neutrino-trapped matter

We investigate the composition and the equation of state of the kaon condensed phase in neutrino-free and neutrino-trapped star matter within the framework of the Brueckner-Hartree-Fock approach with three-body forces. We find that neutrino trapping shifts the onset density of kaon condensation to a larger baryon density, and reduces considerably the kaon abundance. As a consequence, when kaons are allowed, the equation of state of neutrino-trapped star matter becomes stiffer than the one of neutrino free matter. The effects of different three-body forces are compared and discussed. Neutrino trapping turns out to weaken the role played by the symmetry energy in determining the composition of stellar matter, and thus reduces the difference between the results obtained by using different three-body forces.

New parameters of many-body potentials: studying the thermal and mechanical properties of noble metals

New parameters of nearest-neighbor EAM (1N-EAM), n-th neighbor EAM (NN-EAM), and the second-moment approximation to the tight-binding (TB-SMA) potentials are obtained by fitting experimental data at different temperatures. In comparison with the available many-body potentials, our results suggest that the 1N-EAM potential with the new parameters is the best description of atomic interactions in studying the thermal expansion of noble metals. For mechanical properties, it is suggested that the elastic constants should be calculated in the experimental zero-stress states for all three potentials. Furthermore, for NNEAM and TB-SMA potentials, the calculated results approach the experimental data as the range of the atomic interaction increases from the first-neighbor to the sixth-neighbor distance.

Waveguide QED: Many-body bound-state effects in coherent and Fock-state scattering from a two-level system

Strong coupling between a two-level system (TLS) and bosonic modes produces dramatic quantum optics effects. We consider a one-dimensional continuum of bosons coupled to a single localized TLS, a system which may be realized in a variety of plasmonic, photonic, or electronic contexts. We present the exact many-body scattering eigenstate obtained by imposing open boundary conditions. Multiphoton bound states appear in the scattering of two or more photons due to the coupling between the photons and the TLS. Such bound states are shown to have a large effect on scattering of both Fock- and coherent-state wave packets, especially in the intermediate coupling-strength regime. We compare the statistics of the transmitted light with a coherent state having the same mean photon number: as the interaction strength increases, the one-photon probability is suppressed rapidly, and the two- and three-photon probabilities are greatly enhanced due to the many-body bound states. This results in non-Poissonian light. © 2010 The American Physical Society.

Accurate localized resolution of identity approach for linear-scaling hybrid density functionals and for many-body perturbation theory

© 2015 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.A key component in calculations of exchange and correlation energies is the Coulomb operator, which requires the evaluation of two-electron integrals. For localized basis sets, these four-center integrals are most efficiently evaluated with the resolution of identity (RI) technique, which expands basis-function products in an auxiliary basis. In this work we show the practical applicability of a localized RI-variant ('RI-LVL'), which expands products of basis functions only in the subset of those auxiliary basis functions which are located at the same atoms as the basis functions. We demonstrate the accuracy of RI-LVL for Hartree-Fock calculations, for the PBE0 hybrid density functional, as well as for RPA and MP2 perturbation theory. Molecular test sets used include the S22 set of weakly interacting molecules, the G3 test set, as well as the G2-1 and BH76 test sets, and heavy elements including titanium dioxide, copper and gold clusters. Our RI-LVL implementation paves the way for linear-scaling RI-based hybrid functional calculations for large systems and for all-electron many-body perturbation theory with significantly reduced computational and memory cost.

Many-body entanglement in decoherence processes

A pure state decoheres into a mixed state as it entangles with an environment. When an entangled two-mode system is embedded in a thermal environment, however, each mode may not be entangled with its environment by their simple linear interaction. We consider an exactly solvable model to study the dynamics of a total system, which is composed of an entangled two-mode system and a thermal environment. The Markovian interaction with the environment is concerned with an array of infinite number of beam splitters. It is shown that many-body entanglement of the system and the environment may play a crucial role in the process of disentangling the system.

Independent particle descriptions of tunneling using the many-body quantum transport approach

Currents across thin insulators are commonly taken as single electrons moving across classically forbidden regions; this independent particle picture is well-known to describe most tunneling phenomena. Examining quantum transport from a different perspective, i.e., by explicit treatment of electron-electron interactions, we evaluate different single particle approximations with specific application to tunneling in metal-molecule-metal junctions. We find maximizing the overlap of a Slater determinant composed of single-particle states to the many-body current-carrying state is more important than energy minimization for defining single-particle approximations in a system with open boundary conditions. Thus the most suitable single particle effective potential is not one commonly in use by electronic structure methods, such as the Hartree-Fock or Kohn-Sham approximations.

Control-limited perfect state transfer, quantum stochastic resonance, and many-body entangling gate in imperfect qubit registers

We propose a protocol for perfect quantum state transfer that is resilient to a broad class of realistic experimental imperfections, including noise sources that could be modeled either as independent Markovian baths or as certain forms of spatially correlated environments. We highlight interesting connections between the fidelity of state transfer and quantum stochastic resonance effects. The scheme is flexible enough to act as an effective entangling gate for the generation of genuine multipartite entanglement in a control-limited setting. Possible experimental implementations using superconducting qubits are also briefly discussed.

Many-body theory of gamma spectra from positron-atom annihilation

A many-body theory approach to the calculation of gamma spectra of positron annihilation on many-electron atoms is developed. We evaluate the first-order correlation correction to the annihilation vertex and perform numerical calculations for the noble gas atoms. Extrapolation with respect to the maximal orbital momentum of the intermediate electron and positron states is used to achieve convergence. The inclusion of correlation corrections improves agreement with experimental gamma spectra.