Evolution of Surface Tension and Hansen Parameters of a Homologous Series of Imidazolium-Based Ionic Liquids.

In this study, we systematically analyze the surface tension and Hansen solubility parameters (HSPs) of imidazolium-based ionic liquids (ILs) with different anions ([NTf2]-, [PF6]-, [I]-, and [Br]-). These anions are combined with the classical 1-alkyl-3-methyl-substituted imidazolium cations ([CnC1Im]+) and a group of oligoether-functionalized imidazolium cations ([(mPEGn)2Im]+) based on methylated polyethylene glycol (mPEGn). In detail, the influences of the length of the alkyl- and the mPEGn-chain, the anion size, and the water content are investigated experimentally. For [CnC1Im]+-based ILs, the surface tension decreases with increasing alkyl chain length in all cases, but the magnitude of this decrease depends on the size of the anion ([NTf2]- < [PF6]- < [Br]- ≤ [I]-). Molecular dynamics (MD) simulations on [CnC1Im]+-based ILs indicate that these differences are caused by the interplay of charged and uncharged domains, in particular in the different anions, which affects the ability of the alkyl chains of the cation to orient toward the liquid-gas interface. An increase in the mPEGn-chain length of the [(mPEGn)2Im][A] ILs does not significantly influence the surface tension. These changes upon variation of the cation/anion combination do not correlate with the evolution of the HSPs for the two sets of ILs. Finally, our data suggest that significant water contents up to water mole fractions of x(H2O) = 0.25 do not significantly affect the surface tension of the studied binary IL-water mixtures.

Reversible chirality inversion of an AuAgx-cysteine coordination polymer by pH change.

Responsive chiral systems have attracted considerable attention, given their potential for diverse applications in biology, optoelectronics, photonics, and related fields. Here we show the reversible chirality inversion of an AuAgx-cysteine (AuAgx-cys) coordination polymer (CP) by pH changes. The polymer can be obtained by mixing HAuCl4 and AgNO3 with L-cysteine (or D-cysteine) in appropriate proportions in H2O (or other surfactant solutions). Circular dichroism (CD) spectrum is used to record the strong optical activity of the AuAg0.06-L-cys enantiomer (denoted as L0.06), which can be switched to that of the corresponding D0.06 enantiomer by alkalization (final dispersion pH > 13) and can be switched back after neutralization (final dispersion pH <8). Multiple structural changes at different pH values (≈9.6, ≈13) are observed through UV-Vis and CD spectral measurements, as well as other controlled experiments. Exploration of the CP synthesis kinetics suggests that the covalent bond formation is rapid and then the conformation of the CP materials would continuously evolve. The reaction stoichiometry investigation shows that the formation of CP materials with chirality inversion behavior requires the balancing between different coordination and polymerization processes. This study provides insights into the potential of inorganic stereochemistry in developing promising functional materials.

The small molecule activator S3969 stimulates the epithelial sodium channel by interacting with a specific binding pocket in the channel's ß-subunit.

The epithelial sodium channel (ENaC) is essential for mediating sodium absorption in several epithelia. Its impaired function leads to severe disorders, including pseudohypoaldosteronism type 1 and respiratory distress. Therefore, pharmacological ENaC activators have potential therapeutic implications. Previously, a small molecule ENaC activator (S3969) was developed. So far, little is known about molecular mechanisms involved in S3969-mediated ENaC stimulation. Here, we identified an S3969-binding site in human ENaC by combining structure-based simulations with molecular biological methods and electrophysiological measurements of ENaC heterologously expressed in Xenopus laevis oocytes. We confirmed a previous observation that the extracellular loop of ß-ENaC is essential for ENaC stimulation by S3969. Molecular dynamics simulations predicted critical residues in the thumb domain of ß-ENaC (Arg388, Phe391, and Tyr406) that coordinate S3969 within a binding site localized at the ß-γ-subunit interface. Importantly, mutating each of these residues reduced (R388H; R388A) or nearly abolished (F391G; Y406A) the S3969-mediated ENaC activation. Molecular dynamics simulations also suggested that S3969-mediated ENaC stimulation involved a movement of the α5 helix of the thumb domain of ß-ENaC away from the palm domain of γ-ENaC. Consistent with this, the introduction of two cysteine residues (ßR437C - γS298C) to form a disulfide bridge connecting these two domains prevented ENaC stimulation by S3969 unless the disulfide bond was reduced by DTT. Finally, we demonstrated that S3969 stimulated ENaC endogenously expressed in cultured human airway epithelial cells (H441). These new findings may lead to novel (patho-)physiological and therapeutic concepts for disorders associated with altered ENaC function.

Organic Nucleation: Water Rearrangement Reveals the Pathway of Ibuprofen.

The organic nucleation of the pharmaceutical ibuprofen is investigated, as triggered by the protonation of ibuprofen sodium salt at elevated pH. The growth and aggregation of nanoscale solution species by Analytical Ultracentrifugation and Molecular Dynamics (MD) simulations is tracked. Both approaches reveal solvated molecules, oligomers, and prenucleation clusters, their size as well as their hydration at different reaction stages. By combining surface-specific vibrational spectroscopy and MD simulations, water interacting with ibuprofen at the air-water interface during nucleation is probed. The results show the structure of water changes upon ibuprofen protonation in response to the charge neutralization. Remarkably, the water structure continues to evolve despite the saturation of protonated ibuprofen at the hydrophobic interface. This further water rearrangement is associated with the formation of larger aggregates of ibuprofen molecules at a late prenucleation stage. The nucleation of ibuprofen involves ibuprofen protonation and their hydrophobic assembly. The results highlight that these processes are accompanied by substantial water reorganization. The critical role of water is possibly relevant for organic nucleation in aqueous environments in general.

Bottom-to-top modeling of epoxy resins: From atomic models to mesoscale fracture mechanisms.

We outline a coarse-grained model of epoxy resins (bisphenol-F-diglycidyl-ether/3,5-diethyltoluene-2,4-diamine) to describe elastic and plastic deformation, cavitation, and fracture at the µm scale. For this, molecular scale simulation data collected from quantum and molecular mechanics studies are coarsened into an effective interaction potential featuring a single type of beads that mimic 100 nm scale building blocks of the material. Our model allows bridging the time-length scale problem toward experimental tensile testing, thus effectively reproducing the deformation and fracture characteristics observed for strain rates of 10-1 to 10-5 s-1. This paves the way to analyzing viscoelastic deformation, plastic behavior, and yielding characteristics by means of "post-atomistic" simulation models that retain the molecular mechanics of the underlying epoxy resin at length scales of 0.1-10 µm.

Molecular Characterization of Mesoporous Silica (Un)loading by Gemcitabine and Ibuprofen - An Interplay of Salt-Bridges and Hydrogen Bonds.

The molecular mechanisms of mesoporous silica nanomaterial (MSN) loading by gemcitabine and ibuprofen molecules, respectively, are elucidated as functions of pore geometry. Based on a small series of MSN archetypes, we use molecular dynamics simulations to systematically explore molecule-by-molecule loading of the carrier material. Apart from predicting the maximum active pharmaceutical ingredient (API) loading capacity, more detailed statistical analysis of the incorporation energy reveals dedicated profiles stemming from the interplay of guest-MSN salt-bridges/hydrogen bonding in concave and convex domains of the silica surfaces - which outcompete interactions among the drug molecules. Only after full coverage of the silica surface, we find secondary layer growth stabilized by guest-guest interactions exclusively. Based on molecular models, we thus outline a two-step type profile for drug release from MSN networks. Subject to the MSN structure, we find 50-75 % of the API within amorphous domains in the inner regions of the pores - from which drug release is provided at constant dissociation energy. In turn, the remaining 50-25 % of drug molecules are drastically hindered from dissociation.

Self-Assembled Supported Ionic Liquids.

Separation and reuse of the catalytically active metal complexes are persistent issues in homogeneous catalysis. Supported Ionic Liquid Phase (SILP) catalysts, where the catalytic center is dissolved in a thin film of a stable ionic liquid, deposited on a solid support, present a promising alternative. However, the dissolution of the metal center in the film leaves little control over its position and its activity. We present here four novel, task-specific ionic liquids [FPhn ImH R]I (n=1, 2; R=PEG2 , C12 H25 ), designed to self-assemble on a silica surface without any covalent bonding and offering a metal binding site in a controlled distance to the support. Advanced multinuclear solid-state NMR spectroscopic techniques under Magic Angle Spinning, complemented by molecular dynamics (MD) simulations, allow us to determine their molecular conformation when deposited inside SBA-15 as a model silica support. We provide here conceptual proof for a rational design of ionic liquids self-assembling into thin films, opening an avenue for a second, improved generation of SILP catalysts.

Efficient Assessment of 'Instantaneous pK' Values from Molecular Dynamics Simulations.

We present a molecular simulation approach to studying the role of local and momentary molecular environment for potential acid-base reactions. For this, we combine thermodynamic considerations on the pK of ionic species with rapid sampling of energy changes related to (de)protonation. Using dispersed carbonate ions in water as a reference, our approach aims at the fast assessment of the momentary protonation energy, and thus the 'instantaneous pK', of calcium-carbonate ion aggregates. The latter include transient complexes that are elusive to long sampling runs. This motivated the elaboration of approximate, yet particularly fast assessable sampling strategies. Along this line, we were able to characterize instantaneous pK values at a statistical accuracy of 0.4 pK units within sampling runs of only 10 ps duration, whereas statistical errors reduce to 0.1 pK units in 75 ps sampling runs, respectively. This readily enabled the required time resolution for the characterization of [Cax (CO3 )y ]2(x-y) aggregates with x=1,2 and y=1,2,3, respectively. In turn, the analysis of the pH-dependent nature of calcite-water interfaces and dynamically ordered liquid-like oxyanion polymers (dollop) domains is outlined at 10 ps resolution.

Interfaces in reinforced epoxy resins: from molecular scale understanding towards mechanical properties.

CONTEXT: We report on atomic level of detail analyses of polymer composite models featuring epoxy resin interfaces to silica, iron oxide, and cellulose layers. Using "reactive" molecular dynamics simulations to explore epoxy network formation, resin hardening is investigated in an unprejudiced manner. This allows the detailed characterization of salt-bridges and hydrogen bonds at the interfaces. Moreover, our sandwich-type composite systems are subjected to tensile testing along the interface normal. To elucidate the role of relaxation processes, we contrast (i) direct dissociation of the epoxy-metal oxide/cellulose contact layer, (ii) constant strain-rate molecular dynamics studies featuring (visco-)elastic deformation and bond rupture of the epoxy resin, and (iii) extrapolated relaxation dynamics mimicking quasi-static conditions. While the fracture mechanism is clearly identified as interface dissociation of the composite constituents, we still find damaging of the nearby polymer phase. The observed plastic deformation and local cavitation are rationalized from the comparably large stress required for the dissociation of salt-bridges, hydrogen bonds, and van der Waals contacts. Indeed, the delamination of the contact layers of epoxy resins with slabs of silica, magnetite, and cellulose call for a maximum stress of 33, 26, and 21 MPa, respectively, as compared to 84 MPa required for bulk epoxy yielding. METHODS: Molecular dynamics simulations using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code were augmented by a Monte Carlo-type procedure to probe epoxy bond formation (Macromolecules 53(22): 9698-9705). The underlying interaction models are split into conventional Generalized Amber Force Fields (GAFF) for non-reacting moieties and a recently developed reactive molecular mechanics potential enabling epoxy bond formation and cleavage (ACS Polymers Au 1(3): 165-174).

Structure Formation in an Ionic Liquid Wetting Layer: A Combined STM, IRAS, DFT and MD Study of [C2 C1 Im][OTf] on Au(111).

In a solid catalyst with ionic liquid layer (SCILL), ionic liquid (IL) coatings are used to improve the selectivity of noble metal catalysts. To understand the origins of this selectivity control, we performed model studies by surface science methods in ultrahigh vacuum (UHV). We investigated the growth and thermal stability of ultrathin IL films by infrared reflection absorption spectroscopy (IRAS). We combined these experiments with scanning tunneling microscopy (STM) to obtain information on the orientation of the ions, the interactions with the surface, the intermolecular interactions, and the structure formation. Additionally, we performed DFT calculations and molecular dynamics (MD) simulations to interpret the experimental data. We studied the IL 1-ethyl-3-methylimidazolium trifluoromethanesulfonate [C2 C1 Im][OTf] on Au(111) surfaces. We observe a weakly bound multilayer of [C2 C1 Im][OTf], which is stable up to 390 K, while the monolayer desorbs at ∼450 K. [C2 C1 Im][OTf] preferentially adsorbs at the step edges and elbows of the herringbone reconstruction of Au(111). The anion adsorbs via the SO3 group with the molecular axis perpendicular to the surface. At low coverage, the [C2 C1 Im][OTf] crystallizes in a glass-like 2D phase with short-range order. At higher coverage, we observe a phase transition to a 6-membered ring structure with long-range order.

Molecular Simulations and Network Analyses of Surface/Interface Effects in Epoxy Resins: How Bonding Adapts to Boundary Conditions.

In this study, we unravel the atomic structure of a covalent resin near boundaries such as surfaces and composite constituents. For this, a molecular simulation analysis of epoxy resin hardening under various boundary conditions was performed. On the atomic level of detail, molecular dynamics simulations were employed to study crosslinking reactions and self-organization of the polymer network within nm scale slab models. The resulting structures were then coarsened into a graph theoretical description for connectivity analysis of the nodes and combined with characterization of the node-to-node vector orientation. On this basis, we show that the local bonding of epoxy resins near interfaces tends to avoid under-coordinated linker sites. For both epoxy-vacuum surface models and epoxy-silica/epoxy cellulose interfaces, we find almost fully cured polymer networks. These feature a local increase in network linking lateral to the surface/interface, rather than the dangling of unreacted epoxy groups. Consequently, interface tension is low (as compared to the work of separating bulk epoxy), and the reactivity of the resin surface appears negligible.

Multi-Scale Modelling of Plastic Deformation, Damage and Relaxation in Epoxy Resins.

Epoxy resin plasticity and damage was studied from molecular dynamic simulations and interpreted by the help of constitutive modelling. For the latter, we suggested a physically motivated approach that aims at interpolating two well-defined limiting cases; namely, pulling at the vanishing strain rate and very rapid deformation; here, taken as 50% of the speed of sound of the material. In turn, to consider 0.1-10-m/s-scale deformation rates, we employed a simple relaxation model featuring exponential stress decay with a relaxation time of 1.5 ns. As benchmarks, deformation and strain reversal runs were performed by molecular dynamic simulations using two different strain rates. Our analyses show the importance of molecular rearrangements within the epoxy network loops for rationalizing the strain-rate dependence of plasticity and residual stress upon strain reversal. To this end, our constitutive model reasonably reproduced experimental data of elastic and visco-elastic epoxy deformation, along with the maximum stress experienced before fracturing. Moreover, we show the importance of introducing damage elements for mimicking the mechanical behavior of epoxy resins.

Small-Molecular-Weight Additives Modulate Calcification by Interacting with Prenucleation Clusters on the Molecular Level.

Small-molecular-weight (MW) additives can strongly impact amorphous calcium carbonate (ACC), playing an elusive role in biogenic, geologic, and industrial calcification. Here, we present molecular mechanisms by which these additives regulate stability and composition of both CaCO3 solutions and solid ACC. Potent antiscalants inhibit ACC precipitation by interacting with prenucleation clusters (PNCs); they specifically trigger and integrate into PNCs or feed PNC growth actively. Only PNC-interacting additives are traceable in ACC, considerably stabilizing it against crystallization. The selective incorporation of potent additives in PNCs is a reliable chemical label that provides conclusive chemical evidence that ACC is a molecular PNC-derived precipitate. Our results reveal additive-cluster interactions beyond established mechanistic conceptions. They reassess the role of small-MW molecules in crystallization and biomineralization while breaking grounds for new sustainable antiscalants.

Directed Dehydration as Synthetic Tool for Generation of a New Na4 SnS4 Polymorph: Crystal Structure, Na+ Conductivity, and Influence of Sb-Substitution.

We present the convenient synthesis and characterization of the new ternary thiostannate Na4 SnS4 (space group I 4 1 / a c d ) by directed removal of crystal water molecules from Na4 SnS4 ⋅14 H2 O. The compound represents a new kinetically stable polymorph of Na4 SnS4 , which is transformed into the known, thermodynamically stable form (space group P 4 ‾ 2 1 c ) at elevated temperatures. Thermal co-decomposition of mixtures with Na3 SbS4 ⋅9 H2 O generates solid solution products Na4-x Sn1-x Sbx S4 (x=0.01, 0.10) isostructural to the new polymorph (x=0). Incorporation of Sb5+ affects the bonding and local structural situation noticeably evidenced by X-ray diffraction, 119 Sn and 23 Na NMR, and 119 Sn Mössbauer spectroscopy. Electrochemical impedance spectroscopy demonstrates an enormous improvement of the ionic conductivity with increasing Sb content for the solid solution (σ25°C =2×10-3 , 2×10-2 , and 0.1 mS cm-1 for x=0, 0.01, and 0.10), being several orders of magnitude higher than for the known Na4 SnS4 polymorph.

Molecular Simulation Analyses of Polymorphism Control Factors by the Example of Carbamazepine Forms I-IV: A Blueprint for Industrial Drug Formulation?

We outline comparably simple molecular simulation techniques to elucidate the interactions that determine the polymorphism of carbamazepine. Starting from the established GAFF molecular mechanics model, only a small series of tailor-made improvements is needed to tackle the subtle differences in the interaction energies of polymorphs I - IV. On this basis, molecular dynamics simulations provide melting enthalpies at < 1 kcal/mol accuracy (0.2 kcal/mol for forms I-III) as compared to the experiment. Yet, the predicted stability ranking of III > I > II > IV only partially reproduces the experimentally observed III > I > IV > II series. Despite this limitation, we demonstrate how insights from molecular simulation offer the elucidation of possible factors for polymorph control. Apart from characterizing bulk crystals, we outline the evaluation of size-dependent profiles of crystallite formation energy. Contrasting the contributions of bulk, surface and edge terms to the formation energy of nano-scale precipitates, we suggest a multi-step nucleation mechanism leading from amorphous aggregates to crystallites. We argue that carbamazepine aggregates of less than ∼100 molecules adopt a spherical shape to minimize edge/surface energy - overcompensating the loss in bulk energy inherent to non-crystalline ordering in the inner core. In turn, for large crystallites polymorph form III is preferred, whilst suitable spatial confinement to crystallites of 100-500 carbamazepine molecules appears to promote form II.

Molecular Dynamics Simulations of Nitrate/MgO Interfaces and Understanding Metastability of Thermochemical Materials.

We explore the role of molten nitrate interfaces on MgO surface treatment for improving the reversibility of thermochemical energy storage via sorption and desorption of water or CO2. Our molecular dynamics simulations focus on melts of LiNO3, NaNO3, KNO3, and the triple eutectic mixture Li0.38Na0.18K0.44NO3 on the surface of MgO to provide atomic scale details of adsorbed layers and to rationalize interface energies. On this basis, a thermodynamic model is elaborated to characterize the effect of nitrate melts on the dehydration of Mg(OH)2 and to quantitatively explain the difference in dehydration temperatures of intact and LiNO3-doped Mg(OH)2.

On the Role of Amides and Imides for Understanding GaN Syntheses from Ammonia Solution: Molecular Mechanics Models of Ammonia, Amide and Imide Interactions with Gallium Nitride.

A key requisite to characterizing GaN precipitation from ammonia solution from molecular simulations is the availability of reliable molecular mechanics models for the interactions of gallium ions with NH3 , NH2- , and NH2- species, respectively. Here, we present a tailor-made force field which is fully compatible to an earlier developed GaN model, thus bridging the analyses of Ga3+ ions in ammonia solution with the aggregation of [Gax (NH)y (NH2 )z ]+3x-2y-z precursors and the modelling of GaN crystals. For this, quantum mechanical characterization of a series of Ga-coordination clusters is used for parameterization and benchmarking the generalized amber force field (GAFF2) and tailor-made refinements needed to achieve good agreement of both structural features and formation energy, respectively. The perspectives of our models for larger scale molecular dynamics simulations are demonstrated by the analyses of amide and imide defects arrangement during the growth of GaN crystal faces.

Molecular dynamics simulation study of NH4+ and NH2- in liquid ammonia: interaction potentials, structural and dynamical properties.

We provide tailor-made GAFF2-type interaction potentials for modeling ammonium and amide ions in ammonia. Based on harmonic approximation of intra-molecular bond stretching and bending, our force fields nicely reproduce the vibrational modes of NH4+ and NH2-, respectively. Moreover, quantum calculations of pair-wise NH4+/NH2--NH3 interactions were used for inter-molecular force field parameterization, while (NH3)n, [(NH4)(NH3)n]+, and [(NH2)(NH3)n]- complexes with n > 2, respectively, were reserved for benchmarking in terms of both structure and formation energy. Despite the limited reliability of molecular mechanics models for describing dimer complexes (n = 1), we find that GAFF2 reasonably reproduces [(NH4)(NH3)n]+ species for n = 2-4. For the assessment of [(NH2)(NH3)n]- complexes with n = 2-4, we however suggest the introduction of specific van der Waals parameters for amide-ammonia interactions. The application of the (extended) GAFF2 models is demonstrated for the study of ammonium and amide solvation in liquid ammonia at 240 K and 1 atm, respectively. On this basis, we suggest the applicability of our model for both gas phase and liquid states of ammonia.

The nano- and meso-scale structure of amorphous calcium carbonate.

Understanding the underlying processes of biomineralization is crucial to a range of disciplines allowing us to quantify the effects of climate change on marine organisms, decipher the details of paleoclimate records and advance the development of biomimetic materials. Many biological minerals form via intermediate amorphous phases, which are hard to characterize due to their transient nature and a lack of long-range order. Here, using Monte Carlo simulations constrained by X-ray and neutron scattering data together with model building, we demonstrate a method for determining the structure of these intermediates with a study of amorphous calcium carbonate (ACC) which is a precursor in the bio-formation of crystalline calcium carbonates. We find that ACC consists of highly ordered anhydrous nano-domains of approx. 2 nm that can be described as nanocrystalline. These nano-domains are held together by an interstitial net-like matrix of water molecules which generate, on the mesoscale, a heterogeneous and gel-like structure of ACC. We probed the structural stability and dynamics of our model on the nanosecond timescale by molecular dynamics simulations. These simulations revealed a gel-like and glassy nature of ACC due to the water molecules and carbonate ions in the interstitial matrix featuring pronounced orientational and translational flexibility. This allows for viscous mobility with diffusion constants four to five orders of magnitude lower than those observed in solutions. Small and ultra-small angle neutron scattering indicates a hierarchically-ordered organization of ACC across length scales that allow us, based on our nano-domain model, to build a comprehensive picture of ACC formation by cluster assembly from solution. This contribution provides a new atomic-scale understanding of ACC and provides a framework for the general exploration of biomineralization and biomimetic processes.

Deciphering the molecular mechanism of water boiling at heterogeneous interfaces.

Water boiling control evolution of natural geothermal systems is widely exploited in industrial processes due to the unique non-linear thermophysical behavior. Even though the properties of water both in the liquid and gas state have been extensively studied experimentally and by numerical simulations, there is still a fundamental knowledge gap in understanding the mechanism of the heterogeneous nucleate boiling controlling evaporation and condensation. In this study, the molecular mechanism of bubble nucleation at the hydrophilic and hydrophobic solid-water interface was determined by performing unbiased molecular dynamics simulations using the transition path sampling scheme. Analyzing the liquid to vapor transition path, the initiation of small void cavities (vapor bubbles nuclei) and their subsequent merging mechanism, leading to successively growing vacuum domains (vapor phase), has been elucidated. The molecular mechanism and the boiling nucleation sites' location are strongly dependent on the solid surface hydrophobicity and hydrophilicity. Then simulations reveal the impact of the surface functionality on the adsorbed thin water molecules film structuring and the location of high probability nucleation sites. Our findings provide molecular-scale insights into the computational aided design of new novel materials for more efficient heat removal and rationalizing the damage mechanisms.