Accurate calculations of geometries and singlet-triplet energy differences for active-site models of [NiFe] hydrogenase
Delcey, Mickael G. ; Pierloot, Kristine ; Phung, Quan M. ; Vancoillie, Steven ; Lindh, Roland and Ryde, Ulf LU
(2014)
In Physical Chemistry Chemical Physics
16(17).
p.7927-7938
- Abstract
- We have studied the geometry and singlet-triplet energy difference of two mono-nuclear Ni2+ models related to the active site in [NiFe] hydrogenase. Multiconfigurational second-order perturbation theory based on a complete active-space wavefunction with an active space of 12 electrons in 12 orbitals, CASPT2(12,12), reproduces experimental bond lengths to within 1 pm. Calculated singlet-triplet energy differences agree with those obtained from coupled-cluster calculations with single, double and (perturbatively treated) triple excitations (CCSD(T)) to within 12 kJ mol(-1). For a bimetallic model of the active site of [NiFe] hydrogenase, the CASPT2(12,12) results were compared with the results obtained with an extended active space of 22... (More)
- We have studied the geometry and singlet-triplet energy difference of two mono-nuclear Ni2+ models related to the active site in [NiFe] hydrogenase. Multiconfigurational second-order perturbation theory based on a complete active-space wavefunction with an active space of 12 electrons in 12 orbitals, CASPT2(12,12), reproduces experimental bond lengths to within 1 pm. Calculated singlet-triplet energy differences agree with those obtained from coupled-cluster calculations with single, double and (perturbatively treated) triple excitations (CCSD(T)) to within 12 kJ mol(-1). For a bimetallic model of the active site of [NiFe] hydrogenase, the CASPT2(12,12) results were compared with the results obtained with an extended active space of 22 electrons in 22 orbitals. This is so large that we need to use restricted active-space theory (RASPT2). The calculations predict that the singlet state is 48-57 kJ mol(-1) more stable than the triplet state for this model of the Ni-Sl(a) state. However, in the [NiFe] hydrogenase protein, the structure around the Ni ion is far from the square-planar structure preferred by the singlet state. This destabilises the singlet state so that it is only similar to 24 kJ mol(-1) more stable than the triplet state. Finally, we have studied how various density functional theory methods compare to the experimental, CCSD(T), CASPT2, and RASPT2 results. Semi-local functionals predict the best singlet-triplet energy differences, with BP86, TPSS, and PBE giving mean unsigned errors of 12-13 kJ mol(-1) (maximum errors of 25-31 kJ mol(-1)) compared to CCSD(T). For bond lengths, several methods give good results, e. g. TPSS, BP86, and M06, with mean unsigned errors of 2 pm for the bond lengths if relativistic effects are considered. (Less)
Please use this url to cite or link to this publication:
https://lup.lub.lu.se/record/4439478
- author
- Delcey, Mickael G.
; Pierloot, Kristine
; Phung, Quan M.
; Vancoillie, Steven
; Lindh, Roland
and Ryde, Ulf
LU
- organization
- publishing date
- 2014
- type
- Contribution to journal
- publication status
- published
- subject
- in
- Physical Chemistry Chemical Physics
- volume
- 16
- issue
- 17
- pages
- 7927 - 7938
- publisher
- Royal Society of Chemistry
- external identifiers
-
- wos:000334200400033
- scopus:84898730426
- pmid:24647807
- ISSN
- 1463-9084
- DOI
- 10.1039/c4cp00253a
- language
- English
- LU publication?
- yes
- additional info
- The information about affiliations in this record was updated in December 2015. The record was previously connected to the following departments: Theoretical Chemistry (S) (011001039)
- id
- a781bfb7-fe76-4bd5-9c0f-37f61e9967cb (old id 4439478)
- date added to LUP
- 2016-04-01 14:25:20
- date last changed
- 2023-04-06 22:21:09
@article{a781bfb7-fe76-4bd5-9c0f-37f61e9967cb,
abstract = {{We have studied the geometry and singlet-triplet energy difference of two mono-nuclear Ni2+ models related to the active site in [NiFe] hydrogenase. Multiconfigurational second-order perturbation theory based on a complete active-space wavefunction with an active space of 12 electrons in 12 orbitals, CASPT2(12,12), reproduces experimental bond lengths to within 1 pm. Calculated singlet-triplet energy differences agree with those obtained from coupled-cluster calculations with single, double and (perturbatively treated) triple excitations (CCSD(T)) to within 12 kJ mol(-1). For a bimetallic model of the active site of [NiFe] hydrogenase, the CASPT2(12,12) results were compared with the results obtained with an extended active space of 22 electrons in 22 orbitals. This is so large that we need to use restricted active-space theory (RASPT2). The calculations predict that the singlet state is 48-57 kJ mol(-1) more stable than the triplet state for this model of the Ni-Sl(a) state. However, in the [NiFe] hydrogenase protein, the structure around the Ni ion is far from the square-planar structure preferred by the singlet state. This destabilises the singlet state so that it is only similar to 24 kJ mol(-1) more stable than the triplet state. Finally, we have studied how various density functional theory methods compare to the experimental, CCSD(T), CASPT2, and RASPT2 results. Semi-local functionals predict the best singlet-triplet energy differences, with BP86, TPSS, and PBE giving mean unsigned errors of 12-13 kJ mol(-1) (maximum errors of 25-31 kJ mol(-1)) compared to CCSD(T). For bond lengths, several methods give good results, e. g. TPSS, BP86, and M06, with mean unsigned errors of 2 pm for the bond lengths if relativistic effects are considered.}},
author = {{Delcey, Mickael G. and Pierloot, Kristine and Phung, Quan M. and Vancoillie, Steven and Lindh, Roland and Ryde, Ulf}},
issn = {{1463-9084}},
language = {{eng}},
number = {{17}},
pages = {{7927--7938}},
publisher = {{Royal Society of Chemistry}},
series = {{Physical Chemistry Chemical Physics}},
title = {{Accurate calculations of geometries and singlet-triplet energy differences for active-site models of [NiFe] hydrogenase}},
url = {{https://lup.lub.lu.se/search/files/3966301/5266645.pdf}},
doi = {{10.1039/c4cp00253a}},
volume = {{16}},
year = {{2014}},
}