Consistent and reproducible computation of the glass transition temperature from molecular dynamics simulations.

In many fields, from semiconductors for opto-electronic applications to ionic liquids (ILs) for separations, the glass transition temperature (Tg) of a material is a useful gauge for its potential use in practical settings. As a result, there is a great deal of interest in predicting Tg using molecular simulations. However, the uncertainty and variation in the trend shift method, a common approach in simulations to predict Tg, can be high. This is due to the need for human intervention in defining a fitting range for linear fits of density with temperature assumed for the liquid and glass phases across the simulated cooling. The definition of such fitting ranges then defines the estimate for the Tg as the intersection of linear fits. We eliminate this need for human intervention by leveraging the Shapiro-Wilk normality test and proposing an algorithm to define the fitting ranges and, consequently, Tg. Through this integration, we incorporate into our automated methodology that residuals must be normally distributed around zero for any fit, a requirement that must be met for any regression problem. Consequently, fitting ranges for realizing linear fits for each phase are statistically defined rather than visually inferred, obtaining an estimate for Tg without any human intervention. The method is also capable of finding multiple linear regimes across density vs temperature curves. We compare the predictions of our proposed method across multiple IL and semiconductor molecular dynamics simulation results from the literature and compare other proposed methods for automatically detecting Tg from density-temperature data. We believe that our proposed method would allow for more consistent predictions of Tg. We make this methodology available and open source through GitHub.

Molecular Dynamics Simulation of the Influence of External Electric Fields on the Glass Transition Temperature of the Ionic Liquid 1-Ethyl-3-methylimidazolium Bis(trifluoromethylsulfonyl)imide.

We present the results of molecular dynamics simulations of the ionic liquid (IL) 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [C2C1im][NTf2] in the presence of external electric fields (EEFs) of varying strengths to understand the effects of EEFs on the glass transition temperature Tg. We compute Tg with an automated and objective method and observe a depression in Tg when cooling the IL within an EEF above a critical strength. The effect is reversible, and glasses prepared with EEFs recover their original zero-field Tg when heated. By examining the dynamics and structure of the liquid phase, we find that the EEF lowers the activation energy for diffusion, reducing the energetic barrier for movement and consequently Tg. We show that the effect can be leveraged to drive an electrified nonvapor compression refrigeration cycle.