Biomedical subjects

Henrik Koch

Publications and source records attributed to Henrik Koch.

4 recordsLinked to original sources

PubMed2006

Fast noniterative orbital localization for large molecules.

We use Cholesky decomposition of the density matrix in atomic orbital basis to define a new set of occupied molecular orbital coefficients. Analysis of the resulting orbitals ("Cholesky molecular orbitals") demonstrates their localized character inherited from the sparsity of the density matrix. Comparison with the results of traditional iterative localization schemes shows minor differences with respect to a number of suitable measures of locality, particularly the scaling with system size of orbital pair domains used in local correlation methods. The Cholesky procedure for generating orthonormal localized orbitals is noniterative and may be made linear scaling. Although our present implementation scales cubically, the algorithm is significantly faster than any of the conventional localization schemes. In addition, since this approach does not require starting orbitals, it will be useful in local correlation treatments on top of diagonalization-free Hartree-Fock optimization algorithms.

Journal Article↗

PubMed2006

Ab initio potential-energy surface and rovibrational states of the HCN-HCl complex.

A four-dimensional intermolecular potential-energy surface has been calculated for the HCN-HCl complex, with the use of the coupled cluster method with single and double excitations and noniterative inclusion of triples. Data for more than 13,000 geometries were represented by an angular expansion in terms of coupled spherical harmonics; the dependence of the expansion coefficients on the intermolecular distance R was described by the reproducing kernel Hilbert space method. The global minimum with De=1565 cm(-1) and Re=7.47a0 has a linear HCN-HCl hydrogen-bonded structure with HCl as the donor. A secondary hydrogen-bonded equilibrium structure with De=564 cm(-1) and Re=8.21a0 has a T-shaped geometry with HCN as the donor and the acceptor HCl molecule nearly perpendicular to the intermolecular axis. This potential surface was used in a variational approach to compute a series of bound states of the isotopomers HCN-H35Cl, DCN-H35Cl, and HCN-H37Cl for total angular momentum J=0,1,2 and spectroscopic parities e, f. The results could be analyzed in terms of the approximate quantum numbers of a linear polyatomic molecule with two coupled bend modes, plus a quantum number for the intermolecular stretch vibration. They are in good agreement with the recent high resolution spectrum of Larsen et al. [Phys. Chem. Chem. Phys. 7, 1953 (2005)] in the region of 330 cm(-1) corresponding to the HCl libration. The (partly anomalous) effects of isotopic substitutions on the properties of the complex were explained with the aid of the calculations.

Journal Article↗

PubMed2005

Ab initio calculation of optical rotation in (P)-(+)-[4]triangulane.

Optical rotation, the angle through which plane-polarized light rotates when passed through an enantiomerically pure medium, plays a vital role in the determination of the absolute configurations of chiral molecules such as natural products. We describe new quantum mechanical methodology designed to assist in this endeavor by providing high-accuracy computational optical rotatory dispersion data for matching to experimental results. Comparison between theory and experiment for the rigid, helical molecule trispiro[2.0.0.2.1.1]nonane [also known as (P)-(+)-[4]triangulane], recently synthesized with enantiomeric purity, shows that the coupled cluster quantum chemical model provides superb agreement for optical rotation across a wide range of wavelengths (589-365 nm), with errors averaging only 1%.

Journal Article↗

PubMed2004

Polarizability and optical rotation calculated from the approximate coupled cluster singles and doubles CC2 linear response theory using Cholesky decompositions.

A new implementation of the approximate coupled cluster singles and doubles CC2 linear response model using Cholesky decomposition of the two-electron integrals is presented. Significantly reducing storage demands and computational effort without sacrificing accuracy compared to the conventional model, the algorithm is well suited for large-scale applications. Extensive basis set convergence studies are presented for the static and frequency-dependent electric dipole polarizability of benzene and C60, and for the optical rotation of CNOFH2 and (-)-trans-cyclooctene (TCO). The origin-dependence of the optical rotation is calculated and shown to persist for CC2 even at basis set convergence.

Journal Article↗