Effects of amino-acid functionalization and pH value on temperature-dependent water dynamics in silica confinement

2H nuclear magnetic resonance (NMR) field-cycling relaxometry and broadband dielectric spectroscopy (BDS) studies show that water dynamics in silica pores with similar diameters (∼6 nm) strongly depend on the functionalization of the inner surfaces. In all studied confinements, we observe two prominent changes in the temperature dependence of water reorientation. Specifically, the activation energy of Ea ∼ 0.3 eV in the fully liquid state more than triples to Ea ∼ 1.0 eV upon partial crystallization at Tm ∼ 258 K. Furthermore, in the partially crystallized state, the liquid fraction shows a dynamical crossover at ∼185 K, where the common low-temperature behavior of confined water with Ea = 0.4-0.5 eV is established. However, the correlation times of water reorientation are up to two orders of magnitude longer in amino-acid functionalized silica pores than in pristine ones. Comparing the results for different functional groups, NMR and BDS consistently show that the slowdown is strongest for basic lysine followed by neutral alanine and, finally, acidic glutamic acid. Based on this order, one may speculate that the changed dynamics are a consequence of different pH values of water in confinements with different functional groups. Although pH measurements confirm that the pH value strongly depends on the amino-acid functionalization, this speculation must be rejected due to the observation that water with very different pH values does not show diverse reorientation dynamics when enclosed in identical pores.

Structures and Dynamics of Complex Guest Molecules in Confinement, Revealed by Solid-State NMR, Molecular Dynamics, and Calorimetry

This review gives an overview of current trends in the investigation of confined molecules such as water, small and higher alcohols, carbonic acids, ethylene glycol, and non-ionic surfactants, such as polyethylene glycol or Triton-X, as guest molecules in neat and functionalized mesoporous silica materials employing solid-state NMR spectroscopy, supported by calorimetry and molecular dynamics simulations. The combination of steric interactions, hydrogen bonds, and hydrophobic and hydrophilic interactions results in a fascinating phase behavior in the confinement. Combining solid-state NMR and relaxometry, DNP hyperpolarization, molecular dynamics simulations, and general physicochemical techniques, it is possible to monitor these confined molecules and gain deep insights into this phase behavior and the underlying molecular arrangements. In many cases, the competition between hydrogen bonding and electrostatic interactions between polar and non-polar moieties of the guests and the host leads to the formation of ordered structures, despite the cramped surroundings inside the pores.

Dipolar Dispersion Forces in Water-Methanol Mixtures: Enhancement of Water Interactions upon Dilution Drives Self-Association

Mixtures of short-chain alcohols and water produce anomalous thermodynamic and structural quantities, including molecular segregation into water-rich and alcohol-rich components. Herein, we used molecular dynamics simulations with polarizable models to investigate interactions that could drive the self-association of water molecules in mixtures with methanol (MeOH). As water was diluted with MeOH, significant changes in the distribution of molecules and solvation properties occurred, where water exhibited a clear preference for self-association. When common structural quantities were analyzed, it was found that there was a clear reduction in water-water hydrogen bonding and tetrahedral order (both in terms of typical bulk behavior), contrary to the observed water self-association. However, when dipolar dispersion forces between all molecules as a function of system composition were analyzed, it was found that water-water dipolar interactions became significantly stronger with dilution (6-fold stronger interaction in 75% MeOH compared to 0% MeOH). This was only observed for water, where MeOH-MeOH interactions became weaker as the systems were more dilute in MeOH. These forces result from specific dipole orientations, likely occurring to adopt lower energy configurations (i.e., head-to-tail or antiparallel). For water, this may result from lost other interactions (e.g., hydrogen bonding), leading to more rotational freedom between the dipole moments. These intriguing changes in dipolar interactions, which directly result from structural changes, can therefore explain, in part, the driving force for water self-association in MeOH-water mixtures.

A New Study on the Structure, and Phase Transition Temperature of Bulk Silicate Materials by Simulation Method of Molecular Dynamics

In this paper, the structure and phase transition temperature of bulk silicate materials are studied by the simulation method (SM) of molecular dynamics (MD). In this research, all samples are prepared on the same nanoscale material model with the atomic number of 3000 atoms, for which the SM of MD is performed with Beest-Kramer-van Santen and van Santen pair interaction potentials under cyclic boundary conditions. The obtained results show that both the model size (l) and the total energy of the system (Etot) increase slowly in the low temperature (T) region (negative T values) at pressure (P), P = 0 GPa. However, the increase of l determines the Etot value with very large values in the high T region. It is found that l decreases greatly in the high T region with increasing P, and vice versa. In addition, when P increases, the decrease in the Etot value is small in the low T region, but large in the high T region. As a consequence, a change appears in the lengths of the Si-Si, Si-O, and O-O bonds, which are very large in the high T and high P regions, but insignificant in the low T and low P regions. Furthermore, the structural unit number of SiO7 appears at T > 2974 K in the high P region. The obtained results will serve as the basis for future experimental studies to exploit the stored energy used in semiconductor devices.

Structural and Dynamical Properties of Liquids in Confinements: A Review of Molecular Dynamics Simulation Studies

Molecular dynamics (MD) simulations are a powerful tool for detailed studies of altered properties of liquids in confinement, in particular, of changed structures and dynamics. They allow, on one hand, for perfect control and systematic variation of the geometries and interactions inherent in confinement situations and, on the other hand, for type-selective and position-resolved analyses of a huge variety of structural and dynamical parameters. Here, we review MD simulation studies on various types of liquids and confinements. The main focus is confined aqueous systems, but also ionic liquids and polymer and silica melts are discussed. Results for confinements featuring different interactions, sizes, shapes, and rigidity will be presented. Special attention will be given to situations in which the confined liquid and the confining matrix consist of the same type of particles and, hence, disparate liquid-matrix interactions are absent. Findings for the magnitude and the range of wall effects on molecular positions and orientations and on molecular dynamics, including vibrational motion and structural relaxation, are reviewed. Moreover, their dependence on the parameters of the confinement and their relevance to theoretical approaches to the glass transition are addressed.

NMR Relaxometry Accessing the Relaxation Spectrum in Molecular Glass Formers

It is a longstanding question whether universality or specificity characterize the molecular dynamics underlying the glass transition of liquids. In particular, there is an ongoing debate to what degree the shape of dynamical susceptibilities is common to various molecular glass formers. Traditionally, results from dielectric spectroscopy and light scattering have dominated the discussion. Here, we show that nuclear magnetic resonance (NMR), primarily field-cycling relaxometry, has evolved into a valuable method, which provides access to both translational and rotational motions, depending on the probe nucleus. A comparison of ¹H NMR results indicates that translation is more retarded with respect to rotation for liquids with fully established hydrogen-bond networks; however, the effect is not related to the slow Debye process of, for example, monohydroxy alcohols. As for the reorientation dynamics, the NMR susceptibilities of the structural (α) relaxation usually resemble those of light scattering, while the dielectric spectra of especially polar liquids have a different broadening, likely due to contributions from cross correlations between different molecules. Moreover, NMR relaxometry confirms that the excess wing on the high-frequency flank of the α-process is a generic relaxation feature of liquids approaching the glass transition. However, the relevance of this feature generally differs between various methods, possibly because of their different sensitivities to small-amplitude motions. As a major advantage, NMR is isotope specific; hence, it enables selective studies on a particular molecular entity or a particular component of a liquid mixture. Exploiting these possibilities, we show that the characteristic Cole-Davidson shape of the α-relaxation is retained in various ionic liquids and salt solutions, but the width parameter may differ for the components. In contrast, the low-frequency flank of the α-relaxation can be notably broadened for liquids in nanoscopic confinements. This effect also occurs in liquid mixtures with a prominent dynamical disparity in their components.

Confinement effects on glass-forming mixtures: Insights from a combined experimental approach to aqueous ethylene glycol solutions in silica pores

We perform nuclear magnetic resonance, broadband dielectric spectroscopy, and differential scanning calorimetry studies to ascertain the dynamical behaviors of aqueous ethylene glycol (EG) solutions in silica pores over broad temperature ranges. Both translational and rotational motions are analyzed, and the pore diameter (2.4–9.2 nm) and the EG concentration (12–57 mol. %) are varied, leading to fully liquid or partially crystalline systems. It is found that the translational diffusion coefficient strongly decreases when the diameter is reduced, resulting in a slowdown of nearly three orders of magnitude in the narrowest pores, while the confinement effects on the rotational correlation times are moderate. For the fully liquid solutions, we attribute bulk-like and slowed down reorientation processes to the central and interfacial pore regions, respectively. This coexistence is found in all the studied pores, and, hence, the range of the wall effects on the solution dynamics does not exceed ∼1 nm. Compared to the situation in the bulk, the concentration dependence is reduced in confinements, implying that the specific interactions of the molecular species with the silica walls lead to preferential adsorption. On the other hand, bulk-like structural relaxation is not observed in the partially frozen samples, where the liquid is sandwiched between the silica walls and the ice crystallites. Under such circumstances, there is another relaxation process with a weaker temperature dependence, which is observed in various kinds of partially frozen aqueous systems and denoted as the x process.

Effects of partial crystallization on the glassy slowdown of aqueous ethylene glycol solutions

Combining differential scanning calorimetry, nuclear magnetic resonance, and broadband dielectric spectroscopy studies, we ascertain the glass transition of aqueous ethylene glycol (EG) solutions, in particular the effects of partial crystallization on their glassy slowdown. For the completely liquid solutions in the weakly supercooled regime, it is found that the dynamics of the components occur on very similar time scales, rotational and translational motions are coupled, and the structural (α) relaxation monotonously slows down with increasing EG concentration. Upon cooling, partial crystallization strongly alters the glassy dynamics of EG-poor solutions; in particular, it strongly retards the α relaxation of the remaining liquid fraction, causing a non-monotonous concentration dependence, and it results in a crossover from non-Arrhenius to Arrhenius temperature dependence. In the deeply supercooled regime, a recrossing of the respective α-relaxation times results from the Arrhenius behaviors of the partially frozen EG-poor solutions together with the non-Arrhenius behavior of the fully liquid EG-rich solutions. Exploiting the isotope selectivity of nuclear magnetic resonance, we observe different rotational dynamics of the components in this low-temperature range and determine the respective contributions to the ν relaxation decoupling from the α relaxation when the glass transition is approached. The results suggest that the ν process, which is usually regarded as a water process, actually also involves the EG molecules. In addition, we show that various kinds of partially crystalline aqueous systems share a common relaxation process, which is associated with the frozen fraction and differs from that of bulk hexagonal ice.

Structure and dynamics of water confined in cylindrical nanopores with varying hydrophobicity

We report a detailed study of the main structural and dynamical features of water confined in model Lennard-Jones nanopores with tunable hydrophobicity and finite length ([Formula: see text] Å). The generic model of cylindrical confinement used is able to reproduce the wetting features of a large class of technologically and biologically relevant systems spanning from crystalline nanoporous materials, to mesoporous silica and ion channels. The aim of this work is to discuss the influence of parameters such as wall hydrophobicity, temperature, and pore size on the structural and dynamical features of confined water. Our simulation campaign confirmed the existence of a core domain in which water displays bulk-like structural features even in extreme ([Formula: see text] Å) confinement, while dynamical properties were shown to depend non-trivially on the size and hydrophobicity of the pores. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

2H NMR study on temperature-dependent water dynamics in amino-acid functionalized silica nanopores

We prepare various amino-acid functionalized silica pores with diameters of ∼6 nm and study the temperature-dependent reorientation dynamics of water in these confinements. Specifically, we link basic Lys, neutral Ala, and acidic Glu to the inner surfaces and combine ²H nuclear magnetic resonance spin-lattice relaxation and line shape analyses to disentangle the rotational motions of the surfaces groups and the crystalline and liquid water fractions coexisting below partial freezing. Unlike the crystalline phase, the liquid phase shows reorientation dynamics, which strongly depends on the chemistry of the inner surfaces. The water reorientation is slowest for the Lys functionalization, followed by Ala and Glu and, finally, the native silica pores. In total, the rotational correlation times of water at the different surfaces vary by about two orders of magnitude, where this span is largely independent of the temperature in the range ∼200-250 K.

MD simulations of charged binary mixtures reveal a generic relation between high- and low-temperature behavior

Experimental studies of the glassy slowdown in molecular liquids indicate that the high-temperature activation energy E_∞ of glass-forming liquids is directly related to their glass transition temperature T_g. To further investigate such a possible relation between high- and low-temperature dynamics in glass-forming liquids, we analyze the glassy dynamics of binary mixtures using molecular dynamics simulations. We consider a binary mixture of charged Lennard-Jones particles and vary the partial charges of the particles and, thus, the high-temperature activation energy and the glass transition temperature of the system. Based on previous results, we introduce a phenomenological model describing relaxation times over the whole temperature regime from high temperatures to temperatures well inside the supercooled regime. By investigating the dynamics of both particle species on molecular and diffusive length scales along isochoric and isobaric pathways, we find a quadratic charge dependence of both E_∞ and T_g, resulting in an approximately constant ratio of both quantities independent of the underlying observable, the thermodynamic ensemble, and the particle species, and this result is robust against the actual definition of T_g. This generic relation between the activation energy and the glass transition temperature indicates that high-temperature dynamics and the glassy slowdown are related phenomena, and the knowledge of E_∞ may allow us to approximately predict T_g.

Confinement Effects on Glass-Forming Aqueous Dimethyl Sulfoxide Solutions

Combining broadband dielectric spectroscopy and nuclear magnetic resonance studies, we analyze the reorientation dynamics and the translational diffusion associated with the glassy slowdown of the eutectic aqueous dimethyl sulfoxide solution in nano-sized confinements, explicitly, in silica pores with different diameters and in ficoll and lysozyme matrices at different concentrations. We observe that both rotational and diffusive dynamics are slower and more heterogeneous in the confinements than in the bulk but the degree of these effects depends on the properties of the confinement and differs for the components of the solution. For the hard and the soft matrices, the slowdown and the heterogeneity become more prominent when the size of the confinement is reduced. In addition, the dynamics are more retarded for dimethyl sulfoxide than for water, implying specific guest-host interactions. Moreover, we find that the temperature dependence of the reorientation dynamics and of the translational diffusion differs in severe confinements, indicating a breakdown of the Stokes-Einstein-Debye relation. It is discussed to what extent these confinement effects can be rationalized in the framework of core-shell models, which assume bulk-like and slowed-down motions in central and interfacial confinement regions, respectively.

Small Molecules, Non-Covalent Interactions, and Confinement

This review gives an overview of current trends in the investigation of small guest molecules, confined in neat and functionalized mesoporous silica materials by a combination of solid-state NMR and relaxometry with other physico-chemical techniques. The reported guest molecules are water, small alcohols, and carbonic acids, small aromatic and heteroaromatic molecules, ionic liquids, and surfactants. They are taken as characteristic role-models, which are representatives for the typical classes of organic molecules. It is shown that this combination delivers unique insights into the structure, arrangement, dynamics, guest-host interactions, and the binding sites in these confined systems, and is probably the most powerful analytical technique to probe these systems.

Quasielastic neutron scattering studies on couplings of protein and water dynamics in hydrated elastin

Performing quasielastic neutron scattering measurements and analyzing both elastic and quasielasic contributions, we study protein and water dynamics of hydrated elastin. At low temperatures, hydration-independent methyl group rotation dominates the findings. It is characterized by a Gaussian distribution of activation energies centered at about E_m = 0.17 eV. At ∼195 K, coupled protein-water motion sets in. The hydration water shows diffusive motion, which is described by a Gaussian distribution of activation energies with E_m = 0.57 eV. This Arrhenius behavior of water diffusion is consistent with previous results for water reorientation, but at variance with a fragile-to-strong crossover at ∼225 K. The hydration-related elastin backbone motion is localized and can be attributed to the cage rattling motion. We speculate that its onset at ∼195 K is related to a secondary glass transition, which occurs when a β relaxation of the protein has a correlation time of τ_β ∼ 100 s. Moreover, we show that its temperature-dependent amplitude has a crossover at the regular glass transition T_g = 320 K of hydrated elastin, where the α relaxation of the protein obeys τ_α ∼ 100 s. By contrast, we do not observe a protein dynamical transition when water dynamics enters the experimental time window at ∼240 K.

The Influence of Surface Modified Silica Nanoparticles: Properties of Epoxy Nanocomposites

The influence of the surface-modified (CCS) and un-modified (UCS) silica nanoparticles on epoxy nanocomposites were studied. Two different nanocomposites systems were synthesized using tetraethylorthosilicate (TEOS) and 3-(triethoxysilyl) propylamine APTES as a precursor and coupling agent, respectively. In the uncoupled composite system (UCS) the silica particles were solely generated using TEOS as a precursor. The APTES was used as a coupling agent to chemically link the silica (SiO2) particles to the matrix in the coupled composite system (CCS). Both composite systems were fabricated as thin films. The SiO2 epoxy nanocomposites thin films were characterized by Differential Scanning Calorimetry (DSC), RAMAN, Fourier Transform Infra-Red (FT-IR), Scanning Electron Microscopy (SEM) and Thermal Gravimetric (TGA) analysis. RAMAN and FT-IR analysis confirmed the curing of epoxy resin and the generation of the inorganic structural network formation. SEM analysis of these nanocomposites revealed that silica particles were uniformly dispersed in the epoxy matrix. DSC analysis of the nano-composites revealed an increase in glass transition (T g) temperature with the addition of nanofiller. TGA analysis shows enhanced thermal stability of the coupled composite system in comparison to the neat and uncoupled epoxy composite system.

Static field gradient NMR studies of water diffusion in mesoporous silica

NMR diffusometry is used to ascertain the pore-size dependent water diffusion in MCM-41 and SBA-15 silica over broad temperature ranges. Detailed analysis of 1H and 2H NMR stimulated-echo decays reveals that fast water motion through voids between different silica particles impairs such studies in the general case. However, water diffusion inside single pores is probed in the present approach, which applies high static field gradients to enhance the spatial resolution of the experiment and uses excess water in combination with subzero temperatures to embed the silica particles in an ice matrix and, thus, to suppress interparticle water motion. It is found that the diffusion of confined water slows down by almost two orders of magnitude when the pore diameter is reduced from 5.4 nm to 2.1 nm at weak cooling. In the narrower silica pores, the temperature dependence of the self-diffusion coefficient of water is well described by an Arrhenius law with an activation energy of Ea = 0.40 eV. The Arrhenius behavior extends over a broad temperature range of at least 207-270 K, providing evidence against a fragile-to-strong crossover in response to a proposed liquid-liquid phase transition near 225 K. In the wider silica pores, partial crystallization results in a discontinuous temperature dependence. Explicitly, the diffusion coefficients drop when cooling through the pore-size dependent melting temperatures Tm of confined water. This finding can be rationalized by the fact that water can explore the whole pore volumes above Tm, but is restricted to narrow interfacial layers sandwiched between silica walls and ice crystallites below this temperature. Comparing our findings for water diffusion with previous results for water reorientation, we find significantly different temperature dependencies, indicating that the Stokes-Einstein-Debye relation is not obeyed.

Characterization of the non-Arrhenius behavior of supercooled liquids by modeling nonadditive stochastic systems

The characterization of the formation mechanisms of amorphous solids is a large avenue for research, since understanding its non-Arrhenius behavior is challenging to overcome. In this context, we present one path toward modeling the diffusive processes in supercooled liquids near glass transition through a class of nonhomogeneous continuity equations, providing a consistent theoretical basis for the physical interpretation of its non-Arrhenius behavior. More precisely, we obtain the generalized drag and diffusion coefficients that allow us to model a wide range of non-Arrhenius processes. This provides a reliable measurement of the degree of fragility of the system and an estimation of the fragile-to-strong transition in glass-forming liquids, as well as a generalized Stokes-Einstein equation, leading to a better understanding of the classical and quantum effects on the dynamics of nonadditive stochastic systems.

Effects of confinement on supercooled tetrahedral liquids

We use molecular dynamics simulations to ascertain the effects of geometrical restriction on glass-forming tetrahedral liquids. Striving for a broad approach, we study families of waterlike and silicalike liquids, for which we systematically scale the partial charges and, hence, the relevance of the tetrahedral networks. The confined liquids and the confining matrices consist of the same type of particles to avoid disruptive interactions and distorted structures at the interfaces. Spatially resolved analyses show that these neutral confinements still impose static mobility gradients and density correlations on the liquids. We quantify the increasing degree and range of the altered properties upon cooling. For both families of models, common relations describe the confinement effects of all systems with tetrahedral order, while deviations occur for systems with lower polarities and different structures. The observations are rationalized by considering the fact that a pinned wall imprints a static energy landscape to a neighboring liquid. We explore the properties of this landscape based on changes in vibrational motion and structural relaxation and find that typical barrier heights amount to two to three times the activation energy of bulk dynamics. Combining the present and previous results, we predict the evolution of confinement effects down to the glass transition temperature for liquids without fragile-to-strong crossover. In addition, it is found for silicalike liquids that the temperature dependence of dynamic and static correlation lengths from confinement studies is not affected when cooling through fragile-to-strong transitions of the bulk materials, casting doubt on the relevance of these length scales for the glassy slowdown.

Hydrogen Bond Interaction of Ascorbic Acid with Urea: Experimental and Theoretical Study

The interactions of ascorbic acid (AA) with urea were investigated by using the cyclic voltammetry, density functional theory, atoms in molecules and natural bond orbital analyses. The experimental and theoretical results show that the hydrogen bonds are formed between AA and urea, wherein the mainly interaction sites are the hydrogen atoms on enediol of AA and the oxygen atom on carbonyl of urea. The electrochemical behavior of AA was significantly affected by above interactions.

Properties of Hydrogen-Bonded Liquids at Interfaces

Effects of interfaces on hydrogen-bonded liquids play major roles in nature and technology. Despite their importance, a fundamental understanding of these effects is still lacking. In large parts, this shortcoming is due to the high complexity of these systems, leading to an interference of various interactions and effects. Therefore, it is advisable to take gradual approaches, which start from well designed and defined model systems and systematically increase the level of intricacy towards more complex mimetics. Moreover, it is necessary to combine insights from a multitude of methods, in particular, to link novel preparation strategies and comprehensive experimental characterization with inventive computational and theoretical modeling. Such concerted approach was taken by a group of preparative, experimentally, and theoretically working scientists in the framework of Research Unit FOR 1583 funded by the Deutsche Forschungsgemeinschaft (German Research Foundation). This special issue summarizes the outcome of this collaborative research. In this introductory article, we give an overview of the covered topics and the main results of the whole consortium. The following contributions are review articles or original works of individual research projects.

2H NMR Studies on the Dynamics of Pure and Mixed Hydrogen-Bonded Liquids in Confinement

2H NMR is used to ascertain dynamical behaviors of pure and mixed hydrogen-bonded liquids in bulk and in confinement. Detailed comparisons of previous and new results in broad dynamic and temperature ranges reveal that confinement effects differ for various liquids and confinements. For water, molecular reorientation strongly depends on the confinement size, with much slower and less fragile structural relaxation under more severe geometrical restriction. Moreover, a dynamical crossover occurs when a fraction of solid water forms so that the dynamics of the fraction of liquid water becomes even more restricted and, as a consequence, changes from bulk-like to interface-dominated. For glycerol, by contrast, confinement has weak effects on the reorientation dynamics. Mixed hydrogen-bonded liquids show even more complex dynamical behaviors. For aqueous solutions, the temperature dependence of the structural relaxation becomes discontinuous when the concentration changes due to a freezing of water fractions. This tendency for partial crystallization is enhanced rather than reduced by confinement, because different liquid-matrix interactions of the molecular species induce micro-phase segregation, which facilitates ice formation in water-rich regions. In addition, dynamical couplings at solvent-protein interfaces are discussed. It is shown that, on the one hand, solvent dynamics are substantially slowed down at protein surfaces and, on the other hand, protein dynamics significantly depend on the composition and, thus, the viscosity of the solvent. Furthermore, a protein dynamical transition occurs when the amplitude of water-coupled restricted backbone dynamics vanishes upon cooling.