Predicting bond dissociation energy and bond length for bimetallic diatomic molecules: A challenge for electronic structure theory

Junwei Lucas Bao, Xin Zhang, Xuefei Xu, Donald G Truhlar

Research output: Contribution to journalArticle

11 Citations (Scopus)

Abstract

Accurately predicting bond length and bond dissociation energy for bimetallic diatomic molecules that involve metal-metal multiple bonds is a great challenge for electronic structure theory, in part because many of these molecules have inherently multi-configuration wave functions, a characteristic that is variously labeled as strong correlation or multireference character. Although various popular density functionals are widely used in studying metal-metal bonding in catalysis, their accuracy can be questioned, and it is important to see both how well and how poorly a functional can perform. Here we test 50 Kohn-Sham exchange-correlation density functionals for selected 3d and 4d hetero- and homonuclear bimetallic diatomic molecules against experimental bond lengths and bond energies. We found that for the majority of the density functionals, the mean unsigned error in predicting the bond length is larger than 0.08 Å, and for the bond energy, half of the functionals give a mean unsigned error larger than 20 kcal mol-1. This indicates that such highly multireference bimetallic systems are challenging for KS-DFT. However, some exchange-correlation functionals perform significantly better than average for both bond energies and bond lengths, in particular, BLYP, M06-L, N12-SX, OreLYP, RPBE, and revPBE, and are recommended for both kinds of calculations. Other functionals that perform relatively well for bond lengths include MGGA-MS0, MOHLYP, OLYP, PBE, and SOGGA11, and other functionals that perform relatively well for bond energies include GAM, M05, M06, MN15, and τ-HCTHhyb. Although some of these functionals (M05, M06, MN15, N12-SX, and τ-HCTHhyb) contain a nonzero percentage of Hartree-Fock exchange, a broader conclusion is that Hartree-Fock exchange brings in a static correlation error and usually tends to make the results, especially the bond lengths, less accurate. We find some significant differences between all-electron calculations and calculations with effective core potentials. For analysis, the article also presents CASSCF calculations of the percentage contributions of the dominant configurations, and the paper compares orbitals and configurations obtained in DFT calculations to those in CASSCF calculations. The equilibrium bond distance of Rh2 is not available from experiments, and we predict it to be 2.22 Å. The bond energy of VCr is not available from experiments, and we predict it to be 52.9 kcal mol-1.

Original languageEnglish (US)
Pages (from-to)5839-5854
Number of pages16
JournalPhysical Chemistry Chemical Physics
Volume19
Issue number8
DOIs
StatePublished - Jan 1 2017

Fingerprint

Bond length
diatomic molecules
functionals
Electronic structure
dissociation
electronic structure
Molecules
Metals
Discrete Fourier transforms
energy
configurations
metal-metal bonding
Wave functions
Catalysis
Experiments
metals
catalysis
Electrons
wave functions
orbitals

Cite this

Predicting bond dissociation energy and bond length for bimetallic diatomic molecules : A challenge for electronic structure theory. / Bao, Junwei Lucas; Zhang, Xin; Xu, Xuefei; Truhlar, Donald G.

In: Physical Chemistry Chemical Physics, Vol. 19, No. 8, 01.01.2017, p. 5839-5854.

Research output: Contribution to journalArticle

@article{b376b556ad284a5d920843991d2571bf,
title = "Predicting bond dissociation energy and bond length for bimetallic diatomic molecules: A challenge for electronic structure theory",
abstract = "Accurately predicting bond length and bond dissociation energy for bimetallic diatomic molecules that involve metal-metal multiple bonds is a great challenge for electronic structure theory, in part because many of these molecules have inherently multi-configuration wave functions, a characteristic that is variously labeled as strong correlation or multireference character. Although various popular density functionals are widely used in studying metal-metal bonding in catalysis, their accuracy can be questioned, and it is important to see both how well and how poorly a functional can perform. Here we test 50 Kohn-Sham exchange-correlation density functionals for selected 3d and 4d hetero- and homonuclear bimetallic diatomic molecules against experimental bond lengths and bond energies. We found that for the majority of the density functionals, the mean unsigned error in predicting the bond length is larger than 0.08 {\AA}, and for the bond energy, half of the functionals give a mean unsigned error larger than 20 kcal mol-1. This indicates that such highly multireference bimetallic systems are challenging for KS-DFT. However, some exchange-correlation functionals perform significantly better than average for both bond energies and bond lengths, in particular, BLYP, M06-L, N12-SX, OreLYP, RPBE, and revPBE, and are recommended for both kinds of calculations. Other functionals that perform relatively well for bond lengths include MGGA-MS0, MOHLYP, OLYP, PBE, and SOGGA11, and other functionals that perform relatively well for bond energies include GAM, M05, M06, MN15, and τ-HCTHhyb. Although some of these functionals (M05, M06, MN15, N12-SX, and τ-HCTHhyb) contain a nonzero percentage of Hartree-Fock exchange, a broader conclusion is that Hartree-Fock exchange brings in a static correlation error and usually tends to make the results, especially the bond lengths, less accurate. We find some significant differences between all-electron calculations and calculations with effective core potentials. For analysis, the article also presents CASSCF calculations of the percentage contributions of the dominant configurations, and the paper compares orbitals and configurations obtained in DFT calculations to those in CASSCF calculations. The equilibrium bond distance of Rh2 is not available from experiments, and we predict it to be 2.22 {\AA}. The bond energy of VCr is not available from experiments, and we predict it to be 52.9 kcal mol-1.",
author = "Bao, {Junwei Lucas} and Xin Zhang and Xuefei Xu and Truhlar, {Donald G}",
year = "2017",
month = "1",
day = "1",
doi = "10.1039/c6cp08896a",
language = "English (US)",
volume = "19",
pages = "5839--5854",
journal = "Physical Chemistry Chemical Physics",
issn = "1463-9076",
publisher = "Royal Society of Chemistry",
number = "8",

}

TY - JOUR

T1 - Predicting bond dissociation energy and bond length for bimetallic diatomic molecules

T2 - A challenge for electronic structure theory

AU - Bao, Junwei Lucas

AU - Zhang, Xin

AU - Xu, Xuefei

AU - Truhlar, Donald G

PY - 2017/1/1

Y1 - 2017/1/1

N2 - Accurately predicting bond length and bond dissociation energy for bimetallic diatomic molecules that involve metal-metal multiple bonds is a great challenge for electronic structure theory, in part because many of these molecules have inherently multi-configuration wave functions, a characteristic that is variously labeled as strong correlation or multireference character. Although various popular density functionals are widely used in studying metal-metal bonding in catalysis, their accuracy can be questioned, and it is important to see both how well and how poorly a functional can perform. Here we test 50 Kohn-Sham exchange-correlation density functionals for selected 3d and 4d hetero- and homonuclear bimetallic diatomic molecules against experimental bond lengths and bond energies. We found that for the majority of the density functionals, the mean unsigned error in predicting the bond length is larger than 0.08 Å, and for the bond energy, half of the functionals give a mean unsigned error larger than 20 kcal mol-1. This indicates that such highly multireference bimetallic systems are challenging for KS-DFT. However, some exchange-correlation functionals perform significantly better than average for both bond energies and bond lengths, in particular, BLYP, M06-L, N12-SX, OreLYP, RPBE, and revPBE, and are recommended for both kinds of calculations. Other functionals that perform relatively well for bond lengths include MGGA-MS0, MOHLYP, OLYP, PBE, and SOGGA11, and other functionals that perform relatively well for bond energies include GAM, M05, M06, MN15, and τ-HCTHhyb. Although some of these functionals (M05, M06, MN15, N12-SX, and τ-HCTHhyb) contain a nonzero percentage of Hartree-Fock exchange, a broader conclusion is that Hartree-Fock exchange brings in a static correlation error and usually tends to make the results, especially the bond lengths, less accurate. We find some significant differences between all-electron calculations and calculations with effective core potentials. For analysis, the article also presents CASSCF calculations of the percentage contributions of the dominant configurations, and the paper compares orbitals and configurations obtained in DFT calculations to those in CASSCF calculations. The equilibrium bond distance of Rh2 is not available from experiments, and we predict it to be 2.22 Å. The bond energy of VCr is not available from experiments, and we predict it to be 52.9 kcal mol-1.

AB - Accurately predicting bond length and bond dissociation energy for bimetallic diatomic molecules that involve metal-metal multiple bonds is a great challenge for electronic structure theory, in part because many of these molecules have inherently multi-configuration wave functions, a characteristic that is variously labeled as strong correlation or multireference character. Although various popular density functionals are widely used in studying metal-metal bonding in catalysis, their accuracy can be questioned, and it is important to see both how well and how poorly a functional can perform. Here we test 50 Kohn-Sham exchange-correlation density functionals for selected 3d and 4d hetero- and homonuclear bimetallic diatomic molecules against experimental bond lengths and bond energies. We found that for the majority of the density functionals, the mean unsigned error in predicting the bond length is larger than 0.08 Å, and for the bond energy, half of the functionals give a mean unsigned error larger than 20 kcal mol-1. This indicates that such highly multireference bimetallic systems are challenging for KS-DFT. However, some exchange-correlation functionals perform significantly better than average for both bond energies and bond lengths, in particular, BLYP, M06-L, N12-SX, OreLYP, RPBE, and revPBE, and are recommended for both kinds of calculations. Other functionals that perform relatively well for bond lengths include MGGA-MS0, MOHLYP, OLYP, PBE, and SOGGA11, and other functionals that perform relatively well for bond energies include GAM, M05, M06, MN15, and τ-HCTHhyb. Although some of these functionals (M05, M06, MN15, N12-SX, and τ-HCTHhyb) contain a nonzero percentage of Hartree-Fock exchange, a broader conclusion is that Hartree-Fock exchange brings in a static correlation error and usually tends to make the results, especially the bond lengths, less accurate. We find some significant differences between all-electron calculations and calculations with effective core potentials. For analysis, the article also presents CASSCF calculations of the percentage contributions of the dominant configurations, and the paper compares orbitals and configurations obtained in DFT calculations to those in CASSCF calculations. The equilibrium bond distance of Rh2 is not available from experiments, and we predict it to be 2.22 Å. The bond energy of VCr is not available from experiments, and we predict it to be 52.9 kcal mol-1.

UR - http://www.scopus.com/inward/record.url?scp=85027037824&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85027037824&partnerID=8YFLogxK

U2 - 10.1039/c6cp08896a

DO - 10.1039/c6cp08896a

M3 - Article

C2 - 28177019

AN - SCOPUS:85027037824

VL - 19

SP - 5839

EP - 5854

JO - Physical Chemistry Chemical Physics

JF - Physical Chemistry Chemical Physics

SN - 1463-9076

IS - 8

ER -