1
|
Csóka J, Kállay M. Analytic gradients for local density fitting Hartree-Fock and Kohn-Sham methods. J Chem Phys 2023; 158:024110. [PMID: 36641408 DOI: 10.1063/5.0131683] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
We present analytic gradients for local density fitting Hartree-Fock (HF) and hybrid Kohn-Sham (KS) density functional methods. Due to the non-variational nature of the local fitting algorithm, the method of Lagrange multipliers is used to avoid the solution of the coupled perturbed HF and KS equations. We propose efficient algorithms for the solution of the arising Z-vector equations and the gradient calculation that preserve the third-order scaling and low memory requirement of the original local fitting algorithm. In order to demonstrate the speed and accuracy of our implementation, gradient calculations and geometry optimizations are presented for various molecular systems. Our results show that significant speedups can be achieved compared to conventional density fitting calculations without sacrificing accuracy.
Collapse
Affiliation(s)
- József Csóka
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
2
|
Mato J, Duster AW, Guidez EB, Lin H. Adaptive-Partitioning Multilayer Dynamics Simulations: 1. On-the-Fly Switch between Two Quantum Levels of Theory. J Chem Theory Comput 2021; 17:5456-5465. [PMID: 34448578 DOI: 10.1021/acs.jctc.1c00556] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We propose to generalize the previously developed two-layer permuted adaptive-partitioning quantum-mechanics/molecular-mechanics (QM/MM), which reclassifies atoms as QM or MM on-the-fly in dynamics simulations, to multilayer adaptive-partitioning algorithms that enable multiple levels of theory. In this work, we formulate two new algorithms that smoothly interpolate the energy between two QM (Q1 and Q2) levels of theory. The first "permuted adaptive-partitioning" scheme is based on the weighted many-body expansion of the potential, as in the adaptive-partitioning QM/MM. Unconventional and potentially more efficient, the second "interpolated adaptive-partitioning" method employs alchemical QM calculations with Q1/Q2-mixed basis sets, Fock matrices, and overlap matrices. To our knowledge, this is the first time that such alchemical calculations are performed in QM, although they are routinely done in MM. Test calculations on water-cluster models show that both new algorithms indeed yield smooth energy curves when water molecules shift between Q1 and Q2.
Collapse
Affiliation(s)
- Joani Mato
- Department of Chemistry, University of Colorado, Denver, Denver, Colorado 80217, United States
| | - Adam W Duster
- Department of Chemistry, University of Colorado, Denver, Denver, Colorado 80217, United States
| | - Emilie B Guidez
- Department of Chemistry, University of Colorado, Denver, Denver, Colorado 80217, United States
| | - Hai Lin
- Department of Chemistry, University of Colorado, Denver, Denver, Colorado 80217, United States
| |
Collapse
|
3
|
Csóka J, Kállay M. Speeding up Hartree-Fock and Kohn-Sham calculations with first-order corrections. J Chem Phys 2021; 154:164114. [PMID: 33940810 DOI: 10.1063/5.0041276] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
Several approaches are presented to improve the efficiency of Hartree-Fock and Kohn-Sham self-consistent field (SCF) calculations relying on a simple first-order energy correction reminiscent of the scheme used in dual-basis SCF methods. The basic idea is to perform an initial SCF calculation computing approximate Fock-matrices and, in the final iteration step, to use a more complete Fock-matrix builder together with the energy correction to diminish the error. The approximation is tested for conventional and local density fitting (DF) SCF approaches combining various auxiliary basis sets, fitting metrics, and Fock-matrix construction algorithms in the initial and final iterations as well as for seminumerical SCF methods combining integration grids of different qualities. We also report the implementation of the occupied orbital resolution of identity exchange construction algorithm with local DF approximations. Benchmark calculations are presented for total energies, reaction energies, and molecular geometries. Our results show that speedups of up to 80% can be expected utilizing the new approaches without significant loss of accuracy.
Collapse
Affiliation(s)
- József Csóka
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
4
|
Maschio L, Kirtman B. Coupled Perturbation Theory Approach to Dual Basis Sets for Molecules and Solids. 1. General Theory and Application to Molecules. J Chem Theory Comput 2019; 16:340-353. [DOI: 10.1021/acs.jctc.9b00922] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Lorenzo Maschio
- Dipartimento di Chimica, Università di Torino, Via Giuria 5, 10125 Torino, Italy
| | - Bernard Kirtman
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
5
|
Shaw RA, Hill JG. A Linear-Scaling Method for Noncovalent Interactions: An Efficient Combination of Absolutely Localized Molecular Orbitals and a Local Random Phase Approximation Approach. J Chem Theory Comput 2019; 15:5352-5369. [PMID: 31465215 DOI: 10.1021/acs.jctc.9b00615] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A novel method for the accurate and efficient calculation of interaction energies in weakly bound complexes composed of a large number of molecules is presented. The new ALMO+RPAd method circumvents the prohibitive scaling of coupled cluster singles and doubles while still providing similar accuracy across a diverse range of weakly bound chemical systems. Linear-scaling procedures for the Fock build are given utilizing absolutely localized molecular orbitals (ALMOs), resulting in the a priori exclusion of basis set superposition errors. A bespoke data structure and algorithm using density fitting are described, leading to linear scaling for the storage and computation of the two-electron integrals. Electron correlation is included through a new, linear-scaling pairwise local random phase approximation approach, including exchange interactions, and decomposed into purely dispersive excitations (RPAxd). Collectively, these allow meaningful decomposition of the interaction energy into physically distinct contributions: electrostatic, polarization, charge transfer, and dispersion. Comparison with symmetry-adapted perturbation theory shows good qualitative agreement. Tests on various dimers and the S66 benchmark set demonstrate results within 0.5 kcal mol-1 of coupled cluster singles and doubles results. On a large cluster of water molecules, we achieve calculations involving over 3500 orbital and 12,000 auxiliary basis functions in under 10 min on a single CPU core.
Collapse
Affiliation(s)
- Robert A Shaw
- Department of Chemistry , University of Sheffield , Sheffield S3 7HF , United Kingdom
| | - J Grant Hill
- Department of Chemistry , University of Sheffield , Sheffield S3 7HF , United Kingdom
| |
Collapse
|
6
|
Hégely B, Nagy PR, Kállay M. Dual Basis Set Approach for Density Functional and Wave Function Embedding Schemes. J Chem Theory Comput 2018; 14:4600-4615. [DOI: 10.1021/acs.jctc.8b00350] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Affiliation(s)
- Bence Hégely
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Péter R. Nagy
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
7
|
Mao Y, Horn PR, Mardirossian N, Head-Gordon T, Skylaris CK, Head-Gordon M. Approaching the basis set limit for DFT calculations using an environment-adapted minimal basis with perturbation theory: Formulation, proof of concept, and a pilot implementation. J Chem Phys 2016; 145:044109. [DOI: 10.1063/1.4959125] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yuezhi Mao
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Paul R. Horn
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Teresa Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA
- Department of Bioengineering, University of California, Berkeley, California 94720, USA
- Chemical Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton, SO17 1BJ, United Kingdom
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
8
|
Cheng X, Herr JD, Steele RP. Accelerating Ab Initio Path Integral Simulations via Imaginary Multiple-Timestepping. J Chem Theory Comput 2016; 12:1627-38. [DOI: 10.1021/acs.jctc.6b00021] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaolu Cheng
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| | - Jonathan D. Herr
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| | - Ryan P. Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| |
Collapse
|
9
|
König G, Mei Y, Pickard FC, Simmonett AC, Miller BT, Herbert JM, Woodcock HL, Brooks BR, Shao Y. Computation of Hydration Free Energies Using the Multiple Environment Single System Quantum Mechanical/Molecular Mechanical Method. J Chem Theory Comput 2015; 12:332-44. [PMID: 26613419 DOI: 10.1021/acs.jctc.5b00874] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A recently developed MESS-E-QM/MM method (multiple-environment single-system quantum mechanical molecular/mechanical calculations with a Roothaan-step extrapolation) is applied to the computation of hydration free energies for the blind SAMPL4 test set and for 12 small molecules. First, free energy simulations are performed with a classical molecular mechanics force field using fixed-geometry solute molecules and explicit TIP3P solvent, and then the non-Boltzmann-Bennett method is employed to compute the QM/MM correction (QM/MM-NBB) to the molecular mechanical hydration free energies. For the SAMPL4 set, MESS-E-QM/MM-NBB corrections to the hydration free energy can be obtained 2 or 3 orders of magnitude faster than fully converged QM/MM-NBB corrections, and, on average, the hydration free energies predicted with MESS-E-QM/MM-NBB fall within 0.10-0.20 kcal/mol of full-converged QM/MM-NBB results. Out of five density functionals (BLYP, B3LYP, PBE0, M06-2X, and ωB97X-D), the BLYP functional is found to be most compatible with the TIP3P solvent model and yields the most accurate hydration free energies against experimental values for solute molecules included in this study.
Collapse
Affiliation(s)
- Gerhard König
- Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute , 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy and Department of Physics and Institute of Theoretical and Computational Science, East China Normal University , Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, China 200062
| | - Frank C Pickard
- Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute , 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute , 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Benjamin T Miller
- Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute , 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University , Columbus, Ohio 43210, United States
| | - H Lee Woodcock
- Department of Chemistry, University of South Florida , 4202 E. Fowler Avenue, Tampa, Florida 33620, United States
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute , 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Yihan Shao
- Q-Chem Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, United States
| |
Collapse
|
10
|
Steele RP. Multiple-Timestep ab Initio Molecular Dynamics Using an Atomic Basis Set Partitioning. J Phys Chem A 2015; 119:12119-30. [DOI: 10.1021/acs.jpca.5b05850] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ryan P. Steele
- Department of Chemistry and
Henry Eyring Center for Theoretical Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| |
Collapse
|
11
|
Deng J, Gilbert ATB, Gill PMW. MP2[V] – A Simple Approximation to Second-Order Møller–Plesset Perturbation Theory. J Chem Theory Comput 2015; 11:1639-44. [DOI: 10.1021/acs.jctc.5b00147] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jia Deng
- Research
School of Chemistry, The Australian National University, Canberra, Australian Capital Territory 2601, Australia
| | - Andrew T. B. Gilbert
- Research
School of Chemistry, The Australian National University, Canberra, Australian Capital Territory 2601, Australia
| | - Peter M. W. Gill
- Research
School of Chemistry, The Australian National University, Canberra, Australian Capital Territory 2601, Australia
| |
Collapse
|
12
|
Pruitt SR, Brorsen KR, Gordon MS. Ab initio investigation of the aqueous solvation of the nitrate ion. Phys Chem Chem Phys 2015; 17:27027-34. [DOI: 10.1039/c5cp04445f] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The turning point between surface and interior solvation of NO3− is predicted to lie around a cluster size of (H2O)64.
Collapse
Affiliation(s)
- Spencer R. Pruitt
- Argonne Leadership Computing Facility
- Argonne
- USA
- Department of Chemistry
- Iowa State University
| | | | | |
Collapse
|
13
|
Sodt AJ, Mei Y, König G, Tao P, Steele RP, Brooks BR, Shao Y. Multiple environment single system quantum mechanical/molecular mechanical (MESS-QM/MM) calculations. 1. Estimation of polarization energies. J Phys Chem A 2014; 119:1511-23. [PMID: 25321186 PMCID: PMC4353191 DOI: 10.1021/jp5072296] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
In
combined quantum mechanical/molecular mechanical (QM/MM) free
energy calculations, it is often advantageous to have a frozen geometry
for the quantum mechanical (QM) region. For such multiple-environment
single-system (MESS) cases, two schemes are proposed here for estimating
the polarization energy: the first scheme, termed MESS-E, involves
a Roothaan step extrapolation of the self-consistent field (SCF) energy;
whereas the other scheme, termed MESS-H, employs a Newton–Raphson
correction using an approximate inverse electronic Hessian of the
QM region (which is constructed only once). Both schemes are extremely
efficient, because the expensive Fock updates and SCF iterations in
standard QM/MM calculations are completely avoided at each configuration.
They produce reasonably accurate QM/MM polarization energies: MESS-E
can predict the polarization energy within 0.25 kcal/mol in terms
of the mean signed error for two of our test cases, solvated methanol
and solvated β-alanine, using the M06-2X or ωB97X-D functionals;
MESS-H can reproduce the polarization energy within 0.2 kcal/mol for
these two cases and for the oxyluciferin–luciferase complex,
if the approximate inverse electronic Hessians are constructed with
sufficient accuracy.
Collapse
Affiliation(s)
- Alexander J Sodt
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health , 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | | | | | | | | | | | | |
Collapse
|
14
|
Shao Y, Gan Z, Epifanovsky E, Gilbert AT, Wormit M, Kussmann J, Lange AW, Behn A, Deng J, Feng X, Ghosh D, Goldey M, Horn PR, Jacobson LD, Kaliman I, Khaliullin RZ, Kuś T, Landau A, Liu J, Proynov EI, Rhee YM, Richard RM, Rohrdanz MA, Steele RP, Sundstrom EJ, Woodcock HL, Zimmerman PM, Zuev D, Albrecht B, Alguire E, Austin B, Beran GJO, Bernard YA, Berquist E, Brandhorst K, Bravaya KB, Brown ST, Casanova D, Chang CM, Chen Y, Chien SH, Closser KD, Crittenden DL, Diedenhofen M, DiStasio RA, Do H, Dutoi AD, Edgar RG, Fatehi S, Fusti-Molnar L, Ghysels A, Golubeva-Zadorozhnaya A, Gomes J, Hanson-Heine MW, Harbach PH, Hauser AW, Hohenstein EG, Holden ZC, Jagau TC, Ji H, Kaduk B, Khistyaev K, Kim J, Kim J, King RA, Klunzinger P, Kosenkov D, Kowalczyk T, Krauter CM, Lao KU, Laurent AD, Lawler KV, Levchenko SV, Lin CY, Liu F, Livshits E, Lochan RC, Luenser A, Manohar P, Manzer SF, Mao SP, Mardirossian N, Marenich AV, Maurer SA, Mayhall NJ, Neuscamman E, Oana CM, Olivares-Amaya R, O’Neill DP, Parkhill JA, Perrine TM, Peverati R, Prociuk A, Rehn DR, Rosta E, Russ NJ, Sharada SM, Sharma S, Small DW, Sodt A, Stein T, Stück D, Su YC, Thom AJ, Tsuchimochi T, Vanovschi V, Vogt L, Vydrov O, Wang T, Watson MA, Wenzel J, White A, Williams CF, Yang J, Yeganeh S, Yost SR, You ZQ, Zhang IY, Zhang X, Zhao Y, Brooks BR, Chan GK, Chipman DM, Cramer CJ, Goddard WA, Gordon MS, Hehre WJ, Klamt A, Schaefer HF, Schmidt MW, Sherrill CD, Truhlar DG, Warshel A, Xu X, Aspuru-Guzik A, Baer R, Bell AT, Besley NA, Chai JD, Dreuw A, Dunietz BD, Furlani TR, Gwaltney SR, Hsu CP, Jung Y, Kong J, Lambrecht DS, Liang W, Ochsenfeld C, Rassolov VA, Slipchenko LV, Subotnik JE, Van Voorhis T, Herbert JM, Krylov AI, Gill PM, Head-Gordon M. Advances in molecular quantum chemistry contained in the Q-Chem 4 program package. Mol Phys 2014. [DOI: 10.1080/00268976.2014.952696] [Citation(s) in RCA: 1769] [Impact Index Per Article: 176.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
|
15
|
Use of an auxiliary basis set to describe the polarization in the fragment molecular orbital method. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.02.029] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
16
|
Byrd JN, Bartlett RJ, Montgomery JA. At what chain length do unbranched alkanes prefer folded conformations? J Phys Chem A 2014; 118:1706-12. [PMID: 24524689 DOI: 10.1021/jp4121854] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Short unbranched alkanes are known to prefer linear conformations, whereas long unbranched alkanes are folded. It is not known with certainty at what chain length the linear conformation is no longer the global minimum. To clarify this point, we use ab initio and density functional methods to compute the relative energies of the linear and hairpin alkane conformers for increasing chain lengths. Extensive electronic structure calculations are performed to obtain optimized geometries, harmonic frequencies, and accurate single point energies for the selected alkane conformers from octane through octadecane. Benchmark CCSD(T)/cc-pVTZ single point calculations are performed for chains through tetradecane, whereas approximate methods are required for the longer chains up to octadecane. Using frozen natural orbitals to unambiguously truncate the virtual orbital space, we are able to compute composite CCSD FNO(T) single point energies for all the chain lengths. This approximate composite method has significant computational savings compared to full CCSD(T) while retaining ∼0.15 kcal/mol accuracy compared to the benchmark results. More approximate dual-basis resolution-of-the-identity double-hybrid DFT calculations are also performed and shown to have reasonable 0.2-0.4 kcal/mol errors compared with our benchmark values. After including contributions from temperature dependent internal energy shifts, we find the preference for folded conformations to lie between hexadecane and octadecane, in excellent agreement with recent experiments [ Lüttschwager , N. O. ; Wassermann , T. N. ; Mata , R. A. ; Suhm , M. A. Angew. Chem. Int. Ed. 2013 , 52 , 463 ].
Collapse
Affiliation(s)
- Jason N Byrd
- Quantum Theory Project, University of Florida , Gainesville, Florida 32611, United States
| | | | | |
Collapse
|
17
|
Mardirossian N, Lambrecht DS, McCaslin L, Xantheas SS, Head-Gordon M. The Performance of Density Functionals for Sulfate–Water Clusters. J Chem Theory Comput 2013; 9:1368-80. [DOI: 10.1021/ct4000235] [Citation(s) in RCA: 65] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Narbe Mardirossian
- Department of Chemistry, University
of California, Berkeley, California, United States
| | - Daniel S. Lambrecht
- Department of Chemistry, University
of Pittsburgh, Pittsburgh, Pennsylvania, United States
| | - Laura McCaslin
- Department of Chemistry and Biochemistry,
University of Texas—Austin, Texas, United States
| | - Sotiris S. Xantheas
- Physical Sciences Division,
Pacific Northwest National Laboratory, Richland, Washington, United
States
| | - Martin Head-Gordon
- Department of Chemistry, University
of California, Berkeley, California, United States
| |
Collapse
|
18
|
Chwee TS, Lim GS, Fan WY, Sullivan MB. Computational study of molecular properties with dual basis sets. Phys Chem Chem Phys 2013; 15:16566-73. [DOI: 10.1039/c3cp51055g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
19
|
Wu A, Xu X. DCMB that combines divide-and-conquer and mixed-basis set methods for accurate geometry optimizations, total energies, and vibrational frequencies of large molecules. J Comput Chem 2012; 33:1421-32. [DOI: 10.1002/jcc.22973] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2012] [Accepted: 03/02/2012] [Indexed: 11/06/2022]
|
20
|
Deng J, Gilbert ATB, Gill PMW. Communication: Efficient counterpoise corrections by a perturbative approach. J Chem Phys 2011; 135:081105. [DOI: 10.1063/1.3632054] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
|
21
|
Mizukami W, Nakajima T, Hirao K, Yanai T. A dual-level approach to four-component relativistic density-functional theory. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.04.031] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
22
|
Kong L, Valeev EF. Perturbative correction for the basis set incompleteness error of complete-active-space self-consistent field. J Chem Phys 2010; 133:174126. [DOI: 10.1063/1.3499600] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
|
23
|
Marshall MS, Sears JS, Burns LA, Brédas JL, Sherrill CD. An Error and Efficiency Analysis of Approximations to Møller−Plesset Perturbation Theory. J Chem Theory Comput 2010. [DOI: 10.1021/ct100468f] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Michael S. Marshall
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States, and Center for Computational Molecular Science and Technology, Center for Organic Photonics and Electronics, and School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - John S. Sears
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States, and Center for Computational Molecular Science and Technology, Center for Organic Photonics and Electronics, and School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Lori A. Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States, and Center for Computational Molecular Science and Technology, Center for Organic Photonics and Electronics, and School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Jean-Luc Brédas
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States, and Center for Computational Molecular Science and Technology, Center for Organic Photonics and Electronics, and School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - C. David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States, and Center for Computational Molecular Science and Technology, Center for Organic Photonics and Electronics, and School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
24
|
|
25
|
Steele RP, Head-Gordon M, Tully JC. Ab Initio Molecular Dynamics with Dual Basis Set Methods. J Phys Chem A 2010; 114:11853-60. [DOI: 10.1021/jp107342g] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ryan P. Steele
- Department of Chemistry, Yale University, 225 Prospect St., New Haven, Connecticut 06520, United States, Department of Chemistry, University of California-Berkeley, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Department of Chemistry, Yale University, 225 Prospect St., New Haven, Connecticut 06520, United States, Department of Chemistry, University of California-Berkeley, Berkeley, California 94720, United States
| | - John C. Tully
- Department of Chemistry, Yale University, 225 Prospect St., New Haven, Connecticut 06520, United States, Department of Chemistry, University of California-Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
26
|
Steele RP, Head-Gordon M. Dual-basis self-consistent field methods: 6-31G* calculations with a minimal 6-4G primary basis. Mol Phys 2010. [DOI: 10.1080/00268970701519754] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
27
|
Distasio RA, Steele RP, Head-Gordon M. The analytical gradient of dual-basis resolution-of-the-identity second-order Møller–Plesset perturbation theory. Mol Phys 2010. [DOI: 10.1080/00268970701624687] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
28
|
Deng J, Gilbert ATB, Gill PMW. Hartree–Fock perturbative corrections for total and reaction energies. J Chem Phys 2010; 133:044116. [DOI: 10.1063/1.3463800] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
|
29
|
|
30
|
Steele RP, DiStasio Jr RA, Head-Gordon M, Li Y, Galli G. The 1,4-phenylenediisocyanide dimer: gas-phase properties and insights into organic self-assembled monolayers. Phys Chem Chem Phys 2010; 12:82-96. [DOI: 10.1039/b902194a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
31
|
Deng J, Gilbert ATB, Gill PMW. Approaching the Hartree–Fock limit by perturbative methods. J Chem Phys 2009; 130:231101. [DOI: 10.1063/1.3152864] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
|
32
|
Steele RP, DiStasio RA, Head-Gordon M. Non-Covalent Interactions with Dual-Basis Methods: Pairings for Augmented Basis Sets. J Chem Theory Comput 2009; 5:1560-72. [DOI: 10.1021/ct900058p] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ryan P. Steele
- Department of Chemistry, University of California, Berkeley, California 94720
| | - Robert A. DiStasio
- Department of Chemistry, University of California, Berkeley, California 94720
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720
| |
Collapse
|
33
|
|
34
|
Tokura S, Sato T, Tsuneda T, Nakajima T, Hirao K. A dual-level state-specific time-dependent density-functional theory. J Comput Chem 2008; 29:1187-97. [DOI: 10.1002/jcc.20871] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
35
|
Prascher BP, Wilson BR, Wilson AK. Behavior of density functionals with respect to basis set. VI. Truncation of the correlation consistent basis sets. J Chem Phys 2007; 127:124110. [PMID: 17902896 DOI: 10.1063/1.2768602] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Systematic truncation of the correlation consistent basis sets has been investigated in first and second row atoms and molecules to better understand basis set requirements for density functional theory, particularly the need for high angular momentum functions, as well as to understand possible computational cost savings that could be achieved by using reduced basis sets. The truncation scheme employed follows that recently introduced for ab initio methods [B. Mintz et al., J. Chem. Phys. 121, 5629 (2004)]. Properties examined in the current study include geometries, ionization potentials, electron affinities, and dissociation energies. In general, this investigation shows that a degree of truncation of higher angular momentum functions is possible with limited impact upon energetic properties, and does result in useful CPU time savings. However, not all properties investigated have the same level of dependence upon high angular momentum functions, and, thus, careful selection of truncated basis sets should be made.
Collapse
Affiliation(s)
- Brian P Prascher
- Center for Advanced Scientific Computing and Modeling (CASCaM), Department of Chemistry, University of North Texas, Box 305070, Denton, Texas 76203-5070, USA
| | | | | |
Collapse
|