Elucidating the Mechanism of Efficient Eu(III) and Yb(III) Sensitisation from a Re(I) Tetrazolato Triangular Assembly.

The reaction of Re(CO)5Br with deprotonated 1H-(5-(2,2':6',2''-terpyridine)pyrid-2-yl)tetrazole yields a triangular assembly formed by tricarbonyl Re(I) vertices. Photophysical measurements reveal blue-green emission with a maximum at 520 nm, 32% quantum yield, and 2430 ns long-lived excited state decay lifetime in deaerated dichloromethane solution. Coordination of lanthanoid ions to the terpyridine units red-shifts the emission to 570 nm and also reveals efficient (90%) and fast sensitisation to both Eu(III) and Yb(III) at room temperature, with a similar rate constant kET of the order of 107 s-1. Efficient sensitisation of Eu(III) from Re(I) is unprecedented, especially when considering the close proximity in energy between the donor and acceptor excited states. On the other hand, comparative measurements at 77 K reveal that energy transfer to Yb(III) is two orders of magnitude slower than that to Eu(III). A two-step mechanism of sensitisation is therefore proposed, whereby the rate-determining step is a thermally activated energy transfer step between the Re(I) centre and the terpyridine functionality, followed by rapid energy transfer to the respective Ln(III) excited states. At 77 K, the direct Re(I) to Eu(III) energy transfer seems to proceed via a ligand-mediated superexchange Dexter-type mechanism.

New model for aspartic acid species in aqueous calcium carbonate growth environments: challenges and perspectives.

The lack of experimental data on the dynamics of aspartic acid species in water for its range of protonation states and the details of their atomic-level interaction with aqueous calcium carbonate species is a driver for accurate force field development. A classical model that is consistent with the few pieces of experimental data available and with first principles calculations has been developed. The complex dynamics of the aspartate anions relevant to biomineralization and calcium carbonate crystal growth has been explored in water, providing a quantitative description of solvation structure and free energies, including conformational free energy profiles and pairing free energies. The model has been used to probe the structure and dynamics of aqueous calcium aspartate homo- and hetero-chiral clusters, confirming their unlikelihood due to weak and water-mediated interactions. This supports the hypothesis that the formation of such clusters, observed while growing vaterite in the presence of acidic chiral amino acids, is favoured by the presence of the crystal surface.

Solubility-consistent force field simulations for aqueous metal carbonate systems using graphical processing units.

Crystallization of alkaline earth metal carbonates from water is important for biomineralization and environmental geochemistry. Here, large-scale computer simulations are a useful approach to complement experimental studies by providing atomistic insights and even by quantitatively determining the thermodynamics of individual steps. However, this is dependent on the existence of force field models that are sufficiently accurate while being computationally efficient enough to sample complex systems. Here, we introduce a revised force field for aqueous alkaline earth metal carbonates that reproduces both the solubilities of the crystalline anhydrous minerals, as well as the hydration free energies of the ions. The model is also designed to run efficiently on graphical processing units thereby reducing the cost of such simulations. The performance of the revised force field is compared against previous results for important properties relevant to crystallization, including ion-pairing and mineral-water interfacial structure and dynamics. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

Computational workflows for perovskites: case study for lanthanide manganites.

Robust computational workflows are important for explorative computational studies, especially for cases where detailed knowledge of the system structure or other properties is not available. In this work, we propose a computational protocol for appropriate method selection for the study of lattice constants of perovskites using density functional theory, based strictly on open source software. The protocol does not require a starting crystal structure. We validate this protocol using a set of crystal structures of lanthanide manganites, surprisingly finding N12+U to be the best performing method for this class of materials out of the 15 density functional approximations studied. We also highlight that +U values derived from linear response theory are robust and their use leads to improved results. We investigate whether the performance of methods for predicting the bond length of related gas phase diatomics correlates with their performance for bulk structures, showing that care is required when interpreting benchmark results. Finally, using defective LaMnO3 as a case study, we investigate whether the four shortlisted methods (HCTH120, OLYP, N12+U, PBE+U) can computationally reproduce the experimentally determined fraction of MnIV+ at which the orthorhombic to rhombohedral phase transition occurs. The results are mixed, with HCTH120 providing good quantitative agreement with experiment, but failing to capture the spatial distribution of defects linked to the electronic structure of the system.

Structure and Surface Complexation at the Calcite(104)-Water Interface.

Calcite is the most stable polymorph of calcium carbonate (CaCO3) under ambient conditions and is ubiquitous in natural systems. It plays a major role in controlling pH in environmental settings. Electrostatic phenomena at the calcite-water interface and the surface reactivity of calcite in general have important environmental implications. They may strongly impact nutrient and contaminant mobility in soils and other subsurface environments, they control oil recovery from limestone reservoirs, and they may impact the safety of nuclear waste disposal sites. Besides the environmental relevance, the topic is significant for industrial applications and cultural heritage preservation. In this study, the structure of the calcite(104)-water interface is investigated on the basis of a new extensive set of crystal truncation rod data. The results agree with recently reported structures and resolve previous ambiguities with respect to the coordination sphere of surface Ca ions. These structural features are introduced into an electrostatic three-plane surface complexation model, describing ion adsorption and charging at the calcite-water interface. Inner surface potential data for calcite, as measured with a calcite single-crystal electrode, are used as constraints for the model in addition to zeta potential data. Ion adsorption parameters are compared with molecular dynamics simulations. All model parameters, including protonation constants, ion-binding parameters, and Helmholtz capacitances, are within physically and chemically plausible ranges. A PhreeqC version of the model is presented, which we hope will foster application of the model in environmental studies.

Simulating the binding of key organic functional groups to aqueous calcium carbonate species.

The interaction of organic molecules with mineral systems is relevant to a wide variety of scientific problems both in the environment and minerals processing. In this study, the coordination of small organics that contain the two most relevant functional groups for biomineralisation of calcium carbonate, namely carboxylate and ammonium, with the corresponding mineral ions are examined in aqueous solution. Specifically, two force fields have been examined based on rigid-ion or polarisable models, with the latter being within the AMOEBA formalism. Here the parameters for the rigid-ion model are determined to target the accurate reproduction of the hydration structure and solvation thermodynamics, while both force fields are designed to be compatible with the corresponding recently published models for aqueous calcium carbonate. The application of these force fields to ion pairing in aqueous solution is studied in order to quantitatively determine the extent of association.

The atomic structure and dynamics at the CaCO3 vaterite-water interface: A classical molecular dynamics study.

Classical molecular and lattice dynamics were applied to explore the structure and dynamics of water on different surfaces of vaterite, the least abundant calcium carbonate polymorph. Surfaces were generated starting from the three possible structural models for vaterite (monoclinic, hexagonal/trigonal, and triclinic) and pre-screened using their surface energies in an implicit solvent. Surfaces with energies lower than 0.55 J/m2 were then run in explicit water. The majority of these surfaces dissolve in less than 100 ns, highlighting the low stability of this phase in abiotic environments. Three stable surfaces were identified; they exhibited only minor structural changes when in contact with explicit water and did not show any tendency to dissolve during 1 µs molecular dynamics simulations. The computed water density profiles show that all these surfaces have two distinct hydration layers. The water residence time at the various calcium sites was computed to be within 0.7 and 20.5 ns, which suggests that specific Ca ions will be more readily available to bind with organic molecules present in solution. This analysis is a step forward in understanding the structure of this complex mineral and its role in biomineralization, as it provides a solid theoretical background to explore its surface chemistry. In particular, this study provides realistic surface models and predicts the effect of water exchange at the surface active sites on the adsorption of other molecules.

Uncovering the Role of Bicarbonate in Calcium Carbonate Formation at Near-Neutral pH.

Mechanistic pathways relevant to mineralization are not well-understood fundamentally, let alone in the context of their biological and geological environments. Through quantitative analysis of ion association at near-neutral pH, we identify the involvement of HCO3- ions in CaCO3 nucleation. Incorporation of HCO3- ions into the structure of amorphous intermediates is corroborated by solid-state nuclear magnetic resonance spectroscopy, complemented by quantum mechanical calculations and molecular dynamics simulations. We identify the roles of HCO3- ions as being through (i) competition for ion association during the formation of ion pairs and ion clusters prior to nucleation and (ii) incorporation as a significant structural component of amorphous mineral particles. The roles of HCO3- ions as active soluble species and structural constituents in CaCO3 formation are of fundamental importance and provide a basis for a better understanding of physiological and geological mineralization.

Molecular dynamics simulations of liquid-liquid interfaces in an electric field: The water-1,2-dichloroethane interface.

The polarized interface between two immiscible liquids plays a central role in many technological processes. In particular, for electroanalytical and ion extraction applications, an external electric field is typically used to selectively induce the transfer of ionic species across the interfaces. Given that it is experimentally challenging to obtain an atomistic insight into the ion transfer process and the structure of liquid-liquid interfaces, atomistic simulations have often been used to fill this knowledge gap. However, due to the long-range nature of the electrostatic interactions and the use of 3D periodic boundary conditions, the use of external electric fields in molecular dynamics simulations requires special care. Here, we show how the simulation setup affects the dielectric response of the materials and demonstrate how by a careful design of the system it is possible to obtain the correct electric field on both sides of a liquid-liquid interface when using standard 3D Ewald summation methods. In order to prove the robustness of our approach, we ran extensive molecular dynamics simulations with a rigid-ion and polarizable force field of the water/1,2-dichloroethane interface in the presence of weak external electric fields.

Resolving Point Defects in the Hydration Structure of Calcite (10.4) with Three-Dimensional Atomic Force Microscopy.

It seems natural to assume that defects at mineral surfaces critically influence interfacial processes such as the dissolution and growth of minerals in water. The experimental verification of this claim, however, is challenging and requires real-space methods with utmost spatial resolution, such as atomic force microscopy (AFM). While defects at mineral-water interfaces have been resolved in 2D AFM images before, the perturbation of the surrounding hydration structure has not yet been analyzed experimentally. In this Letter, we demonstrate that point defects on the most stable and naturally abundant calcite (10.4) surface can be resolved using high-resolution 3D AFM-even within the fifth hydration layer. Our analysis of the hydration structure surrounding the point defect shows a perturbation of the hydration with a lateral extent of approximately one unit cell. These experimental results are corroborated by molecular dynamics simulations.

Investigating Intracellular Localisation and Cytotoxicity Trends for Neutral and Cationic Iridium Tetrazolato Complexes in Live Cells.

A family of five neutral cyclometalated iridium(III) tetrazolato complexes and their methylated cationic analogues have been synthesised and characterised. The complexes are distinguished by variations of the substituents or degree of π conjugation on either the phenylpyridine or tetrazolato ligands. The photophysical properties of these species have been evaluated in organic and aqueous media, revealing predominantly a solvatochromic emission originating from mixed metal-to-ligand and ligand-to-ligand charge transfer excited states of triplet multiplicity. These emissions are characterised by typically long excited-state lifetimes (∼hundreds of ns), and quantum yields around 5-10 % in aqueous media. Methylation of the complexes caused a systematic red-shift of the emission profiles. The behaviour and the effects of the different complexes were then examined in cells. The neutral species localised mostly in the endoplasmic reticulum and lipid droplets, whereas the majority of the cationic complexes localised in the mitochondria. The amount of complexes found within cells does not depend on lipophilicity, which potentially suggests diverse uptake mechanisms. Methylated analogues were found to be more cytotoxic compared to the neutral species, a behaviour that might to be linked to a combination of uptake and intracellular localisation.

Uncovering the Atomistic Mechanism for Calcite Step Growth.

Determining a complete atomic-level picture of how minerals grow from aqueous solution remains a challenge as macroscopic rates can be a convolution of many reactions. For the case of calcite (CaCO3 ) in water, computer simulations have been used to map the complex thermodynamic landscape leading to growth of the two distinct steps, acute and obtuse, on the basal surface. The carbonate ion is found to preferentially adsorb at the upper edge of acute steps while Ca2+ only adsorbs after CO32- . In contrast to the conventional picture, ion pairs prefer to bind at the upper edge of the step with only one ion, at most, coordinated to the step and lower terrace. Migration of the first carbonate ion to a growth site is found to be rate-limiting for kink nucleation, with this process having a lower activation energy on the obtuse step.

A Supramolecular Ice Growth Inhibitor.

Safranine O, a synthetic dye, was found to inhibit growth of ice at millimolar concentrations with an activity comparable to that of highly evolved antifreeze glycoproteins. Safranine inhibits growth of ice crystals along the crystallographic a-axis, resulting in bipyramidal needles extended along the <0001> directions as well as and plane-specific thermal hysteresis (TH) activity. The interaction of safranine with ice is reversible, distinct from the previously reported behavior of antifreeze proteins. Spectroscopy and molecular dynamics indicate that safranine forms aggregates in aqueous solution at micromolar concentrations. Metadynamics simulations and aggregation theory suggested that as many as 30 safranine molecules were preorganized in stacks at the concentrations where ice growth inhibition was observed. The simulations and single-crystal X-ray structure of safranine revealed regularly spaced amino and methyl substituents in the aggregates, akin to the ice-binding site of antifreeze proteins. Collectively, these observations suggest an unusual link between supramolecular assemblies of small molecules and functional proteins.

Resorcinol Crystallization from the Melt: A New Ambient Phase and New "Riddles".

Structures of the α and ß phases of resorcinol, a major commodity chemical in the pharmaceutical, agrichemical, and polymer industries, were the first polymorphic pair of molecular crystals solved by X-ray analysis. It was recently stated that "no additional phases can be found under atmospheric conditions" (Druzbicki, K. et al. J. Phys. Chem. B 2015, 119, 1681). Herein is described the growth and structure of a new ambient pressure phase, ε, through a combination of optical and X-ray crystallography and by computational crystal structure prediction algorithms. α-Resorcinol has long been a model for mechanistic crystal growth studies from both solution and vapor because prisms extended along the polar axis grow much faster in one direction than in the opposite direction. Research has focused on identifying the absolute sense of the fast direction-the so-called "resorcinol riddle"-with the aim of identifying how solvent controls crystal growth. Here, the growth velocity dissymmetry in the melt is analyzed for the ß phase. The ε phase only grows from the melt, concomitant with the ß phase, as polycrystalline, radially growing spherulites. If the radii are polar, then the sense of the polar axis is an essential feature of the form. Here, this determination is made for spherulites of ß resorcinol (ε, point symmetry 222, does not have a polar axis) with additives that stereoselectively modify growth velocities. Both ß and ε have the additional feature that individual radial lamellae may adopt helicoidal morphologies. We correlate the appearance of twisting in ß and ε with the symmetry of twist-inducing additives.

Can Point Defects in Surfaces in Solution be Atomically Resolved by Atomic Force Microscopy?

While the atomic force microscope (AFM) is able to image mineral surfaces in solution with atomic resolution, so far, it has been a matter of debate whether imaging point defects is also possible under these conditions. The difficulties stem from the limited knowledge of what types of defects may be stable in the presence of an AFM tip, as well as from the complicated imaging mechanism involving interactions between hydration layers over the surface and around the tip apex. Here, we present atomistic molecular dynamics and free energy calculations of the AFM imaging of vacancies and ionic substitutions in the calcite (101[over ¯]4) surface in water, using a new silica AFM tip model. Our results indicate that both calcium and carbonate vacancies, as well as a magnesium substitution, could be resolved in an AFM experiment, albeit with different imaging mechanisms.

Entropy Drives Calcium Carbonate Ion Association.

The understanding of the molecular mechanisms underlying the early stages of crystallisation is still incomplete. In the case of calcium carbonate, experimental and computational evidence suggests that phase separation relies on so-called pre-nucleation clusters (PNCs). A thorough thermodynamic analysis of the enthalpic and entropic contributions to the overall free energy of PNC formation derived from three independent methods demonstrates that solute clustering is driven by entropy. This can be quantitatively rationalised by the release of water molecules from ion hydration layers, explaining why ion association is not limited to simple ion pairing. The key role of water release in this process suggests that PNC formation should be a common phenomenon in aqueous solutions.

Photophysical and photochemical trends in tricarbonyl rhenium(I) N-heterocyclic carbene complexes.

A family of tricarbonyl Re(I) complexes of the formulation fac-[Re(CO)3(NHC)L] has been synthesized and characterized, both spectroscopically and structurally. The NHC ligand represents a bidentate N-heterocyclic carbene species where the central imidazole ring is substituted at the N3 atom by a butyl, a phenyl, or a mesityl group and substituted at the N1 atom by a pyridyl, a pyrimidyl, or a quinoxyl group. On the other hand, the ancillary L ligand alternates between chloro and bromo. For the majority of the complexes, the photophysical properties suggest emission from the lowest triplet metal-to-ligand charge transfer states, which are found partially mixed with triplet ligand-to-ligand charge transfer character. The nature and relative energy of the emitting states appear to be mainly influenced by the identity of the substituent on the N3 atom of the imidazole ring; thus, the pyridyl complexes have blue-shifted emission in comparison to the more electron deficient pyrimidyl complexes. The quinoxyl complexes show an unexpected blue-shifted emission, possibly occurring from ligand-centered excited states. No significant variations were found upon changing the substituent on the imidazole N3 atom and/or the ancillary ligand. The photochemical properties of the complexes have also been investigated, with only the complexes bound to the pyridyl-substituted NHC ligands showing photoinduced CO dissociation upon excitation at 370 nm, as demonstrated by the change in the IR and NMR spectra as well as a red shift in the emission profile after photolysis. Temperature-dependent photochemical experiments show that CO dissociation occurs at temperatures as low as 233 K, suggesting that the Re-C bond cleaves from excited states of metal-to-ligand charge transfer nature rather than thermally activated ligand field excited states. A photochemical mechanism that takes into account the reactivity of the complexes bound to the pyridyl-NHC ligand as well as the stability of those bound to the pyrimidyl- and quinoxyl-NHC ligands is proposed.

Proton-induced reversible modulation of the luminescent output of rhenium(I), iridium(III), and ruthenium(II) tetrazolate complexes.

One of the distinct features of metal-tetrazolate complexes is the possibility of performing electrophilic additions onto the imine-type nitrogens of the coordinated five-membered ring. These reactions, in particular, provide a useful tool for varying the main structural and electronic properties of the starting tetrazolate complexes. In this paper, we demonstrate how the use of a simple protonation-deprotonation protocol enables us to reversibly change, to a significant extent, the light-emission output and performance of a series of Re(I)-tetrazolate-based phosphors of the general formulation fac-[Re(N(∧)N)(CO)3L], where N(∧)N denotes diimine-type ligands such as 2,2'-bipyridine (bpy) or 1,10-phenanthroline (phen) and L represents a series of different 5-aryl tetrazolates. Indeed, upon addition of triflic acid to these neutral Re(I) complexes, a consistent blue shift (Δλmax ca. 50 nm) of the emission maximum is observed and the protonated species also display increased quantum yield values (4-13 times greater than the starting compounds) and longer decay lifetimes. This alteration can be reversed to the initial condition by further treating the protonated Re(I) complex with a base such as triethylamine. Interestingly, the reversible modulation of luminescent features by the same protonation-deprotonation mechanism appears as a quite general characteristic of photoactive metal tetrazolate complexes, even for compounds in which the 2-pyridyl tetrazolate ligands coordinate the metal center with a bidentate mode, such as the corresponding Ir(III) cyclometalates [Ir(C(∧)N)2L] and the Ru(II) polypyridyl derivatives [Ru(bpy)2L](+). In these cases, the protonation of the starting materials leads to red-shifted and more intense emissions for the Ir(III) complexes, while almost complete quenching is observed in the case of the Ru(II) analogues.

Redefined ion association constants have consequences for calcium phosphate nucleation and biomineralization.

Calcium orthophosphates (CaPs), as hydroxyapatite (HAP) in bones and teeth are the most important biomineral for humankind. While clusters in CaP nucleation have long been known, their speciation and mechanistic pathways to HAP remain debated. Evidently, mineral nucleation begins with two ions interacting in solution, fundamentally underlying solute clustering. Here, we explore CaP ion association using potentiometric methods and computer simulations. Our results agree with literature association constants for Ca2+ and H2PO4-, and Ca2+ and HPO42-, but not for Ca2+ and PO43- ions, which previously has been strongly overestimated by two orders of magnitude. Our data suggests that the discrepancy is due to a subtle, premature phase separation that can occur at low ion activity products, especially at higher pH. We provide an important revision of long used literature constants, where association of Ca2+ and PO43- actually becomes negligible below pH 9.0, in contrast to previous values. Instead, [CaHPO4]0 dominates the aqueous CaP speciation between pH ~6-10. Consequently, calcium hydrogen phosphate association is critical in cluster-based precipitation in the near-neutral pH regime, e.g., in biomineralization. The revised thermodynamics reveal significant and thus far unexplored multi-anion association in computer simulations, constituting a kinetic trap that further complicates aqueous calcium phosphate speciation.

Crosslinking Rapidly Cured Epoxy Resin Thermosets: Experimental and Computational Modeling and Simulation Study.

The power of computational modeling and simulation for establishing clear links between materials' intrinsic properties and their atomic structure has more and more increased the demand for reliable and reproducible protocols. Despite this increased demand, no one approach can provide reliable and reproducible outcomes to predict the properties of novel materials, particularly rapidly cured epoxy-resins with additives. This study introduces the first computational modeling and simulation protocol for crosslinking rapidly cured epoxy resin thermosets based on solvate ionic liquid (SIL). The protocol combines several modeling approaches, including quantum mechanics (QMs) and molecular dynamics (MDs). Furthermore, it insightfully provides a wide range of thermo-mechanical, chemical, and mechano-chemical properties, which agree with experimental data.