Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userUniversal iso-density polarizable continuum model for molecular solvents
156
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Implicit electron-density solvation models based on joint density-functional theory offer a computationally efficient solution to the problem of calculating thermodynamic quantities of solvated systems from firstprinciples quantum mechanics. However, despite much recent interest in such models, to date the applicability of such models in the plane-wave context to non-aqueous solvents has been limited because the determination of the model parameters requires fitting to a large database of experimental solvation energies for each new solvent considered. This work presents an alternate approach which allows development of new iso-density models for a large class of protic and aprotic solvents from only simple, single-molecule ab initio calculations and readily available bulk thermodynamic data.
rate research
Read More
A new empirical potential for efficient, large scale molecular dynamics simulation of water is presented. The HIPPO (Hydrogen-like Intermolecular Polarizable POtential) force field is based upon the model electron density of a hydrogen-like atom. This framework is used to derive and parameterize individual terms describing charge penetration damped permanent electrostatics, damped polarization, charge transfer, anisotropic Pauli repulsion, and damped dispersion interactions. Initial parameter values were fit to Symmetry Adapted Perturbation Theory (SAPT) energy components for ten water dimer configurations, as well as the radial and angular dependence of the canonical dimer. The SAPT-based parameters were then systematically refined to extend the treatment to water bulk phases. The final HIPPO water model provides a balanced representation of a wide variety of properties of gas phase clusters, liquid water and ice polymorphs, across a range of temperatures and pressures. This water potential yields a rationalization of water structure, dynamics and thermodynamics explicitly correlated with an ab initio energy decomposition, while providing a level of accuracy comparable or superior to previous polarizable atomic multipole force fields. The HIPPO water model serves as a cornerstone around which similarly detailed physics-based models can be developed for additional molecular species.
A common way to evaluate electronic integrals for polyatomic molecules is to use Beckes partitioning scheme [J. Chem. Phys.88, 2547 (1988)] in conjunction with overlapping grids centered at each atomic site. The Becke scheme was designed for integrands that fall off rapidly at large distances, such as those approximating bound electronic states. When applied to states in the electronic continuum, however, Becke scheme exhibits slow convergence and it is highly redundant. Here, we present a modified version of Becke scheme that is applicable to functions of the electronic continuum, such as those involved in molecular photoionization and electron-molecule scattering, and which ensures convergence and efficiency comparable to those realized in the calculation of bound states. In this modified scheme, the atomic weights already present in Beckes partition are smoothly switched off within a range of few bond lengths from their respective nuclei, and complemented by an asymptotically unitary weight. The atomic integrals are evaluated on small spherical grids, centered on each atom, with size commensurate to the support of the corresponding atomic weight. The residual integral of the interstitial and long-range region is evaluated with a central master grid. The accuracy of the method is demonstrated by evaluating integrals involving integrands containing Gaussian Type Orbitals and Yukawa potentials, on the atomic sites, as well as spherical Bessel functions centered on the master grid. These functions are representative of those encountered in realistic electron-scattering and photoionization calculations in polyatomic molecules.
We construct a density-functional formalism adapted to uniform external magnetic fields that is intermediate between conventional Density Functional Theory and Current-Density Functional Theory (CDFT). In the intermediate theory, which we term LDFT, the basic variables are the density, the canonical momentum, and the paramagnetic contribution to the magnetic moment. Both a constrained-search formulation and a convex formulation in terms of Legendre--Fenchel transformations are constructed. Many theoretical issues in CDFT find simplified analogues in LDFT. We prove results concerning $N$-representability, Hohenberg--Kohn-like mappings, existence of minimizers in the constrained-search expression, and a restricted analogue to gauge invariance. The issue of additivity of the energy over non-interacting subsystems, which is qualitatively different in LDFT and CDFT, is also discussed.
Lithium based deep eutectic solvents (DESs) are excellent candidates for eco-friendly electrolytes in lithium ion batteries. While some of these DES have shown promising results, a clear mechanism of lithium ion transport in DESs is not yet established. This work reports the study on the solvation and transport of lithium in a DES made from lithium perchlorate and acetamide using Molecular Dynamics (MD) simulation and neutron scattering techniques. Based on hydrogen bonding (H-bonding) of acetamide with neighbouring molecules/ions, two states are largely prevalent: 1) acetamide molecules which are H-bonded to lithium ions (~ 36 %) and 2) acetamide molecules that are entirely free (~ 58%). Analysing their stochastic dynamics independently, it is observed that the long-range diffusion of the former is significantly slower than the latter one. This is also validated from the neutron scattering experiment on the same DES system. Further, the analysis the lithium dynamics shows that the diffusion of acetamide molecules in the first category is strongly coupled to that of lithium ions. On an average the lithium ions are H-bonded to ~ 3.2 acetamide molecules in their first solvation. These observations are further bolstered through the analysis of the H-bond correlation function between acetamide and lithium ions, which show that ~ 90% of lithium ionic transport is achieved by vehicular motion where the ions diffuse along with its first solvation shell. The findings of this work are an important advancement in understanding solvation and transport of lithium ion in DES.
The energies of molecular excited states arise as solutions to the electronic Schr{o}dinger equation and are often compared to experiment. At the same time, nuclear quantum motion is known to be important and to induce a red-shift of excited state energies. However, it is thus far unclear whether incorporating nuclear quantum motion in molecular excited state calculations leads to a systematic improvement of their predictive accuracy, making further investigation necessary. Here we present such an investigation by employing two first-principles methods for capturing the effect of quantum fluctuations on excited state energies, which we apply to the Thiel set of organic molecules. We show that accounting for zero-point motion leads to much improved agreement with experiment, compared to `static calculations which only account for electronic effects, and the magnitude of the red-shift can become as large as 1.36 eV. Moreover, we show that the effect of nuclear quantum motion on excited state energies largely depends on the molecular size, with smaller molecules exhibiting larger red-shifts. Our methodology also makes it possible to analyze the contribution of individual vibrational normal modes to the red-shift of excited state energies, and in several molecules we identify a limited number of modes dominating this effect. Overall, our study provides a foundation for systematically quantifying the shift of excited state energies due to nuclear quantum motion, and for understanding this effect at a microscopic level.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
