By coupling an enhanced sampling algorithm with an orbital-localized variant of Car-Parrinello molecular dynamics, the so-called atomic centered density matrix propagation model, we reconstruct the free energy profiles along reaction pathways using different density functional approximations (DFAs) ranging from locals to hybrids. In particular, we compare the computed free energy barrier height of proton transfer (PT) reactions to those obtained by a more traditional static approach, based on the intrinsic reaction coordinate (IRC), for two case systems, namely malonaldehyde and formic acid dimer. The obtained results show that both the IRC profiles and the potentials of mean force, derived from biased dynamic trajectories, are very sensitive to the density functional approximation applied. More precisely, we observe that, with the notable exception of M06-L, local density functionals always strongly underestimate the reaction barrier heights. More generally, we find that also the shape of the free energy profile is very sensitive to the density functional choice, thus highlighting the effect, often neglected, that the choice of DFA has also in the case of dynamics simulations.

Free Energy Profiles of Proton Transfer Reactions: Density Functional Benchmark from Biased Ab Initio Dynamics / Brémond, Éric; Savarese, Marika; Rega, Nadia; Ciofini, Ilaria; Adamo, Carlo. - In: JOURNAL OF CHEMICAL THEORY AND COMPUTATION. - ISSN 1549-9626. - 18:3(2022), pp. 1501-1511. [10.1021/acs.jctc.1c01002]

Free Energy Profiles of Proton Transfer Reactions: Density Functional Benchmark from Biased Ab Initio Dynamics

Savarese, Marika;Rega, Nadia;
2022

Abstract

By coupling an enhanced sampling algorithm with an orbital-localized variant of Car-Parrinello molecular dynamics, the so-called atomic centered density matrix propagation model, we reconstruct the free energy profiles along reaction pathways using different density functional approximations (DFAs) ranging from locals to hybrids. In particular, we compare the computed free energy barrier height of proton transfer (PT) reactions to those obtained by a more traditional static approach, based on the intrinsic reaction coordinate (IRC), for two case systems, namely malonaldehyde and formic acid dimer. The obtained results show that both the IRC profiles and the potentials of mean force, derived from biased dynamic trajectories, are very sensitive to the density functional approximation applied. More precisely, we observe that, with the notable exception of M06-L, local density functionals always strongly underestimate the reaction barrier heights. More generally, we find that also the shape of the free energy profile is very sensitive to the density functional choice, thus highlighting the effect, often neglected, that the choice of DFA has also in the case of dynamics simulations.
2022
Free Energy Profiles of Proton Transfer Reactions: Density Functional Benchmark from Biased Ab Initio Dynamics / Brémond, Éric; Savarese, Marika; Rega, Nadia; Ciofini, Ilaria; Adamo, Carlo. - In: JOURNAL OF CHEMICAL THEORY AND COMPUTATION. - ISSN 1549-9626. - 18:3(2022), pp. 1501-1511. [10.1021/acs.jctc.1c01002]
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11588/915598
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 6
  • ???jsp.display-item.citation.isi??? 6
social impact