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.

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.

Atomistic simulations of calcium aluminosilicate interfaced with liquid water

The dissolution behavior of calcium aluminosilicate based glass fibers, such as stone wool fibers, is an important consideration in mineral wool applications for both the longevity of the mineral wool products in humid environments and limiting the health impacts of released and inhaled fibers from the mineral wool product. Balancing these factors requires a molecular-level understanding of calcium aluminosilicate glass dissolution mechanisms, details that are challenging to resolve with experiment alone. Molecular dynamics simulations are a powerful tool capable of providing complementary atomistic insights regarding dissolution; however, they require force fields capable of describing not-only the calcium aluminosilicate surface structure but also the interactions relevant to dissolution phenomena. Here, a new force field capable of describing amorphous calcium aluminosilicate surfaces interfaced with liquid water is developed by fitting parameters to experimental and first principles simulation data of the relevant oxide-water interfaces, including ab initio molecular dynamics simulations performed for this work for the wüstite and periclase interfaces. Simulations of a calcium aluminosilicate surface interfaced with liquid water were used to test this new force field, suggesting moderate ingress of water into the porous glass interface. This design of the force field opens a new avenue for the further study of calcium and network-modifier dissolution phenomena in calcium aluminosilicate glasses and stone wool fibers at liquid water interfaces.

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 LaMnO₃ 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 Mn^IV+ 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.

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.

Revisiting the Conformational Isomerism of Dihaloethanes: A Hybrid Computational and Experimental Laboratory for the Undergraduate Curriculum

The conformational isomerism of disubstituted ethanes is a well-known concept that is part of every chemistry curriculum. Due to the species' simplicity, studying the (free) energy difference between the gauche and anti isomers has been the testing ground of experimental and computational techniques, such as Raman and IR spectroscopy, quantum chemistry, and atomistic simulations. While students normally receive formal training in spectroscopic techniques during their early undergraduate years, computational methods often receive less attention. In this work, we revisit the conformational isomerism of 1,2-dichloroethane and 1,2-dibromoethane and design a hybrid computational and experimental laboratory for our undergraduate chemistry curriculum with a focus on introducing computational techniques as a complementary research tool to experimentation. We show how commonly available Raman spectrometers and atomistic simulations performed on desktop computers can be combined to study the conformational isomerism of disubstituted ethanes while discussing the advantages and limitations of the different approaches.

Obtaining Consistent Free Energies for Ion Binding at Surfaces from Solution: Pathways versus Alchemy for Determining Kink Site Stability

Ion incorporation or removal from a solid at the interface with solution is a fundamental part of crystal growth. Despite this, there have been few quantitative determinations of the thermodynamics for such processes from atomistic molecular dynamics due to the associated technical challenges. In this study, we compute the free energies for ion removal from kink sites at the interface between NaCl and water as an illustrative example. To examine the influence of the free energy technique used, we compare methods that follow an explicit pathway for dissolution with those that focus on the thermodynamics of the initial and final states using metadynamics and free energy perturbation, respectively. While the initial results of the two approaches are found to be completely different, it is demonstrated that the thermodynamics can be reconciled with appropriate corrections for the standard states, thus illustrating the need for caution in interpreting raw free energy curves for ion binding as widely found in the literature. In addition, a new efficient approach is introduced to correct for the system size dependence of kink site energies both due to the periodic interaction of charges in an inhomogeneous dielectric system and due to the dipolar interactions between pairs of kinks along a row. Ultimately, it is shown that with suitable care, both classes of free energy techniques are capable of producing kink site stabilities that are consistent with the solubility of the underlying bulk solid. However, the precise values for individual kink sites exhibit a small systematic offset, which can be ascribed to the contribution of the interfacial potential to the pathway-based results. For the case of NaCl, the free energies of the kink sites relative to a 1 M aqueous solution for Na⁺ and Cl^- are found to be surprisingly different and of opposite sign, despite the ions having very similar hydration free energies.

Calculations of the Evolution of the Ca L23 Fine Structure in Amorphous Calcium Carbonate

Amorphous calcium carbonate (ACC) has been found in many different organisms. Biogenic ACC is frequently a precursor in the formation of calcite and aragonite. The process of structural transformation is therefore of great interest in the study of crystallization pathways in biomineralization. Changes in the prepeak/main peak (L₂'/L₂) intensity ratio of the Ca L₂₃-edge X-ray absorption spectroscopy (XAS) of Ca-rich particles in skeleton-building cells of sea urchin larva revealed that ACC precipitates through a continuum of states rather than through abrupt phase transitions involving two distinct phases as formerly believed. Using an atomic multiplet code, we show that only a tetragonal or "umbrella-like" distortion of the Ca coordination polyhedron can give rise to the observed continuum of states. We also show on the basis of the structures obtained from previous molecular dynamics simulations of hydrated nanoparticles that the Ca L₂₃-edge is not sensitive to atomic arrangements in the early stages of the transformation process.

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.

A Universal Force Field for Materials, Periodic GFN-FF: Implementation and Examination

In this study, the adaption of the recently published molecular GFN-FF for periodic boundary conditions (pGFN-FF) is described through the use of neighbor lists combined with appropriate charge sums to handle any dimensionality from 1D polymers to 2D surfaces and 3D solids. Numerical integration over the Brillouin zone for the calculation of π bond orders of periodic fragments is also included. Aside from adapting the GFN-FF method to handle periodicity, improvements to the method are proposed in regard to the calculation of topological charges through the inclusion of a screened Coulomb term that leads to more physical charges and avoids a number of pathological cases. Short-range damping of three-body dispersion is also included to avoid collapse of some structures. Analytic second derivatives are also formulated with respect to both Cartesian and strain variables, including prescreening of terms to accelerate the dispersion/coordination number contribution to the Hessian. The modified pGFN-FF scheme is then applied to a wide range of different materials in order to examine how well this universal model performs.

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

Calcite is the most stable polymorph of calcium carbonate (CaCO₃) 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.

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 HCO₃⁻ ions in CaCO₃ nucleation. Incorporation of HCO₃⁻ 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 HCO₃⁻ 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 HCO₃⁻ ions as active soluble species and structural constituents in CaCO₃ formation are of fundamental importance and provide a basis for a better understanding of physiological and geological mineralization.

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/m² 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.

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.

Ion Pairing and Multiple Ion Binding in Calcium Carbonate Solutions Based on a Polarizable AMOEBA Force Field and Ab Initio Molecular Dynamics

The speciation of calcium carbonate in water is important to the geochemistry of the world's oceans and has ignited significant debate regarding the mechanism by which nucleation occurs. Here, it is vital to be able to quantify the thermodynamics of ion pairing versus higher order association processes in order to distinguish between possible pathways. Given that it is experimentally challenging to quantify such species, here we determine the thermodynamics for ion pairing and multiple binding of calcium carbonate species using bias-enhanced molecular dynamics. In order to examine the uncertainties underlying these results, we derived a new polarizable force field for both calcium carbonate and bicarbonate in water based on the AMOEBA model to compare against our earlier rigid ion model, both of which are further benchmarked against ab initio molecular dynamics for the ion pair. Both force fields consistently indicate that the association of calcium carbonate ion pairs to form larger species is stable, though with an equilibrium constant that is lower than for ion pairing itself.

Wetting Properties of the CO2-Water-Calcite System via Molecular Simulations: Shape and Size Effects

Assessment of the risks and environmental impacts of carbon geosequestration requires knowledge about the wetting behavior of mineral surfaces in the presence of CO₂ and the pore fluids. In this context, the interfacial tension (IFT) between CO₂ and the aqueous fluid and the contact angle, θ, with the pore mineral surfaces are the two key parameters that control the capillary pressure in the pores of the candidate host rock. Knowledge of these two parameters and their dependence on the local conditions of pressure, temperature, and salinity is essential for the correct prediction of structural and residual trapping. We have performed classical molecular dynamics simulations to predict the CO₂-water IFT and the CO₂-water-calcite contact angle. The IFT results are consistent with previous simulations, where simple point charge water models have been shown to underestimate the water surface tension, thus affecting the simulated IFT values. When combined with the EPM2 CO₂ model, the SPC/Fw water model indeed underestimates the IFT in the low-pressure region at all temperatures studied. On the other hand, at high pressure and low temperature, the IFT is overestimated by ∼5 mN/m. Literature data regarding the CO₂/water/calcite contact angle on calcite are contradictory. Using our new set of force field parameters, we performed NVT simulations at 323 K and 20 MPa to calculate the contact angle of a water droplet on the calcite {10.4} surface in a CO₂ atmosphere. We performed simulations for both spherical and cylindrical droplet configurations for different initial radii to study the size dependence of the water contact angle on calcite in the presence of CO₂. Our results suggest that the contact angle of a cylindrical droplet, is independent of droplet size, for droplets with a radius of 50 Å or more. On the contrary, spherical droplets make a contact angle that is strongly influenced by their size. At the largest size explored in this study, both spherical and cylindrical droplets converge to the same contact angle, 38°, indicating that calcite is strongly wetted by water.

Simulation of Calcium Phosphate Prenucleation Clusters in Aqueous Solution: Association beyond Ion Pairing

Classical molecular dynamics simulations and free energy methods have been used to obtain a better understanding of the molecular processes occurring prior to the first nucleation event for calcium phosphate biominerals. The association constants for the formation of negatively charged complexes containing calcium and phosphate ions in aqueous solution have been computed, and these results suggest that the previously proposed calcium phosphate building unit, [Ca(HPO4)3]4-, should only be present in small amounts under normal experimental conditions. However, the presence of an activation barrier for the removal of an HPO4 2- ion from this complex indicates that this species could be kinetically trapped. Aggregation pathways involving CaHPO4, [Ca(HPO4)2]2-, and [Ca(HPO4)3]4- complexes have been explored with the finding that dimerization is favorable up to a Ca/HPO4 ratio of 1:2.

Homochirality in biomineral suprastructures induced by assembly of single-enantiomer amino acids from a nonracemic mixture

Since Pasteur first successfully separated right-handed and left-handed tartrate crystals in 1848, the understanding of how homochirality is achieved from enantiomeric mixtures has long been incomplete. Here, we report on a chirality dominance effect where organized, three-dimensional homochiral suprastructures of the biomineral calcium carbonate (vaterite) can be induced from a mixed nonracemic amino acid system. Right-handed (counterclockwise) homochiral vaterite helicoids are induced when the amino acid L-Asp is in the majority, whereas left-handed (clockwise) homochiral morphology is induced when D-Asp is in the majority. Unexpectedly, the Asp that incorporates into the homochiral vaterite helicoids maintains the same enantiomer ratio as that of the initial growth solution, thus showing chirality transfer without chirality amplification. Changes in the degree of chirality of the vaterite helicoids are postulated to result from the extent of majority enantiomer assembly on the mineral surface. These mechanistic insights potentially have major implications for high-level advanced materials synthesis.

Anomalous heat transport in binary hard-sphere gases

Equilibrium and nonequilibrium molecular dynamics (MD) are used to investigate the thermal conductivity of binary hard-sphere fluids. It is found that the thermal conductivity of a mixture can not only lie outside the series and parallel bounds set by their pure component values, but can lie beyond even the pure component fluid values. The MD simulations verify that revised Enskog theory can accurately predict nonequilibrium thermal conductivities at low densities and this theory is applied to explore the model parameter space. Only certain mass and size ratios are found to exhibit conductivity enhancements above the parallel bounds and dehancement below the series bounds. The anomalous dehancement is experimentally accessible in helium-hydrogen gas mixtures and a review of the literature confirms the existence of mixture thermal conductivity below the series bound and even below the pure fluid values, in accordance with the predictions of revised Enskog theory. The results reported here may reignite the debate in the nanofluid literature on the possible existence of anomalous thermal conductivities outside the series and parallel bounds as this Rapid Communication demonstrates they are a fundamental feature of even simple fluids.

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 (CaCO₃ ) 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 Ca²⁺ only adsorbs after CO₃^2- . 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.

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.

Structure, Energetics, and Dynamics of Screw Dislocations in Even n-Alkane Crystals

Spiral hillocks on n-alkane crystal surfaces were observed immediately after Frank recognized the importance of screw dislocations for crystal growth, yet their structures and energies in molecular crystals remain ill-defined. To illustrate the structural chemistry of screw dislocations that are responsible for plasticity in organic crystals and upon which the organic electronics and pharmaceutical industries depend, molecular dynamics was used to examine heterochiral dislocation pairs with Burgers vectors along [001] in n-hexane, n-octane, and n-decane crystals. The cores were anisotropic and elongated in the (110) slip plane, with significant local changes in molecular position, orientation, conformation, and energy. This detailed atomic level picture produced a distribution of strain consistent with linear elastic theory, giving confidence in the simulations. Dislocations with doubled Burgers vectors split into pairs with elementary displacements. These results suggest a pathway to understanding the mechanical properties and failure associated with elastic and plastic deformation in soft crystals.

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.

Negatively charged Ir(iii) cyclometalated complexes containing a chelating bis-tetrazolato ligand: synthesis, photophysics and the study of reactivity with electrophiles

The synthesis, the reactivity toward electrophiles and use for Ir(iii) based soft salts of new anionic Ir(iii) complexes containing a bis-tetrazolato ligand are described herein.

Photophysical and photochemical studies of tricarbonyl rhenium(i) N-heterocyclic carbene complexes containing azide and triazolate ligands

Rhenium NHC complexes bound to azide anions readily react with alkynes to form the corresponding triazolate complexes, a new class of photochemically active species.

The Use of Anisotropic Potentials in Modeling Water and Free Energies of Hydration

We propose a novel, anisotropic rigid-body intermolecular potential model to predict the properties of water and the hydration free energies of neutral organic solutes. The electrostatic interactions of water and the solutes are modeled using atomic multipole moments up to hexadecapole; these are obtained from distributed multipole analysis of the quantum mechanically computed charge densities and include average polarization effects in solution. The repulsion-dispersion water-water interactions are modeled with a three-site, exp-6 model fitted to the experimental liquid water density and oxygen-oxygen radial distribution function at ambient conditions. The proposed water model reproduces well several water properties not used in its parametrization, including vapor-liquid coexistence densities, the maximum in liquid water density at atmospheric pressure, the structure of ordered ice polymorphs, and the liquid water heat capacity. The model is used to compute the hydration free energy of 10 neutral organic solutes using explicit-solvent free energy perturbation. The solute-solute repulsion-dispersion intermolecular potential is obtained from previous parametrizations on organic crystal structures. In order to calculate the free energies of hydration, water-solute repulsion-dispersion interactions are modeled using Lorenz-Berthelot combining rules. The root-mean-square error of the predicted hydration free energies is 1.5 kcal mol(-1), which is comparable to the error found using a continuum mean-field quantum mechanical approach parametrized using experimental free energy of hydration data. The results are also contrasted with explicit-solvent hydration free energies obtained with an atomic charge representation of the solute's charge density computed at the same level of theory used to compute the distributed multipoles. Replacing the multipole description of the solute's charge density with an atomic charge model changes the free energy of hydration by as much as 3 kcal mol(-1) and provides an estimate for the effect of the modeling quality of the intermolecular electrostatic forces in free energy of solvation calculations.

Rhenium tetrazolato complexes coordinated to thioalkyl-functionalised phenanthroline ligands: synthesis, photophysical characterisation, and incubation in live HeLa cells

Three new complexes of formulation fac-[Re(CO)3(diim)L], where diim is either 1,10-phenanthroline or 1,10-phenanthroline functionalised at position 5 by a thioalkyl chain, and L is either a chloro or aryltetrazolato ancillary ligand, were synthesised and photophysically characterised. The complexes exhibit phosphorescent emission with maxima around 600 nm, originating from triplet metal-to-ligand charge transfer states with partially mixed ligand-to-ligand charge transfer character. The emission is relatively long-lived, within the 200-400 ns range, and with quantum yields of 2-4%. The complexes were trialed as cellular markers in live HeLa cells, along with two previously reported rhenium tetrazolato complexes bound to unsubstituted 1,10-phenanthroline. All five complexes exhibit good cellular uptake and non-specific perinuclear localisation. Upon excitation at 405 nm, the emission from the rhenium complexes could be clearly distinguished from autofluorescence, as demonstrated by spectral detection within the live cells. Four of the complexes did not appear to be toxic, however prolonged excitation could result in membrane blebbing. No major sign of photobleaching was detected upon multiple imaging on the same cell sample.

Microscopic evidence for liquid-liquid separation in supersaturated CaCO3 solutions

Recent experimental observations of the onset of calcium carbonate (CaCO3) mineralization suggest the emergence of a population of clusters that are stable rather than unstable as predicted by classical nucleation theory. This study uses molecular dynamics simulations to probe the structure, dynamics, and energetics of hydrated CaCO3 clusters and lattice gas simulations to explore the behavior of cluster populations before nucleation. Our results predict formation of a dense liquid phase through liquid-liquid separation within the concentration range in which clusters are observed. Coalescence and solidification of nanoscale droplets results in formation of a solid phase, the structure of which is consistent with amorphous CaCO3. The presence of a liquid-liquid binodal enables a diverse set of experimental observations to be reconciled within the context of established phase-separation mechanisms.

Enhanced deep-blue emission from Pt(II) complexes bound to 2-pyridyltetrazolate and an ortho-xylene-linked bis(NHC)cyclophane

The coordination of 2-pyridyltetrazolate and ortho-xylene-linked bis(NHC)cyclophane to Pt(II) yielded a novel complex characterised by enhanced pure deep-blue emission, whose intensity can be modulated via methylation of the tetrazole ring.

Self-assembling of zinc phthalocyanines on ZnO (1010) surface through multiple time scales

We adopt a hierarchic combination of theoretical methods to study the assembling of zinc phthalocyanines (ZnPcs) on a ZnO (1010) surface through multiple time scales. Atomistic simulations, such as model potential molecular dynamics and metadynamics, are used to study the energetics and short time evolution (up to ∼100 ns) of small ZnPc aggregates. The stability and the lifetime of large clusters is then studied by means of an atomistically informed coarse-grained model using classical molecular dynamics. Finally, the macroscopic time scale clustering phenomenon is studied by Metropolis Monte Carlo algorithms as a function of temperature and surface coverage. We provide evidence that at room temperature the aggregation is likely to occur at sufficiently high coverage, and we characterize the nature, morphology, and lifetime of ZnPc's clusters. We identify the molecular stripes oriented along [010] crystallographic directions as the most energetically stable aggregates.