Research method
Density Functional Theory (DFT)
Density functional theory approximates the electronic energy of a molecule or extended solid from its electron density. Chemists use it to rank isomers, estimate barriers, assign redox character, and generate descriptors that a statistical model can then learn. The number that comes out is a computed energy or geometry at a chosen functional and basis, not a measured rate, and 'chemical accuracy' (about 1 kcal mol⁻¹) is a target some hybrid workflows hit only after training on DFT itself.
Researchers reach for DFT when an experiment cannot easily see a transition state, a spin-state gap, or which ligand orbital is being reduced. It answers 'which structure or path is lower at this level of theory?' Its main limitation in this library is domain: octahedral first-row complexes, SNAr barriers, Au(100)–water, or a gas-phase heme model do not automatically transfer to a different metal, solvent, or applied potential.
Evidence
What the evidence shows
Drawn from 28 studies in this library. Each finding starts with a plain-language takeaway, then the denser detail. Supports means evidence for a finding; Challenges means evidence against a stated position; Qualifies marks scope with a short note on each study’s contribution. Challenged positions are labeled — they are not findings.
ML trained on DFT descriptors can beat raw DFT on a held-out SNAr set and still depend on being able to compute the transition state. A Gaussian-process model reached R² = 0.93 and MAE = 0.77 kcal mol⁻¹ versus raw DFT MAE 2.93; selectivity top-1 accuracy was 86%. Acid/base catalysis and other reaction classes were not in the model; TS geometry remains the bottleneck.
Neural nets can learn DFT spin-state gaps for a defined ligand set. On octahedral Cr–Ni complexes, test RMSE was 3.1 kcal mol⁻¹ and ground state was right in 528/538 test cases; CSD complexes were worse (MUE 10 kcal mol⁻¹ overall, 5 inside a reliability cutoff), and cyclams and early metals were outliers around 30 kcal mol⁻¹. Predictions target the B3LYP-family DFT the network was trained on.
DFT is often a supporting assignment tool, not the primary result. Periodic-DFT on a Zn phenanthroline–maleate crystal gave a 3.45 eV gap; gas-phase heme-model pathways accompanied FT-ICR epoxidation rates with estimated ±30% rate error; ligand-centered versus metal-centered redox in cobalt PY4/PY3PZ complexes was parsed with DFT after electrocatalysis.
- A new Zn phenanthroline–maleate crystal
- Gas-phase rates of heme-model olefin epoxidation
- Pyrazine reservoirs that lower cobalt HER overpotential
Study Role Design N Population Outcome A new Zn phenanthroline–maleate crystal Supports OtherCrystal structure, periodic DFT, and secondary MIC assays of a Zn–phenanthroline–maleate complex Coordination-chemistry characterization — no cohort N Zn(II) phenanthroline/maleate coordination compound Structure, electronics, and weak antibacterial activity vs S. mutans Gas-phase rates of heme-model olefin epoxidation Supports OtherGas-phase FT-ICR bimolecular rate constants for [FeIV(O)(TPFPP+·)]+ with olefins plus DFT Physical-organic kinetics panel of olefins — not a sample-N study Gas-phase iron-oxo porphyrin cation and olefin substrates Epoxidation rate constants and efficiency versus olefin ionization potential Pyrazine reservoirs that lower cobalt HER overpotential Supports OtherCobalt PY4/PY3PZ complexes for aqueous electrocatalytic and photoredox H2 evolution with DFT Molecular catalysis study — no sample N Cobalt redox-active ligand complexes in aqueous catalytic assays Overpotential and H2 evolution as a function of pendant pyrazine redox reservoirs Learned potentials trained on DFT can extend the timescale of a mechanism study. A PaiNN ensemble on PBE-D3 Au(100)–water data plus metadynamics gave an associative ORR path with an O₂-to-hydroxyl barrier of about 0.3 eV — a landscape from the neural-network potential, not a measured Tafel slope, and without explicit cations or applied potential.
Open questions
Tensions and limits
Some items are genuine disagreements on the same question. Others mark different assays, populations, or outcomes — limits on how far one study travels — not a forced fight between papers.
These papers do not agree on what DFT is for. One trains a statistical model to correct DFT barriers; one trains a net to emulate DFT spin states; others use a single-point or periodic calculation to rationalise a crystal, a gas-phase ion, or a ligand redox event. Quoting a kcal mol⁻¹ figure without saying whether it is raw DFT, a GP correction, or an NNP barrier mixes those jobs.
- ML plus transition states predict SNAr barriers
- Neural nets predict TM spin states and bonds
- NN potentials map gold–water ORR paths
- Pyrazine reservoirs that lower cobalt HER overpotential
Study Role Design N Population Outcome ML plus transition states predict SNAr barriers Supports Computational / modellingDFT/physical-organic descriptors with GPR for SNAr barriers and selectivity ML on reaction barriers/selectivity; learning curves emphasize 50–150 training samples — no single cohort N Nucleophilic aromatic substitution reactions (including patent-selectivity cases) Predicted activation free energies and regio-/chemoselectivity accuracy Neural nets predict TM spin states and bonds Supports Computational / modellingANN trained on DFT octahedral Cr–Ni complexes to predict spin states and bond lengths N=2690 · 2,690 DFT geometry optimizations after filtering spin-contaminated cases Homoleptic and heteroleptic octahedral transition-metal complexes Predicted spin-state energetics and metal–ligand bond lengths NN potentials map gold–water ORR paths Supports Computational / modellingPaiNN active-learning NNPs enabling path-CV metadynamics of ORR on Au(100)–water 2.5 ns production metadynamics after AIMD active learning — no sample N Au(100)–water models with O2/OH adsorbates Associative ORR pathway and ~0.3 eV O2-to-hydroxyl barrier Pyrazine reservoirs that lower cobalt HER overpotential Supports OtherCobalt PY4/PY3PZ complexes for aqueous electrocatalytic and photoredox H2 evolution with DFT Molecular catalysis study — no sample N Cobalt redox-active ligand complexes in aqueous catalytic assays Overpotential and H2 evolution as a function of pendant pyrazine redox reservoirs Gas-phase ions, aqueous electrocatalysis, and a crystal gap are different DFT worlds. Naked [Feᴵⱽ(O)(porphyrin)]⁺ epoxidation in FT-ICR is not P450 Compound I in water; the Zn crystal gap was not tested in a device; cobalt HER overpotential is an electrochemical measurement with DFT as interpretation.
- Gas-phase rates of heme-model olefin epoxidation
- A new Zn phenanthroline–maleate crystal
- Pyrazine reservoirs that lower cobalt HER overpotential
Study Role Design N Population Outcome Gas-phase rates of heme-model olefin epoxidation Supports OtherGas-phase FT-ICR bimolecular rate constants for [FeIV(O)(TPFPP+·)]+ with olefins plus DFT Physical-organic kinetics panel of olefins — not a sample-N study Gas-phase iron-oxo porphyrin cation and olefin substrates Epoxidation rate constants and efficiency versus olefin ionization potential A new Zn phenanthroline–maleate crystal Supports OtherCrystal structure, periodic DFT, and secondary MIC assays of a Zn–phenanthroline–maleate complex Coordination-chemistry characterization — no cohort N Zn(II) phenanthroline/maleate coordination compound Structure, electronics, and weak antibacterial activity vs S. mutans Pyrazine reservoirs that lower cobalt HER overpotential Supports OtherCobalt PY4/PY3PZ complexes for aqueous electrocatalytic and photoredox H2 evolution with DFT Molecular catalysis study — no sample N Cobalt redox-active ligand complexes in aqueous catalytic assays Overpotential and H2 evolution as a function of pendant pyrazine redox reservoirs
Common misconceptions
If a model reaches 0.77 kcal mol⁻¹ MAE, DFT has been replaced and no longer needs to be computed.
The SNAr GP still needs reliable transition-state geometries, which the authors call the bottleneck, and it does not cover acid/base catalysis or other reaction classes.
A neural net that is 98% correct on spin states will be 98% correct on any transition-metal complex.
Training was octahedral first-row Cr–Ni with a specified ligand set. CSD transfer was poorer, and cyclams and early metals were ~30 kcal mol⁻¹ outliers. The target is the DFT functional used in training, not experiment.
A 0.3 eV ORR barrier from metadynamics is a measured activation energy on a gold electrode in a real cell.
It comes from a neural-network potential trained on PBE-D3 Au(100)–water without explicit cations or applied potential; facet dependence is left for later work.
Exam-style questions
Short-answer questions that ask you to explain or compare, not recall.
The SNAr paper reports raw DFT MAE 2.93 kcal mol⁻¹ and GP MAE 0.77. What is the student supposed to conclude about 'DFT accuracy'?
On that set, uncorrected DFT barrier errors are larger than the hybrid model's errors. The improvement is a statistical correction using physical-organic and TS descriptors, not a new functional, and it still requires computing the TS.
Why is 528/538 correct ground states on a test set compatible with 10 kcal mol⁻¹ MUE on CSD complexes?
The test set is drawn from the same octahedral Cr–Ni design as training. CSD structures include chemistries outside that box; a Euclidean reliability cutoff halves the CSD MUE, which is an admission that out-of-domain complexes should not be trusted at the test-set figure.
A colleague treats the gold ORR 0.3 eV barrier as proof of the experimental mechanism in alkaline fuel cells. What is missing from the model?
Explicit electrolyte cations, applied potential, and other facets. The path is associative *OOH then two *OH in pure water on Au(100) at 350 K in a learned PBE-D3 potential.
When should DFT be cited as the result versus as an interpretation of an experiment in this library?
As the result: SNAr barrier models, TM spin-state networks, NNP metadynamics. As interpretation: Zn crystal electronics, heme-model paths next to FT-ICR rates, ligand versus metal redox after measuring HER overpotential. Mixing those citations makes a computed gap look like a measured MIC or a Faradaic efficiency.
The studies
28 studies in this library bear on Density Functional Theory (DFT), ordered by citations. The first 8 are shown.
- Single-atom Co–N–C catalyzes nitroarene azo coupling
A self-supporting Co–N–C single-atom catalyst with a CoN4C8-1-2O2 site hydrogenatively couples nitroarenes to azo compounds.
- Halogen bonds are not purely electrostatic
Relativistic DFT shows halogen and hydrogen bonds share HOMO–LUMO covalency; halogen bonds have weaker electrostatics but stronger orbital mixing.
- DFT: graphene SACs break CO2 scaling relations
Single metal atoms in graphene vacancies are predicted to reduce CO2 to methanol or methane at milder limiting potentials than metal terraces by weakening *CO more than *CHO.
- Charged amino acids absorb past 250 nm
A 67-residue protein with no aromatic side chains still absorbs from 250 to 800 nm because Lys/Glu charge-transfer transitions create Protein Charge Transfer Spectra.
- FeMn single-atom nanozyme dual-readout HER2 test
An etched FeMn–N–C nanozyme boosts peroxidase-like activity 2.03-fold and reports HER2 by electrochemistry and 808 nm photothermal readout.
- Neural nets predict TM spin states and bonds
Graph-based ANNs trained on 2690 DFT geometries of Cr–Ni octahedral complexes predict spin splitting and metal–ligand distances without 3D coordinates.
- Unsymmetrical ligands still make one Pd2L4 cage isomer
Steric and geometric ligand design plus solvent polarity steer unsymmetrical dipyridyl ligands into a single cis or trans Pd2L4 cage rather than a statistical isomer mix.
- Bismuth oxyhalide films as photoelectrodes
AACVD BiOX films show halide-tuned bandgaps; untreated BiOBr gives about 0.38 mA cm−2 photoanodic current without a sacrificial donor.
Show 20 more studiesShow fewer studies
- ML plus transition states predict SNAr barriers
A Gaussian-process model that includes SNAr transition-state features predicts experimental activation energies with 0.77 kcal mol−1 MAE, beating raw DFT.
- Ultrathin porous CoP nanosheets drive HER
Phosphidation of Co3O4 yields sub-1.1 nm porous CoP sheets with exposed {200} facets and high HER mass activity.
- Defect hydroxides in UiO-66 shuttle protons
Simulations of missing-linker UiO-66 put charge-balancing OH on under-coordinated Zr and show those sites swap between hydroxide and water by rapid proton transfer with extra-framework water.
- Soft polymers wet UiO-66 without interfacial voids
MD shows flexible PVDF and PEG fill UiO-66 surface pockets; rigid PIM-1 and PS leave microvoids and fail as 70 wt% films.
- Modelled peptides inhibit SARS-CoV-2 main protease
Simulations prefer a neutral His41/Cys145 dyad and guided four designed peptides that bind the Mpro dimer and inhibit substrate cleavage with IC50 values of 3–5 μM.
- PQ and water oxy-trifluoromethylate enynes
Sunlight-excited phenanthrenequinone turns Langlois’ CF3SO2Na into CF3· and uses water as the oxygen atom to build CF3 benzofurans, benzothiophenes and indoles.
- DASA photoswitches work in chloroform if N-methyl
First-generation donor–acceptor Stenhouse adducts with a small N-methyl substituent switch well in chloroform; kinetics split photoisomerisation from cyclisation.
- Gas-phase rates of heme-model olefin epoxidation
FT-ICR and DFT show iron(IV)–oxo porphyrin cation radicals epoxidize olefins with rates that track substrate ionization energy.
- Pyrazine reservoirs that lower cobalt HER overpotential
Moving a redox-active pyrazine on a cobalt polypyridine ligand cuts hydrogen-evolution overpotential by about 200 mV versus pyridine analogues in water.
- Carbon-coated nickel nanoparticles that make CO from CO2
N-doped carbon plus a carbon coat lets metallic Ni nanoparticles reach about 94% CO faradaic efficiency, suppressing hydrogen evolution.
- Why gold carbenes hit phenol para-C–H
DFT plus deuterium controls show (PhO)3PAu carbenes add at phenol para-C, then two waters shuttle the proton, beating O–H insertion.
- ML and high-throughput screens find organic H2 photocatalysts
Screening 572 then 96 organic molecules under identical HER conditions, plus ML, uncovered unexpected molecular photocatalysts rivaling conjugated polymers.
- AIMD maps GFP's excited-state proton shuttle
TD-DFT AIMD and wavelet spectra show GFP ESPT is concerted once low-frequency chromophore modes planarize the H-bond wire.
- Solid Rh crystals isomerize butene without dissolving
A crystalline Rh σ-alkane complex swaps NBA for alkenes in the solid state and the porous ethene polymorph isomerizes 1-butene with high TOF.
- NN potentials map gold–water ORR paths
Equivariant neural-network potentials drive nanosecond metadynamics showing associative ORR on Au(100) in which *OOH is reduced to two *OH, with a ~0.3 eV barrier.
- On–off allosteric control of a bifunctional catalyst
A Pt(II) weak-link tweezer and a sulfonate regulator switch a squaramide–amine co-catalyst fully on or off in situ.
- He tagging maps copper cluster ion shapes
Helium nanodroplet mass spectrometry of 63Cu clusters shows magic He shells that match computed binding sites, including twisted-X Cu5+.
- Gold chloride cluster anions grow as zigzag chains
LDI-MS detects Au_nCl_{n+3}– anions whose lowest isomers are zigzag chains starting from AuCl4–, growing exergonically by about 41–44 kcal mol–1.
- Wide-band TA maps an Fe(II) four-state cascade
Femtosecond TA from 370–1200 nm resolves 1PALCT → 3PALCT → 3MC → 5MC decay in an iron amido chromophore, including a ~100 fs NIR CT signature.
- A new Zn phenanthroline–maleate crystal
Slow evaporation yields [Zn(phen)(maleate)(H2O)]·H2O, a distorted square-pyramidal complex with a 3.45 eV DFT gap and H-bond-dominated packing.
Learn alongside