van der Waals corrections to density functional theory calculations: Methane,ethane, ethylene,benzene, formaldehyde,ammonia, water,PBE, and CPMD

Parameters are developed for a practical application of the empirical van der Waals (vdW) correction infrastructure available in the CPMD density functional theory (DFT) code. The binding energy, geometry, and potential energy surface (PES) are examined for methane, ethane, ethylene, formaldehyde, ammonia, three benzene dimer geometries, and three benzene–water geometries. The vdW corrected results compare favorably with MP2 and CCSD(T) calculations near the complete basis set limits, and with experimental results where they are available.     

水、甲醇和乙醇液体微结构性质的Car-Parrinello 分子动力学模拟

采用Car-Parrinello 分子动力学(CPMD)方法分别研究了水、甲醇和乙醇的液体微结构性质.研究结果显示:在水、甲醇和乙醇三个体系中O…O径向分布函数曲线的第一个峰位置分别为0.278、0.276 和0.275nm; O…H径向分布函数曲线的第一个峰位置分别为0.178、0.176和0.177 nm.表明基团(氢原子、甲基、乙基)的差异对O…O第一个峰的位置影响很小.但基团的差异对径向分布函数峰高的影响却很显著,由水到乙醇第一个峰的高度逐渐变高.空间分布函数表明氧原子和氢原子在溶剂分子周围有取向地分布,这与径向分布函数所表现出尖锐的第一个峰相一致.氢键分布分析显示,水、甲醇和乙醇的平均氢键数分别为3.62、1.99 和1.87,表明水形成了网状氢键结构,而甲醇、乙醇形成链状氢键结构.     

CPMD/GULP QM/MM interface for modeling periodic solids: Implementation and its application in the study of Y‐zeolite supported Rhn clusters

We report here the development of hybrid quantum mechanics/molecular mechanics (QM/MM) interface between the plane‐wave density functional theory based CPMD code and the empirical force‐field based GULP code for modeling periodic solids and surfaces. The hybrid QM/MM interface is based on the electrostatic coupling between QM and MM regions. The interface is designed for carrying out full relaxation of all the QM and MM atoms during geometry optimizations and molecular dynamics simulations, including the boundary atoms. Both Born–Oppenheimer and Car–Parrinello molecular dynamics schemes are enabled for the QM part during the QM/MM calculations. This interface has the advantage of parallelization of both the programs such that the QM and MM force evaluations can be carried out in parallel to model large systems. The interface program is first validated for total energy conservation and parallel scaling performance is benchmarked. Oxygen vacancy in α‐cristobalite is then studied in detail and the results are compared with a fully QM calculation and experimental data. Subsequently, we use our implementation to investigate the structure of rhodium cluster (Rh_n; n = 2 to 6) formed from Rh(C₂H₄)₂ complex adsorbed within a cavity of Y‐zeolite in a reducible atmosphere of H₂ gas. © 2016 Wiley Periodicals, Inc.     

Be2+在水、甲醇和乙醇中结构性质的从头算分子动力学模拟

采用Car-Parrinello分子动力学(CPMD)方法分别研究了Be²⁺在水、甲醇和乙醇中的溶剂结构性质, 并对Be²⁺的第一溶剂壳结构的实验及理论结果进行了比较. 所得第一溶剂壳结构与已报道的实验和理论结果较为一致. 对径向分布函数、配位数以及角度分布等进行了详细的分析. 结果表明: 在水、甲醇和乙醇中, Be²⁺第一溶剂壳为稳定理想的四面体结构. 在本文的模拟时间尺度内,没有观察到第一溶剂壳中的分子与第二溶剂壳中的分子进行交换, 进一步证明Be²⁺第一溶剂壳为稳定的四配位结构. 根据计算得到的空间分布函数, Be²⁺在溶剂分子的等高面上主要集中分布在溶剂分子接受氢键的方向. 根据氧原子在Be²⁺周围的分布, 壳层分子主要集中分布在Be²⁺周围的四个区域, 进一步证实溶剂壳为四面体对称.     

Naphthazarin Derivatives in the Light of Intra- and Intermolecular Forces

Our long-term investigations have been devoted the characterization of intramolecular hydrogen bonds in cyclic compounds. Our previous work covers naphthazarin, the parent compound of two systems discussed in the current work: 2,3-dimethylnaphthazarin (1) and 2,3-dimethoxy-6-methylnaphthazarin (2). Intramolecular hydrogen bonds and substituent effects in these compounds were analyzed on the basis of Density Functional Theory (DFT), Møller–Plesset second-order perturbation theory (MP2), Coupled Clusters with Singles and Doubles (CCSD) and Car-Parrinello Molecular Dynamics (CPMD). The simulations were carried out in the gas and crystalline phases. The nuclear quantum effects were incorporated a posteriori using the snapshots taken from ab initio trajectories. Further, they were used to solve a vibrational Schrödinger equation. The proton reaction path was studied using B3LYP, ωB97XD and PBE functionals with a 6-311++G(2d,2p) basis set. Two energy minima (deep and shallow) were found, indicating that the proton transfer phenomena could occur in the electronic ground state. Next, the electronic structure and topology were examined in the molecular and proton transferred (PT) forms. The Atoms In Molecules (AIM) theory was employed for this purpose. It was found that the hydrogen bond is stronger in the proton transferred (PT) forms. In order to estimate the dimers’ stabilization and forces responsible for it, the Symmetry-Adapted Perturbation Theory (SAPT) was applied. The energy decomposition revealed that dispersion is the primary factor stabilizing the dimeric forms and crystal structure of both compounds. The CPMD results showed that the proton transfer phenomena occurred in both studied compounds, as well as in both phases. In the case of compound 2, the proton transfer events are more frequent in the solid state, indicating an influence of the environmental effects on the bridged proton dynamics. Finally, the vibrational signatures were computed for both compounds using the CPMD trajectories. The Fourier transformation of the autocorrelation function of atomic velocity was applied to obtain the power spectra. The IR spectra show very broad absorption regions between 700 cm−1–1700 cm−1 and 2300 cm−1–3400 cm−1 in the gas phase and 600 cm−1–1800 cm−1 and 2200 cm−1–3400 cm−1 in the solid state for compound 1. The absorption regions for compound 2 were found as follows: 700 cm−1–1700 cm−1 and 2300 cm−1–3300 cm−1 for the gas phase and one broad absorption region in the solid state between 700 cm−1 and 3100 cm−1. The obtained spectroscopic features confirmed a strong mobility of the bridged protons. The inclusion of nuclear quantum effects showed a stronger delocalization of the bridged protons.     

Intermolecular Interactions and Spectroscopic Signatures of the Hydrogen-Bonded System—n-Octanol in Experimental and Theoretical Studies

n-Octanol is the object of experimental and theoretical study of spectroscopic signatures and intermolecular interactions. The FTIR measurements were carried out at 293 K for n-octanol and its deuterated form. Special attention was paid to the vibrational features associated with the O-H stretching and the isotope effect. Density Functional Theory (DFT) in its classical formulations was applied to develop static models describing intermolecular hydrogen bond (HB) and isotope effect in the gas phase and using solvent reaction field reproduced by Polarizable Continuum Model (PCM). The Atoms in Molecules (AIM) theory enabled electronic structure and molecular topology study. The Symmetry-Adapted Perturbation Theory (SAPT) was used for energy decomposition in the dimers of n-octanol. Finally, time-evolution methods, namely classical molecular dynamics (MD) and Car-Parrinello Molecular Dynamics (CPMD) were employed to shed light onto dynamical nature of liquid n-octanol with emphasis put on metric and vibrational features. As a reference, CPMD gas phase results were applied. Nuclear quantum effects were included using Path Integral Molecular Dynamics (PIMD) and a posteriori method by solving vibrational Schrödinger equation. The latter applied procedure allowed to study the deuterium isotope effect.     

Revealing Intra- and Intermolecular Interactions Determining Physico-Chemical Features of Selected Quinolone Carboxylic Acid Derivatives

The intra- and intermolecular interactions of selected quinolone carboxylic acid derivatives were studied in monomers, dimers and crystals. The investigated compounds are well-recognized as medicines or as bases for further studies in drug design. We employed density functional theory (DFT) in its classical formulation to develop gas-phase and solvent reaction field (PCM) models describing geometric, energetic and electronic structure parameters for monomers and dimers. The electronic structure was investigated based on the atoms in molecules (AIM) and natural bond orbital (NBO) theories. Special attention was devoted to the intramolecular hydrogen bonds (HB) present in the investigated compounds. The characterization of energy components was performed using symmetry-adapted perturbation theory (SAPT). Finally, the time-evolution methods of Car–Parrinello molecular dynamics (CPMD) and path integral molecular dynamics (PIMD) were employed to describe the hydrogen bond dynamics as well as the spectroscopic signatures. The vibrational features of the O-H stretching were studied using Fourier transformation of the autocorrelation function of atomic velocity. The inclusion of quantum nuclear effects provided an accurate depiction of the bridged proton delocalization. The CPMD and PIMD simulations were carried out in the gas and crystalline phases. It was found that the polar environment enhances the strength of the intramolecular hydrogen bonds. The SAPT analysis revealed that the dispersive forces are decisive factors in the intermolecular interactions. In the electronic ground state, the proton-transfer phenomena are not favourable. The CPMD results showed generally that the bridged proton is localized at the donor side, with possible proton-sharing events in the solid-phase simulation of stronger hydrogen bridges. However, the PIMD enabled the quantitative estimation of the quantum effects inclusion—the proton position was moved towards the bridge midpoint, but no qualitative changes were detected. It was found that the interatomic distance between the donor and acceptor atoms was shortened and that the bridged proton was strongly delocalized.