Molecular Crystals Research Articles

Wigner molecular crystals from multielectron moiré artificial atoms.

Semiconductor moiré superlattices provide a versatile platform to engineer quantum solids composed of artificial atoms on moiré sites. Previous studies have mostly focused on the simplest correlated quantum solid-the Fermi-Hubbard model-in which intra-atom interactions are simplified to a single onsite repulsion energy U. Here we report the experimental observation of Wigner molecular crystals emerging from multielectron artificial atoms in twisted bilayer tungsten disulfide moiré superlattices. Using scanning tunneling microscopy, we demonstrate that Wigner molecules appear in multielectron artificial atoms when Coulomb interactions dominate. The array of Wigner molecules observed in a moiré superlattice comprises a crystalline phase of electrons: the Wigner molecular crystal, which is shown to be highly tunable through mechanical strain, moiré period, and carrier charge type.

Surface Engineering of the Mechanical Properties of Molecular Crystals via an Atomistic Model for Computing the Facet Stress Response of Solids.

Invited for the cover of this issue are Mubarak Almehairbi, Vikram C. Joshi, Changquan Calvin Sun and Sharmarke Mohamed. The image depicts the digital exploration of the mechanical properties of crystals on specific facets that may be of interest for materials applications by "dialing-in" their stress response. Read the full text of the article at 10.1002/chem.202400779.

Exfoliatable Layered 2D Honeycomb Crystals of Host‐guest Complexes Networked by CH‐π Hydrogen Bonds

AbstractStudies of graphene show that robust chemical bonds such as covalent bonds with trigonal‐planar atoms afford layered atomic 2D crystals possessing unique properties. Although layered molecular crystals are of interest to diversify elements and structures of 2D materials, the structural diversity of molecules as well as weak intermolecular interactions inevitably makes the design to be one‐off and individual. We herein report a versatile method to assemble layered molecular crystals. By developing a D3‐symmetry host at vertices to form a honeycomb layer, a diverse range of layered 2D host–guest crystals were obtained. Substituents on the host, elements/structures of the guest, the stereochemistry of the host and types of intercalants were diversified, which should allow for 6×32×3×2 combinations for structural diversification.

Exfoliatable Layered 2D Honeycomb Crystals of Host-guest Complexes Networked by CH-π Hydrogen Bonds.

Studies of graphene show that robust chemical bonds such as covalent bonds with trigonal-planar atoms afford layered atomic 2D crystals possessing unique properties. Although layered molecular crystals are of interest to diversify elements and structures of 2D materials, the structural diversity of molecules as well as weak intermolecular interactions inevitably makes the design to be one-off and individual. We herein report a versatile method to assemble layered molecular crystals. By developing a D3-symmetry host at vertices to form a honeycomb layer, a diverse range of layered 2D host-guest crystals were obtained. Substituents on the host, elements/structures of the guest, the stereochemistry of the host and types of intercalants were diversified, which should allow for 6×32×3×2 combinations for structural diversification.

Machine Learning Classification of Local Environments in Molecular Crystals.

Identifying local structural motifs and packing patterns of molecular solids is a challenging task for both simulation and experiment. We demonstrate two novel approaches to characterize local environments in different polymorphs of molecular crystals using learning models that employ either flexibly learned or handcrafted molecular representations. In the first case, we follow our earlier work on graph learning in molecular crystals, deploying an atomistic graph convolutional network combined with molecule-wise aggregation to enable per-molecule environmental classification. For the second model, we develop a new set of descriptors based on symmetry functions combined with a point-vector representation of the molecules, encoding information about the positions and relative orientations of the molecule. We demonstrate very high classification accuracy for both approaches on urea and nicotinamide crystal polymorphs and practical applications to the analysis of dynamical trajectory data for nanocrystals and solid-solid interfaces. Both architectures are applicable to a wide range of molecules and diverse topologies, providing an essential step in the exploration of complex condensed matter phenomena.

Dual responsive behaviors in fluorescence of molecular crystals based on naphthalene pyridyl derivatives

Three compounds based on naphthalene and pyridine structures with different nitrogen positions (NP-2, NP-3, and NP-4) were synthesized. Except for displaying aggregation-induced emission (AIE) properties, the solution and crystals of these compounds show reversible fluorescence changes under alternative acid and base treatments. Experimental results, including X-ray diffraction, optical spectroscopy measurements, and cross-polarized microscope characterization, suggested that compounds might undergo reversible crystal-to-crystal transition upon external acid and base stimuli. It also revealed that the position of the heteroatom nitrogen in the pyridine unit significantly affected the intra- and intermolecular interactions, leading to restricted-intermolecular motions (RIMs) that caused the AIE effect. Besides, computational calculations indicated significant changes in molecular geometries and electronic interactions between the initial states and the protonated states. An information encryption system involving activation and operation codes was established based on the multiple stimuli-responsive properties. Our work not only presents a facile way to develop multiple responsive properties in molecular crystal systems but also explores their potential as highly emissive functional dyes for future applications.

3D Nanocrystallography and the Imperfect Molecular Lattice.

Crystallographic analysis relies on the scattering of quanta from arrays of atoms that populate a repeating lattice. While large crystals built of lattices that appear ideal are sought after by crystallographers, imperfections are the norm for molecular crystals. Additionally, advanced X-ray and electron diffraction techniques, used for crystallography, have opened the possibility of interrogating micro- and nanoscale crystals, with edges only millions or even thousands of molecules long. These crystals exist in a size regime that approximates the lower bounds for traditional models of crystal nonuniformity and imperfection. Accordingly, data generated by diffraction from both X-rays and electrons show increased complexity and are more challenging to conventionally model. New approaches in serial crystallography and spatially resolved electron diffraction mapping are changing this paradigm by better accounting for variability within and between crystals. The intersection of these methods presents an opportunity for a more comprehensive understanding of the structure and properties of nanocrystalline materials.

Insight into material behavior via surface free energy calculations for common energetic materials

AbstractThe Gibbs Free energy is a driving force for equilibrium crystal shapes and the formation of crystal facets in molecular crystals. Orientation dependence of interfacial properties is linked to surface free energy (SFE). Prediction of orientation‐dependent properties such as thermal stability, mechanical response, and compatibility with binders require a systematic approach to the quantification of SFE. In molecular crystals, entropy has a much larger contribution among all ordered crystalline materials. In this paper, we extend our previously developed method to quantify SFE and entropy of β‐HMX to other common energetic materials–TATB, α‐RDX, and PETN. Two complimentary approaches, Nonequilibrium Thermodynamic Integration (NETI) and Steered Molecular Dynamics (SMD) methods are used to obtain insight into interfacial phenomena along with surface free energy estimates. We discuss the relevance of surface free energy and the importance of surface entropy for facetted molecular crystals in understanding crystal properties, activation of slip planes, and potential pathways for fracture. These values allow us to predict theoretical crystal shape using Wulff Construction, better understand the effect of hydrogen bonding on SFE, and the diversity of bonding environment in energetic crystals. In particular, in crystals with low stacking fault energy, the SMD values can be inconclusive due to the triggering of slip plane motions. In cases where SMD simulations lead to large deformations and high uncertainty, the NETI approach can still provide SFE estimates.

Guest Editorial: Special Issue on Multi-Functional Molecular Liquids, Solids, and Crystals for Materials Engineering & Manufacturing — Part 2

Guest Editorial: Special Issue on Multi-Functional Molecular Liquids, Solids, and Crystals for Materials Engineering & Manufacturing — Part 2

Social Self-Sorting of Quasi-Racemates: A Unique Approach for Dual-Pore Molecular Crystals.

Despite recent advances in porous organic molecular crystals, the engineering of dual-pore systems within the intermolecular voids remains a significant challenge. In this study, we have achieved the crystallization-induced social self-sorting of "quasi-racemic" dialdehydes into a macrocyclic imine. X-ray crystallographic analysis unambiguously characterizes the resulting structure as incorporating two quasi-racemate pairs with four diamine molecules. Notably, different alkyl substituents on the quasi-racemates afford two types of one-dimensional pores within the macrocyclic imine crystal. The different adsorption properties of these pores were substantiated through adsorption experiments. An intriguing helical arrangement of guest molecules was observed within one of the pores. This study provides pioneering evidence that the social self-sorting of quasi-racemates offers a new methodology for creating dual-functional supramolecular materials.

First principles study on interactions in inorganic molecular crystals at zero dimensions

Inorganic molecular crystals (IMCs) are attracting significant scientific interest due to their distinctive anisotropic crystal structure. The molecular units in the form of cages are the fundamental unit of these materials and most of their properties are dependent on the geometry, their arrangement and the interactions present in between these cages. This work is carried out to shed light on the structural arrangement and the interaction properties of three IMCs including antimony oxide (Sb2O3), phosphorous tri-selenium (P4Se3) and phosphorous trioxide (P4O6) in zero dimensions. The reported calculations were carried out using density Functional Theory (DFT) to investigate the interactions between the cages based on ETS-NOCV, NBO, Wiberg bond indices and QTAIM analysis. The results revealed weak VdWs association along the inter-cages and significant covalent bonding within the cage. The investigation shows that Sb2O3 and P4Se3 IMCs have potential applications in high-pressure and stress-bearing materials.

Two-dimensional crystallization of precise side-chain giant molecules with constant building blocks ratio

This paper describes the two-dimensional (2D) crystallization of side-chain giant molecules (SCGMs) having the fixed number ratio of two types of functionalized building blocks with 1:1. SCGMs are regarded to a certain degree as size-amplified versions of universal synthetic polymers that are prepared by precisely connecting molecular nanoparticles (MNPs) to polymer chains. Our previous experimental results have shown that individually changing the number ratio of one of the building blocks would significantly affect the structural parameters of molecular crystals, mainly the longitudinal thickness. Herein, we systematically discussed the 2D crystallization of two categories of POSS-based SCGMs with constant building blocks ratio. Solution self-assembly of such SCGMs resulted in well-defined 2D nanosheets with similar structural configurations, showing great tolerance in the number of building blocks. Subsequently, we proposed a “sandwiched-type” mode to illustrate the bilayer 2D supramolecular framework, in which crystalline blocks pack in head-to-head manner. In addition, sequence isomers provide an ideal platform to further verify that the thickness of 2D nanocrystal is linked with the tethering density determined by the number ratio of building blocks. We believe that this study could be a better complement to our previous work and then make us to be more in-depth understand the crystallization mechanism of SCGMs.

Surfactant-Polymer Complexation and Competition on Drug Nanocrystal Surfaces Control Crystallinity.

Nanosizing drug crystals has emerged as a successful approach to enabling oral bioavailability, as increasing drug crystal surface area improves dissolution kinetics and effective solubility. Recently, bottom-up methods have been developed to directly assemble nanosized crystals by leveraging polymer and surfactant excipients during crystallization to control crystal size, morphology, and structure. However, while significant research has investigated how polymers and other single additives inhibit or promote crystallization in pharmaceutical systems, there is little work studying the mechanistic interactions of multiple excipients on drug crystal structure and the extent of crystallinity, which can influence formulation performance. This study explores how the structure and crystallinity of a model hydrophobic drug crystal, fenofibrate, change as a result of competitive interfacial chemisorption between common nonionic surfactants (polysorbate 80 and sorbitan monooleate) and a surface-active polymer excipient (methylcellulose). Classical molecular dynamics simulations highlight how key intermolecular interactions, including surfactant-polymer complexation and surfactant screening of the crystal surface, modify the resulting crystal structure. In parallel, experiments generating drug nanocrystals in hydrogel thin films validate that drug crystallinity increases with an increasing weight fraction of surfactant. Simulation results reveal a connection between accelerated dynamics in the bulk crystal and the experimentally measured extent of crystallinity. To our knowledge, these are the first simulations that directly characterize structural changes in a drug crystal as a result of excipient surface composition and relate the experimental extent of crystallinity to structural changes in the molecular crystal. Our approach provides a mechanistic understanding of crystallinity in nanocrystallization, which can expand the range of orally deliverable small molecule therapies.

Comprehensive investigation of 2-methylquinolinium-5-sulphosalicylate monohydrate: Synthesis, characterization, multifaceted analysis of molecular structure, optical, and nonlinear optical properties

We introduce the synthesis and characterization of a novel molecular salt crystal, 2-methylquinolinium-5-sulphosalicylate monohydrate (MQSSA), obtained through an equimolar reaction between donors (D) and acceptors (A) employing crystal engineering principles. The resulting proton transfer complex, MQSSA, crystallized into a monoclinic crystal system with space group P21/C, as revealed by SCXRD analysis. Significant peaks observed in PXRD analysis confirmed the crystallinity and purity of the material. FT-IR and FT-Raman spectroscopy were employed to analyze various functional groups present in the molecule. UV–Vis-NIR spectroscopy provided insights into absorbance transmittance, lower cutoff wavelength, and optical band gap, suggesting it's suitability for diverse optical applications, while characteristic emission was identified using PL measurements. Thermal stability was assessed via TG-DTA and TG-DSC techniques. Nonlinear absorption coefficient (β), nonlinear refractive index (n2), and third-order nonlinear susceptibility (χ3) were quantified using the Z-scan technique. Hirshfeld surface analysis and fingerprint plots were utilized to gain insight into intermolecular interactions and atomic arrangements within the MQSSA crystal, emphasizing the roles of crystal structure in H...H and O...H/H...O interactions. Additionally, quantum computational methods were employed to evaluate various molecular-level electronic characteristics, including molecular orbitals, molecular electrostatic potential maps, and first and second-order nonlinear optical (NLO) polarizability, providing insights into the molecular-level NLO response properties.

Solar overall water-splitting by a spin-hybrid all-organic semiconductor

Direct solar-to-hydrogen conversion from pure water using all-organic heterogeneous catalysts remains elusive. The challenges are twofold: (i) full-band low-frequent photons in the solar spectrum cannot be harnessed into a unified S1 excited state for water-splitting based on the common Kasha-allowed S0 → S1 excitation; (ii) the H+ → H2 evolution suffers the high overpotential on pristine organic surfaces. Here, we report an organic molecular crystal nanobelt through the self-assembly of spin-one open-shell perylene diimide diradical anions (:PDI2-) and their tautomeric spin-zero closed-shell quinoid isomers (PDI2-). The self-assembled :PDI2-/PDI2- crystal nanobelt alters the spin-dependent excitation evolution, leading to spin-allowed S0S1 → 1(TT) → T1 + T1 singlet fission under visible-light (420 nm~700 nm) and a spin-forbidden S0 → T1 transition under near-infrared (700 nm~1100 nm) within spin-hybrid chromophores. With a triplet-triplet annihilation upconversion, a newly formed S1 excited state on the diradical-quinoid hybrid induces the H+ reduction through a favorable hydrophilic diradical-mediated electron transfer, which enables simultaneous H2 and O2 production from pure water with an average apparent quantum yield over 1.5% under the visible to near-infrared solar spectrum.

Porphyrin/Fullerene Porous Molecular Cocrystal Featuring a Robust One-Dimensional Channel.

Microporous molecular crystals are promising materials because of their designable porosity as well as their resistance to chemical and other stimuli. Here, we developed microporous molecular cocrystals by taking advantage of the specific interactions between porphyrins and fullerene molecules. Single-crystal X-ray diffraction analysis revealed that one nickel(II) porphyrin interacts with two fullerene molecules to form a two-dimensional honeycomb network with an eclipsed stacking mode, providing one-dimensional void channels. After the pores were activated by heat treatment or mechanical grinding, the prepared cocrystal can incorporate gas and solvent molecules reversibly while maintaining its single-crystallinity. Also, it retained its single-crystallinity in the presence of water, acid-base, or high pressure. These findings in this study expand the availability of molecular crystals based on intermolecular interactions as porous materials, which are expected to work under conditions that have not been applicable to other molecule-based porous materials.

Electronic Transport Properties of Molecular Clusters Sb4O6, P4Se3, and P4O6.

Inorganic molecular crystal (IMC) is a trending class of materials in which structural units comprise molecular cages or clusters bonded via van der Waal forces. The structure-property relationship in IMCs is less known due to the unusual assembly of molecular clusters in these materials. In this paper, the density functional theory-calculated electronic transport properties of the molecular clusters of antimony oxide (Sb4O6), phosphorus triselenide (P4Se3), and phosphorus trioxide (P4O6) are described in detail. The calculated values of highest occupied molecular orbital-lowest unoccupied molecular orbital gaps appeared as 5.487, 2.296, and 4.425 eV for Sb4O6, P4Se3, and P4O6, respectively. The work was carried out to explore the charge transport mechanism in IMCs in order to disclose their potential in practical applications. The calculations involved charge-transfer integral based on Marcus theory to compute the electronic coupling (V), reorganization energies (λ), and hopping rate (k) in the structures. The hopping rate for Sb4O6, P4Se3, and P4O6 is found as 8.49 × 10-12, 1.28 × 10-14, and 2.51 × 10-20 s-1, respectively. The transport properties of Sb4O6 are found better, which predicts the application of the relevant IMC for device grade applications. The findings of this study are important for future application of the IMCs in electronic and optoelectronic applications.

Suppression of Disorder in Benzamide and Thiobenzamide Crystals by Fluorine Substitution.

Disorder is a common feature of molecular crystals that complicates determination of structures and can potentially affect electric and mechanical properties. Suppression of disorder is observed in otherwise severely disordered benzamide and thiobenzamide crystals by substituting hydrogen with fluorine in the ortho-position of the phenyl ring. Fluorine occupancies of 20-30% are sufficient to suppress disorder without changing the packing motif. Crystal structure prediction calculations reveal a much denser lattice energy landscape for benzamide compared to 2-fluorobenzamide, suggesting that fluorine substitution makes disorder less likely.

Contributions beyond direct random-phase approximation in the binding energy of solid ethane, ethylene, and acetylene

The random-phase approximation (RPA) includes a subset of higher than second-order correlation-energy contributions, but stays in the same complexity class as the second-order Møller–Plesset perturbation theory (MP2) in both Gaussian-orbital and plane-wave codes. This makes RPA a promising ab initio electronic structure approach for the binding energies of molecular crystals. Still, some issues stand out in practical applications of RPA. Notably, compact clusters of nonpolar molecules are poorly described, and the interaction energies strongly depend on the reference single-determinant state. Using the many-body expansion of the binding energy of a crystal, we investigate those issues and the effect of beyond-RPA corrections. We find the beneficial effect of quartic-scaling exchange and non-ring coupled-cluster doubles corrections. The nonadditive interactions in compact trimers of molecules are improved by using the self-consistent Hartree–Fock orbitals instead of the usual Kohn–Sham states, but this kind of orbital input also leads to underestimated dimer energies. Overall, a substantial improvement over the RPA with a renormalized singles approach is possible at a modest quartic-scaling cost, which encourages further research into additional RPA corrections.

Predicting the Color Polymorphism of ROY from a Time-Dependent Optimally Tuned Screened Range-Separated Hybrid Functional.

Polymorphism is a well-known property of molecular crystals, which allows the same molecule to form solids with several crystalline structures that can differ significantly in physical properties. Polymorphs that possess different optical absorption properties in the visible range may exhibit different perceived colors, a phenomenon known as color polymorphism. One striking example of color polymorphism is given by 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecarbonitrile, known as ROY for its red-orange-yellow colors. First-principles prediction of color polymorphism may help in polymorph assignment and design but has proven to be challenging. Here, we predict the absorption spectra and simulate the colors of 12 ROY polymorphs using the general, nonempirical method of time-dependent (TD) optimally tuned screened range-separated hybrid (OT-SRSH) functional. For 5 ROY polymorphs with known experimental absorption spectra, we show that the TD-OT-SRSH approach predicts absorption spectra in quantitative agreement with experiment. For all polymorphs, we show that an accurate simulation of the colors is obtained, paving the way to a fully predictive, low-cost calculation of color polymorphism.