Abstract:The accurate modeling of dissociative chemisorption of molecules on metal surfaces presents an exciting scientific challenge to theorists, and is practically relevant to modeling heterogeneously catalyzed reactive processes in computational catalysis. The first important scientific challenge in the field is that accurate barriers for dissociative chemisorption are not yet available from first principles methods. For systems that are not prone to charge transfer (for which the difference between the work function of the surface and the electron affinity of the molecule is less than 7 eV) this problem can be circumvented: chemically accurate barrier heights can be extracted with a semi-empirical version of density functional theory (DFT). However, a second important challenge is posed by systems that are prone to (full or partial) electron transfer from the surface to the molecule. For these systems the Born-Oppenheimer approximation breaks down, and currently no method of established accuracy exists for modeling the resulting effect of non-adiabatic energy dissipation on the dissociative chemisorption reaction. Because two problems exist for this class of reactions, a semi-empirical approach to computing barrier heights, which would demand that computed and experimental dissociative chemisorption probabilities match, is unlikely to work. This Perspective presents a vision on how these two problems may be solved. We suggest an approach in which parameterized density functionals are used as in the previous semi-empirical approach to DFT, but in which the parameters are based on calculations with first principles electronic structure methods. We also suggest that the diffusion Monte-Carlo (DMC) and the random phase approximation (RPA) probably are the best two first principles electronic structure methods to pursue in the framework of the approach that we call first-principles based DFT (FPB-DFT) – providing DMC and the RPA with a steppingstone towards benchmarking and future applications in computational catalysis. Probably the FPB density functional is best based on screened hybrid exchange in combination with non-local van der Waals correlation.There is a second paragraph we cannot at this stage include, but we feel it should be there, Please see the uploaded paper.

I. Dissociation free energies of drug-receptor systems via non-equilibrium alchemical simulations: a theoretical framework.

Monitoring receptor-ligand interactions between surfaces by thermal fluctuations.

Perturbation Potentials to Overcome Order/Disorder Transitions in Alchemical Binding Free Energy Calculations

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

Understanding Aldose Reductase-Inhibitors interactions with free energy simulation

Multistate Method to Efficiently Account for Tautomerism and Protonation in Alchemical Free-Energy Calculations

Binding Selectivity Analysis from Alchemical Receptor Hopping and Swapping Free Energy Calculations

Addressing Suboptimal Poses in Nonequilibrium Alchemical Calculations

On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities

A novel method for calculating relative free energy of similar molecules in two environments

Alchemical Transfer Approach to Absolute Binding Free Energy Estimation

Equilibrium and Non-equilibrium Ensemble Methods for Accurate, Precise and Reproducible Absolute Binding Free Energy Calculations

Alchemical Binding Free Energy Calculations in AMBER20: Advances and Best Practices for Drug Discovery

Alchemical Transformations and Beyond: Recent Advances and Real-World Applications of Free Energy Calculations in Drug Discovery

Validation of the Alchemical Transfer Method for the Estimation of Relative Binding Affinities of Molecular Series

Can Quantum Chemistry Improve the understanding of Protein-Ligand Interactions? Implications for Structure Based Drug Discovery

What to Make of Zero: Resolving the Statistical Noise from Conformational Reorganization in Alchemical Binding Free Energy Estimates with Metadynamics Sampling

Best-of-both-worlds computational approaches to difficult-to-model dissociation reactions on metal surfaces

Combining Alchemical Transformation with Physical Pathway to Accurately Compute Absolute Binding Free Energy

Recent theoretical and computational advances for modeling protein–ligand binding affinities