Self-intermediate scattering function analysis of supercooled water confined in hydrophilic silica nanopores

Role of anisotropy in understanding the molecular grounds for density scaling in dynamics of glass-forming liquids

Molecular Dynamics (MD) simulations of glass-forming liquids play a pivotal role in uncovering the molecular nature of the liquid vitrification process. In particular, much focus was given to elucidating the interplay between the character of intermolecular potential and molecular dynamics behaviour. This has been tried to achieve by simulating the spherical particles interacting via isotropic potential. However, when simulation and experimental data are analysed in the same way by using the density scaling approaches, serious inconsistency is revealed between them. Similar scaling exponent values are determined by analysing the relaxation times and pVT data obtained from computer simulations. In contrast, these values differ significantly when the same analysis is carried out in the case of experimental data. As discussed thoroughly herein, the coherence between results of simulation and experiment can be achieved if anisotropy of intermolecular interactions is introduced to MD simulations. In practice, it has been realized in two different ways: (1) by using the anisotropic potential of the Gay-Berne type or (2) by replacing the spherical particles with quasi-real polyatomic anisotropic molecules interacting through isotropic Lenard-Jones potential. In particular, the last strategy has the potential to be used to explore the relationship between molecular architecture and molecular dynamics behaviour. Finally, we hope that the results presented in this review will also encourage others to explore how 'anisotropy' affects remaining aspects related to liquid-glass transition, like heterogeneity, glass transition temperature, glass forming ability, etc.

Mesoporous Matrices as a Promising New Generation of Carriers for Multipolymorphic Active Pharmaceutical Ingredient Aripiprazole

The enhancement of the properties (i.e., poor solubility and low bioavailability) of currently available active pharmaceutical ingredients (APIs) is one of the major goals of modern pharmaceutical sciences. Among different strategies, a novel and innovative route to reach this milestone seems to be the application of nanotechnology, especially the incorporation of APIs into porous membranes composed of pores of nanometric size and made of nontoxic materials. Therefore, in this work, taking the antipsychotic API aripiprazole (APZ) infiltrated into various types of mesoporous matrices (anodic aluminum oxide, native, and silanized silica) characterized by similar pore diameters (d = 8-10 nm) as an example, we showed the advantage of incorporated systems in comparison to the bulk substance considering the crystallization kinetics, molecular dynamics, and physical stability. Calorimetric investigations supported by the temperature-dependent X-ray diffraction measurements revealed that in the bulk system the recrystallization of polymorph III, which next is converted to the mixture of forms IV and I, is visible, while in the case of confined samples polymorphic forms I and III of APZ are produced upon heating of the molten API with different rates. Importantly, the two-step crystallization observed in thermograms obtained for the API infiltrated into native silica templates may suggest crystal formation by the interfacial and core molecules. Furthermore, dielectric studies enabled us to conclude that there is no trace of crystallization of spatially restricted API during one month of storage at T = 298 K. Finally, we found that in contrast to the crystalline and amorphous bulk samples, all examined confined systems show a logarithmic increase in API dissolution over time (very close to a prolonged release effect) without any sign of precipitation. Our data demonstrated that mesoporous matrices appear to be interesting candidates as carriers for unstable amorphous APIs, like APZ. In addition to protecting them against crystallization, they can provide the desired prolonged release effect, which may increase the drug concentration in the blood (resulting in higher bioavailability). We believe that the "nanostructirization" in terms of the application of porous membranes as a novel generation of drug carriers might open unique perspectives in the further development of drugs characterized by prolonged release.

Liquid-Liquid Flow at Nanoscale: Slip and Hydrodynamic Boundary Conditions

Nonequilibrium molecular dynamics (NEMD) simulations have been performed to describe the flow of a fluid nanolayer confined by another fluid. The results show that the behavior of liquids can still be described by the Navier-Stokes equation at the nanoscale, i.e., when only few molecular layers are involved. NEMD furthermore gives additional knowledge on flow. Indeed, while a very small slip is evidenced for a solid-liquid interface as, e.g., in lubrication, the slip lengths are significantly larger at the liquid-liquid interface, as encountered, e.g., in droplet coalescence. The slip lengths of the two fluids are linked. The increase in hydrodynamic slip for liquid-liquid interfaces is attributed to the enhancement of fluid diffusion, which reduces friction.

Studying the Intermolecular Interactions, Structural Dynamics, and Non-Equilibrium Kinetics of Cilnidipine Infiltrated into Alumina and Silica Pores

In the present study, the behavior of the calcium channel blocker cilnidipine (CLN) infiltrated into silica (SiO₂) and anodic aluminum oxide (AAO) porous membranes characterized by a similar pore size (d = 8 nm and d = 10 nm, respectively) as well as the bulk sample has been investigated using differential scanning calorimetry, broadband dielectric spectroscopy (BDS), and Fourier-transform infrared spectroscopy (FTIR) techniques. The obtained data suggested the existence of two sets of CLN molecules in both confined systems (core and interfacial). They also revealed the lack of substantial differences in inter- and intramolecular dynamics of nanospatially restricted samples independently of the applied porous membranes. Moreover, the annealing experiments (isothermal time-dependent measurements) performed on the confined CLN clearly indicated that the whole equilibration process under confinement is governed by structural relaxation. It was also found that the β_anneal parameters obtained from BDS and FTIR data upon equilibration of both confined samples are comparable (within 10%) to each other, while the equilibration constants are significantly different. This finding strongly emphasizes that there is a close connection between the inter- and intramolecular dynamics under nanospatial restriction.

New Insights from Nonequilibrium Kinetics Studies on Highly Polar S-Methoxy-PC Infiltrated into Pores

Herein, the annealing of highly polar (S)-(-)-4-methoxymethyl-1,3-dioxolan-2-one (S-methoxy-PC) within alumina and silica porous membranes characterized by different pore diameters was studied by means of dielectric spectroscopy. We found a significant slowing down of the structural dynamics of confined S-methoxy-PC with annealing time below and, surprisingly, also above the glass transition temperatures of the interfacial layer, T_{g,interfacial}. Furthermore, unexpectedly, a change in the slope of temperature dependencies of the characteristic time scale of this process τ_anneal, at T_{g,interfacial} for all confined samples, was reported. By modeling τ_anneal(T), we noted that the observed enormous variation of τ_anneal results from a decrease of the pore radius due to the vitrification of the interfacial molecules. This indicates that the enhanced dynamics of confined materials upon cooling is mainly controlled by the interfacial molecules.

Dynamic Heterogeneity at the Interface of an Intrinsically Disordered Peptide

It is believed that water around an intrinsically disordered protein or peptide (IDP) in an aqueous environment plays an important role in guiding its conformational properties and aggregation behavior. However, despite its importance, only a handful of studies exploring the correlation between the conformational motions of an IDP and the microscopic properties of water at its surface are reported. Attempts have been made in this work to study the dynamic properties of water present in the vicinity of α-synuclein, an IDP associated with Parkinson's disease (PD). Room temperature molecular dynamics (MD) simulations of eight α-synuclein_1-95 peptides with a wide range of initial conformations have been carried out in aqueous media. The calculations revealed that due to solid-like caging motions, the translational and rotational mobility of water molecules near the surfaces of the peptide repeat unit segments R1 to R7 are significantly restricted. A small degree of dynamic heterogeneity in the hydration environment around the repeat units has been observed with water near the hydrophobic R6 unit exhibiting relatively more restricted diffusivity. The time scales involving the overall structural relaxations of peptide-water and water-water hydrogen bonds near the peptide have been found to be correlated with the time scale of diffusion of the interfacial water molecules. We believe that the relatively more hindered dynamic environment near R6 can help create water-mediated contacts centered around R6 between peptide monomers at a higher concentration, thereby enhancing the early stages of peptide aggregation.

Density Scaling of Translational and Rotational Molecular Dynamics in a Simple Ellipsoidal Model near the Glass Transition

In this paper, we show that a simple anisotropic model of supercooled liquid properly reflects some density scaling properties observed for experimental data, contrary to many previous results obtained from isotropic models. We employ a well-known Gay–Berne model earlier parametrized to achieve a supercooling and glass transition at zero pressure to find the point of glass transition and explore volumetric and dynamic properties in the supercooled liquid state at elevated pressure. We focus on dynamic scaling properties of the anisotropic model of supercooled liquid to gain a better insight into the grounds for the density scaling idea that bears hallmarks of universality, as follows from plenty of experimental data collected near the glass transition for different dynamic quantities. As a result, the most appropriate values of the scaling exponent γ are established as invariants for a given anisotropy aspect ratio to successfully scale both the translational and rotational relaxation times considered as single variable functions of density^γ/temperature. These scaling exponent values are determined based on the density scaling criterion and differ from those obtained in other ways, such as the virial–potential energy correlation and the equation of state derived from the effective short-range intermolecular potential, which is qualitatively in accordance with the results yielded from experimental data analyses. Our findings strongly suggest that there is a deep need to employ anisotropic models in the study of glass transition and supercooled liquids instead of the isotropic ones very commonly exploited in molecular dynamics simulations of supercooled liquids over the last decades.

Viscoelastic response of fluid trapped between two dissimilar van der Waals surfaces

Employing grand canonical Monte-Carlo and molecular dynamics simulations, the viscoelastic response of trapped fluid under molecularly thin confinement by walls having different wall-fluid interaction strengths, is investigated. With increase in slit asymmetry, given by the ratio of interaction strengths of the wall having strong wall-fluid interaction to that of the wall with weak wall-fluid interaction, a crossover in effective density of the fluid film, from rarer (R) to denser (D) than the bulk density is observed. Upon increasing asymmetry further, the dense fluid (F) layers undergo bond-orientational (S) ordering. The variation of viscoelastic relaxation time with scaled asymmetry shows a universal behavior, independent of slit width, with two distinct regimes. Below a critical value of asymmetry, the viscoelastic relaxation time is a slowly varying function of asymmetry, comparable with the structural relaxation time. Beyond the critical asymmetry, on the other hand, viscoelastic response time shows a sharp increase upon increasing asymmetry, deviating markedly from the structural relaxation time. Interestingly the critical asymmetry value is found to correlate with R to D crossover. The microscopic origin of the two-regime universal behavior of viscoelastic response time is found to stem from the fact that below critical asymmetry, the overall viscoelastic behaviour of the slit is dominated by that of the fast relaxing layer close to the weakly attracting surface, while above the critical asymmetry, the relaxation behaviour is guided by the dense fluid layer adjacent to the strongly attracting wall. In vicinity of fluid to ordering transition, the loss and storage moduli merge for low frequencies as in gel-like mechanical behaviour. The storage modulus takes over the loss modulus in the phase co-existence region even before the long ranged order sets in. Our findings bear important implications for fluid transport in hetero-structured geometry in nanotechnology.

Local structure and molecular dynamics of highly polar propylene carbonate derivative infiltrated within alumina and silica porous templates

Herein, we examined the effect of finite size and wettability on the structural dynamics and the molecular arrangement of the propylene carbonate derivative, (S)-(-)-4-methoxymethyl-1,3-dioxolan-2-one (assigned as s-methoxy-PC), incorporated into alumina and silica porous templates of pore diameters d = 4 nm-10 nm using Raman and broadband dielectric spectroscopy, differential scanning calorimetry, and x-ray diffraction. It was demonstrated that only subtle changes in the molecular organization and short-range order of confined s-methoxy-PC molecules were detected. Yet, a significant deviation of the structural dynamics and depression of the glass transition temperatures, T_g, was found for all confined samples with respect to the bulk material. Interestingly, these changes correlate with neither the finite size effects nor the interfacial energy but seem to vary with wettability, generally. Nevertheless, for s-methoxy-PC infiltrated into native (more hydrophilic) and modified (more hydrophobic) silica templates of the same nanochannel size (d = 4 nm), a change in the dynamics and T_g was negligible despite a significant variation in wettability. These results indicated that although wettability might be a suitable variable to predict alteration of the structural dynamics and depression of the glass transition temperature, other factors, i.e., surface roughness and the density packing, might also have a strong contribution to the observed confinement effects.

Are hydrogen supramolecular structures being suppressed upon nanoscale confinement? The case of monohydroxy alcohols

In this paper, the molecular dynamics, H-bonding pattern and wettability of the primary and secondary monohydroxyalcohols, 2-ethyl-1-hexanol (2E1H), 2-ethyl-1-butanol (2E1B) and 5-methyl-3-heptanol (5M3H) infiltrated into native and functionalized silica and alumina pores having pore diameters, d = 4 nm and d = 10 nm, have been studied with the use of Broadband Dielectric (BDS) and Fourier Transform InfraRed (FTIR) spectroscopies, as well as contact angle measurements. We found significant differences in the behavior of alcohols forming chain- (2E1H, 2E1B) or micelle-like (5M3H) supramolecular structures despite of their similarities in the wettability and interfacial energy. It turned out that nanoassociates as well as H-bonds are more or less affected by the confinement dependently on the chemical structure and alcohol order. Moreover, a peculiar behavior of the self-assemblies at the interface was noted in the latter material (5M3H). Finally, it was found that irrespectively to the sample, type of pores, functionalization, the temperature evolution of Debye relaxation times, τ_D, of the confined systems deviates from the bulk behavior always at similar τ_D due to vitrification of the interfacial layer. This finding is a clear indication that unexpectedly dynamics (mobility) of the supramolecular structures close to the hydrophilic and hydrophobic surfaces is similar in each system.

Tracer Diffusion in Tightly-Meshed Homogeneous Polymer Networks: A Brownian Dynamics Simulation Study

We report Brownian dynamics simulations of tracer diffusion in regularly crosslinked polymer networks in order to elucidate the transport of a tracer particle in polymer networks. The average mesh size of homogeneous polymer networks is varied by assuming different degrees of crosslinking or swelling, and the size of a tracer particle is comparable to the average mesh size. Simulation results show subdiffusion of a tracer particle at intermediate time scales and normal diffusion at long times. In particular, the duration of subdiffusion is significantly prolonged as the average mesh size decreases with increasing degree of crosslinking, for which long-time diffusion occurs via the hopping processes of a tracer particle after undergoing rattling motions within a cage of the network mesh for an extended period of time. On the other hand, the cage dynamics and hopping process are less pronounced as the mesh size decreases with increasing polymer volume fractions. The interpretation is provided in terms of fluctuations in network mesh size: at higher polymer volume fractions, the network fluctuations are large enough to allow for collective, structural changes of network meshes, so that a tracer particle can escape from the cage, whereas, at lower volume fractions, the fluctuations are so small that a tracer particle remains trapped within the cage for a significant period of time before making infrequent jumps out of the cage. This work suggests that fluctuation in mesh size, as well as average mesh size itself, plays an important role in determining the dynamics of molecules and nanoparticles that are embedded in tightly meshed polymer networks.

High-pressure experiments as a novel perspective to study the molecular dynamics of glass-forming materials confined at the nanoscale

Herein, we report the pioneering high-pressure dielectric studies on the dynamics of a model van der Waals glass-forming liquid bisphenol-A diglycidyl ether (DGEBA) infiltrated into anodic aluminum oxide (AAO) templates of the mean pore sizes, d = 150 and d = 18 nm. It was found that although the shape of the structural relaxation process varies with the confinement, it remains constant under varying thermodynamic conditions for a given pore diameter. Consequently, the time-temperature-pressure (TTP) rule satisfied for the majority of bulk substances is also obeyed for the spatially restricted liquid. We have also shown for the first time that there is a decoupling between the core and interfacial mobility at elevated pressure. Moreover, it was noted that the structural dynamics of the former fraction of molecules becomes systematically shorter with respect to the bulk DGEBA during the compression. The enhanced structural dynamics of the core material, as well as the varying pressure coefficients of the glass transition temperature of the interfacial and core molecules, have been discussed in the context of a distinct evolution in their free volume/density packing with respect to the bulk DGEBA, and a change in the interfacial tension, which may lead to the enhanced wettability of the liquid adsorbed onto the pore walls under different thermodynamic conditions. The performed high-pressure measurements offer novel perspectives to explore the combination of two different effects, compression and confinement, which might be a breakthrough in the study of the glass transition phenomenon and the behavior of soft materials confined at the nanoscale.