Simple coupled-cluster singles and doubles method with perturbative inclusion of triples and explicitly correlated geminals: The CCSD(T)R12 model

To approach the complete basis set limit of the "gold-standard" coupled-cluster singles and doubles plus perturbative triples [CCSD(T)] method, we extend the recently proposed perturbative explicitly correlated coupled-cluster singles and doubles method, CCSD(2)(R12) [E. F. Valeev, Phys. Chem. Chem. Phys. 8, 106 (2008)], to account for the effect of connected three-electron correlations. The natural choice of the zeroth-order Hamiltonian produces a perturbation expansion with rigorously separable second-order energy corrections due to the explicitly correlated geminals and conventional triple and higher excitations. The resulting CCSD(T)(R12) energy is defined as a sum of the standard CCSD(T) energy and an amplitude-dependent geminal correction. The method is technically very simple: Its implementation requires no modification of the standard CCSD(T) program and the formal cost of the geminal correction is small. We investigate the performance of the open-shell version of the CCSD(T)(R12) method as a possible replacement of the standard complete-basis-set CCSD(T) energies in the high accuracy extrapolated ab initio thermochemistry model of Stanton et al. [J. Chem. Phys. 121, 11599 (2004)]. Correlation contributions to the heat of formation computed with the new method in an aug-cc-pCVXZ basis set have mean absolute basis set errors of 2.8 and 1.0 kJmol when X is T and Q, respectively. The corresponding errors of the standard CCSD(T) method are 9.1, 4.0, and 2.1 kJmol when X=T, Q, and 5. Simple two-point basis set extrapolations of standard CCSD(T) energies perform better than the explicitly correlated method for absolute correlation energies and atomization energies, but no such advantage found when computing heats of formation. A simple Schwenke-type two-point extrapolation of the CCSD(T)(R12)aug-cc-pCVXZ energies with X=T,Q yields the most accurate heats of formation found in this work, in error on average by 0.5 kJmol and at most by 1.7 kJmol.     

Basis set convergence of the coupled-cluster correction, δ(MP2)(CCSD(T)): best practices for benchmarking non-covalent interactions and the attendant revision of the S22, NBC10, HBC6, and HSG databases

In benchmark-quality studies of non-covalent interactions, it is common to estimate interaction energies at the complete basis set (CBS) coupled-cluster through perturbative triples [CCSD(T)] level of theory by adding to CBS second-order perturbation theory (MP2) a "coupled-cluster correction," δ(MP2)(CCSD(T)), evaluated in a modest basis set. This work illustrates that commonly used basis sets such as 6-31G*(0.25) can yield large, even wrongly signed, errors for δ(MP2)(CCSD(T)) that vary significantly by binding motif. Double-ζ basis sets show more reliable results when used with explicitly correlated methods to form a δ(MP2-F12)(CCSD(T(*))-F12) correction, yielding a mean absolute deviation of 0.11 kcal mol(-1) for the S22 test set. Examining the coupled-cluster correction for basis sets up to sextuple-ζ in quality reveals that δ(MP2)(CCSD(T)) converges monotonically only beyond a turning point at triple-ζ or quadruple-ζ quality. In consequence, CBS extrapolation of δ(MP2)(CCSD(T)) corrections before the turning point, generally CBS (aug-cc-pVDZ,aug-cc-pVTZ), are found to be unreliable and often inferior to aug-cc-pVTZ alone, especially for hydrogen-bonding systems. Using the findings of this paper, we revise some recent benchmarks for non-covalent interactions, namely the S22, NBC10, HBC6, and HSG test sets. The maximum differences in the revised benchmarks are 0.080, 0.060, 0.257, and 0.102 kcal mol(-1), respectively.     

On the accuracy of explicitly correlated coupled-cluster interaction energies--have orbital results been beaten yet?

The basis set convergence of weak interaction energies for dimers of noble gases helium through krypton is studied for six variants of the explicitly correlated, frozen geminal coupled-cluster singles, doubles, and noniterative triples [CCSD(T)-F12] approach: the CCSD(T)-F12a, CCSD(T)-F12b, and CCSD(T)(F12*) methods with scaled and unscaled triples. These dimers were chosen because CCSD(T) complete-basis-set (CBS) limit benchmarks are available for them to a particularly high precision. The dependence of interaction energies on the auxiliary basis sets has been investigated and it was found that the default resolution-of-identity sets cc-pVXZ/JKFIT are far from adequate in this case. Overall, employing the explicitly correlated approach clearly speeds up the basis set convergence of CCSD(T) interaction energies, however, quite surprisingly, the improvement is not as large as the one achieved by a simple addition of bond functions to the orbital basis set. Bond functions substantially improve the CCSD(T)-F12 interaction energies as well. For small and moderate bases with bond functions, the accuracy delivered by the CCSD(T)-F12 approach cannot be matched by conventional CCSD(T). However, the latter method in the largest available bases still delivers the CBS limit to a better precision than CCSD(T)-F12 in the largest bases available for that approach. Our calculations suggest that the primary reason for the limited accuracy of the large-basis CCSD(T)-F12 treatment are the approximations made at the CCSD-F12 level and the non-explicitly correlated treatment of triples. In contrast, the explicitly correlated second-order Mo?ller-Plesset perturbation theory (MP2-F12) approach is able to pinpoint the complete-basis-set limit MP2 interaction energies of rare gas dimers to a better precision than conventional MP2. Finally, we report and analyze an unexpected failure of the CCSD(T)-F12 method to deliver the core-core and core-valence correlation corrections to interaction energies consistently and accurately.     

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. Ko?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.     

Quintuple-zeta quality coupled-cluster correlation energies with triple-zeta basis sets

The explicitly-correlated coupled-cluster method CCSD(T)(R12) is extended to include F12 geminal basis functions that decay exponentially with the interelectronic distance and reproduce the form of the average Coulomb hole more accurately than linear-r(12). Equations derived using the Ansatz 2 strong orthogonality projector are presented. The convergence of the correlation energy with orbital basis set for the new CCSD(T)(F12) method is studied and found to be rapid, 98% of the basis set limit correlation energy is typically recovered using triple-zeta orbital basis sets. The performance for reaction enthalpies is assessed via a test set of 15 reactions involving 23 molecules. The title statement is found to hold equally true for total and relative correlation energies.     

Partially linearized, fully size-extensive, and reduced multireference coupled-cluster methods. I. Formalism and mutual relationship

We describe a fully size-extensive alternative of the reduced multireference (RMR) coupled-cluster (CC) method with singles (S) and doubles (D) that generates a subset of higher-than-pair cluster amplitudes, using linearized CC equations from the full CC chain, projected onto the corresponding higher-than-doubly excited configurations. This approach is referred to as partially linearized (pl) MR CCSD method and characterized by the acronym plMR CCSD. In contrast to a similar CCSDT-1 method [Y. S. Lee et al., J. Chem. Phys. 81, 5906 (1984)] this approach also considers higher than triples (currently up to hexuples), while focusing only on a small subset of such amplitudes, referred to as the primary ones. These amplitudes are selected using similar criteria as in RMR CCSD. An extension considering secondary triples via the standard (T)-type corrections, resulting in the plMR CCSD(T) method, is also considered. The relationship of RMR and plMR CCSD and CCSD(T) approaches is discussed, and their performance and characteristics are the subject of the subsequent Part II of this paper.     

Ab initio MO study on the global minimum structure of C3H3 anion: propynl-1-yl versus allenyl anions

The geometrical structures of the C₃H⁻₃ anion are surveyed at the coupled-cluster doubles (CCD) level of theory with the aug-cc-pVDZ basis set. To clarify the CCD geometries, the stable two isomers -- propynl-l-yl 1 and allenyl 2 anions -- are further optimized at the coupled-cluster singles, doubles (triples) (CCSD(T)) level of theory both with the aug-cc-pVDZ and aug-cc-pVTZ basis sets. The final energies are calculated at the CCSD(T) and the complete active space self-consistent field (CASSCF) multi-reference internally contracted CI (MRCI) levels of theory with the aug-cc-pVTZ basis set. At the MRCI level of theory including both the corrections due to the cluster energies (MRCI+Q) and the zero-point vibrational energies, the allenyl anion 2 is about 1.3 kcal mol⁻¹ lower in energy than the propynl-l-yl anion 1. These results contrast with the previous theoretical estimates, where the propynl-l-yl anion 1 is 2-3 kcal mol⁻¹ lower in energy than the allenyl anion 2. The activation energies of the intramolecular hydrogen transfer in the 1 → 2 conversion reactions are 63.5 kcal mol⁻¹ at the MRCI+Q level of theory with the aug-cc-pVTZ basis set including the zero-point energy corrections. The adiabatic electron affinity of the planer propargyl (H₂CCCH) radical, which is the global minimum of the C₃H₃ radical, is calculated to be 0.976 eV (after correction for the zero-point energy changes) at the CCSD(T) level of theory with the aug-cc-pVTZ basis set. The present electron affinity is in fairly good agreement with the experimental one (0.893 eV) observed by Oakes and Ellison.     

Approximations to complete basis set-extrapolated, highly correlated non-covalent interaction energies

The first-principles calculation of non-covalent (particularly dispersion) interactions between molecules is a considerable challenge. In this work we studied the binding energies for ten small non-covalently bonded dimers with several combinations of correlation methods (MP2, coupled-cluster single double, coupled-cluster single double (triple) (CCSD(T))), correlation-consistent basis sets (aug-cc-pVXZ, X = D, T, Q), two-point complete basis set energy extrapolations, and counterpoise corrections. For this work, complete basis set results were estimated from averaged counterpoise and non-counterpoise-corrected CCSD(T) binding energies obtained from extrapolations with aug-cc-pVQZ and aug-cc-pVTZ basis sets. It is demonstrated that, in almost all cases, binding energies converge more rapidly to the basis set limit by averaging the counterpoise and non-counterpoise corrected values than by using either counterpoise or non-counterpoise methods alone. Examination of the effect of basis set size and electron correlation shows that the triples contribution to the CCSD(T) binding energies is fairly constant with the basis set size, with a slight underestimation with CCSD(T)∕aug-cc-pVDZ compared to the value at the (estimated) complete basis set limit, and that contributions to the binding energies obtained by MP2 generally overestimate the analogous CCSD(T) contributions. Taking these factors together, we conclude that the binding energies for non-covalently bonded systems can be accurately determined using a composite method that combines CCSD(T)∕aug-cc-pVDZ with energy corrections obtained using basis set extrapolated MP2 (utilizing aug-cc-pVQZ and aug-cc-pVTZ basis sets), if all of the components are obtained by averaging the counterpoise and non-counterpoise energies. With such an approach, binding energies for the set of ten dimers are predicted with a mean absolute deviation of 0.02 kcal/mol, a maximum absolute deviation of 0.05 kcal/mol, and a mean percent absolute deviation of only 1.7%, relative to the (estimated) complete basis set CCSD(T) results. Use of this composite approach to an additional set of eight dimers gave binding energies to within 1% of previously published high-level data. It is also shown that binding within parallel and parallel-crossed conformations of naphthalene dimer is predicted by the composite approach to be 9% greater than that previously reported in the literature. The ability of some recently developed dispersion-corrected density-functional theory methods to predict the binding energies of the set of ten small dimers was also examined.     

Theoretical reference values for the AE6 and BH6 test sets from explicitly correlated coupled-cluster theory

We propose a new computational protocol to obtain highly accurate theoretical reference data. This protocol employs the explicitly correlated coupled-cluster method with iterative single and double excitations as well as perturbative triple excitations, CCSD(T)(F12), using quadruple-z\zeta basis sets. Higher excitations are accounted for by conventional CCSDT(Q) calculations using double-z\zeta basis sets, while core/core-valence correlation effects are estimated by conventional CCSD(T) calculations using quadruple-z\zeta basis sets. Finally, scalar-relativistic effects are accounted for by conventional CCSD(T) calculations using triple-z\zeta basis sets. In the present article, this protocol is applied to the popular test sets AE6 and BH6. An error analysis shows that the new reference values obtained by our computational protocol have an uncertainty of less than 1 kcal/mol (chemical accuracy). Furthermore, concerning the atomization energies, a cancellation of the basis set incompleteness error in the CCSD(T)(F12) perturbative triples contribution with the corresponding error in the contribution from higher excitations is observed. This error cancellation is diminished by the CCSD(T*)(F12) method. Thus, we recommend the use of the CCSD(T*)(F12) method only for small- and medium-sized basis sets, while the CCSD(T)(F12) approach is preferred for high-accuracy calculations in large basis sets.     

Coupled-cluster dynamic polarizabilities including triple excitations

Dynamic polarizabilities for open- and closed-shell molecules were obtained by using coupled-cluster (CC) linear response theory with full treatment of singles, doubles, and triples (CCSDT-LR) with large basis sets utilizing the NWChem software suite. By using four approximate CC methods in conjunction with augmented cc-pVNZ basis sets, we are able to evaluate the convergence in both many-electron and one-electron spaces. For systems with primarily dynamic correlation, the results for CC3 and CCSDT are almost indistinguishable. For systems with significant static correlation, the CC3 tends to overestimate the triples contribution, while the PS(T) approximation [J. Chem. Phys. 127, 164105 (2007)] produces mixed results that are heavily dependent on the accuracies provided by noniterative approaches used to correct the equation-of-motion CCSD excitation energies. Our results for open-shell systems show that the choice of reference (restricted open-shell Hartree-Fock versus unrestricted Hartree-Fock) can have a significant impact on the accuracy of polarizabilities. A simple extrapolation based on pentuple-zeta CCSD calculations and triple-zeta CCSDT calculations reproduces experimental results with good precision in most cases.     

Extension of the renormalized coupled-cluster methods exploiting left eigenstates of the similarity-transformed hamiltonian to open-shell systems: a benchmark study

The recently formulated completely renormalized coupled-cluster method with singles, doubles, and noniterative triples, exploiting the biorthogonal form of the method of moments of coupled-cluster equations (Piecuch, P.; W?och, M. J. Chem. Phys. 2005, 123, 224105; Piecuch, P.; W?och, M.; Gour, J. R.; Kinal, A. Chem. Phys. Lett. 2006, 418, 467), termed CR-CC(2,3), is extended to open-shell systems. Test calculations for bond breaking in the OH radical and the F2+ ion and singlet-triplet gaps in the CH2, HHeH, and (HFH)- biradical systems indicate that the CR-CC(2,3) approach employing the restricted open-shell Hartree--Fock (ROHF) reference is significantly more accurate than the widely used CCSD(T) method and other noniterative triples coupled-cluster approximations without making the calculations substantially more expensive. A few molecular examples, including the activation energies of the C2H4 + H --> C2H5 forward and reverse reactions and the triplet states of the CH2 and H2Si2O2 biradicals, are used to show that the dependence of the ROHF-based CR-CC(2,3) energies on the method of canonicalization of the ROHF orbitals is, for all practical purposes, negligible.     

Accurate Interaction Energies of CO2 with the 20 Naturally Occurring Amino Acids

We have performed a series of highly accurate calculations between CO₂ and the 20 naturally occurring amino acids for the investigation of the attractive noncovalent interactions. Different nucleophilic groups present in the amino acid structures were considered (α-NH₂, COOH, side groups), and the stronger binding sites were identified. A database of accurate reference interactions energies was compiled as computed by explicitly-correlated coupled-cluster singles-and-doubles, together with perturbative triples extrapolated to the complete-basis-set limit. The CCSD(F12)(T)/CBS reference values were used for comparing a variety of popular density functionals with different basis sets. Our results show that most density functionals with the triple-zeta basis set def2-TZVPP align with the CCSD(F12)(T)/CBS reference values, but errors range from 0.1 kcal/mol up to 1.0 kcal/mol.     

Renormalized coupled-cluster methods exploiting left eigenstates of the similarity-transformed Hamiltonian

Completely renormalized (CR) coupled-cluster (CC) approaches, such as CR-CCSD(T), in which one corrects the standard CC singles and doubles (CCSD) energy for the effects of triply (T) and other higher-than-doubly excited clusters [K. Kowalski and P. Piecuch, J. Chem. Phys. 113, 18 (2000)], are reformulated in terms of the left eigenstates Phimid R:L of the similarity-transformed Hamiltonian of CC theory. The resulting CR-CCSD(T)(L) or CR-CC(2,3) and other CR-CC(L) methods are derived from the new biorthogonal form of the method of moments of CC equations (MMCC) in which, in analogy to the original MMCC theory, one focuses on the noniterative corrections to standard CC energies that recover the exact, full configuration-interaction energies. One of the advantages of the biorthogonal MMCC theory, which will be further analyzed and extended to excited states in a separate paper, is a rigorous size extensivity of the basic ground-state CR-CC(L) approximations that result from it, which was slightly violated by the original CR-CCSD(T) and CR-CCSD(TQ) approaches. This includes the CR-CCSD(T)(L) or CR-CC(2,3) method discussed in this paper, in which one corrects the CCSD energy by the relatively inexpensive noniterative correction due to triples. Test calculations for bond breaking in HF, F(2), and H(2)O indicate that the noniterative CR-CCSD(T)(L) or CR-CC(2,3) approximation is very competitive with the standard CCSD(T) theory for nondegenerate closed-shell states, while being practically as accurate as the full CC approach with singles, doubles, and triples in the bond-breaking region. Calculations of the activation enthalpy for the thermal isomerizations of cyclopropane involving the trimethylene biradical as a transition state show that the noniterative CR-CCSD(T)(L) approximation is capable of providing activation enthalpies which perfectly agree with experiment.     

Implementation of the CCSD(T)-F12 method using cusp conditions

The explicitly-correlated coupled-cluster singles and doubles with perturbative triples method (CCSD(T)-F12) is implemented using the cusp conditions. Numerical tests for a set of 16 molecules have shown agreement of correlation energies within 1 mE(h) between the cusp-condition and fully-optimized CCSD(T)-F12 methods. Benchmark calculations on 13 chemical reactions with the cusp-condition CCSD(T)-F12 method reproduce experimental enthalpies within 2 kJ mol(-1). It is also shown that regular unitary-invariant ansatz cannot exactly satisfy singlet and triplet cusp conditions in open-shell situations. We present an extended ansatz which can handle both conditions exactly.     

Aliphatic C-H/pi interactions: Methane-benzene, methane-phenol, and methane-indole complexes

Noncovalent C-H/pi interactions are prevalent in biochemistry and are important in molecular recognition. In this work, we present potential energy curves for methane-benzene, methane-phenol, and methane-indole complexes as prototypes for interactions between C-H bonds and the aromatic components of phenylalanine, tyrosine, and tryptophan. Second-order perturbation theory (MP2) is used in conjunction with the aug-cc-pVDZ and aug-cc-pVTZ basis sets to determine the counterpoise-corrected interaction energy for selected complex configurations. Using corrections for higher-order electron correlation determined with coupled-cluster theory through perturbative triples [CCSD(T)] in the aug-cc-pVDZ basis set, we estimate, through an additive approximation, results at the very accurate CCSD(T)/aug-cc-pVTZ level of theory. Symmetry-adapted perturbation theory (SAPT) is employed to determine the physically significant components of the total interaction energy for each complex.     

Hybrid coupled cluster methods: combining active space coupled cluster methods with coupled cluster singles, doubles, and perturbative triples

Based on the coupled-cluster singles, doubles, and a hybrid treatment of triples (CCSD(T)-h) method developed by us [J. Shen, E. Xu, Z. Kou, and S. Li, J. Chem. Phys. 132, 114115 (2010); and ibid. 133, 234106 (2010); and ibid. 134, 044134 (2011)], we developed and implemented a new hybrid coupled cluster (CC) method, named CCSD(T)q-h, by combining CC singles and doubles, and active triples and quadruples (CCSDtq) with CCSD(T) to deal with the electronic structures of molecules with significant multireference character. These two hybrid CC methods can be solved with non-canonical and canonical MOs. With canonical MOs, the CCSD(T)-like equations in these two methods can be solved directly without iteration so that the storage of all triple excitation amplitudes can be avoided. A practical procedure to divide canonical MOs into active and inactive subsets is proposed. Numerical calculations demonstrated that CCSD(T)-h with canonical MOs can well reproduce the corresponding results obtained with non-canonical MOs. For three atom exchange reactions, we found that CCSD(T)-h can offer a significant improvement over the popular CCSD(T) method in describing the reaction barriers. For the bond-breaking processes in F(2) and H(2)O, our calculations demonstrated that CCSD(T)q-h is a good approximation to CCSDTQ over the entire bond dissociation processes.     

Full potential energy curve for N2 by the reduced multireference coupled-cluster method

Relying on a 56-dimensional reference space and using up to the correlation-consistent, polarized, valence-quadruple-zeta (cc-pVQZ) basis sets, the reduced multireference (RMR) coupled-cluster method with singles and doubles (CCSD), as well as its perturbatively corrected version for secondary triples [RMR CCSD(T)], is employed to generate the full potential energy curves for the nitrogen molecule. The resulting potentials are then compared to the recently published accurate analytic potential based on an extensive experimental data analysis [R. J. Le Roy et al., J. Chem. Phys. 125, 164310 (2006)], and the vibrational term values of these potentials are compared over the entire well. A comparison with single-reference CCSD and CCSD(T) results, as well as with earlier obtained eight-reference RMR CC results, is also made. Excellent performance of RMR CCSD, and its systematic improvement with the increasing dimension of the reference space employed, is demonstrated. For the first 19 vibrationally excited levels, which are based on experimentally observed bands, we find an absolute average deviation of 8 cm(-1) from the computed RMR CCSD/cc-pVQZ values. The perturbative correction for triples increases this deviation to 126 cm(-1), but only to 61 cm(-1) when extrapolated to the basis set limit. Both RMR CCSD and RMR CCSD(T) potentials perform well when compared to the experiment-based analytic potential in the entire range of internuclear separations.     

Ab initio potential energy and dipole moment surfaces of (H2O)2

A full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) are reported for the water dimer, (H2O)2. The CCSD(T)-PES is a very precise fit to 19,805 ab initio energies obtained with the coupled-cluster (CCSD(T)) method, using an aug-cc-pVTZ basis. The standard counterpoise correction was applied to approximately eliminate basis set superposition errors. The fit is based on an approach that incorporates the permutational symmetry of identical atoms [Huang, X.; Braams, B.; Bowman, J. M. J. Chem.Phys. 2005, 122, 044308]. The DMS is a fit to the dipole moment obtained with M?ller-Plesset (MP2) theory, using an aug-cc-pVTZ basis. The PES has an RMS fitting error of 31 cm(-1) for energies below 20,000 cm(-1), relative to the global minimum. This surface can describe various internal floppy motions, including various monomer inversions, and isomerization pathways. Ten characteristic stationary points have been located on the surface, four of which are transition structures and the rest are higher order saddle points. Their geometrical and vibrational properties are presented and compared with available previous theoretical work. The CCSD(T)-PES and MP2-DMS dissociate correctly (and symmetrically) to two H2O monomers, with D(e) = 1665.7 cm(-1) (19.93 kJ/mol). Accurate quantum calculations of the zero-point energy of the dimer (using diffusion Monte Carlo) and the monomers (using a vibrational configuration interaction approach) are reported, and these together with D(e) give a value of D0 of 1042 cm(-1) (12.44 kJ/mol). A best estimated value is 1130 cm(-1) (13.5 kJ/mol).     

Partially linearized, fully size-extensive, and reduced multireference coupled-cluster methods. II. Applications and performance

The partially linearized (pl), fully size-extensive multireference (MR) coupled-cluster (CC) method, fully accounting for singles (S) and doubles (D) and approximately for a subset of primary higher than doubles, referred to as plMR CCSD, as well as its plMR CCSD(T) version corrected for secondary triples, as described in Part I of this paper [X. Li and J. Paldus, J. Chem. Phys. 128, 144118 (2008)], are applied to the problem of bond breaking in the HF, F2, H2O, and N2 molecules, as well as to the H4 model, using basis sets of a DZ or a cc-pVDZ quality that enable a comparison with the full configuration interaction (FCI) exact energies for a given ab initio model. A comparison of the performance of the plMR CCSD/CCSD(T) approaches with those of the reduced MR (RMR) CCSD/CCSD(T) methods, as well as with the standard single reference (SR) CCSD and CCSD(T) methods, is made in each case. For the H4 model and N2 we also compare our results with the completely renormalized (CR) CC(2,3) method [P. Piecuch and M. W?och, J. Chem. Phys. 123, 224105 (2005)]. An important role of a proper choice of the model space for the MR-type methods is also addressed. The advantages and shortcomings of all these methods are pointed out and discussed, as well as their size-extensivity characteristics, in which case we distinguish supersystems involving noninteracting SR and MR subsystems from those involving only MR-type subsystems. Although the plMR-type approaches render fully size-extensive results, while the RMR CCSD may slightly violate this property, the latter method yields invariably superior results to the plMR CCSD ones and is more easy to apply in highly demanding cases, such as the triple-bond breaking in the nitrogen molecule.