The hydration of Cu2+: Can the Jahn-Teller effect be detected in liquid solution?

Hydration numbers of biologically relevant divalent metal cations from ab initio molecular dynamics and continuum solvation methods

Hydration and, in particular, the coordination number of a metal ion is of paramount importance as it defines many of its (bio)physicochemical properties. It is not only essential for understanding its behavior in aqueous solutions but also determines the metal ion reference state and its binding energy to (bio)molecules. In this paper, for divalent metal cations Ca2+, Cd2+, Cu2+, Fe2+, Hg2+, Mg2+, Ni2+, Pb2+, and Zn2+, we compare two approaches for predicting hydration numbers: (1) a mixed explicit/continuum DFT-D3//COSMO-RS solvation model and (2) density functional theory based ab initio molecular dynamics. The former approach is employed to calculate the Gibbs free energy change for the sequential hydration reactions, starting from [M(H2O)2]2+ aqua complexes to [M(H2O)9]2+, allowing explicit water molecules to bind in the first or second coordination sphere and determining the most stable [M(H2O)n]2+ structure. In the latter approach, the hydration number is obtained by integrating the ion-water radial distribution function. With a couple of exceptions, the metal ion hydration numbers predicted by the two approaches are in mutual agreement, as well as in agreement with the experimental data.

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.

A Combined Wave Function and Density Functional Approach for K-Edge X-ray Absorption Near-Edge Spectroscopy: A Case Study of Hydrated First-Row Transition Metal Ions

The prediction of X-ray absorption spectra (XAS) of transition metal complexes has important and broad application areas in chemistry and biology. In this letter, we have investigated the predictive ability of multiconfiguration pair-density functional theory (MC-PDFT) for X-ray absorption spectra by calculating the metal K pre-edge features of aquated 3d transition metal ions in common oxidation states. MC-PDFT results were compared with experimentally measured spectra as well as analyzed against results from restricted active-space second-order perturbation theory (RASPT2) and time-dependent density functional theory (TDDFT). As expected, TDDFT performs well for excited states that can be accurately represented by singly excited configurations but fails for excited states where higher order excitations become important. On the other hand, both RASPT2 and MC-PDFT provide quantitatively accurate results for all excited states irrespective of their character. While core-level spectroscopy with RASPT2 is accurate, it is computationally expensive. Our results show that MC-PDFT performs equally well with significantly lower computational cost and is an encouraging alternate approach for X-ray spectroscopies.

A Comparison of Bonded and Nonbonded Zinc(II) Force Fields with NMR Data

Classical molecular dynamics (MD) simulations are widely used to inspect the behavior of zinc(II)-proteins at the atomic level, hence the need to properly model the zinc(II) ion and the interaction with its ligands. Different approaches have been developed to represent zinc(II) sites, with the bonded and nonbonded models being the most used. In the present work, we tested the well-known zinc AMBER force field (ZAFF) and a recently developed nonbonded force field (NBFF) to assess how accurately they reproduce the dynamic behavior of zinc(II)-proteins. For this, we selected as benchmark six zinc-fingers. This superfamily is extremely heterogenous in terms of architecture, binding mode, function, and reactivity. From repeated MD simulations, we computed the order parameter (S2) of all backbone N-H bond vectors in each system. These data were superimposed to heteronuclear Overhauser effect measurements taken by NMR spectroscopy. This provides a quantitative estimate of the accuracy of the FFs in reproducing protein dynamics, leveraging the information about the protein backbone mobility contained in the NMR data. The correlation between the MD-computed S2 and the experimental data indicated that both tested FFs reproduce well the dynamic behavior of zinc(II)-proteins, with comparable accuracy. Thus, along with ZAFF, NBFF represents a useful tool to simulate metalloproteins with the advantage of being extensible to diverse systems such as those bearing dinuclear metal sites.

Near-Quantitative Predictions of the First-Shell Coordination Structure of Hydrated First-Row Transition Metal Ions Using K-Edge X-ray Absorption Near-Edge Spectroscopy

The solvation structure of transition metal ions is important for applications in geochemistry, biochemistry, energy storage, and environmental chemistry. We study the X-ray absorption pre-edge and near-edge spectra at the K-edge of a nearly complete series of hydrated first-row transition metal ions with d orbital occupancy from d² to d¹⁰. We optimize all of the structures at the density functional theory (DFT) level with explicit solvation and then compute the pre-edge X-ray absorption spectra with time-dependent DFT (TDDFT) and restricted active space second-order perturbation theory (RASPT2). TDDFT provides accurate results for spectra that are dominated by single excitations, while RASPT2 correctly distinguishes between singly and doubly excited states with quantitative accuracy compared with experiment. We analyze the pre-edge features for each metal ion to reveal the impact of the variations in d orbital occupancy on the first-shell coordination environment. We also report the lowest-energy ligand field d-d transitions using complete active space second-order perturbation theory.

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.

Resolving Configurational Disorder for Impurities in a Low-Entropy Phase

Hematite (α-Fe₂O₃) exerts a strong control over the transport of minor but critical metals in the environment and is used in multiple industrial applications; the photocatalysis community has explored the properties of hematite nanoparticles over a wide range of transition metal dopants. Nonetheless, simplistic assumptions are used to rationalize the local coordination environment of impurities in hematite. Here, we use ab initio molecular dynamics (AIMD)-guided structural analysis to model the extended X-ray absorption fine structure (EXAFS) of Cu²⁺- and Zn²⁺-doped hematite nanoparticles. Specific defect-impurity associations were identified, and the local coordination environments of Cu and Zn both displayed considerable configurational disorder that, in aggregate, approached Jahn-Teller-like distortion for Cu but, in contrast, maintained hematite-like symmetry for Zn. This study highlights the role of defects in accommodating impurities in a nominally low-entropy phase and the limits to traditional shell-by-shell fitting of EXAFS for dopants/impurities in unprecedented bonding environments.

A spectroscopic study of solvent effects on the formation of Cu(ii)-chloride complexes in aqueous solution

A combination of electronic (UV-vis) and X-ray absorption (EXAFS, XANES) spectroscopies has been used to investigate the formation of copper(ii)/chloride complexes in concentrated aqueous solutions. It is established that lowering the water activity by the addition of Mg(ClO₄)₂ at a constant Cl^-/Cu(ii) ratio results in the replacement of water molecules by Cl^- ions in the primary coordination shell of Cu(ii). This behavior closely parallels the effect of increasing the Cl^-/Cu(ii) ratio and demonstrates that full understanding of the stoichiometry and structures of the complexes formed in concentrated metal-ion chloride solutions requires explicit consideration of the role of the solvent.

Data science assisted investigation of catalytically active copper hydrate in zeolites for direct oxidation of methane to methanol using H2O2

Dozens of Cu zeolites with MOR, FAU, BEA, FER, CHA and MFI frameworks are tested for direct oxidation of CH₄ to CH₃OH using H₂O₂ as oxidant. To investigate the active structures of the Cu zeolites, 15 structural variables, which describe the features of the zeolite framework and reflect the composition, the surface area and the local structure of the Cu zeolite active site, are collected from the Database of Zeolite Structures of the International Zeolite Association (IZA). Also analytical studies based on inductively coupled plasma-optical emission spectrometry (ICP-OES), X-ray fluorescence (XRF), N₂ adsorption specific surface area measurement and X-ray absorption fine structure (XAFS) spectral measurement are performed. The relationships between catalytic activity and the structural variables are subsequently revealed by data science techniques, specifically, classification using unsupervised and supervised machine learning and data visualization using pairwise correlation. Based on the unveiled relationships and a detailed analysis of the XAFS spectra, the local structures of the Cu zeolites with high activity are proposed.

Development of 64Cu-loaded Perfluoropentane Nanodroplet: A Potential Tumor Theragnostic Nano-carrier and Dual-Modality PET-Ultrasound Imaging Agents

The purpose of this study was to develop and preliminarily evaluate phospholipid-shelled nanodroplets (NDs) encapsulating perfluoropentane (PFP) and radioactive ⁶⁴Cu as a hybrid positron emission tomography (PET)-ultrasound (US) probe. PFP NDs were fabricated by mixing liquid-phase PFP with a phospholipid solution. The ⁶⁴Cu was encapsulated into the NDs in a size-controlled manner by exploiting the hydrophobicity of ⁶⁴Cu-diacetyl-bis(N⁴-methylthiosemicarbazone) (⁶⁴Cu-ATSM) using a vial mixer and an extruder. The fabricated ⁶⁴Cu-loaded PFP NDs (⁶⁴Cu-PFP NDs) were evaluated using in vitro/in vivo PET-computed tomography (PET-CT), US imaging and transmission electron microscopy. In the in vitro PET images, the ⁶⁴Cu-PFP NDs were observed as a hot spot in the lower section of the test tube. In the acquired US images, the mean region of interest brightness values of ⁶⁴Cu-PFP NDs were revealed by their strong echo image. In a tumor-bearing mouse animal model, tumor uptake of the ⁶⁴Cu-PFP NDs was low, that is, approximately 65%, compared with that of only free ⁶⁴Cu, as determined by PET-delayed imaging analysis. The dual-function concept of the NDs is expected to contribute to the prognosis and effectiveness of therapy by fusing the science and technology of nuclear medicine and US.

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 Å.

A model for the effect of ion pairing on an outer sphere electron transfer

Ion pairing can strongly affect the rates of electron transfer reactions. To explain this effect, we propose a model Hamiltonian that describes the interactions between the pairing ion and the reactant, solvent and inner sphere reorganization, and bond breaking. Explicit expressions for the energies of the initial and final states, and for the energy of activation are derived in the weak adiabatic limit. The model is applied to the reduction of Cu(ii) in the presence of chloride ions. For this purpose, the pertinent system parameters are obtained from density functional theory. Our model explains why the chloride ion enhances the rate of the first electron transfer in copper deposition.

[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.

Composition-driven Cu-speciation and reducibility in Cu-CHA zeolite catalysts: a multivariate XAS/FTIR approach to complexity

The small pore Cu-CHA zeolite is attracting increasing attention as a versatile platform to design novel single-site catalysts for deNO _x applications and for the direct conversion of methane to methanol. Understanding at the atomic scale how the catalyst composition influences the Cu-species formed during thermal activation is a key step to unveil the relevant composition-activity relationships. Herein, we explore by in situ XAS the impact of Cu-CHA catalyst composition on temperature-dependent Cu-speciation and reducibility. Advanced multivariate analysis of in situ XANES in combination with DFT-assisted simulation of XANES spectra and multi-component EXAFS fits as well as in situ FTIR spectroscopy of adsorbed N₂ allow us to obtain unprecedented quantitative structural information on the complex dynamics during the speciation of Cu-sites inside the framework of the CHA zeolite.

Structure and Stability of a Copper(II) Lactate Complex in Alkaline Solution: a Case Study by Energy-Dispersive X-ray Absorption Spectroscopy

Energy-dispersive X-ray absorption spectroscopy was applied, aimed at solving the problem of the structure and stability of a copper(II) lactate complex in alkaline solution, used as a precursor for the electrodeposition of Cu₂O. The application of multiple scattering calculations to the simulation of the X-ray absorption near-edge structure part of the spectra allowed an accurate resolution of the structure: the copper(II) cation is surrounded by four lactate ions in a distorted tetrahedral environment, with the lactate anions acting as monodentate ligands. This results in an atomic arrangement where copper is surrounded by four oxygen atoms located at quite a short distance (ca. 1.87 Å) and four oxygen atoms located quite far apart (ca. 3.1-3.2 Å). The complex was finally found to be stable in a wide range of applied potentials.

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

Identifying Coordination Geometries of Metal Aquaions in Water: Application to the Case of Lanthanoid and Actinoid Hydrates

The angular distribution function (ADF) associated with the oxygen-metal ion-oxygen angle (OMO) of several trivalent lanthanoid and actinoid aquaions has been used to identify the most probable coordination geometry of these aquaions in aqueous solutions. The ADFs extracted from the molecular dynamics trajectories have been compared with continuous distribution functions corresponding to the geometry of a reference polyhedron pattern. The procedure incorporates specific quantum-mechanical information on the aquaion under study. The new method is applied to the analysis of four M(H₂O)_n³⁺ aquaions in water, M = Lu and Cf for n = 8, and M = La and Ac for n = 9. For those that are 8-coordinated, the square antiprism (SA) coordination geometry is preferred. For the 9-fold coordination, the simulation ADFs are more similar to the continuous ADF derived from a Gyro-elongated-SA rather than to the usually proposed trigonal tricapped prism. Advantages of these continuous distributions with respect to the usually employed discrete distributions are emphasized as well as further applications are suggested.

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.

Origins of optical absorption characteristics of Cu2+complexes in aqueous solutions

Spectrometry is combined with first-principles simulations of solvated Cu²⁺complexes to elucidate the roles of ligand coordination, d-orbital fine structure, complex geometry, and solvation environment in determining optical absorption profiles.

Revisiting the nature of Cu sites in the activated Cu-SSZ-13 catalyst for SCR reaction

Cu-SSZ-13 is a highly active NH3-SCR catalyst for the abatement of harmful nitrogen oxides (NO x , x = 1, 2) from the exhausts of lean-burn engines. The study of Cu-speciation occurring upon thermal dehydration is a key step for the understanding of the enhanced catalytic properties of this material and for identifying the SCR active sites and their redox capability. Herein, we combined FTIR, X-ray absorption (XAS) and emission (XES) spectroscopies with DFT computational analysis to elucidate the nature and location of the most abundant Cu sites in the activated catalyst. Different Cu species have been found to be dominant as a function of the dehydration temperature and conditions. Data analysis revealed that the dehydration process of Cu cations is essentially completed at 250 °C, with the formation of dehydrated [CuOH]+ species hosted in close proximity to 1-Al sites in both d6r and 8r units of the SSZ-13 matrix. These species persist at higher temperatures only if a certain amount of O2 is present in the gas feed, while under inert conditions they undergo virtually total "self-reduction" as a consequence of an OH extra-ligand loss, resulting in bi-coordinated bare Cu+ cations. Synchrotron characterization supported by computational analysis allowed an unprecedented quantitative refinement of the local environment and structural parameters of these Cu(ii) and Cu(i) species.

An Angular Overlap Model for Cu(II) Ion in the AMOEBA Polarizable Force Field

An extensible polarizable force field for transition metal ion was developed based on AMOEBA and the angular overlap model (AOM) with consistent treatment of electrostatics for all atoms. Parameters were obtained by fitting molecular mechanics (MM) energies to various ab initio gas-phase calculations. The results of parameterization were presented for copper (II) ion ligated to water and model fragments of amino acid residues involved in the copper binding sites of type 1 copper proteins. Molecular dynamics (MD) simulations were performed on aqueous copper (II) ion at various temperatures, as well as plastocyanin (1AG6) and azurin (1DYZ). Results demonstrated that the AMOEBA-AOM significantly improves the accuracy of classical MM in a number of test cases when compared to ab initio calculations. The Jahn-Teller distortion for hexa-aqua copper (II) complex was handled automatically without specifically designating axial and in-plane ligands. Analyses of MD trajectories resulted in a 6-coordination first solvation shell for aqueous copper (II) ion and a 1.8ns average residence time of water molecules. The ensemble average geometries of 1AG6 and 1DYZ copper binding sites were in general agreement with X-ray and previous computational studies.

The structure of water in solutions containing di- and trivalent cations by empirical potential structure refinement

Empirical potential structure refinement (EPSR) has been used to build experimentally consistent models of a range of electrolyte solutions containing di- and trivalent cations: Cu(ClO4)2 at concentrations of 0.5 and 2.0 m, and Cr(NO3)3, YCl3 and LaCl3 at a concentration of 1.0 m. The resulting models are used to investigate the perturbation of these electrolytes on the pair distribution and triplet angle correlations between solvent water molecules, compared with those found in the pure solvent. The results elucidate the differences that derive from the reflected range of highly structured local cation environments and provide a complementary viewpoint on the hydration shell geometries.