A combined QM/MM molecular dynamics simulations study of nitrate anion (NO3-) in aqueous solution

The structural and dynamical properties of NO3- in dilute aqueous solution have been investigated by means of two combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations, namely HF/MM and B3LYP/MM, in which the ion and its surrounding water molecules were treated at HF and B3LYP levels of accuracy, respectively, using the DZV+ basis set. On the basis of both HF and B3LYP methods, a well-defined first hydration shell of NO3- is obtainable, but the shell is quite flexible and the hydrogen-bond interactions between NO3- and water are rather weak. With respect to the detailed analysis of the geometrical arrangement and vibrations of NO3-, the experimentally observed solvent-induced symmetry breaking of the ion is well reflected. In addition, the dynamical information, i.e., the bond distortions and shifts in the corresponding bending and stretching frequencies as well as the mean residence time of water molecules surrounding the NO3- ion, clearly indicates the "structure-breaking" ability of this ion in aqueous solution. From a methodical point of view it seems that both the HF and B3LYP methods are not too different in describing this hydrated ion by means of a QM/MM simulation. However, the detailed analysis of the dynamics properties indicates a better suitability of the HF method compared to the B3LYP-DFT approach.     

The electronic structure of Keggin anion (PW12O40)3− and catalytic properties

The electronic structure of the Keggin anion (PW₁₂O₄₀)^3? has been calculated by the DV–Xα method. The nature of the chemical bonds of the Keggin anion and the relationship between electronic structure and catalytic behavior are discussed. At the same time, the difference in electronic structure and the relationship in catalytic behavior between (PW₁₂O₄₀)^3? and (PMo₁₂O₄₀)^3? are compared. © 1992 John Wiley & Sons, Inc.     

Rational addition of capping groups to the phosphomolybdate Keggin anion [PMo12O40](3-) by mild, non-aqueous reductive aggregation

Controlled reductive assembly of capped Keggin anions [PMo(12)O(40)(ML(m))(n)](3-) has been achieved by reduction of [PMo(12)O(40)](3-) with sodium-mercury amalgam in the presence of metal halides, as exemplified by the rational syntheses of mono-capped [PMo(12)O(40){Co(MeCN)(2)}](3-) and bi-capped [PMo(12)O(40)(VO)(2)](3-) and [PMo(12)O(40)Sb(2)](3-).     

Keggin polyoxoanions in aqueous solution: ion pairing and its effect on dynamic properties by molecular dynamics simulations

The dynamics of Keggin polyoxoanions in aqueous solution in the presence of monovalent cations is analyzed through molecular dynamics simulations. Together with structural information yielding the radial distribution functions of Li(+), Na(+), and K(+) with three polyoxometalates (POMs) bearing 3-, 4-, and 5- charges, the diffusion coefficient of these POMs is calculated. We found that the effect of the microscopic molecular details of the solvent is a key aspect to interpreting the structural and dynamic data because a competition between electrostatic interactions between the ions and the stability of the solvation shell is established. Furthermore, we show that solvent-shared structures weakly bound to the POM anion play a crucial role in the determination of the dynamic properties of the anion. The nature of these ion pairs, structurally characterized for the first time, is consistent with experimental data available.     

Fabrication of bulk-modified carbon paste electrode containing alpha-PW12O40(3-) polyanion supported on modified silica gel: Preparation, electrochemistry and electrocatalysis

Alpha-PW(12)O(40)(3-) (PW(12)) supported on the surface of silica gel derivatized by 3-aminopropyl(triethoxy)silane (devoted briefly as SiNH(3)PW(12)) was synthesized and used as bulk modifier to fabricate a renewable three-dimensional chemically modified electrode. The electrochemical behavior of the modified electrode was characterized by cyclic voltammetry. There is an ionic bonding character between PW(12) and the surface amino groups of modified silica, which greatly improves the stability of SiNH(3)PW(12)-modified carbon paste electrode due to insolubility of silica gel in water. The SiNH(3)PW(12) bulk-modified carbon paste electrode not only maintains the electrochemical activity of PW(12), but also exhibits remarkable advantages of renewability, as well as simple preparation and inexpensive material. The modified electrode offers an excellent and stable electrocatalytic response for the reduction of IO(3)(-) and hydrogen peroxide. The SiNH(3)PW(12)-CPE is successfully applied as an electrochemical detector to monitor IO(3)(-) in flow injection analysis (FIA). The catalytic peak current was found to be linear with the IO(3)(-) concentration in the range 5x10(-6) to 1x10(-3)molL(-1). The detection limit of the proposed method was found to be 3.1x10(-6)molL(-1) for IO(3)(-) determination.     

Keggin结构多金属氧酸盐α-[SiW12O40]4-的从头算理论研究

采用从头算Hartree-Fock方法对具有Keggin结构的多金属氧酸盐离子α-[SiW12O40]4-进行优化. 以基态几何为基础, 进行单激发组态相互作用(CIS)计算. 结果表明, WO6单元形成扭曲的八面体, 而SiO4部分仍具有准四面体结构; 阴离子的4个三金属簇W3O13可容纳3～4个负电荷, 从微观结构上揭示了杂多阴离子拥有大量负电荷仍能稳定的原因.     

Core-crosslinked block copolymer nanorods as templates for grafting [SiMo(12)O(40)](4-) Keggin ions

Core-crosslinked PB-P2VP block copolymer nanorods are used as templates for the synthesis of Keggin-type heteropolyoxometalate (POM) nanostructures by grafting [SiMo(12)O(40)](4-) Keggin ions on the template.     

Oxidations of hypophosphorus and arsenious acids by 12-tungstocobaltate(III) anion in aqueous solution

Summary The kinetics of oxidation of hypophosphorus and arsenious acids by 12-tungstocobaltate(III) anion have been studied in aqueous hydrochloric acid at constant ionic strength (I=2.0 M NaCl). The reactions obey the second-order rate law d[oxidant]/dt=2k [oxidant] [reductant]. Variation of [H⁺] in the range 0.10–1.50 M has no effect on the rates. Possible mechanistic interpretations of these observations are suggested.     

Molecular dynamics and Monte Carlo simulations for the structure of the aqueous trimethylammonium chloride solution in the 0.2-1 molar range

Molecular dynamics and Monte Carlo simulations have been performed for characterizing the structure of the 0.2 and 1 molar aqueous trimethylammonium chloride solutions. Atomic charges were derived through the CHELPG and RESP fits to the molecular electrostatic potentials calculated for the cation in water at the IEF-PCM/B3LYP level using the 6-31G* and 6-311++G** basis sets. Maxima and minima of the calculated radial distribution functions were not sensitive to the four atomic charge sets. Simulated solution structures suggest non-negligible solute-solute interactions and remarkable inhomogeneity at both concentrations. This means that equilibrium ratios, derived for conformers/tautomers by means of ab initio calculations with the IEF-PCM continuum dielectric solvent model, should be corrected for free energy changes following solute association when compared to experimental data obtained for the 0.1-1 molar aqueous solutions.     

Photoelectron spectroscopy of free multiply charged Keggin anions alpha-[PM12O40]3- (M = Mo, W) in the gas phase

Two polyoxometalate Keggin-type anions, alpha-PM12O40(3-) (M = Mo, W), were transferred to the gas phase by electrospray; their electronic structure and stability were probed by photoelectron spectroscopy. These triply charged anions were found to be highly stable in the gas phase with large adiabatic electron detachment energies of 1.7 and 2.1 eV for M = Mo and W, respectively. The magnitude of the repulsive Coulomb barrier was measured as approximately 3.4 eV for both anions, providing an experimental estimate for the intramolecular Coulomb repulsion present in these highly charged anions. Density functional theory calculations were carried out and compared with the experimental data, providing insight into the electronic structure and valence molecular orbitals of the two Keggin anions. The calculations indicated that the highest occupied molecular orbital and other frontier orbitals for PM12O40(3-) are localized primarily on the mu2-oxo bridging ligands of the polyoxometalate framework, consistent with the reactivity on the mu2-oxo sites observed in solution. It was shown that the HOMO of PW12O40(3-) is stabilized relative to that of PMo12O40(3-) by approximately 0.35 eV. The experimental adiabatic electron detachment energies of PM12O40(3-) (i.e., the electron affinities of PM12O40(2-)) are combined with recent calculations on the proton affinity of PM12O40(3-) to yield O-H bond dissociation energies in PM12O39(OH)2- as approximately 5.1 eV.     

Molecular dynamics simulations of end-to-end contact formation in hydrocarbon chains in water and aqueous urea solution

We probe the urea-denaturation mechanism using molecular dynamics simulations of an elementary "folding" event, namely, the formation of end-to-end contact in the linear hydrocarbon chain (HC) CH(3)(CH(2))(18)CH(3). Electrostatic effects are examined using a model HC in which one end of the chain is positively charged (+0.2e) and the other contains a negative charge (-0.2e). For these systems multiple transitions between "folded" (conformations in which the chain ends are in contact) and "unfolded" (end-to-end contact is broken) can be observed during 4 ns molecular dynamics simulations. In water and 6 M aqueous urea solution HC and the charged HC fluctuate between collapsed globular conformations and a set of expanded structures. The collapsed conformation adopted by the HC in water is slightly destablized in 6 M urea. In contrast, the end-to-end contact is disrupted in the charged HC only in aqueous urea solution. Despite the presence of a large hydrophobic patch, on length scales on the order of approximately 8-10 A "denaturation" (transition to the expanded unfolded state) occurs by a direct interaction of urea with charges on the chain ends. The contiguous patch of hydrophobic moieties leads to "mild dewetting", which becomes more pronounced in the charged HC in 6 M aqueous urea solution. Our simulations establish that the urea denaturation mechanism is most likely electrostatic in origin.     

Molecular dynamics simulations of sorption of organic compounds at the clay mineral/aqueous solution interface

The adsorption of trichloroethene, C₂HCl₃, on clay mineral surfaces in the presence of water has been modeled as an example describing a general program that uses molecular dynamics simulations to study the sorption of organic materials at the clay mineral/aqueous solution interface. Surfaces of the clay minerals kaolinite and pyrophyllite were hydrated at different water levels corresponding to partial and complete monolayers of water. In agreement with experimental trends, water was found to outcompete C₂HCl₃ for clay surface sites. The simulations suggest that at least three distinct mechanisms coexist for C₂HCl₃ on clay minerals in the environment. The most stable interaction of C₂HCl₃ with clay surfaces is by full molecular contact, coplanar with the basal surface. This kind of interaction is suppressed by increasing water loads. A second less stable and more reversible interaction involves adsorption through single-atom contact between one Cl atom and the surface. In a third mechanism, adsorbed C₂HCl₃ never contacts the clay directly but sorbs onto the first water layer. To test the efficacy of existing force field parameters of organic compounds in solid state simulations, molecular dynamics simulations of several representative organic crystals were also performed and compared with the experimental crystal structures. These investigations show that, in general, in condensed-phase studies, parameter evaluations are realistic only when thermal motion effects are included in the simulations. For chlorohydrocarbons in particular, further explorations are needed of atomic point charge assignments. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 144–153, 1998     

A DFT study of the electronic spectrum of the alpha-keggin anion [Co(II)W(12)O(40)](6-)

Multiplet splittings for several excited configurations of [Co(II)W(12)O(40)](6-) were calculated using DFT methods. In agreement with the experimental interpretation of the spectrum the calculations found that the first strong band corresponds to Co d-d transitions, but it is worth noting that superposed to these transitions there are charge transfer transitions from cobalt to tungsten. The calculations also showed the importance of Jahn-Teller distortions in the excited states. With the exception of the consequences derived from a smaller splitting of d cobalt orbitals the d-d spectrum of [CoCl(4)](2-) is similar to that of the more complex Keggin anion. Finally, the energy of the bielectronic transition (4)A(2) --> (4)T(1)(P) was estimated via an approximate procedure based on ligand field theory.     

Molecular dynamics simulations of La2O3 thin films on SiO2

Classical molecular dynamics is used to investigate the equilibrium state of the surface region and interface of heteroepitaxial La2O3 thin films.Due to the lattice mismatch,heteroepitaxial thin films are subject to very large stress.For this reason the behavior of La2O3 thin films at SiO2interface becomes an important concern.Our result indicates that La2O3 can uniformly wet SiO2 surface.The properties of the simulated films are analyzed and the lack of any discernible crystalline phase in epitaxial La2O3 on SiO2 indicates that the lattice mismatch between SiO2 and La2O3 is sufficiently large to prevent the formation of even short-range orders in La2O3 film.     

The self-assembly mechanism of the Lindqvist anion [W(6)O(19)](2-) in aqueous solution: a density functional theory study

The formation mechanism is always a fundamental and confused issue for polyoxometalate chemistry. Two formation mechanisms (M1 and M2) of the Lindqvist anion [W(6)O(19)](2-) have been adopted to investigate it's self-assembly reaction pathways at a density functional theory (DFT) level. The potential energy surfaces reveal that both the mechanisms are thermodynamically favorable and overall barrierless at room temperature, but M2 is slightly dominant to M1. The formation of the pentanuclear species [W(5)O(16)](2-) and [W(5)O(15)(OH)](-) are recognized as the rate-determining steps in the whole assembly polymerization processes. These two steps involve the highest energy barriers with 30.48 kcal mol(-1) and 28.90 kcal mol(-1), respectively, for M1 and M2. [W(4)O(13)](2-) and [W(4)O(12)(OH)](-) are proved to be the most stable building blocks. In addition, DFT results reveal that the formation of [W(3)O(10)](2-) experiences a lower barrier along the chain channel.     

Mechanistic analysis of the electrocatalytic properties of dissolved alpha and beta isomers of [SiW12O40]4- and solid [Ru(bipy)3]2[alpha-SiW12O40] on the reduction of nitrite in acidic aqueous media

Voltammetric studies on the reduction of alpha and beta isomers of the Keggin polyoxometalate anion [SiW12O40]4- reveal a series of electrochemically reversible processes in acidic aqueous media. In the presence of NO2-, catalytic current is detected in the potential region of the [SiW12O40]4-/5- process. Electronic spectroscopy and simulation of voltammetric data undertaken at variable [NO2-] and [H+] allow the following mechanism to be postulated, [SiW12O40]4- + e- <-->[SiW12O40]5-, H+ + HNO2 <--> NO+ + H2O, NO+ + [SiW12O40]5- --> NO + [SiW12O40]4-. The second-order rate constant for the rate-determining step is faster for the alpha isomer than for the beta one. This may be attributed to the different reversible potentials of -0.144 V (alpha isomer) and -0.036 V vs Ag/AgCl (beta isomer) and, hence, smaller driving force for an assumed outer sphere electron-transfer reaction in the case of beta isomer. A stable, water-insoluble, thin-film [Ru(bipy)3]2[alpha-SiW12O40] chemically modified electrode was generated electrochemically via ion-exchange of [Ru(bipy)3]2+ with Bu4N+ in the [Bu4N]4[alpha-SiW12O40] solid. The first reduction process with this modified electrode gives rise to the reaction [Ru(bipy)3]2[alpha-SiW12O40](solid) + H+(soln) + e- <--> H[Ru(bipy)3]2[alpha-SiW12O40](solid). The need to transfer a proton from the solution to the solid phase for charge neutralization purposes introduces a hydrogen-ion concentration dependence into this reaction, which is not found in the solution-phase study. Nevertheless, the voltammetric catalytic activity with respect to nitrite reduction is retained with the chemically modified electrode. However, nitrite catalysis with the [Ru(bipy)3]2[alpha-SiW12O40]-modified electrode is now independent of concentration of H+, rather than exhibiting a first-order dependence, and full mechanistic details for this process are unknown.