Performance of Electronic Structure Methods for the Description of Hückel-Möbius Interconversions in Extended π-Systems

2019 ◽  
Author(s):  
Tatiana Woller ◽  
Ambar Banerjee ◽  
Nitai Sylvetsky ◽  
Xavier Deraet ◽  
Frank De Proft ◽  
...  
Keyword(s):  
Electronic Structure ◽  
Basis Set ◽  
Dft Methods ◽  
Molecular Switching ◽  
Wide Range ◽  
Electronic Structure Methods ◽  
Explicitly Correlated ◽  
Relative Errors ◽  
The Stability ◽  
Basis Set Expansion

<p>Expanded porphyrins provide a versatile route to molecular switching devices due to their ability to shift between several π-conjugation topologies encoding distinct properties. Taking into account its size and huge conformational flexibility, DFT remains the workhorse for modeling such extended macrocycles. Nevertheless, the stability of Hückel and Möbius conformers depends on a complex interplay of different factors, such as hydrogen bonding, p···p stacking, steric effects, ring strain and electron delocalization. As a consequence, the selection of an exchange-correlation functional for describing the energy profile of topological switches is very difficult. For these reasons, we have examined the performance of a variety of wavefunction methods and density functionals for describing the thermochemistry and kinetics of topology interconversions across a wide range of macrocycles. Especially for hexa- and heptaphyrins, the Möbius structures have a pronouncedly stronger degree of static correlation than the Hückel and figure-eight structures, and as a result the relative energies of singly-twisted structures are a challenging test for electronic structure methods. Comparison of limited orbital space full CI calculations with CCSD(T) calculations within the same active spaces shows that post-CCSD(T) correlation contributions to relative energies are very minor. At the same time, relative energies are weakly sensitive to further basis set expansion, as proven by the minor energy differences between MP2/cc-pVDZ and explicitly correlated MP2-F12/cc-pVDZ-F12 calculations. Hence, our CCSD(T) reference values are reasonably well-converged in both 1-particle and n-particle spaces. While conventional MP2 and MP3 yield very poor results, SCS-MP2 and particularly SOS-MP2 and SCS-MP3 agree to better than 1 kcal mol<sup>-1</sup> with the CCSD(T) relative energies. Regarding DFT methods, only M06-2X provides relative errors close to chemical accuracy with a RMSD of 1.2 kcal mol<sup>-1</sup>. While the original DSD-PBEP86 double hybrid performs fairly poorly for these extended p-systems, the errors drop down to 2 kcal mol<sup>-1</sup> for the revised revDSD-PBEP86-NL, again showing that same-spin MP2-like correlation has a detrimental impact on performance like the SOS-MP2 results. </p>

Performance of Electronic Structure Methods for the Description of Hückel-Möbius Interconversions in Extended π-Systems

2019 ◽  
Author(s):  
Tatiana Woller ◽  
Ambar Banerjee ◽  
Nitai Sylvetsky ◽  
Xavier Deraet ◽  
Frank De Proft ◽  
...  
Keyword(s):  
Electronic Structure ◽  
Basis Set ◽  
Dft Methods ◽  
Molecular Switching ◽  
Wide Range ◽  
Electronic Structure Methods ◽  
Explicitly Correlated ◽  
Relative Errors ◽  
The Stability ◽  
Basis Set Expansion

<p>Expanded porphyrins provide a versatile route to molecular switching devices due to their ability to shift between several π-conjugation topologies encoding distinct properties. Taking into account its size and huge conformational flexibility, DFT remains the workhorse for modeling such extended macrocycles. Nevertheless, the stability of Hückel and Möbius conformers depends on a complex interplay of different factors, such as hydrogen bonding, p···p stacking, steric effects, ring strain and electron delocalization. As a consequence, the selection of an exchange-correlation functional for describing the energy profile of topological switches is very difficult. For these reasons, we have examined the performance of a variety of wavefunction methods and density functionals for describing the thermochemistry and kinetics of topology interconversions across a wide range of macrocycles. Especially for hexa- and heptaphyrins, the Möbius structures have a pronouncedly stronger degree of static correlation than the Hückel and figure-eight structures, and as a result the relative energies of singly-twisted structures are a challenging test for electronic structure methods. Comparison of limited orbital space full CI calculations with CCSD(T) calculations within the same active spaces shows that post-CCSD(T) correlation contributions to relative energies are very minor. At the same time, relative energies are weakly sensitive to further basis set expansion, as proven by the minor energy differences between MP2/cc-pVDZ and explicitly correlated MP2-F12/cc-pVDZ-F12 calculations. Hence, our CCSD(T) reference values are reasonably well-converged in both 1-particle and n-particle spaces. While conventional MP2 and MP3 yield very poor results, SCS-MP2 and particularly SOS-MP2 and SCS-MP3 agree to better than 1 kcal mol<sup>-1</sup> with the CCSD(T) relative energies. Regarding DFT methods, only M06-2X provides relative errors close to chemical accuracy with a RMSD of 1.2 kcal mol<sup>-1</sup>. While the original DSD-PBEP86 double hybrid performs fairly poorly for these extended p-systems, the errors drop down to 2 kcal mol<sup>-1</sup> for the revised revDSD-PBEP86-NL, again showing that same-spin MP2-like correlation has a detrimental impact on performance like the SOS-MP2 results. </p>

scholarly journals Assessment of basis sets for F12 explicitly-correlated molecular electronic-structure methods

2009 ◽  
Vol 107 (8-12) ◽  
pp. 963-975 ◽  
Author(s):  
Florian A. Bischoff ◽  
Sandra Wolfsegger ◽  
David P. Tew ◽  
Wim Klopper
Keyword(s):  
Electronic Structure ◽  
Basis Sets ◽  
Molecular Electronic ◽  
Electronic Structure Methods ◽  
Explicitly Correlated ◽  
Molecular Electronic Structure

scholarly journals Unveiling the complex pattern of intermolecular interactions responsible for the stability of the DNA duplex

2021 ◽  
Vol 12 (38) ◽  
pp. 12785-12793
Author(s):  
Ahmet Altun ◽  
Miquel Garcia-Ratés ◽  
Frank Neese ◽  
Giovanni Bistoni
Keyword(s):  
Electronic Structure ◽  
Intermolecular Interactions ◽  
Complex Pattern ◽  
Duplex Structure ◽  
Dna Duplex ◽  
Intrinsic Stability ◽  
Human Dna ◽  
Electronic Structure Methods ◽  
The Stability ◽  
Insight Into

Advanced electronic structure methods provide a new insight into the intermolecular interactions responsible for the intrinsic stability of the duplex structure of human DNA.

Divergence of Many-Body Perturbation Theory for Noncovalent Interactions of Large Molecules

2019 ◽  
Author(s):  
Brian Nguyen ◽  
Guo P Chen ◽  
Matthew M. Agee ◽  
Asbjörn M. Burow ◽  
Matthew Tang ◽  
...  
Keyword(s):  
Perturbation Theory ◽  
System Size ◽  
State Density ◽  
Second Order ◽  
Binding Energies ◽  
Basis Set ◽  
Many Body ◽  
Many Body Perturbation Theory ◽  
Relative Errors ◽  
Ground State Density

Prompted by recent reports of large errors in noncovalent interaction (NI) energies obtained from many-body perturbation theory (MBPT), we compare the performance of second-order Møller–Plesset MBPT (MP2), spin-scaled MP2, dispersion-corrected semilocal density functional approximations (DFA), and the post-Kohn–Sham random phase approximation (RPA) for predicting binding energies of supramolecular complexes contained in the S66, L7, and S30L benchmarks. All binding energies are extrapolated to the basis set limit, corrected for basis set superposition errors, and compared to reference results of the domain-based local pair-natural orbital coupled-cluster (DLPNO-CCSD(T)) or better quality. Our results confirm that MP2 severely overestimates binding energies of large complexes, producing relative errors of over 100% for several benchmark compounds. RPA relative errors consistently range between 5-10%, significantly less than reported previously using smaller basis sets, whereas spin-scaled MP2 methods show limitations similar to MP2, albeit less pronounced, and empirically dispersion-corrected DFAs perform almost as well as RPA. Regression analysis reveals a systematic increase of relative MP2 binding energy errors with the system size at a rate of approximately 1‰ per valence electron, whereas the RPA and dispersion-corrected DFA relative errors are virtually independent of the system size. These observations are corroborated by a comparison of computed rotational constants of organic molecules to gas-phase spectroscopy data contained in the ROT34 benchmark. To analyze these results, an asymptotic adiabatic connection symmetry-adapted perturbation theory (AC-SAPT) is developed which uses monomers at full coupling whose ground-state density is constrained to the ground-state density of the complex. Using the fluctuation–dissipation theorem, we obtain a nonperturbative “screened second-order” expression for the dispersion energy in terms of monomer quantities which is exact for non-overlapping subsystems and free of induction terms; a first-order RPA-like approximation to the Hartree, exchange, and correlation kernel recovers the macroscopic Lifshitz limit. The AC-SAPT expansion of the interaction energy is obtained from Taylor expansion of the coupling strength integrand. Explicit expressions for the convergence radius of the AC-SAPT series are derived within RPA and MBPT and numerically evaluated. Whereas the AC-SAPT expansion is always convergent for nondegenerate monomers when RPA is used, it is found to spuriously diverge for second-order MBPT, except for the smallest and least polarizable monomers. The divergence of the AC-SAPT series within MBPT is numerically confirmed within RPA; prior numerical results on the convergence of the SAPT expansion for MBPT methods are revisited and support this conclusion once sufficiently high orders are included. The cause of the failure of MBPT methods for NIs of large systems is missing or incomplete “electrodynamic” screening of the Coulomb interaction due to induced particle–hole pairs between electrons in different monomers, leaving the effective interaction too strong for AC-SAPT to converge. Hence, MBPT cannot be considered reliable for quantitative predictions of NIs, even in moderately polarizable molecules with a few tens of atoms. The failure to accurately account for electrodynamic polarization makes MBPT qualitatively unsuitable for applications such as NIs of nanostructures, macromolecules, and soft materials; more robust non-perturbative approaches such as RPA or coupled cluster methods should be used instead whenever possible.<br>

Ideal–Gas Thermochemical Properties for Alkanolamine and Related Species Involved in Carbon Capture Applications

2020 ◽  
Author(s):  
Nayyereh hatefi ◽  
William Smith
Keyword(s):  
Heat Capacity ◽  
Electronic Structure ◽  
Co2 Capture ◽  
Carbon Capture ◽  
Ideal Gas ◽  
Thermochemical Properties ◽  
Temperature Range ◽  
Comparison Results ◽  
Electronic Structure Methods ◽  
Protonated Forms

<div>Ideal{gas thermochemical properties (enthalpy, entropy, Gibbs energy, and heat capacity, Cp) of 49 alkanolamines potentially suitable for CO2 capture applications and their carbamate and protonated forms were calculated using two high{order electronic structure methods, G4 and G3B3 (or G3//B3LYP). We also calculate for comparison results from the commonly used B3LYP/aug-cc-pVTZ method. This data is useful for the construction of molecular{based thermodynamic models of CO2 capture processes involving these species. The Cp data for each species over the temperature range 200 K{1500 K is presented as functions of temperature in the form of NASA seven-term polynomial expressions, permitting the set of thermochemical properties to be calculated over this temperature range. The accuracy of the G3B3 and G4 results is estimated to be 1 kcal/mol and the B3LYP/aug-cc-pVTZ results are of nferior quality..</div>

Sign in / Sign up

Export Citation Format

Share Document