ترغب بنشر مسار تعليمي؟ اضغط هنا

Meta-local density functionals: a new rung on Jacobs ladder

88   0   0.0 ( 0 )
 نشر من قبل Susi Lehtola
 تاريخ النشر 2020
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

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.



قيم البحث

اقرأ أيضاً

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 properti es 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.
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 bo nd 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.
Kernel ridge regression is used to approximate the kinetic energy of non-interacting fermions in a one-dimensional box as a functional of their density. The properties of different kernels and methods of cross-validation are explored, and highly accu rate energies are achieved. Accurate {em constrained optimal densities} are found via a modified Euler-Lagrange constrained minimization of the total energy. A projected gradient descent algorithm is derived using local principal component analysis. Additionally, a sparse grid representation of the density can be used without degrading the performance of the methods. The implications for machine-learned density functional approximations are discussed.
The stationary nonempirical simulations of Na+(H2O)n clusters with n in a range of 28 to 51 carried out at the density functional level with a hybrid B3LYP functional and the Born-Oppenheimer molecular dynamics modeling of the size selected clusters reveal the interrelated structural and energetic peculiarities of sodium hydration structures. Surface, bulk, and transient configurations of the clusters are distinguished with the different location of the sodium nucleus (close to either the spatial center of the structure or one of its side faces) and its consistently changing coordination number (which typically equals five or six). The <rNaO> mean Na-O distances for the first-shell water molecules are found to depend both on the spatial character of the structure and the local coordination of sodium. The <rNaO> values are compared to different experimental estimates, and the virtual discrepancy of the latter is explained based on the results of the cluster simulations. Different coordination neighborhoods of sodium are predicted depending on its local fraction in the actual specimens.
Chirped-Pulse millimetre-Wave (CPmmW) rotational spectroscopy provides a new class of information about photolysis transition state(s). Measured intensities in rotational spectra determine species-isomer-vibrational populations, provided that rotatio nal populations can be thermalized. The formation and detection of S0 vinylidene is discussed in the limits of low and high initial rotational excitation. CPmmW spectra of 193 nm photolysis of Vinyl Cyanide (Acrylonitrile) contain J=0-1 transitions in more than 20 vibrational levels of HCN, HNC, but no transitions in vinylidene or highly excited local-bender vibrational levels of acetylene. Reasons for the non-observation of the vinylidene co-product of HCN are discussed.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا