Probing the performance of DFT in the structural characterization of [FeFe] hydrogenase models.

In this work, a series of DFT and DFT-D methods is combined with double-ζ basis sets to benchmark their performance in predicting the structures of five newly synthesized hexacarbonyl diiron complexes with a bridging ligand featuring a µ-S2C3 motif in a ring-containing unit functionalized with aromatic groups. Such complexes have been considered as [FeFe] hydrogenase catalytic site models with potential for eco-friendly energetic applications. According to this assessment, r2SCAN is identified as the density functional recommended for the reliable description of the molecular and crystal structures of the herein studied models. However, the butterfly (µ-S)2Fe2 core of the models demonstrates a minor deformation of its optimized geometry obtained from both molecular and periodic calculations. The FeFe bond length is slightly underestimated while the FeS bonds tend to be too long. Adding the D3(BJ) correction to r2SCAN does not lead to any improvement in the calculated structures.

Photoactuators Based on Plastically Flexible α-Cyanostilbene Molecular Crystals Driven by the Solid-State [2+2] Cycloaddition Reaction.

Organic photomechanical molecular crystals are promising candidates for photoactuators, which have potential applications as smart materials in various fields. However, it is still challenging to fabricate photomechanical molecular crystals with flexibility because most of the molecular crystals are brittle and the mechanism of flexible crystals remains controversial. Here, a plastically flexible α-cyanostilbene crystal has been synthesized that can undergo solid-state [2+2] cycloaddition reaction under violet or UV irradiation and exhibits excellent photomechanical bending properties. A hook-shaped crystal can lift 0.7 mg object upward by 1.5 cm, which proves its potential for application as photoactuators. When complex with the agarose polymer, the molecules will be in the form of macroscopic crystals, which can drive the composite films to exhibit excellent photomechanical bending performance. Upon irradiation with UV light, the composite film can quickly lift 18.0 mg object upward by 0.3 cm. The results of this work may facilitate the application of macroscale crystals as photoactuators.

Stereo-Active Lone Pairs Induced Second Harmonic Generation Responses and Electrocatalytic Activity in Hybrid Material.

The lone pair electrons in the electronic structure of molecules have been a prominent research focus in chemistry for more than a century. Stable s2 lone pair electrons significantly influence material properties, including thermoelectric properties, nonlinear optical properties, ferroelectricity, and electro(photo)catalysis. While major advances have been achieved in understanding the influence of lone pair electrons on material characteristics, research on this effect in organic-inorganic hybrid materials is in its initial stage. In this work, we successfully obtained a novel organic-inorganic hybrid multifunctional material incorporating Ge with 4s2 lone pair electrons, (MeHDabco)2[GeBr3]4-H2O (MeHDabco=N-methyl-1,4-diazabicyclo[2.2.2]octane) (1). Driven by the stereochemically active lone pair electrons on the Ge2+, 1 crystallizes in the noncentrosymmetric space group P21 at room temperature and exhibits good second harmonic generation (SHG) responses. Interestingly, 1 also shows electrocatalytic activity for the hydrogen evolution reaction (HER) due to the existence of lone pair electrons on Ge2+ cations. The electrochemical experiment combined with the density functional theory (DFT) calculations revealed that the lone pair electrons act as both an active site for proton adsorption and facilitate the ionization of water. This work not only emphasizes the important role of lone pair electrons in material properties and functions but also provides new insight for designing novel Ge-based multifunctional hybrid materials.

Benzodifuranone Crystals with Solid-State Emission Arising from the Introduction of Bulky Substituents.

Benzodifuranone derivatives, structural analogues of widely studied diketopyrrolopyrrole, have been reported to show photoluminescence exclusively in solution states. Here, upon introducing bulky aromatic groups into benzodifuranone dyes, crystals with unprecedented solid-state emission were obtained. A crystallographic analysis revealed that the emissive properties should most likely be attributed to the absence of stacking between the dye scaffolds. In addition to the solid-state emission, the compound showed responsivity to external stimuli, i. e., luminescent mechanochromism and thermosalient effect.

Transition Behaviors of Isostructural Hydrogen-Bonded Frameworks Composed of Naphthalene, Quinoxaline, and Pyrazinopyrazine Derivatives.

A series of isostructural reticular frameworks with systematic differences on chemical structures allows us to disclose correlations between specific structural factors and properties, providing insights for designing novel porous materials. However, even slight differences in the molecular structure often lead to non-isostructural polymorphic frameworks particularly in the case of hydrogen-bonded organic frameworks (HOFs) because the structures of HOFs are based on a subtle balance of reversible interactions. In this study, we found that three simple analogues of tetracarboxylic acids with naphthalene, quinoxaline, and pyrazinopyrazine cores (NT, QX, and PP, respectively) yielded isostructural solvated HOFs (NT-1, QX-1, and PP-1, respectively), where hydrogen-bonded sql-networked sheets were slip-stacked with closely similar manners. More importantly, these isostructural HOFs underwent structural transformations in different manners upon removal of the guest solvents. Comparison of the crystal structures of the HOFs before and after the transformation revealed that intermolecular interactions of the core significantly affected on rearrangements of hydrogen bonds in the transformation. The results suggest the potential to control the properties and functions of isostructural HOFs by elements in the core.

Finite Temperature String with Order Parameter as Collective Variables for Molecular Crystal: A Case of Polymorphic Transformation of TNT under External Electric Field.

An external electric field is an effective tool to induce the polymorphic transformation of molecular crystals, which is important practically in the chemical, material, and energy storage industries. However, the understanding of this mechanism is poor at the molecular level. In this work, two types of order parameters (OPs) were constructed for the molecular crystal based on the intermolecular distance, bond orientation, and molecular orientation. Using the K-means clustering algorithm for the sampling of OPs based on the Euclidean distance and density weight, the polymorphic transformation of TNT was investigated using a finite temperature string (FTS) under external electric fields. The potential of mean force (PMF) was obtained, and the essence of the polymorphic transformation between o-TNT and m-TNT was revealed, which verified the effectiveness of the FTS method based on K-means clustering to OPs. The differences in PMFs between the o-TNT and transition state were decreased under external electric fields in comparison with those in no field. The fields parallel to the c-axis obviously affected the difference in PMF, and the relationship between the changes in PMFs and field strengths was found. Although the external electric field did not promote the convergence, the time of the polymorphic transformation was reduced under the external electric field in comparison to its absence. Moreover, under the external electric field, the polymorphic transformation from o-TNT to m-TNT occurred while that from m-TNT to o-TNT was prevented, which was explained by the dipole moment of molecule, relative permittivity, chemical potential difference, nucleation work and nucleation rate. This confirmed that the polymorphic transformation orientation of the molecular crystal could be controlled by the external electric field. This work provides an effective way to explore the polymorphic transformation of the molecular crystals at a molecular level, and it is useful to control the production process and improve the performance of energetic materials by using the external electric fields.

Molecular Crystals Constructed by Polar Molecular Cages: A Promising System for Exploring High-performance Infrared Nonlinear Optical Crystals.

Polar molecular crystals, with their densely stacked polar nonlinear optical (NLO) active units, are favored for their large second harmonic generation (SHG) responses and birefringence. However, their potential for practical applications as Infrared (IR) NLO materials has historically been underappreciated due to the weak inter-molecular interaction forces that may compromise their physicochemical properties. In this study, we propose molecular crystals with polar molecular cages as a treasure-house for the development of superior IR NLO materials and a representative system, binary chalcogenide molecular crystals, composed of [P4 Sn ] (n=3-9) polar molecular cages, is introduced. These crystals may not only achieve wide band gap, large SHG response, and birefringence in a single structure, but also exhibit favorable physicochemical properties. We subsequently obtained a polar molecular crystal, α-P4 S5 , which demonstrated exceptional IR optical properties, including a strong SHG response (1.1×AGS), wide band gap (3.02 eV), large birefringence (0.134@2050 nm), and a broad transmission range (0.41-14.7 µm). Moreover, it showed excellent water resistance and hardness. These findings highlight the potential of polar molecular crystals as a promising platform for the development of high-performance IR NLO materials.

Two-Dimensional Nonbenzenoid Heteroacene Crystals Synthesized via In-Situ Embedding of Ladder Bipyrazinylenes on Au(111).

Precisely introducing topological defects is an important strategy in nanographene crystal engineering because defects can tune π-electronic structures and control molecular assemblies. The synergistic control of the synthesis and assembly of nanographenes by embedding the topological defects to afford two-dimensional (2D) crystals on surfaces is still a great challenge. By in-situ embedding ladder bipyrazinylene (LBPy) into acene, the narrowest nanographene with zigzag edges, we have achieved the precise preparation of 2D nonbenzenoid heteroacene crystals on Au(111). Through intramolecular electrocyclization of o-diisocyanides and Au adatom-directed [2+2] cycloaddition, the nonbenzenoid heteroacene products are produced with high chemoselectivity, and lead to the molecular 2D assembly via LBPy-derived interlocking hydrogen bonds. Using bond-resolved scanning tunneling microscopy, we determined the atomic structures of the nonbenzenoid heteroacene product and diverse organometallic intermediates. The tunneling spectroscopy measurements revealed the electronic structure of the nonbenzenoid heteroacene, which is supported by density functional theory (DFT) calculations. The observed distinct organometallic intermediates during progression annealing combined with DFT calculations demonstrated that LBPy formation proceeds via electrocyclization of o-diisocyanides, trapping of heteroarynes by Au adatoms, and stepwise elimination of Au adatoms.

Discriminating and understanding molecular crystal polymorphism.

Polymorph discrimination for a molecular crystal has long been a challenging task, which, nonetheless, is a major concern in the pharmaceutical industry. In this work, we have investigated polymorph discrimination on three different molecular crystals, tetrolic acid, oxalic acid, and oxalyl dihydrazide, covering both packing polymorphism and conformational polymorphism. To gain more understanding, we have performed energy decomposition analysis based on many-body expansion, and have compared the results from the XO-PBC method, that is, the eXtended ONIOM method (XO) with the periodic boundary condition (PBC), with those from some commonly used dispersion corrected density functional theory (DFT-D) methods. It is shown here that, with the XYG3 doubly hybrid functional chosen as the target high level to capture the intra- and short-range intermolecular interactions, and the periodic PBE as the basic low level to take long range interactions into account, the XO-PBC(XYG3:PBE) method not only obtains the correct experimental stability orderings, but also predicts reasonable polymorph energy ranges for all three cases. Our results have demonstrated the usefulness of the present theoretical methods, in particular XO-PBC, while highlighted the importance of a better treatment of different kinds of interactions to be beneficial to polymorph control.

The Factor beyond Schmidt's Criteria Impacting the Photo-Induced [2+2] Cycloaddition Reactivity and Photoactuation of Molecular Crystals Based on Cyclic Chalcone Analogues.

Generally, the potential reactive "olefin pairs" in the molecular crystals satisfying Schmidt's criteria could undergo topological [2+2] cycloaddition. In this study, another factor that affects the photodimerization reactivity of chalcone analogues was found. The cyclic chalcone analogues of (E)-2-(2,4-dichlorobenzylidene)-2,3-dihydro-1H-inden-1-one (BIO), (E)-2-(naphthalen-2-ylmethylene)-2,3-dihydro-1H-inden-1-one (NIO), (Z)-2-(2,4-dichlorobenzylidene)benzofuran-3(2H)-one (BFO), and (Z)-2-(2,4-dichlorobenzylidene)benzo[b]thiophen-3(2H)-one (BTO) have been synthesized. While the geometrical parameters for the molecular packing of the above four compounds did not exceed Schmidt's criteria, [2+2] cycloaddition did not occur in the crystals of BIO and BTO. The single crystal structures and Hirshfeld surface analyses revealed that interactions of C=O⋅⋅⋅H (CH2 ) existed between adjacent molecules in the crystal of BIO. Therefore, the carbonyl and methylene groups linked with one carbon atom in carbon-carbon double bond were tightly confined in the lattice, acting as a tweezer to inhibit free movement of the double bond and suppressing [2+2] cycloaddition. In the crystal of BTO, similar interactions of Cl⋅⋅⋅S and C=O⋅⋅⋅H (C6 H4 ) prevented free movement of the double bond. In contrast, the intermolecular interaction of C=O⋅⋅⋅H only exists around the carbonyl group in the crystals of BFO and NIO, leaving the C=C double bonds to move freely and allowing the occurrence of [2+2] cycloaddition. Driven by photodimerization, the needle-like crystals of BFO and NIO displayed evident photo-induced bending behavior. This work demonstrates that the intermolecular interactions around carbon-carbon double bond affect the [2+2] cycloaddition reactivity except for Schmidt's criteria. These findings provide valuable insights into the design of photomechanical molecular crystalline materials.

The Ambiguous Origin of Thermochromism in Molecular Crystals of Dichalcogenides: Chalcogen Bonds versus Dynamic Se-Se/Te-Te Bonds.

We report thermochromism in crystals of diphenyl diselenide (dpdSe) and diphenyl ditelluride (dpdTe), which is at variance with the commonly known mechanisms of thermochromism in molecular crystals. Variable temperature neutron diffraction studies indicated no conformational change, tautomerization or phase transition between 100 K and 295 K. High-pressure crystallography studies indicated no associated piezochromism in dpdSe and dpdTe crystals. The evolution of the crystal structures and their electronic band structure with pressure and temperature reveal the contributions of intramolecular and intermolecular factors towards the origin of thermochromism-especially the intermolecular Se⋅⋅⋅Se and Te⋅⋅⋅Te chalcogen bonds and torsional modes of vibrations around the dynamic Se-Se and Te-Te bonds. Further, a co-crystal of dpdSe with iodine (dpdSe-I2 ) and an alloy crystal of dpdSe and dpdTe implied a predominantly intramolecular origin of the observed thermochromism associated with vibronic coupling.

An Ultrathin Inorganic Molecular Crystal Interfacial Layer for Stable Zn Anode.

Zn metal anode suffers from dendrite growth and side reactions during cycling, significantly deteriorating the lifespan of aqueous Zn metal batteries. Herein, we introduced an ultrathin and ultra-flat Sb2 O3 molecular crystal layer to stabilize Zn anode. The in situ optical and atomic force microscopes observations show that such a 10 nm Sb2 O3 thin layer could ensure uniform under-layer Zn deposition with suppressed tip growth effect, while the traditional WO3 layer undergoes an uncontrolled up-layer Zn deposition. The superior regulation capability is attributed to the good electronic-blocking ability and low Zn affinity of the molecular crystal layer, free of dangling bonds. Electrochemical tests exhibit Sb2 O3 layer can significantly improve the cycle life of Zn anode from 72 h to 2800 h, in contrast to the 900 h of much thicker WO3 even in 100 nm. This research opens up the application of inorganic molecular crystals as the interfacial layer of Zn anode.

Inorganic Low k Cage-molecular Crystals.

For the interlayer dielectric in microelectronics, light element compounds are preferably accepted due to less electronic polarization. Here, the nontrivial dielectric nature of the Sb4O6 cage-molecular crystal, known as α-antimony trioxide (α-Sb2O3), is reported. The gas-phase synthesized α-Sb2O3 nanoflakes are of high crystal quality, from which the abnormal local admittance responses were revealed by scanning microwave impedance microscopy (sMIM). The remarkably low dielectric constant (k), 2.0∼2.5, is corroborated by the analysis of the thickness-dependent sMIM-capacitance signal. In light of the theoretical calculations, the ultralow molecular density and the significantly suppressed ionic polarization are both crucial to the highly reduced k. Combining with the excellent optical band gap, thermal stability, and breakdown strength, α-Sb2O3 is a promising low k dielectric.

Molecular Forcefield Methods for Describing Energetic Molecular Crystals: A Review.

Energetic molecular crystals are widely applied for military and civilian purposes, and molecular forcefields (FF) are indispensable for treating the microscopic issues therein. This article reviews the three types of molecular FFs that are applied widely for describing energetic crystals-classic FFs, consistent FFs, and reactive FFs (ReaxFF). The basic principle of each type of FF is briefed and compared, with the application introduced, predicting polymorph, morphology, thermodynamics, vibration spectra, thermal property, mechanics, and reactivity. Finally, the advantages and disadvantages of these FFs are summarized, and some directions of future development are suggested.

Fast and accurate quantum Monte Carlo for molecular crystals.

Computer simulation plays a central role in modern-day materials science. The utility of a given computational approach depends largely on the balance it provides between accuracy and computational cost. Molecular crystals are a class of materials of great technological importance which are challenging for even the most sophisticated ab initio electronic structure theories to accurately describe. This is partly because they are held together by a balance of weak intermolecular forces but also because the primitive cells of molecular crystals are often substantially larger than those of atomic solids. Here, we demonstrate that diffusion quantum Monte Carlo (DMC) delivers subchemical accuracy for a diverse set of molecular crystals at a surprisingly moderate computational cost. As such, we anticipate that DMC can play an important role in understanding and predicting the properties of a large number of molecular crystals, including those built from relatively large molecules which are far beyond reach of other high-accuracy methods.

High Li-Ion Conductivity in Li{N(SO2F)2}(NCCH2CH2CN)2 Molecular Crystal.

There is an urgent need to develop solid electrolytes based on organic molecular crystals for application in energy devices. However, the quest for molecular crystals with high Li-ion conductivity is still in its infancy. In this study, the high Li-ion conductivity of a Li{N(SO2F)2}(NCCH2CH2CN)2 molecular crystal is reported. The crystal shows a Li-ion conductivity of 1 × 10-4 S cm-1 at 30 °C and 1 × 10-5 S cm-1 at -20 °C, with a low activation energy of 28 kJ mol-1. The conductivity at 30 °C is one of the highest values attainable by molecular crystals, whereas that at -20 °C is approximately 2 orders of magnitude higher than previously reported values. Furthermore, the all-solid-state Li-battery fabricated using this solid electrolyte demonstrates stable cycling, thereby maintaining 90% of the initial capacity after 100 charge-discharge cycles. The finding of high Li-ion conductivity in molecular crystals paves the way for their application in all-solid-state Li-batteries.

Effects of MEH-PPV Molecular Ordering in the Emitting Layer on the Luminescence Efficiency of Organic Light-Emitting Diodes.

We investigated the effects of molecular ordering on the electro-optical characteristics of organic light-emitting diodes (OLEDs) with an emission layer (EML) of poly[2-methoxy-5-(2-ethylhexyloxy)-1,4-phenylenevinylene] (MEH-PPV). The EML was fabricated by a solution process which can make molecules ordered. The performance of the OLED devices with the molecular ordering method was compared to that obtained through fabrication by a conventional spin coating method. The turn-on voltage and the luminance of the conventional OLEDs were 5 V and 34.75 cd/m2, whereas those of the proposed OLEDs were 4.5 V and 120.3 cd/m2, respectively. The underlying mechanism of the higher efficiency with ordered molecules was observed by analyzing the properties of the EML layer using AFM, SE, XRD, and an LCR meter. We confirmed that the electrical properties of the organic thin film can be improved by controlling the molecular ordering of the EML, which plays an important role in the electrical characteristics of the OLED.

Vibrational Dynamics of Crystalline 4-Phenylbenzaldehyde from INS Spectra and Periodic DFT Calculations.

The present work emphasizes the value of periodic density functional theory (DFT) calculations in the assessment of the vibrational spectra of molecular crystals. Periodic calculations provide a nearly one-to-one match between the calculated and observed bands in the inelastic neutron scattering (INS) spectrum of crystalline 4-phenylbenzaldehyde, thus validating their assignment and correcting previous reports based on single molecule calculations. The calculations allow the unambiguous assignment of the phenyl torsional mode at ca. 118-128 cm-1, from which a phenyl torsional barrier of ca. 4000 cm-1 is derived, and the identification of the collective mode involving the antitranslational motion of CH···O bonded pairs, a hallmark vibrational mode of systems where C-H···O contacts are an important feature.

In Situ Assessment of Intrinsic Strength of X-I⋯OA-Type Halogen Bonds in Molecular Crystals with Periodic Local Vibrational Mode Theory.

Periodic local vibrational modes were calculated with the rev-vdW-DF2 density functional to quantify the intrinsic strength of the X-I⋯OA-type halogen bonding (X = I or Cl; OA: carbonyl, ether and N-oxide groups) in 32 model systems originating from 20 molecular crystals. We found that the halogen bonding between the donor dihalogen X-I and the wide collection of acceptor molecules OA features considerable variations of the local stretching force constants (0.1-0.8 mdyn/Å) for I⋯O halogen bonds, demonstrating its powerful tunability in bond strength. Strong correlations between bond length and local stretching force constant were observed in crystals for both the donor X-I bonds and I⋯O halogen bonds, extending for the first time the generalized Badger's rule to crystals. It is demonstrated that the halogen atom X controlling the electrostatic attraction between the σ -hole on atom I and the acceptor atom O dominates the intrinsic strength of I⋯O halogen bonds. Different oxygen-containing acceptor molecules OA and even subtle changes induced by substituents can tweak the n → σ ∗ (X-I) charge transfer character, which is the second important factor determining the I⋯O bond strength. In addition, the presence of the second halogen bond with atom X of the donor X-I bond in crystals can substantially weaken the target I⋯O halogen bond. In summary, this study performing the in situ measurement of halogen bonding strength in crystalline structures demonstrates the vast potential of the periodic local vibrational mode theory for characterizing and understanding non-covalent interactions in materials.

Toward Developing a Predictive Approach To Assess Electron Beam Instability during Transmission Electron Microscopy of Drug Molecules.

During drug development control of polymorphism, particle properties and impurities are critical for ensuring a good quality, reproducible, and safe medicine. A wide variety of analytical techniques are employed in demonstrating the regulators control over the drug substance and product manufacturing, storage, and supply. Transmission electron microscopy (TEM) offers the opportunity to analyze in detail pharmaceutical systems at a length scale and limit of detection not readily achieved by many traditional techniques. However, the use of TEM as a characterization tool for drug development is uncommon due to possible damage caused by the electron beam. This work outlines the development of a model, using molecular descriptors, to predict the electron beam stability of active pharmaceutical ingredients (API). For a given set of conditions and a particular imaging or analytical mode, the total number of electrons per unit area, which causes observable damage to a sample in the TEM, can be defined as the critical fluence ( CF). Here the CF of 20 poorly water-soluble APIs were measured using selected area electron diffraction. Principal component analysis was used to select the most influential molecular descriptors on CF, which were shown to be descriptors involving the degree of conjugation, the number of hydrogen bond donors and acceptors, and the number of rotatable bonds. These were used to generate several multiple linear regression models. The model that provided the best fit to the measured CF data included the ratio of the number of conjugated carbons to nonconjugated carbons, the ratio of the number of hydrogen bond donors to acceptors, and the ratio of the number of hydrogen bond acceptors to donors. Using this model, the CF of the majority of the compounds was predicted within ±2 e-/Å2. Molecules with no hydrogen bond acceptors did not fit the model accurately possibly due to the limited sample size or the influence of other parameters not included in this model, such as intermolecular bond energies. The model presented can be used to support pharmaceutical development by quickly assessing the stability of other poorly soluble drugs in TEM. Provided that the model suggests that the API is relatively stable to electron irradiation, TEM offers the prospect of determining the presence of crystalline material at low levels at length scales and limits of detection unobtainable by other techniques. This is particularly so for amorphous solid dispersions.