Search

With the goal of accelerating the design and discovery of metal–organic frameworks (MOFs) for electronic, optoelectronic, and energy storage applications, we present a dataset of predicted electronic structure properties for thousands of MOFs carried out using multiple density functional approximations. Compared to more accurate hybrid functionals, we find that the widely used PBE generalized gradient approximation (GGA) functional severely underpredicts MOF band gaps in a largely systematic manner for semi-conductors and insulators without magnetic character. However, an even larger and less predictable disparity in the band gap prediction is present for MOFs with open-shell 3d transition metal cations. With regards to partial atomic charges, we find that different density functional approximations predict similar charges overall, although hybrid functionals tend to shift electron density away from the metal centers and onto the ligand environments compared to the GGA point of reference. Much more significant differences in partial atomic charges are observed when comparing different charge partitioning schemes. We conclude by using the dataset of computed MOF properties to train machine-learning models that can rapidly predict MOF band gaps for all four density functional approximations considered in this work, paving the way for future high-throughput screening studies. To encourage exploration and reuse of the theoretical calculations presented in this work, the curated data is made publicly available via an interactive and user-friendly web application on the Materials Project. [ABSTRACT FROM AUTHOR]

Multiconfiguration pair-density functional (MC-PDFT) theory provides an economical way to calculate the ground-state and excited-state energetics of strongly correlated systems. The energy is calculated from the kinetic energy, density, and on-top pair-density of a multiconfiguration wave function as the sum of kinetic energy, classical Coulomb energy, and on-top density functional energy. We have usually found good results with the translated Perdew–Burke–Ernzerhof (tPBE) on-top density functional, and in this article, we examine whether the results can be systematically improved by introducing scaling constants into the exchange and correlation terms. We find that only a small improvement is possible for electronic excitation energies and that no improvement is possible for bond energies. [ABSTRACT FROM AUTHOR]

This work examines the question of how a Hubbard U correction to a local exchange–correlation functional compares with adding Hartree–Fock exchange to a local functional for both solid-state and molecular properties. We compute a solid-state property, namely the band gap, and thermochemical molecular properties, in particular, main-group bond energies, transition metal–ligand bond energies, and barrier heights, to elucidate whether the DFT+U method mimics hybrid DFT. We find that a calculation with a Hubbard U correction may or may not mimic a hybrid functional—depending on the atom, the subshell, and the property to which it is applied. For band gaps, we find that adding a Hubbard U correction to the valence d orbitals of transition metals increases the band gap, which thereby gets closer to the experimental value, while adding a Hubbard U correction to valence s or p orbitals of main-group elements need not always increase the band gap. For molecular thermochemistry, we find that adding a Hubbard U correction to a local density functional need not have the same effect as adding Hartree–Fock exchange to a local density functional. For example when compared to a DFT calculation with a local exchange-correlation functional, hybrid DFT increases the barrier height in all cases, but DFT+U does not always increase the barrier height. For the band gaps of transition metal monoxides, the Hubbard-corrected results lowered the mean errors significantly and were comparable to what could be achieved with a much more expensive hybrid functional, but for reaction barrier heights and bond energies of molecules, the Hubbard correction was found to lower the mean error by only approximately a kcal/mol. As part of the analysis, we also compare VASP and Gaussian 09 calculations for the same density functional. [ABSTRACT FROM AUTHOR]

Enzymatic haem and non-haem high-valent iron-oxo species are known to activate strong C-H bonds, yet duplicating this reactivity in a synthetic system remains a formidable challenge. Although instability of the terminal iron-oxo moiety is perhaps the foremost obstacle, steric and electronic factors also limit the activity of previously reported mononuclear iron(IV)-oxo compounds. In particular, although nature's non-haem iron(IV)-oxo compounds possess high-spin S = 2 ground states, this electronic configuration has proved difficult to achieve in a molecular species. These challenges may be mitigated within metal-organic frameworks that feature site-isolated iron centres in a constrained, weak-field ligand environment. Here, we show that the metal-organic framework Fe2(dobdc) (dobdc4− = 2,5-dioxido-1,4-benzenedicarboxylate) and its magnesium-diluted analogue, Fe0.1Mg1.9(dobdc), are able to activate the C-H bonds of ethane and convert it into ethanol and acetaldehyde using nitrous oxide as the terminal oxidant. Electronic structure calculations indicate that the active oxidant is likely to be a high-spin S = 2 iron(IV)-oxo species. [ABSTRACT FROM AUTHOR]

We applied the solvation models SM8, SM8AD, and SMD in combination with the Minnesota M06-2X density functional to predict vacuum-water transfer free energies (Task 1) and tautomeric ratios in aqueous solution (Task 2) for the SAMPL2 test set. The bulk-electrostatic contribution to the free energy of solvation is treated as follows: SM8 employs the generalized Born model with the Coulomb field approximation, SM8AD employs the generalized Born approximation with asymmetric descreening, and SMD solves the nonhomogeneous Poisson equation. The non-bulk-electrostatic contribution arising from short-range interactions between the solute and solvent molecules in the first solvation shell is treated as a sum of terms that are products of geometry-dependent atomic surface tensions and solvent-accessible surface areas of the individual atoms of the solute. On average, three models tested in the present work perform similarly. In particular, we achieved mean unsigned errors of 1.3 (SM8), 2.0 (SM8AD), and 2.6 kcal/mol (SMD) for the aqueous free energies of 30 out of 31 compounds with known reference data involved in Task 1 and mean unsigned errors of 2.7 (SM8), 1.8 (SM8AD), and 2.4 kcal/mol (SMD) in the free energy differences (tautomeric ratios) for 21 tautomeric pairs in aqueous solution involved in Task 2. [ABSTRACT FROM AUTHOR]

We present two new hybrid meta exchange- correlation functionals, called M06 and M06-2X. The M06 functional is parametrized including both transition metals and nonmetals, whereas the M06-2X functional is a high-nonlocality functional with double the amount of nonlocal exchange (2X), and it is parametrized only for nonmetals.The functionals, along with the previously published M06-L local functional and the M06-HF full-Hartree–Fock functionals, constitute the M06 suite of complementary functionals. We assess these four functionals by comparing their performance to that of 12 other functionals and Hartree–Fock theory for 403 energetic data in 29 diverse databases, including ten databases for thermochemistry, four databases for kinetics, eight databases for noncovalent interactions, three databases for transition metal bonding, one database for metal atom excitation energies, and three databases for molecular excitation energies. We also illustrate the performance of these 17 methods for three databases containing 40 bond lengths and for databases containing 38 vibrational frequencies and 15 vibrational zero point energies. We recommend the M06-2X functional for applications involving main-group thermochemistry, kinetics, noncovalent interactions, and electronic excitation energies to valence and Rydberg states. We recommend the M06 functional for application in organometallic and inorganometallic chemistry and for noncovalent interactions. [ABSTRACT FROM AUTHOR]

Hybrid density functional theory, which is a combined Hartree-Fock and density functional method, provides a simple but effective way to incorporate nonlocal exchange effects and static and dynamical correlation energy into an orbital-based theory with affordable computational cost for many important problems of gas-phase chemistry. The inclusion of a reaction field representing an implicit solvent in a self-consistent hybrid density functional calculation provides an effective and efficient way to extend this approach to problems of liquid-phase chemistry. In previous work, we have parameterized several models based on this approach, and in the present article, we present several new parameterizations based on implicit solvation models SM5.43 and SM5.43R. In particular, we extend the applicability of these solvation models to several combinations of the MPWXhybrid-density functional with various one-electron basis sets, where MPWXdenotes a combination of Barone and Adamo’s modified version of Perdew and Wang’s exchange functional, Perdew and Wang’s correlation functional, and a percentageXof exact Hartree-Fock exchange. SM5.43R parameter optimizations are presented for the MPWX/MIDI!, MPWX/MIDI!6D, and MPWX/6-31+G(d,p) combinations withX=0 (i.e., pure density functional theory), 25, 42.8, and 60.6, and for MPWX/6-31G(d) and MPWX/6-31+G(d), withX=0, 42.8, and 60.6; this constitutes a total of 18 new parameter sets. [Note that parameter optimizations using MPW25/6-31G(d) and MPW25/6-31+G(d) were carried out in a previous SM5.43R parameterization.] For each of the five basis sets, we found no significant loss in the accuracy of the model when parameters averaged over the four values ofXare used instead of the parameters optimized for a specific value ofX. Therefore for each of the five basis sets used here, the SM5.43R and SM5.43 models are defined to have a single parameter set that can be used for any value ofXbetween 0 and 60.6. The new models yield accurate free energies of solvation for a broad range of solutes in both water and organic solvents. On the average, the mean-unsigned errors, as compared with those from experiment, of the free energies of solvation of neutral solutes range from 0.50 to 0.55 kcal/mol and those for ions range from 4.5 to 4.9 kcal/mol. Since the SM5.43R model computes the free energy of solvation as a sum of bulk-electrostatic and non-bulk-electrostatic contributions, it may be used for detailed analysis of the physical effects underlying a calculation of the free energy of solvation. Several calculations illustrating the partitioning of these contributions for a variety of solutes inn-hexadecane, 1-octanol, and water are presented. [ABSTRACT FROM AUTHOR]

Equilibrium free-energy cycles relating oxidation and reduction potentials in solution to ionization potentials and electron affinities in the gas phase are constructed and the utilities of various levels of theory for computing particular free-energy changes within these cycles are discussed within the context of several examples. Emphasis is placed on the use of quantum-mechanical continuum solvation models to compute free energies of solvation. Key systems discussed include quinones, substituted anilines, substituted phenols, and reductive dechlorination reactions. [ABSTRACT FROM AUTHOR]

Beginning with the MIDI! basis set (also called MIDIX), we introduce the MIDIX+, MIDIY, and MIDIY+ basis sets. Using correlated ab initio and hybrid density functional theory, we compare their performance to that of several existing basis sets for electronic structure calculations. The new basis sets are tested with databases of 358 energies of reactions, 44 barrier heights, 31 electron affinities, 18 geometries, and 29 dipole moments. The MIDI!, MIDIX+, MIDIY, and MIDIY+ basis sets are shown to be cost-efficient methods for calculating relative energies, geometries, and dipole moments. The MIDIX+ basis is shown to be particularly efficient for calculating electron affinities of large molecules. [ABSTRACT FROM AUTHOR]

Molecular interactions with both oxides and metals are essential for heterogenous catalysis, leading to remarkable synergistic impacts on activity and selectivity. Here, we show that the direct link between the two phases (and not merely being together) is required to selectively hydrogenate CO2 to methanol on catalysts containing Cu and ZrO2. Materials consisting of isolated Cu particles or atomically dispersed Cu–O–Zr sites only catalyze the reverse water-gas shift reaction. In contrast, a metal organic framework structure (UiO-66) with Cu nanoparticles occupying missing-linker defects maximizes the fraction of metallic Cu interfaced to ZrO2 nodes leading to a material with high adsorption capacity for CO2 and high activity and selectivity for low-temperature methanol synthesis. Molecular interactions with both oxides and metals are essential for heterogenous catalysis, leading to remarkable impacts on activity. Here the authors show that a direct link between Cu and ZrO2 in a metal organic framework is required to hydrogenate CO2 to methanol. [ABSTRACT FROM AUTHOR]

Summarizes matters involved at the Symposium on High-Performance Computer Applications in the Behavioral Sciences. Overview of the activities in scientific computation; Inclusion of international symposia and a graduate degree program in the activities; Emergence of scientific computation as an important field in academia and industry.

A highly selective Ir catalyst supported on the metal–organic framework (MOF) UiO-67 for the catalytic borylation of methane has recently been synthesized. The high chemoselectivity of the catalyst toward monoborylated methane (CH3Bpin, Bpin = pinacolborane) instead of diborylated methane (CH2Bpin2) was speculated to be caused by the steric confinement of MOF UiO-67. In this study, we applied quantum mechanical methods to determine: (1) the steric effect of the UiO-67 framework in promoting the chemoselectivity of the Ir catalyst toward CH3Bpin and (2) the borylation mechanisms over the Ir catalyst supported on UiO-67. Our results show that UiO-67 framework sterically obstructs the diffusion of the larger CH2Bpin2 molecule within the MOF while allowing the smaller CH3Bpin molecule to pass through with little energy penalty. The diffusion of CH2Bpin2 from the tetrahedral pore to the tetragonal pyramidal pore within modified UiO-67 with coordinated Ir(Bpin)3 complex has an estimated barrier of 24.7 kcal/mol and is 14.2 kcal/mol higher than the diffusion of CH3Bpin. The electronic and steric effects of the support at the Ir catalytic center are much smaller than this confinement effect on diffusion, and the catalytic center behaves similarly to the homogeneous Ir catalyst. We determined an overall free energy of activation of 34.6 kcal/mol for the CH4 borylation reaction using the Ir(III) catalyst. We also determined that the turnover-determining step for the catalytic methane borylation is the isomerization of seven-coordinated Ir(V) complex instead of the commonly assumed C–H bond activation by oxidative addition. [ABSTRACT FROM AUTHOR]

Elucidating atmospheric oxidation mechanisms is necessary for estimating the lifetimes of atmospheric species and understanding secondary organic aerosol formation and atmospheric oxidation capacity. We report an unexpectedly fast mechanistic pathway for the unimolecular reactions of large stabilized Criegee intermediates, which involves the formation of bicyclic structures from large Criegee intermediates containing an aldehyde group. The barrier heights of the mechanistic pathways are unexpectedly low – about 2–3 kcal/mol – and are at least 10 kcal/mol lower than those of hydrogen shift processes in large syn Criegee intermediates; and the calculated rate constants show that the mechanistic pathways are 105-109 times faster than those of the corresponding hydrogen shift processes. The present findings indicate that analogous low-energy pathways can now also be expected in other large Criegee intermediates and that oxidative capacity of some Criegee intermediates is smaller than would be predicted by existing models. Criegee intermediates have received much attention in atmospheric chemistry because of their importance in ozonolysis mechanisms. Here, using quantum mechanical kinetics, the authors reveal an unexpectedly fast mechanistic pathway for unimolecular reactions of large stabilized Criegee intermediates. [ABSTRACT FROM AUTHOR]

The continuously increasing number and size of lithium-based batteries developed for large-scale applications raise serious environmental concerns. Herein, we address the issues related to electrolyte toxicity and safety by proposing a "water-in-ionomer" type of electrolyte which replaces organic solvents by water and expensive and toxic fluorinated lithium salts by a non-fluorinated, inexpensive and non-toxic superabsorbing ionomer, lithium polyacrylate. Interestingly, the electrochemical stability window of this electrolyte is extended greatly, even for high water contents. Particularly, the gel with 50 wt% ionomer exhibits an electrochemical stability window of 2.6 V vs. platinum and a conductivity of 6.5 mS cm−1 at 20 °C. Structural investigations suggest that the electrolytes locally self-organize and most likely switch local structures with the change of water content, leading to a 50% gel with good conductivity and elastic properties. A LiTi2(PO4)3/LiMn2O4 lithium-ion cell incorporating this electrolyte provided an average discharge voltage > 1.5 V and a specific energy of 77 Wh kg−1, while for an alternative cell chemistry, i.e., TiO2/LiMn2O4, a further enhanced average output voltage of 2.1 V and an initial specific energy of 124.2 Wh kg−1 are achieved. The flammability and toxicity of the currently used electrolytes are the concerns that must be addressed. Here the authors show a non-fluorinated and non-toxic ionomeric aqueous gel electrolyte called water-in-ionomer that allows an enlargement of electrochemical stability window and design of environmentally friendly battery cell chemistries. [ABSTRACT FROM AUTHOR]

The article offers the author's insights on a study conducted by Claire Adjiman and colleagues which developed a method that allows identification of solvent to improve the kinetics of a reaction. He says that the focus on free energies is one key advantage in modern reaction dynamics, and describes the transition state theory. He states the definition of an inverse catalyst design as described by the Purdue University Center for Catalyst Design.

The article presents information on the journal "Theoretical Chemistry Accounts: Theory, Computation, & Modeling" (TCA). The journal is opening a new editorial office in China, where chemistry professor Yundong Wu will be associate editor. It is reported that TCA is the oldest journal in the world dedicated to the topic of theoretical chemistry, and it has a long and distinguished history of publishing some of the most important papers in the field.

A correction to the article "The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06 Functionals and 12 Other Functionals," that was published in a previous issue is presented.

The article reflects on the developments of the 2006 issue of the periodical "Theoretical Chemistry Accounts." The recent issue presented 33 articles which cover areas related to theoretical chemistry by modern theoretical chemists. Senior advisers suggested contributors including Biman Bagchi, Vincenzo Barone, Charles Brooks and other colleagues.

Kohn-Sham density functional theory is widely used in chemistry, but no functional can accurately predict the whole range of chemical properties, although recent progress by some doubly hybrid functionals comes close. Here, we optimized a singly hybrid functional called CF22D with higher across-the-board accuracy for chemistry than most of the existing non-doubly hybrid functionals by using a flexible functional form that combines a global hybrid meta-nonseparable gradient approximation that depends on density and occupied orbitals with a damped dispersion term that depends on geometry. We optimized this energy functional by using a large database and performance-triggered iterative supervised training. We combined several databases to create a very large, combined database whose use demonstrated the good performance of CF22D on barrier heights, isomerization energies, thermochemistry, noncovalent interactions, radical and nonradical chemistry, small and large systems, simple and complex systems and transition-metal chemistry., (© 2022. The Author(s).)