Cu2+ in liquid ammonia-The impact of solvent flexibility and electron correlation in ab initio quantum mechanical charge field molecular dynamics

The impact of solvent flexibility and electron correlation on the simulation results of Cu²⁺ in liquid ammonia has been investigated via an ab initio quantum mechanical charge field molecular dynamics (QMCF MD) simulation approach. To achieve this, three different simulation systems were considered in this study, namely Cu²⁺ in rigid and flexible ammonia at Hartree-Fock (HF) level of theory, as well as resolution of identity second order Møller-Plesset (MP2) perturbation theory in the rigid body case. In all cases, a stable octahedral [Cu(NH₃ )₆ ]²⁺ complex subject to dynamic Jahn-Teller distortions without the occurrence of ligand exchange was observed. The Cu²⁺  - NH₃ distance in the first shell agrees well with the experimental and other theoretical data. In all three cases, the structural data shows that the rigid-body ammonia model in conjunction with the HF level of theory provides accurate data for the first solvation shell, while at the same time, the computational demand and thus the achievable simulation time are much more beneficial. The vibrational analysis of the Cu²⁺  - NH₃ interaction yields similar force constants in the three investigated systems indicating that there is no distinct difference on the dynamical properties of the first solvation shell. In addition to the QMCF MD simulations, a number of natural bond orbital (NBO) analyses were carried out, confirming the strong electrostatic character of the Cu²⁺  - NH₃ interaction.

Hydration Properties of HnPO4n-3 (n = 0-3) From Ab Initio Molecular Dynamics Simulations

For a comprehensive and detailed microscopic understanding of the hydration properties of primary aqueous phosphorus species of valence states V (viz., H₃PO₄, H₂PO₄^-, HPO₄^2-, and PO₄^3-), a series of extensive ab initio molecular dynamics simulations is conducted at ambient temperature. In each of these cases, the spatially resolved, three-dimensional hydration shells are computed, allowing for a direct microscopic visual understanding of the hydration shells around the species. Since these species are excellent agents for the formation of hydrogen bonds (H-bonds) in water, which determine a wide range of their structural, dynamic, and spectroscopic features, a detailed analysis of the qualitative and quantitative aspects of the H-bonds, including their lifetime calculations, is performed. Vibrational density of states (VDOS) is calculated for each of the species in solute phases, resolved for each H-bonding site, and compared against the gas-phase normal modes of H₃PO₄ for the purpose of understanding the signatures of the peaks in VDOS plots and, in particular, the effects of solvation and H-bonding mechanisms. The results are well in line with available experimental data and other recent computer-aided studies in the literature.

The Jahn-Teller effect in mixed aqueous solution: the solvation of Cu2+ in 18.6% aqueous ammonia obtained from ab initio quantum mechanical charge field molecular dynamics

The solvation structure and dynamics of Cu2+ in 18.6 % aqueous ammonia have been investigated using an ab initio quantum mechanical charge field molecular dynamics (QMCF MD) simulation approach at the Hartree–Fock (HF) level of theory applying the LANL2DZ ECP and Dunning DZP basis sets for Cu2+, ammonia and water, respectively. During a simulation time of 20 ps, only NH3 molecules are observed within the first solvation shell of Cu2+, resulting in the formation of an octahedral [Cu(NH3)6]2+ complex. While no exchange of these ligands with the second solvation shell are observed along the simulation, the monitoring of the associated N-Ntrans distances highlight the dynamics of the associated Jahn-Teller distortions, showing on average 2 elongated axial (2.19 Å) and 4 equatorial Cu–N bonds (2.39 Å). The observed structural properties are found in excellent agreement with experimental studies. In addition, an NBO analysis was carried out, confirming the strong electrostatic character of the Cu2+–NH3 interaction.

Remarkable conformational flexibility of aqueous 18-crown-6 and its strontium(II) complex-ab initio molecular dynamics simulations

Ab initio QMCF-MD simulations of aqueous 18-crown-6 (18C6) and strontium(II)-18-crown-6 (18C6-Sr) were performed to gather insight into their hydration properties. Strongly different characteristics were found for the two solutes in terms of structure and dynamics such as H-bonding. They, however, have in common that their backbone shows high flexibility in aqueous medium, adopting structures significantly differing from idealized gas phase geometries. In particular, planar oxyethylene units stable in the picosecond range occurred in 18C6, while the strontium complex readily exhibits a bent structure. Detailed analysis of this high flexibility was done via two dimensional root mean square deviation plots as well as the evolution of dihedral angles and angles within the simulation trajectory. The vibrational spectra obtained from the QMCF-MD simulations are in excellent agreement with experimental data and show a pronounced blueshift upon complexation of the strontium(II) ion in 18C6.

Hydration of guanidinium depends on its local environment

Hydration of gaseous guanidinium (Gdm+) with up to 100 water molecules attached was investigated using infrared photodissociation spectroscopy in the hydrogen stretch region between 2900 and 3800 cm-1. Comparisons to IR spectra of low-energy computed structures indicate that at small cluster size, water interacts strongly with Gdm+ with three inner shell water molecules each accepting two hydrogen bonds from adjacent NH2 groups in Gdm+. Comparisons to results for tetramethylammonium (TMA+) and Na+ enable structural information for larger clusters to be obtained. The similarity in the bonded OH region for Gdm(H2O)20+vs. Gdm(H2O)100+ and the similarity in the bonded OH regions between Gdm+ and TMA+ but not Na+ for clusters with <50 water molecules indicate that Gdm+ does not significantly affect the hydrogen-bonding network of water molecules at large size. These results indicate that the hydration around Gdm+ changes for clusters with more than about eight water molecules to one in which inner shell water molecules only accept a single H-bond from Gdm+. More effective H-bonding drives this change in inner-shell water molecule binding to other water molecules. These results show that hydration of Gdm+ depends on its local environment, and that Gdm+ will interact with water even more strongly in an environment where water is partially excluded, such as the surface of a protein. This enhanced hydration in a limited solvation environment may provide new insights into the effectiveness of Gdm+ as a protein denaturant.

Structure and dynamics of chromatographically relevant Fe(III)-chelates

Immobilized metal ion affinity chromatography (IMAC) is an important chromatographic technique for biomolecules. In order to get a detailed understanding of the hydration of immobilized Fe(III), complexes of Fe(III) with methyl substituted iminodiacetate ([Fe(MSIDA)(H2O)3](+)) as well as with methyl substituted nitrilotriacetate ([Fe(MSNTA)(H2O)2]) were simulated in aqueous solutions with the quantum mechanical charge field molecular dynamics (QMCF MD) approach. The simulations were carried out at the Hartree-Fock (HF) level of theory, since cluster calculations at the HF, MP2, and B3LYP levels of theory showed that this method results in a good compromise between computational effort and accuracy. None of the coordinating water molecules were exchanged during the simulation period of 15 ps. The Fe-OH2O bond distances as well as the Fe-OH2O stretching motions differed among the coordinating water molecules, indicating different bond strengths. For the water molecules in the second hydration layer, mean residence times of 2.7 and 1.9 ps were obtained for [Fe(MSIDA)(H2O)3](+) and [Fe(MSNTA)(H2O)2], respectively. Furthermore, infrared measurements were carried out to characterize the most prominent bond features of aqueous Fe(III)-NTA and to discuss these results in conjunction with the computationally derived frequencies.

Perspectives for hybrid ab initio/molecular mechanical simulations of solutions: from complex chemistry to proton-transfer reactions and interfaces

As a consequence of the ongoing development of enhanced computational resources, theoretical chemistry has become an increasingly valuable field for the investigation of a variety of chemical systems. Simulations employing a hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) technique have been shown to be a particularly promising approach, whenever ultrafast (i.e., picosecond) dynamical properties are to be studied, which are in many cases difficult to access via experimental techniques. Details of the quantum mechanical charge field (QMCF) ansatz, an advanced QM/MM protocol, are discussed and simulation results for various systems ranging from simple ionic hydrates to solvated organic molecules and coordination complexes in solution are presented. A particularly challenging application is the description of proton-transfer reactions in chemical simulations, which is a prerequisite to study acidified and basic systems. The methodical requirements for a combination of the QMCF methodology with a dissociative potential model for the description of the solvent are discussed. Furthermore, the possible extension of QM/MM approaches to solid/liquid interfaces is outlined.

Combined Ab Initio Computational and Infrared Spectroscopic Study of the cis- and trans-Bis(glycinato)copper(II) Complexes in Aqueous Environment

The cis- and trans-bis(glycinato)copper(II) complexes in aqueous solution have been investigated by means of a combined theoretical and experimental approach. The conducted quantum mechanical charge field molecular dynamics (QMCF-MD) studies, being the first quantum mechanical simulations of organometallic complexes by this method, yielded accurate structural details of the investigated isomers as well as novel dynamic data, which has successfully been confirmed and extended by subsequent mid-infrared measurements. The spectroscopic results, critically assessed by adjacent multivariate data analysis, indicate an isomeric stability at ambient conditions, vanishing at elevated temperatures.