Picosecond dynamic heterogeneity, hopping, and Johari-Goldstein relaxation in glass-forming liquids

Glassy dynamics in a liquid of anisotropic molecules: Bifurcation of relaxation spectrum

In experimental and theoretical studies of glass transition phenomena, one often finds a sharp crossover in dynamical properties at a temperature Tcr. A bifurcation of a relaxation spectrum is also observed at a temperature TB≈Tcr; both lie significantly above the glass transition temperature. In order to better understand these phenomena, we introduce a new model of glass-forming liquids, a binary mixture of prolate and oblate ellipsoids. This model system exhibits sharp thermodynamic and dynamic anomalies, such as the specific heat jump during heating and a sharp variation in the thermal expansion coefficient around a temperature identified as the glass transition temperature, Tg. The same temperature is obtained from the fit of the calculated relaxation times to the Vogel-Fulcher-Tammann (VFT) form. As the temperature is lowered, the calculated single peak rotational relaxation spectrum splits into two peaks at TB above the estimated Tg. Similar bifurcation is also observed in the distribution of short-to-intermediate time translational diffusion. Interrogation of the two peaks reveals a lower extent of dynamic heterogeneity in the population of the faster mode. We observe an unexpected appearance of a sharp peak in the product of rotational relaxation time τ2 and diffusion constant D at a temperature Tcr, close to TB, but above the glass transition temperature. Additionally, we coarse-grain the system into cubic boxes, each containing, on average, ∼62 particles, to study the average dynamical properties. Clear evidence of large-scale sudden changes in the diffusion coefficient and rotational correlation time signals first-order transitions between low and high-mobility domains.

Local vs. cooperative: Unraveling glass transition mechanisms with SEER

Which phenomenon slows down the dynamics in supercooled liquids and turns them into glasses is a long-standing question of condensed matter. Most popular theories posit that as the temperature decreases, many events must occur in a coordinated fashion on a growing length scale for relaxation to occur. Instead, other approaches consider that local barriers associated with the elementary rearrangement of a few particles or "excitations" govern the dynamics. To resolve this conundrum, our central result is to introduce an algorithm, Systematic Excitation ExtRaction, which can systematically extract hundreds of excitations and their energy from any given configuration. We also provide a measurement of the activation energy, characterizing the liquid dynamics, based on fast quenching and reheating. We use these two methods in a popular liquid model of polydisperse particles. Such polydisperse models are known to capture the hallmarks of the glass transition and can be equilibrated efficiently up to millisecond time scales. The analysis reveals that cooperative effects do not control the fragility of such liquids: the change of energy of local barriers determines the change of activation energy. More generally, these methods can now be used to measure the degree of cooperativity of any liquid model.

Solid-that-Flows Picture of Glass-Forming Liquids

This perspective article reviews arguments that glass-forming liquids are different from those of standard liquid-state theory, which typically have a viscosity in the mPa·s range and relaxation times on the order of picoseconds. These numbers grow dramatically and become 1012 - 1015 times larger for liquids cooled toward the glass transition. This translates into a qualitative difference, and below the "solidity length" which is roughly one micron at the glass transition, a glass-forming liquid behaves much like a solid. Recent numerical evidence for the solidity of ultraviscous liquids is reviewed, and experimental consequences are discussed in relation to dynamic heterogeneity, frequency-dependent linear-response functions, and the temperature dependence of the average relaxation time.

Front propagation in ultrastable glasses is dynamically heterogeneous

Upon heating, ultrastable glassy films transform into liquids via a propagating equilibration front, resembling the heterogeneous melting of crystals. A microscopic understanding of this robust phenomenology is, however, lacking because experimental resolution is limited. We simulate the heterogeneous transformation kinetics of ultrastable configurations prepared using the swap Monte Carlo algorithm, thus allowing a direct comparison with experiments. We resolve the liquid-glass interface both in space and in time as well as the underlying particle motion responsible for its propagation. We perform a detailed statistical analysis of the interface geometry and kinetics over a broad range of temperatures. We show that the dynamic heterogeneity of the bulk liquid is passed on to the front that propagates heterogeneously in space and intermittently in time. This observation allows us to relate the averaged front velocity to the equilibrium diffusion coefficient of the liquid. We suggest that an experimental characterization of the interface geometry during the heterogeneous devitrification of ultrastable glassy films could provide direct experimental access to the long-sought characteristic length scale of dynamic heterogeneity in bulk supercooled liquids.

Fractional Coupling of Primary and Johari-Goldstein Relaxations in a Model Polymer

A polymer model exhibiting heterogeneous Johari−Goldstein (JG) secondary relaxation is studied by extensive molecular-dynamics simulations of states with different temperature and pressure. Time−temperature−pressure superposition of the primary (segmental) relaxation is evidenced. The time scales of the primary and the JG relaxations are found to be highly correlated according to a power law. The finding agrees with key predictions of the Coupling Model (CM) accounting for the decay in a correlation function due to the relaxation and diffusion of interacting systems. Nonetheless, the exponent of the power law, even if it is found in the range predicted by CM (0<ξ<1), deviates from the expected one. It is suggested that the deviation could depend on the particular relaxation process involved in the correlation function and the heterogeneity of the JG process.

Excitations follow (or lead?) density scaling in propylene carbonate

Structural excitations that enable interbasin (IB) barrier crossings on a potential energy landscape are thought to play a facilitating role in the relaxation of liquids. Here, we show that the population of these excitations exhibits the same density scaling observed for α relaxation in propylene carbonate, even though they are heavily influenced by intramolecular modes. We also find that IB crossing modes exhibit a Gr[Formula: see text]neisen parameter ( γ G) that is approximately equivalent to the density scaling parameter γ TS. These observations suggest that the well-documented relationship between γ G and γ TS may be a direct result of the pressure dependence of the frequency of unstable (relaxation) modes associated with IB motion.

Dynamical Manifestations of Supercooling in Amyloid Hydration

The effect of extreme temperature on amyloidogenic species remains sparsely explored. In a recent study (J. Phys. Chem. Lett., 2019, 10, (10)), we employed exhaustive molecular dynamics simulations to explore the cold thermal response of a putative small amyloid oligomer and to elicit the role of solvent modulation. Herein, we investigate the dynamical response of the hydration waters of the oligomer within the supercooled states. Using NMR-based formalism, we delineate the entropic response in terms of the side-chain conformational entropy that corroborates the weakening of the hydrophobic core with lowering of temperature. The translational dynamics of the protein and hydration waters reveal the coupling of protein dynamical fluctuations with solvent dynamics under supercooled conditions. Probing the translational motion as a space-time correlation indicates glassy dynamics exhibited by hydration waters in the supercooled regime. Caging of the water molecules with lowering of temperature and the resultant hopping dynamics are reflected in the longer β-relaxation timescales of translational motion. Furthermore, we utilized mode-coupling theory (MCT) and derived the ideal glass transition temperature from translational and rotational dynamics, around ∼196 and 209 K, respectively. Interestingly, rotational motion in the supercooled regime deviates from the MCT law, exhibits Arrhenius motion, and marks a fragile-to-strong crossover at 227 K. The low-frequency vibrational modes also coincide with the dynamical transition. This exposition lends dynamical insights into the hydration coupling of an amyloid aggregate under cryogenic conditions.

Breakdown of the Stokes-Einstein relation in supercooled water: the jump-diffusion perspective

Although water is the most ubiquitous liquid it shows many thermodynamic and dynamic anomalies. Some of the anomalies further intensify in the supercooled regime. While many experimental and theoretical studies have focused on the thermodynamic anomalies of supercooled water, fewer studies explored the dynamical anomalies very extensively. This is due to the intricacy of the experimental measurement of the dynamical properties of supercooled water. Violation of the Stokes-Einstein relation (SER), an important relation connecting the diffusion of particles with the viscosity of the medium, is one of the major dynamical anomalies. In absence of experimentally measured viscosity, researchers used to check the validity of SER indirectly using average translational relaxation time or α-relaxation time. Very recently, the viscosity of supercooled water was accurately measured at a wide range of temperatures and pressures. This allowed direct verification of the SER at different temperature-pressure thermodynamic state points. An increasing breakdown of the SER was observed with decreasing temperature. Increasing pressure reduces the extent of breakdown. Although some well-known theories explained the above breakdown, a detailed molecular mechanism was still elusive. Recently, a translational jump-diffusion (TJD) approach has been able to quantitatively explain the breakdown of the SER in pure supercooled water and an aqueous solution of methanol. The objective of this article is to present a detailed and state-of-the-art analysis of the past and present works on the breakdown of SER in supercooled water with a specific focus on the new TJD approach for explaining the breakdown of the SER.

Mutual Information in Molecular and Macromolecular Systems

The relaxation properties of viscous liquids close to their glass transition (GT) have been widely characterised by the statistical tool of time correlation functions. However, the strong influence of ubiquitous non-linearities calls for new, alternative tools of analysis. In this respect, information theory-based observables and, more specifically, mutual information (MI) are gaining increasing interest. Here, we report on novel, deeper insight provided by MI-based analysis of molecular dynamics simulations of molecular and macromolecular glass-formers on two distinct aspects of transport and relaxation close to GT, namely dynamical heterogeneity (DH) and secondary Johari–Goldstein (JG) relaxation processes. In a model molecular liquid with significant DH, MI reveals two populations of particles organised in clusters having either filamentous or compact globular structures that exhibit different mobility and relaxation properties. In a model polymer melt, MI provides clearer evidence of JG secondary relaxation and sharper insight into its DH. It is found that both DH and MI between the orientation and the displacement of the bonds reach (local) maxima at the time scales of the primary and JG secondary relaxation. This suggests that, in (macro)molecular systems, the mechanistic explanation of both DH and relaxation must involve rotation/translation coupling.

Nonlinear measurements of kinetics and generalized dynamical modes. I. Extracting the one-dimensional Green's function from a time series

Often, a single correlation function is used to measure the kinetics of a complex system. In contrast, a large set of k-vector modes and their correlation functions are commonly defined for motion in free space. This set can be transformed to the van Hove correlation function, which is the Green's function for molecular diffusion. Here, these ideas are generalized to other observables. A set of correlation functions of nonlinear functions of an observable is used to extract the corresponding Green's function. Although this paper focuses on nonlinear correlation functions of an equilibrium time series, the results are directly connected to other types of nonlinear kinetics, including perturbation-response experiments with strong fields. Generalized modes are defined as the orthogonal polynomials associated with the equilibrium distribution. A matrix of mode-correlation functions can be transformed to the complete, single-time-interval (1D) Green's function. Diagonalizing this matrix finds the eigendecays. To understand the advantages and limitation of this approach, Green's functions are calculated for a number of models of complex dynamics within a Gaussian probability distribution. Examples of non-diffusive motion, rate heterogeneity, and range heterogeneity are examined. General arguments are made that a full set of nonlinear 1D measurements is necessary to extract all the information available in a time series. However, when a process is neither dynamically Gaussian nor Markovian, they are not sufficient. In those cases, additional multidimensional measurements are needed.

Dynamic heterogeneity, cooperative motion, and Johari-Goldstein [Formula: see text]-relaxation in a metallic glass-forming material exhibiting a fragile-to-strong transition

We investigate the Johari-Goldstein (JG) [Formula: see text]-relaxation process in a model metallic glass-forming (GF) material ([Formula: see text]), previously studied extensively by both frequency-dependent mechanical measurements and simulation studies devoted to equilibrium properties, by molecular dynamics simulations based on validated and optimized interatomic potentials with the primary aim of better understanding the nature of this universal relaxation process from a dynamic heterogeneity (DH) perspective. The present relatively low temperature and long-time simulations reveal a direct correspondence between the JG [Formula: see text]-relaxation time [Formula: see text] and the lifetime of the mobile particle clusters [Formula: see text], defined as in previous DH studies, a relationship dual to the corresponding previously observed relationship between the [Formula: see text]-relaxation time [Formula: see text] and the lifetime of immobile particle clusters [Formula: see text]. Moreover, we find that the average diffusion coefficient D nearly coincides with [Formula: see text] of the smaller atomic species (Al) and that the 'hopping time' associated with D coincides with [Formula: see text] to within numerical uncertainty, both trends being in accord with experimental studies. This indicates that the JG [Formula: see text]-relaxation is dominated by the smaller atomic species and the observation of a direct relation between this relaxation process and rate of molecular diffusion in GF materials at low temperatures where the JG [Formula: see text]-relaxation becomes the prevalent mode of structural relaxation. As an unanticipated aspect of our study, we find that [Formula: see text] exhibits fragile-to-strong (FS) glass formation, as found in many other metallic GF liquids, but this fact does not greatly alter the geometrical nature of DH in this material and the relation of DH to dynamical properties. On the other hand, the temperature dependence of the DH and dynamical properties, such as the structural relaxation time, can be significantly altered from 'ordinary' GF liquids.

Experimental evidence of mosaic structure in strongly supercooled molecular liquids

When a liquid is cooled to produce a glass its dynamics, dominated by the structural relaxation, become very slow, and at the glass-transition temperature T_g its characteristic relaxation time is about 100 s. At slightly elevated temperatures (~1.2 T_g) however, a second process known as the Johari-Goldstein relaxation, β_JG, decouples from the structural one and remains much faster than it down to T_g. While it is known that the β_JG-process is strongly coupled to the structural relaxation, its dedicated role in the glass-transition remains under debate. Here we use an experimental technique that permits us to investigate the spatial and temporal properties of the β_JG relaxation, and give evidence that the molecules participating in it are highly mobile and spatially connected in a system-spanning, percolating cluster. This correlation of structural and dynamical properties provides strong experimental support for a picture, drawn from theoretical studies, of an intermittent mosaic structure in the deeply supercooled liquid phase.

The Johari-Goldstein relaxation is a precursor of the glass transition in liquids. Caporaletti et al. use time-dependent interferometry data to substantiate its suggested structural appearance as a globally percolating, fluctuating mosaic.

The polarizability response of a glass-forming liquid reveals intrabasin motion and interbasin transitions on a potential energy landscape

Potential energy landscape (PEL) concepts have been useful in conceptualizing the effects of intermolecular interactions on dynamic and thermodynamic properties of liquids and glasses. "Basins", or regions of reduced potential energy associated with locally preferred molecular packing are important PEL features. The molecular configurations at the bottom of these basins are referred to as inherent structures (ISs). Experimental methods for directly characterizing PEL features such as these are rare, largely relegating PEL concepts to theory and simulation studies, and impeding their exploration in real systems. Recently, we showed that quasielastic neutron scattering (QENS) data from propylene carbonate (PC) exhibit signatures of picosecond timescale motion that are consistent with intrabasin motion and interbasin transitions [Cicerone et al., J. Chem. Phys., 2017, 146, 054502]. Here we present optically-heterodyne-detected optical Kerr effect (OHD-OKE) spectroscopy studies on PC. The data exhibit signatures of motion within and transitions between basins that agree quantitatively with and extend the QENS results. We show that the librational component of the OKE response corresponds to intrabasin dynamics, and the enigmatic intermediate OKE response corresponds to interbasin transition events. The OKE data extend the measurement range of these parameters and reveal their utility in characterizing PEL features of real systems.

Coincident Correlation between Vibrational Dynamics and Primary Relaxation of Polymers with Strong or Weak Johari-Goldstein Relaxation

The correlation between the vibrational dynamics, as sensed by the Debye-Waller factor, and the primary relaxation in the presence of secondary Johari-Goldstein (JG) relaxation, has been investigated through molecular dynamics simulations. Two melts of polymer chains with different bond length, resulting in rather different strength of the JG relaxation are studied. We focus on the bond-orientation correlation function, exhibiting higher JG sensitivity with respect to alternatives provided by torsional autocorrelation function and intermediate scattering function. We find that, even if changing the bond length alters both the strength and the relaxation time of the JG relaxation, it leaves unaffected the correlation between the vibrational dynamics and the primary relaxation. The finding is in harmony with previous studies reporting that numerical models not showing secondary relaxations exhibit striking agreement with experimental data of polymers also where the presence of JG relaxation is known.

Understanding the Origin of the Breakdown of the Stokes-Einstein Relation in Supercooled Water at Different Temperature-Pressure Conditions

A recent experiment has measured the viscosity of water down to approximately 244 K and up to 300 MPa. The correct viscosity and translational diffusivity data at various temperature-pressure (T-P) state points allowed for checking the validity of the Stokes-Einstein (SE) relation, which accounts for the coupling between translational self-diffusion and medium viscosity. The diffusion-viscosity decoupling increases with decreasing temperature, but the increasing pressure reduces the extent of the decoupling. Earlier simulation studies explained the breakdown of the SE relation in terms of the location of the Widom line, emanating from the liquid-liquid critical point (LLCP). Although these studies made a significant contribution to the current understanding of the above phenomena, a detailed molecular picture is still lacking. Recently, our group has explained the diffusion-viscosity decoupling from a jump-diffusion perspective. The jump-diffusion coefficient, emanating from the jump translation of water molecules, is calculated using a quantitative approach for different temperatures at ambient pressure. It has been observed that jump-diffusion is the key factor for diffusion-viscosity decoupling in supercooled water. The same method is adopted in the present work to estimate the jump-diffusion coefficient for different T-P state points and, thereby, explains the role of jump-diffusion for the different extents of the SE relation breakdown at different pressures. The residual diffusion coefficient, the other component of the total diffusion that originates from small step displacement and that is calculated by subtracting the jump-diffusion coefficient from the total diffusion, is seen to be fairly coupled to the viscosity at the entire range of temperature and pressure. Furthermore, we have calculated the average number of H-bonds per water molecule and the tetrahedral order for different T-P state points and investigated an approximate correlation between the average local structure and the contribution of the jump-diffusion to the total diffusion of water. This study, therefore, puts forward a new perspective for explaining the SE relation breakdown in supercooled water under different pressure conditions.

Universal nature of dynamic heterogeneity in glass-forming liquids: A comparative study of metallic and polymeric glass-forming liquids

Glass-formation is a ubiquitous phenomenon that is often observed in a broad class of materials ranging from biological matter to commonly encountered synthetic polymer, as well as metallic and inorganic glass-forming (GF) materials. Despite the many regularities in the dynamical properties of GF materials, the structural origin of the universal dynamical properties of these materials has not yet been identified. Recent simulations of coarse-grained polymeric GF liquids have indicated the coexistence of clusters of mobile and immobile particles that appear to be directly linked, respectively, to the rate of molecular diffusion and structural relaxation. The present work examines the extent to which these distinct types of "dynamic heterogeneity" (DH) arise in metallic GF liquids (Cu-Zr, Ni-Nb, and Pd-Si alloys) having a vastly different molecular structure and chemistry. We first identified mobile and immobile particles and their transient clusters and found the DH in the metallic alloys to be remarkably similar in form to polymeric GF liquids, confirming the "universality" of the DH phenomenon. Furthermore, the lifetime of the mobile particle clusters was found to be directly related to the rate of diffusion in these materials, while the lifetime of immobile particles was found to be proportional to the structural relaxation time, providing some insight into the origin of decoupling in GF liquids. An examination of particles having a locally preferred atomic packing, and clusters of such particles, suggests that there is no one-to-one relation between these populations of particles so that an understanding of the origin of DH in terms of static fluid structure remains elusive.

A microscopic look at the Johari-Goldstein relaxation in a hydrogen-bonded glass-former

Understanding the glass transition requires getting the picture of the dynamical processes that intervene in it. Glass-forming liquids show a characteristic decoupling of relaxation processes when they are cooled down towards the glassy state. The faster (βJG) process is still under scrutiny, and its full explanation necessitates information at the microscopic scale. To this aim, nuclear γ-resonance time-domain interferometry (TDI) has been utilized to investigate 5-methyl-2-hexanol, a hydrogen-bonded liquid with a pronounced βJG process as measured by dielectric spectroscopy. TDI probes in fact the center-of-mass, molecular dynamics at scattering-vectors corresponding to both inter- and intra-molecular distances. Our measurements demonstrate that, in the undercooled liquid phase, the βJG relaxation can be visualized as a spatially-restricted rearrangement of molecules within the cage of their closest neighbours accompanied by larger excursions which reach out at least the inter-molecular scale and are related to cage-breaking events. In-cage rattling and cage-breaking processes therefore coexist in the βJG relaxation.

Structural rearrangements governing Johari-Goldstein relaxations in metallic glasses

The Johari-Goldstein secondary (β) relaxations are an intrinsic feature of supercooled liquids and glasses. They are crucial to many properties of glassy materials, but the underlying mechanisms are still not established. In a model metallic glass, we study the atomic rearrangements by molecular dynamics simulations at time scales of up to microseconds. We find that the distributions of single-particle displacements exhibit multiple peaks, whose positions quantitatively match the pair distribution function. These are identified as the structural signature of cooperative string-like excitations. Furthermore, the most probable time of the string-like motions coincides with the β-relaxation time as probed by dynamical mechanical simulations over a wide temperature range and is consistent with a theoretical model. Our results provide insights into the long-standing puzzle regarding the structural origin of β relaxations in glassy metallic materials.

Enhanced Stabilization in Dried Silk Fibroin Matrices

Preliminary studies have shown that silk fibroin can protect biomacromolecules from thermal degradation, but a deeper understanding of underlying mechanisms needed to fully leverage the stabilizing potential of this matrix has not been realized. In this study, we investigate stabilization of plasma C-reactive protein (CRP), a diagnostic indicator of infection or inflammation, to gain insight into stabilizing mechanisms of silk. We observed that the addition of antiplasticizing excipients that suppress β-relaxation amplitudes in silk matrices resulted in enhanced stability of plasma CRP. These observations are consistent with those made in sugar-glass-based protein-stabilizing matrices and suggest fundamental insight into mechanisms as well as practical strategies to employ with silk protein matrices for enhanced stabilization utility.

Quantifying the sensitivity of the network structure and properties from simultaneous measurements during photopolymerization

We present a method that combines experimental and computational approaches to assess a comprehensive set of structural and functional evolution during a network formation process via photopolymerization. Our work uses the simultaneous measurement of the degree of conversion, polymerization stress, the change in reaction temperature, and shrinkage strain in situ. These measurements are combined with the theory of viscoelastic materials to deduce the relaxation time and frequency-dependent modulus of the polymerizing network. The relaxation time and degree of conversion are used to demonstrate the effect of processing parameters (e.g. curing protocol adjusted by the light intensity) in creating different network structures for the same initial resin. We describe experimental trends using effective medium calculations on a cross-linked polymer network model. In particular, we show that the effect of curing conditions on the spatial heterogeneity in crosslink density can be quantified using multiparametric measurements and modeling. Collectively, the present method is a way to examine holistically the complex structural and functional evolution of the network formation process.

Effect of adding nanometre-sized heterogeneities on the structural dynamics and the excess wing of a molecular glass former

We present the relaxation dynamics of glass-forming glycerol mixed with 1.1 nm sized polyhedral oligomeric silsesquioxane (POSS) molecules using dielectric spectroscopy (DS) and two different neutron scattering (NS) techniques. Both, the reorientational dynamics as measured by DS and the density fluctuations detected by NS reveal a broadening of the α relaxation when POSS molecules are added. Moreover, we find a significant slowing down of the α-relaxation time. These effects are in accord with the heterogeneity scenario considered for the dynamics of glasses and supercooled liquids. The addition of POSS also affects the excess wing in glycerol arising from a secondary relaxation process, which seems to exhibit a dramatic increase in relative strength compared to the α relaxation.

A robust comparison of dynamical scenarios in a glass-forming liquid

We use Bayesian inference methods to provide fresh insights into the sub-nanosecond dynamics of glycerol, a prototypical glass-forming liquid.

Silk-based stabilization of biomacromolecules

Silk fibroin is a high molecular weight amphiphilic protein that self-assembles into robust biomaterials with remarkable properties including stabilization of biologicals and tunable release kinetics correlated to processing conditions. Cells, antibiotics,monoclonal antibodies and peptides, among other biologics, have been encapsulated in silk using various processing approaches and material formats. The mechanistic basis for the entrapment and stabilization features, along with insights into the modulation of release of the entrained compounds from silks will be reviewed with a focus on stabilization of bioactive molecules.

Stabilization of proteins in solid form

Immunogenicity of aggregated or otherwise degraded protein delivered from depots or other biopharmaceutical products is an increasing concern, and the ability to deliver stable, active protein is of central importance. We review characterization approaches for solid protein dosage forms with respect to metrics that are intended to be predictive of protein stability against aggregation and other degradation processes. Each of these approaches is ultimately motivated by hypothetical connections between protein stability and the material property being measured. We critically evaluate correlations between these properties and stability outcomes, and use these evaluations to revise the currently standing hypotheses. Based on this we provide simple physical principles that are necessary (and possibly sufficient) for generating solid delivery vehicles with stable protein loads. Essentially, proteins should be strongly coupled (typically through H-bonds) to the bulk regions of a phase-homogeneous matrix with suppressed β relaxation. We also provide a framework for reliable characterization of solid protein forms with respect to stability.

Using the fluorescence red edge effect to assess the long-term stability of lyophilized protein formulations

Nanosecond relaxation processes in sugar matrices are causally linked through diffusional processes to protein stability in lyophilized formulations. Long-term protein degradation rates track mean-squared displacement (⟨u(2)⟩) of hydrogen atoms in sugar glasses, a parameter describing dynamics on a time scale of picoseconds to nanoseconds. However, measurements of ⟨u(2)⟩ are usually performed by neutron scattering, which is not conducive to rapid formulation screening in early development. Here, we present a benchtop technique to derive a ⟨u(2)⟩ surrogate based on the fluorescence red edge effect. Glycerol, lyophilized trehalose, and lyophilized sucrose were used as model systems. Samples containing 10(-6) mole fraction of rhodamine 6G, a fluorophore, were excited at either 532 nm (main peak) or 566 nm (red edge), and the ⟨u(2)⟩ surrogate was determined based the corresponding Stokes shifts. Results showed reasonable agreement between ⟨u(2)⟩ from neutron scattering and the surrogate from fluorescence, although deviations were observed at very low temperatures. We discuss the sources of the deviations and suggest technique improvements to ameliorate these. We expect that this method will be a valuable tool to evaluate lyophilized sugar matrices with respect to their ability to protect proteins from diffusion-limited degradation processes during long-term storage. Additionally, the method may have broader applications in amorphous pharmaceutical solids.