Center-of-mass problem in truncated configuration interaction and coupled-cluster calculations

The problem of center-of-mass (CM) contaminations in ab initio nuclear structure calculations using configuration interaction (CI) and coupled-cluster (CC) approaches is analyzed. A rigorous and quantitative scheme for diagnosing the CM contamination of intrinsic observables is proposed and applied to ground-state calculations for He-4 and O-16. The CI and CC calculations for O-16 based on model spaces defined via a truncation of the single-particle basis lead to sizable CM contaminations, while the importance-truncated no-core shell model based on the $N_{\max}\hbar\Omega$ space is virtually free of CM contaminations.Comment: 6 pages, 1 figure, 2 tables; v2: minor modifications, Phys. Lett. B in prin

Breaking Bonds of Open-Shell Species with the Restricted Open-Shell Size Extensive Left Eigenstate Completely Renormalized Coupled-Cluster Method

The recently developed restricted open-shell, size extensive, left eigenstate, completely renormalized (CR), coupled-cluster (CC) singles (S), doubles (D), and noniterative triples (T) approach, termed CR-CC(2,3) and abbreviated in this paper as ROCCL, is compared with the unrestricted CCSD(T) [UCCSD(T)] and multireference second-order perturbation theory (MRMP2) methods to assess the accuracy of the calculated potential energy surfaces (PESs) of eight single bond-breaking reactions of open-shell species that consist of C, H, Si, and Cl; these types of reactions are interesting because they account for part of the gas-phase chemistry in the silicon carbide chemical vapor deposition. The full configuration interaction (FCI) and multireference configuration interaction with Davidson quadruples correction [MRCI(Q)] methods are used as benchmark methods to evaluate the accuracy of the ROCCL, UCCSD(T), and MRMP2 PESs. The ROCCL PESs are found to be in reasonable agreement with the corresponding FCI or MRCI(Q) PESs in the entire region R = 1−3Re for all of the studied bond-breaking reactions. The ROCCL PESs have smaller nonparallelity error (NPE) than the UCCSD(T) ones and are comparable to those obtained with MRMP2. Both the ROCCL and UCCSD(T) PESs have significantly smaller reaction energy errors (REE) than the MRMP2 ones. Finally, an efficient strategy is proposed to estimate the ROCCL/cc-pVTZ PESs using an additivity approximation for basis set effects and correlation corrections

Ab initio coupled-cluster and configuration interaction calculations for 16-O using V_UCOM

Using the ground-state energy of 16-O obtained with the realistic V_UCOM interaction as a test case, we present a comprehensive comparison of different configuration interaction (CI) and coupled-cluster (CC) methods, analyzing the intrinsic advantages and limitations of each of the approaches. In particular, we use the importance-truncated (IT) CI and no-core shell model (NCSM) schemes with up to 4-particle-4-hole (4p4h) excitations as well as the size extensive CC methods with a complete treatment of one- and two-body clusters (CCSD) and a non-iterative treatment of connected three-body clusters via the completely renormalized correction to the CCSD energy defining the CR-CC(2,3) approach. We discuss the impact of the center-of-mass contaminations, the choice of the single-particle basis, and size-extensivity on the resulting energies. When the IT-CI and IT-NCSM methods include the 4p4h excitations and when the CC calculations include the 1p1h, 2p2h, and 3p3h clusters, as in the CR-CC(2,3) approach, we observe an excellent agreement among the different methodologies. This shows that despite their individual limitations, the IT-CI, IT-NCSM, and CC methods can provide precise and consistent ab initio nuclear structure predictions. Furthermore, the IT-CI, IT-NCSM, and CC ground-state energy values obtained with 16-O are in good agreement with the experimental value, proving that the V_UCOM two-body interaction allows for a realistic description of binding energies for heavier nuclei and that all of the methods used in this study account for most of the relevant particle correlation effects.Comment: 20 pages, 4 figures, 1 table (v2: extended version in response to referees' comments

Effects of Anacetrapib in Patients with Atherosclerotic Vascular Disease

BACKGROUND: Patients with atherosclerotic vascular disease remain at high risk for cardiovascular events despite effective statin-based treatment of low-density lipoprotein (LDL) cholesterol levels. The inhibition of cholesteryl ester transfer protein (CETP) by anacetrapib reduces LDL cholesterol levels and increases high-density lipoprotein (HDL) cholesterol levels. However, trials of other CETP inhibitors have shown neutral or adverse effects on cardiovascular outcomes. METHODS: We conducted a randomized, double-blind, placebo-controlled trial involving 30,449 adults with atherosclerotic vascular disease who were receiving intensive atorvastatin therapy and who had a mean LDL cholesterol level of 61 mg per deciliter (1.58 mmol per liter), a mean non-HDL cholesterol level of 92 mg per deciliter (2.38 mmol per liter), and a mean HDL cholesterol level of 40 mg per deciliter (1.03 mmol per liter). The patients were assigned to receive either 100 mg of anacetrapib once daily (15,225 patients) or matching placebo (15,224 patients). The primary outcome was the first major coronary event, a composite of coronary death, myocardial infarction, or coronary revascularization. RESULTS: During the median follow-up period of 4.1 years, the primary outcome occurred in significantly fewer patients in the anacetrapib group than in the placebo group (1640 of 15,225 patients [10.8%] vs. 1803 of 15,224 patients [11.8%]; rate ratio, 0.91; 95% confidence interval, 0.85 to 0.97; P=0.004). The relative difference in risk was similar across multiple prespecified subgroups. At the trial midpoint, the mean level of HDL cholesterol was higher by 43 mg per deciliter (1.12 mmol per liter) in the anacetrapib group than in the placebo group (a relative difference of 104%), and the mean level of non-HDL cholesterol was lower by 17 mg per deciliter (0.44 mmol per liter), a relative difference of -18%. There were no significant between-group differences in the risk of death, cancer, or other serious adverse events. CONCLUSIONS: Among patients with atherosclerotic vascular disease who were receiving intensive statin therapy, the use of anacetrapib resulted in a lower incidence of major coronary events than the use of placebo. (Funded by Merck and others; Current Controlled Trials number, ISRCTN48678192 ; ClinicalTrials.gov number, NCT01252953 ; and EudraCT number, 2010-023467-18 .)

Comparison of the completely renormalized equation-of-motion coupled-cluster and Quantum Monte Carlo results for the low-lying electronic states of methylene

The left-eigenstate completely renormalized (CR) equation-of-motion (EOM) coupled-cluster (CC) method with singles, doubles, and non-iterative triples, abbreviated as CR-EOMCC(2,3) [M. Woch et al., Mol. Phys. 104, 2149 (2006); P. Piecuch et al., Int. J. Quantum Chem. 109, 3268 (2009)], and the companion ground-state CR-CC(2,3) methodology [P. Piecuch and M. Woch, J. Chem. Phys. 123, 224105 (2005); P. Piecuch et al., Chem. Phys. Lett. 418, 467 (2006)] are used to determine the total electronic and adiabatic excitation energies corresponding to the ground and lowest three excited states of methylene. The emphasis is on comparing the CR-CC(2,3)/CR-EOMCC(2,3) results obtained with the large correlation-consistent basis sets of the aug-cc-pCV xZ (x = T, Q, 5) quality and the corresponding complete basis set (CBS) limits with the recently published variational and diffusion Quantum Monte Carlo (QMC) data [P. Zimmerman et al., J. Chem. Phys. 131, 124103 (2009)]. It is demonstrated that the CBS CR-CC(2,3)/CR-EOMCC(2,3) results are in very good agreement with the best QMC, i.e. diffusion MC (DMC) data, with errors in the total and adiabatic excitation energies of all calculated states on the order of a few millihartree and less than 0.1 eV, respectively, even for the challenging, strongly multi-reference C1A1 state for which the basic EOMCC approach with singles and doubles completely fails. The agreement between the CBS CR-CC(2,3)/CR-EOMCC(2,3) and variational MC (VMC) results for the total energies is not as good as in the DMC case, but the excitation energies resulting from the CBS CR-CC(2,3)/CR-EOMCC(2,3) and VMC calculations agree very well. © 2010 Taylor & Francis

Left-eigenstate completely renormalized equation-of-motion coupled-cluster methods: Review of key concepts, extension to excited states of open-shell systems, and comparison with electron-attached and ionized approaches

The recently proposed left-eigenstate completely renormalized (CR) coupled-cluster (CC) method with singles, doubles, and noniterative triples, termed CR-CC(2, 3) [Piecuch and Wloch, J Chem Phys, 2005, 123, 224105; Piecuch et al. Chem Phys Lett, 2006, 418, 467] and the companion CR-EOMCC(2, 3) methodology, which has been previously applied to singlet excited states of closed-shell molecular systems [Wloch et al. Mol Phys, 2006, 104, 2149] and in which relatively inexpensive noniterative corrections due to triple excitations derived from the biorthogonal method of moments of CC equations (MMCC) are added to the CC singles and doubles (CCSD) or equation-of-motion (EOM) CCSD energies, have been extended to excited states of open-shell species. The resulting highly efficient computer codes for the open-shell CR-EOMCC(2, 3) approach exploiting the recursively generated intermediates and fast matrix multiplication routines have been developed and interfaced with the GAMESS package, enabling CR-EOMCC(2, 3) calculations for singlet as well as nonsinglet ground and excited states of closed- and open-shell systems using the restricted Hartree-Fock or restricted open-shell Hartree-Fock references. A number of important mathematical and algorithmic details related to formal aspects and computer implementation of the CR-EOMCC(2, 3) method have been discussed, in addition to overviewing the key concepts behind the CR-EOMCC(2, 3) and biorthogonal MMCC methodologies for ground and excited states, and the numerical results involving low-lying states of the CH, CNC, C2N, N 3, and NCO species, including states dominated by two-electron transitions, have been presented. The results of the CR-EOMCC(2, 3) calculations have been compared with other CC/EOMCC approaches, including the EOMCCSD and EOMCC singles, doubles, and triples methods, and their full and active-space valence counterparts based on the electron-attached and ionized EOMCC methodologies, and the predecessor of CR-EOMCC(2, 3) termed CR-EOMCCSD(T) [Kowalski and Piecuch, J Chem Phys, 2004, 120, 1715]

Dealing with chemical reaction pathways and electronic excitations in molecular systems via renormalized and active-space coupled-cluster methods

© 2015 AIP Publishing LLC. Coupled-cluster (CC) theory has become the de facto standard for high-accuracy molecular calculations, but the widely used CC and equation-of-motion (EOM) CC approaches, such as CCSD(T) and EOMCCSD, have difficulties with capturing stronger electron correlations that characterize multi-reference molecular problems. This presentation demonstrates that many of these difficulties can be addressed by exploiting the completely renormalized (CR) CC and EOMCC approaches, such as CR-CC(2,3), CR-EOMCCSD(T), and CR-EOMCC(2,3), and their local correlation counterparts applicable to systems with hundreds of atoms, and the active-space CC/EOMCC approaches, such as CCSDt and EOMCCSDt, and their extensions to valence systems via the electron-attached and ionized formalisms