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 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.

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.

Universal Behavior of Dielectric Responses of Glass Formers: Role of Dipole-Dipole Interactions.

From an exhaustive examination of the molecular dynamics in practically all van der Waals molecular glass formers ever probed by dielectric spectroscopy, we found that the width of the α-loss peak at or near the glass transition temperature T_{g} is strongly anticorrelated with the polarity of the molecule. The larger the dielectric relaxation strength Δε(T_{g}) of the system, the narrower is the α-loss peak. This remarkable property is explained by the contribution from the dipole-dipole interaction potential V_{dd}(r)=-Dr^{-6} to the attractive part of the intermolecular potential, making the resultant potential more harmonic, and the effect increases rapidly with the dipole moment µ and Δε(T_{g}) in view of the relation, D∝(µ^{4}/kT_{g})∝kT_{g}[Δε(T_{g})]^{2}. Since the novel correlation discovered encompasses practically all van der Waals molecular glass formers studied by dielectric spectroscopy, it impacts the large dielectric research community as well as those engaged in solving the glass transition problem.

Effects of dynamic heterogeneity and density scaling of molecular dynamics on the relationship among thermodynamic coefficients at the glass transition.

In this paper, we define and experimentally verify thermodynamic characteristics of the liquid-glass transition, taking into account a kinetic origin of the process. Using the density scaling law and the four-point measure of the dynamic heterogeneity of molecular dynamics of glass forming liquids, we investigate contributions of enthalpy, temperature, and density fluctuations to spatially heterogeneous molecular dynamics at the liquid-glass transition, finding an equation for the pressure coefficient of the glass transition temperature, dTg/dp. This equation combined with our previous formula for dTg/dp, derived solely from the density scaling criterion, implies a relationship among thermodynamic coefficients at Tg. Since this relationship and both the equations for dTg/dp are very well validated using experimental data at Tg, they are promising alternatives to the classical Prigogine-Defay ratio and both the Ehrenfest equations in case of the liquid-glass transition.

The complex, non-monotonic thermal response of the volumetric space of simple liquids.

In this paper, an intricate effect of the isothermal compression of simple liquids on their volumetric response is reported for α,ω-halogenoalkanes as examples. We apply an accurate experimental technique, scanning transitiometry, to directly measure isobaric thermal volume expansivities αp of the liquids in a wide density range. To thoroughly analyze the observed intersection of the experimental isothermal pressure dependences of αp, we develop a class of equations of state derived in the density scaling regime for molecular dynamics, finding successful temperature parameterizations of an isothermal equation of state (EOS) intrinsically adapted to describe volumetric data in an extremely wide density range. The EOS based analyses of the scanning transitiometry data as a function of temperature T and pressure p undoubtedly show that the previously considered crossing point of the isothermal dependences αp(p) is in general represented by a non-linear and non-monotonic curve in the (T-p) phase diagram.

On the scaling behavior of electric conductivity in [C4mim][NTf2].

In this work we examine, for the first time, the molar conductivity behavior of the deeply supercooled room temperature ionic liquid [C4mim][NTf2] in the temperature, pressure and volume thermodynamic space in terms of density scaling (TV(γ))(-1) combined with the equation of state (EOS). The exponent γσ determined from the Avramov model analysis is compared with the coefficient obtained from the viscosity studies carried out at moderate temperatures. Therefore, the experimental results presented herein provide the answer to the long-standing question regarding the validity of thermodynamic scaling of ionic liquids over a wide temperature range, i.e. from the normal liquid state to the glass transition point. Finally, we investigate the relationship between the dynamic and thermodynamic properties of [C4mim][NTf2] represented by scaling exponent γ and Grüneisen constant γG, respectively.

Equation of state in the generalized density scaling regime studied from ambient to ultra-high pressure conditions.

In this paper, based on the effective intermolecular potential with well separated density and configuration contributions and the definition of the isothermal bulk modulus, we derive two similar equations of state dedicated to describe volumetric data of supercooled liquids studied in the extremely wide pressure range related to the density range, which is extremely wide in comparison with the experimental range reached so far in pressure-volume-temperature measurements of glass-forming liquids. Both the equations comply with the generalized density scaling law of molecular dynamics versus h(ρ)/T at different densities ρ and temperatures T, where the scaling exponent can be in general only a density function γ(ρ) = d ln h/d ln ρ as recently argued by the theory of isomorphs. We successfully verify these equations of state by using data obtained from molecular dynamics simulations of the Kob-Andersen binary Lennard-Jones liquid. As a very important result, we find that the one-parameter density function h(ρ) analytically formulated in the case of this prototypical model of supercooled liquid, which implies the one-parameter density function γ(ρ), is able to scale the structural relaxation times with the value of this function parameter determined by fitting the volumetric simulation data to the equations of state. We also show that these equations of state properly describe the pressure dependences of the isothermal bulk modulus and the configurational isothermal bulk modulus in the extremely wide pressure range investigated by the computer simulations. Moreover, we discuss the possible forms of the density functions h(ρ) and γ(ρ) for real glass formers, which are suggested to be different from those valid for the model of supercooled liquid based on the Lennard-Jones intermolecular potential.

General rules prospected for the liquid fragility in various material groups and different thermodynamic conditions.

The fragility parameter has been acknowledged as one of the most important characteristics of glass-forming liquids. We show that the mystery of the dramatic change in molecular dynamics of systems approaching the glass transition can be better understood by the high pressure study of fragility parameters defined in different thermodynamic conditions. We formulate and experimentally confirm a few rules obeyed by the fragility parameters, which are also rationalized by the density scaling law and its modification suggested for associated liquids. In this way, we successfully explore and gain a new insight into the pressure effect on molecular dynamics of van der Waals liquids, polymer melts, ionic liquids, and hydrogen-bonded systems near the glass transition.

Effect of temperature and density fluctuations on the spatially heterogeneous dynamics of glass-forming Van der Waals liquids under high pressure.

In this Letter, we show how temperature and density fluctuations affect the spatially heterogeneous dynamics at ambient and elevated pressures. By using high-pressure experimental data for van der Waals liquids, we examine contributions of the temperature and density fluctuations to the dynamics heterogeneity. We show that the dynamic heterogeneity decreases significantly with increasing pressure at a constant structural relaxation time (isochronal condition), while the broadening of the relaxation spectrum remains constant. This observation questions the relationship between spectral broadening and dynamic heterogeneity.

Relaxation dynamics and crystallization study of sildenafil in the liquid and glassy states.

In this paper, the physical stability and molecular dynamics of amorphous sildenafil are investigated in both the liquid and glassy states. We have established that the amorphous sildenafil is resistant to recrystallization at temperatures below the glass transition temperature Tg during the experimental period of its storage (i.e., above 6 months), however, it easily undergoes cold crystallization at T > Tg. To determine the crystallization mechanism, the isothermal and non-isothermal studies of the cold crystallization kinetics of the drug are performed by using the broadband dielectric spectroscopy (BDS) and the differential scanning calorimetry (DSC), respectively. The cold crystallization mechanism has been found to be similar in both the isothermal and non-isothermal cases. This mechanism has been analyzed from the point of view of the molecular mobility of sildenafil investigated in the supercooled liquid and glassy states by using the BDS measurements in the wide temperature range. This analysis has been enriched with a new approach based on a recently reported measure of dynamic heterogeneity given by a four-point dynamic susceptibility function. No tendency to recrystallization of glassy sildenafil at T < Tg is also discussed in relation to molecular dynamics of sildenafil in the glassy state. The relatively small molecular mobility reflected in one secondary relaxation as well as the predicted large time scale of structural relaxation of glassy sildenafil suggests that amorphous sildenafil should not recrystallize during its long-term storage at room temperature.

Enhancement of amorphous celecoxib stability by mixing it with octaacetylmaltose: the molecular dynamics study.

In this paper, we present a novel way of stabilization of amorphous celecoxib (CEL) against recrystallization by preparing binary amorphous celecoxib-octaacetylmaltose (CEL-acMAL) systems by quench-cooling of the molten phase. As far as we know this is the first application of carbohydrate derivatives with acetate groups to enhance the stability of an amorphous drug. We found that CEL in the amorphous mixture with acMAL is characterized by a much better solubility than pure CEL. We report very promising results of the long-term measurements of stability of the CEL-acMAL binary amorphous system with small amount of stabilizer during its storage at room temperature. Moreover, we examined the effect of adding acMAL on molecular dynamics of CEL in the wide temperature range in both the supercooled liquid and glassy states. We found that the molecular mobility of the mixture of CEL with 10 wt % acMAL in the glassy state is much more limited than that in the case of pure CEL, which correlates with the better stability of the amorphous binary system. By dielectric measurements and theoretical calculations within the framework of density functional theory (DFT), we studied the role of acMAL in enhancing the stability of amorphous CEL in mixtures and postulated which interactions between CEL and acMAL molecules can be responsible for preventing devitrification.

Mechanism of mutarotation in supercooled liquid phase: Studies on L-sorbose.

We have studied mutarotation in anhydrous supercooled L-sorbose by means of dielectric spectroscopy. The phenomenon observed in L-sorbose is much faster than in the structurally similar D-fructose. The kinetics of this process has been determined by applying 1st order kinetics model. Activation energy equal to 68 kJ/mol was obtained from temperature dependence of rate constants. To understand differences in mutarotation rate between D-fructose and L-sorbose, quantum mechanical calculations were performed to study mechanism of this phenomenon. The possible impact of water absorbed from air on the mutarotation in supercooled liquid state has been checked. It turned out that the process is probably intermolecular and the water molecules or other carbohydrate molecules assist in the proton transfer process. Finally we have shown that the rate constant can be alternatively determined from frequency of the maximum of peak, obtained by performing Fourier transform of kinetic curve.

High pressure study of molecular dynamics of protic ionic liquid lidocaine hydrochloride.

In this paper, we investigate the effect of pressure on the molecular dynamics of protic ionic liquid lidocaine hydrochloride, a commonly used pharmaceutical, by means of dielectric spectroscopy and pressure-temperature-volume methods. We observed that near T(g) the pressure dependence of conductivity relaxation times reveals a peculiar behavior, which can be treated as a manifestation of decoupling between ion migration and structural relaxation times. Moreover, we discuss the validity of thermodynamic scaling in lidocaine HCl. We also employed the temperature-volume Avramov model to determine the value of pressure coefficient of glass transition temperature, dT(g)/dP|(P = 0.1). Finally, we investigate the role of thermal and density fluctuations in controlling of molecular dynamics of the examined compound.

Heart Transplantation as a Result of Pediatric Inflammatory Multisystem Syndrome in an Adolescent.

Several reviews have shown that COVID-19 in children is a relatively mild disease. However, a rare complication affecting children and adolescents after COVID-19 has been identified. Pediatric inflammatory multisystem syndrome temporally associated with COVID-19 (PIMS-TS), which in some cases manifests itself as a hyperinflammatory syndrome with a multiorgan failure, may lead to death. We report a case of a 17-year-old patient who was admitted to the hospital with cardiogenic shock of unknown etiology. The disease was life-threatening, thus necessitating mechanical ventilation, circulatory support, and extracorporeal therapy due to renal and liver dysfunction. The patient tested negative for SARS-CoV-2 Reverse Transcription Polymerase Chain Reaction. Other infectious causes of illness were excluded. However, the patient had a positive IgG antibody test result and high levels of interleukin-6, which helped to diagnose PIMS-TS. Intravenous immunoglobulin and steroid therapies were initiated, unfortunately, with poor outcome. The patient's critical condition, particularly end-stage heart failure, led to mechanical circulatory support implantation and finally orthotopic heart transplantation. After the surgery, the patient's condition improved gradually. PIMS-TS manifests itself with different clinical images and as a state of varying severity, ultimately causing multiorgan dysfunction with shock resembling toxic shock syndrome. Ultimately, myocardial complications of PIMS-TS necessitated heart transplantation in the described patient.

Outdoor exposure in children from Buenos Aires Province, Argentina.

OBJECTIVES: To evaluate myopia risk factors, mainly outdoor exposure and reading habits, in a country with low prevalence of myopia (Buenos Aires Province, Argentina). METHODS: Consecutive children interviewed in a clinical private practice setting were autorefracted under cycloplegia with cyclopentolate 1%. Their parents consented to fill a questionnaire about schooling, tutorial classes, outdoor exposure, reading habits, and cellphone use, both on weekdays and weekends. The Spanish questionnaire was based on past English questionnaires of myopia clinical trials. The spherical equivalent of the right eye was used for the refractive distribution. The average daily hours spent for each activity were calculated. RESULTS: This study involved 115 children aged 10.48 ±â€¯3.65 years (range 5-18 years), with 56.5% being girls. Children had 8 h of schooling per day in 62.6% of cases, and only 14.8 % had tutorial classes after school. There were 38.3% myopes (<-0.50 D), 24.3% hyperopes (>+2.00 D) and the rest were emmetropes. The mean time that these children spent outdoors per day was 3.94 ±â€¯1.45 h (27.60 ±â€¯10.16 h per week). The total mean time spent reading and writing per day was 1.50 ±â€¯0.98 h, and that spent using cellphones and tablets was 2.43 ±â€¯1.66 h. CONCLUSION: In an environment with low myopia prevalence, children spend almoast 4 h per day outdoors, much more than the usual recommendation of 2 h a day for myopia prevention.

Communication: Relationships between Intermolecular potential, thermodynamics, and dynamic scaling in viscous systems.

In this communication, we provide a recipe for a consistent relation between dynamic scaling and thermodynamic properties well-grounded by the same intermolecular generalized Lennard-Jones potential, which is derived by using an essentially modified Avramov model within the framework of the "thermodynamic scaling" idea. This relation is experimentally verified very well for supercooled van der Waals liquids, and consequently, it can be a good basis for a proper universal description of molecular dynamics and thermodynamics of viscous systems.

Study of molecular dynamics of the pharmaceutically important protic ionic liquid verapamil hydrochloride. II. Test of entropic models.

In this work we analyzed the structural relaxation times as a function of both temperature and pressure in terms of the entropic models by using dielectric and PVT measurements data presented in our previous research on the ionic liquid verapamil hydrochloride [Z. Wojnarowska, M. Paluch, A. Grzybowski, et al., J. Chem. Phys. 131, 104505 (2009)]. Two different approaches were used to analyze the tau(alpha)(T,P) dependence: the modified Avramov model as well as the pressure extended Adam-Gibbs model in the forms proposed by Casalini (AG(C)) and Schwartz (AG(S)). In every case a satisfactory description of the structural relaxation times was achieved. Additionally, using both mentioned models the pressure dependence of the fragility m(P) and the glass transition temperature T(g) were determined. We also compared the value of dT(g)/dP|(P=0) calculated on the basis of the considered entropic models with the experimental value evaluated in our recent work. Consequently, we were able to estimate which of the examined models in the best way relates the dynamic to the thermodynamic parameters.

Virial-potential-energy correlation and its relation to density scaling for quasireal model systems.

In this paper, we examine the virial- and the potential-energy correlation for quasireal model systems. This correlation constitutes the framework of the theory of the isomorph in the liquid phase diagram commonly examined using simple liquids. Interestingly, our results show that for the systems characterized by structural anisotropy and flexible bonds, the instantaneous values of total virial and total potential energy are entirely uncorrelated. It is due to the presence of the intramolecular interactions because the contributions to the virial and potential energy resulting from the intermolecular interactions still exhibit strong linear dependence. Interestingly, in contrast to the results reported for simple liquids, the slope of the mentioned linear dependence is different than the values of the density scaling exponent. However, our findings show that for quasireal materials, the slope of dependence between the virial and potential energy (resulting from the intermolecular interactions) strongly depends on the interval of intermolecular distances that are taken into account. Consequently, the value of the slope of the discussed relationship, which enables satisfactory density scaling, can be obtained. Interestingly, this conclusion is supported by the results obtained for analogous systems without intermolecular attraction, for which the value the slope of the virial-potential-energy correlation is independent of considered intermolecular distances, directly corresponds to the exponent of the intermolecular repulsion, and finally leads to accurate density scaling.