The impact of stoichiometry on the initial steps of crystal formation: Stability and lifetime of charged triple-ion complexes.

Minerals form in natural systems from solutions with varying ratios of their lattice ions, yet non-stoichiometric conditions have generally been overlooked in investigations of new formation (nucleation) of ionic crystals. Here, we investigated the influence of cation:anion ratio in the solution on the initial steps of nucleation by studying positively and negatively charged triple ion complexes and subsequent particle size evolution. Our model systems are carbonates and sulfates of calcium and barium, as it was recently shown that solution stoichiometry affects the timing and rate of their nucleation. Molecular dynamics (MD) simulations and dynamic light scattering (DLS) flow experiments show that nucleation correlates with the stability and lifetime of the initial complexes, which were significantly impacted by the cation:anion stoichiometry and ion type. Specifically, B a S O 4 2 2 - ${{\rm B}{\rm a}{\left({{\rm S}{\rm O}}_{4}\right)}_{2}^{2-}}$ was found to have higher association constants and its lifetime was twofold longer than B a 2 S O 4 2 + ${{{\rm B}{\rm a}}_{2}{{\rm S}{\rm O}}_{4}^{2+}}$ . Similar trends were observed for B a C O 3 ${{{\rm B}{\rm a}{\rm C}{\rm O}}_{3}}$ and C a S O 4 ${{{\rm C}{\rm a}{\rm S}{\rm O}}_{4}}$ . Contrastingly, for C a C O 3 ${{{\rm C}{\rm a}{\rm C}{\rm O}}_{3}}$ , C a C O 3 2 2 - ${{\rm C}{\rm a}{\left({{\rm C}{\rm O}}_{3}\right)}_{2}^{2-}}$ was found to have lower association constants and its lifetime was shorter than C a 2 C O 3 2 + ${{{\rm C}{\rm a}}_{2}{{\rm C}{\rm O}}_{3}^{2+}}$ . These trends in stability and lifetime follow the same asymmetrical behaviour as observed experimentally for particle formation using techniques like DLS. This suggests a causal relationship between the stability and lifetime of the initial charged complexes and the nucleation under non-stoichiometric conditions.

A Database of Solution Additives Promoting Mg2+ Dehydration and the Onset of MgCO3 Nucleation.

Formed via aqueous carbonation of Mg2+ ions, the crystallization of magnesite (MgCO3) is a promising route to carbon capture and reuse, albeit limited by the slow precipitation of MgCO3. Although magnesite is naturally abundant, forming at low temperature conditions, its industrial production is an energy-intensive process due to the temperatures required to prevent the formation of hydrated phases. The principal difficulty in aqueous conditions arises from the very strong Mg2+···H2O interaction, with high barriers to Mg2+ dehydration. Using atomistic simulations, we have investigated the influence of 30 additive anions (X n-, n = 1-3), ranging from simple halides to more complex molecules, on the first two steps of MgCO3 aggregation from solution, as follows: Mg2+ dehydration and subsequent prenucleative Mg2+···CO3 2- pairing. We have computed the thermodynamic stabilities of solvent shared ion pairs (Mg2+···H2O···X n-) and contact ion pairs (Mg2+···X n-) to reveal the propensity of solution additives to inhibit or promote Mg2+···CO3 2- formation. We have determined the stabilization of undercoordinated hydrated Mg2+ states with a vacant coordination site to which CO3 2- can bind, subsequently initiating MgCO3 nucleation or Mg2+ incorporation into the crystal lattice. Extensive molecular dynamics simulations of electrolyte solutions containing Na2CO3 with different sources of Mg2+ (i.e., MgCl2, MgSO4, and Mg(CH3COO)2) further show that the degree of dehydration of Mg2+ and the structure of prenucleation MgCO3 clusters change depending on the counterion identity. Through a fundamental understanding of the role of solution additives in the mechanism of Mg2+ dehydration, our results help to rationalize previously reported experimental observation of the effect of solvation environments on the growth of magnesite. This understanding may contribute to identifying the solution composition and conditions that could promote the low-temperature CO2 conversion into MgCO3 at industrially relevant scales.

Density functional theory based molecular dynamics study of solution composition effects on the solvation shell of metal ions.

We present an ab initio molecular dynamics study of the alkali metal ions Li+, Na+, K+ and Cs+, and of the alkaline earth metal ions Mg2+ and Ca2+ in both pure water and electrolyte solutions containing the counterions Cl- and SO42-. Simulations were conducted using different density functional theory methods (PBE, BLYP and revPBE), with and without the inclusion of dispersion interactions (-D3). Analysis of the ion-water structure and interaction strength, water exchange between the first and second hydration shell, and hydrogen bond network and low-frequency reorientation dynamics around the metal ions have been used to characterise the influence of solution composition on the ionic solvation shell. Counterions affect the properties of the hydration shell not only when they are directly coordinated to the metal ion, but also when they are at the second coordination shell. Chloride ions reduce the sodium hydration shell and expand the calcium hydration shell by stabilizing under-coordinated hydrated Na(H2O)5+ complexes and over-coordinated Ca(H2O)72+. The same behaviour is observed in CaSO4(aq), where Ca2+ and SO42- form almost exclusively solvent-shared ion pairs. Water exchange between the first and second hydration shell around Ca2+ in CaSO4(aq) is drastically decelerated compared with the simulations of the hydrated metal ion (single Ca2+, no counterions). Velocity autocorrelation function analysis, used to probe the strength of the local ion-water interaction, shows a smoother decay of Mg2+ in MgCl2(aq), which is a clear indication of a looser inter-hexahedral vibration in the presence of chloride ions located in the second coordination shell of Mg2+. The hydrogen bond statistics and orientational dynamics in the ionic solvation shell show that the influence on the water-water network cannot only be ascribed to the specific cation-water interaction, but also to the subtle interplay between the level of hydration of the ions, and the interactions between ions, especially those of opposite charge. As many reactive processes involving solvated metal ions occur in environments that are far from pure water but rich in ions, this computational study shows how the solution composition can result in significant differences in behaviour and function of the ionic solvation shell.

Mechanism of Photoinduced Triplet Intermolecular Hydrogen Transfer between Cycloxydim and Chlorothalonil.

The possible reaction mechanisms for the experimentally observed hydrogen transfer between the herbicide cycloxydim (CD) and the triplet fungicide chlorothalonil (CT) were identified with density functional theory (DFT) and time-dependent density function theory (TDDFT) computations. Excited energy transfer (EET) calculations indicate that reactants for intermolecular hydrogen transfer were formed via energy transfer from triplet CT to ground state CD. Three possible reaction pathways after EET were identified, and hydrogen transfer from the hydroxyl group on the cyclohexane ring of CD to CT exhibited the lowest energy barrier. Natural population analysis (NPA) along the reaction pathways has confirmed that the pathways involved either electron transfer induced proton transfer or coupled electron-proton transfer, leading to different potential energy profiles. Electrostatic potential (ESP) study substantiated the reaction mechanisms in different pathways. This study suggests an explanation for the accelerated photodegradation of CD by CT and provides a pipeline for future studies of photoinduced intermolecular hydrogen transfer.

Can small hydrophobic gold nanoparticles inhibit ß2-microglobulin fibrillation?

Inorganic nanoparticles stabilized by a shell of organic ligands can enhance or suppress the natural propensity of proteins to form fibrils. Functionalization facilitates targeted delivery of the nanoparticles to various cell types, bioimaging, drug delivery and other therapeutic and diagnostic applications. In this study, we provide a computational model of the effect of a prototypical thiol-protected gold nanoparticle, Au25L18(-) (L = S(CH2)2Ph) on the ß2-microglobulin natural fibrillation propensity. To reveal the molecular basis of the protein-nanoparticle association process, we performed various simulations at multiple levels (Classical Molecular Dynamics and Brownian Dynamics) that cover multiple length- and timescales. The results provide a model of the ensemble of structures constituting the protein-gold nanoparticle complexes, and insights into the driving forces for the binding of ß2-microglobulin to hydrophobic small size gold nanoparticles. We have found that the small nanoparticles can bind the protein to form persistent complexes. This binding of nanoparticles is able to block the active sites of domains from binding to another protein, thus leading to potential inhibition of the fibrillation activity. A comparison with the binding patches identified for the interaction of the protein with a known inhibitor of fibrillation, supports our conclusion.

Proposed alteration of images of molecular orbitals obtained using a scanning tunneling microscope as a probe of electron correlation.

Scanning tunneling spectroscopy (STS) allows us to image single molecules decoupled from the supporting substrate. The obtained images are routinely interpreted as the square moduli of molecular orbitals, dressed by the mean-field electron-electron interaction. Here we demonstrate that the effect of electron correlation beyond the mean field qualitatively alters the uncorrelated STS images. Our evidence is based on the ab initio many-body calculation of STS images of planar molecules with metal centers. We find that many-body correlations alter significantly the image spectral weight close to the metal center of the molecules. This change is large enough to be accessed experimentally, surviving to molecule-substrate interactions.

Peptide synthesis of gold nanoparticles: the early steps of gold reduction investigated by density functional theory.

Gold nanoparticles can be synthesized by reducing chloroaurate(III) ions in the presence of peptides. Here, such reduction for serine and tyrosine is studied by density functional theory including solvent effects. We find that the formation of chloroaurate complexes of these amino acids is thermodynamically viable and facilitates the reduction of Au(III), to a greater degree for tyrosine as found in experiments. Our results also suggest a rationale for the behavior of tyrosine-intercalated peptides.

Visualizing electron correlation by means of ab initio scanning tunneling spectroscopy images of single molecules.

Scanning tunneling microscopy (STM) has been a fundamental tool to characterize many-body effects in condensed matter systems, from extended solids to quantum dots. STM of molecules decoupled from the supporting conductive substrate has the potential to extend STM characterization of many-body effects to the molecular world as well. In this paper, we describe a many-body tunneling theory for molecules decoupled from the STM substrate, and we report on the use of standard quantum chemical methods to calculate the quantities necessary to provide the "correlated" STM molecular image. The developed approach has been applied to 18 different molecules to explore the effects of their chemical nature and of their substituents, as well as to verify the possible contribution by transition metal centers. Whereas the bulk of calculations has been performed with the configuration interaction method with single and double excitations (CISD), because of the computational cost some tests have been also performed with the more accurate coupled cluster with single and double excitations (CCSD) method to quantify the importance of the computational level on many-body STM images. We have found that correlation induces a remarkable squeezing of the images, and that correlated images are not derived from Hartree-Fock HOMO or LUMO alone, but include contributions from other orbitals as well. Although correlation effects are too small to be resolved by present STM experiments for the studied molecules, our results provide hints for seeking out other species with larger, and possibly experimentally detectable, correlation effects.

Structure of the gas-phase glycine tripeptide.

The structural preferences of the neutral gas-phase glycine tripeptide have been investigated using a variety of strategies including a hierarchy of electronic structure theory (encompassing HF/3-21G single-point energy calculation and geometry optimisation, B3LYP/6-31G(d) geometry optimisation and MP2/6-31+G(d) single-point energy calculation and/or geometry optimisation). The structures and relative stabilities of the 20 most stable conformers identified were verified by M05-2X and mPW2-PLYP-D calculations. The most stable conformer located has a folded gamma-turn structure, with an NH...N interaction between the N-terminal nitrogen and the amide hydrogen of glycine (2) and an NH...O interaction between the amide hydrogen of glycine (3) and the carboxyl oxygen of glycine (1). The results show a clear preference for folded over extended structures.