Neutron Instruments for Research in Coordination Chemistry.

Neutron diffraction and spectroscopy offer unique insight into structures and properties of solids and molecular materials. All neutron instruments located at the various neutron sources are distinct, even if their designs are based on similar principles, and thus, they are usually less familiar to the community than commercial X-ray diffractometers and optical spectrometers. Major neutron instruments in the USA, which are open to scientists around the world, and examples of their use in coordination chemistry research are presented here, along with a list of similar instruments at main neutron facilities in other countries. The reader may easily and quickly find from this minireview an appropriate neutron instrument for research. The instruments include single-crystal and powder diffractometers to determine structures, inelastic neutron scattering (INS) spectrometers to probe magnetic and vibrational excitations, and quasielastic neutron scattering (QENS) spectrometers to study molecular dynamics such as methyl rotation on ligands. Key and unique features of the diffraction and neutron spectroscopy that are relevant to inorganic chemistry are reviewed.

Zr-Porphyrin Metal-Organic Framework as nanoreactor for boosting the formation of hydrogen clathrates.

We report the first experimental evidence for rapid formation of hydrogen clathrates under mild pressure and temperature conditions within the cavities of a zirconium-metalloporphyrin framework, specifically PCN-222. PCN-222 has been selected for its 1D mesoporous channels, high water-stability, and proper hydrophilic behavior. Firstly, we optimize a microwave (MW)-assisted method for the synthesis of nanosized PCN-222 particles with precise structure control (exceptional homogeneity in morphology and crystalline phase purity), taking advantage of MW in terms of rapid/homogeneous heating, time and energy savings, as well as potential scalability of the synthetic method. Second, we explore the relevance of the large mesoporous 1D open channels within the PCN-222 to promote the nucleation and growth of confined hydrogen clathrates. Experimental results show that PCN-222 drives the nucleation process at a lower pressure than the bulk system (1.35 kbar vs 2 kbar), with fast kinetics (minutes), using pure water, and with a nearly complete water-to-hydrate conversion. Unfortunately, PCN-222 cannot withstand these high pressures, which lead to a significant alteration of the mesoporous structure while the microporous network remains mainly unchanged.

Direct Conversion of Methane to Ethylene and Acetylene over an Iron-Based Metal-Organic Framework.

Conversion of methane (CH4) to ethylene (C2H4) and/or acetylene (C2H2) enables routes to a wide range of products directly from natural gas. However, high reaction temperatures and pressures are often required to activate and convert CH4 controllably, and separating C2+ products from unreacted CH4 can be challenging. Here, we report the direct conversion of CH4 to C2H4 and C2H2 driven by non-thermal plasma under ambient (25 °C and 1 atm) and flow conditions over a metal-organic framework material, MFM-300(Fe). The selectivity for the formation of C2H4 and C2H2 reaches 96% with a high time yield of 334 µmol gcat-1 h-1. At a conversion of 10%, the selectivity to C2+ hydrocarbons and time yield exceed 98% and 2056 µmol gcat-1 h-1, respectively, representing a new benchmark for conversion of CH4. In situ neutron powder diffraction, inelastic neutron scattering and solid-state nuclear magnetic resonance, electron paramagnetic resonance (EPR), and diffuse reflectance infrared Fourier transform spectroscopies, coupled with modeling studies, reveal the crucial role of Fe-O(H)-Fe sites in activating CH4 and stabilizing reaction intermediates via the formation of an Fe-O(CH3)-Fe adduct. In addition, a cascade fixed-bed system has been developed to achieve online separation of C2H4 and C2H2 from unreacted CH4 for direct use. Integrating the processes of CH4 activation, conversion, and product separation within one system opens a new avenue for natural gas utility, bridging the gap between fundamental studies and practical applications in this area.

Vibrations and Phase Stability in Mixed Valence Antimony Oxide.

α-Sb2O4 (cervantite) and ß-Sb2O4 (clinocervantite) are mixed valence compounds with equal proportions of SbIII and SbV as represented in the formula SbIIISbVO4. Their structure and properties can be difficult to calculate owing to the SbIII lone-pair electrons. Here, we present a study of the lattice dynamics and vibrational properties using a combination of inelastic neutron scattering, Mössbauer spectroscopy, nuclear inelastic scattering, and density functional theory (DFT) calculations. DFT calculations that account for lone-pair electrons match the experimental densities of phonon states. Mössbauer spectroscopy reveals the ß phase to be significantly harder than the α phase. Calculations with O vacancies reveal the possibility for nonstoichiometric proportions of SbIII and SbV in both phases. An open question is what drives the stability of the α phase over the ß phase, as the latter shows pronounced kinetic stability and lower symmetry despite being in the high-temperature phase. Since the vibrational entropy difference is small, it is unlikely to stabilize the α phase. Our results suggest that the α phase is more stable only because the material is not fully stoichiometric.

Neutron scattering study of polyamorphic THF·17(H2O) - toward a generalized picture of amorphous states and structures derived from clathrate hydrates.

From crystalline tetrahydrofuran clathrate hydrate, THF-CH (THF·17H2O, cubic structure II), three distinct polyamorphs can be derived. First, THF-CH undergoes pressure-induced amorphization when pressurized to 1.3 GPa in the temperature range 77-140 K to a form which, in analogy to pure ice, may be called high-density amorphous (HDA). Second, HDA can be converted to a densified form, VHDA, upon heat-cycling at 1.8 GPa to 180 K. Decompression of VHDA to atmospheric pressure below 130 K produces the third form, recovered amorphous (RA). Results from neutron scattering experiments and molecular dynamics simulations provide a generalized picture of the structure of amorphous THF hydrates with respect to crystalline THF-CH and liquid THF·17H2O solution (∼2.5 M). Although fully amorphous, HDA is heterogeneous with two length scales for water-water correlations (less dense local water structure) and guest-water correlations (denser THF hydration structure). The hydration structure of THF is influenced by guest-host hydrogen bonding. THF molecules maintain a quasiregular array, reminiscent of the crystalline state, and their hydration structure (out to 5 Å) constitutes ∼23H2O. The local water structure in HDA is reminiscent of pure HDA-ice featuring 5-coordinated H2O. In VHDA, the hydration structure of HDA is maintained but the local water structure is densified and resembles pure VHDA-ice with 6-coordinated H2O. The hydration structure of THF in RA constitutes ∼18 H2O molecules and the water structure corresponds to a strictly 4-coordinated network, as in the liquid. Both VHDA and RA can be considered as homogeneous.

Modulation of Uptake and Reactivity of Nitrogen Dioxide in Metal-Organic Framework Materials.

We report the modulation of reactivity of nitrogen dioxide (NO2 ) in a charged metal-organic framework (MOF) material, MFM-305-CH3 in which unbound N-centres are methylated and the cationic charge counter-balanced by Cl- ions in the pores. Uptake of NO2 into MFM-305-CH3 leads to reaction between NO2 and Cl- to give nitrosyl chloride (NOCl) and NO3 - anions. A high dynamic uptake of 6.58 mmol g-1 at 298 K is observed for MFM-305-CH3 as measured using a flow of 500 ppm NO2 in He. In contrast, the analogous neutral material, MFM-305, shows a much lower uptake of 2.38 mmol g-1 . The binding domains and reactivity of adsorbed NO2 molecules within MFM-305-CH3 and MFM-305 have been probed using in situ synchrotron X-ray diffraction, inelastic neutron scattering and by electron paramagnetic resonance, high-field solid-state nuclear magnetic resonance and UV/Vis spectroscopies. The design of charged porous sorbents provides a new platform to control the reactivity of corrosive air pollutants.

Synthesis of Nitro Compounds from Nitrogen Dioxide Captured in a Metal-Organic Framework.

Increasing levels of air pollution are driving the need for the development of new processes that take "waste-to-chemicals". Herein, we report the capture and conversion under ambient conditions of a major air pollutant, NO2, using a robust metal-organic framework (MOF) material, Zr-bptc (H4bptc = 3,3',5,5'-biphenyltetracarboxylic acid), comprising {Zr6(µ3-O)4(µ3-OH)4(COO)12} clusters linked by 4-connected bptc4- ligands in an ftw topology. At 298 K, Zr-bptc shows exceptional stability and adsorption of NO2 at both low (4.9 mmol g-1 at 10 mbar) and high pressures (13.8 mmol g-1 at 1.0 bar), as measured by isotherm experiments. Dynamic breakthrough experiments have confirmed the selective retention of NO2 by Zr-bptc at low concentrations under both dry and wet conditions. The immobilized NO2 can be readily transformed into valuable nitro compounds relevant to construction, agrochemical, and pharmaceutical industries. In situ crystallographic and spectroscopic studies reveal strong binding interactions of NO2 to the {Zr6(µ3-O)4(µ3-OH)4(COO)12} cluster node. This study paves a circular pathway to enable the integration of nitrogen-based air pollutants into the production of fine chemicals.

Vibrational Spectroscopy of Hexahalo Complexes.

Halogenated inorganic complexes Ax[MHaly] (A = alkali metal or alkaline earth, M = transition or main group metal, x = 1-3, and y = 2-9) are an archetypal class of compounds that provide entry points to large areas of inorganic and physical chemistry. All of the hexahalo complexes adopt an octahedral, Oh, symmetry (or nearly so). Consequently, one of the bending modes is forbidden in both the infrared and Raman spectra. In the solid state, many of the complexes crystallize in the cubic space group Fm3̅m, which preserves the octahedral symmetry. Even for those that are not cubic, the octahedral symmetry of the [MHal6]n- ion is largely retained and, to a good approximation, so are the selection rules. In the present work, we show that by using the additional information provided by neutron vibrational spectroscopy, in combination with conventional optical spectroscopies, we can generate complete and unambiguous assignments for all the modes. Comparison of the experimental and calculated transition energies for the systems where periodic-density functional theory was possible (i.e., those for which the crystal structure is known) shows that the agreement is almost quantitative. We also provide a linear relationship that enables the prediction of the forbidden mode.

Comprehensive Studies of Magnetic Transitions and Spin-Phonon Couplings in the Tetrahedral Cobalt Complex Co(AsPh3)2I2.

A combination of inelastic neutron scattering (INS), far-IR magneto-spectroscopy (FIRMS), and Raman magneto-spectroscopy (RaMS) has been used to comprehensively probe magnetic excitations in Co(AsPh3)2I2 (1), a reported single-molecule magnet (SMM). With applied field, the magnetic zero-field splitting (ZFS) peak (2D') shifts to higher energies in each spectroscopy. INS placed the ZFS peak at 54 cm-1, as revealed by both variable-temperature (VT) and variable-magnetic-field data, giving results that agree well with those from both far-IR and Raman studies. Both FIRMS and RaMS also reveal the presence of multiple spin-phonon couplings as avoided crossings with neighboring phonons. Here, phonons refer to both intramolecular and lattice vibrations. The results constitute a rare case in which the spin-phonon couplings are observed with both Raman-active (g modes) and far-IR-active phonons (u modes; space group P21/c, no. 14, Z = 4 for 1). These couplings are fit using a simple avoided crossing model with coupling constants of ca. 1-2 cm-1. The combined spectroscopies accurately determine the magnetic excited level and the interaction of the magnetic excitation with phonon modes. Density functional theory (DFT) phonon calculations compare well with INS, allowing for the assignment of the modes and their symmetries. Electronic calculations elucidate the nature of ZFS in the complex. Features of different techniques to determine ZFS and other spin-Hamiltonian parameters in transition-metal complexes are summarized.

Joint neutron/molecular dynamics vibrational spectroscopy reveals softening of HIV-1 protease upon binding of a tight inhibitor.

Biomacromolecules are inherently dynamic, and their dynamics are interwoven into function. The fast collective vibrational dynamics in proteins occurs in the low picosecond timescale corresponding to frequencies of ∼5-50 cm-1. This sub-to-low THz frequency regime covers the low-amplitude collective breathing motions of a whole protein and vibrations of the constituent secondary structure elements, such as α-helices, ß-sheets and loops. We have used inelastic neutron scattering experiments in combination with molecular dynamics simulations to demonstrate the vibrational dynamics softening of HIV-1 protease, a target of HIV/AIDS antivirals, upon binding of a tight clinical inhibitor darunavir. Changes in the vibrational density of states of matching structural elements in the two monomers of the homodimeric protein are not identical, indicating asymmetric effects of darunavir on the vibrational dynamics. Three of the 11 major secondary structure elements contribute over 40% to the overall changes in the vibrational density of states upon darunavir binding. Molecular dynamics simulations informed by experiments allowed us to estimate that the altered vibrational dynamics of the protease would contribute -3.6 kcal mol-1 at 300 K, or 25%, to the free energy of darunavir binding. As HIV-1 protease drug resistance remains a concern, our results open a new avenue to help establish a direct quantitative link between protein vibrational dynamics and drug resistance.

A multi-faceted structural, thermodynamic, and spectroscopic approach for investigating ethanol dehydration over transition phase aluminas.

Understanding the role that the surface of a material plays in the mediation of a chemical reaction at the atomic level is paramount to the optimization and improvement of catalytic materials. While this area of research has matured over several decades, few techniques are sensitive enough to directly examine and differentiate the behavior of molecular adsorbates during the course of the chemical reaction with a substrate. In this study, a combined approach which involves structural characterization techniques, volumetric adsorption, temperature programmed desorption, and inelastic neutron scattering (INS) was used to investigate the mechanism of ethanol dehydration on the surface of transition phase aluminas. The alumina samples employed were extensively characterized using X-ray diffraction, solid-state 27Al nuclear magnetic resonance spectroscopy, and thermogravimetric analysis with differential scanning calorimetry. A high-precision volumetric adsorption apparatus was used to characterize the surface area and to controllably dose ethanol onto the surface of the aluminas. A modified temperature programmed desorption (TPD) method which samples the molecular composition of the vapor at discrete temperatures in a closed cell is described. INS results were used to confirm adsorption of ethanol on γ- and θ-alumina and show the reaction of ethanol and subsequent formation of ethylene as a function of temperature. The TPD and INS results affirm that the dehydration reaction and subsequent formation of ethylene on both γ- and θ-aluminas occur rapidly at 300 °C, though ethanol is still observed on θ-alumina indicating fewer active sites. These results demonstrate the value of a multi-faceted characterization approach, featuring INS, towards providing a detailed understanding of the ethanol dehydration mechanism on θ-alumina and further provide the basis for extending this approach to other systems in heterogeneous catalysis and areas where molecule-substrate interactions are poorly understood.

Inelastic Neutron Scattering from Thin Film Biaxially Oriented Polyethylene Terephthalate.

Recent interest in emerging processes for polymer manufacturing and bio-based chemistries for direct chemical recycling/upcycling has motivated new research focused on a deeper understanding of atomic-scale polymer properties and how they influence macroscopic phenomena. Uncovering the fundamental properties of polymers that give rise to macroscopic behavior could enable new pathways for improved recyclability or utilization of alternative "greener" polymer analogues. In this study, the neutron vibrational spectrum was measured for a film of biaxially oriented polyethylene terephthalate (BoPET) using inelastic neutron scattering (INS), to investigate the relationship between the structure and dynamics of a widely used polymer. Compared to conventional spectroscopic techniques, the use of INS is advantageous for polymeric materials due to the absence of selection rules (i.e., all transitions are allowed), broad-band energy range, and considerable sensitivity to hydrogen modes. In order to distinguish the vibrational modes caused by trans and gauche rotational isomerism, the normal modes of vibration were calculated from a density functional theory-optimized structure of crystalline PET (cPET), representative of the all-trans state, and compared with INS from "highly crystalline" PET powder. Although in- and out-of-plane wagging of hydrogens on the ring structure exhibit significant contribution to both BoPET and cPET spectra, the wagging, rocking, and twisting modes of hydrogen on the ethylene glycol group are, in most cases, conformation-specific. These results were further rationalized by investigating the role of hyperconjugation in stabilizing both conformations using the natural bond order method. Through comparison of experimental and calculated INS results, this work provides the fundamental basis for discovering the role of structure and dynamics in shaping the macroscopic properties of PET and polymer analogues.

Structural and Dynamic Analysis of Sulphur Dioxide Adsorption in a Series of Zirconium-Based Metal-Organic Frameworks.

We report reversible high capacity adsorption of SO2 in robust Zr-based metal-organic framework (MOF) materials. Zr-bptc (H4 bptc=biphenyl-3,3',5,5'-tetracarboxylic acid) shows a high SO2 uptake of 6.2 mmol g-1 at 0.1 bar and 298 K, reflecting excellent capture capability and removal of SO2 at low concentration (2500 ppm). Dynamic breakthrough experiments confirm that the introduction of amine, atomically-dispersed CuII or heteroatomic sulphur sites into the pores enhance the capture of SO2 at low concentrations. The captured SO2 can be converted quantitatively to a pharmaceutical intermediate, aryl N-aminosulfonamide, thus converting waste to chemical values. In situ X-ray diffraction, infrared micro-spectroscopy and inelastic neutron scattering enable the visualisation of the binding domains of adsorbed SO2 molecules and host-guest binding dynamics in these materials at the atomic level. Refinement of the pore environment plays a critical role in designing efficient sorbent materials.

Exceptional Packing Density of Ammonia in a Dual-Functionalized Metal-Organic Framework.

We report the reversible adsorption of ammonia (NH3) up to 9.9 mmol g-1 in a robust Al-based metal-organic framework, MFM-303(Al), which is functionalized with free carboxylic acid and hydroxyl groups. The unique pore environment decorated with these acidic sites results in an exceptional packing density of NH3 at 293 K (0.801 g cm-3) comparable to that of solid NH3 at 193 K (0.817 g cm-3). In situ synchrotron X-ray diffraction and inelastic neutron scattering reveal the critical role of free -COOH and -OH groups in immobilizing NH3 molecules. Breakthrough experiments confirm the excellent performance of MFM-303(Al) for the capture of NH3 at low concentrations under both dry and wet conditions.

High Ammonia Adsorption in MFM-300 Materials: Dynamics and Charge Transfer in Host-Guest Binding.

Ammonia (NH3) is a promising energy resource owing to its high hydrogen density. However, its widespread application is restricted by the lack of efficient and corrosion-resistant storage materials. Here, we report high NH3 adsorption in a series of robust metal-organic framework (MOF) materials, MFM-300(M) (M = Fe, V, Cr, In). MFM-300(M) (M = Fe, VIII, Cr) show fully reversible capacity for >20 cycles, reaching capacities of 16.1, 15.6, and 14.0 mmol g-1, respectively, at 273 K and 1 bar. Under the same conditions, MFM-300(VIV) exhibits the highest uptake among this series of MOFs of 17.3 mmol g-1. In situ neutron powder diffraction, single-crystal X-ray diffraction, and electron paramagnetic resonance spectroscopy confirm that the redox-active V center enables host-guest charge transfer, with VIV being reduced to VIII and NH3 being oxidized to hydrazine (N2H4). A combination of in situ inelastic neutron scattering and DFT modeling has revealed the binding dynamics of adsorbed NH3 within these MOFs to afford a comprehensive insight into the application of MOF materials to the adsorption and conversion of NH3.

Quantitative production of butenes from biomass-derived γ-valerolactone catalysed by hetero-atomic MFI zeolite.

The efficient production of light olefins from renewable biomass is a vital and challenging target to achieve future sustainable chemical processes. Here we report a hetero-atomic MFI-type zeolite (NbAlS-1), over which aqueous solutions of γ-valerolactone (GVL), obtained from biomass-derived carbohydrates, can be quantitatively converted into butenes with a yield of >99% at ambient pressure under continuous flow conditions. NbAlS-1 incorporates simultaneously niobium(V) and aluminium(III) centres into the framework and thus has a desirable distribution of Lewis and Brønsted acid sites with optimal strength. Synchrotron X-ray diffraction and absorption spectroscopy show that there is cooperativity between Nb(V) and the Brønsted acid sites on the confined adsorption of GVL, whereas the catalytic mechanism for the conversion of the confined GVL into butenes is revealed by in situ inelastic neutron scattering, coupled with modelling. This study offers a prospect for the sustainable production of butene as a platform chemical for the manufacture of renewable materials.

Applying Unconventional Spectroscopies to the Single-Molecule Magnets, Co(PPh3 )2 X2 (X=Cl, Br, I): Unveiling Magnetic Transitions and Spin-Phonon Coupling.

Large separation of magnetic levels and slow relaxation in metal complexes are desirable properties of single-molecule magnets (SMMs). Spin-phonon coupling (interactions of magnetic levels with phonons) is ubiquitous, leading to magnetic relaxation and loss of memory in SMMs and quantum coherence in qubits. Direct observation of magnetic transitions and spin-phonon coupling in molecules is challenging. We have found that far-IR magnetic spectra (FIRMS) of Co(PPh3 )2 X2 (Co-X; X=Cl, Br, I) reveal rarely observed spin-phonon coupling as avoided crossings between magnetic and u-symmetry phonon transitions. Inelastic neutron scattering (INS) gives phonon spectra. Calculations using VASP and phonopy programs gave phonon symmetries and movies. Magnetic transitions among zero-field split (ZFS) levels of the S=3/2 electronic ground state were probed by INS, high-frequency and -field EPR (HFEPR), FIRMS, and frequency-domain FT terahertz EPR (FD-FT THz-EPR), giving magnetic excitation spectra and determining ZFS parameters (D, E) and g values. Ligand-field theory (LFT) was used to analyze earlier electronic absorption spectra and give calculated ZFS parameters matching those from the experiments. DFT calculations also gave spin densities in Co-X, showing that the larger Co(II) spin density in a molecule, the larger its ZFS magnitude. The current work reveals dynamics of magnetic and phonon excitations in SMMs. Studies of such couplings in the future would help to understand how spin-phonon coupling may lead to magnetic relaxation and develop guidance to control such coupling.

La2Zr2O7 Nanoparticle-Mediated Synthesis of Porous Al-Doped Li7La3Zr2O12 Garnet.

In this work, we modified the reaction pathway to quickly (minutes) incorporate lithium and stabilize the ionic conducting garnet phase by decoupling the formation of a La-Zr-O network from the addition of lithium. To do this, we synthesized La2Zr2O7 (LZO) nanoparticles to which LiNO3 was added. This method is a departure from typical solid-state synthesis methods that require high-energy milling to promote mixing and intimate particle-particle contact and from sol-gel syntheses as a unique porous microstructure is obtained. We show that the reaction time is limited by the rate of nitrate decomposition and that this method produces a porous high-Li-ion-conducting cubic phase, within an hour, that may be used as a starting structure for a composite electrolyte.

Theory-Guided Inelastic Neutron Scattering of Crystalline Alkaline Aluminate Salts Bearing Principal Motifs of Solution-State Species.

Aluminate salts precipitated from caustic alkaline solutions exhibit a correlation between the anionic speciation and the identity of the alkali cation in the precipitate, with the aluminate ions occurring either in monomeric (Al(OH)4-) or dimeric (Al2O(OH)62-) forms. The origin of this correlation is poorly understood as are the roles that oligomeric aluminate species play in determining the solution structure, prenucleation clusters, and precipitation pathways. Characterization of aluminate solution speciation with vibrational spectroscopy results in spectra that are difficult to interpret because the ions access a diverse and dynamic configurational space. To investigate the Al(OH)4- and Al2O(OH)62- anions within a well-defined crystal lattice, inelastic neutron scattering (INS) and Raman spectroscopic data were collected and simulated by density functional theory for K2[Al2O(OH)6], Rb2[Al2O(OH)6], and Cs[Al(OH) 4]·2H2O. These structures capture archetypal solution aluminate species: the first two salts contain dimeric Al2O(OH)62- anions, while the third contains the monomeric Al(OH)4- anion. Comparisons were made to the INS and Raman spectra of sodium aluminate solutions frozen in a glassy state. In contrast to solution systems, the crystal lattice of the salts results in well-defined vibrations and associated resolved bands in the INS spectra. The use of a theory-guided analysis of the INS of this solid alkaline aluminate series revealed that differences were related to the nature of the hydrogen-bonding network and showed that INS is a sensitive probe of the degree of completeness and strength of the bond network in hydrogen-bonded materials. Results suggest that the ionic size may explain cation-specific differences in crystallization pathways in alkaline aluminate salts.

Davis Computational Spectroscopy Workflow-From Structure to Spectra.

We describe an automated workflow that connects a series of atomic simulation tools to investigate the relationship between atomic structure, lattice dynamics, materials properties, and inelastic neutron scattering (INS) spectra. Starting from the atomic simulation environment (ASE) as an interface, we demonstrate the use of a selection of calculators, including density functional theory (DFT) and density functional tight binding (DFTB), to optimize the structures and calculate interatomic force constants. We present the use of our workflow to compute the phonon frequencies and eigenvectors, which are required to accurately simulate the INS spectra in crystalline solids like diamond and graphite as well as molecular solids like rubrene. We have also implemented a machine-learning force field based on Chebyshev polynomials called the Chebyshev interaction model for efficient simulation (ChIMES) to improve the accuracy of the DFTB simulations. We then explore the transferability of our DFTB/ChIMES models by comparing simulations derived from different training sets. We show that DFTB/ChIMES demonstrates ∼100× reduction in computational expense while retaining most of the accuracy of DFT as well as yielding high accuracy for different materials outside of our training sets. The DFTB/ChIMES method within the workflow expands the possibilities to use simulations to accurately predict materials properties of increasingly complex structures that would be unfeasible with ab initio methods.