skip to main content

Title: Moving beyond the constraints of chemistry via crystal structure discovery with isotropic multiwell pair potentials

The rigid constraints of chemistry—dictated by quantum mechanics and the discrete nature of the atom—limit the set of observable atomic crystal structures. What structures are possible in the absence of these constraints? Here, we systematically crystallize one-component systems of particles interacting with isotropic multiwell pair potentials. We investigate two tunable families of pairwise interaction potentials. Our simulations self-assemble a multitude of crystal structures ranging from basic lattices to complex networks. Sixteen of the structures have natural analogs spanning all coordination numbers found in inorganic chemistry. Fifteen more are hitherto unknown and occupy the space between covalent and metallic coordination environments. The discovered crystal structures constitute targets for self-assembly and expand our understanding of what a crystal structure can look like.

more » « less
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Proceedings of the National Academy of Sciences
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Page Range / eLocation ID:
Article No. e2024034118
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Investigations into the reactivity, properties, and applications of osmium(IV) tetraaryl complexes have been hampered by their low yielding syntheses from volatile and toxic OsO4 (typically ≤34%). Here we show that known air-stable M(aryl)4 compounds (M = Os, Ru; aryl = 2-tolyl, 2,5-xylyl) can be prepared in ≤73% yields using new, less hazardous (Oct4N)2[MX6] precursors (M = Os, Ru; X = Cl, Br). This approach also facilitates the preparation of Os(mesityl)4 (Os3) for the first time, a complex comprising bulky 2,6-dimethyl substituted aryl ligands, albeit in low yield (5%). To better understand these yield extremes, we track, by synthesizing two additional new complexes with different 2-substituted σ-aryl ligands, a clear relationship between the yields of Os(aryl)4 and ligand steric bulk. Single-crystal X-ray structures of these compounds indicate that the observed yield trend reflects the ease of accommodating aryl substituents into an open pocket that lies directly opposite each M-aryl coordination site. We perform variable-temperature 1H NMR studies of Os3, utilize a "tetrahedricity" metric to assess geometric distortion in Ru(aryl)4 and Os(aryl)4 materials, and calculate cone angle and percentage buried volume metrics to further illustrate and help quantify -aryl ligand steric properties. Solution cyclic voltammograms of Os(aryl)4 show that the potentials of their reversible 1−/0 and 0/1+ redox features can be fine-tuned by varying aryl substituents, and that Os3 exhibits an additional 1+/2+ redox event not previously observed in this class of compounds. Taken together, this work helps to advance the potential application of these relatively underexplored organometallic complexes in established and emerging areas of molecular materials science, such as extended molecular frameworks and self-assembled monolayers, where analogous tetraphenylmethane and silane species (M = C, Si) have been frequently targeted. 
    more » « less
  2. Abstract

    Iron oxide is commonly found in natural or industrial glass compositions and can exist as ferrous (Fe2+) or ferric (Fe3+) species, with their ratios depending on glass composition, temperature, pressure and the redox reactions during the glass forming process. The iron redox ratio plays an important role on silicate glass structures and consequently various properties. This work aims to study the effect of iron oxide, and particularly the iron redox ratio, on the structures of borosilicate and boroaluminosilicate glasses using molecular dynamics simulations with newly developed iron potential parameters that are compatible with the borosilicate potentials. The results provide detailed cation coordination states of both iron species and the effect of redox ratio on boron coordination and other structural features. Particularly, competition for charge compensating modifier cations (such as Na+) among the fourfold‐coordinated cations such as B3+, Al3+, and Fe3+is investigated by calculating the cation–cation pair distribution functions and coordination preferential ratios. The results show that the trivalent ferric ions, with a shorter Fe–O bond distance and better defined first coordiation shell with mainly four‐fold coordination, act as a glass former whereas the divalent ferrous ions mainly play the role of glass modifier. The ferrous/ferric ratio (Fe2+/Fe3+) was found to affect the glass chemistry and hence glass properties by regulating the amount of four‐coordinated boron, the fraction of non‐briding oxygen and other features. The results are compared with available experimental data to gain insights of the complex structures and charge compensation schemes of the glass system.

    more » « less
  3. Abstract

    The coordination of traits from individual organs to whole plants is under strong selection because of environmental constraints on resource acquisition and use. However, the tight coordination of traits may provide underlying mechanisms of how locally adapted plant populations can become maladapted because of climate change.

    To better understand local adaptation in intraspecific trait coordination, we studied trait variability in the widely distributed foundation tree species,Populus fremontiiusing a common garden near the mid‐elevational point of this species distribution. We examined 28 traits encompassing four spectra: phenology, leaf economic spectrum (LES), whole‐tree architecture (Corner's Rule) and wood economic spectrum (WES).

    Based on adaptive syndrome theory, we hypothesized that trait expression would be coordinated among and within trait spectra, reflecting local adaptation to either exposure to freeze‐thaw conditions in genotypes sourced from high‐elevation populations or exposure to extreme thermal stress in genotypes sourced from low‐elevation populations.

    High‐elevation genotypes expressed traits within the phenology and WES that limit frost exposure and tissue damage. Specifically, genotypes sourced from high elevations had later mean budburst, earlier mean budset, higher wood densities, higher bark fractions and smaller xylem vessels than their low‐elevation counterparts. Conversely, genotypes sourced from low elevations expressed traits within the LES that prioritized hydraulic efficiency and canopy thermal regulation to cope with extreme heat exposure, including 40% smaller leaf areas, 67% higher stomatal densities and 34% higher mean theoretical maximum stomatal conductance. Low‐elevation genotypes also expressed a lower stomatal control over leaf water potentials that subsequently dropped to pressures that could induce hydraulic failure.

    Synthesis. Our results suggest thatPopulus fremontiiexpresses a high degree of coordination across multiple trait spectra to adapt to local climate constraints on photosynthetic gas exchange, growth and survival. These results, therefore, increase our mechanistic understanding of local adaptation and the potential effects of climate change that in turn, improves our capacity to identify genotypes that are best suited for future restoration efforts.

    more » « less
  4. null (Ed.)
    Abstract This paper introduces the use of topological data analysis (TDA) as an unsupervised machine learning tool to uncover classification criteria in complex inorganic crystal chemistries. Using the apatite chemistry as a template, we track through the use of persistent homology the topological connectivity of input crystal chemistry descriptors on defining similarity between different stoichiometries of apatites. It is shown that TDA automatically identifies a hierarchical classification scheme within apatites based on the commonality of the number of discrete coordination polyhedra that constitute the structural building units common among the compounds. This information is presented in the form of a visualization scheme of a barcode of homology classifications, where the persistence of similarity between compounds is tracked. Unlike traditional perspectives of structure maps, this new “Materials Barcode” schema serves as an automated exploratory machine learning tool that can uncover structural associations from crystal chemistry databases, as well as to achieve a more nuanced insight into what defines similarity among homologous compounds. 
    more » « less
  5. The digital alchemy framework is an extended ensemble simulation technique that incorporates particle attributes as thermodynamic variables, enabling the inverse design of colloidal particles for desired behavior. Here, we extend the digital alchemy framework for the inverse design of patchy spheres that self-assemble into target crystal structures. To constrain the potentials to non-trivial solutions, we conduct digital alchemy simulations with constant second virial coefficient. We optimize the size, range, and strength of patchy interactions in model triblock Janus spheres to self-assemble the 2D kagome and snub square lattices and the 3D pyrochlore lattice, and demonstrate self-assembly of all three target structures with the designed models. The particles designed for the kagome and snub square lattices assemble into high quality clusters of their target structures, while competition from similar polymorphs lower the yield of the pyrochlore assemblies. We find that the alchemically designed potentials do not always match physical intuition, illustrating the ability of the method to find nontrivial solutions to the optimization problem. We identify a window of second virial coefficients that result in self-assembly of the target structures, analogous to the crystallization slot in protein crystallization. 
    more » « less