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.

The origin of the density scaling exponent for polyatomic molecules and the estimation of its value from the liquid structure.

In this article, we unravel the problem of interpreting the density scaling exponent for the polyatomic molecules representing the real van der Waals liquids. Our studies show that the density scaling exponent is a weighted average of the exponents of the repulsive terms of all interatomic interactions that occur between molecules, where the potential energy of a given interaction represents its weight. It implies that potential energy is a key quantity required to calculate the density scaling exponent value for real molecules. Finally, we use the well-known method for potential energy estimation and show that the density scaling exponent could be successfully predicted from the liquid structure for fair representatives of the real systems.

The impact of the length of alkyl chain on the behavior of benzyl alcohol homologues - the interplay between dispersive and hydrogen bond interactions.

In this work, we examined the effect of the length of alkyl chain attached to the benzene ring on the self-assembling phenomena for a series of phenyl alcohol (PhA) derivatives, from phenylmethanol (benzyl alcohol) to 7-phenyl-1-heptanol, by means of X-Ray Diffraction (XRD), Differential Scanning Calorimetry (DSC), Fourier Transform Infrared (FTIR) spectroscopy, and Broadband Dielectric Spectroscopy (BDS) methods. XRD data in the reciprocal and real spaces showed a gradual increase in the local order with the elongation of the alkyl chain. However, the position and full width at half maximum of the main diffraction peak exhibited a non-systematic behavior. To better understand this fact, PhAs were subjected to FTIR spectroscopic studies. These investigations revealed that the association degree and the activation energy of dissociation increase as the alkyl chain length grows. On the other hand, BDS data showed a non-monotonic variation in the Kirkwood correlation factor with increasing length of the alkyl chain, indicating a competition between interactions of the non-polar and polar parts of the molecules in the studied PhAs. Finally, it was also found that the molar surface entropy for PhAs increases with the number of methylene groups, approaching values reported for alkanes, which indicates suppression of the surface order for PhAs with a long alkyl chain. This variability of the various parameters as a function of the length of the side chain shows that the interplay between soft interactions has a strong impact on the local structure and intra and intermolecular dynamics of the studied PhAs.

Correlation between configurational entropy, excess entropy, and ion dynamics in imidazolium-based ionic liquids: Test of the Adam-Gibbs model.

The Adam-Gibbs (AG) model, linking thermodynamics with molecular dynamics of glass-forming liquids, plays a crucial role in the studies of the glass transition phenomenon. We employ this approach to investigate the relationship between ion dynamics and thermodynamics in three imidazolium-based ionic liquids in the current work. We show that the AG relation, -log10σdc ∝ (TSc)-1 (where σdc, T, and Sc denote the dc-conductivity, absolute temperature, and configurational entropy, respectively), does not work when the whole supercooled liquid state is considered. Meanwhile, a linear relationship between -log10σdc and (TSe)-1 (where Se denotes the excess entropy) was observed in the entire supercooled range. On the other hand, the generalized AG model log10σdc ∝ (TSc α)-1 with an additional free parameter α successfully describes the relation between σdc and Sc. The determined α values being less than unity indicate that the configurational entropy is insufficient to govern the ion dynamics. Meanwhile, we found a systematical decrease in α with the elongation of the alkyl chain attached to the imidazolium ring.

The relation between molecular dynamics and configurational entropy in room temperature ionic liquids: Test of Adam-Gibbs model.

In this communication, the Adam-Gibbs model connecting molecular dynamics with configurational entropy is tested for the first time for ionic liquids. For this purpose, we investigate simultaneously the shear viscosity η and configurational entropy Sc of an aprotic ionic liquid: 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (BMIm TFSI). Comparing the Sc data obtained by the combination of Vogel-Fulcher-Tammann and Adam-Gibbs equations to the Sc points determined directly from the calorimetric experiment, good agreement is found in the entire supercooled liquid region. These results indicate the validity of the Adam-Gibbs model in materials with electrostatic interactions being dominated. These important findings not only generalize the applications of the Adam-Gibbs theory but also provide an opportunity to gain insight into the relationship between thermodynamics and molecular dynamics in ionic liquids.

Discharge of the Nanopore Confinement Effect on the Glass Transition Dynamics via Viscous Flow.

Using dielectric spectroscopy, we demonstrate that confinement-induced changes in the glass transition dynamics, as observed for polymethylphenylsiloxane in alumina nanopores, reveal a pronounced nonequilibrium nature. Our results indicate that glass formers confined to nanopores are able to recover their bulklike mobility. We found that the characteristic time constant of such an equilibration process correlates with an extremely slow viscous flow rate in cylindrical channels of nanometer size. Thus, all the way to equilibrium, confinement effects seen in faster segmental dynamics are released through the viscous flow which eventually helps to eliminate surplus volume gained by nanoconstrained polymers upon cooling.

Experimental Evidence for a State-Point-Independent Density-Scaling Exponent in Ionic Liquids.

This Letter addresses a fundamental issue of condensed-matter physics, which is the validity of the density-scaling concept. For this purpose, the ambient and high-pressure conductivity measurements of two selected ionic liquids (ILs), with the different contribution of H-bonding interactions, were performed in the dynamic range of 13 orders of magnitude and corresponding to the density changes as large as 20%. All experimental data obtained within one compound are shown to superimpose each other when plotted as a function of ρ^{γ}/T. These results clearly show that for studied ILs the scaling exponent is a state-point-independent parameter that is in odds with the recent findings for van der Waals liquid [Sanz et al., Phys. Rev. Lett. 122, 055501 (2019)PRLTAO0031-900710.1103/PhysRevLett.122.055501].

Dramatic slowing down of the conformational equilibrium in the silyl derivative of glucose in the vicinity of the glass transition temperature.

The vitrification process is usually preceded by a significant change (around 6-8 decades) in the viscosity, structural relaxation times, or diffusion that occurs in a relatively small range of temperatures in fragile liquids. Along with this phenomenon, conformations of the molecules vary as well. In fact, this process is studied in bulk polymers and high molecular weight materials deposited in the form of thin films. On the other hand, spatial rearrangement of small glass formers in the supercooled liquid state has not been intensively investigated, so far. Herein, data obtained from measurements carried out using various experimental techniques on supercooled 1,2,3,4,6-penta-O-(trimethylsilyl)-d-glucopyranose (S-GLU) have revealed that rotations of silyl moieties along with the deformation in the saccharide ring are significantly slowed down in the vicinity of the glass transition temperature (Tg). These intramolecular reorganizations affect the structural relaxation time, atomic pair distribution function, integrated intensity, as well as a number of bands and signals observed, respectively, in the Raman and NMR spectra. Data reported herein offer a better understanding of the conformational variation and time scale of this process in the complex and flexible molecules around the Tg.

The effect of molecular architecture on the physical properties of supercooled liquids studied by MD simulations: Density scaling and its relation to the equation of state.

Theoretical concepts in condensed matter physics are typically verified and also developed by exploiting computer simulations mostly in simple models. Predictions based on these usually isotropic models are often at odds with measurement results obtained for real materials. One of the examples is an intriguing problem within the density scaling idea that has attracted attention in recent decades due to its hallmarks of universality, i.e., the fact that the difference between the density scaling exponent and the exponent of the equation of state is observed for real materials, whereas it has not been reported for the model system. In this paper, we use new model molecules of simple but anisotropic architecture to study the effect of molecular anisotropy on the dynamic and thermodynamic properties of the system. We identify the applicable range of intermolecular interactions for a given physical process, and then we explain the reason for observed differences between the behavior of the model and real systems. It demonstrates that the new model systems open broad perspectives for simulation and theoretical research, for example, into unifying concepts in the glass transition physics.

Formation of new polymorphs and control of crystallization in molecular glass-formers by electric field.

Static electric fields were observed to modify the crystallization behavior in a simple supercooled liquid, leading to a new crystal polymorph that could not be obtained in the absence of a field, even under high-pressure conditions. Using different thermal protocols and field amplitudes in the range from 40 to 200 kV cm-1, changes in both nucleation and crystal growth rates of 4-vinyl-propylene carbonate (vinyl-PC) are revealed. Remarkably, all field-induced changes in the crystallization behaviour were found to be fully reversible and do not affect dynamics of the tested liquid. Because vinyl-PC is a simple polar molecule, these field induced features are expected to occur in many other materials having permanent dipole moments. Our results highlight the important role of an external electric field as an additional control variable to influence the crystallization tendency of molecular glass-formers, and provide new opportunities in pharmaceutical science or organic electronics.

The effect of electrostatic interactions on the formation of pharmaceutical eutectics.

Over the past decade, the formation of pharmaceutical eutectics has become a very attractive strategy to increase the bioavailability of active pharmaceutical ingredients (APIs). A great advantage of a eutectic phase, which can be obtained by simple physical mixing of solid materials, is the possibility to obtain a material with desired physicochemical properties only by varying the molar ratio of the parent components. In this work, we have investigated the ability of two protic ionic liquids (PILs), which are hydrochloride salts of lidocaine and prilocaine, as well as their non-ionic counterparts, to form eutectic mixtures. To gain an insight into the calorimetric properties of the formed dipolar and ionic mixtures, differential scanning calorimetry was employed. The mechanism of formation of deep eutectic mixtures on the molecular level was investigated by ab initio quantum mechanics calculations. The effect of electrostatic interactions on the eutectic transition, glass forming ability and the physical stability of pharmaceutical eutectics was also revealed.

Studying molecular dynamics of the slow, structural, and secondary relaxation processes in series of substituted ibuprofens.

In this paper, the molecular dynamics of a series of ester derivatives of ibuprofen (IBU), in which the hydrogen atom from the hydroxyl group was substituted by the methyl, isopropyl, hexyl, and benzyl moieties, has been investigated using Broadband dielectric (BD), Nuclear magnetic resonance (NMR), and Raman spectroscopies. We found that except for benzyl IBU (Ben-IBU), an additional process (slow mode, SM) appears in dielectric spectra in all examined compounds. It is worth noting that this relaxation process was observed for the first time in non-modified IBU (a Debye relaxation). According to suggestions by Affouard and Correia [J. Phys. Chem. B. 114, 11397 (2010)] as well as further studies by Adrjanowicz et al. [J. Chem. Phys. 139, 111103 (2013)] on Met-IBU, it was attributed to synperiplanar-antiperiplanar conformational changes within the molecule. Herein, we have shown that with an increasing molecular weight of the substituent, the relaxation times of the SM become longer and its activation energy significantly increases. Moreover, this new relaxation mode was found to be broader than a simple Debye relaxation in Iso-IBU and Hex-IBU. Additional complementary NMR studies indicated that either there is a significant slowdown of the rotation around the O=C-O-R moiety or this kind of movement is completely suppressed in the case of Ben-IBU. Therefore, the SM is not observed in the dielectric loss spectra of this compound. Finally, we carried out isothermal experiments on the samples which have a different thermal history. Interestingly, it turned out that the relaxation times of the structural processes are slightly shorter with respect to those obtained from temperature dependent measurements. This effect was the most prominent in the case of Hex-IBU, while for Ben-IBU, it was not observed at all. Additional time-dependent measurements revealed the ongoing equilibration manifested by the continuous shift of the structural process, until it finally reached its equilibrium position. Further Raman investigations showed that this effect may be related to the rotational/conformational equilibration of the long hexyl chains. Our results are the first ones demonstrating that the structural process is sensitive to the conformational equilibration occurring in the specific highly viscous systems.

A New Method To Identify Physically Stable Concentration of Amorphous Solid Dispersions (I): Case of Flutamide + Kollidon VA64.

In this paper, a novel approach to determine stable concentration in API-polymer systems is presented. As a model, binary amorphous mixtures flutamide (FL) drug with a copolymer Kollidon VA64 (PVP/VA) have been used. It is worthwhile to note that finding an effective method to achieve this goal is a matter of great importance because physical stability of the amorphous pharmaceuticals is the key issue that is investigated worldwide. Due to the fact that molecular dynamics was found to be the crucial factor affecting physical stability of disordered pharmaceuticals, we examined it for both neat FL and its PVP/VA mixtures by means of broadband dielectric spectroscopy (BDS). Thorough investigation of the impact of polymeric additive on the molecular mobility of disordered FL reveals unusual, previously unreported behavior. Namely, simultaneously with the beginning of the recrystallization process, we observe some transformation from unstable supersaturated concentration of investigated mixture to the different, unknown concentration of FL-PVP/VA. Observed, during BDS experiment, transformation enables us to determine the limiting, highly physically stable concentration of FL in PVP/VA polymer (saturated solution), which is equivalent to FL + 41% wt. of PVP/VA. The described high physical stability of this unveiled system has been confirmed by means of long-term XRD measurements. According to our knowledge, this is the first time when such a behavior has been observed by means of BDS.

Molecular Factors Governing the Liquid and Glassy States Recrystallization of Celecoxib in Binary Mixtures with Excipients of Different Molecular Weights.

Transformation of poorly water-soluble crystalline pharmaceuticals to the amorphous form is one of the most promising strategies to improve their oral bioavailability. Unfortunately, the amorphous drugs are usually thermodynamically unstable and may quickly return to their crystalline form. A very promising way to enhance the physical stability of amorphous drugs is to prepare amorphous compositions of APIs with certain excipients which can be characterized by significantly different molecular weights, such as polymers, acetate saccharides, and other APIs. By using different experimental techniques (broadband dielectric spectroscopy, differential scanning calorimetry, X-ray diffraction) we compare the effect of adding the large molecular weight polymer-polyvinylpyrrolidone (PVP K30)-and the small molecular weight excipient-octaacetylmaltose (acMAL)-on molecular dynamics as well as the tendency to recrystallization of the amorphous celecoxib (CEL) in the amorphous solid dispersions: CEL-PVP and CEL-acMAL. The physical stability investigations of the binary systems were performed in both the supercooled liquid and glassy states. We found that acMAL is a better inhibitor of recrystallization of amorphous CEL than PVP K30 deep in the glassy state (T < Tg). In contrast, PVP K30 is a better crystallization inhibitor of CEL than acMAL in the supercooled liquid state (at T > Tg). We discuss molecular factors governing the recrystallization of amorphous CEL in examined solid dispersions.

Studying the Crystallization of Various Polymorphic Forms of Nifedipine from Binary Mixtures with the Use of Different Experimental Techniques.

In this paper the crystal growth of nifedipine from pure system and from binary mixtures composed of active substance (API) and two acetylated disaccharides, maltose and sucrose (NIF-acMAL, NIF-acSUC, 5:1 weight ratio), was investigated. Optical snapshots supported by X-ray diffraction (XRD) and Fourier transform infrared spectroscopy (FTIR) measurements showed that mainly ß and α forms of nifedipine grow up in all investigated samples. They also revealed that the morphology of growing crystals strongly depends on the presence of modified carbohydrates and temperature conditions. Interestingly, it was found that the activation barrier for the crystal growth of the ß polymorph is not affected by acetylated saccharides while the one estimated for the α form changes significantly from 48.5 kJ/mol (pure API) up to 122 kJ/mol (NIF-acMAL system). Moreover, the relationship between the crystal growth rate and structural relaxation times for pure NIF and solid dispersions were analyzed. It turned out that there is a clear decoupling between the crystal growth rate and structural dynamics in both NIF-acMAL and NIF-acSUC binary mixtures. This is in line with recent reports indicating the decoupling phenomenon to be a universal feature of soft matter in the close vicinity of the glass transition temperature.

Atorvastatin as a Promising Crystallization Inhibitor of Amorphous Probucol: Dielectric Studies at Ambient and Elevated Pressure.

The aim of this article was to check the physical stability of the amorphous form of probucol at both standard storage and manufacturing conditions. Our studies clearly show that disordered form of the examined, cholesterol lowering, agent stored at ambient pressure does not reveal any tendency toward recrystallization. The physical stability of neat probucol stored at ambient pressure has been investigated (i) at room temperature by means of X-ray diffraction technique (XRD) as well as (ii) at T = 333 K by means of broadband dielectric spectroscopy (BDS). Due to the fact that compression is an important stage of drugs manufacturing we additionally performed physical stability tests of amorphous probucol at elevated pressure. The recrystallization tendency of the examined pharmaceutical has been tracked online from the initial and further up to a few hours after compression by means of the high pressure BDS technique. These experiments indicate that even very small pressure applied during the sample compression immediately induce its recrystallization. Since, the sensitivity on pressure eliminates probucol from the group of physically stable amorphous APIs, its stabilization is required. Taking into account that there are many scientific reports describing the positive effect of coadministration of probucol with the drug atorvastatin, we used the latter as probucol's crystallization inhibitor.

Comparison of high pressure and nanoscale confinement effects on crystallization of the molecular glass-forming liquid, dimethyl phthalate.

High pressure and nanoscopic confinement are two different strategies commonly employed to modify the physicochemical properties of various materials. Both strategies act mostly by changing the molecular packing. In this work, we performed a comparative study on the effect of compression and confined geometry on crystallization of a molecular liquid. Dielectric spectroscopy was employed to investigate the crystallization of the van der Waals liquid, dimethyl phthalate, in nanoporous alumina of different pore sizes as well as on increased pressure (up to 200 MPa). The analysis of the crystallization kinetics under varying thermodynamic conditions revealed that both strategies affect the crystallization behavior of the sample in very distinct ways. Compression shifts the maximum crystallization rate towards a higher temperature and broadens it. As a result, it is more challenging to avoid crystallization upon cooling the liquid at high pressure. In contrast, when the same material is incorporated into nanopores, crystallization significantly slows down and the maximum rate shifts towards a lower temperature with decreasing pore size. Finally, we show that crystallization in nanoporous alumina is accompanied by pre-crystallization effects upon which a shift of the α-relaxation peak is observed. An equilibration process prior to the initiation of crystallization was detected for the confined material both above and below the glass transition temperature of the interfacial layer, while not in the bulk.

A study on the progress of mutarotation above and below the Tg and the relationship between constant rates and structural relaxation times.

Comprehensive FTIR studies on the progress of mutarotation in d-fructose mixed with maltitol have been carried out over a wide range of temperatures, both above and below the glass transition temperature Tg. In addition to the analysis of single bands, we have developed a completely new approach considering the full spectral range to follow the overall progress of the reaction. We have found that at the calorimetric Tg, there is a clear change in the temperature dependence of constant rates. The activation barrier for mutarotation changes from around 59 kJ mol-1 (the supercooled state) to around 249 kJ mol-1 (the glassy state). This dramatic variation in the activation barrier is consistent with the change in the mechanism of this specific chemical conversion, as theoretically considered by Wlodarczyk et al. [Phys. Chem. Chem. Phys., 2014, 16, 4694-4698]. Alternatively, it can also be connected to the change in the viscosity of the sample. Additionally, we investigated the relationship between constant rates (k) of mutarotation, structural relaxation times (τα), and dc conductivity (σdc) above and below the glass transition temperature. It was found that there was a linear correlation between all these quantities; they scale with various exponents changing at Tg. Our results also indicate that a single activation barrier might not be sufficient to describe the mutarotation process.

How is charge transport different in ionic liquids? The effect of high pressure.

Modern ionic liquids (ILs) are considered green solvents for the future applications due to their inherited advantages and remarkable transport properties. One of the ubiquitous properties of ILs is their intrinsic ionic conductivity. However, understanding of the super-Arrhenius behavior of the ionic conductivity process at elevated pressure still remains elusive and crucial in glass science. In this work, we investigate the ion transport properties of 1-butyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide: [C4mim][NTf2], 1-butylimidazolium bis[(trifluoromethyl)-sulfonyl]imide: [C4Him][NTf2] and 1-butylimidazolium hydrogen sulfate: [C4Him][HSO4] ILs in the supercooled liquid state using dielectric spectroscopy at ambient and high pressure. We present the experimental data in the dynamic window of the conductivity formalism to examine the charge transport properties. The frequency-dependent ionic conductivity data have been analyzed using the time-temperature superposition principle. In the Arrhenius diagram, the thermal evolution of the dc-conductivity reveals similar temperature dependence for both protic and aprotic ILs thus making it difficult to distinguish the ion dynamics. However, our results demonstrate the key role of high pressure that unambiguously separates the charge transport properties of protic ILs from aprotic ones through the apparent activation volume parameter. We also highlight that the activation volume can be employed to assess the information connecting the ability of ionic systems to form H-bond networks and the impact of proton transfer involved in the conduction process.

Nonlinear dielectric features of highly polar glass formers: Derivatives of propylene carbonate.

We have measured the nonlinear dielectric behavior of several highly polar propylene carbonate (PC) derivatives in the vicinity of their glass transition temperatures. Focus is on the effects of a large static electric field on the frequency dependent permittivity and on the cubic susceptibility measured using sinusoidal fields of high amplitude. The case of vinyl-PC shows dielectric saturation as well as an electro-rheological effect, i.e., a field induced increase of dielectric relaxation times, whose magnitude changes linearly with the apparent activation energy. The extent of this shift of the loss profile caused by the field correlates strongly with the peak magnitude of the cubic susceptibility, |χ3|, underlining the notion of a link between the |χ3| "hump" and this electro-rheological behavior. Further support for this picture emerges from the observation that the most polar of these liquids, (S)-(-)-methoxy-PC with εs ≈ 250, lacks both the electro-rheological effect in ε″(ω) and the "hump" typically observed in |χ3(ω)|. The absence of any sensitivity of the dynamics to an electric field is contrary to the expectation that the electro-rheological effect correlates with the field induced entropy change, which is extraordinarily high for this liquid. The results suggest that the dependence of the relaxation time on the electric field is not directly linked to the entropy change.