Solid-state Density Functional Theory Research Articles

Donepezil hydrochloride fingerprint spectral and formation mechanism analysis.

Donepezil hydrochloride is a widely used medication for treating Alzheimer's disease. This study utilized terahertz time-domain spectroscopy to analyze the fingerprint spectra of donepezil hydrochloride, identifying five characteristic absorption peaks at 1.65, 2.44, 2.56, 3.31, and 3.75THz. The vibrational spectrum of the donepezil hydrochloride crystal was further examined using solid-state density functional theory. Based on simulation calculations, the characteristic peaks were identified and analyzed in detail, focusing on long-range ordering and weak interaction networks. The results demonstrate that terahertz spectroscopy is an effective tool for studying intermolecular interactions in drug crystals.

Low-Frequency Vibrational Spectra of 4-Fluorophenol Studied by THz Spectroscopy and Solid-State Density Functional Theory.

Hydrogen bonds serve as important intermolecular interactions organizing the spatial arrangement of molecular crystals. Determining the hydrogen-bond orientations remains a challenging task. In the previous XRD study, the authors assumed a single crystal structure of 4-fluorophenol in which molecules form hexamer-ring clusters via anticlockwise hydrogen bonding. However, the existence of the other structure, which adopts clockwise hydrogen bonding, remains uncertain and warrants further exploration. We unveil distinctive fingerprint information associated with both structures by using high-resolution terahertz spectroscopy. The distinct structures result in different intramolecular geometries of 4-fluorophenol regarding the O-H bond configurations, which are manifested by a noticeable peak splitting in the 20-200 cm-1 frequency range. This result illustrates the sensitivity of THz spectroscopy to hydrogen-bond conformational polymorphism. Our findings represent a significant advancement in utilizing high-resolution THz spectroscopy to resolve hydrogen-bond orientations in molecular crystals with possible broad applicability across diverse hydrogen-bonded organic and inorganic framework materials.

Surface Tailoring of MoS2 Nanosheets with Substituted Aromatic Diazonium Salts for Gas Sensing: A DFT Study.

Two-dimensional (2D) nanomaterials are useful for building gas sensors owing to their desirable electronic and optical properties. However, they usually suffer from selectivity, because they cannot discriminate between gas molecules. Functionalization with organic molecules can be used to tailor their surfaces to recognize a specific family of compounds. In this study, solid-state density functional theory (DFT) was used to elucidate the functionalization of MoS2 with substituted aromatic diazonium salts (R = -H, - CH3, -CO2H, -CHO, -OCH3, and -NO2). Results showed that chemical reaction with diazonium salts is favored to their physical adsorption (E ads = -0.04 to -0.38 eV vs E rxn = -1.47 to -2.20 eV), where organic cations have a preference to attach atop of sulfur atoms. Chemical functionalization induced a small variation in the bandgap energy not exceeding 0.04 eV; thus, the optical properties were well preserved. In the presence of ammonia, the substituted MoS 2 /2(a-f) responded to the target analyte through a change in the interaction energy, varying from -0.08 to -0.83 eV, where the best interaction energy was obtained for MoS 2 /2c, bearing the carboxylic acid group. In the presence of other gases such as CO2, SO2, and H2S, the interaction energy is lower (-0.14 to -0.35 eV), indicating good selectivity of the nanomaterials. Furthermore, the interaction increased in the presence of humidity, which was more realistic than that in the presence of neat NH3. This interaction was confirmed by computing the partial charges. Recovery times estimated from the interaction energies ranged from 0.31 s to several minutes, depending on the interacting molecules. Phenylcarboxyl-modified MoS2 nanosheets show great potential as candidates for the development of chemoresistive gas sensors that are specifically designed for detecting ammonia.

Structure of vanadia-titania catalysts doped with alkali metals according to solid-state NMR, ESR spectroscopies and DFT

In this work, V/M/TiO2 catalysts (M = Li, Na, K, Rb and Cs) were investigated with Solid-state nuclear magnetic resonance (SSNMR), electron paramagnetic resonance (ESR) and first-principles calculations. The total vanadium content corresponds to half monolayer (4 V atoms·nm−2 of TiO2 support surface), with V/M atomic ratio being 4/0.6. Static 51V and MAS 1H, 7Li, 51V, 23Na, 133Cs NMR spectra were acquired. Surface moieties of vanadium oxide on (001) anatase surface were modeled using density functional theory (DFT); their 51V NMR parameters were calculated with the Gauge-Including Projected Augmented Wave (GIPAW) method and compared to experimental data obtained with 51V SSNMR spectroscopy. It was found that mainly strongly bound vanadium sites (V3) are formed on anatase samples, located on TiO2 terminal oxygen atoms. Weakly bound vanadium sites (V1, V2) are formed on bridged anatase oxygen atoms. Supported alkali metals (Li, Na, K, Rb and Cs) on TiO2 were found to interact with protons located on the terminal and bridged TiO2 oxygen atoms. The relative ratio of weakly bound vanadium sites (V1, V2) increased on alkali-containing samples. Calculations performed showed lithium cation to prefer V-O-Ti or V-O-V bonds over VO, unsystematically deshielding vanadium nuclei. The presence of V3+ is likely to cause a loss of intensity in 51V NMR spectra.

Zirconium-Based Metal-Organic Frameworks with Free Hydroxy Groups for Enhanced Perfluorooctanoic Acid Uptake in Water.

Perfluorooctanoic acid (PFOA) is a highly recalcitrant organic pollutant, and its bioaccumulation severely endangers human health. While various methods are developed for PFOA removal, the targeted design of adsorbents with high efficiency and reusability remains largely unexplored. Here the rational design and synthesis of two novel zirconium-based metal‒organic frameworks (MOFs) bearing free ortho-hydroxy sites, namely noninterpenetrated PCN-1001 and twofold interpenetrated PCN-1002,are presented. Single crystal analysis of the pure ligand reveals that intramolecular hydrogen bonding plays a pivotal role in directing the formation of MOFs with free hydroxy groups. Furthermore, the transformation from PCN-1001 to PCN-1002 is realized. Compared to PCN-1001, PCN-1002 displays higher chemical stability due to interpenetration, thereby demonstrating an exceptional PFOA adsorption capacity of up to 632mgg-1 (1.53mmolg-1), which is comparable to the reported record values. Moreover, PCN-1002 shows rapid kinetics, high selectivity, and long-life cycles in PFOA removal tests. Solid-state nuclear magnetic resonance results and density functional theory calculations reveal that multiple hydrogen bonds between the free ortho-hydroxy sites and PFOA, along with Lewis acid-base interaction, work collaboratively to enhance PFOA adsorption.

Ion Dynamics at the Intermediate Charging State of the Sodium Vanadium Fluorophosphate Cathode.

Na super ionic conductor (NASICON)-type polyanionic vanadium fluorophosphate Na3V2O2(PO4)2F (NVOPF) is a promising cathode material for high-energy sodium-ion batteries. The dynamic diffusion and exchange of sodium ions in the lattice of NVOPF are crucial for its electrochemical performance. However, standard characterizations are mostly focused on the as-synthesized material without cycling, which is different from the actual battery operation conditions. In this work, we investigated the hopping processes of sodium in NVOPF at the intermediate charging state with 23Na solid-state nuclear magnetic resonance (ssNMR) and density functional theory (DFT) calculations. Our experimental characterizations revealed six distinct sodium coordination sites in the intermediate structure and determined the exchange rates among these sites at variable temperatures. The theoretical calculations showed that these dynamic processes correspond to different ion transport pathways in the crystalline lattice. Our combined experimental and theoretical study uncovered the underlying mechanisms of the ion transport in cycled NVOPF and these understandings may help the optimization of cathode materials for sodium-ion batteries.

A Blyholder mechanism in the chemisorption of N2O on Ni(111) – studied with Auger-photoelectron coincidence spectroscopy

In heterogeneous catalysis the surface-adsorbate bond strength is critical for the function of the system. Here we study a series consisting of multilayer, bilayer and monolayer N2O on Ni(111) and employ Auger-photoelectron coincidence spectroscopy (APECS) to study the interaction between the molecule and the substrate directly. We observe intensity in the nitrogen Auger spectra that arise from the interaction between molecule and surface (not observed in free molecules) whereas the oxygen spectra are thickness-independent. Since the two nitrogen atoms of N2O are chemically inequivalent we can assign the intensity present in the bilayer and monolayer cases to orbitals centered on the terminal nitrogen which is closest to the Ni(111) surface. Using ab initio, molecular dynamics and solid-state density functional theory calculations we infer a Blyholder model of the surface bond as consisting of donation from the terminal nitrogen lone-pair valence orbital with back-donation from the metal into the unoccupied orbitals on that nitrogen. This coincidence technique can readily be used to study substrate–adsorbate interactions directly with chemical and orbital specificity — this opens up prospects to study fundamental steps of molecular adsorption and heterogeneous catalysis with unprecedented detail.

Where do the Fluorine Atoms Go in Inorganic-Oxide Fluorinations? A Fluorooxoborate Illustration under Terahertz Light.

The substitution of fluorine atoms for oxygen atoms/hydroxyl groups has emerged as a promising strategy to enhance the physical and chemical properties of oxides/hydroxides in fluorine chemistry. However, distinguishing fluorine from oxygen/hydroxyl in the reaction products poses a significant challenge in existing characterization methods. In this study, we illustrate that terahertz (THz) spectroscopy provides a powerful tool for addressing this challenge. To this end, we investigated two fluorination reactions of boric acid, utilizing MHF2 (M=Na, C(NH2)3) as fluorine reagents. Through an interplay between THz spectroscopy and solid-state density functional theory, we have conclusively demonstrated that fluorine atoms exclusively bind with the sp3-boron but not with the sp2-boron in the reaction products of Na[B(OH)3][B3O3F2(OH)2] (NaBOFH) and [C(NH2)3]2B3O3F4OH (GBF2). Based on this evidence, we have proposed a reaction pathway for the fluorinations under investigation, a process previously hindered due to structural ambiguity. This work represents a step forward in gaining a deeper understanding of the precise structures and reaction mechanisms involved in the fluorination of oxides/hydroxides, illuminated by the insights provided by THz spectroscopy.

Machine-Learned Potentials by Active Learning from Organic Crystal Structure Prediction Landscapes.

A primary challenge in organic molecular crystal structure prediction (CSP) is accurately ranking the energies of potential structures. While high-level solid-state density functional theory (DFT) methods allow for mostly reliable discrimination of the low-energy structures, their high computational cost is problematic because of the need to evaluate tens to hundreds of thousands of trial crystal structures to fully explore typical crystal energy landscapes. Consequently, lower-cost but less accurate empirical force fields are often used, sometimes as the first stage of a hierarchical scheme involving multiple stages of increasingly accurate energy calculations. Machine-learned interatomic potentials (MLIPs), trained to reproduce the results of ab initio methods with computational costs close to those of force fields, can improve the efficiency of the CSP by reducing or eliminating the need for costly DFT calculations. Here, we investigate active learning methods for training MLIPs with CSP datasets. The combination of active learning with the well-developed sampling methods from CSP yields potentials in a highly automated workflow that are relevant over a wide range of the crystal packing space. To demonstrate these potentials, we illustrate efficiently reranking large, diverse crystal structure landscapes to near-DFT accuracy from force field-based CSP, improving the reliability of the final energy ranking. Furthermore, we demonstrate how these potentials can be extended to more accurately model structures far from lattice energy minima through additional on-the-fly training within Monte Carlo simulations.

An NMR crystallographic characterisation of solid (+)-usnic acid.

We use a combination of one- and two-dimensional solid-state nuclear magnetic resonance (NMR) spectroscopy and density functional theory (DFT) calculations to obtain a full assignment of the 1H and 13C signals for solid (+)-usnic acid, which contains two molecules in the asymmetric unit. By combining through-space 1H-1H correlation data with computation it is possible to assign signals not just to the same molecules (relative assignment) but to assign the signals to specific crystallographic molecules (absolute assignment). Variable-temperature measurements reveal that there is some variation in many of the 13C chemical shifts with temperature, likely arising from varying populations of different tautomeric forms of the molecule. The NMR spectrum of crystalline (+)-usnic acid is then compared with that of ground Usnea dasopoga lichen (the source material of the usnic acid). The abundance of usnic acid is so great in the lichen that this natural product can be observed directly in the NMR spectrum without further purification. This natural sample of usnic acid appears to have the same crystalline form as that in the pure commercial sample.

1H and 13C chemical shift-structure effects in anhydrous β-caffeine and four caffeine-diacid cocrystals probed by solid-state NMR experiments and DFT calculations.

By using density functional theory (DFT) calculations, we refined the H atom positions in the structures of β-caffeine (C), α-oxalic acid (OA; (COOH)2), α-(COOH)2·2H2O, β-malonic acid (MA), β-glutaric acid (GA), and I-maleic acid (ME), along with their corresponding cocrystals of 2 : 1 (2C-OA, 2C-MA) or 1 : 1 (C-GA, C-ME) stoichiometry. The corresponding 13C/1H chemical shifts obtained by gauge including projector augmented wave (GIPAW) calculations agreed overall very well with results from magic-angle-spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy experiments. Chemical-shift/structure trends of the precursors and cocrystals were examined, where good linear correlations resulted for all COO1H sites against the H⋯O and/or H⋯N H-bond distance, whereas a general correlation was neither found for the aliphatic/caffeine-stemming 1H sites nor any 13C chemical shift against either the intermolecular hydrogen- or tetrel-bond distance, except for the 13COOH sites of the 2C-OA, 2C-MA, and C-GA cocrystals, which are involved in a strong COOH⋯N bond with caffeine that is responsible for the main supramolecular stabilization of the cocrystal. We provide the first complete 13C NMR spectral assignment of the structurally disordered anhydrous β-caffeine polymorph. The results are discussed in relation to previous literature on the disordered α-caffeine polymorph and the ordered hydrated counterpart, along with recommendations for NMR experimentation that will secure sufficient 13C signal-resolution for reliable resonance/site assignments.

Electrolyte Mixtures and Organic Cathode Materials for Rechargeable Aluminum Batteries

Rechargeable aluminum (Al) metal batteries are an emerging energy storage technology ideal for use at a global scale: Al metal is energy dense, low cost, inherently safe, earth abundant, and highly recyclable. Despite such great promise, their technological progress has been hindered by the few electrolytes able to reversibly electrodeposit Al metal at room temperature, coupled with the limited number of positive electrode materials that are compatible with them. In this context, recent progress will be discussed in the development of electrolyte mixtures and organic cathode materials for rechargeable aluminum batteries. Chloroaluminate ionic liquid electrolytes and ionic liquid analogues prepared using mixtures of organic cations and/or neutral solvent species will be presented, which exhibit improved electrochemical properties for aluminum electrodeposition and in aluminum-graphite batteries, over temperatures ranging from ambient conditions down to -60 °C. Different organic structures will also be presented as organic cathode materials for aluminum batteries, where the effects of molecular structure on bulk energy storage properties will be discussed. Molecular-scale understanding of their ionic and electronic charge storage mechanisms will be elucidated through a combination of solid-state nuclear magnetic resonance (NMR) spectroscopy and density functional theory (DFT) calculations. Lastly, by using different electrolytes in combination with an organic cathode, it will be shown that electrolyte speciation can affect the local environments of charge-compensating ions in cycled organic electrodes. Overall, the results and analyses are aimed at developing next-generation rechargeable aluminum metal batteries for diverse energy storage applications.

Predicting 51V nuclear magnetic resonance observables in molecular crystals.

Solid-state nuclear magnetic resonance (NMR) spectroscopy and quantum chemical density functional theory (DFT) calculations are widely used to characterize vanadium centers in biological and pharmaceutically relevant compounds. Several techniques have been recently developed to improve the accuracy of predicted NMR parameters obtained from DFT. Fragment-based and planewave-corrected methods employing hybrid density functionals are particularly effective tools for solid-state applications. A recent benchmark study involving molecular crystal compounds found that fragment-based NMR calculations using hybrid density functionals improve the accuracy of predicted 51V chemical shieldings by 20% relative to traditional planewave methods. This work extends the previous study, including a careful analysis of 51V chemical shift anisotropy, electric field gradient calculations, and a more extensive test set. The accuracy of planewave-corrected techniques and recently developed fragment-based methods using electrostatic embedding based on the polarized continuum model (PCM) are found to be highly competitive with previous methods. Planewave-corrected methods achieve a 34% improvement in the errors of predicted 51V chemical shieldings relative to planewave. Additionally, planewave-corrected and fragment-based calculations were performed using PCM embedding, improving the accuracy of predicted 51V chemical shielding (CS) tensor principal values by 30% and values by 15% relative to traditional planewave methods. The performance of these methods is further examined using a redox-active oxovandium complex and a common 51V NMR reference compound.

Charge-clustering induced fast ion conduction in 2LiX-GaF3: A strategy for electrolyte design.

2LiX-GaF3 (X = Cl, Br, I) electrolytes offer favorable features for solid-state batteries: mechanical pliability and high conductivities. However, understanding the origin of fast ion transport in 2LiX-GaF3 has been challenging. The ionic conductivity order of 2LiCl-GaF3 (3.20 mS/cm) > 2LiBr-GaF3 (0.84 mS/cm) > 2LiI-GaF3 (0.03 mS/cm) contradicts binary LiCl (10-12 S/cm) < LiBr (10-10 S/cm) < LiI (10-7 S/cm). Using multinuclear 7Li, 71Ga, 19F solid-state nuclear magnetic resonance and density functional theory simulations, we found that Ga(F,X)n polyanions boost Li+-ion transport by weakening Li+-X- interactions via charge clustering. In 2LiBr-GaF3 and 2LiI-GaF3, Ga-X coordination is reduced with decreased F participation, compared to 2LiCl-GaF3. These insights will inform electrolyte design based on charge clustering, applicable to various ion conductors. This strategy could prove effective for producing highly conductive multivalent cation conductors such as Ca2+ and Mg2+, as charge clustering of carboxylates in proteins is found to decrease their binding to Ca2+ and Mg2+.

Synergistic effect of surface oxygen vacancies and hydroxyl groups on Cu-doped TiO2 photocatalyst for hydrogen evolution

The intricate nature of the titanium dioxide (TiO2) photocatalytic system, influenced by Cu doping, often yields conflicting outcomes. This study addresses this complexity by examining the impact of Cu ion doping in nanostructured TiO2 to augment its photocatalytic efficacy in hydrogen generation from ethanol-water mixtures. The improved performance is primarily attributed to the generation of oxygen vacancies (Ovac) within the TiO2 lattice. By adjusting the thermal processing parameters, the chemical state of Cu is modulated to optimize photocatalytic performance. Solid-state nuclear magnetic resonance (SSNMR) and density functional theory (DFT) calculations confirm that the presence of Ovac and hydroxyl groups (OH−) surface species not only facilitates local spatial charge separation but also enhances light absorption, thereby improving photocatalytic hydrogen evolution reaction (HER) rate up to 20.2 mmol h−1 g−1 under UV light. Furthermore, a proposed mechanism involves the participation of Ovac and OH− species, supported by Mott-Schottky (M − S) plots and ultraviolet photoelectron spectroscopy (UPS), revealing that Cu–TiO2 exhibits a lower flat band potential (Efb) and smaller work function compared to TiO2. This research underscores the pivotal role of Cu ions and Ovac:OH− in advancing the photocatalytic activity of TiO2-based crystals, imparting noteworthy implications for sustainable hydrogen generation.

Diamine-appended metal-organic framework for carbon capture from wet flue gas: Characteristics and mechanism

Due to the superior CO2 adsorption capacity and selectivity, amine-functional solid adsorbents have received increasing attention. Nevertheless, amine-functionalized polyporous materials have a large performance deviation when considering carbon removal capacity under wet flue gas. Therefore, a variety of metal–organic frameworks (MOFs) with amine functionality are synthesized in this work by appending diamine with a propylene linker to Mg2(dobpdc) (dobpdc4- = 4,4′-dioxido-3,3′-biphenyldicarboxylate) and screened under the simulated flue gas environment with a high water content (7 v%). The results indicate that ndmpn-grafted Mg2(dobpdc) (ndmpn = N,N-dimethyl-1,3-propanediamine) has a high CO2 adsorption capacity of 1.72 mmol/g at 60 °C with an extremely low regeneration heat (2.26 GJ/ton CO2 at the desorption temperature of 145 °C). Adsorption property of ndmpn-Mg2(dobpdc) at high humidity is even more excellent than that of another well-reported diaminopropene-appended MOF, dmpn-Mg2(dobpdc) (dmpn = 2,2-dimethyl-1,3-propanediamine). Herein, van der Waals (vdW)-corrected density functional theory (DFT) calculations are used to quantitatively analyze binding strengths and water effects on it between diamine-appended Mg2(dobpdc) and CO2. It is revealed that the addition of water makes ammonium carbamate chains produced in ndmpn-Mg2(dobpdc) more stable. These chemisorption structures under humid conditions have also been disclosed by using 15N solid-state nuclear magnetic resonance (NMR) spectroscopy experiments and DFT calculation. To conclude, we unveil the binary adsorption characteristics and mechanism of ndmpn-Mg2(dobpdc). It enables a thorough understanding of this compound and highlights its potential as a promising adsorbent under humid conditions for post-combustion applications.

One-Year Water-Stable and Porous Bi(III) Halide Semiconductor with Broad-Spectrum Antibacterial Performance.

Hybrid metal halide semiconductors are a unique family of materials with immense potential for numerous applications. For this to materialize, environmental stability and toxicity deficiencies must be simultaneously addressed. We report here a porous, visible light semiconductor, namely, (DHS)Bi2I8 (DHS = [2.2.2] cryptand), which consists of nontoxic, earth-abundant elements, and is water-stable for more than a year. Gas- and vapor-sorption studies revealed that it can selectively and reversibly adsorb H2O and D2O at room temperature (RT) while remaining impervious to N2 and CO2. Solid-state NMR measurements and density functional theory (DFT) calculations verified the incorporation of H2O and D2O in the molecular cages, validating the porous nature. In addition to porosity, the material exhibits broad band-edge light emission centered at 600 nm with a full width at half-maximum (fwhm) of 99 nm, which is maintained after 6 months of immersion in H2O. Moreover, (DHS)Bi2I8 exhibits bacteriocidal action against three Gram-positive and three Gram-negative bacteria, including antibiotic-resistant strains. This performance, coupled with the recorded water stability and porous nature, renders it suitable for a plethora of applications, from solid-state batteries to water purification and disinfection.

Probing low-frequency terahertz calculation: A comparison between periodic boundary conditions and cluster models

Terahertz spectroscopy performs as a powerful tool for revealing valuable information on molecular interactions. In recent years, low-frequency terahertz simulation through periodic boundary conditions has proved to be a great success with regard to long-range force and complete noncovalent interaction rather than cluster models or unit cell models. However, quantum chemistry simulation with a limited number of molecules also enables many in-depth and comprehensive analytical methods, therefore, encouraging proving the outcome validity of the cluster models based on quantum chemistry calculation. In this article, we compared the computational results of both periodic boundary conditions and cluster models through solid-state density functional theory, revealing the necessity of a periodic structure in calculating low-frequency terahertz bands, which is instrumental for the selection of computational models for terahertz spectra.

Terahertz spectroscopy analysis of L-Phenylalanine and its fluorinated derivatives

3-Fluoro-L-Phenylalanine (3-F-L-Phe) and 4-Fluoro-L-Phenylalanine (4-F-L-Phe) are fluorinated derivatives of L-Phenylalanine (L-Phe) that hold significant medical applications. To investigate the impact of introducing fluorine atoms on the structure of L-Phe, this study employed terahertz (THz) spectroscopy to distinguish three crystal structure differences. A comparative analysis of the absorption peaks of 3-F-L-Phe and 4-F-L-Phe with that of L-Phe revealed varying degrees of red shift and blue shift. The dielectric loss analysis demonstrated that the substitution of fluorine atoms in different positions of the benzene ring changes the crystal structure and electric dipole moment. By combining lattice network and vibration mode comparative analysis, the study concluded that the differences in the spectrum could be attributed to variations in the hydrogen bond network and changes in the unit cell volume. These findings highlight the effectiveness of using terahertz time domain spectroscopy (THz-TDS) and solid-state density functional theory (DFT) as a new method for the qualitative analysis of fluorinated amino acid derivatives.

Low-Frequency Raman Spectroscopy of Pure and Cocrystallized Mycophenolic Acid.

The aqueous solubility of solid-state pharmaceuticals can often be enhanced by cocrystallization with a coformer to create a binary cocrystal with preferred physical properties. Greater understanding of the internal and external forces that dictate molecular structure and intermolecular packing arrangements enables more efficient design of new cocrystals. Low-frequency (sub-200 cm-1) Raman spectroscopy experiments and solid-state density functional theory simulations have been utilized together to investigate the crystal lattice vibrations of mycophenolic acid, an immunosuppressive drug, in its pure form and as a cocrystal with 2,2'-dipyridylamine. The lattice vibrations primarily consist of large-amplitude translations and rotations of the crystal components, thereby providing insights into the critical intermolecular forces governing cohesion of the molecular solids. The simulations reveal that despite mycophenolic acid having a significantly unfavorable conformation in the cocrystal as compared to the pure solid, the cocrystal exhibits greater thermodynamic stability over a wide temperature range. The energetic penalty due to the conformational strain is more than compensated for by the strong intermolecular forces between the drug and 2,2'-dipyridylamine. Quantifying the balance of internal and external energy factors in cocrystal formation indicates a path forward in the development of future mycophenolic acid cocrystals.