Polymorphic phase transitions: Macroscopic theory and molecular simulation

Tale of Two Polymorphs: Investigating the Structural Differences and Dynamic Relationship between Nirmatrelvir Solid Forms (Paxlovid)

Two anhydrous polymorphs of the novel antiviral medicine nirmatrelvir were discovered during the development of Paxlovid, Pfizer's oral Covid-19 treatment. A comprehensive experimental and computational approach was necessary to distinguish the two closely related polymorphs, herein identified as Forms 1 and 4. This approach paired experimental methods, including powder X-ray diffraction and single-crystal X-ray diffraction, solid-state experimental methods, thermal analysis, solid-state nuclear magnetic resonance and Raman spectroscopy with computational investigations comprising crystal structure prediction, Gibbs free energy calculations, and molecular dynamics simulations of the polymorphic transition. Forms 1 and 4 were ultimately determined to be enantiotropically related polymorphs with Form 1 being the stable form above the transition temperature of ∼17 °C and designated as the nominated form for drug development. The work described in this paper shows the importance of using highly specialized orthogonal approaches to elucidate the subtle differences in structure and properties of similar solid-state forms. This synergistic approach allowed for unprecedented speed in bringing Paxlovid to patients in record time amidst the pandemic.

Ultrasound treatment of crystalline oil-in-water emulsions stabilized by sodium caseinate: Impact on emulsion stability through altered crystallization behavior in the oil globules

Partial coalescence is a key factor contributing to the instability of crystalline oil-in-water emulsions in products like dressings and sauces, reducing shelf life. The intrinsic characteristics of semi-crystalline droplets, including solid fat content, fat crystal arrangement, and polymorphism, play a pivotal role in influencing partial coalescence, challenging prevention efforts even with emulsifiers like amphiphilic proteins. High-intensity ultrasound (HIU) has emerged as an efficient and cost-effective technology for manipulating bulk fat crystallization, thereby enhancing physical properties. This study specifically investigates the impact of HIU treatment on fat crystallization on protein-stabilized crystalline emulsions, utilizing palm olein stearin (POSt) as the lipid phase and sodium caseinate (NaCas) as the surfactant under various HIU powers (100, 150, 200, 300, and 400 W). Results show that increasing HIU power maintained the interfacial potential (-20 mV) provided by NaCas in the emulsions without significant differences. Higher HIU power induced the most stable polymorphic form (β) in the emulsions. Engagingly, the emulsions at 200 W exhibited better storage stability and slower partial coalescence kinetics. Semi-crystalline globules had more uniform and integral crystal clusters that were distributed tangentially near the droplet boundary, perhaps attributed to intermediate subcooling (40.4 °C) at 200 W. The acoustic energy of HIU significantly translates into thermal effects, influencing subcooling degrees as a dominant factor affecting crystallisation in the emulsions. This study establishes ultrasonic crystallization as a novel strategy for modifying the stability of emulsions containing fat crystals.

Bifurcated Polymorphic Transition and Thermochromic Fluorescence of a Molecular Crystal Involving Three-Dimensional Supramolecular Gear Rotation

The ability of molecules to move and rearrange in the solid state accounts for the polymorphic transition and stimuli-responsive properties of molecular crystals. However, how the crystal structure determines the molecular motion ability remains poorly understood. Here, we report that a three-dimensional (3D) supramolecular gear network in the green-emissive polymorph 1G of a dialkylamino-substituted anthracene-pentiptycene π-system (1) enables an unusual bifurcated polymorphic transition into a yellow-emissive polymorph (1Y) and a new green-emissive polymorph (1G*) via 3D correlated supramolecular rotation. The 90° forward correlated rotation causes the molecular conformation between the octyl and the anthracene units to change from syn to anti, the ladder-like supramolecular columns to constrict, and the gear network to disengage. This cooperative molecular motion is marked by the gradual formation of an intermediate state (1I) across the entire crystal from 170 to 230 °C, which then undergoes bifurcated (forward or backward rotation) and irreversible transitions to form polymorphs 1Y and 1G* at 230-235 °C. Notably, 1G* is similar to 1G but lacks gear engagement, preventing its transformation into 1Y. Nevertheless, 1G can be restored by grinding 1Y or 1G* or fuming with dichloromethane (DCM) vapor. This work illustrates the correlation between the crystal structure and solid-state molecular motion behavior and demonstrates how a 3D molecular gear system efficiently transmits thermal energy to drive the polymorphic transition and induce fluorochromism through significant conformational and packing changes.

Direct in situ measurement of polymorphic transition temperatures under thermo-mechanochemical conditions

For the first time, we directly measured the onset and completion temperatures of polymorphic transitions under thermo-mechanochemical conditions by simultaneous in situ synchrotron powder X-ray diffraction and temperature monitoring. We determined the thermo-mechanochemical polymorphic transition temperature in 1-adamantyl-1-diamantyl ether to be 31 °C lower than the transition temperature determined by DSC. Our findings highlight the uniqueness of thermo-mechanochemical conditions, with potential applications in polymorph screening.

Effect of organic acids on the solid-state polymorphic phase transformation of piracetam

Metastable polymorphs are frequently used in oral solid dosage forms to enhance the absorption of poorly water-soluble drug compounds. However, the solid phase transformation from the metastable polymorph to the thermodynamically stable polymorph during manufacturing or storage poses a major challenge for product development and quality control. Here, we report that low-content organic acids can exhibit distinct effects on the solid-state polymorphic phase transformation of piracetam (PCM), a nootropic drug used for memory enhancement. The addition of 1 mol% citric acid (CA) and tricarballylic acid (TA) can significantly inhibit the phase transformation of PCM Form I to Form II, while glutaric acid (GA) and adipic acid (AA) produce a minor effect. A molecular simulation shows that organic acid molecules can adsorb on the crystal surface of PCM Form I, thus slowing the movement of molecules from the metastable form to the stable form. Our study provides deeper insights into the mechanisms of solid-state polymorphic phase transformation of drugs in the presence of additives and facilitates opportunities for controlling the stability of metastable pharmaceuticals.

Investigating the function and design of molecular materials through terahertz vibrational spectroscopy

Terahertz spectroscopy has proved to be an essential tool for the study of condensed phase materials. Terahertz spectroscopy probes the low-frequency vibrational dynamics of atoms and molecules, usually in the condensed phase. These nuclear dynamics, which typically involve displacements of entire molecules, have been linked to bulk phenomena ranging from phase transformations to semiconducting efficiency. The terahertz region of the electromagnetic spectrum has historically been referred to as the 'terahertz gap', but this is a misnomer, as there exist a multitude of methods for accessing terahertz frequencies, and now there are cost-effective instruments that have made terahertz studies much more user-friendly. This Review highlights some of the most exciting applications of terahertz vibrational spectroscopy so far, and provides an in-depth overview of the methods of this technique and its utility to the study of the chemical sciences.

How can machine learning and multiscale modeling benefit ocular drug development?

The eyes possess sophisticated physiological structures, diverse disease targets, limited drug delivery space, distinctive barriers, and complicated biomechanical processes, requiring a more in-depth understanding of the interactions between drug delivery systems and biological systems for ocular formulation development. However, the tiny size of the eyes makes sampling difficult and invasive studies costly and ethically constrained. Developing ocular formulations following conventional trial-and-error formulation and manufacturing process screening procedures is inefficient. Along with the popularity of computational pharmaceutics, non-invasive in silico modeling & simulation offer new opportunities for the paradigm shift of ocular formulation development. The current work first systematically reviews the theoretical underpinnings, advanced applications, and unique advantages of data-driven machine learning and multiscale simulation approaches represented by molecular simulation, mathematical modeling, and pharmacokinetic (PK)/pharmacodynamic (PD) modeling for ocular drug development. Following this, a new computer-driven framework for rational pharmaceutical formulation design is proposed, inspired by the potential of in silico explorations in understanding drug delivery details and facilitating drug formulation design. Lastly, to promote the paradigm shift, integrated in silico methodologies were highlighted, and discussions on data challenges, model practicality, personalized modeling, regulatory science, interdisciplinary collaboration, and talent training were conducted in detail with a view to achieving more efficient objective-oriented pharmaceutical formulation design.

Polymorph control by designed ultrasound application strategy: The role of molecular self-assembly

Molecular self-assembly plays a vital role in the nucleation process and sometimes determines the nucleation outcomes. In this study, ultrasound technology was applied to control polymorph nucleation. For the first time, different ultrasonic application methods based on the nucleation mechanisms have been proposed. For PZA-water and DHB-toluene systems that the molecular self-assembly in solution resembles the synthon in crystal structure, ultrasound pretreatment strategy was conducted to break the original molecular interactions to alter the nucleated form. When the solute molecular self-associates can't give sufficient information to predict the nucleated polymorph like INA-ethanol system, the method of introducing continuous ultrasonic irradiation in the nucleation stage was applied. The induction of ultrasound during nucleation process can break the original interactions firstly by shear forces and accelerate the occurrence of nucleation to avoid the reorientation and rearrangement of solute molecules. These strategies were proved to be effective in polymorph control and have a degree of applicability.

Molecular Mechanism of Organic Crystal Nucleation: A Perspective of Solution Chemistry and Polymorphism

Crystal nucleation determining the formation and assembly pathway of first organic materials is the central science of various scientific disciplines such as chemical, geochemical, biological, and synthetic materials. However, our current understanding of the molecular mechanisms of nucleation remains limited. Over the past decades, the advancements of new experimental and computational techniques have renewed numerous interests in detailed molecular mechanisms of crystal nucleation, especially structure evolution and solution chemistry. These efforts bifurcate into two categories: (modified) classical nucleation theory (CNT) and non-classical nucleation mechanisms. In this review, we briefly introduce the two nucleation mechanisms and summarize current molecular understandings of crystal nucleation that are specifically applied in polymorphic crystallization systems of small organic molecules. Many important aspects of crystal nucleation including molecular association, solvation, aromatic interactions, and hierarchy in intermolecular interactions were examined and discussed for a series of organic molecular systems. The new understandings relating to molecular self-assembly in nucleating systems have suggested more complex multiple nucleation pathways that are associated with the formation and evolution of molecular aggregates in solution.

Exploring the influence factors and improvement strategies of drug polymorphic transformation combined kinetic and thermodynamic perspectives during the formation of nanosuspensions

Nanosuspensions can effectively increase saturation solubility and improve the bioavailability of poorly water-soluble drugs attributed to high loading and surface-to-volume ratio. Wet media milling has been regarded as a scalable method to prepare nanosuspensions because of its simple operation and easy scale-up. In recent years, besides particle aggregation and Ostwald ripening, polymorphic transformation induced by processing has become a critical factor leading to the instability of nanosuspensions. Therefore, this review aims to discuss the influence factors comprehensively and put forward the corresponding improvement strategies of polymorphic transformation during the formation of nanosuspensions. In addition, this review also demonstrates the implication of molecular simulation in polymorphic transformation. The competition between shear-induced amorphization and thermally activated crystallization is the global mechanism of polymorphic transformation during media milling. The factors affecting the polymorphic transformation and corresponding improvement strategies are summarized from formulation and process parameters perspectives during the formation of nanosuspensions. The development of analytical techniques has promoted the qualitative and quantitative characterization of polymorphic transformation, and some techniques can in-situ monitor dynamic transformation. The microhydrodynamic model can be referenced to study the stress intensities by analyzing formulation and process parameters during wet media milling. Molecular simulation can be used to explore the possible polymorphic transformation based on the crystal structure and energy. This review is helpful to improve the stability of nanosuspensions by regulating polymorphic transformation, providing quality assurance for nanosuspension-based products.

Recent advances in drug polymorphs: Aspects of pharmaceutical properties and selective crystallization

Drug polymorphism, an established term used to describe the phenomenon that a drug can exist in different crystalline phases, has attracted great interests in pharmaceutical field in consideration of its important role in affecting the pharmaceutical performance of oral formulations. This paper presents an overview of recent advances in the research on polymorphic drug systems including understandings on nucleation, crystal growth, dissolution, mechanical properties, polymorphic transformation, etc. Moreover, new strategies and mechanisms in the control of polymorphic forms are also highlighted in this review. Furthermore, challenges and trends in the development of polymorphic drugs are briefly discussed, aiming at developing effective and efficient pharmaceutical formulations containing the polymorphic drugs.

The rise of an exciton in solid ammonia

We trace a polymorphic phase change in solid ammonia films through the emergence of a Frenkel exciton at 194.4 nm, for deposition temperatures of 48 K, 50 K and 52 K. Observations on a timescale of hours give unparalleled access to the individual processes of nucleation and the phase change itself. The excitonic transition is forbidden in the low temperature phase, but greater flexing of the solid state structure in the higher temperature phase makes the transition allowed, as the nano-crystals approach ∼30 unit cells through nucleation. We find activation energies of 21.7 ± 0.6 kJ mol^-1 for nucleation and 22.8 ± 0.6 kJ mol^-1 for the phase change, corresponding to the breaking of two to three hydrogen bonds.

Crystal Structure of Colloidally Prepared Metastable Ag2Se Nanocrystals

Structural polymorphism is known for many bulk materials; however, on the nanoscale metastable polymorphs tend to form more readily than in the bulk, and with more structural variety. One such metastable polymorph observed for colloidal Ag₂Se nanocrystals has traditionally been referred to as the "tetragonal" phase. While there are reports on the chemistry and properties of this metastable polymorph, its crystal structure, and therefore electronic structure, has yet to be determined. We report that an anti-PbCl₂-like structure type (space group P2₁/n) more accurately describes the powder X-ray diffraction and X-ray total scattering patterns of colloidal Ag₂Se nanocrystals prepared by several different methods. Density functional theory (DFT) calculations indicate that this anti-PbCl₂-like Ag₂Se polymorph is a dynamically stable, narrow-band-gap semiconductor. The anti-PbCl₂-like structure of Ag₂Se is a low-lying metastable polymorph at 5-25 meV/atom above the ground state, depending on the exchange-correlation functional used.

The mechanism driving a solid-solid phase transition in a biomacromolecular crystal

Solid-solid phase transitions (SSPTs) occur between distinguishable crystalline forms. Because of their importance in application and theory in materials science and condensed-matter physics, SSPTs have been studied most extensively in metallic alloys, inorganic salts and small organic molecular crystals, but much less so in biomacromolecular crystals. In general, the mechanisms of SSPTs at the atomic and molecular levels are not well understood. Here, the ordered molecular rearrangements in biomacromolecular crystals of the adenine riboswitch aptamer are described using real-time serial crystallography and solution atomic force microscopy. Large, ligand-induced conformational changes drive the initial phase transition from the apo unit cell (AUC) to the trans unit cell 1 (TUC1). During this transition, coaxial stacking of P1 duplexes becomes the dominant packing interface, whereas P2-P2 interactions are almost completely disrupted, resulting in 'floating' layers of molecules. The coupling points in TUC1 and their local conformational flexibility allow the molecules to reorganize to achieve the more densely packed and energetically favorable bound unit cell (BUC). This study thus reveals the interplay between the conformational changes and the crystal phases - the underlying mechanism that drives the phase transition. Using polarized video microscopy to monitor SSPTs in small crystals at high ligand concentration, the time window during which the major conformational changes take place was identified, and the in crystallo kinetics have been simulated. Together, these results provide the spatiotemporal information necessary for informing time-resolved crystallography experiments. Moreover, this study illustrates a practical approach to characterization of SSPTs in transparent crystals.

A combined approach to characterize ligand-induced solid-solid phase transitions in biomacromolecular crystals

Solid-solid phase transitions (SSPTs) are widespread naturally occurring phenomena. Understanding the molecular mechanisms and kinetics of SSPTs in various crystalline materials, however, has been challenging due to technical limitations. In particular, SSPTs in biomacromolecular crystals, which may involve large-scale changes and particularly complex sets of interactions, are largely unexplored, yet may have important implications for time-resolved crystallography and for developing synthetic biomaterials. The adenine riboswitch (riboA) is an RNA control element that uses ligand-induced conformational changes to regulate gene expression. Crystals of riboA, upon the addition of a ligand, undergo an SSPT from monoclinic to triclinic to orthorhombic. Here, solution atomic force microscopy (AFM) and polarized video microscopy (PVM) are used to characterize the multiple transition states throughout the SSPT in both the forward and the reverse directions. This contribution describes detailed protocols for growing crystals directly on mica or glass surfaces for AFM and PVM characterization, respectively, as well as methods for image processing and phase-transition kinetics analysis.

Impact of Crystal Size and Morphology on Switchability Characteristics in Pillared-Layer Metal-Organic Framework DUT-8(Ni)

Variation of the crystallite size in flexible porous coordination polymers can significantly influence or even drastically change the flexibility characteristics. The impact of crystal morphology, however, on the dynamic properties of flexible metal-organic frameworks (MOFs) is poorly investigated so far. In the present work, we systematically modulated the particle size of a model gate pressure MOF (DUT-8(Ni), Ni2(2,6-ndc)2(dabco), 2,6-ndc-2,6-naphthalenedicarboxylate, dabco-1,4-diazabicyclo[2.2.2]octane) and investigated the influence of the aspect ratio, length, and width of anisotropically shaped crystals on the gate opening characteristics. DUT-8 is a member of the pillared-layer MOF family, showing reversible structural transition, i.e., upon nitrogen physisorption at 77 K. The framework crystalizes as rod-like shaped crystals in conventional synthesis. To understand which particular crystal surfaces dominate the phenomena observed, crystals similar in size and differing in morphology were involved in a systematic study. The analysis of the data shows that the width of the rods (corresponding to the crystallographic directions along the layer) represents a critical parameter governing the dynamic properties upon adsorption of nitrogen at 77 K. This observation is related to the anisotropy of the channel-like pore system and the nucleation mechanism of the solid-solid phase transition triggered by gas adsorption.

Formulation of ionic liquid APIs via spray drying processes to enable conversion into single and two-phase solid forms

Ionic liquid (IL) forms of drugs are increasingly being explored to address problems presented by poorly water-soluble drugs and solid-state stability. However, before ILs of active pharmaceutical ingredients (APIs) can be routinely incorporated into oral solid dosage forms (OSDs), challenges surrounding their ease of handling and manufacture must be addressed. To this end a framework for transforming API-ILs into solid forms at high loadings based on spray encapsulation using an immiscible polymer has recently been demonstrated. The current work demonstrates that this framework can be applied to a broad range of newly synthesized low glass transition temperature (T_g) API-ILs. Furthermore, the work explores a second novel approach to solidification of API-ILs based on polymer-API-IL miscibility that, to the best of our knowledge, has not been previously demonstrated. Modulated differential scanning calorimetry (mDSC) and attenuated total reflectance Fourier transform infrared spectroscopy showed that it was possible to produce spray dried solid materials, at acceptable loadings and yields for OSD applications in the form of both two-phase phase encapsulated systems and single phase amorphous solid dispersions (ASDs). This was achieved by the appropriate selection of an API-IL insoluble polymer (ethyl cellulose) for phase separated systems, or a miscible polymer with an exceptionally high T_g (the polysaccharide, maltodextrin) for the ASDs. Both approaches successfully overcame the T_g suppression associated with room temperature ILs. This work represents the first step to understanding the fundamental critical physical attributes of these systems to facilitate a more mechanistic methodology for their design.

Dependence of phase transition uniformity on crystal sizes characterized using birefringence

Solid-solid phase transitions (SSPTs) have been widely observed in crystals of organic or inorganic small-molecules. Although SSPTs in macromolecular crystals have been reported, the majority involve local atomic changes, such as those induced by changes in hydration. SSPTs driven by large conformational changes, however, can be more difficult to characterize since they often significantly disrupt lattice packing interactions. Such drastic changes make the cooperativity of molecular motion at the atomic level less easily achieved and more dependent on intrinsic properties of the crystal that define lattice order. Here, we investigate the effect of crystal size on the uniformity of SSPT in thin plate-like crystals of the adenine riboswitch aptamer RNA (riboA) by monitoring changes in crystal birefringence upon the diffusion of adenine ligand. The birefringence intensity is directly related to molecular order and the concurrent changes to polarizability of molecules that results from structural changes throughout the phase transition. The riboA crystals were loosely grouped into three categories (small, medium, and large) based on the surface area of the crystal plates. The time width of transition increased as a function of crystal size, ranging from ∼13 s for small crystals to ∼40 s for the largest crystal. Whereas the transitions in small crystals (<10 μm2) were mostly uniform throughout, the medium and large crystals exhibited large variations in the time and width of the transition peak depending on the region of the crystal being analyzed. Our study provides insight into the spatiotemporal behavior of phase transitions in crystals of biological molecules and is of general interest to time-resolved crystallographic studies, where the kinetics of conformational changes may be governed by the kinetics of an associated SSPT.

Synchronous RNA conformational changes trigger ordered phase transitions in crystals

Time-resolved studies of biomacromolecular crystals have been limited to systems involving only minute conformational changes within the same lattice. Ligand-induced changes greater than several angstroms, however, are likely to result in solid-solid phase transitions, which require a detailed understanding of the mechanistic interplay between conformational and lattice transitions. Here we report the synchronous behavior of the adenine riboswitch aptamer RNA in crystal during ligand-triggered isothermal phase transitions. Direct visualization using polarized video microscopy and atomic force microscopy shows that the RNA molecules undergo cooperative rearrangements that maintain lattice order, whose cell parameters change distinctly as a function of time. The bulk lattice order throughout the transition is further supported by time-resolved diffraction data from crystals using an X-ray free electron laser. The synchronous molecular rearrangements in crystal provide the physical basis for studying large conformational changes using time-resolved crystallography and micro/nanocrystals.

Time-resolved crystallography (TRX) is used for monitoring only small conformational changes of biomacromolecules within the same lattice. Here, the authors report the interplay between synchronous molecular rearrangements and lattice phase transitions in RNA crystals, providing the basis for the investigation of large conformational changes using TRX.

Hydrogen-Bonded Crystalline Molecular Machines with Ultrafast Rotation and Displacive Phase Transitions

Two new crystalline rotors 1 and 2 assembled through N-H⋅⋅⋅N hydrogen bonds by using halogenated carbazole as stators and 1,4-diaza[2.2.2]bicyclooctane (DABCO) as the rotator, are described. The dynamic characterization through ¹ H T₁ relaxometry experiments indicate very low rotational activation barriers (E_a ) of 0.67 kcal mol^-1 for 1 and 0.26 kcal mol^-1 for 2, indicating that DABCO can reach a THz frequency at room temperature in the latter. These E_a values are supported by solid-state density functional theory computations. Interestingly, both supramolecular rotors show a phase transition between 298 and 250 K, revealed by differential scanning calorimetry and single-crystal X-ray diffraction. The subtle changes in the crystalline environment of these rotors that can alter the motion of an almost barrierless DABCO are discussed here.

Solubility prediction for a soluble organic molecule via chemical potentials from density of states

While the solubility of a substance is a fundamental property of widespread significance, its prediction from first principles (starting from only the knowledge of the molecular structure of the solute and solvent) remains a challenge. Recently, we proposed a robust and efficient method to predict the solubility from the density of states of a solute-solvent system using classical molecular simulation. The efficiency, and indeed the generality, of the method has now been enhanced by extending it to calculate solution chemical potentials (rather than probability distributions as done previously), from which solubility may be accessed. The method has been employed to predict the chemical potential of Form 1 of urea in both water and methanol for a range of concentrations at ambient conditions and for two charge models. The chemical potential calculations were validated by thermodynamic integration with the two sets of values being in excellent agreement. The solubility determined from the chemical potentials for urea in water ranged from 0.46 to 0.50 mol kg^-1, while that for urea in methanol ranged from 0.62 to 0.85 mol kg^-1, over the temperature range 298-328 K. In common with other recent studies of solubility prediction from molecular simulation, the predicted solubilities differ markedly from experimental values, reflecting limitations of current forcefields.

Imaging of dehydration in particulate matter using Raman line-focus microscopy

Crystalline solids can incorporate water molecules into their crystal lattice causing a dramatic impact on their properties. This explains the increasing interest in understanding the dehydration pathways of these solids. However, the classical thermal analytical techniques cannot spatially resolve the dehydration pathway of organic hydrates at the single particle level. We have developed a new method for imaging the dehydration of organic hydrates using Raman line-focus microscopy during heating of a particle. Based on this approach, we propose a new metastable intermediate of theophylline monohydrate during the three-step dehydration process of this system and further, we visualize the complex nature of the three-step dehydration pathway of nitrofurantoin monohydrate to its stable anhydrous form. A Raman line-focus mapping option was applied for fast simultaneous mapping of differently sized and shaped particles of nitrofurantoin monohydrate, revealing the appearance of multiple solid-state forms and the non-uniformity of this particle system during the complex dehydration process. This method provides an in-depth understanding of phase transformations and can be used to explain practical industrial challenges related to variations in the quality of particulate materials.

Mapping the polymorphic transformation gateway vibration in crystalline 1,2,4,5-tetrabromobenzene

The thermosalient behavior of 1,2,4,5-tetrabromobenzene (TBB) is related to a temperature-induced polymorphic structural change. The mechanism behind the phase transition has been investigated in this work using low-frequency (10-250 cm-1) Raman spectroscopy and solid-state density functional theory simulations. Careful adjustments of the probing laser power permitted thermal control of the polymorph populations and enabled high-quality Raman vibrational spectra to be obtained for both the β (low temperature) and γ (high temperature) forms of TBB. Numerous well-defined vibrational features appear in the Raman spectra of both polymorphs which could be assigned to specific motions of the solid-state TBB molecules. It was discovered that the lowest-frequency vibration at 15.5 cm-1 in β-TBB at 291 K is a rotational mode that functions as a gateway for inducing the polymorphic phase transition to γ-TBB, and serves as the initiating step in the storage of mechanical strain for subsequent macroscopic release. Computationally mapping the potential energy surface along this vibrational coordinate reveals that the two TBB polymorphs are separated by a 2.40 kJ mol-1 barrier and that γ-TBB exhibits an enhanced cohesion energy that stabilizes its structure.