Reaching precise proton affinities in non-Born-Oppenheimer calculations

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Nikkanen, Luukas, Lehtola, Susi
Formato: Preprint
Publicado: 2025
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866910136839176192
author Nikkanen, Luukas
Lehtola, Susi
author_facet Nikkanen, Luukas
Lehtola, Susi
contents An attractive way to model nuclear quantum effects is to describe select nuclei quantum mechanically at the same level as the electrons. This non-Born-Oppenheimer (non-BO) method is known by many names including the nuclear-electronic orbital (NEO) and the multicomponent method. Two basis sets are typically used for such calculations: a nuclear basis set and an electronic basis set. In this work, we investigate the convergence of non-BO proton affinities (PAs) with respect to the protonic and electronic basis sets. PAs are a sensitive measure of the proton and electron densities. We demonstrate that most protonic basis sets are sufficient for non-BO density-functional calculations of PAs, resulting in convergence to within 0.1 kcal/mol of the complete protonic basis set limit. This indicates that the truncation error is dominated by the electronic basis, and that smaller protonic basis sets could be developed. We show that non-BO calculations should use uncontracted electronic basis sets on the quantum protons. The contraction coefficients in typical electronic basis sets have been derived for point nuclear charge distributions, and uncontracting the electronic basis set on the quantized proton leads to significantly faster convergence to the electronic basis set limit. Uncontraction leads to results of one $ζ$-level higher quality with negligible additional computational cost in multiple diffuse basis set families: Jensen's polarization consistent aug-pc-X basis sets, Dunning's correlation-consistent aug-cc-pVXZ basis sets, as well as the Karlsruhe def2-XZPD basis sets. In specific, the aug-pc-3 electronic basis set already affords PAs converged beyond 0.1 kcal/mol when uncontracted on the quantum proton.
format Preprint
id arxiv_https___arxiv_org_abs_2503_03966
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Reaching precise proton affinities in non-Born-Oppenheimer calculations
Nikkanen, Luukas
Lehtola, Susi
Chemical Physics
Computational Physics
An attractive way to model nuclear quantum effects is to describe select nuclei quantum mechanically at the same level as the electrons. This non-Born-Oppenheimer (non-BO) method is known by many names including the nuclear-electronic orbital (NEO) and the multicomponent method. Two basis sets are typically used for such calculations: a nuclear basis set and an electronic basis set. In this work, we investigate the convergence of non-BO proton affinities (PAs) with respect to the protonic and electronic basis sets. PAs are a sensitive measure of the proton and electron densities. We demonstrate that most protonic basis sets are sufficient for non-BO density-functional calculations of PAs, resulting in convergence to within 0.1 kcal/mol of the complete protonic basis set limit. This indicates that the truncation error is dominated by the electronic basis, and that smaller protonic basis sets could be developed. We show that non-BO calculations should use uncontracted electronic basis sets on the quantum protons. The contraction coefficients in typical electronic basis sets have been derived for point nuclear charge distributions, and uncontracting the electronic basis set on the quantized proton leads to significantly faster convergence to the electronic basis set limit. Uncontraction leads to results of one $ζ$-level higher quality with negligible additional computational cost in multiple diffuse basis set families: Jensen's polarization consistent aug-pc-X basis sets, Dunning's correlation-consistent aug-cc-pVXZ basis sets, as well as the Karlsruhe def2-XZPD basis sets. In specific, the aug-pc-3 electronic basis set already affords PAs converged beyond 0.1 kcal/mol when uncontracted on the quantum proton.
title Reaching precise proton affinities in non-Born-Oppenheimer calculations
topic Chemical Physics
Computational Physics
url https://arxiv.org/abs/2503.03966