1
|
Rastetter U, Jacobi von Wangelin A, Herrmann C. Redox-active ligands as a challenge for electronic structure methods. J Comput Chem 2023; 44:468-479. [PMID: 36326153 DOI: 10.1002/jcc.27013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2022] [Revised: 07/27/2022] [Accepted: 08/19/2022] [Indexed: 11/06/2022]
Abstract
To improve the catalytic activity of 3d transition metal catalysts, redox-active ligands are a promising tool. These ligands influence the oxidation state of the metal center as well as the ground spin-state and make the experimental determination of both properties challenging. Therefore, first-principles calculations, in particular employing density functional theory with a proper choice of exchange-correlation (xc) functional, are crucial. Common xc functionals were tested on a simple class of metal complexes: homoleptic, octahedral tris(diimine) iron(II) complexes. The spin-state energy splittings for most of these complexes showed the expected linear dependence on the amount of exact exchange included in the xc functionals. Even though varying redox-activity affects the electronic structure of the complexes considerably, the sensitivity of the spin-state energetics to the exact exchange admixture is surprisingly small. For iron(II) complexes with highly redox-active ligands and for a broad range of ligands in the reduced tris(diimine) iron(I) complexes, self-consistent field convergence to local minima was observed, which differ from the global minimum in the redox state of the ligand. This may also result in convergence to a molecular structure that corresponds to an energetically higher-lying local minimum. One criterion to detect such behavior is a change in the sign of the slope for the dependence of the spin-state energy splittings on the amount of exact exchange. We discuss possible protocols for dealing with such artifacts in cases in which a large number of calculations makes checking by hand unfeasible.
Collapse
Affiliation(s)
- Ursula Rastetter
- Department of Chemistry, University of Hamburg, Hamburg, Germany
| | | | - Carmen Herrmann
- Department of Chemistry, University of Hamburg, Hamburg, Germany
| |
Collapse
|
2
|
Reimann M, Kaupp M. Spin-State Splittings in 3d Transition-Metal Complexes Revisited: Benchmarking Approximate Methods for Adiabatic Spin-State Energy Differences in Fe(II) Complexes. J Chem Theory Comput 2022; 18:7442-7456. [PMID: 36417564 DOI: 10.1021/acs.jctc.2c00924] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
The CASPT2+δMRCI composite approach reported in a companion paper has been extended and used to provide high-quality reference data for a series of adiabatic spin gaps (defined as ΔE = Equintet - Esinglet) of [FeIIL6]2+ complexes (L = CNH, CO, NCH, NH3, H2O), either at nonrelativistic level or including scalar relativistic effects. These highly accurate data have been used to evaluate the performance of various more approximate methods. Coupled-cluster theory with singles, doubles, and perturbative triples, CCSD(T), is found to agree well with the new reference data for Werner-type complexes but exhibits larger underestimates by up to 70 kJ/mol for the π-acceptor ligands, due to appreciable static correlation in the low-spin states of these systems. Widely used domain-based local CCSD(T) calculations, DLPNO-CCSD(T), are shown to depend very sensitively on the cutoff values used to construct the localized domains, and standard values are not sufficient. A large number of density functional approximations have been evaluated against the new reference data. The B2PLYP double hybrid gives the smallest deviations, but several functionals from different rungs of the usual ladder hierarchy give mean absolute deviations below 20 kJ/mol. This includes the B97-D semilocal functional, the PBE0* global hybrid with 15% exact-exchange admixture, as well as the local hybrids LH07s-SVWN and LH07t-SVWN. Several further functionals achieve mean absolute errors below 30 kJ/mol (M06L-D4, SSB-D, B97-1-D4, LC-ωPBE-D4, LH12ct-SsirPW92-D4, LH12ct-SsifPW92-D4, LH14t-calPBE-D4, LHJ-HFcal-D4, and several further double hybrids) and thereby also still overall outperform CCSD(T) or uncorrected CASPT2. While exact-exchange admixture is a crucial factor in favoring high-spin states, the present evaluations confirm that other aspects can be important as well. A number of the better-performing functionals underestimate the spin gaps for the π-acceptor ligands but overestimate them for L = NH3, H2O. In contrast to a previous suggestion, non-self-consistent density functional theory (DFT) computations on top of Hartree-Fock orbitals are not a promising path to produce accurate spin gaps in such complexes.
Collapse
Affiliation(s)
- Marc Reimann
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
3
|
Bajaj A, Duan C, Nandy A, Taylor MG, Kulik HJ. Molecular orbital projectors in non-empirical jmDFT recover exact conditions in transition-metal chemistry. J Chem Phys 2022; 156:184112. [DOI: 10.1063/5.0089460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Low-cost, non-empirical corrections to semi-local density functional theory are essential for accurately modeling transition-metal chemistry. Here, we demonstrate the judiciously modified density functional theory (jmDFT) approach with non-empirical U and J parameters obtained directly from frontier orbital energetics on a series of transition-metal complexes. We curate a set of nine representative Ti(III) and V(IV) d1 transition-metal complexes and evaluate their flat-plane errors along the fractional spin and charge lines. We demonstrate that while jmDFT improves upon both DFT+U and semi-local DFT with the standard atomic orbital projectors (AOPs), it does so inefficiently. We rationalize these inefficiencies by quantifying hybridization in the relevant frontier orbitals. To overcome these limitations, we introduce a procedure for computing a molecular orbital projector (MOP) basis for use with jmDFT. We demonstrate this single set of d1 MOPs to be suitable for nearly eliminating all energetic delocalization error and static correlation error. In all cases, MOP jmDFT outperforms AOP jmDFT, and it eliminates most flat-plane errors at non-empirical values. Unlike DFT+U or hybrid functionals, jmDFT nearly eliminates energetic delocalization error and static correlation error within a non-empirical framework.
Collapse
Affiliation(s)
- Akash Bajaj
- Massachusetts Institute of Technology, United States of America
| | - Chenru Duan
- Massachusetts Institute of Technology, United States of America
| | - Aditya Nandy
- Massachusetts Institute of Technology, United States of America
| | | | - Heather J. Kulik
- Dept of Chemical Engineering, Massachusetts Institute of Technology, United States of America
| |
Collapse
|
4
|
On the Investigation of Effective Factors on Electronic Structure Properties of Transition Metal Complexes: Robust Modeling Using GPR Approach. INTERNATIONAL JOURNAL OF CHEMICAL ENGINEERING 2022. [DOI: 10.1155/2022/8264297] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Materials discovery is usually done using high-throughput computational screening. The use of costly and complex direct density functional theory (DFT) simulation methods has been commonly used to determine subtle trends in spin-state ordering and inorganic bonding of inorganic materials and, in general, to predict the electronic structure properties of transition metal complexes. A Gaussian process regression (GPR) framework consisting of four kernel functions is introduced for spin-state splitting estimation through inorganic chemistry-appropriate empirical inputs. To this end, the present study reviewed an extensive range of data values from earlier works. According to statistical analysis, the GPR model showed very good performance. The coefficients of determination were calculated to be 0.986 for the exponential and Matern kernel functions, suggesting the highest predictive power of these methods. Moreover, the sensitivity of output to inputs was measured. Artificial intelligence (AI) helped accurately predict the target values through various input ranges.
Collapse
|
5
|
Harper DR, Kulik HJ. Computational Scaling Relationships Predict Experimental Activity and Rate-Limiting Behavior in Homogeneous Water Oxidation. Inorg Chem 2022; 61:2186-2197. [PMID: 35037756 DOI: 10.1021/acs.inorgchem.1c03376] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
While computational screening with first-principles density functional theory (DFT) is essential for evaluating candidate catalysts, limitations in accuracy typically prevent the prediction of experimentally relevant activities. Exemplary of these challenges are homogeneous water oxidation catalysts (WOCs) where differences in experimental conditions or small changes in ligand structure can alter rate constants by over an order of magnitude. Here, we compute mechanistically relevant electronic and energetic properties for 19 mononuclear Ru transition-metal complexes (TMCs) from three experimental water oxidation catalysis studies. We discover that 15 of these TMCs have experimental activities that correlate with a single property, the ionization potential of the Ru(II)-O2 catalytic intermediate. This scaling parameter allows the quantitative understanding of activity trends and provides insight into the rate-limiting behavior. We use this approach to rationalize differences in activity with different experimental conditions, and we qualitatively analyze the source of distinct behavior for different electronic states in the other four catalysts. Comparison to closely related single-atom catalysts and modified WOCs enables rationalization of the source of rate enhancement in these WOCs.
Collapse
Affiliation(s)
- Daniel R Harper
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
6
|
Duan C, Chen S, Taylor MG, Liu F, Kulik HJ. Machine learning to tame divergent density functional approximations: a new path to consensus materials design principles. Chem Sci 2021; 12:13021-13036. [PMID: 34745533 PMCID: PMC8513898 DOI: 10.1039/d1sc03701c] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2021] [Accepted: 09/01/2021] [Indexed: 01/17/2023] Open
Abstract
Virtual high-throughput screening (VHTS) with density functional theory (DFT) and machine-learning (ML)-acceleration is essential in rapid materials discovery. By necessity, efficient DFT-based workflows are carried out with a single density functional approximation (DFA). Nevertheless, properties evaluated with different DFAs can be expected to disagree for cases with challenging electronic structure (e.g., open-shell transition-metal complexes, TMCs) for which rapid screening is most needed and accurate benchmarks are often unavailable. To quantify the effect of DFA bias, we introduce an approach to rapidly obtain property predictions from 23 representative DFAs spanning multiple families, “rungs” (e.g., semi-local to double hybrid) and basis sets on over 2000 TMCs. Although computed property values (e.g., spin state splitting and frontier orbital gap) differ by DFA, high linear correlations persist across all DFAs. We train independent ML models for each DFA and observe convergent trends in feature importance, providing DFA-invariant, universal design rules. We devise a strategy to train artificial neural network (ANN) models informed by all 23 DFAs and use them to predict properties (e.g., spin-splitting energy) of over 187k TMCs. By requiring consensus of the ANN-predicted DFA properties, we improve correspondence of computational lead compounds with literature-mined, experimental compounds over the typically employed single-DFA approach. Machine learning (ML)-based feature analysis reveals universal design rules regardless of density functional choices. Using the consensus among multiple functionals, we identify robust lead complexes in ML-accelerated chemical discovery.![]()
Collapse
Affiliation(s)
- Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology Cambridge MA 02139 USA +1-617-253-4584.,Department of Chemistry, Massachusetts Institute of Technology Cambridge MA 02139 USA
| | - Shuxin Chen
- Department of Chemical Engineering, Massachusetts Institute of Technology Cambridge MA 02139 USA +1-617-253-4584.,Department of Chemistry, Massachusetts Institute of Technology Cambridge MA 02139 USA
| | - Michael G Taylor
- Department of Chemical Engineering, Massachusetts Institute of Technology Cambridge MA 02139 USA +1-617-253-4584
| | - Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology Cambridge MA 02139 USA +1-617-253-4584
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology Cambridge MA 02139 USA +1-617-253-4584
| |
Collapse
|
7
|
Nandy A, Duan C, Taylor MG, Liu F, Steeves AH, Kulik HJ. Computational Discovery of Transition-metal Complexes: From High-throughput Screening to Machine Learning. Chem Rev 2021; 121:9927-10000. [PMID: 34260198 DOI: 10.1021/acs.chemrev.1c00347] [Citation(s) in RCA: 70] [Impact Index Per Article: 23.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Transition-metal complexes are attractive targets for the design of catalysts and functional materials. The behavior of the metal-organic bond, while very tunable for achieving target properties, is challenging to predict and necessitates searching a wide and complex space to identify needles in haystacks for target applications. This review will focus on the techniques that make high-throughput search of transition-metal chemical space feasible for the discovery of complexes with desirable properties. The review will cover the development, promise, and limitations of "traditional" computational chemistry (i.e., force field, semiempirical, and density functional theory methods) as it pertains to data generation for inorganic molecular discovery. The review will also discuss the opportunities and limitations in leveraging experimental data sources. We will focus on how advances in statistical modeling, artificial intelligence, multiobjective optimization, and automation accelerate discovery of lead compounds and design rules. The overall objective of this review is to showcase how bringing together advances from diverse areas of computational chemistry and computer science have enabled the rapid uncovering of structure-property relationships in transition-metal chemistry. We aim to highlight how unique considerations in motifs of metal-organic bonding (e.g., variable spin and oxidation state, and bonding strength/nature) set them and their discovery apart from more commonly considered organic molecules. We will also highlight how uncertainty and relative data scarcity in transition-metal chemistry motivate specific developments in machine learning representations, model training, and in computational chemistry. Finally, we will conclude with an outlook of areas of opportunity for the accelerated discovery of transition-metal complexes.
Collapse
Affiliation(s)
- Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Michael G Taylor
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Adam H Steeves
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
8
|
Vennelakanti V, Nandy A, Kulik HJ. The Effect of Hartree-Fock Exchange on Scaling Relations and Reaction Energetics for C–H Activation Catalysts. Top Catal 2021. [DOI: 10.1007/s11244-021-01482-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
|
9
|
Duan C, Liu F, Nandy A, Kulik HJ. Putting Density Functional Theory to the Test in Machine-Learning-Accelerated Materials Discovery. J Phys Chem Lett 2021; 12:4628-4637. [PMID: 33973793 DOI: 10.1021/acs.jpclett.1c00631] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
Accelerated discovery with machine learning (ML) has begun to provide the advances in efficiency needed to overcome the combinatorial challenge of computational materials design. Nevertheless, ML-accelerated discovery both inherits the biases of training data derived from density functional theory (DFT) and leads to many attempted calculations that are doomed to fail. Many compelling functional materials and catalytic processes involve strained chemical bonds, open-shell radicals and diradicals, or metal-organic bonds to open-shell transition-metal centers. Although promising targets, these materials present unique challenges for electronic structure methods and combinatorial challenges for their discovery. In this Perspective, we describe the advances needed in accuracy, efficiency, and approach beyond what is typical in conventional DFT-based ML workflows. These challenges have begun to be addressed through ML models trained to predict the results of multiple methods or the differences between them, enabling quantitative sensitivity analysis. For DFT to be trusted for a given data point in a high-throughput screen, it must pass a series of tests. ML models that predict the likelihood of calculation success and detect the presence of strong correlation will enable rapid diagnoses and adaptation strategies. These "decision engines" represent the first steps toward autonomous workflows that avoid the need for expert determination of the robustness of DFT-based materials discoveries.
Collapse
Affiliation(s)
- Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
10
|
Mariano LA, Vlaisavljevich B, Poloni R. Improved Spin-State Energy Differences of Fe(II) Molecular and Crystalline Complexes via the Hubbard U-Corrected Density. J Chem Theory Comput 2021; 17:2807-2816. [PMID: 33831303 DOI: 10.1021/acs.jctc.1c00034] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We recently showed that the DFT+U approach with a linear-response U yields adiabatic energy differences biased toward high spin [Mariano et al. J. Chem. Theory Comput. 2020, 16, 6755-6762]. Such bias is removed here by employing a density-corrected DFT approach where the PBE functional is evaluated on the Hubbard U-corrected density. The adiabatic energy differences of six Fe(II) molecular complexes computed using this approach, named PBE[U] here, are in excellent agreement with coupled cluster-corrected CASPT2 values for both weak- and strong-field ligands resulting in a mean absolute error (MAE) of 0.44 eV, smaller than that of the recently proposed Hartree-Fock density-corrected DFT (1.22 eV) and any other tested functional, including the best performer TPSSh (0.49 eV). We take advantage of the computational efficiency of this approach and compute the adiabatic energy differences of five molecular crystals using PBE[U] with periodic boundary conditions. The results show, again, an excellent agreement (MAE = 0.07 eV) with experimentally extracted values and a superior performance compared with the best performers M06-L (MAE = 0.08 eV) and TPSSh (MAE = 0.31 eV) computed on molecular fragments.
Collapse
Affiliation(s)
- Lorenzo A Mariano
- University Grenoble Alpes, CNRS, Grenoble-INP, SIMaP, F-38042 Grenoble, France
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| | - Roberta Poloni
- University Grenoble Alpes, CNRS, Grenoble-INP, SIMaP, F-38042 Grenoble, France
| |
Collapse
|
11
|
Bajaj A, Kulik HJ. Molecular DFT+U: A Transferable, Low-Cost Approach to Eliminate Delocalization Error. J Phys Chem Lett 2021; 12:3633-3640. [PMID: 33826346 DOI: 10.1021/acs.jpclett.1c00796] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
While density functional theory (DFT) is widely applied for its combination of cost and accuracy, corrections (e.g., DFT+U) that improve it are often needed to tackle correlated transition-metal chemistry. In principle, the functional form of DFT+U, consisting of a set of localized atomic orbitals (AOs) and a quadratic energy penalty for deviation from integer occupations of those AOs, enables the recovery of the exact conditions of piecewise linearity and the derivative discontinuity. Nevertheless, for practical transition-metal complexes, where both atomic states and ligand orbitals participate in bonding, standard DFT+U can fail to eliminate delocalization error (DE). Here, we show that by introducing an alternative valence-state (i.e., molecular orbital or MO) basis to the DFT+U approach, we recover exact conditions in cases for which standard DFT+U corrections have no error-reducing effect. This MO-based DFT+U also eliminates DE where standard AO-based DFT+U is already successful. We demonstrate the transferability of our approach on representative transition-metal complexes with a range of ligand field strengths, electron configurations (i.e., from Sc to Zn), and spin states.
Collapse
Affiliation(s)
- Akash Bajaj
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
12
|
Yokogawa D, Suda K. Electrostatic Potential Fitting Method Using Constrained Spatial Electron Density Expanded with Preorthogonal Natural Atomic Orbitals. J Phys Chem A 2020; 124:9665-9673. [DOI: 10.1021/acs.jpca.0c07425] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Daisuke Yokogawa
- Graduate School of Arts and Science, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Kayo Suda
- Graduate School of Arts and Science, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| |
Collapse
|
13
|
Taylor MG, Yang T, Lin S, Nandy A, Janet JP, Duan C, Kulik HJ. Seeing Is Believing: Experimental Spin States from Machine Learning Model Structure Predictions. J Phys Chem A 2020; 124:3286-3299. [PMID: 32223165 PMCID: PMC7311053 DOI: 10.1021/acs.jpca.0c01458] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
![]()
Determination of ground-state spins
of open-shell transition-metal
complexes is critical to understanding catalytic and materials properties
but also challenging with approximate electronic structure methods.
As an alternative approach, we demonstrate how structure alone can
be used to guide assignment of ground-state spin from experimentally
determined crystal structures of transition-metal complexes. We first
identify the limits of distance-based heuristics from distributions
of metal–ligand bond lengths of over 2000 unique mononuclear
Fe(II)/Fe(III) transition-metal complexes. To overcome these limits,
we employ artificial neural networks (ANNs) to predict spin-state-dependent
metal–ligand bond lengths and classify experimental ground-state
spins based on agreement of experimental structures with the ANN predictions.
Although the ANN is trained on hybrid density functional theory data,
we exploit the method-insensitivity of geometric properties to enable
assignment of ground states for the majority (ca. 80–90%) of
structures. We demonstrate the utility of the ANN by data-mining the
literature for spin-crossover (SCO) complexes, which have experimentally
observed temperature-dependent geometric structure changes, by correctly
assigning almost all (>95%) spin states in the 46 Fe(II) SCO complex
set. This approach represents a promising complement to more conventional
energy-based spin-state assignment from electronic structure theory
at the low cost of a machine learning model.
Collapse
Affiliation(s)
- Michael G Taylor
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Tzuhsiung Yang
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Sean Lin
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Jon Paul Janet
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
14
|
Nandy A, Chu DBK, Harper DR, Duan C, Arunachalam N, Cytter Y, Kulik HJ. Large-scale comparison of 3d and 4d transition metal complexes illuminates the reduced effect of exchange on second-row spin-state energetics. Phys Chem Chem Phys 2020; 22:19326-19341. [DOI: 10.1039/d0cp02977g] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
The origin of distinct 3d vs. 4d transition metal complex sensitivity to exchange is explored over a large data set.
Collapse
Affiliation(s)
- Aditya Nandy
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
- Department of Chemistry
| | - Daniel B. K. Chu
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
| | - Daniel R. Harper
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
- Department of Chemistry
| | - Chenru Duan
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
- Department of Chemistry
| | - Naveen Arunachalam
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
| | - Yael Cytter
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
| | - Heather J. Kulik
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
| |
Collapse
|
15
|
Farcaș AA, Bende A. Improving the Light-Induced Spin Transition Efficiency in Ni(II)-Based Macrocyclic-Ligand Complexes. Molecules 2019; 24:molecules24234249. [PMID: 31766599 PMCID: PMC6930591 DOI: 10.3390/molecules24234249] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2019] [Revised: 11/15/2019] [Accepted: 11/20/2019] [Indexed: 01/29/2023] Open
Abstract
The structural stability and photoabsorption properties of Ni(II)-based metal-organic complexes with octahedral coordination having different planar ligand ring structures were investigated employing density functional theory (DFT) and its time-dependent extension (TD-DFT) considering the M06 exchange-correlation functional and the Def2-TZVP basis set. The results showed that the molecular composition of different planar cyclic ligand structures had significant influences on the structural stability and photoabsorption properties of metal-organic complexes. Only those planar ligands that contained aromatic rings met the basic criteria (thermal stability, structural reversibility, and appropriate excitation frequency domain) for light-induced excited spin state trapping, but their spin transition efficiencies were very different. While, in all three aromatic cases, the singlet electronic excitations induced charge distribution that could help in the singlet-to-triplet spin transition, and triplet excitations, which could assist in the backward (triplet-to-singlet) spin transition, was found only for one complex.
Collapse
Affiliation(s)
- Alex-Adrian Farcaș
- National Institute for Research and Development of Isotopic and Molecular Technologies, Donat Street, No. 67-103, Ro-400293 Cluj-Napoca, Romania;
- Faculty of Physics, “Babeş-Bolyai” University, Mihail Kogalniceanu Street No. 1, Ro-400084 Cluj-Napoca, Romania
| | - Attila Bende
- National Institute for Research and Development of Isotopic and Molecular Technologies, Donat Street, No. 67-103, Ro-400293 Cluj-Napoca, Romania;
- Correspondence:
| |
Collapse
|
16
|
Liu F, Kulik HJ. Impact of Approximate DFT Density Delocalization Error on Potential Energy Surfaces in Transition Metal Chemistry. J Chem Theory Comput 2019; 16:264-277. [DOI: 10.1021/acs.jctc.9b00842] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
17
|
Moltved KA, Kepp KP. Performance of Density Functional Theory for Transition Metal Oxygen Bonds. Chemphyschem 2019; 20:3210-3220. [DOI: 10.1002/cphc.201900862] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2019] [Revised: 10/01/2019] [Indexed: 11/09/2022]
Affiliation(s)
- Klaus A. Moltved
- Technical University of DenmarkDTU Chemistry, Building 206, 2800 Kgs. Lyngby DK – Denmark
| | - Kasper P. Kepp
- Technical University of DenmarkDTU Chemistry, Building 206, 2800 Kgs. Lyngby DK – Denmark
| |
Collapse
|
18
|
Pandharkar R, Hermes MR, Cramer CJ, Gagliardi L. Spin-State Ordering in Metal-Based Compounds Using the Localized Active Space Self-Consistent Field Method. J Phys Chem Lett 2019; 10:5507-5513. [PMID: 31429583 DOI: 10.1021/acs.jpclett.9b02077] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Quantitatively accurate calculations for spin-state ordering in transition-metal complexes typically demand a robust multiconfigurational treatment. The poor scaling of such methods with increasing size makes them impractical for large, strongly correlated systems. Density matrix embedding theory (DMET) is a fragmentation approach that can be used to specifically address this challenge. The single-determinantal bath framework of DMET is applicable in many situations, but it has been shown to perform poorly for molecules characterized by strong correlation when a multiconfigurational self-consistent field solver is used. To ameliorate this problem, the localized active space self-consistent field (LASSCF) method was recently described. In this work, LASSCF is applied to predict spin-state energetics in mono- and di-iron systems, and we show that the model offers an accuracy equivalent to that of CASSCF but at a substantially lower computational cost. Performance as a function of basis set and active space is also examined.
Collapse
Affiliation(s)
- Riddhish Pandharkar
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute , University of Minnesota , Minneapolis , Minnesota 55455 , United States
| | - Matthew R Hermes
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute , University of Minnesota , Minneapolis , Minnesota 55455 , United States
| | - Christopher J Cramer
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute , University of Minnesota , Minneapolis , Minnesota 55455 , United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute , University of Minnesota , Minneapolis , Minnesota 55455 , United States
| |
Collapse
|
19
|
Zhao Q, Kulik HJ. Stable Surfaces That Bind Too Tightly: Can Range-Separated Hybrids or DFT+U Improve Paradoxical Descriptions of Surface Chemistry? J Phys Chem Lett 2019; 10:5090-5098. [PMID: 31411023 PMCID: PMC6748670 DOI: 10.1021/acs.jpclett.9b01650] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Accepted: 08/14/2019] [Indexed: 05/25/2023]
Abstract
Approximate, semilocal density functional theory (DFT) suffers from delocalization error that can lead to a paradoxical model of catalytic surfaces that both overbind adsorbates yet are also too stable. We investigate the effect of two widely applied approaches for delocalization error correction, (i) affordable DFT+U (i.e., semilocal DFT augmented with a Hubbard U) and (ii) hybrid functionals with an admixture of Hartree-Fock (HF) exchange, on surface and adsorbate energies across a range of rutile transition metal oxides widely studied for their promise as water-splitting catalysts. We observe strongly row- and period-dependent trends with DFT+U, which increases surface formation energies only in early transition metals (e.g., Ti and V) and decreases adsorbate energies only in later transition metals (e.g., Ir and Pt). Both global and local hybrids destabilize surfaces and reduce adsorbate binding across the periodic table, in agreement with higher-level reference calculations. Density analysis reveals why hybrid functionals correct both quantities, whereas DFT+U does not. We recommend local, range-separated hybrids for the accurate modeling of catalysis in transition metal oxides at only a modest increase in computational cost over semilocal DFT.
Collapse
Affiliation(s)
- Qing Zhao
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, Cambridge, Massachusetts 02139, United States
- Department
of Mechanical Engineering, Massachusetts
Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
20
|
Kulik HJ. Making machine learning a useful tool in the accelerated discovery of transition metal complexes. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1439] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Heather J. Kulik
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge Massachusetts
| |
Collapse
|
21
|
Francés‐Monerris A, Gros PC, Assfeld X, Monari A, Pastore M. Toward Luminescent Iron Complexes: Unravelling the Photophysics by Computing Potential Energy Surfaces. CHEMPHOTOCHEM 2019. [DOI: 10.1002/cptc.201900100] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Antonio Francés‐Monerris
- Laboratoire de Physique et Chimie Théoriques (LPCT)Université de Lorraine, CNRS 54000 Nancy France
| | - Philippe C. Gros
- Laboratoire Lorrain de Chimie Moléculaire (L2CM)Université de Lorraine, CNRS 54000 Nancy France
| | - Xavier Assfeld
- Laboratoire de Physique et Chimie Théoriques (LPCT)Université de Lorraine, CNRS 54000 Nancy France
| | - Antonio Monari
- Laboratoire de Physique et Chimie Théoriques (LPCT)Université de Lorraine, CNRS 54000 Nancy France
| | - Mariachiara Pastore
- Laboratoire de Physique et Chimie Théoriques (LPCT)Université de Lorraine, CNRS 54000 Nancy France
| |
Collapse
|
22
|
Katari M, Carmichael D, Jacquemin D, Frison G. Structure of Electronically Reduced N-Donor Bidentate Ligands and Their Heteroleptic Four-Coordinate Zinc Complexes: A Survey of Density Functional Theory Results. Inorg Chem 2019; 58:7169-7179. [PMID: 31117621 DOI: 10.1021/acs.inorgchem.8b03549] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The role of Hartree-Fock exchange in describing the structural changes occurring upon reduction of bipyridine-based ligands and their complexes is investigated within the framework of density functional theory (DFT) calculations. A set of four free ligands in their neutral and radical anionic forms, and two of their zinc complexes in their dicationic and monocationic radical forms, is used to compare a large panel of pure, conventional, and long-range corrected hybrid DFT functionals; coupled cluster single and double calculations are used alongside experimental results as benchmarks. Particular attention has been devoted to the magnitude of the change, upon reduction, of the Δ-parameter, which measures the difference between the Cpy-Cpy and the C-N bond lengths in bipyridine ligand and is known to experimentally correlate with the charge of the ligands. Our results indicate that the structural changes significantly depend on the amount of exact exchange included in the functional. A progressive evolution is observed for the free ligands, whereas two distinct sets of results are obtained for the complexes. Functionals with a small degree of HF exchange, e.g., B3LYP, do not adequately describe geometric changes for the considered species, and, quite surprisingly, the same holds for the CC2 method. The best agreement to experimental and CCSD values is obtained with functionals that include a significant but not excessive part of exact exchange, e.g., CAM-B3LYP, M06-2X, and ωB97X-D. The calculated localization of the added electron after reduction, which depends on the self-interaction error, is used to rationalize these outcomes. Static correlation is also shown to play a role in the accurate description of the electronic structure.
Collapse
Affiliation(s)
| | - Duncan Carmichael
- LCM, CNRS, Ecole Polytechnique , IP Paris , F-91128 Palaiseau , France
| | - Denis Jacquemin
- University of Nantes , CNRS, CEISAM (UMR 6230), 2 chemin de la Houssinière , 44322 Nantes , Cedex 03 , France
| | - Gilles Frison
- LCM, CNRS, Ecole Polytechnique , IP Paris , F-91128 Palaiseau , France
| |
Collapse
|
23
|
Bajaj A, Liu F, Kulik HJ. Non-empirical, low-cost recovery of exact conditions with model-Hamiltonian inspired expressions in jmDFT. J Chem Phys 2019; 150:154115. [PMID: 31005112 DOI: 10.1063/1.5091563] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Density functional theory (DFT) is widely applied to both molecules and materials, but well known energetic delocalization and static correlation errors in practical exchange-correlation approximations limit quantitative accuracy. Common methods that correct energetic delocalization errors, such as the Hubbard U correction in DFT+U or Hartree-Fock exchange in global hybrids, do so at the cost of worsening static correlation errors. We recently introduced an alternate approach [Bajaj et al., J. Chem. Phys. 147, 191101 (2017)] known as judiciously modified DFT (jmDFT), wherein the deviation from exact behavior of semilocal functionals over both fractional spin and charge, i.e., the so-called flat plane, was used to motivate functional forms of second order analytic corrections. In this work, we introduce fully nonempirical expressions for all four coefficients in a DFT+U+J-inspired form of jmDFT, where all coefficients are obtained only from energies and eigenvalues of the integer-electron systems. We show good agreement for U and J coefficients obtained nonempirically as compared with the results of numerical fitting in a jmDFT U+J/J' correction. Incorporating the fully nonempirical jmDFT correction reduces and even eliminates the fractional spin error at the same time as eliminating the energetic delocalization error. We show that this approach extends beyond s-electron systems to higher angular momentum cases including p- and d-electrons. Finally, we diagnose some shortcomings of the current jmDFT approach that limit its ability to improve upon DFT results for cases such as weakly bound anions due to poor underlying semilocal functional behavior.
Collapse
Affiliation(s)
- Akash Bajaj
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
24
|
Liu F, Yang T, Yang J, Xu E, Bajaj A, Kulik HJ. Bridging the Homogeneous-Heterogeneous Divide: Modeling Spin for Reactivity in Single Atom Catalysis. Front Chem 2019; 7:219. [PMID: 31041303 PMCID: PMC6476907 DOI: 10.3389/fchem.2019.00219] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Accepted: 03/20/2019] [Indexed: 12/03/2022] Open
Abstract
Single atom catalysts (SACs) are emergent catalytic materials that have the promise of merging the scalability of heterogeneous catalysts with the high activity and atom economy of homogeneous catalysts. Computational, first-principles modeling can provide essential insight into SAC mechanism and active site configuration, where the sub-nm-scale environment can challenge even the highest-resolution experimental spectroscopic techniques. Nevertheless, the very properties that make SACs attractive in catalysis, such as localized d electrons of the isolated transition metal center, make them challenging to study with conventional computational modeling using density functional theory (DFT). For example, Fe/N-doped graphitic SACs have exhibited spin-state dependent reactivity that remains poorly understood. However, spin-state ordering in DFT is very sensitive to the nature of the functional approximation chosen. In this work, we develop accurate benchmarks from correlated wavefunction theory (WFT) for relevant octahedral complexes. We use those benchmarks to evaluate optimal DFT functional choice for predicting spin state ordering in small octahedral complexes as well as models of pyridinic and pyrrolic nitrogen environments expected in larger SACs. Using these guidelines, we determine Fe/N-doped graphene SAC model properties and reactivity as well as their sensitivities to DFT functional choice. Finally, we conclude with broad recommendations for computational modeling of open-shell transition metal single-atom catalysts.
Collapse
Affiliation(s)
- Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Tzuhsiung Yang
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Jing Yang
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Eve Xu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Akash Bajaj
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States.,Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States
| |
Collapse
|
25
|
Duan C, Janet JP, Liu F, Nandy A, Kulik HJ. Learning from Failure: Predicting Electronic Structure Calculation Outcomes with Machine Learning Models. J Chem Theory Comput 2019; 15:2331-2345. [DOI: 10.1021/acs.jctc.9b00057] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
|
26
|
Janet JP, Liu F, Nandy A, Duan C, Yang T, Lin S, Kulik HJ. Designing in the Face of Uncertainty: Exploiting Electronic Structure and Machine Learning Models for Discovery in Inorganic Chemistry. Inorg Chem 2019; 58:10592-10606. [PMID: 30834738 DOI: 10.1021/acs.inorgchem.9b00109] [Citation(s) in RCA: 59] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Recent transformative advances in computing power and algorithms have made computational chemistry central to the discovery and design of new molecules and materials. First-principles simulations are increasingly accurate and applicable to large systems with the speed needed for high-throughput computational screening. Despite these strides, the combinatorial challenges associated with the vastness of chemical space mean that more than just fast and accurate computational tools are needed for accelerated chemical discovery. In transition-metal chemistry and catalysis, unique challenges arise. The variable spin, oxidation state, and coordination environments favored by elements with well-localized d or f electrons provide great opportunity for tailoring properties in catalytic or functional (e.g., magnetic) materials but also add layers of uncertainty to any design strategy. We outline five key mandates for realizing computationally driven accelerated discovery in inorganic chemistry: (i) fully automated simulation of new compounds, (ii) knowledge of prediction sensitivity or accuracy, (iii) faster-than-fast property prediction methods, (iv) maps for rapid chemical space traversal, and (v) a means to reveal design rules on the kilocompound scale. Through case studies in open-shell transition-metal chemistry, we describe how advances in methodology and software in each of these areas bring about new chemical insights. We conclude with our outlook on the next steps in this process toward realizing fully autonomous discovery in inorganic chemistry using computational chemistry.
Collapse
Affiliation(s)
- Jon Paul Janet
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Fang Liu
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Aditya Nandy
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States.,Department of Chemistry , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Chenru Duan
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States.,Department of Chemistry , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Tzuhsiung Yang
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Sean Lin
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Heather J Kulik
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| |
Collapse
|
27
|
Nandy A, Duan C, Janet JP, Gugler S, Kulik HJ. Strategies and Software for Machine Learning Accelerated Discovery in Transition Metal Chemistry. Ind Eng Chem Res 2018. [DOI: 10.1021/acs.iecr.8b04015] [Citation(s) in RCA: 76] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Jon Paul Janet
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Stefan Gugler
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
28
|
Siig OS, Kepp KP. Iron(II) and Iron(III) Spin Crossover: Toward an Optimal Density Functional. J Phys Chem A 2018; 122:4208-4217. [PMID: 29630380 DOI: 10.1021/acs.jpca.8b02027] [Citation(s) in RCA: 63] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Spin crossover (SCO) plays a major role in biochemistry, catalysis, materials, and emerging technologies such as molecular electronics and sensors, and thus accurate prediction and design of SCO systems is of high priority. However, the main tool for this purpose, density functional theory (DFT), is very sensitive to applied methodology. The most abundant SCO systems are Fe(II) and Fe(III) systems. Even with average good agreement, a functional may be significantly more accurate for Fe(II) or Fe(III) systems, preventing balanced study of SCO candidates of both types. The present work investigates DFT's performance for well-known Fe(II) and Fe(III) SCO complexes, using various design types and customized versions of GGA, hybrid, meta-GGA, meta-hybrid, double-hybrid, and long-range-corrected hybrid functionals. We explore the limits of DFT performance and identify proficient Fe(II)-Fe(III)-balanced functionals. We identify and quantify remarkable differences in the DFT description of Fe(II) and Fe(III) systems. Most functionals become more accurate once Hartree-Fock exchange is adjusted to 10-17%, regardless of the type of functionals involved. However, this typically introduces a clear Fe(II)-Fe(III) bias. The most accurate functionals measured by mean absolute errors <10 kJ/mol are CAMB3LYP-17, B3LYP*, and B97-15 with 15-17% Hartree-Fock exchange, closely followed by CAMB3LYP and CAMB3LYP-15, OPBE, rPBE-10, and B3P86-15. While GGA functionals display a small Fe(II)-Fe(III) bias, they are generally inaccurate, except the O exchange functional. Hybrid functionals (including B2PLYP double hybrids and meta hybrids) tend to favor HS too much in Fe(II) vs Fe(III), which is important in many studies where the oxidation state of iron can vary, e.g. rational SCO design and studies of catalytic processes involving iron. The only functional with a combined bias <5 kJ/mol and a decent MAE (15 kJ/mol) is our customized PBE0-12 functional. Alternatively one has to sacrifice Fe(II)-Fe(III) balance to use the best functionals for each group separately. We also investigated the precision (measured as the standard deviation of errors) and show that the target accuracy for iron SCO is 10 kJ/mol for accuracy and 5 kJ/mol for precision, and DFT is probably not going to break this limit in the near future. Importantly, all four types of functional behavior (accurate/precise, accurate/imprecise, inaccurate/precise, inaccurate/imprecise) are observed. More generally, our work illustrates the importance not only of overall accuracy but also of balanced accuracy for systems likely to occur in context.
Collapse
Affiliation(s)
- Oliver S Siig
- DTU Chemistry , Technical University of Denmark , Building 206 , 2800 Kgs. Lyngby , Denmark
| | - Kasper P Kepp
- DTU Chemistry , Technical University of Denmark , Building 206 , 2800 Kgs. Lyngby , Denmark
| |
Collapse
|
29
|
Song S, Kim MC, Sim E, Benali A, Heinonen O, Burke K. Benchmarks and Reliable DFT Results for Spin Gaps of Small Ligand Fe(II) Complexes. J Chem Theory Comput 2018; 14:2304-2311. [DOI: 10.1021/acs.jctc.7b01196] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- Suhwan Song
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Min-Cheol Kim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Eunji Sim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | | | | | - Kieron Burke
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| |
Collapse
|
30
|
Janet JP, Chan L, Kulik HJ. Accelerating Chemical Discovery with Machine Learning: Simulated Evolution of Spin Crossover Complexes with an Artificial Neural Network. J Phys Chem Lett 2018; 9:1064-1071. [PMID: 29425453 DOI: 10.1021/acs.jpclett.8b00170] [Citation(s) in RCA: 100] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Machine learning (ML) has emerged as a powerful complement to simulation for materials discovery by reducing time for evaluation of energies and properties at accuracy competitive with first-principles methods. We use genetic algorithm (GA) optimization to discover unconventional spin-crossover complexes in combination with efficient scoring from an artificial neural network (ANN) that predicts spin-state splitting of inorganic complexes. We explore a compound space of over 5600 candidate materials derived from eight metal/oxidation state combinations and a 32-ligand pool. We introduce a strategy for error-aware ML-driven discovery by limiting how far the GA travels away from the nearest ANN training points while maximizing property (i.e., spin-splitting) fitness, leading to discovery of 80% of the leads from full chemical space enumeration. Over a 51-complex subset, average unsigned errors (4.5 kcal/mol) are close to the ANN's baseline 3 kcal/mol error. By obtaining leads from the trained ANN within seconds rather than days from a DFT-driven GA, this strategy demonstrates the power of ML for accelerating inorganic material discovery.
Collapse
Affiliation(s)
- Jon Paul Janet
- Department of Chemical Engineering, Massachusetts Institute of Technology , Cambridge, Massachusetts 02139, United States
| | - Lydia Chan
- Department of Chemical Engineering, Massachusetts Institute of Technology , Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology , Cambridge, Massachusetts 02139, United States
| |
Collapse
|
31
|
Wilbraham L, Adamo C, Ciofini I. Communication: Evaluating non-empirical double hybrid functionals for spin-state energetics in transition-metal complexes. J Chem Phys 2018; 148:041103. [DOI: 10.1063/1.5019641] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Liam Wilbraham
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris (IRCP), 11 rue P. et M. Curie, F-75005 Paris, France
| | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris (IRCP), 11 rue P. et M. Curie, F-75005 Paris, France
| | - Ilaria Ciofini
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris (IRCP), 11 rue P. et M. Curie, F-75005 Paris, France
| |
Collapse
|
32
|
Zhao Q, Kulik HJ. Where Does the Density Localize in the Solid State? Divergent Behavior for Hybrids and DFT+U. J Chem Theory Comput 2018; 14:670-683. [PMID: 29298057 DOI: 10.1021/acs.jctc.7b01061] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Approximate density functional theory (DFT) is widely used in chemistry and physics, despite delocalization errors that affect energetic and density properties. DFT+U (i.e., semilocal DFT augmented with a Hubbard U correction) and global hybrid functionals are two commonly employed practical methods to address delocalization error. Recent work demonstrated that in transition-metal complexes both methods localize density away from the metal and onto surrounding ligands, regardless of metal or ligand identity. In this work, we compare density localization trends with DFT+U and global hybrids on a diverse set of 34 transition-metal-containing solids with varying magnetic state, electron configuration and valence shell, and coordinating-atom orbital diffuseness (i.e., O, S, Se). We also study open-framework solids in which the metal is coordinated by molecular ligands, i.e., MCO3, M(OH)2, M(NCNH)2, K3M(CN)6 (M = V-Ni). As in transition-metal complexes, incorporation of Hartree-Fock exchange consistently localizes density away from the metal, but DFT+U exhibits diverging behavior, localizing density (i) onto the metal in low-spin and late transition metals and (ii) away from the metal in other cases in agreement with hybrids. To isolate the effect of the crystal environment, we extract molecular analogues from open-framework transition-metal solids and observe consistent localization of the density away from the metal in all cases with both DFT+U and hybrid exchange. These observations highlight the limited applicability of trends established for functional tuning on transition-metal complexes even to equivalent coordination environments in the solid state.
Collapse
Affiliation(s)
- Qing Zhao
- Department of Chemical Engineering and ‡Department of Mechanical Engineering, Massachusetts Institute of Technology , Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering and ‡Department of Mechanical Engineering, Massachusetts Institute of Technology , Cambridge, Massachusetts 02139, United States
| |
Collapse
|
33
|
Prokopiou G, Kronik L. Spin-State Energetics of Fe Complexes from an Optimally Tuned Range-Separated Hybrid Functional. Chemistry 2017; 24:5173-5182. [DOI: 10.1002/chem.201704014] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2017] [Revised: 10/05/2017] [Indexed: 01/20/2023]
Affiliation(s)
- Georgia Prokopiou
- Department of Materials and Interfaces; Weizmann Institute of Science; Rehovoth 76100 Israel
| | - Leeor Kronik
- Department of Materials and Interfaces; Weizmann Institute of Science; Rehovoth 76100 Israel
| |
Collapse
|
34
|
Janet JP, Kulik HJ. Resolving Transition Metal Chemical Space: Feature Selection for Machine Learning and Structure-Property Relationships. J Phys Chem A 2017; 121:8939-8954. [PMID: 29095620 DOI: 10.1021/acs.jpca.7b08750] [Citation(s) in RCA: 121] [Impact Index Per Article: 17.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Machine learning (ML) of quantum mechanical properties shows promise for accelerating chemical discovery. For transition metal chemistry where accurate calculations are computationally costly and available training data sets are small, the molecular representation becomes a critical ingredient in ML model predictive accuracy. We introduce a series of revised autocorrelation functions (RACs) that encode relationships of the heuristic atomic properties (e.g., size, connectivity, and electronegativity) on a molecular graph. We alter the starting point, scope, and nature of the quantities evaluated in standard ACs to make these RACs amenable to inorganic chemistry. On an organic molecule set, we first demonstrate superior standard AC performance to other presently available topological descriptors for ML model training, with mean unsigned errors (MUEs) for atomization energies on set-aside test molecules as low as 6 kcal/mol. For inorganic chemistry, our RACs yield 1 kcal/mol ML MUEs on set-aside test molecules in spin-state splitting in comparison to 15-20× higher errors for feature sets that encode whole-molecule structural information. Systematic feature selection methods including univariate filtering, recursive feature elimination, and direct optimization (e.g., random forest and LASSO) are compared. Random-forest- or LASSO-selected subsets 4-5× smaller than the full RAC set produce sub- to 1 kcal/mol spin-splitting MUEs, with good transferability to metal-ligand bond length prediction (0.004-5 Å MUE) and redox potential on a smaller data set (0.2-0.3 eV MUE). Evaluation of feature selection results across property sets reveals the relative importance of local, electronic descriptors (e.g., electronegativity, atomic number) in spin-splitting and distal, steric effects in redox potential and bond lengths.
Collapse
Affiliation(s)
- Jon Paul Janet
- Department of Chemical Engineering, Massachusetts Institute of Technology , Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology , Cambridge, Massachusetts 02139, United States
| |
Collapse
|
35
|
Gani TZH, Kulik HJ. Unifying Exchange Sensitivity in Transition-Metal Spin-State Ordering and Catalysis through Bond Valence Metrics. J Chem Theory Comput 2017; 13:5443-5457. [DOI: 10.1021/acs.jctc.7b00848] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Terry Z. H. Gani
- Department
of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department
of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
36
|
Janet JP, Kulik HJ. Predicting electronic structure properties of transition metal complexes with neural networks. Chem Sci 2017; 8:5137-5152. [PMID: 30155224 PMCID: PMC6100542 DOI: 10.1039/c7sc01247k] [Citation(s) in RCA: 109] [Impact Index Per Article: 15.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2017] [Accepted: 05/09/2017] [Indexed: 12/24/2022] Open
Abstract
High-throughput computational screening has emerged as a critical component of materials discovery. Direct density functional theory (DFT) simulation of inorganic materials and molecular transition metal complexes is often used to describe subtle trends in inorganic bonding and spin-state ordering, but these calculations are computationally costly and properties are sensitive to the exchange-correlation functional employed. To begin to overcome these challenges, we trained artificial neural networks (ANNs) to predict quantum-mechanically-derived properties, including spin-state ordering, sensitivity to Hartree-Fock exchange, and spin-state specific bond lengths in transition metal complexes. Our ANN is trained on a small set of inorganic-chemistry-appropriate empirical inputs that are both maximally transferable and do not require precise three-dimensional structural information for prediction. Using these descriptors, our ANN predicts spin-state splittings of single-site transition metal complexes (i.e., Cr-Ni) at arbitrary amounts of Hartree-Fock exchange to within 3 kcal mol-1 accuracy of DFT calculations. Our exchange-sensitivity ANN enables improved predictions on a diverse test set of experimentally-characterized transition metal complexes by extrapolation from semi-local DFT to hybrid DFT. The ANN also outperforms other machine learning models (i.e., support vector regression and kernel ridge regression), demonstrating particularly improved performance in transferability, as measured by prediction errors on the diverse test set. We establish the value of new uncertainty quantification tools to estimate ANN prediction uncertainty in computational chemistry, and we provide additional heuristics for identification of when a compound of interest is likely to be poorly predicted by the ANN. The ANNs developed in this work provide a strategy for screening transition metal complexes both with direct ANN prediction and with improved structure generation for validation with first principles simulation.
Collapse
Affiliation(s)
- Jon Paul Janet
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , MA 02139 , USA . ; Tel: +1-617-253-4584
| | - Heather J Kulik
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , MA 02139 , USA . ; Tel: +1-617-253-4584
| |
Collapse
|
37
|
Proppe J, Reiher M. Reliable Estimation of Prediction Uncertainty for Physicochemical Property Models. J Chem Theory Comput 2017; 13:3297-3317. [PMID: 28581746 DOI: 10.1021/acs.jctc.7b00235] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
One of the major challenges in computational science is to determine the uncertainty of a virtual measurement, that is the prediction of an observable based on calculations. As highly accurate first-principles calculations are in general unfeasible for most physical systems, one usually resorts to parameteric property models of observables, which require calibration by incorporating reference data. The resulting predictions and their uncertainties are sensitive to systematic errors such as inconsistent reference data, parametric model assumptions, or inadequate computational methods. Here, we discuss the calibration of property models in the light of bootstrapping, a sampling method that can be employed for identifying systematic errors and for reliable estimation of the prediction uncertainty. We apply bootstrapping to assess a linear property model linking the 57Fe Mössbauer isomer shift to the contact electron density at the iron nucleus for a diverse set of 44 molecular iron compounds. The contact electron density is calculated with 12 density functionals across Jacob's ladder (PWLDA, BP86, BLYP, PW91, PBE, M06-L, TPSS, B3LYP, B3PW91, PBE0, M06, TPSSh). We provide systematic-error diagnostics and reliable, locally resolved uncertainties for isomer-shift predictions. Pure and hybrid density functionals yield average prediction uncertainties of 0.06-0.08 mm s-1 and 0.04-0.05 mm s-1, respectively, the latter being close to the average experimental uncertainty of 0.02 mm s-1. Furthermore, we show that both model parameters and prediction uncertainty depend significantly on the composition and number of reference data points. Accordingly, we suggest that rankings of density functionals based on performance measures (e.g., the squared coefficient of correlation, r2, or the root-mean-square error, RMSE) should not be inferred from a single data set. This study presents the first statistically rigorous calibration analysis for theoretical Mössbauer spectroscopy, which is of general applicability for physicochemical property models and not restricted to isomer-shift predictions. We provide the statistically meaningful reference data set MIS39 and a new calibration of the isomer shift based on the PBE0 functional.
Collapse
Affiliation(s)
- Jonny Proppe
- Laboratorium für Physikalische Chemie, ETH Zürich , Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich , Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
38
|
Flores-Leonar MM, Moreno-Esparza R, Ugalde-Saldívar VM, Amador-Bedolla C. Correlating Properties in Iron(III) Complexes: A DFT Description of Structure, Redox Potential and Spin Crossover Phenomena. ChemistrySelect 2017. [DOI: 10.1002/slct.201700547] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Martha M. Flores-Leonar
- Facultad de Química (UNAM); Edificio B; Av. Universidad 3000, Coyoacán Ciudad de México 04510 México
| | - Rafael Moreno-Esparza
- Facultad de Química (UNAM); Edificio B; Av. Universidad 3000, Coyoacán Ciudad de México 04510 México
| | - Víctor M. Ugalde-Saldívar
- Facultad de Química (UNAM); Edificio B; Av. Universidad 3000, Coyoacán Ciudad de México 04510 México
| | - Carlos Amador-Bedolla
- Facultad de Química (UNAM); Edificio B; Av. Universidad 3000, Coyoacán Ciudad de México 04510 México
| |
Collapse
|
39
|
Wilbraham L, Verma P, Truhlar DG, Gagliardi L, Ciofini I. Multiconfiguration Pair-Density Functional Theory Predicts Spin-State Ordering in Iron Complexes with the Same Accuracy as Complete Active Space Second-Order Perturbation Theory at a Significantly Reduced Computational Cost. J Phys Chem Lett 2017; 8:2026-2030. [PMID: 28436662 DOI: 10.1021/acs.jpclett.7b00570] [Citation(s) in RCA: 45] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The spin-state orderings in nine Fe(II) and Fe(III) complexes with ligands of diverse ligand-field strength were investigated with multiconfiguration pair-density functional theory (MC-PDFT). The performance of this method was compared to that of complete active space second-order perturbation theory (CASPT2) and Kohn-Sham density functional theory. We also investigated the dependence of CASPT2 and MC-PDFT results on the size of the active-space. MC-PDFT reproduces the CASPT2 spin-state ordering, the dependence on the ligand field strength, and the dependence on active space at a computational cost that is significantly reduced as compared to CASPT2.
Collapse
Affiliation(s)
- Liam Wilbraham
- PSL Research University , Institut de Recherche de Chimie Paris IRCP, CNRS - Chimie ParisTech, 11 rue Pierre et Marie Curie, F-75005 Paris, France
| | - Pragya Verma
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota , Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota , Minneapolis, Minnesota 55455-0431, United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota , Minneapolis, Minnesota 55455-0431, United States
| | - Ilaria Ciofini
- PSL Research University , Institut de Recherche de Chimie Paris IRCP, CNRS - Chimie ParisTech, 11 rue Pierre et Marie Curie, F-75005 Paris, France
| |
Collapse
|
40
|
Janet JP, Gani TZH, Steeves AH, Ioannidis EI, Kulik HJ. Leveraging Cheminformatics Strategies for Inorganic Discovery: Application to Redox Potential Design. Ind Eng Chem Res 2017. [DOI: 10.1021/acs.iecr.7b00808] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jon Paul Janet
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Terry Z. H. Gani
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Adam H. Steeves
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Efthymios I. Ioannidis
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
41
|
Ashley DC, Jakubikova E. Ironing out the photochemical and spin-crossover behavior of Fe(II) coordination compounds with computational chemistry. Coord Chem Rev 2017. [DOI: 10.1016/j.ccr.2017.02.005] [Citation(s) in RCA: 65] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
|