Triples Correction Research Articles

Fock-space multi-reference coupled-cluster response with the effect of triples on dipole moment of ClO and SF radicals#

Dipole moment calculations of SF and ClO radicals have been carried out using the recently developed partial triples correction to Fock-space multi-reference coupled cluster method. Theoretical calculation of the doublet SF and ClO radicals is useful due to their importance in atmospheric chemistry. The dipole moments of these radicals are extremely sensitive to correlation effects. A brief insight to the way the triples correction has been implemented is presented. We compare the results obtained from our analytic response treatment with that of restricted open Hartree-Fock (ROHF) calculations. Results are presented for both relaxed and non-relaxed approach in the ROHF method. Results suggest the importance of triples corrections. The effects of orbital relaxation are also analysed from the results.

Resummation methods

AbstractResummation methods can significantly improve the accuracy of ab initio electronic structure computations without increasing the computational cost. For perturbation theories, resummation methods can be designed by constructing approximants to model the known singularity structure of the theory in the complex plane of the perturbation parameter. Quadratic approximants for the fourth‐order Møller–Plesset perturbation theory (MP4) greatly improve the accuracy for the ground‐state energy and provide information about singularity positions that can be used to select an optimal summation method. The Coupled cluster theories CCSD (coupled clusters with single and double excitations), CCSDT (with triple excitations), CCSDTQ (with quadruple excitations), and CCSD(T) (with a triples correction from perturbation theory) can be resummed using approximants that model the empirically observed convergence patterns of the Hartree–Fock (HF), CCSD, CCSD(T) and HF, CCSD, CCSDT, CCSDTQ sequences. Coupling‐constant perturbation theories of molecular vibration and of atoms in external fields, and semiclassical perturbation theories also benefit from appropriate approximants. © 2011 John Wiley & Sons, Ltd.This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods

Initiating molecular growth in the interstellar medium via dimeric complexes of observed ions and molecules

A feasible initiation step for particle growth in the interstellar medium (ISM) is simulated by means of ab quantum chemistry methods. The systems studied are dimer ions formed by pairing nitrogen containing small molecules known to exist in the ISM with ions of unsaturated hydrocarbons or vice versa. Complexation energies, structures of ensuing complexes and electronic excitation spectra of the encounter complexes are estimated using various quantum chemistry methods. Moller-Plesset perturbation theory (MP2, Z-averaged perturbation theory (ZAP2), coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)), and density functional theory (DFT) methods (B3LYP) were employed along with the correlation consistent cc-pVTZ and aug-cc-pVTZ basis sets. Two types of complexes are predicted. One type of complex has electrostatic binding with moderate (7-20 kcal per mol) binding energies, that are nonetheless significantly stronger than typical van der Waals interactions between molecules of this size. The other type of complex develops strong covalent bonds between the fragments. Cyclic isomers of the nitrogen containing complexes are produced very easily by ion-molecule reactions. Some of these complexes show intense ultraviolet visible spectra for electronic transitions with large oscillator strengths at the B3LYP, omegaB97, and equations of motion coupled cluster (EOM-CCSD) levels. The open shell nitrogen containing carbonaceous complexes especially exhibit a large oscillator strength electronic transition in the visible region of the electromagnetic spectrum.

Ab initiocomputation of the energies of circular quantum dots

We perform coupled-cluster and diffusion Monte Carlo calculations of the energies of circular quantum dots up to 20 electrons. The coupled-cluster calculations include triples corrections and a renormalized Coulomb interaction defined for a given number of low-lying oscillator shells. Using such a renormalized Coulomb interaction brings the coupled-cluster calculations with triples correlations in excellent agreement with the diffusion Monte Carlo calculations. This opens up perspectives for doing ab initio calculations for much larger systems of electrons.

Coupled cluster and DFT calculations of 14N nuclear quadrupole coupling constants

AbstractThe dependence of 14N quadrupole coupling constants calculated using coupled cluster theory on the level of approximation is examined for a series of small molecules. For HCN, HNC, CH3CN, and CH3NC, we use the coupled cluster singles‐and‐doubles with a noniterative perturbative triples correction—CCSD(T)—approach, and we analyze the basis set dependence of the results. For aziridine, diazirine, and cyclopropyl cyanide, we use the CCSD(T) approach, but smaller basis sets, and for the largest studied molecules—quinuclidine and hexamine—we present CCSD results. The differences between computed and experimental values for the best basis sets used are ≈ 5% at the CCSD level and decrease noticeably at the CCSD(T) level. The ‐ N≡C bonds are an exception—in this case the quadrupole coupling constants are very small, hence the differences between theory and experiment become larger (up to 9%). We also consider the performance of density functional theory, comparing the results for different density functionals with the coupled cluster values of the same constants. Most of the functionals provide results systematically improved with respect to the Hartree–Fock values, with 14N coupling constants in ‐ N≡C bonds being again an exception. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011

Structures and properties of neutral gallium clusters: A theoretical investigation

A systematic and unbiased structure search based on a genetic algorithm in combination with density functional theory (DFT) procedures has been carried out to locate low-energy isomers of Ga(n) up to n = 25. For the smaller clusters up to n = 8 results are checked by coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)) employing a quadruple zeta type basis set. The CCSD(T) calculations confirm a (3)Π(u) ground state for the dimer. Ga(3) has a doublet ground state 0.2 eV below two quartet states, whereas two isoenergetic triplet states are predicted for Ga(4) with D(4h) and a rhombus structure (D(2h)). Three low-lying isomers with doublet electronic states are found for Ga(5): a W-structure (C(2v)), a planar envelope (C(s)) at 0.015 eV, and a non-planar envelope (C(1)) 0.086 eV above the ground state. A triplet state for a trigonal prism (D(3h)) and a singlet for an open prism (C(2v)) are computed with virtually identical energy for Ga(6). The global minimum for Ga(7) is a capped trigonal prism (C(s)) and that for Ga(8) a distorted cube in D(2h). DFT provides a fair agreement with CCSD(T), deviations in dissociation energies are up to 0.2 eV for n ≤ 8. The structures for Ga(n) are mostly irregular for n ≥ 9, those for Ga(12) to Ga(17) can be derived from the truncated decahedron with D(5h) symmetry though highly distorted by Jahn-Teller effects, for example. For Ga(18) to Ga(23) we find stacks of five- and six-membered rings as global minima, e.g., 5-1-5-1-6 for Ga(18). Ga(24) and Ga(25) consist of layers with packing sequence ABCBA similar to those found for clusters of aluminum. The most important feature of computed cohesive energies is a rapid increase with n: for Ga(25) it reaches 2.46 eV, the experimental bulk value is 2.84 eV. Particularly stable clusters for Ga(n) are seen for n = 7, 14, and 20.

SU-E-T-184: Usefulness of Triple Channel Correction for Gafchromic EBT2 Film on Patient Specific Quality Assurance of IMRT and IMAT

Purpose: Triple channel correction acquisition (TCCA) is a model‐based method for obtaining the characteristic curve of Gafchromic EBT2 film (EBT2). In this study the TCCA method was applied to EBT2 film dosimetry of intensity modulated radiation therapy(IMRT) and intensity modulated arc therapy (IMAT). Method and Materials:EBT2 films were calibrated by a field‐by‐field method in 13 steps. All films were scanned by using flat‐bed scanners at one day after the irradiation. Seven characteristic curves obtained by the TCCA method were compared with those of a single channel acquisition (SCA) method. For clinical applications, the planar dose distributions of IMRT and IMAT plans were compared with dose distributions measured with EBT2 films. The film analyses were performed by applying both TCCA and SCA. We chose two IMRT plans and two Varian IMAT (RapidArc) plans, which were created by the BrainLAB i‐plan and Eclipse Varian treatment planning systems for IMRT and RapidArc, respectively. Results: The constancy of the characteristic curves obtained by TCCA was better than that of SCA. The most significant difference between the characteristic curves of TCCA and SCA was observed in a low dose range between 25 cGy and 100 cGy. The results of clinical cases demonstrated that the dose difference between the measured and planned doses in high dose gradient regions was smaller with TCCA than that with SCA. Conclusions: We showed that characteristic curves obtained by the TCCA method have a good constancy compared with those of SCA. The potential advantage of TCCA was seen in high dose gradient regions in clinically relevant IMRT and IMAT cases. Hence, EBT2 film dosimetry with the TCCA method can be used for high precision patient specific QA.

A universal state-selective approach to multireference coupled-cluster non-iterative corrections

A new form of the asymmetric energy functional for multireference coupled cluster (MRCC) theories is discussed from the point of view of an energy expansion in a quasidegenerate situation. The resulting expansion for the exact electronic energy can be used to define the non-iterative corrections to approximate MRCC approaches. In particular, we show that in the proposed framework the essential part of dynamic correlation can be encapsulated in the so-called correlation Hamiltonian, which in analogy to the effective Hamiltonian, is defined in the model space (M(0)). The proper parametrization of the exact/trial wavefunctions leads to the cancellation of the overlap-type numerators and to a connected form of the correlation Hamiltonian and size-extensive energies. Within this parametrization, when the trial wavefunctions are determined without invoking a specific form of the MRCC sufficiency conditions, the ensuing correction can be universally applied to any type of the approximate MRCC method. The analogies with other MRCC triples corrections to MRCC theories with singles and doubles (MRCCSD) are outlined. In particular, we discuss the approach, which in analogy to the Λ-Mk-MRCCSD(T) method [F. A. Evangelista, E. Prochnow, J. Gauss, H. F. Schaefer III, J. Chem. Phys. 132, 074107 (2010)], introduces an approximate form of the triply-excited clusters into the effective and correlation Hamiltonians. Since the discussed corrections can be calculated as a sum of independent reference-related contributions, possible parallel algorithms are also outlined.

Small clusters of aluminum and tin: Highly correlated calculations and validation of density functional procedures

We present results of molecular electronic structure treatments of multireference configuration interaction (MRCI) type for clusters Al(n) and Sn(n) in the range up to n = 4, and of coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)) type in the range up to n = 10. Basis sets of quadruple zeta size are employed, computed energy differences, such as cohesive energies, E(coh), or dissociation energies for the removal of a single atom, D(e), differ from the complete basis set limit by only a few 0.01 eV. MRCI and CCSD(T) results are then compared to those obtained from density functional theory (DFT) treatments, which show that all computational procedures agree with the general features of D(e) and E(coh). The best agreement of DFT with CCSD(T) is found for the meta-GGA (generalized gradient approximation) TPSS (Tao, Perdew, Staroverov, Scuseria) for which D(e) differs from CCSD(T) by at most 0.15 eV for Al(n) and 0.21 eV for Sn(n). The GGA PBE (Perdew, Burke, Ernzerhof) is slightly poorer with maximum deviations of 0.23 and 0.24 eV, whereas hybrid functionals are not competitive with GGA and meta-GGA functionals. A general conclusion is that errors of D(e) and/or energy differences of isomers computed with DFT procedures may easily reach 0.2 eV and errors for cohesive energies E(coh) 0.1 eV.

Determination of Shear Viscosity of Molecular Nitrogen (N2): Molecular Dynamic Hard Rotor Methodology and the Results

Determination of shear viscosity of molecular nitrogen (N(2)) by molecular dynamics (MD) in the high density range needs explicit incorporation of the rotational motion and therefore precise knowledge of angular dependence of N(2)-N(2) intermolecular potential. Newly designed Couette flow nonequilibrium molecular dynamic (NEMD) simulation procedure employs corrugated moving boundary, coupling the moving walls to translational and rotational motion exactly. Low density data on nitrogen viscosity show good agreement between MD data and experiment, confirming the radial dependence of the potential derived from quantum mechanical (QM) high precision calculations (coupled-cluster singles-and-doubles with a perturbative triples corrections [CCSD(T)]). Additionally, the angular dependence of the potential is verified using shear viscosity data for high density region, obtained from newly developed molecular dynamics (MD) simulations. It was demonstrated that the corrugated wall flow simulations provide results that are independent of the details of wall potential, fulfilling a basic requirement for application of MD simulations. These results are in good agreement with the equilibrium molecular dynamics (EMD) viscosity, derived from the Green-Kubo formula. Derived analytical dependence of the shear viscosity on the density and temperature shows that the MD data are in good agreement with experiment. Thus, MD simulations indicate that the CCSD(T) potential angular form is sufficiently precise for determination of the viscosity in a wide range of temperature and pressure.

Accurate potential energy curves for the group 12 dimers Zn2, Cd2, and Hg2

Potential energy curves of the electronic ground states of the group 12 dimers Zn2 and Cd2 were computed at the CCSD(T) level of theory, including full triple corrections \(\Updelta\)T in the coupled-cluster procedure, and spin-orbit (SO) contributions from four-component coupled-cluster calculations, extrapolated to the complete basis set (CBS) limit. For Hg2, the potential energy curve published recently (Pahl et al. in J Chem Phys 132:114301, 2010] is complemented in this work by non-relativistic calculations to quantify and discuss relativistic effects. We obtain very accurate fits of our CBS/CCSD(T) and CBS/CCSD(T)+\(\Updelta\)T data points to an analytically simple and computationally efficient extended Lennard Jones form. For the CBS/CCSD(T)+\(\Updelta\)T+SO curves, we obtain dissociation energies of De = 226 cm−1 and De = 319 cm−1 for Zn2 and Cd2 respectively, in very good agreement with recent theoretical calculations and experimental data. We also present equilibrium distances and rotational and vibrational spectroscopic constants to compare with available theoretical and experimental data. The results obtained for non-relativistically treated Hg2 continue nicely the trends with increasing atom number preset by Zn2 and Cd2, confirming that indeed, relativistic effects account for the known peculiarities for the mercury dimer.

Quantum chemical assessment of the binding energy of CuO+

We present a detailed theoretical investigation on the dissociation energy of CuO(+), carried out by means of coupled cluster theory, the multireference averaged coupled pair functional (MR-ACPF) approach, diffusion quantum Monte Carlo (DMC), and density functional theory (DFT). At the respective extrapolated basis set limits, most post-Hartree-Fock approaches agree within a narrow error margin on a D(e) value of 26.0 kcal mol(-1) [coupled-cluster singles and doubles level augmented by perturbative triples corrections, CCSD(T)], 25.8 kcal mol(-1) (CCSDTQ via the high accuracy extrapolated ab initio thermochemistry protocol), and 25.6 kcal mol(-1) (DMC), which is encouraging in view of the disaccording data published thus far. The configuration-interaction based MR-ACPF expansion, which includes single and double excitations only, gives a slightly lower value of 24.1 kcal mol(-1), indicating that large basis sets and triple excitation patterns are necessary ingredients for a quantitative assessment. Our best estimate for D(0) at the CCSD(T) level is 25.3 kcal mol(-1), which is somewhat lower than the latest experimental value (D(0) = 31.1 ± 2.8 kcal mol(-1)[semicolon] reported by the Armentrout group) [Int. J. Mass Spectrom. 182/183, 99 (1999)]. These highly correlated methods are, however, computationally very demanding, and the results are therefore supplemented with those of more affordable DFT calculations. If used in combination with moderately-sized basis sets, the M05 and M06 hybrid functionals turn out to be promising candidates for studies on much larger systems containing a [CuO](+) core.

F+ and F− Affinities of Simple NxFy and OxFy Compounds

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for the neutral and ionic N(x)F(y) and O(x)F(y) systems using coupled cluster theory with single and double excitations and including a perturbative triples correction (CCSD(T)) method with correlation consistent basis sets extrapolated to the complete basis set (CBS) limit. To achieve near chemical accuracy (±1 kcal/mol), three corrections to the electronic energy were added to the frozen core CCSD(T)/CBS binding energies: corrections for core-valence, scalar relativistic, and first order atomic spin-orbit effects. Vibrational zero point energies were computed at the CCSD(T) level of theory where possible. The calculated heats of formation are in good agreement with the available experimental values, except for FOOF because of the neglect of higher order correlation corrections. The F(+) affinity in the N(x)F(y) series increases from N(2) to N(2)F(4) by 63 kcal/mol, while that in the O(2)F(y) series decreases by 18 kcal/mol from O(2) to O(2)F(2). Neither N(2) nor N(2)F(4) is predicted to bind F(-), and N(2)F(2) is a very weak Lewis acid with an F(-) affinity of about 10 kcal/mol for either the cis or trans isomer. The low F(-) affinities of the nitrogen fluorides explain why, in spite of the fact that many stable nitrogen fluoride cations are known, no nitrogen fluoride anions have been isolated so far. For example, the F(-) affinity of NF is predicted to be only 12.5 kcal/mol which explains the numerous experimental failures to prepare NF(2)(-) salts from the well-known strong acid HNF(2). The F(-) affinity of O(2) is predicted to have a small positive value and increases for O(2)F(2) by 23 kcal/mol, indicating that the O(2)F(3)(-) anion might be marginally stable at subambient temperatures. The calculated adiabatic ionization potentials and electron affinities are in good agreement with experiment considering that many of the experimental values are for vertical processes.

Is HO3 minimum cis or trans? An analytic full-dimensional ab initio isomerization path

The minimum energy path for isomerization of HO(3) has been explored in detail using accurate high-level ab initio methods and techniques for extrapolation to the complete basis set limit. In agreement with other reports, the best estimates from both valence-only and all-electron single-reference methods here utilized predict the minimum of the cis-HO(3) isomer to be deeper than the trans-HO(3) one. They also show that the energy varies by less than 1 kcal mol(-1) or so over the full isomerization path. A similar result is found from valence-only multireference configuration interaction calculations with the size-extensive Davidson correction and a correlation consistent triple-zeta basis, which predict the energy difference between the two isomers to be of only Δ = -0.1 kcal mol(-1). However, single-point multireference calculations carried out at the optimum triple-zeta geometry with basis sets of the correlation consistent family but cardinal numbers up to X = 6 lead upon a dual-level extrapolation to the complete basis set limit of Δ = (0.12 ± 0.05) kcal mol(-1). In turn, extrapolations with the all-electron single-reference coupled-cluster method including the perturbative triples correction yield values of Δ = -0.19 and -0.03 kcal mol(-1) when done from triple-quadruple and quadruple-quintuple zeta pairs with two basis sets of increasing quality, namely cc-cpVXZ and aug-cc-pVXZ. Yet, if added a value of 0.25 kcal mol(-1) that accounts for the effect of triple and perturbative quadruple excitations with the VTZ basis set, one obtains a coupled cluster estimate of Δ = (0.14 ± 0.08) kcal mol(-1). It is then shown for the first time from systematic ab initio calculations that the trans-HO(3) isomer is more stable than the cis one, in agreement with the available experimental evidence. Inclusion of the best reported zero-point energy difference (0.382 kcal mol(-1)) from multireference configuration interaction calculations enhances further the relative stability to ΔE(ZPE) = (0.51 ± 0.08) kcal mol(-1). A scheme is also suggested to model the full-dimensional isomerization potential-energy surface using a quadratic expansion that is parametrically represented by a Fourier analysis in the torsion angle. The method illustrated at the raw and complete basis-set limit coupled-cluster levels can provide a valuable tool for a future analysis of the available (incomplete thus far) experimental rovibrational data.

The CCSD(T) model with Cholesky decomposition of orbital energy denominators

AbstractA new implementation of the coupled cluster singles and doubles with approximate triples correction method [CCSD(T)] using Cholesky decomposition of the orbital energy denominators is described. The new algorithm reduces the scaling of CCSD(T) from N7 to N6, where N is the number of orbitals. The Cholesky decomposition is carried out using simple analytical expressions that allow us to evaluate a priori the order in which the decomposition should be carried out and to obtain the relevant parts of the vectors whenever needed in the calculation. Several benchmarks have been carried out comparing the performance of the conventional and Cholesky CCSD(T) implementations. The Cholesky implementation shows a speed‐up factor larger than O2/V, where O is the number of occupied and V the number of virtual orbitals, and in general at most 5 vectors are needed to get a precision of μEh. We demonstrate that the Cholesky algorithm is better suited for studying large systems. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011

The ‘tailored’ CCSD(T) description of the automerization of cyclobutadiene

An alternative route to extend the CCSD(T) approach to multireference problems is presented. The well-known defect of the CCSD(T) model in describing the non-dynamic electron correlation effects is remedied by ‘tailoring’ the underlying coupled-cluster singles and doubles (CCSD) approach and applying the perturbative triples correction to it. The TCCSD(T) approach suggested in the paper has the same computational demands as the CCSD(T) method, though being mostly free from its drawbacks pertinent to multireference (quasidegenerate) situations. To test the approach we calculate the potential energy surface for the automerization of cyclobutadiene where the transition state exhibits a strong multireference character.

Communication: Rotational g-factor and spin-rotation constant of CH+

The rotational g-factor and spin-rotation constants of the methylidynium ion CH(+) have been calculated for the first time with a large multiconfigurational self-consistent field wave function and at the coupled-cluster singles and doubles level augmented by a perturbative triples correction. The results for an equilibrium internuclear distance as well as for the v=0, J=1 vibration-rotational state are presented.

Explicitly correlated coupled-cluster theory using cusp conditions. II. Treatment of connected triple excitations

The coupled-cluster singles and doubles method augmented with single Slater-type correlation factors (CCSD-F12) determined by the cusp conditions (also denoted as SP ansatz) yields results close to the basis set limit with only small overhead compared to conventional CCSD. Quantitative calculations on many-electron systems, however, require to include the effect of connected triple excitations at least. In this contribution, the recently proposed [A. Köhn, J. Chem. Phys. 130, 131101 (2009)] extended SP ansatz and its application to the noniterative triples correction CCSD(T) is reviewed. The approach allows to include explicit correlation into connected triple excitations without introducing additional unknown parameters. The explicit expressions are presented and analyzed, and possible simplifications to arrive at a computationally efficient scheme are suggested. Numerical tests based on an implementation obtained by an automated approach are presented. Using a partial wave expansion for the neon atom, we can show that the proposed ansatz indeed leads to the expected (L(max)+1)(-7) convergence of the noniterative triples correction, where L(max) is the maximum angular momentum in the orbital expansion. Further results are reported for a test set of 29 molecules, employing Peterson's F12-optimized basis sets. We find that the customary approach of using the conventional noniterative triples correction on top of a CCSD-F12 calculation leads to significant basis set errors. This, however, is not always directly visible for total CCSD(T) energies due to fortuitous error compensation. The new approach offers a thoroughly explicitly correlated CCSD(T)-F12 method with improved basis set convergence of the triples contributions to both total and relative energies.

Ab initiocoupled-cluster approach to nuclear structure with modern nucleon-nucleon interactions

We perform coupled-cluster calculations for the doubly magic nuclei 4He, 16O, 40Ca and 48Ca, for neutron-rich isotopes of oxygen and fluorine, and employ "bare" and secondary renormalized nucleon-nucleon interactions. For the nucleon-nucleon interaction from chiral effective field theory at order next-to-next-to-next-to leading order, we find that the coupled-cluster approximation including triples corrections binds nuclei within 0.4 MeV per nucleon compared to data. We employ interactions from a resolution-scale dependent similarity renormalization group transformations and assess the validity of power counting estimates in medium-mass nuclei. We find that the missing contributions due to three-nucleon forces are consistent with these estimates. For the unitary correlator model potential, we find a slow convergence with respect to increasing the size of the model space. For the G-matrix approach, we find a weak dependence of ground-state energies on the starting energy combined with a rather slow convergence with respect to increasing model spaces. We also analyze the center-of-mass problem and present a practical and efficient solution.

Efficient evaluation of triple excitations in symmetry-adapted perturbation theory via second-order Møller–Plesset perturbation theory natural orbitals

An accurate description of dispersion interactions is required for reliable theoretical studies of many noncovalent complexes. This can be obtained with the wave function-based formulation of symmetry-adapted perturbation theory (SAPT) provided that the contribution of triple excitations to dispersion is included. Unfortunately, this triples dispersion correction limits the applicability of SAPT due to its O(N(7)) scaling. The efficiency of the evaluation of this correction can be greatly improved by removing virtual orbitals from the computation. The error incurred from truncating the virtual space is reduced if second-order Mo̸ller-Plesset perturbation theory (MP2) natural orbitals are used in place of the canonical Hartree-Fock molecular orbitals that are typically used. This approximation is further improved if the triples correction to dispersion is scaled to account for the smaller virtual space. If virtual MP2 natural orbitals are removed according to their occupation numbers, in practice, roughly half of the virtual orbitals can be removed (with the aug-cc-pVDZ basis set) with negligible errors if the remaining triples dispersion contribution is scaled. This typically leads to speedups of 15-20 times for the cases considered here. By combining the truncated virtual space with the frozen core approximation, the triples correction can be evaluated approximately 50 times faster than the canonical computation. These approximations cause less than 1% error (or at most 0.02 kcal mol(-1)) for the cases considered. Truncation of greater fractions of the virtual space is possible for larger basis sets (leading to speedups of over 40 times before additional speedups from the frozen core approximation).