Garnet Kin-Lic Chan

The radical character of the acenes: A density matrix renormalization group study

We present a detailed investigation of the acene series using high-level wave function theory. Our ab initio density matrix renormalization group algorithm has enabled us to carry out complete active space calculations on the acenes from napthalene to dodecacene correlating the full pi-valence space. While we find that the ground state is a singlet for all chain lengths, examination of several measures of radical character, including the natural orbitals, effective number of unpaired electrons, and various correlation functions, suggests that the longer acene ground states are polyradical in nature.

The Density Matrix Renormalization Group in Quantum Chemistry

The density matrix renormalization group is a method that is useful for describing molecules that have strongly correlated electrons. Here we provide a pedagogical overview of the basic challenges of strong correlation, how the density matrix renormalization group works, a survey of its existing applications to molecular problems, and some thoughts on the future of the method.

Orbital optimization in the density matrix renormalization group, with applications to polyenes and β-carotene

In previous work we have shown that the density matrix renormalization group (DMRG) enables near-exact calculations in active spaces much larger than are possible with traditional complete active space algorithms. Here, we implement orbital optimization with the DMRG to further allow the self-consistent improvement of the active orbitals, as is done in the complete active space self-consistent field (CASSCF) method. We use our resulting DMRG-CASSCF method to study the low-lying excited states of the all-trans polyenes up to C24H26 as well as beta-carotene, correlating with near-exact accuracy the optimized complete pi-valence space with up to 24 active electrons and orbitals, and analyze our results in the light of the recent discovery from resonance Raman experiments of new optically dark states in the spectrum.

Canonical transformation theory for multireference problems

We propose a theory to describe dynamic correlations in bonding situations where there is also significant nondynamic character. We call this the canonical transformation (CT) theory. When combined with a suitable description of nondynamic correlation, such as given by a complete-active-space self-consistent Field (CASSCF) or density matrix renormalization group wave function, it provides a theory to describe bonding situations across the entire potential energy surface with quantitative accuracy for both dynamic and nondynamic correlation. The canonical transformation theory uses a unitary exponential ansatz, is size consistent, and has a computational cost of the same order as a single-reference coupled cluster theory with the same level of excitations. Calculations using the CASSCF based CT method with single and double operators for the potential energy curves for water and nitrogen molecules, the BeH(2) insertion reaction, and hydrogen fluoride and boron hydride bond breaking, consistently yield quantitative accuracies typical of equilibrium region coupled cluster theory, but across all geometries, and better than obtained with multireference perturbation theory.

State-of-the-art density matrix renormalization group and coupled cluster theory studies of the nitrogen binding curve

We study the nitrogen binding curve with the density matrix renormalization group (DMRG) and single-reference and multireference coupled cluster (CC) theory. Our DMRG calculations use up to 4000 states and our single-reference CC calculations include up to full connected hextuple excitations. Using the DMRG, we compute an all-electron benchmark nitrogen binding curve, at the polarized, valence double-zeta level (28 basis functions), with an estimated accuracy of 0.03 mEh. We also assess the performance of more approximate DMRG and CC theories across the nitrogen curve. We provide an analysis of the relative strengths and merits of the DMRG and CC theory under different correlation conditions.

Ab initio determination of the crystalline benzene lattice energy to sub-kilojoule/mole accuracy

Working out how to pack benzene in silico Many organic compounds crystallize in several different energetically similar packing arrangements, or polymorphs. This complicates processes such as drug formulation that rely on reproducible crystallization. Yang et al. have now achieved the long-standing goal of calculating a crystal packing arrangement from first principles to an accuracy that can distinguish polymorphs (see the Perspective by Price). They used benzene as a prototypical test case and applied quantum chemical methods that improve estimates of multibody interactions. The results bode well for future applications of theory to optimization of crystallization protocols. Science, this issue p. 640; see also p. 619 Theoretical calculations of molecular packing in crystals attain sufficient accuracy to distinguish polymorphs. [Also see Perspective by Price] Computation of lattice energies to an accuracy sufficient to distinguish polymorphs is a fundamental bottleneck in crystal structure prediction. For the lattice energy of the prototypical benzene crystal, we combined the quantum chemical advances of the last decade to attain sub-kilojoule per mole accuracy, an order-of-magnitude improvement in certainty over prior calculations that necessitates revision of the experimental extrapolation to 0 kelvin. Our computations reveal the nature of binding by improving on previously inaccessible or inaccurate multibody and many-electron contributions and provide revised estimates of the effects of temperature, vibrations, and relaxation. Our demonstration raises prospects for definitive first-principles resolution of competing polymorphs in molecular crystal structure prediction.

Multireference correlation in long molecules with the quadratic scaling density matrix renormalization group

We have devised a local ab initio density matrix renormalization group algorithm to describe multireference correlations in large systems. For long molecules that are extended in one of their spatial dimensions, we can obtain an exact characterization of correlation, in the given basis, with a cost that scales only quadratically with the size of the system. The reduced scaling is achieved solely through integral screening and without the construction of correlation domains. We demonstrate the scaling, convergence, and robustness of the algorithm in polyenes and hydrogen chains. We converge to exact correlation energies (in the sense of full configuration interaction, with 1-10 microE(h) precision) in all cases and correlate up to 100 electrons in 100 active orbitals. We further use our algorithm to obtain exact energies for the metal-insulator transition in hydrogen chains and compare and contrast our results with those from conventional quantum chemical methods.

Multireference quantum chemistry through a joint density matrix renormalization group and canonical transformation theory

We describe the joint application of the density matrix renormalization group and canonical transformation theory to multireference quantum chemistry. The density matrix renormalization group provides the ability to describe static correlation in large active spaces, while the canonical transformation theory provides a high-order description of the dynamic correlation effects. We demonstrate the joint theory in two benchmark systems designed to test the dynamic and static correlation capabilities of the methods, namely, (i) total correlation energies in long polyenes and (ii) the isomerization curve of the [Cu(2)O(2)](2+) core. The largest complete active spaces and atomic orbital basis sets treated by the joint DMRG-CT theory in these systems correspond to a (24e,24o) active space and 268 atomic orbitals in the polyenes and a (28e,32o) active space and 278 atomic orbitals in [Cu(2)O(2)](2+).

Tensor factorizations of local second-order Møller-Plesset theory.

Efficient electronic structure methods can be built around efficient tensor representations of the wavefunction. Here we first describe a general view of tensor factorization for the compact representation of electronic wavefunctions. Next, we use this language to construct a low-complexity representation of the doubles amplitudes in local second-order Møller-Plesset perturbation theory. We introduce two approximations--the direct orbital-specific virtual approximation and the full orbital-specific virtual approximation. In these approximations, each occupied orbital is associated with a small set of correlating virtual orbitals. Conceptually, the representation lies between the projected atomic orbital representation in Pulay-Saebø local correlation theories and pair natural orbital correlation theories. We have tested the orbital-specific virtual approximations on a variety of systems and properties including total energies, reaction energies, and potential energy curves. Compared to the Pulay-Saebø ansatz, we find that these approximations exhibit favorable accuracy and computational times while yielding smooth potential energy curves.

The orbital-specific-virtual local coupled cluster singles and doubles method

We extend the orbital-specific-virtual tensor factorization, introduced for local Møller-Plesset perturbation theory in Ref. [J. Yang, Y. Kurashige, F. R. Manby and G. K. L. Chan, J. Chem. Phys. 134, 044123 (2011)], to local coupled cluster singles and doubles theory (OSV-LCCSD). The method is implemented by modifying an efficient projected-atomic-orbital local coupled cluster program (PAO-LCCSD) described recently, [H.-J. Werner and M. Schütz, J. Chem. Phys. 135, 144116 (2011)]. By comparison of both methods we find that the compact representation of the amplitudes in the OSV approach affords various advantages, including smaller computational time requirements (for comparable accuracy), as well as a more systematic control of the error through a single energy threshold. Overall, the OSV-LCCSD approach together with an MP2 correction yields small domain errors in practical calculations. The applicability of the OSV-LCCSD is demonstrated for molecules with up to 73 atoms and realistic basis sets (up to 2334 basis functions).