Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFirst-principles many-body non-additive polarization energies from monomer and dimer calculations only : A case study on water
138
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The many-body polarization energy is the major source of non-additivity in strongly polar systems such as water. This non-additivity is often considerable and must be included, if only in an average manner, to correctly describe the physical properties of the system. Models for the polarization energy are usually parameterized using experimental data, or theoretical estimates of the many-body effects. Here we show how many-body polarization models can be developed for water complexes using data for the monomer and dimer only using ideas recently developed in the field of intermolecular perturbation theory and state-of-the-art approaches for calculating distributed molecular properties based on the iterated stockholder atoms (ISA) algorithm. We show how these models can be calculated, and validate their accuracy in describing the many-body non-additive energies of a range of water clusters. We further investigate their sensitivity to the details of the polarization damping models used. We show how our very best polarization models yield many-body energies that agree with those computed with coupled-cluster methods, but at a fraction of the computational cost.
rate research
Read More
During compression of a water dimer calculated with high-precision first-principles methods, the trends of H-bond and O-H bond lengths show quantum effect of the electronic structure. We found that the H-bond length keeps decreasing, while the O-H bond length increases up to the stable point and decreases beyond it when the water dimer is further compressed. The remarkable properties are different from those observed in most previous researches which can be understood and extrapolated through classical simulation. The observations can be explained by the decrease in orbital overlap and change in the exchange repulsion interaction between water monomers. The dominant interaction between water monomers changes from electrostatic interaction to exchange repulsion at the turning point of the O-H bond length when the O...O distance is decreased. These findings highlight the quantum effect on the hydrogen bond in extreme conditions and play an important role in the recognition of the hydrogen bond structure and mechanism.
The treatment of atomic anions with Kohn-Sham density functional theory (DFT) has long been controversial since the highest occupied molecular orbital (HOMO) energy, $E_{HOMO}$, is often calculated to be positive with most approximate density functionals. We assess the accuracy of orbital energies and electron affinities for all three rows of elements in the periodic table (H-Ar) using a variety of theoretical approaches and customized basis sets. Among all of the theoretical methods studied here, we find that a non-empirically tuned range-separated approach (constructed to satisfy DFT-Koopmans theorem for the anionic electron system) provides the best accuracy for a variety of basis sets - even for small basis sets where most functionals typically fail. Previous approaches to solve this conundrum of positive $E_{HOMO}$ values have utilized non-self-consistent methods; however electronic properties, such as electronic couplings/gradients (which require a self-consistent potential and energy), become ill-defined with these approaches. In contrast, the non-empirically tuned range-separated procedure used here yields well-defined electronic couplings/gradients and correct $E_{HOMO}$ values since both the potential and resulting electronic energy are computed self-consistently. Orbital energies and electron affinities are further analyzed in the context of the electronic energy as a function of electronic number (including fractional numbers of electrons) to provide a stringent assessment of self-interaction errors for these complex anion systems.
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.
The homogeneous electron gas (HEG) is a key ingredient in the construction of most exchange-correlation functionals of density-functional theory. Often, the energy of the HEG is parameterized as a function of its spin density $n$, leading to the local density approximation (LDA) for inhomogeneous systems. However, the connection between the electron density and kinetic energy density of the HEG can be used to generalize the LDA by evaluating it on a weighted geometric average of the local spin density and the spin density of a HEG that has the local kinetic energy density of the inhomogeneous system, with a mixing ratio $x$. This leads to a new family of functionals that we term meta-local density approximations (meta-LDAs), which are still exact for the HEG, which are derived only from properties of the HEG, and which form a new rung of Jacobs ladder of density functionals. The first functional of this ladder, the local $tau$ approximation (LTA) of Ernzerhof and Scuseria that corresponds to $x=1$ is unfortunately not stable enough to be used in self-consistent field calculations, because it leads to divergent potentials as we show in this work. However, a geometric averaging of the LDA and LTA densities with smaller values of $x$ not only leads to numerical stability of the resulting functional, but also yields more accurate exchange energies in atomic calculations than the LDA, the LTA, or the tLDA functional ($x=1/4$) of Eich and Hellgren. We choose $x=0.50$ as it gives the best total energy in self-consistent exchange-only calculations for the argon atom. Atomization energy benchmarks confirm that the choice $x=0.50$ also yields improved energetics in combination with correlation functionals in molecules, almost eliminating the well-known overbinding of the LDA and reducing its error by two thirds.
Reliable first-principles calculations of electrochemical processes require accurate prediction of the interfacial capacitance, a challenge for current computationally-efficient continuum solvation methodologies. We develop a model for the double layer of a metallic electrode that reproduces the features of the experimental capacitance of Ag(100) in a non-adsorbing, aqueous electrolyte, including a broad hump in the capacitance near the Potential of Zero Charge (PZC), and a dip in the capacitance under conditions of low ionic strength. Using this model, we identify the necessary characteristics of a solvation model suitable for first-principles electrochemistry of metal surfaces in non-adsorbing, aqueous electrolytes: dielectric and ionic nonlinearity, and a dielectric-only region at the interface. The dielectric nonlinearity, caused by the saturation of dipole rotational response in water, creates the capacitance hump, while ionic nonlinearity, caused by the compactness of the diffuse layer, generates the capacitance dip seen at low ionic strength. We show that none of the previously developed solvation models simultaneously meet all these criteria. We design the Nonlinear Electrochemical Soft-Sphere solvation model (NESS) which both captures the capacitance features observed experimentally, and serves as a general-purpose continuum solvation model.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
