Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAb initio molecular dynamics study of collective excitations in liquid H$_2$O and D$_2$O: Effect of dispersion corrections
86
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The collective dynamics in liquid water is an active research topic experimentally, theoretically and via simulations. Here, ab initio molecular dynamics simulations are reported in heavy and ordinary water at temperature 323.15 K, or 50$^circ$C. The simulations in heavy water were performed both with and without dispersion corrections. We found that the dispersion correction (DFT-D3) changes the relaxation of density-density time correlation functions from a slow, typical of a supercooled state, to exponential decay behaviour of regular liquids. This implies an essential reduction of the melting point of ice in simulations with DFT-D3. Analysis of longitudinal (L) and transverse (T) current spectral functions allowed us to estimate the dispersions of acoustic and optic collective excitations and to observe the L-T mixing effect. The dispersion correction shifts the L and T optic (O) modes to lower frequencies and provides by almost thirty per cent smaller gap between the longest-wavelength LO and TO excitations, which can be a consequence of a larger effective high-frequency dielectric permittivity in simulations with dispersion corrections. Simulation in ordinary water with the dispersion correction results in frequencies of optic excitations higher than in D$_2$O, and in a long-wavelength LO-TO gap of 24 ps$^{-1}$ (127 cm$^{-1}$).
rate research
Read More
The impact of the inner structure and thermal history of planets on their observable features, such as luminosity or magnetic field, crucially depends on the poorly known heat and charge transport properties of their internal layers. The thermal and electric conductivities of different phases of water (liquid, solid, and super-ionic) occurring in the interior of ice giant planets, such as Uranus or Neptune, are evaluated from equilibrium ab initio molecular dynamics, leveraging recent progresses in the theory and data analysis of transport in extended systems. The implications of our findings on the evolution models of the ice giants are briefly discussed
We present the results of first-principles molecular-dynamics simulations of molten silicates, based on the density functional formalism. In particular, the structural properties of a calcium aluminosilicate $ [$ CaO-Al$_2$O$_3$-SiO$_2$ $ ]$ melt are compared to those of a silica melt. The local structures of the two melts are in good agreement with the experimental understanding of these systems. In the calcium aluminosilicate melt, the number of non-bridging oxygens found is in excess of the number obtained from a simple stoichiometric prediction. In addition, the aluminum avoidance principle, which states that links between AlO$_4$ tetrahedra are absent or rare, is found to be violated. Defects such as 2-fold rings and 5-fold coordinated silicon atoms are found in comparable proportions in both liquids. However, in the calcium aluminosilicate melt, a larger proportion of oxygen atoms are 3-fold coordinated. In addition, 5-fold coordinated aluminum atoms are observed. Finally evidence of creation and anihilation of non-bridging oxygens is observed, with these oxygens being mostly connected to Si tetrahedra.
Bi-based cuprate superconductors are important materials for both fundamental research and applications. As in other cuprates, the superconducting phase in the Bi compounds lies close to an antiferromagnetic phase. Our density functional theory calculations based on the strongly-constrained-and-appropriately-normed (SCAN) exchange correlation functional in Bi$_2$Sr$_2$CaCu$_2$O$_{8+delta}$ reveal the persistence of magnetic moments on the copper ions for oxygen concentrations ranging from the pristine phase to the optimally hole-doped compound. We also find the existence of ferrimagnetic solutions in the heavily doped compounds, which are expected to suppress superconductivity.
From a detailed study, including polarization dependence, of the normal state angle-resolved photoemission spectra for Bi$_2$Sr$_2$CaCu$_2$O$_8$, we find only one CuO$_2$ band related feature. All other spectral features can be ascribed either to umklapps from the superlattice or to ``shadow bands. Even though the dispersion of the peaks looks like band theory, the lineshape is anomalously broad and no evidence is found for bilayer splitting. We argue that the ``dip feature in the spectrum below $T_c$ arises not from bilayer splitting, but rather from many body effects.
We present a general computational scheme based on molecular dynamics (m.d.) simulation for calculating the chemical potential of adsorbed molecules in thermal equilibrium on the surface of a material. The scheme is based on the calculation of the mean force in m.d. simulations in which the height of a chosen molecule above the surface is constrained, and subsequent integration of the mean force to obtain the potential of mean force and hence the chemical potential. The scheme is valid at any coverage and temperature, so that in principle it allows the calculation of the chemical potential as a function of coverage and temperature. It avoids all statistical mechanical approximations, except for the use of classical statistical mechanics for the nuclei, and assumes nothing in advance about the adsorption sites. From the chemical potential, the absolute desorption rate of the molecules can be computed, provided the equilibration rate on the surface is faster than the desorption rate. We apply the theory by {em ab initio} m.d. simulation to the case of H$_2$O on MgO (001) in the low-coverage limit, using the Perdew-Burke-Ernzerhof (PBE) form of exchange-correlation. The calculations yield an {em ab initio} value of the Polanyi-Wigner frequency prefactor, which is more than two orders of magnitude greater than the value of $10^{13}$ s$^{-1}$ often assumed in the past. Provisional comparison with experiment suggests that the PBE adsorption energy may be too low, but the extension of the calculations to higher coverages is needed before firm conclusions can be drawn. The possibility of including quantum nuclear effects by using path-integral simulations is noted.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
