Anion–arene adducts: C–H hydrogen bonding, anion– interaction, and carbon bonding motifs

This article summarizes experimental and theoretical evidence for the existence of four distinct binding modes for complexes of anions with charge-neutral arenes. These include C–H hydrogen bonding and three motifs involving the arene– system—the noncovalent anion– interaction, weakly covalent interaction, and strongly covalent interaction

Computational Study of Copper(II) Complexation and Hydrolysis in Aqueous Solutions Using Mixed Cluster/Continuum Models

We use density functional theory (B3LYP) and the COSMO continuum solvent model to characterize the structure and stability of the hydrated Cu(II) complexes [Cu(MeNH_2)(H_2O)_(n−1)]^(2+) and [Cu(OH)_x(H_2O)_(n−x)]^(2−x) (x = 1−3) as a function of metal coordination number (4−6) and cluster size (n = 4−8, 18). The small clusters with n = 4−8 are found to be the most stable in the nearly square-planar four-coordinate configuration, except for [Cu(OH)_3(H_2O)]^−, which is three-coordinate. In the presence of the two full hydration shells (n = 18), however, the five-coordinate square-pyramidal geometry is the most favorable for Cu(MeNH_2)^(2+) (5, 6) and Cu(OH)^+ (5, 4, 6), and the four-coordinate geometry is the most stable for Cu(OH)_2 (4, 5) and Cu(OH)_3^− (4). (Other possible coordination numbers for these complexes in the aqueous phase are shown in parentheses.) A small energetic difference between these structures (0.23−2.65 kcal/mol) suggests that complexes with different coordination numbers may coexist in solution. Using two full hydration shells around the Cu^(2+) ion (18 ligands) gives Gibbs free energies of aqueous reactions that are in excellent agreement with experiment. The mean unsigned error is 0.7 kcal/mol for the three consecutive hydrolysis steps of Cu^(2+) and the complexation of Cu^(2+) with methylamine. Conversely, calculations for the complexes with only one coordination shell (four equatorial ligands) lead to a mean unsigned error that is >6.0 kcal/mol. Thus, the explicit treatment of the first and the second shells is critical for the accurate prediction of structural and thermodynamic properties of Cu(II) species in aqueous solution

PAMAM Dendrimers Undergo pH Responsive Conformational Changes without Swelling

Atomistic molecular dynamics (MD) simulations of a G4-NH2 PAMAM dendrimer were carried out in aqueous solution using explicit water molecules and counterions (with the Dreiding III force field optimized using quantum mechanics). Our simulations predict that the radius of gyration (R_g) of the dendrimer changes little with pH from 21.1 Å at pH 10 (uncharged PAMAM) to 22.1 Å at pH 5 (charged with 126 protons), which agrees quantitatively with recent small angle neutron scattering (SANS) experiments (from 21.4 Å at pH 10 to 21.5 Å at pH 5). Even so we predict a dramatic change in the conformation. The ion pairing in the low pH form leads to a locally compact dense shell with an internal surface area only 37% of the high pH form with a dense core. This transformation from “dense core” at high pH to “dense shell” at low pH could facilitate the encapsulation and release of guest molecules (e.g., drugs) using pH as the trigger, making dendrimers a unique drug delivery vehicle

Evaluation of B3LYP, X3LYP, and M06-class density functionals for predicting the binding energies of neutral, protonated, and deprotonated water clusters

In this paper we assess the accuracy of the B3LYP, X3LYP, and newly developed M06-L, M06-2X, and M06 functionals to predict the binding energies of neutral and charged water clusters including (H_2O)_n, n = 2−8, 20), H_3O+(H_2O_)n, n = 1−6, and OH−(H_2O)_n, n = 1−6. We also compare the predicted energies of two ion hydration and neutralization reactions on the basis of the calculated binding energies. In all cases, we use as benchmarks calculated binding energies of water clusters extrapolated to the complete basis set limit of the second-order Møller−Plesset perturbation theory with the effects of higher order correlation estimated at the coupled-cluster theory with single, double, and perturbative triple excitations in the aug-cc-pVDZ basis set. We rank the accuracy of the functionals on the basis of the mean unsigned error (MUE) between calculated benchmark and density functional theory energies. The corresponding MUE (kcal/mol) for each functional is listed in parentheses. We find that M06-L (0.73) and M06 (0.84) give the most accurate binding energies using very extended basis sets such as aug-cc-pV5Z. For more affordable basis sets, the best methods for predicting the binding energies of water clusters are M06-L/aug-cc-pVTZ (1.24), B3LYP/6-311++G(2d,2p) (1.29), and M06/aug-cc-PVTZ (1.33). M06-L/aug-cc-pVTZ also gives more accurate energies for the neutralization reactions (1.38), whereas B3LYP/6-311++G(2d,2p) gives more accurate energies for the ion hydration reactions (1.69)

pK_a Calculations of Aliphatic Amines, Diamines, and Aminoamides via Density Functional Theory with a Poisson−Boltzmann Continuum Solvent Model

In order to make reliable predictions of the acid−base properties of macroligands with a large number of ionizable sites such as dendrimers, one needs to develop and validate computational methods that accurately estimate the acidity constants (pK_a) of their chemical building blocks. In this article, we couple density functional theory (B3LYP) with a Poisson−Boltzmann continuum solvent model to calculate the aqueous pK_a of aliphatic amines, diamines, and aminoamides, which are building blocks for several classes of dendrimers. No empirical correction terms were employed in the calculations except for the free energy of solvation of the proton (H^+) adjusted to give the best match with experimental data. The use of solution-phase optimized geometries gives calculated pK_a values in excellent agreement with experimental measurements. The mean absolute error is 0.5 pK_a unit compared to experimental measurements. We find that geometry optimization in solution is essential for making accurate pK_a predictions for systems possessing intramolecular hydrogen bonds

Addition of H_2O and O_2 to Acetone and Dimethylsulfoxide Ligated Uranyl(V) Dioxocations

Gas-phase complexes of the formula [UO_2(lig)]^+ (lig = acetone (aco) or dimethylsulfoxide (dmso)) were generated by electrospray ionization (ESI) and studied by tandem ion-trap mass spectrometry to determine the general effect of ligand charge donation on the reactivity of UO_2^+ with respect to water and dioxygen. The original hypothesis that addition of O_2 is enhanced by strong σ-donor ligands bound to UO_2^+ is supported by results from competitive collision-induced dissociation (CID) experiments, which show near exclusive loss of H_2O from [UO_2(dmso)(H_2O)(O_2)]^+, whereas both H_2O and O_2 are eliminated from the corresponding [UO_2(aco)(H_2O)(O_2)]^+ species. Ligand-addition reaction rates were investigated by monitoring precursor and product ion intensities as a function of ion storage time in the ion-trap mass spectrometer: these experiments suggest that the association of dioxygen to the UO_2^+ complex is enhanced when the more basic dmso ligand was coordinated to the metal complex. Conversely, addition of H_2O is favored for the analogous complex ion that contains an aco ligand. Experimental rate measurements are supported by density function theory calculations of relative energies, which show stronger bonds between UO_2^+ and O_2 when dmso is the coordinating ligand, whereas bonds to H_2O are stronger for the aco complex

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_2O)_n]^(2+) cluster structures and energies) is augmented by single-point quantum mechanices (QM) energies from [Cu(H_2O)_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_2O)_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

ab initio informed inelastic neutron scattering for time-resolved local dynamics in molten MgCl2

Ion dynamics that drive the transport and thermophysical properties of molten salts are poorly understood due to challenges in precisely quantifying the spatial and temporal fluctuations of specific ions in highly disordered systems. While the Van Hove correlation function (VHF) obtained from inelastic neutron scattering (INS) probes these dynamics directly, its interpretation is limited by the inherent species-averaging of experiments, which obscures analysis of key ion transport and solvation mechanisms. Here, ab initio molecular dynamics (AIMD) is used to model the VHF, unravel its partial contributions, and elucidate its underlying ionic transport mechanisms. Slow decorrelation is revealed for oppositely charged ions (Mg2+ and Cl-) caused by ion exchange across the solvation shell between adjoining ionocovalent complexes. Furthermore, transport coefficients are accurately recovered and connections between macroscopic properties and ion dynamics are revealed. This study demonstrates the potential of ab initio-informed VHF to resolve long-standing challenges in uncovering relationships between picosecond-scale ion dynamics, mechanisms, and emergent physical properties of molten salts

Experimental and Computational Analysis of the Solvent-Dependent O2/Li+-O2− Redox Couple: Standard Potentials, Coupling Strength, and Implications for Lithium-Oxygen Batteries

Understanding and controlling the kinetics of O[subscript 2] reduction in the presence of Li+-containing aprotic solvents, to either Li+-O[subscript 2]− by one-electron reduction or Li[subscript 2]O[subscript 2] by two-electron reduction, is instrumental to enhance the discharge voltage and capacity of aprotic Li-O[subscript 2] batteries. Standard potentials of O[subscript 2]/Li+-O[subscript 2]− and O[subscript 2]/O[subscript 2]− were experimentally measured and computed using a mixed cluster-continuum model of ion solvation. Increasing combined solvation of Li+ and O[subscript 2]− was found to lower the coupling of Li+-O[subscript 2]− and the difference between O[subscript 2]/Li+-O[subscript 2]− and O[subscript 2]/O[subscript 2]− potentials. The solvation energy of Li+ trended with donor number (DN), and varied greater than that of O[subscript 2]− ions, which correlated with acceptor number (AN), explaining a previously reported correlation between Li+-O[subscript 2]− solubility and DN. These results highlight the importance of the interplay between ion–solvent and ion–ion interactions for manipulating the energetics of intermediate species produced in aprotic metal–oxygen batteries.National Science Foundation (U.S.) (NSF award no. ECS-0335765)China Clean Energy Research Center-Clean Vehicles Consortium (CERC-CVC)United States. Dept. of Energy (Award number DEPI0000012)National Science Foundation (U.S.) (Award number DMR-0819762)Robert Bosch GmbH (Bosch Energy Research Network (BERN) Grant)Skolkovo Institute of Science and Technology (Skoltech-MIT Center for Electochemical Energy Storage