Generalized spin-ratio scaled MP2 method for accurate prediction of intermolecular interactions for neutral and ionic species

Research output: Contribution to journalArticleResearchpeer-review

Abstract

The accurate calculation of intermolecular interactions is important to our understanding of properties in large molecular systems. The high computational cost of the current “gold standard” method, coupled cluster with singles and doubles and perturbative triples (CCSD(T), limits its application to small- to medium-sized systems. Second-order Møller-Plesset perturbation (MP2) theory is a cheaper alternative for larger systems, although at the expense of its decreased accuracy, especially when treating van der Waals complexes. In this study, a new modification of the spin-component scaled MP2 method was proposed for a wide range of intermolecular complexes including two well-known datasets, S22 and S66, and a large dataset of ionic liquids consisting of 174 single ion pairs, IL174. It was found that the spin ratio, ϵΔs=EINT OS/EINT SS, calculated as the ratio of the opposite-spin component to the same-spin component of the interaction correlation energy fell in the range of 0.1 and 1.6, in contrast to the range of 3-4 usually observed for the ratio of absolute correlation energy, ϵs=EOS/ESS, in individual molecules. Scaled coefficients were found to become negative when the spin ratio fell in close proximity to 1.0, and therefore, the studied intermolecular complexes were divided into two groups: (1) complexes with ϵΔs< 1 and (2) complexes with ϵΔs≥ 1. A separate set of coefficients was obtained for both groups. Exclusion of counterpoise correction during scaling was found to produce superior results due to decreased error. Among a series of Dunning’s basis sets, cc-pVTZ and cc-pVQZ were found to be the best performing ones, with a mean absolute error of 1.4 kJ mol−1 and maximum errors below 6.2 kJ mol−1. The new modification, spin-ratio scaled second-order Møller-Plesset perturbation, treats both dispersion-driven and hydrogen-bonded complexes equally well, thus validating its robustness with respect to the interaction type ranging from ionic to neutral species at minimal computational cost.

Original languageEnglish
Article number064108
Number of pages14
JournalJournal of Chemical Physics
Volume146
Issue number6
DOIs
Publication statusPublished - 14 Feb 2017

Cite this

@article{289be201fa0b4a4a962360fb199c1386,
title = "Generalized spin-ratio scaled MP2 method for accurate prediction of intermolecular interactions for neutral and ionic species",
abstract = "The accurate calculation of intermolecular interactions is important to our understanding of properties in large molecular systems. The high computational cost of the current “gold standard” method, coupled cluster with singles and doubles and perturbative triples (CCSD(T), limits its application to small- to medium-sized systems. Second-order M{\o}ller-Plesset perturbation (MP2) theory is a cheaper alternative for larger systems, although at the expense of its decreased accuracy, especially when treating van der Waals complexes. In this study, a new modification of the spin-component scaled MP2 method was proposed for a wide range of intermolecular complexes including two well-known datasets, S22 and S66, and a large dataset of ionic liquids consisting of 174 single ion pairs, IL174. It was found that the spin ratio, ϵΔs=EINT OS/EINT SS, calculated as the ratio of the opposite-spin component to the same-spin component of the interaction correlation energy fell in the range of 0.1 and 1.6, in contrast to the range of 3-4 usually observed for the ratio of absolute correlation energy, ϵs=EOS/ESS, in individual molecules. Scaled coefficients were found to become negative when the spin ratio fell in close proximity to 1.0, and therefore, the studied intermolecular complexes were divided into two groups: (1) complexes with ϵΔs< 1 and (2) complexes with ϵΔs≥ 1. A separate set of coefficients was obtained for both groups. Exclusion of counterpoise correction during scaling was found to produce superior results due to decreased error. Among a series of Dunning’s basis sets, cc-pVTZ and cc-pVQZ were found to be the best performing ones, with a mean absolute error of 1.4 kJ mol−1 and maximum errors below 6.2 kJ mol−1. The new modification, spin-ratio scaled second-order M{\o}ller-Plesset perturbation, treats both dispersion-driven and hydrogen-bonded complexes equally well, thus validating its robustness with respect to the interaction type ranging from ionic to neutral species at minimal computational cost.",
author = "Samuel Tan and {Barrera Acevedo}, Santiago and Izgorodina, {Ekaterina I}",
year = "2017",
month = "2",
day = "14",
doi = "10.1063/1.4975326",
language = "English",
volume = "146",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "AIP",
number = "6",

}

Generalized spin-ratio scaled MP2 method for accurate prediction of intermolecular interactions for neutral and ionic species. / Tan, Samuel; Barrera Acevedo, Santiago; Izgorodina, Ekaterina I.

In: Journal of Chemical Physics, Vol. 146, No. 6, 064108, 14.02.2017.

Research output: Contribution to journalArticleResearchpeer-review

TY - JOUR

T1 - Generalized spin-ratio scaled MP2 method for accurate prediction of intermolecular interactions for neutral and ionic species

AU - Tan, Samuel

AU - Barrera Acevedo, Santiago

AU - Izgorodina, Ekaterina I

PY - 2017/2/14

Y1 - 2017/2/14

N2 - The accurate calculation of intermolecular interactions is important to our understanding of properties in large molecular systems. The high computational cost of the current “gold standard” method, coupled cluster with singles and doubles and perturbative triples (CCSD(T), limits its application to small- to medium-sized systems. Second-order Møller-Plesset perturbation (MP2) theory is a cheaper alternative for larger systems, although at the expense of its decreased accuracy, especially when treating van der Waals complexes. In this study, a new modification of the spin-component scaled MP2 method was proposed for a wide range of intermolecular complexes including two well-known datasets, S22 and S66, and a large dataset of ionic liquids consisting of 174 single ion pairs, IL174. It was found that the spin ratio, ϵΔs=EINT OS/EINT SS, calculated as the ratio of the opposite-spin component to the same-spin component of the interaction correlation energy fell in the range of 0.1 and 1.6, in contrast to the range of 3-4 usually observed for the ratio of absolute correlation energy, ϵs=EOS/ESS, in individual molecules. Scaled coefficients were found to become negative when the spin ratio fell in close proximity to 1.0, and therefore, the studied intermolecular complexes were divided into two groups: (1) complexes with ϵΔs< 1 and (2) complexes with ϵΔs≥ 1. A separate set of coefficients was obtained for both groups. Exclusion of counterpoise correction during scaling was found to produce superior results due to decreased error. Among a series of Dunning’s basis sets, cc-pVTZ and cc-pVQZ were found to be the best performing ones, with a mean absolute error of 1.4 kJ mol−1 and maximum errors below 6.2 kJ mol−1. The new modification, spin-ratio scaled second-order Møller-Plesset perturbation, treats both dispersion-driven and hydrogen-bonded complexes equally well, thus validating its robustness with respect to the interaction type ranging from ionic to neutral species at minimal computational cost.

AB - The accurate calculation of intermolecular interactions is important to our understanding of properties in large molecular systems. The high computational cost of the current “gold standard” method, coupled cluster with singles and doubles and perturbative triples (CCSD(T), limits its application to small- to medium-sized systems. Second-order Møller-Plesset perturbation (MP2) theory is a cheaper alternative for larger systems, although at the expense of its decreased accuracy, especially when treating van der Waals complexes. In this study, a new modification of the spin-component scaled MP2 method was proposed for a wide range of intermolecular complexes including two well-known datasets, S22 and S66, and a large dataset of ionic liquids consisting of 174 single ion pairs, IL174. It was found that the spin ratio, ϵΔs=EINT OS/EINT SS, calculated as the ratio of the opposite-spin component to the same-spin component of the interaction correlation energy fell in the range of 0.1 and 1.6, in contrast to the range of 3-4 usually observed for the ratio of absolute correlation energy, ϵs=EOS/ESS, in individual molecules. Scaled coefficients were found to become negative when the spin ratio fell in close proximity to 1.0, and therefore, the studied intermolecular complexes were divided into two groups: (1) complexes with ϵΔs< 1 and (2) complexes with ϵΔs≥ 1. A separate set of coefficients was obtained for both groups. Exclusion of counterpoise correction during scaling was found to produce superior results due to decreased error. Among a series of Dunning’s basis sets, cc-pVTZ and cc-pVQZ were found to be the best performing ones, with a mean absolute error of 1.4 kJ mol−1 and maximum errors below 6.2 kJ mol−1. The new modification, spin-ratio scaled second-order Møller-Plesset perturbation, treats both dispersion-driven and hydrogen-bonded complexes equally well, thus validating its robustness with respect to the interaction type ranging from ionic to neutral species at minimal computational cost.

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

U2 - 10.1063/1.4975326

DO - 10.1063/1.4975326

M3 - Article

VL - 146

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 6

M1 - 064108

ER -