The goal of this paper is to examine the performance of the conventional and renormalized single-reference coupled-cluster (CC) methods in calculations of the potential energy surface of the water molecule. A comparison with the results of the internally contracted multi-reference configuration interaction calculations including the quasi-degenerate Davidson correction (MRCI(Q)) and the spectroscopically accurate potential energy surface of water resulting from the use of the energy switching (ES) approach indicates that the relatively inexpensive completely renormalized (CR) CC methods with singles (S), doubles (D), and a non-iterative treatment of triples (T) or triples and quadruples (TQ), such as CR-CCSD(T), CR-CCSD(TQ), and the recently developed rigorously size extensive extension of CR-CCSD(T), termed CR-CC(2,3), provide substantial improvements in the results of conventional CCSD(T) and CCSD(TQ) calculations at larger internuclear separations. It is shown that the CR-CC(2,3) results corrected for the effect of quadruply excited clusters through the CR-CC(2,3)+Q approach can compete with the highly accurate MRCI(Q) data. The excellent agreement between the CR-CC(2,3)+Q and MRCI(Q) results suggests ways of improving the global potential energy surface of water resulting from the use of the ES approach in the regions of intermediate bond stretches and intermediate energies connecting the region of the global minimum with the asymptotic regions.

VL - 120 IS - 1-3 ER - TY - JOUR T1 - Active-Space Symmetry-Adapted-Cluster Configuration-Interaction and Equation-of-Motion Cluster Methods for High Accuracy Calculations of Potential Energy Surfaces of Radicals JF - The Journal of Chemical Physics Y1 - 2007 A1 - Y. Ohtsuka A1 - Piotr Piecuch A1 - J. R Gour A1 - M. Ehara A1 - H. Nakatsuji KW - configuration interactions KW - coupled cluster calculations KW - electron attachment KW - excited states KW - free radicals KW - ground states KW - ionisation KW - organic compounds KW - oxygen compounds KW - potential energy surfaces AB -The electron-attached (EA) and ionized (IP) symmetry-adapted-cluster configuration-interaction (SAC-CI) methods and their equation-of-motion coupled-cluster (EOMCC) analogs provide an elegant framework for studying open-shell systems. As shown in this study, these schemes require the presence of higher-order excitations, such as the four-particle-three-hole (4p-3h) or four-hole–three-particle (4h-3p) terms, in the electron attaching or ionizing operator R in order to produce accurate ground- and excited-state potential energy surfaces of radicals along bond breaking coordinates. The full inclusion of the 4p-3h/4h-3p excitations in the EA/IP SAC-CI and EOMCC methods leads to schemes which are far too expensive for calculations involving larger radicals and realistic basis sets. In order to reduce the large costs of such schemes without sacrificing accuracy, the active-space EA/IP EOMCC methodology [ J. R. Gour et al., J. Chem. Phys. 123, 134113 (2005) ] is extended to the EA/IP SAC-CI approaches with 4p-3h/4h-3p excitations. The resulting methods, which use a physically motivated set of active orbitals to pick out the most important 3p-2h/3h-2p and 4p-3h/4h-3p excitations, represent practical computational approaches for high-accuracy calculations of potential energy surfaces of radicals. To illustrate the potential offered by the active-space EA/IP SAC-CI approaches with up to 4p-3h/4h-3p excitations, the results of benchmark calculations for the potential energy surfaces of the low-lying doublet states of CH and OH are presented and compared with other SAC-CI and EOMCC methods, and full CI results.

VL - 126 IS - 16 ER - TY - JOUR T1 - Ab Initio Coupled-Cluster Calculations for Nuclei Using Methods of Quantum Chemistry JF - European Physics Journal A Y1 - 2006 A1 - Piotr Piecuch A1 - M. Wloch A1 - D. J. Dean A1 - J. R Gour A1 - M. {Hjorth-Jensen} A1 - T. Papenbrock A1 - K. Kowalski AB -We report preliminary large scale ab initio calculations of ground and excited states of 16O using quantum chemistry inspired coupled cluster methods and realistic two-body interactions. By using the renormalized Hamiltonians obtained with a no-core G-matrix approach, we obtain the virtually converged results at the level of two-body interactions. Due to the polynomial scaling with the system size that characterizes coupled cluster methods, we can probe large model spaces with up to seven major oscillator shells, for which standard non-truncated shell-model calculations are not possible.

VL - 25 IS - 1 ER - TY - JOUR T1 - Ab initio determination of solid-state nanostructure JF - Nature Y1 - 2006 A1 - P. Juhás A1 - D. M. Cherba A1 - P. M. Duxbury A1 - W. F. Punch A1 - S. J. L. Billinge AB -Advances in materials science and molecular biology followed rapidly from the ability to characterize atomic structure using single crystals1, 2, 3, 4. Structure determination is more difficult if single crystals are not available5. Many complex inorganic materials that are of interest in nanotechnology have no periodic long-range order and so their structures cannot be solved using crystallographic methods6. Here we demonstrate that ab initio structure solution of these nanostructured materials is feasible using diffraction data in combination with distance geometry methods. Precise, sub-ångström resolution distance data are experimentally available from the atomic pair distribution function (PDF)6, 7. Current PDF analysis consists of structure refinement from reasonable initial structure guesses6, 7 and it is not clear, a priori, that sufficient information exists in the PDF to obtain a unique structural solution. Here we present and validate two algorithms for structure reconstruction from precise unassigned interatomic distances for a range of clusters. We then apply the algorithms to find a unique, ab initio, structural solution for C60 from PDF data alone. This opens the door to sub-ångström resolution structure solution of nanomaterials, even when crystallographic methods fail.

VL - 440 ER - TY - JOUR T1 - Ab initio studies of the electronic structure of defects in {PbTe} JF - Physical Review B Y1 - 2006 A1 - Salameh Ahmad A1 - S. D. Mahanti A1 - Khang Hoang A1 - M. G. Kanatzidis AB -Understanding the detailed electronic structure of deep defect states in narrow band-gap semiconductors has been a challenging problem. Recently, self-consistent ab initio calculations within density functional theory using supercell models have been successful in tackling this problem. In this paper, we carry out such calculations in {PbTe,} a well-known narrow band-gap semiconductor, for a large class of defects: cationic and anionic substitutional impurities of different valence, and cationic and anionic vacancies. For the cationic defects, we study the chemical trends in the position of defect levels by looking at series of compounds {RPb2n−1Te2n,} where R is vacancy or monovalent, divalent, or trivalent atom. Similarly, for anionic defects, we study compounds {MPb2nTe2n−1,} where M is vacancy, S, Se or I. We find that the density of states near the top of the valence band and the bottom of the conduction band get significantly modified for most of these defects. This suggests that the transport properties of {PbTe} in the presence of impurities may not always be interpreted by simple carrier doping (from bound impurity states in the gap) concepts, confirming such ideas developed from qualitative and semiquantitative arguments.

VL - 74 IS - 15 N1 - Copyright {(C)} 2010 The American Physical Society; Please report any problems to prola@aps.org ER - TY - JOUR T1 - Automated derivation and parallel computer implementation of renormalized and active-space coupled-cluster methods JF - International Journal of Quantum Chemistry Y1 - 2006 A1 - Piotr Piecuch A1 - So Hirata A1 - K. Kowalski A1 - P.-D. Fan A1 - Theresa L. Windus AB -Our recent efforts that have led to an automated derivation and computer implementation of the renormalized and active-space coupled-cluster {(CC)} methods with Tensor Contraction Engine {(TCE)} are summarized. The {TCE-generated} renormalized and active-space {CC} computer codes are parallel and applicable to closed- and open-shell references, enabling accurate calculations of potential energy surfaces along bond-breaking coordinates and excited states displaying a significant multi-reference character. The effectiveness of the new codes in describing electronic quasi-degeneracies is illustrated by the renormalized {CC} calculations of the potential energy curve of {HCl} and the active-space {CC} calculations for the low-lying excited states of the Be3 system. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006

VL - 106 IS - 1 ER - TY - JOUR T1 - Ab Initio Coupled-Cluster Study of 16O JF - Physical Review Letters Y1 - 2005 A1 - M. Włoch A1 - D. J. Dean A1 - J. R Gour A1 - M. {Hjorth-Jensen} A1 - K. Kowalski A1 - T. Papenbrock A1 - P. Piecuch AB -We report converged results for the ground and excited states and matter density of 16O using realistic two-body nucleon-nucleon interactions and coupled-cluster methods and algorithms developed in quantum chemistry. Most of the binding is obtained with the coupled-cluster singles and doubles approach. Additional binding due to three-body clusters (triples) is minimal. The coupled-cluster method with singles and doubles provides a good description of the matter density, charge radius, charge form factor, and excited states of a one-particle, one-hole nature, but it cannot describe the first-excited 0+ state. Incorporation of triples has no effect on the latter finding.

VL - 94 IS - 21 ER - TY - JOUR T1 - Ab Initio Coupled-Cluster Study of O JF - Physical Review Letters Y1 - 2005 A1 - Piotr Piecuch A1 - M. Wloch A1 - D. J. Dean A1 - J. R Gour A1 - M. {Hjorth-Jensen} A1 - K. Kowalski A1 - T. Papenbrock AB -We report converged results for the ground and excited states and matter density of 16O using realistic two-body nucleon-nucleon interactions and coupled-cluster methods and algorithms developed in quantum chemistry. Most of the binding is obtained with the coupled-cluster singles and doubles approach. Additional binding due to three-body clusters (triples) is minimal. The coupled-cluster method with singles and doubles provides a good description of the matter density, charge radius, charge form factor, and excited states of a one-particle, one-hole nature, but it cannot describe the first-excited 0+ state. Incorporation of triples has no effect on the latter finding.

VL - 94 ER - TY - JOUR T1 - Active-Space Coupled-Cluster Study of Electronic States of Be 3 JF - The Journal of Chemical Physics Y1 - 2005 A1 - Piotr Piecuch A1 - K. Kowalski A1 - S. Hirata A1 - M. Wloch A1 - T.L. Windus AB -The active-space coupled-cluster (CC) and equation-of-motion (EOM) CC methods with all single and double excitations with triple excitations defined via active orbitals (CCSDt, EOMCCSDt), as implemented with TENSOR CONTRACTION ENGINE, are applied to the challenging Be3 system, which is characterized by a large number of low-lying excited states dominated by two-electron transitions and significant high-order correlation effects in the ground electronic state. It is demonstrated that the CCSDt and EOMCCSDt methods provide an excellent description of complicated electronic quasidegeneracies present in the Be3 cluster. Different strategies for defining triple excitations within the CCSDt∕EOMCCSDt approach are discussed.

VL - 123 IS - 7 ER - TY - JOUR T1 - Active-Space Equation-of=Motion Coupled-Cluster Methods for Excited States of Radicals and Other Open-Shell Systems: EA-EOMCCSDt and IP-EOMCCSDt JF - Journal of Chemical Physics Y1 - 2005 A1 - Piotr Piecuch A1 - J. R Gour A1 - M. Wloch VL - 123 IS - 13 ER - TY - JOUR T1 - Active-space equation-of-motion coupled-cluster methods for excited states of radicals and other open-shell systems: {EA-EOMCCSDt} and {IP-EOMCCSDt} JF - The Journal of Chemical Physics Y1 - 2005 A1 - Jeffrey R. Gour A1 - Piotr Piecuch A1 - Marta Wloch KW - coupled cluster calculations KW - free radicals KW - molecular electronic states AB -The previously developed active-space coupled-cluster (CC) and equation-of-motion (EOM) CC methods are extended to radicals and other open-shell systems by combining them with the electron attached (EA) and ionized (IP) EOMCC approaches. As illustrated by the calculations for the CH and OH radicals, the resulting EA-EOMCCSDt and IP-EOMCCSDt theories are capable of providing a highly accurate description of the electronic spectra of radical systems, including excited states displaying a manifestly multideterminantal nature, with the low costs that are not much greater that those characterizing the standard EOMCC singles and doubles method.

VL - 123 IS - 13 ER -