Supersoft Norbornene-Based Thermoplastic Elastomers with High Strength and Upper Service Temperature.

With over 6 million tons produced annually, thermoplastic elastomers (TPEs) have become ubiquitous in modern society, due to their unique combination of elasticity, toughness, and reprocessability. Nevertheless, industrial TPEs display a tradeoff between softness and strength, along with low upper service temperatures, typically ≤100 °C. This limits their utility, such as in bio-interfacial applications where supersoft deformation is required in tandem with strength, in addition to applications that require thermal stability (e.g., encapsulation of electronics, seals/joints for aeronautics, protective clothing for firefighting, and biomedical devices that can be subjected to steam sterilization). Thus, combining softness, strength, and high thermal resistance into a single versatile TPE has remained an unmet opportunity. Through de novo design and synthesis of novel norbornene-based ABA triblock copolymers, this gap is filled. Ring-opening metathesis polymerization is employed to prepare TPEs with an unprecedented combination of properties, including skin-like moduli (<100 kPa), strength competitive with commercial TPEs (>5 MPa), and upper service temperatures akin to high-performance plastics (≈260 °C). Furthermore, the materials are elastic, tough, reprocessable, and shelf stable (≥2 months) without incorporation of plasticizer. Structure-property relationships identified herein inform development of next-generation TPEs that are both biologically soft yet thermomechanically durable.

Unlocking Direct Lithium Extraction in Harsh Conditions through Thiol-Functionalized Metal-Organic Framework Subnanofluidic Membranes.

Metal-organic framework (MOF) membranes with high ion selectivity are highly desirable for direct lithium-ion (Li+) separation from industrial brines. However, very few MOF membranes can efficiently separate Li+ from brines of high Mg2+/Li+ concentration ratios and keep stable in ultrahigh Mg2+-concentrated brines. This work reports a type of MOF-channel membranes (MOFCMs) by growing UiO-66-(SH)2 into the nanochannels of polymer substrates to improve the efficiency of MOF membranes for challenging Li+ extraction. The resulting membranes demonstrate excellent monovalent metal ion selectivity over divalent metal ions, with Li+/Mg2+ selectivity up to 103 since Mg2+ should overcome a higher energy barrier than Li+ when transported through the MOF pores, as confirmed by molecular dynamics simulations. Under dual-ion diffusion, as the Mg2+/Li+ mole ratio of the feed solution increases from 0.2 to 30, the membrane Li+/Mg2+ selectivity decreases from 1516 to 19, corresponding to the purity of lithium products between 99.9 and 95.0%. Further research on multi-ion diffusion that involves Mg2+ and three monovalent metal ions (K+, Na+, and Li+, referred to as M+) in the feed solutions shows a significant improvement in Li+/Mg2+ separation efficiency. The Li+/Mg2+ selectivity can go up to 1114 when the Mg2+/M+ molar concentration ratio is 1:1, and it remains at 19 when the ratio is 30:1. The membrane selectivity is also stable for 30 days in a highly concentrated solution with a high Mg2+/Li+ concentration ratio. These results indicate the feasibility of the MOFCMs for direct lithium extraction from brines with Mg2+ concentrations up to 3.5 M. This study provides an alternative strategy for designing efficient MOF membranes in extracting valuable minerals in the future.

Ion and Water Dynamics in the Transition from Dry to Wet Conditions in Salt-Doped PEG.

The influence of the water content on ion and water transport mechanisms in polymer membranes under low to moderate hydration conditions remains poorly understood. In this study, we combine ion and water diffusivity (PFG-NMR) measurements with atomistic molecular dynamics simulations to better understand transport processes in hydrated salt-doped poly(ethylene glycol). Above the water percolation threshold, the experimental and simulated diffusivities are in good agreement with the free volume transport models. At low hydration levels, unlike dry systems, ion diffusion cannot be described by polymer segmental dynamics alone. We rationalize such observations by the interplay between ion-water and ion-polymer solvation of cations and between ion-water and cation-anion interactions for anions. Further, we demonstrate that a two-state model combining ion-water solvation and free volume transport can describe water dynamics across the entire hydration range of interest. Our findings provide a more encompassing analysis of ion and water transport in hydrated polyelectrolytes, specifically in the low hydration regime.

Structural elucidation of polydopamine facilitated by ionic liquid solvation.

Minimal understanding of the formation mechanism and structure of polydopamine (pDA) and its natural analogue, eumelanin, impedes the practical application of these versatile polymers and limits our knowledge of the origin of melanoma. The lack of conclusive structural evidence stems from the insolubility of these materials, which has spawned significantly diverse suggestions of pDA's structure in the literature. We discovered that pDA is soluble in certain ionic liquids. Using these ionic liquids (ILs) as solvents, we present an experimental methodology to solvate pDA, enabling us to identify pDA's chemical structure. The resolved pDA structure consists of self-assembled supramolecular aggregates that contribute to the increasing complexity of the polymer. The underlying molecular energetics of pDA solvation and a macroscopic picture of the disruption of the aggregates using IL solvents have been investigated, along with studies of the aggregation mechanism in water.

Accounting for Ion Pairing Effects on Sulfate Salt Sorption in Cation Exchange Membranes.

Ion exchange membranes (IEMs) are frequently used in water treatment and electrochemical applications, with their ion separation properties largely governed by equilibrium ion partitioning between a membrane and contiguous solution. Despite an expansive literature on IEMs, the influence of electrolyte association (i.e., ion pairing) on ion sorption remains relatively unexplored. In this study, salt sorption in two commercial cation exchange membranes equilibrated with 0.01-1.0 M MgSO4 and Na2SO4 is investigated experimentally and theoretically. Association measurements of salt solutions using conductometric experiments and the Stokes-Einstein approximation show significant concentrations of ion pairs in MgSO4 and Na2SO4 relative to those in simple electrolytes (i.e., NaCl), which is consistent with prior studies of sulfate salts. The Manning/Donnan model, developed and validated for halide salts in previous studies, substantially underpredicts sulfate sorption measurements, presumably due to ion pairing effects not accounted for in this established theory. These findings suggest that ion pairing can enhance salt sorption in IEMs due to partitioning of reduced valence species. By reformulating the Donnan and Manning models, a theoretical framework for predicting salt sorption in IEMs that explicitly considers electrolyte association is developed. Remarkably, theoretical predictions of sulfate sorption are improved by over an order of magnitude by accounting for ion speciation. In some cases, good quantitative agreement is observed between theoretical and experimental values for external salt concentrations between 0.1 and 1.0 M using no adjustable parameters.

An artificial sodium-selective subnanochannel.

Single-ion selectivity with high precision has long been pursued for fundamental bioinspired engineering and applications such as in ion separation and energy conversion. However, it remains a challenge to develop artificial ion channels to achieve single-ion selectivity comparable to their biological analogs, especially for high Na+/K+ selectivity. Here, we report an artificial sodium channel by subnanoconfinement of 4'-aminobenzo-15-crown-5 ethers (15C5s) into ~6-Å-sized metal-organic framework subnanochannel (MOFSNC). The resulting 15C5-MOFSNC shows an unprecedented Na+/K+ selectivity of tens to 102 and Na+/Li+ selectivity of 103 under multicomponent permeation conditions, comparable to biological sodium channels. A co-ion-responsive single-file transport mechanism in 15C-MOFSNC is proposed for the preferential transport of Na+ over K+ due to the synergetic effects of size exclusion, charge selectivity, local hydrophobicity, and preferential binding with functional groups. This study provides an alternative strategy for developing potential single-ion selective channels and membranes for many applications.

Construction of angstrom-scale ion channels with versatile pore configurations and sizes by metal-organic frameworks.

Controllable fabrication of angstrom-size channels has been long desired to mimic biological ion channels for the fundamental study of ion transport. Here we report a strategy for fabricating angstrom-scale ion channels with one-dimensional (1D) to three-dimensional (3D) pore structures by the growth of metal-organic frameworks (MOFs) into nanochannels. The 1D MIL-53 channels of flexible pore sizes around 5.2 × 8.9 Å can transport cations rapidly, with one to two orders of magnitude higher conductivities and mobilities than MOF channels of hybrid pore configurations and sizes, including Al-TCPP with 1D ~8 Å channels connected by 2D ~6 Å interlayers, and 3D UiO-66 channels of ~6 Å windows and 9 - 12 Å cavities. Furthermore, the 3D MOF channels exhibit better ion sieving properties than those of 1D and 2D MOF channels. Theoretical simulations reveal that ion transport through 2D and 3D MOF channels should undergo multiple dehydration-rehydration processes, resulting in higher energy barriers than pure 1D channels. These findings offer a platform for studying ion transport properties at angstrom-scale confinement and provide guidelines for improving the efficiency of ionic separations and nanofluidics.

Multimorphic Materials: Spatially Tailoring Mechanical Properties via Selective Initiation of Interpenetrating Polymer Networks.

Access to multimaterial polymers with spatially localized properties and robust interfaces is anticipated to enable new capabilities in soft robotics, such as smooth actuation for advanced medical and manufacturing technologies. Here, orthogonal initiation is used to create interpenetrating polymer networks (IPNs) with spatial control over morphology and mechanical properties. Base catalyzes the formation of a stiff and strong polyurethane, while blue LEDs initiate the formation of a soft and elastic polyacrylate. IPN morphology is controlled by when the LED is turned "on", with large phase separation occurring for short time delays (≈1-2 min) and a mixed morphology for longer time delays (>5 min), which is supported by dynamic mechanical analysis, small angle X-ray scattering, and atomic force microscopy. Through tailoring morphology, tensile moduli and fracture toughness can be tuned across ≈1-2 orders of magnitude. Moreover, a simple spring model is used to explain the observed mechanical behavior. Photopatterning produces "multimorphic" materials, where morphology is spatially localized with fine precision (<100 µm), while maintaining a uniform chemical composition throughout to mitigate interfacial failure. As a final demonstration, the fabrication of hinges represents a possible use case for multimorphic materials in soft robotics.

Surface Modification of Ultrafiltration Membranes with 1,4-Benzoquinone and Polyetheramines to Improve Fouling Resistance.

Membrane fouling remains a key challenge for membrane separations. Hydrophilic membrane surface modification can mitigate irreversible foulant deposition, thereby improving fouling resistance. We report new hydrophilic membrane coatings based on 1,4-benzoquinone and various commercially available polyetheramines. These coatings, prepared from 1,4-benzoquinone and Jeffamine EDR 148, poly(benzoquinone-Jeffamine EDR 148) (p(BQ-EDR 148)), were used to modify polysulfone (PS) ultrafiltration membranes. In fouling experiments using an oil/water emulsion, membranes exhibited comparable fouling resistance to that of polydopamine (pDA)-modified membranes. Based on contact angle measurements, p(BQ-EDR 148) and pDA-modified membranes have similar levels of hydrophilicity, and both exhibited higher threshold flux values than those of their unmodified analogues. Based on their similar threshold flux values, p(BQ-EDR 148)-modified (76 LMH) and pDA-modified membranes (74 LMH) should have similar fouling resistance. Moreover, the mean pore size of p(BQ-EDR 148)-modified membranes can be tuned, while keeping the pure water permeance constant, by changing the deposition time and molar ratio of benzoquinone to EDR 148 in the modification solution.

The Donnan potential revealed.

Selective transport of solutes across a membrane is critical for many biological, water treatment and energy conversion and storage systems. When a charged membrane is equilibrated with an electrolyte, an unequal distribution of ions arises between phases, generating the so-called Donnan electrical potential at the solution/membrane interface. The Donnan potential results in the partial exclusion of co-ion, providing the basis of permselectivity. Although there are well-established ways to indirectly estimate the Donnan potential, it has been widely reported that it cannot be measured directly. Here we report the first direct measurement of the Donnan potential of an ion exchange membrane equilibrated with salt solutions. Our results highlight the dependence of the Donnan potential on external salt concentration and counter-ion valence, and show a reasonable agreement with current theoretical models of IEMs, which incorporate ion activity coefficients. By directly measuring the Donnan potential, we eliminate ambiguities that arise from limitations inherent in current models.

Impact of Ion-Ion Correlated Motion on Salt Transport in Solvated Ion Exchange Membranes.

The influence of dynamical ion-ion correlations and ion pairing on salt transport in ion exchange membranes remain poorly understood. In this study, we use the framework of Onsager transport coefficients within atomistic molecular dynamics simulations to study the impact of ion-ion correlated motion on salt transport in hydrated polystyrene sulfonate membranes and compare with the results from aqueous salt solutions. At sufficiently high salt concentrations, cation-anion dynamical correlations exert a significant influence on both salt diffusivities and conductivities. Anion-anion distinct correlations, arising from the imbalance between the concentration of free (mobile) cations and anions, and the retarding effect of the fixed charge groups on cations, proves to be an additional important feature for polymer membranes. Our results demonstrate that dynamical correlations should become an important consideration in experimental measurements of salt diffusivities and conductivities for non-dilute salt solutions in polymer membranes.

Polymeric multimaterials by photochemical patterning of crystallinity.

An organized combination of stiff and elastic domains within a single material can synergistically tailor bulk mechanical properties. However, synthetic methods to achieve such sophisticated architectures remain elusive. We report a rapid, facile, and environmentally benign method to pattern strong and stiff semicrystalline phases within soft and elastic matrices using stereo-controlled ring-opening metathesis polymerization of an industrial monomer, cis-cyclooctene. Dual polymerization catalysis dictates polyolefin backbone chemistry, which enables patterning of compositionally uniform materials with seamless stiff and elastic interfaces. Visible light-induced activation of a metathesis catalyst results in the formation of semicrystalline trans polyoctenamer rubber, outcompeting the formation of cis polyoctenamer rubber, which occurs at room temperature. This bottom-up approach provides a method for manufacturing polymeric materials with promising applications in soft optoelectronics and robotics.

Free Volume and Permeability of Mixed Matrix Membranes Made from a Terbutil-M-terphenyl Polyamide and a Porous Polymer Network.

A set of thermally rearranged mixed matrix membranes (TR-MMMs) was manufactured and tested for gas separation. These membranes were obtained through the thermal treatment of a precursor MMM with a microporous polymer network and an o-hydroxypolyamide,(HPA) created through a reaction of 2,2-bis(3-amino-4-hydroxyphenyl)-hexafluoropropane (APAF) and 5'-terbutil-m-terfenilo-3,3″-dicarboxylic acid dichloride (tBTmCl). This HPA was blended with different percentages of a porous polymer network (PPN) filler, which produced gas separation MMMs with enhanced gas permeability but with decreased selectivity. The thermal treatment of these MMMs gave membranes with excellent gas separation properties that did not show the selectivity decreasing trend. It was observed that the use of the PPN load brought about a small decrease in the initial mass losses, which were lower for increasing PPN loads. Regarding the glass transition temperature, it was observed that the use of the filler translated to a slightly lower Tg value. When these MMMs and TR-MMMs were compared with the analogous materials created from the isomeric 5'-terbutil-m-terfenilo-4,4″-dicarboxylic acid dichloride (tBTpCl), the permeability was lower for that of tBTmCl, compared with the one from tBTpCl, although selectivity was quite similar. This fact could be attributed to a lower rigidity as roughly confirmed by the segmental length of the polymer chain as studied by WAXS. A model for FFV calculation was proposed and its predictions compared with those evaluated from density measurements assuming a matrix-filler interaction or ideal independence. It turns out that permeability as a function of FFV for TR-MMMs follows an interaction trend, while those not thermally treated follow the non-interaction trend until relatively high PPN loads were reached.

Gas Permeability, Fractional Free Volume and Molecular Kinetic Diameters: The Effect of Thermal Rearrangement on ortho-hydroxy Polyamide Membranes Loaded with a Porous Polymer Network.

Mixed-matrix membranes (MMMs) consisting of an ortho-hydroxy polyamide (HPA) matrix, and variable loads of a porous polymer network (PPN) were thermally treated to induce the transformation of HPA to polybenzoxazole (ß-TR-PBO). Two different HPAs were synthesized to be used as a matrix, 6FCl-APAF and tBTpCl-APAF, while the PPN used as a filler was prepared by reacting triptycene and trifluoroacetophenone. The permeability of He, H2, N2, O2, CH4 and CO2 gases through these MMMs are analyzed as a function of the fraction of free volume (FFV) of the membrane and the kinetic diameter of the gas, allowing for the evaluation of the free volume. Thermal rearrangement entails an increase in the FFV. Both before and after thermal rearrangement, the free volume increases with the PPN content very similarly for both polymeric matrices. It is shown that there is a portion of free volume that is inaccessible to permeation (occluded volume), probably due to it being trapped within the filler. In fact, permeability and selectivity change below what could be expected according to densities, when the fraction of occluded volume increases. A higher filler load increases the percentage of inaccessible or trapped free volume, probably due to the increasing agglomeration of the filler. On the other hand, the phenomenon is slightly affected by thermal rearrangement. The fraction of trapped free volume seems to be lower for membranes in which the tBTpCl-APAF is used as a matrix than for those with a 6FCl-APAF matrix, possibly because tBTpCl-APAF could approach the PPN better. The application of an effective medium theory for permeability allowed us to extrapolate for a 100% filler, giving the same value for both thermally rearranged and non-rearranged MMMs. The pure filler could also be extrapolated by assuming the same tendency as in the Robeson's plots for MMMs with low filler content.

Hydrogen Recovery by Mixed Matrix Membranes Made from 6FCl-APAF HPA with Different Contents of a Porous Polymer Network and Their Thermal Rearrangement.

Mixed matrix membranes (MMMs) consisting of a blend of a hydroxypolyamide (HPA) matrix and variable loads of a porous polymer network (PPN) were thermally treated to induce the transformation of HPA to polybenzoxazole (ß-TR-PBO). Here, the HPA matrix was a hydroxypolyamide having two hexafluoropropyilidene moieties, 6FCl-APAF, while the PPN was prepared by reacting triptycene (TRP) and trifluoroacetophenone (TFAP) in a superacid solution. The most probable size of the PPN particles was 75 nm with quite large distributions. The resulting membranes were analyzed by SEM and AFM. Up to 30% PPN loads, both SEM and AFM images confirmed quite planar surfaces, at low scale, with limited roughness. Membranes with high hydrogen permeability and good selectivity for the gas pairs H2/CH4 and H2/N2 were obtained. For H2/CO2, selectivity almost vanished after thermal rearrangement. In all cases, their hydrogen permeability increased with increasing loads of PPN until around 30% PPN with ulterior fairly abrupt decreasing of permeability for all gases studied. Thermal rearrangement of the MMMs resulted in higher permeabilities but lower selectivities. For all the membranes and gas pairs studied, the balance of permeability vs. selectivity surpassed the 1991 Robeson's upper bound, and approached or even exceeded the 2008 line, for MMMs having 30% PPN loads. In all cases, the HPA-MMMs before thermal rearrangement provided good selectivity versus permeability compromise, similar to their thermally rearranged counterparts but in the zone of high selectivity. For H2/CH4, H2/N2, these nonthermally rearranged MMMs approach the 2008 Robeson's upper bound while H2/CO2 gives selective transport favoring H2 on the 1991 Robeson's bound. Thus, attending to the energy cost of thermal rearrangement, it could be avoided in some cases especially when high selectivity is the target rather than high permeability.

Selective Separation of Lithium Chloride by Organogels Containing Strapped Calix[4]pyrroles.

Reported herein are two functionalized crown ether strapped calix[4]pyrroles, H1 and H2. As inferred from competitive salt binding experiments carried out in nitrobenzene-d5 and acetonitrile-d3, these hosts capture LiCl selectively over four other test salts, viz. NaCl, KCl, MgCl2, and CaCl2. Support for the selectivity came from density functional theory (DFT) calculations carried out in a solvent continuum. These theoretical analyses revealed a higher innate affinity for LiCl in the case of H1, but a greater selectivity relative to NaCl in the case of H2, recapitulating that observed experimentally. Receptors H1 and H2 were outfitted with methacrylate handles and subject to copolymerization with acrylate monomers and cross-linkers to yield gels, G1 and G2, respectively. These two gels were found to adsorb lithium chloride preferentially from an acetonitrile solution containing a mixture of LiCl, NaCl, KCl, MgCl2, and CaCl2 and then release the lithium chloride in methanol. The gels could then be recycled for reuse in the selective adsorption of LiCl. As such, the present study highlights the use of solvent polarity switching to drive separations with potential applications in lithium purification and recycling.

Engineering Li/Na selectivity in 12-Crown-4-functionalized polymer membranes.

Lithium is widely used in contemporary energy applications, but its isolation from natural reserves is plagued by time-consuming and costly processes. While polymer membranes could, in principle, circumvent these challenges by efficiently extracting lithium from aqueous solutions, they usually exhibit poor ion-specific selectivity. Toward this end, we have incorporated host-guest interactions into a tunable polynorbornene network by copolymerizing 1) 12-crown-4 ligands to impart ion selectivity, 2) poly(ethylene oxide) side chains to control water content, and 3) a crosslinker to form robust solids at room temperature. Single salt transport measurements indicate these materials exhibit unprecedented reverse permeability selectivity (∼2.3) for LiCl over NaCl-the highest documented to date for a dense, water-swollen polymer. As demonstrated by molecular dynamics simulations, this behavior originates from the ability of 12-crown-4 to bind Na+ ions more strongly than Li+ in an aqueous environment, which reduces Na+ mobility (relative to Li+) and offsets the increase in Na+ solubility due to binding with crown ethers. Under mixed salt conditions, 12-crown-4 functionalized membranes showed identical solubility selectivity relative to single salt conditions; however, the permeability and diffusivity selectivity of LiCl over NaCl decreased, presumably due to flux coupling. These results reveal insights for designing advanced membranes with solute-specific selectivity by utilizing host-guest interactions.

Gas Separation by Mixed Matrix Membranes with Porous Organic Polymer Inclusions within o-Hydroxypolyamides Containing m-Terphenyl Moieties.

A hydroxypolyamide (HPA) manufactured from 2,2-bis(3-amino-4-hydroxy phenyl)-hexafluoropropane (APAF) diamine and 5'-terbutyl-m-terphenyl-4,4''-dicarboxylic acid chloride (tBTpCl), and a copolyimide produced by stochiometric copolymerization of APAF and 4,4'-(hexafluoroisopropylidene) diamine (6FpDA), using the same diacid chloride, were obtained and used as polymeric matrixes in mixed matrix membranes (MMMs) loaded with 20% (w/w) of two porous polymer networks (triptycene-isatin, PPN-1, and triptycene-trifluoroacetophenone, PPN-2). These MMMs, and also the thermally rearranged membranes (TR-MMMs) that underwent a thermal treatment process to convert the o-hydroxypolyamide moieties to polybenzoxazole ones, were characterized, and their gas separation properties evaluated for H2, N2, O2, CH4, and CO2. Both TR process and the addition of PPN increased permeability with minor decreases in selectivity for all gases tested. Excellent results were obtained, in terms of the permeability versus selectivity compromise, for H2/CH4 and H2/N2 separations with membranes approaching the 2008 Robeson's trade-off line. The best gas separation properties were obtained when PPN-2 was used. Finally, gas permeation was characterized in terms of chain intersegmental distance and fraction of free volume of the membrane along with the kinetic diameters of the permeated gases. The intersegmental distance increased after TR and/or the addition of PPN-2. Permeability followed an exponential dependence with free volume and a quadratic function of the kinetic diameter of the gas.

Origins of Lithium/Sodium Reverse Permeability Selectivity in 12-Crown-4-Functionalized Polymer Membranes.

Direct lithium extraction via membrane separations has been fundamentally limited by lack of monovalent ion selectivity exhibited by conventional polymeric membranes, particularly between sodium and lithium ions. Recently, a 12-Crown-4-functionalized polynorbornene membrane was shown to have the largest lithium/sodium permeability selectivity observed in a fully aqueous system to date. Using atomistic molecular dynamics simulations, we reveal that this selectivity is due to strong interactions between sodium ions and 12-Crown-4 moieties, which reduce sodium ion diffusivity while leaving lithium ion mobility relatively unaffected. Moreover, the ion diffusivities in the membrane, when scaled by their respective solution diffusivities and free ion fractions, can be collapsed to an almost universal relationship depending on solvent volume fraction.

Prediction of lattice energy of benzene crystals: A robust theoretical approach.

We present an inexpensive and robust theoretical approach based on the fragment molecular orbital methodology and the spin-ratio scaled second-order Møller-Plesset perturbation theory to predict the lattice energy of benzene crystals within 2 kJ⋅mol-1 . Inspired by the Harrison method to estimate the Madelung constant, the proposed approach calculates the lattice energy as a sum of two- and three-body interaction energies between a reference molecule and the surrounding molecules arranged in a sphere. The lattice energy converges rapidly at a radius of 13 Å. Adding the corrections to account for a higher correlated level of theory and basis set superposition for the Hartree Fock (HF) level produced a lattice energy of -57.5 kJ⋅mol-1 for the benzene crystal structure at 138 K. This estimate is within 1.6 kJ⋅mol-1 off the best theoretical prediction of -55.9 kJ⋅mol-1 . We applied this approach to calculate lattice energies of the crystal structures of phase I and phase II-polymorphs of benzene-observed at a higher temperature of 295 K. The stability of these polymorphs was correctly predicted, with phase II being energetically preferred by 3.7 kJ⋅mol-1 over phase I. The proposed approach gives a tremendous potential to predict stability of other molecular crystal polymorphs.