Hydrated germanium (II): Irregular structural and dynamical properties revealed by a quantum mechanical charge field molecular dynamics study

Structural and dynamical properties of Ge (II) in aqueous solution have been investigated using the novel ab initio quantum mechanical charge field (QMCF) molecular dynamics (MD) formalism. The first and second hydration shells were treated by ab initio quantum mechanics at restricted Hartree–Fock (RHF) level using the cc‐pVDZ‐PP basis set for Ge (II) and Dunning double‐ζ plus polarization basis sets for O and H. Besides ligand exchange processes and mean ligand residence times to observe dynamics, tilt‐ and theta‐angle distributions along with an advanced structural parameter, namely radial and angular distribution functions (RAD) for different regions were also evaluated. The combined radial and angular distribution depicted through surface plot and contour map is presented to provide a detailed insight into the density distribution of water molecules around the Ge²⁺ ion. A strongly distorted hydration structure with two trigonal pyramidal substructures within the first hydration shell is observed, which demonstrates the lone‐pair influence and provides a new basis for the interpretation of the catalytic and pharmacological properties of germanium coordination compounds. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

An ab initio quantum mechanical charge field molecular dynamics simulation of a dilute aqueous HCl solution

An ab initio quantum mechanical charge field (QMCF) molecular dynamics simulation has been performed to study the structural and dynamical properties of a dilute aqueous HCl solution. The solute molecule HCl and its surrounding water molecules were treated at Hartree‐Fock level in conjunction with Dunning double‐ζ plus polarization function basis sets. The simulation predicts an average H? Cl bond distance of 1.28 Å, which is in good agreement with the experimental value. The H_HCl···O_w and Cl_HCl···H_w distances of 1.84 and 3.51 Å were found for the first hydration shell. At the hydrogen site of HCl, a single water molecule is the most preferred coordination, whereas an average coordination number of 12 water molecules of the full first shell was observed for the chloride site. The hydrogen bonding at the hydrogen site of HCl is weakened by proton transfer reactions and an associated lability of ligand binding. Two proton transfer processes were observed in the QMCF MD simulation, demonstrating acid dissociation of HCl. A weak structure‐making/breaking effect of HCl in water is recognized from the mean residence times of 2.1 and 0.8 ps for ligands in the neighborhood of Cl and H sites of HCl, respectively. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Quantum mechanical charge field molecular dynamics simulation of the TiO2+ ion in aqueous solution

Structural and dynamical properties of the TiO(2+) ion in aqueous solution have been investigated by using the new ab initio quantum mechanical charge field (QMCF) molecular dynamics (MD) formalism, which does not require any other potential functions except those for solvent-solvent interactions. Both first and second hydration shell have been treated at Hartree-Fock (HF) quantum mechanical level. A Ti-O bond distance of 1.5 A was observed for the [Ti=O](2+) ion. The first hydration shell of the ion shows a varying coordination number ranging from 5 to 7, five being the dominant one and representing one axial and four equatorial water molecules directly coordinated to Ti, which are located at 2.3 A and 2.1 A, respectively. The flexibility in the coordination number reflects the fast exchange processes, which occur only at the oxo atom, where water ligands are weakly bound through hydrogen bonds. Considering the first shell hydration, the composition of the TiO(2+) hydrate can be characterized as [(H(2)O)(0.7)(H(2)O)(4) (eq)(H(2)O)(ax)](2+). The second shell consists in average of 12 water molecules located at a mean distance of 4.4 A. Several other structural parameters such as radial and angular distribution functions and coordination number distributions were analyzed to fully characterize the hydration structure of the TiO(2+) ion in aqueous solution. For the dynamics of the TiO(2+) ion, different sets of dynamical parameters such as Ti=O, Ti-O(eq), and Ti-O(ax) stretching frequencies and ligands' mean residence times were evaluated. During the simulation time of 15 ps, 3 water exchange processes in the first shell were observed at the oxo atom, corresponding to a mean residence time of 3.6 ps. The ligands' mean residence time for the second shell was determined as 3.5 ps.     

Ab Initio quantum mechanical charge field study of hydrated bicarbonate ion: Structural and dynamical properties

The ab initio quantum mechanical charge field molecular dynamics (QMCF MD) formalism was applied to simulate the bicarbonate ion, HCO₃^?, in aqueous solution. The difference in coordination numbers obtained by summation over atoms (6.6) and for the solvent‐accessible surface (5.4) indicates the sharing of some water molecules between the individual atomic hydration shells. It also proved the importance to consider the hydration of the chemically different atoms individually for the evaluation of structural and dynamical properties of the ion. The orientation of water molecules in the hydration shell was visualized by the θ–tilt surface plot. The mean residence time in the surroundings of the HCO₃^? ion classify it generally as a structure‐breaking ion, but the analysis of the individual ion‐water hydrogen bonds revealed a more complex behavior of the different coordination sites. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Temperature dependence of structure and dynamics of the hydrated Ca2+ ion according to ab initio quantum mechanical charge field and classical molecular dynamics

Simulations using ab initio quantum mechanical charge field molecular dynamics (QMCF MD) and classical molecular dynamics using two‐body and three‐body potentials were performed to investigate the hydration of the Ca²⁺ ion at different temperatures. Results from the simulations demonstrate significant effects of temperature on solution dynamics and the corresponding composition and structure of hydrated Ca²⁺. Substantial increase in ligand exchange events was observed in going from 273.15 K to 368.15 K, resulting in a redistribution of coordination numbers to lower values. The effect of temperature is also visible in a red‐shift of the ion‐oxygen stretching frequencies, reflecting weakened ligand binding. Even the moderate increase from ambient to body temperature leads to significant changes in the properties of Ca²⁺ in aqueous environment. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Sulfur dioxide in water: structure and dynamics studied by an ab initio quantum mechanical charge field molecular dynamics simulation

An ab initio Quantum Mechanical Charge Field Molecular Dynamics Simulation (QMCF MD) was performed to investigate structure and dynamics behavior of hydrated sulfur dioxide (SO(2)) at the Hartree-Fock level of theory employing Dunning DZP basis sets for solute and solvent molecules. The intramolecular structural characteristics of SO(2), such as S═O bond lengths and O═S═O bond angle, are in good agreement with the data available from a number of different experiments. The structural features of the hydrated SO(2) were primarily evaluated in the form of S-O(wat) and O(SO(2))-H(wat) radial distribution functions (RDFs) which gave mean distances of 2.9 and 2.2 ?, respectively. The dynamical behavior characterizes the solute molecule to have structure making properties in aqueous solution or water aerosols, where the hydrated SO(2) can easily get oxidized to form a number of sulfur(VI) species, which are believed to play an important role in the atmospheric processes.     

Structure and dynamics of phosphate ion in aqueous solution: an ab initio QMCF MD study

A simulation of phosphate in aqueous solution was carried out employing the new QMCF MD approach which offers the possibility to investigate composite systems with the accuracy of a QMMM method but without the time consuming creation of solute-solvent potential functions. The data of the simulations give a clear picture of the hydration shells of the phosphate anion. The first shell consists of 13 water molecules and each oxygen of the phosphate forms in average three hydrogens bonds to different solvent molecules. Several structural parameters such as radial distribution functions and coordination number distributions allow to fully characterize the embedding of the highly charged phosphate ion in the solvent water. The dynamics of the hydration structure of phosphate are described by mean residence times of the solvent molecules in the first hydration shell and the water exchange rate.     

Structure and ultrafast dynamics of liquid water: a quantum mechanics/molecular mechanics molecular dynamics simulations study

A quantum mechanics/molecular mechanics molecular dynamics simulation was performed for liquid water to investigate structural and dynamical properties of this peculiar liquid. The most important region containing a central reference molecule and all nearest surrounding molecules (first coordination shell) was treated by Hartree-Fock (HF), post-Hartree-Fock [second-order Moller-Plesset perturbation theory (MP2)], and hybrid density functional B3LYP [Becke's three parameter functional (B3) with the correlation functional of Lee, Yang, and Parr (LYP)] methods. In addition, another HF-level simulation (2HF) included the full second coordination shell. Site to site interactions between oxygen-oxygen, oxygen-hydrogen, and hydrogen-hydrogen atoms of all ab initio methods were compared to experimental data. The absence of a second peak and the appearance of a shoulder instead in the gO-O graph obtained from the 2HF simulation is notable, as this feature has been observed so far only for pressurized or heated water. Dynamical data show that the 2HF procedure compensates some of the deficiency of the HF one-shell simulation, reducing the difference between correlated (MP2) and HF results. B3LYP apparently leads to too rigid structures and thus to an artificial slow down of the dynamics.     

Quantum mechanical/molecular mechanical simulations of the Tl(III) ion in water

Classical molecular dynamics (MD) and combined quantum mechanical/molecular mechanical (QM/MM) MD simulations have been performed to investigate the structural and dynamical properties of the Tl(III) ion in water. A six-coordinate hydration structure with a maximum probability of the Tl-O distance at 2.21 A was observed, which is in good agreement with X-ray data. The librational and vibrational spectra of water molecules in the first hydration shell are blue-shifted compared with those of pure liquid water, and the Tl-O stretching force constant was evaluated as 148 Nm(-1). Both structural and dynamical properties show a distortion of the first solvation shell structure. The second shell ligands' mean residence time was determined as 12.8 ps. The Tl(III) ion can be classified as "structure forming" ion; the calculated hydration energy of -986 +/- 9 kcal mol agrees well with the experimental value of -986 kcal mol.     

Making and breaking of small water clusters: A combined quantum chemical and molecular dynamics approach

We present a combined quantum chemical and molecular dynamics study of cyclic and noncyclic water n-mers ([(H₂O]_n, n = 2–6) at four different temperatures and showcase that the dynamics of small water clusters can reproduce the known properties of bulk water reasonably well. We investigate the making and breaking of the water clusters by computing the hydrogen bond strengths, average lifetimes, and relative stabilities, which are important to understand the complex solution dynamics. We compare the behavior of water clusters in the gas phase and in the solution phase as well as the variation in the properties as a function of cluster size and highlight the notably more interesting cluster dynamics of the water trimer when compared to the other water clusters. © 2019 Wiley Periodicals, Inc.     

Structural and dynamical properties of hydrogen fluoride in aqueous solution: an ab initio quantum mechanical charge field molecular dynamics simulation

The novel ab initio quantum mechanical charge field (QMCF) molecular dynamics simulation at the Hartree-Fock level has been employed to investigate hydration structure and dynamics of hydrogen fluoride in aqueous solution. The average H-F bond length of 0.93 A obtained from the QMCF MD simulation is in good agreement with the experimental data. The HHF...Ow distance of 1.62 A was evaluated for the first hydration shell, and 2.00 A was observed for the FHF...Hw distance. The stability of hydrogen bonding is more pronounced in the hydrogen site of hydrogen fluoride, with a single water molecule in this part of the first hydration shell. A wide range of coordination numbers between 3 and 9 with an average value of 5.6 was obtained for the fluorine site. The force constants of 819.1 and 5.9 N/m were obtained for the HHF-FHF and HHF...Ow interactions, respectively, proving the stability of the nondissociated form of hydrogen fluoride in aqueous solution. The mean residence times of 2.1 and 2.5 ps were determined for ligand exchange processes in the neighborhood of fluorine and hydrogen atoms of hydrogen fluoride, respectively, indicating a weak structure-making effect of hydrogen fluoride in water. The corresponding H-bond lifetimes attribute this effect to the H atom site of HF.     

Structure,thermodynamics and diffusion in asymmetric binary mixtures: a molecular dynamics simulation study

Structural and thermodynamic properties as well as diffusion coefficients of binary fluid mixtures with asymmetry in mass, size, charge and their combinations have been studied using classical molecular dynamics simulations. The fluid mixture is modelled as spherical particles interacting via the Weeks–Chandler–Andersen and Coulomb potential. The diameter, charge and mass of the fluid particles are in the range 6–60 Å, 1–10e and 1—500 amu, respectively. Systematic variations in pair-correlation functions, thermodynamic properties as well as the self-diffusion coefficient are found with the size, charge and mass ratio of the particles. The self-diffusion coefficient for systems having more than one type of asymmetry is calculated and expressed in terms of diffusion coefficients of systems with only one type of asymmetry.     

An experimental and computational study on intramolecular charge transfer: a tetrathiafulvalene-fused dipyridophenazine molecule

To study the electronic interactions in donor-acceptor (D-A) ensembles, D and A fragments are coupled in a single molecule. Specifically, a tetrathiafulvalene (TTF)-fused dipyrido[3,2-a:2',3'-c]phenazine (dppz) compound having inherent redox centers has been synthesized and structurally characterized. Its electronic absorption, fluorescence emission, photoinduced intramolecular charge transfer, and electrochemical behavior have been investigated. The observed electronic properties are explained on the basis of density functional theory.     

Performance assessment of semiempirical molecular orbital methods in describing halogen bonding: quantum mechanical and quantum mechanical/molecular mechanical-molecular dynamics study

The performance of semiempirical molecular-orbital methods--MNDO, MNDO-d, AM1, RM1, PM3 and PM6--in describing halogen bonding was evaluated, and the results were compared with molecular mechanical (MM) and quantum mechanical (QM) data. Three types of performance were assessed: (1) geometrical optimizations and binding energy calculations for 27 halogen-containing molecules complexed with various Lewis bases (Two of the tested methods, AM1 and RM1, gave results that agree with the QM data.); (2) charge distribution calculations for halobenzene molecules, determined by calculating the solvation free energies of the molecules relative to benzene in explicit and implicit generalized Born (GB) solvents (None of the methods gave results that agree with the experimental data.); and (3) appropriateness of the semiempirical methods in the hybrid quantum-mechanical/molecular-mechanical (QM/MM) scheme, investigated by studying the molecular inhibition of CK2 protein by eight halobenzimidazole and -benzotriazole derivatives using hybrid QM/MM molecular-dynamics (MD) simulations with the inhibitor described at the QM level by the AM1 method and the rest of the system described at the MM level. The pure MM approach with inclusion of an extra point of positive charge on the halogen atom approach gave better results than the hybrid QM/MM approach involving the AM1 method. Also, in comparison with the pure MM-GBSA (generalized Born surface area) binding energies and experimental data, the calculated QM/MM-GBSA binding energies of the inhibitors were improved by replacing the G(GB,QM/MM) solvation term with the corresponding G(GB,MM) term.     

Diffusion of methanol in zeolite NaY: a molecular dynamics study

Molecular dynamics simulations were performed in order to obtain a detailed understanding of the self-diffusion mechanisms of methanol in the zeolite NaY system. We derived a new force-field term to describe the interactions between the methanol molecules and the extraframework cations. From the simulations, we show that diffusive behavior in the high-temperature range consists of a combination of both short- and long-range motions at low and intermediate loadings. This type of motion is characterized by an activation energy that decreases as the loading increases. At low loadings, we also observe short-range diffusive behavior based on a surface-mediated mechanism. The short-range behavior corresponds to motion only on the length scale of an FAU supercage, whereas the long-range behavior involves intercage diffusion. For the saturation loading corresponding to 96 methanol molecules per unit cell, only short-range motions within the same supercage predominate. Finally, the preferential arrangement of the adsorbate molecules around the extraframework cations are examined and compared with those previously deduced from experimental data.     

Vibrational dynamics of hydrogen‐bonded HCN complexes with OH and NH acids: Computational DFT systematic study

Vibrational properties (band position, infrared [IR], and Raman intensities) of C?N stretching mode were studied in 65 gas phase hydrogen‐bonded 1:1 complexes of HCN with OH acids and NH acids using density functional theory (DFT) calculations at the B3LYP‐6‐311++G(d,p) level. Furthermore, general characteristics of the hydrogen bonds and vibrational changes in acids OH/NH stretching bands were also considered. Experimentally observed blue shift of the C?N stretching band promoted by hydrogen bonding, which shortens the triple bond length, is very well reproduced and quantitatively depends on the hydrogen bond length. Both IR and Raman ν(C?N) band intensities are enhanced, also in good agreement with the experimental results. IR intensity increase is a direct function of the hydrogen bond energy. However, the predicted Raman intensity raise is a more complex function, depending simultaneously on characteristics of both the hydrogen bond (C?N bond length) and the H‐donating acid (polarizability). With these two parameters, ν (C?N) Raman intensities of the complexes are explained with a mean error of ±2.4%. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Limiting diffusion coefficients of ionic liquids in water and methanol: a combined experimental and molecular dynamics study

Mutual diffusion coefficients D(12) of the ionic liquids 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide ([C(2)MIM][NTf(2)]) and [C(4)MIM][NTf(2)] in highly diluted solutions of water and methanol have been measured at different temperatures between 288 K and 313 K using the Taylor dispersion technique. Tracer diffusion coefficients of the two cations [C(2)MIM](+) and [C(4)MIM](+) as well as the anion [NTf(2)](-) in these solutions have been obtained by molecular dynamics (MD) simulations. For our simulations we used well established force fields for the solvents water and methanol and a recently developed force field for imidazolium-based ionic liquid [C(n)MIM][NTf(2)]. Mutual diffusion coefficients D(12) have been calculated from the tracer diffusion coefficients using the Nernst-Hartley equation strictly valid only at low ionic concentration. The agreement between the diffusion coefficients reported in the literature, the experimental data obtained in this work and the MD results is excellent.     

Ab initio quantum mechanical charge field (QMCF) molecular dynamics: a QM/MM – MD procedure for accurate simulations of ions and complexes

A new formalism for quantum mechanical / molecular mechanical (QM/MM) dynamics of chemical species in solution has been developed, which does not require the construction of any other potential functions except those for solvent–solvent interactions, maintains all the advantages of large simulation boxes and ensures the accuracy of ab initio quantum mechanics for all forces acting in the chemically most relevant region. Interactions between solute and more distant solvent molecules are incorporated by a dynamically adjusted force field corresponding to the actual molecular configuration of the simulated system and charges derived from the electron distribution in the solvate. The new formalism has been tested with some examples of hydrated ions, for which accurate conventional ab initio QM/MM simulations have been previously performed, and the comparison shows equivalence and in some aspects superiority of the new method. As this simulation procedure does not require any tedious construction of two-and three-body interaction potentials inherent to conventional QM/MM approaches, it opens the straightforward access to ab initio molecular dynamics simulations of any kind of solutes, such as metal complexes and other composite species in solution.     

A nonadditive methanol force field: bulk liquid and liquid-vapor interfacial properties via molecular dynamics simulations using a fluctuating charge model

We study the bulk and interfacial properties of methanol via molecular dynamics simulations using a CHARMM (Chemistry at HARvard Molecular Mechanics) fluctuating charge force field. We discuss the parametrization of the electrostatic model as part of the ongoing CHARMM development for polarizable protein force fields. The bulk liquid properties are in agreement with available experimental data and competitive with existing fixed-charge and polarizable force fields. The liquid density and vaporization enthalpy are determined to be 0.809 g/cm3 and 8.9 kcal/mol compared to the experimental values of 0.787 g/cm3 and 8.94 kcal/mol, respectively. The liquid structure as indicated by radial distribution functions is in keeping with the most recent neutron diffraction results; the force field shows a slightly more ordered liquid, necessarily arising from the enhanced condensed phase electrostatics (as evidenced by an induced liquid phase dipole moment of 0.7 D), although the average coordination with two neighboring molecules is consistent with the experimental diffraction study as well as with recent density functional molecular dynamics calculations. The predicted surface tension of 19.66+/-1.03 dyn/cm is slightly lower than the experimental value of 22.6 dyn/cm, but still competitive with classical force fields. The interface demonstrates the preferential molecular orientation of molecules as observed via nonlinear optical spectroscopic methods. Finally, via canonical molecular dynamics simulations, we assess the model's ability to reproduce the vapor-liquid equilibrium from 298 to 423 K, the simulation data then used to obtain estimates of the model's critical temperature and density. The model predicts a critical temperature of 470.1 K and critical density of 0.312 g/cm3 compared to the experimental values of 512.65 K and 0.279 g/cm3, respectively. The model underestimates the critical temperature by 8% and overestimates the critical density by 10%, and in this sense is roughly equivalent to the underlying fixed-charge CHARMM22 force field.     

Structure and dynamics of water surrounding the poly(methacrylic acid): a molecular dynamics study

All-atom molecular dynamics simulations are used to study a single chain of poly(methacrylic acid) in aqueous solutions at various degrees of charge density. Through a combination of analysis on the radial distribution functions of water and snapshots of the equilibrated structure, we observe that local arrangements of water molecules, surrounding the functional groups of COO- and COOH in the chain, behave differently and correlated well to the resulting chain conformation behavior. In general, due to strong attractive interactions between water and charged COO- via the formation of hydrogen bonds, water molecules tend to form shell-like layers around the COO- groups. Furthermore, water molecules often act as a bridging agent between two neighboring COO- groups. These bridged water molecules are observed to stabilize the rodlike chain conformation that the highly charged chain reveals, as they significantly limit torsional and bending degrees of the backbone monomers. In addition, they display different dynamic properties from the bulk water. Both the resulting oxygen and hydrogen spectra are greatly shifted due to the presence of strong H-bonded interactions.