New Insights into the Jahn-Teller Effect through ab initio Quantum-Mechanical/Molecular-Mechanical Molecular Dynamics Simulations of CuII in Water

Structural, Thermodynamic, and Spectroscopic Evolution in the Hydration of Copper(II) Ions, Cu2+(H2O)2-8

Gas-phase clusters of the hydrated Cu(II) cation with 2-8 water molecules were investigated using ab initio quantum chemistry. Isomer structures, energies, and vibrational spectra were computed across this size range, yielding a qualitative picture of this ion as an intact Cu2+ hydrate that also partially oxidizes the surrounding water network at equilibrium. At sufficient cluster sizes, these ion hydrates also become thermodynamically preferred over competitive Cu(II) hydroxide hydrates. Competitive coordination environments were found to exist at some cluster sizes, due to both hydrogen-bonding and d-orbital chemical effects, and the dominant coordination number was found in some cases to be temperature-dependent. Clear spectral signatures of the ion's coordination environment were computed to exist at each cluster size, which should make experimental verification of these computational predictions straightforward. Through comparison to recent studies of hydrated CuOH+, the effective charge on the metal center was shown to converge to approximately +1.5 in both cases, despite qualitatively different behavior of their radical spin densities. Therefore, nominally Cu(II) ions exhibit considerable electronic, chemical, and structural flexibility. The electronic origins of this flexibility─including key roles played by the water network itself─are investigated in this work and should provide a conceptual foundation for future studies of copper-based, water-oxidation catalysts.

Deep purification of copper from Cu(II)-EDTA acidic wastewater by Fe(III) replacement/diethyldithiocarbamate precipitation

Cu(II)-EDTA is a highly stable typical metal-organic complex in a wide pH range (3.0-12.0) and it is difficult to deeply purify Cu(II) by conventional precipitation methods. In this study, Fe(III) replacement/diethyldithiocarbamate (DDTC) precipitation combined process is proposed as a promising strategy to achieve the deep purification of Cu(II) from Cu(II)-EDTA acidic wastewater. The underlying mechanism has also been systematically elucidated by chemical equilibriums, experiments, and density functional theory (DFT) calculations, laying a foundation for the development and application. Chemical equilibriums show that Fe(III) replacement favors the stoichiometric release of Cu(II) from Cu(II)-EDTA and the formation of Fe(III)-EDTA complex under acidic conditions. Experimentally, Cu(II) is removed (over 99.99%) and deeply purified (under 0.008 mg/L) under the optimal conditions, which is lower than the most stringent discharge standards of copper ions in electroplating effluent (<0.5 mg/L, China). DFT calculations reveal that DDTC could further precipitate the released free copper ions via the carbon disulfide (-C(=S)-S) chelating group while exhibiting a slight effect on the Fe(III) in Fe(III)-EDTA. Considering these results, the electronic structures of Cu(II) and Fe(III), as well as their interaction with EDTA and DDTC ligands, are discussed to understand the mechanism of Fe(III)/DDTC process. By introducing a low dosage of Fe(III), the DDTC could efficiently purify Cu(II) from the Cu(II)-EDTA acid wastewater and realize the near-zero discharge of metal pollutants in metal-organic complex wastewater. It is believed that the main findings may benefit the water pollution reduction and comprehensive recycling of metal resources.

Effects of MN4-Type Coordination Structure in Metallophthalocyanine for Bio-Inspired Oxidative Desulfurization Performance

Oxidative desulfurization (ODS) is the promising new method for super deep desulfurization of fuel oil. The oxidative desulfurization performance of the metal-N₄-chelates metallophthalocyanines (MPcs) is related to the chemical properties of conjugate structures and the central metal ions. Herein, a biomimetic catalytic system composed of metallophthalocyanines (MPcR₄, M = Mn(II), Fe(II), Co(II), Ni(II), Cu(II), Zn(II); R = -H, -COOH, -NO₂, -NH₂) and molecular O₂ was performed to study the influence of MN₄-type coordination structure in metallophthalocyanines for the degradation of dibenzothiophene (DBT) in model oil containing n-octane. The results reveal that the conjugate structures and the center metal ions of metallophthalocyanines played key roles in oxidative desulfurization performance. The inductive effect of different R substituents strongly affected the electron cloud distribution of the conjugate structures and the catalytic performance. Moreover, the catalytic activity of MPcs, which is related to the d electronic configuration and ligand-field effects, does not sequentially increase with the increase in the d electron number of central metal ions.

Symmetry Breaking in Solution-Phase [Cu(tsc)2(H2O)2]2+: Emergent Asymmetry in Cu-S Distances and in Covalence

The structure of aqueous Cu(II)-bis-thiosemicarbazide, [Cu(tsc)₂]²⁺, is reported following EXAFS and MXAN analyses of the copper K-edge X-ray absorption (XAS) spectrum. The rising K-edge feature at 8987.1 eV is higher energy than those of crystalline models, implying unique electronic and structural solution states. EXAFS analysis (k = 2-13 Å^-1; 2 × Cu-N = 2.02 ± 0.01 Å; 2 × Cu-S = 2.27 ± 0.01 Å; Cu-O_ax = 2.41 ± 0.04 Å) could not resolve 5- versus 6-coordinate models. However, MXAN fits converged to an asymmetric broken symmetry 6-coordinate model with cis-disposed TSC ligands (Cu-O_ax = 2.07 and 2.54 Å; Cu-N = 1.94 Å, 1.98 Å; Cu-S = 2.20 Å, 2.41 Å). Transition dipole integral evaluation of the sulfur K-edge XAS 1s → 3p valence transition feature at 2470.7 eV yielded a Cu-S covalence of 0.66 e^-, indicating Cu^1.34+. The high Cu-S covalence and short Cu-S bond in aqueous [Cu(tsc)₂(H₂O)₂]²⁺ again contradict the need for a protein rack to explain the unique structure of the blue copper active site. MXAN models of dissolved Cu(II) complex ions have invariably featured broken centrosymmetry. The potential energy ground state for dissolved Cu(II) evidently includes the extended solvation field, providing a target for improved physical theory. A revised solvation model for aqueous Cu(II), |[Cu(H₂O)₅]·14H₂O|²⁺, is presented.

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.

EXAFS Study on the Coordination Chemistry of the Solvated Copper(II) Ion in a Series of Oxygen Donor Solvents

The structures of the solvated copper(II) ion in water and nine organic oxygen donor solvents with similar electron-pair donor ability, but with different space-demanding properties at coordination, have been studied by EXAFS. N,N′-Dimethylpropyleneurea and N,N,N′,N′-tetramethylurea are sufficiently space demanding at coordination to make the axial positions not accessible, resulting in square-planar copper(II) solvate complexes with an intense green color. The mean Cu–O bond distances in these two solvate complexes are 1.939(3) and 1.935(3) Å, respectively. The best fits of the remaining solvates, which are light blue in different hues, are obtained with a Jahn–Teller distorted-octahedral model consisting of four strongly bound solvent molecules in the equatorial positions at 1.96(2) Å and two in the axial positions but with different Cu–O_ax bond distances: ca. 2.15 and 2.32 Å. This is in agreement with observations in solid-state structures of compounds containing hexaaquacopper(II) complexes crystallizing in noncentrosymmetric space groups and all reported crystal structures containing a [Cu(H₂O)₅(O-ligand)] complex with Jahn–Teller distortion. Such a structure is in agreement with previous EPR and EXAFS studies proving the hydrated copper(II) ion to be a noncentrosymmetric complex in aqueous solution. The refinements of the EXAFS data of the solids [Cu(H₂O)₆](ClO₄)₂, [Cu(H₂O)₆](BrO₃)₂, [Cu(H₂O)₆]SiF₆, Cu(NO₃)₂·2.5H₂O, and CuSO₄·5H₂O gave Cu–O bond distances significantly different from those reported in the crystallographic studies but similar to the configuration and bond distances in the hydrated copper(II) ion in aqueous solution. This may depend on whether the orientation of the axial positions is random in one or three dimensions, giving a mean structure of the solid with symmetry higher than that of the individual complexes. This study presents the very first experimental data from the new X-ray absorption spectroscopy beamline Balder at the MAX IV synchrotron radiation facility in Lund, Sweden, as well as the utilized properties of the beamline.

The coordination chemistry of the solvated copper(II) ion has been studied in 10 solvents, including water. The copper(II) ion has a noncentrosymmetric Jahn−Teller distorted-octahedral geometry with the axial Cu−O bond distances differing by ca. 0.2 Å.

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.

Binding of Copper Ions with Octapeptide Region in Prion Protein: Simulations with Charge Transfer Model

Copper ions are important cofactors of many metalloproteins. The binding dynamics of proteins to the copper ion is important for biological functions but is less understood at the microscopic level. What are the key factors determining the recognition and the stabilization of the copper ion during the binding? Our work investigates the binding dynamics of the copper ion with a simple system (the N-terminus of PrP) using simulation methods. To precisely characterize the protein?ion interaction, we build up an effective copper?peptide force field based on quantum chemistry calculations. In our model, the effects of charge transfer, protonation/deprotonation, and induced polarization are considered. With this force field, we successfully characterize the local structures and the complex interactions of the octapeptide around the copper ion. Furthermore, using an enhanced sampling method, the binding/unbinding processes of the copper ion with the octapeptide are simulated. Free-energy landscapes are generated in consequence, and multiple binding pathways are characterized. It is observed that various native ligands contribute differently to the binding processes. Some residues are related to the capture of the ion (behaving like ?arm?s), and some others contribute to the stabilization of the coordination structure (acting like ?core?s). These different interactions induce various pathways. Besides, a nonnative binding ligand is determined, and it has essential contributions and modulations to the binding pathways. With all these results, the picture of copper?octapeptide binding is outlined. These features are believed to happen in many ion?peptide interactions, such as the cooperative stabilization between the coordinations with neighboring backbone nitrogens and an auxiliary intermediate coordination with the neighboring oxygen from the N-terminal direction. We believe that our studies are valuable to understand the complicated ion?peptide binding processes.

Cellulose Nanofiber-Graphene Oxide Biohybrids: Disclosing the Self-Assembly and Copper-Ion Adsorption Using Advanced Microscopy and ReaxFF Simulations

The self-assembly of nanocellulose and graphene oxide into highly porous biohybrid materials has inspired the design and synthesis of multifunctional membranes for removing water pollutants. The mechanisms of self-assembly, metal ion capture, and cluster formation on the biohybrids at the nano- and molecular scales are quite complex. Their elucidation requires evidence from the synergistic combination of experimental data and computational models. The AFM-based microscopy studies of (2,2,6,6-tetramethylpiperidine-1-oxylradical)-mediated oxidized cellulose nanofibers (TOCNFs), graphene oxide (GO), and their biohybrid membranes provide strong, direct evidence of self-assembly; small GO nanoparticles first attach and accumulate along a single TOCNF fiber, while the long, flexible TOCNF filaments wrap around the flat, wide GO planes, thus forming an amorphous and porous biohybrid network. The layered structure of the TOCNFs and GO membrane, derived from the self-assembly and its surface properties before and after the adsorption of Cu(II), is investigated by advanced microscopy techniques and is further clarified by the ReaxFF molecular dynamics (MD) simulations. The dynamics of the Cu(II)-ion capture by the TOCNF and GO membranes in solution and the ion cluster formation during drying are confirmed by the MD simulations. The results of this multidisciplinary investigation move the research one step forward by disclosing specific aspects of the self-assembly behavior of biospecies and suggesting effective design strategies to control the pore size and robust materials for industrial applications.

[Cu(aq)]2+ is structurally plastic and the axially elongated octahedron goes missing

High resolution (k = 18 Å^-1 or k = 17 Å^-1) copper K-edge EXAFS and MXAN (Minuit X-ray Absorption Near Edge) analyses have been used to investigate the structure of dissolved [Cu(aq)]²⁺ in 1,3-propanediol (1,3-P) or 1,5-pentanediol (1,5-P) aqueous frozen glasses. EXAFS analysis invariably found a single axially asymmetric 6-coordinate (CN6) site, with 4×O_eq = 1.97 Å, O_ax1 = 2.22 Å, and O_ax2 = 2.34 Å, plus a second-shell of 4×O_water = 3.6 Å. However, MXAN analysis revealed that [Cu(aq)]²⁺ occupies both square pyramidal (CN5) and axially asymmetric CN6 structures. The square pyramid included 4×H₂O = 1.95 Å and 1×H₂O = 2.23 Å. The CN6 sites included either a capped, near perfect, square pyramid with 5×H₂O = 1.94 ± 0.04 Å and H₂O_ax = 2.22 Å (in 1,3-P) or a split axial configuration with 4×H₂O = 1.94, H₂O_ax1 = 2.14 Å, and H₂O_ax2 = 2.28 Å (in 1,5-P). The CN6 sites also included an 8-H₂O second-shell near 3.7 Å, which was undetectable about the strictly pyramidal sites. Equatorial angles averaging 94° ± 5° indicated significant departures from tetragonal planarity. MXAN assessment of the solution structure of [Cu(aq)]²⁺ in 1,5-P prior to freezing revealed the same structures as previously found in aqueous 1M HClO₄, which have become axially compressed in the frozen glasses. [Cu(aq)]²⁺ in liquid and frozen solutions is dominated by a 5-coordinate square pyramid, but with split axial CN6 appearing in the frozen glasses. Among these phases, the Cu-O axial distances vary across 1 Å, and the equatorial angles depart significantly from the square plane. Although all these structures remove the d_x²-y² , d_z² degeneracy, no structure can be described as a Jahn-Teller (JT) axially elongated octahedron. The JT-octahedral description for dissolved [Cu(aq)]²⁺ should thus be abandoned in favor of square pyramidal [Cu(H₂O)₅]²⁺. The revised ligand environments have bearing on questions of the Cu(i)/Cu(ii) self-exchange rate and on the mechanism for ligand exchange with bulk water. The plasticity of dissolved Cu(ii) complex ions falsifies the foundational assumption of the rack-induced bonding theory of blue copper proteins and obviates any need for a thermodynamically implausible protein constraint.

Molecular Facts on the Structure and Dynamics of Electrolyte Species in Cu-Cl Cycle for Hydrogen Generation: An Insight from Molecular Dynamic Simulations

The Cu complex, which is the key chemical species in well-known Cu-Cl hybrid thermochemical cycles and also in numerous metal hydrometallurgical and sedimentary deposit processes, displays a wide variety of structural and dynamical characteristics that are further complicated by the presence of multiple oxidation states of Cu ions with different coordination chemistries, therefore they are difficult to explore from experiments alone. In this article, an attempt has been made to understand the coordination behavior of the Cu complex using MD simulations. The study provides compelling evidence of the experimentally observed multiple stoichiometries of Cu ions, i.e., 1:6:0, 1:5:1, and 1:4:2 for Cu⁺:H₂O:Cl^- and 1:6:0 for Cu²⁺:H₂O:Cl^-. The presence of the anionic Cu complex, [Cu⁺Cl₂]^-·2H₂O, [Cu⁺Cl₂]^-·3H₂O, [Cu²⁺Cl₃]^-·H₂O, and [Cu²⁺Cl₃]^-·2H₂O, was captured in the presence of excess chloride ions. Furthermore, the probability distribution profiles have been estimated to determine the most possible complex in the considered systems. The results establish structural and dynamical reformation of the Cu complex with change in the salt concentration or variation in the solvent medium in which they are dissolved. Moreover, the structure and kinetics of the Cu ions in the Cu-Cl electrolyzer have been explored over a large range of the electric field by extending the simulated systems for varied strengths of the electric fields. It has been observed that with an increase in the strength of the electric field, the water molecules lose their coordination strength with central Cu ions, which, on the other hand, results in a significant change in the structure of the captured complex. The diffusion dynamics of the ions is altered while applying the electric field, which is furthermore modified while increasing the strength of electric field beyond a critical limit. In fact, the diffusion mechanism of the ions was seen to be transformed from Brownian-like to linear motion and then to hopping diffusion with the increasing strength of the electric field. To the best of our knowledge, this is the first time when the multiple oxidation states of the Cu ion are explored using MD simulations, and the coexisting pictures of the multiple coordinations and the solvent effects have been clearly revealed. Also to date, the present article is the first one to report the insights of the structure and the dynamics of the ions in the Cu-Cl electrolyzer over a wide range of the electric field. The present studies will be very helpful in understanding the mechanism involved in numerous metal hydrometallurgical and sedimentary deposit processes and to comprehend the analogies involved in the electrode reactions of the Cu-Cl cycle for hydrogen generation.

Copper Oxidation/Reduction in Water and Protein: Studies with DFTB3/MM and VALBOND Molecular Dynamics Simulations

We apply two recently developed computational methods, DFTB3 and VALBOND, to study copper oxidation/reduction processes in solution and protein. The properties of interest include the coordination structure of copper in different oxidation states in water or in a protein (plastocyanin) active site, the reduction potential of the copper ion in different environments, and the environmental response to copper oxidation. The DFTB3/MM and VALBOND simulation results are compared to DFT/MM simulations and experimental results whenever possible. For a copper ion in aqueous solution, DFTB3/MM results are generally close to B3LYP/MM with a medium basis, including both solvation structure and reduction potential for Cu(II); for Cu(I), however, DFTB3/MM finds a two-water coordination, similar to previous Born-Oppenheimer molecular dynamics simulations using BLYP and HSE, whereas B3LYP/MM leads to a tetrahedron coordination. For a tetraammonia copper complex in aqueous solution, VALBOND and DFTB3/MM are consistent in terms of both structural and dynamical properties of solvent near copper for both oxidation states. For copper reduction in plastocyanin, DFTB3/MM simulations capture the key properties of the active site, and the computed reduction potential and reorganization energy are in fair agreement with experiment, especially when the periodic boundary condition is used. Overall, the study supports the value of VALBOND and DFTB3(/MM) for the analysis of fundamental copper redox chemistry in water and protein, and the results also help highlight areas where further improvements in these methods are desirable.

A high-resolution XAS study of aqueous Cu(II) in liquid and frozen solutions: pyramidal, polymorphic, and non-centrosymmetric

High-resolution EXAFS (k = 18 Å(-1)) and MXAN XAS analyses show that axially elongated square pyramidal [Cu(H2O)5](2+) dominates the structure of Cu(II) in aqueous solution, rather than 6-coordinate JT-octahedral [Cu(H2O)6](2+). Freezing produced a shoulder at 8989.6 eV on the rising XAS edge and an altered EXAFS spectrum, while 1s → 3d transitions remained invariant in energy position and intensity. Core square pyramidal [Cu(H2O)5](2+) also dominates frozen solution. Solvation shells were found at ∼3.6 Å (EXAFS) or ∼3.8 Å (MXAN) in both liquid and frozen phases. However, MXAN analysis revealed that about half the time in liquid solution, [Cu(H2O)5](2+) associates with an axially non-bonding 2.9 Å water molecule. This distant water apparently organizes the solvation shell. When the 2.9 Å water molecule is absent, the second shell is undetectable to MXAN. The two structural arrangements may represent energetic minima of fluxional dissolved aqueous [Cu(H2O)5](2+). The 2.9 Å trans-axial water resolves an apparent conflict of the [Cu(H2O)5](2+) core model with a dissociational exchange mechanism. In frozen solution, [Cu(H2O)5](2+) is associated with either a 3.0 Å axial non-bonded water molecule or an axial ClO4(-) at 3.2 Å. Both structures are again of approximately equal presence. When the axial ClO4(-) is present, Cu(II) is ∼0.5 Å above the mean O4 plane. This study establishes [Cu(H2O)5](2+) as the dominant core structure for Cu(II) in water solution, and is the first to both empirically resolve multiple extended solution structures for fluxional [Cu(H2O)5](2+) and to provide direct evidence for second shell dynamics.

Characterization of structure and dynamics of an aqueous scandium(III) ion by an extended ab initio QM/MM molecular dynamics simulation

Hydration structure and dynamics of an aqueous Sc(III) solution were characterized by means of an extended ab initio quantum mechanical/molecular dynamical (QM/MM) molecular dynamics simulation at Hartree-Fock level. A monocapped trigonal prismatic structure composed of seven water molecules surrounding scandium(III) ion was proposed by the QM/MM simulation including the quantum mechanical effects for the first and second hydration shells. The mean Sc(III)-O bond length of 2.14 Å was identified for six prism water molecules with one capping water located at around 2.26 Å, reproducing well the X-ray diffraction data. The Sc(III)-O stretching frequency of 432 cm(-1) corresponding to a force constant of 130 N m(-1), evaluated from the enlarged QM/MM simulation, is in good agreement with the experimentally determined value of 430 cm(-1) (128 N m(-1)). Various water exchange processes in the second hydration shell of the hydrated Sc(III) ion predict a mean ligand residence time of 7.3 ps.

Mononuclear-dinuclear equilibrium of grafted copper complexes confined in the nanochannels of MCM-41 silica

Following the structural concept of copper-containing proteins in which dinuclear copper centers are connected by hydroxide bridging ligands, a bidentate copper(II) complex has been incorporated into nano-confined MCM-41 silica by a multistep sequential grafting technique. Characterization by a combination of EPR spectroscopy, X-ray photoelectron spectroscopy (XPS), UV/Vis spectroscopy, IR spectroscopy , and solid-state (13)C and (29)Si cross-polarization magic-angle spinning (CP-MAS) NMR suggests that dinuclear Cu complexes are bridged by hydroxide and other counterions (chloride or perchlorate ions), similar to the situation for EPR-undetectable [Cu(II)···Cu(II)] dimer analogues in biological systems. More importantly, a dynamic mononuclear-dinuclear equilibrium between different coordination modes of copper is observed, which strongly depends on the nature of the counterions (Cl(-) or ClO(4)(-)) in the copper precursor and the pore size of the silica matrix (the so-called confinement effect). A proton-transfer mechanism within the hydrogen-bonding network is suggested to explain the dynamic nature of the dinuclear copper complex supported on the MCM-41 silica.

Development and validation of a ReaxFF reactive force field for Cu cation/water interactions and copper metal/metal oxide/metal hydroxide condensed phases

To enable large-scale reactive dynamic simulations of copper oxide/water and copper ion/water interactions we have extended the ReaxFF reactive force field framework to Cu/O/H interactions. To this end, we employed a multistage force field development strategy, where the initial training set (containing metal/metal oxide/metal hydroxide condensed phase data and [Cu(H(2)O)(n)](2+) cluster structures and energies) is augmented by single-point quantum mechanices (QM) energies from [Cu(H(2)O)(n)](2+) clusters abstracted from a ReaxFF molecular dynamics simulation. This provides a convenient strategy to both enrich the training set and to validate the final force field. To further validate the force field description we performed molecular dynamics simulations on Cu(2+)/water systems. We found good agreement between our results and earlier experimental and QM-based molecular dynamics work for the average Cu/water coordination, Jahn-Teller distortion, and inversion in [Cu(H(2)O)(6)](2+) clusters and first- and second-shell O-Cu-O angular distributions, indicating that this force field gives a satisfactory description of the Cu-cation/water interactions. We believe that this force field provides a computationally convenient method for studying the solution and surface chemistry of metal cations and metal oxides and, as such, has applications for studying protein/metal cation complexes, pH-dependent crystal growth/dissolution, and surface catalysis.

Temperature effects on the structure and dynamics of the Jahn-Teller distorted Cr(2+) ion in aqueous solution: a hybrid QM/MM molecular dynamics simulation

The influences of an elevated temperature on the structure and dynamics of the Jahn-Teller distorted [Cr(H(2)O)(6)](2+) complex have been studied using an ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation, showing that the increased temperature affects the lifetime distortions of the hydrated Cr(2+) ion by decreasing the inversion time to 0.5-2 ps.