Showing papers by "Christof Hättig published in 2006"

Distributed memory parallel implementation of energies and gradients for second-order Møller–Plesset perturbation theory with the resolution-of-the-identity approximation

Abstract: We present a parallel implementation of second-order Moller–Plesset perturbation theory with the resolution-of-the-identity approximation (RI-MP2). The implementation is based on a recent improved sequential implementation of RI-MP2 within the Turbomole program package and employs the message passing interface (MPI) standard for communication between distributed memory nodes. The parallel implementation extends the applicability of canonical MP2 to considerably larger systems. Examples are presented for full geometry optimizations with up to 60 atoms and 3300 basis functions and MP2 energy calculations with more than 200 atoms and 7000 basis functions.

Benchmarking two-photon absorption with CC3 quadratic response theory, and comparison with density-functional response theory.

Abstract: We present a detailed study of the effects of electron correlation on two-photon absorption calculated by coupled cluster quadratic response theory. The hierarchy of coupled cluster models CCS, CC2, CCSD, and CC3 has been used to investigate the effects of electron correlation on the two-photon absorption cross sections of formaldehyde (CH2O), diacetylene (C4H2), and water (H2O). In particular, the effects of triple excitations on two-photon transition cross sections are determined for the first time. In addition, we present a detailed comparison of the coupled cluster results with those obtained from Hartree-Fock and density-functional response theories. We have investigated the local-density approximation, the pure Becke-Lee-Yang-Parr (BLYP) functional, the hybrid Becke-3-parameter-Lee-Yang-Parr (B3LYP), and the Coulomb-attenuated B3LYP (CAM-B3LYP) functionals. Our results show that the CAM-B3LYP functional, when used in conjuction with a one-particle basis-set containing diffuse functions, has much promise; however, care must still be exercised for diffuse Rydberg-type states.

Photophysics of organic photostabilizers. Ab initio study of the excited-state deactivation mechanisms of 2-(2'-hydroxyphenyl)benzotriazole.

Abstract: Excited-state reaction paths and the corresponding energy profiles of 2-(2‘-hydroxyphenyl)benzotriazole (TIN−H) have been determined with the CC2 (simplified singles-and-doubles coupled-cluster) ab initio method. Hydrogen transfer along the intramolecular hydrogen bond, torsion of the aromatic rings and pyramidization of the central nitrogen atom are identified as the most relevant photochemical reaction coordinates. The keto-type planar S1 state reached by barrierless intramolecular hydrogen transfer is found to be unstable with respect to torsion. The latter mode, together with a moderate pyramidization of the central nitrogen atom, provides barrierless access to a S1−S0 conical intersection. Only the π-type orbitals of the aromatic rings are involved in the open-shell structures. The S1−S0 conical intersection, which occurs for perpendicular geometry of the aromatic rings, is a pure biradical. From the conical intersection, a barrierless reaction path steers the system back to the enol-type minimum of ...

Inclusion of the (T) triples correction into the linear‐r12 corrected coupled‐cluster model CCSD(R12)

Abstract: The recently developed singles-and-doubles linear-r12 corrected coupled-cluster model CCSD(R12) has been extended by inclusion of triples corrections in exactly the same manner as in the popular CCSD(T) model. The resulting model is denoted CCSD(T)(R12). An orthonormal auxiliary basis set is used for the resolution-of-the-identity approximation to calculate the three- and four-electron integrals needed in the linear-r12 Ansatz. CCSD(T)(R12) reaction energies are calculated for a test set of 15 chemical reactions involving 23 small molecules. The CCSD(T)(R12), CCSD(R12), and MP2-R12 results are compared with experimental data obtained by correcting experimental enthalpies of reaction for zero-point vibrational energies and finite temperature contributions. Conventional MP2, CCSD, and CCSD(T) results are presented for comparison. An estimate of the basis set limit is obtained by an aug-cc-pV(Q + d)Z/aug-cc-pV(5 + d)Z CCSD(T) extrapolation for a subset of 12 reactions. All CCSD(T) results are compared with the extrapolated values, and the errors are analyzed for the CCSD(T) and CCSD(T)(R12) models in various basis sets. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006

Intramolecular charge-transfer mechanism in quinolidines: the role of the amino twist angle.

Abstract: Quantum-chemical calculations with the approximate coupled-cluster singles-and-doubles model CC2 have been carried out for 1-tert-butyl-6-cyano-1,2,3,4-tetrahydroquinoline (NTC6). For this molecule dual fluorescence was experimentally observed, raising the discussion about the importance of the amino twist angle for this process. The calculations suggest that both the ground state and the normal fluorescent state are significantly twisted by 30°−40° and that the molecule is flexible enough to move into an even stronger twisted conformation (60°−70°) in its intramolecular charge-transfer (ICT) state which is responsible for the anomalous fluorescence band. Such a conformation both minimizes the total energy in the S1 state and maximizes the dipole moment. The barrier from the normal fluorescent state to the ICT state region is very small. Comparison to the situation in the 1-methyl-derivative NMC6 suggests that a large alkyl substituent makes the preferably planar normal fluorescent state energetically unf...

Coupled-cluster response theory with linear-r12 corrections: the CC2-R12 model for excitation energies.

Abstract: Coupled-cluster response theory for vertical excitation energies within the second-order approximate coupled-cluster singles-and-doubles model CC2, including linear-r12 corrections, is derived and implemented for Ansatze 1 and 2 of R12 theory. An orthonormal auxiliary basis set is used for the resolution-of-the-identity approximation in order to calculate the three- and four-electron integrals needed in R12 theory. The basis set convergence is investigated for a selected set of atoms and small molecules and it is found that in many cases the convergence is not improved. An analysis of the different contributions to excitation energies shows that the present scheme for the construction of the R12 pair functions leads in response theory to an unbalanced description of ground- and excited-state wave functions and needs to be generalized to carry the high accuracy of R12 methods over to response theory.

High-order correlation effects on dynamic hyperpolarizabilities and their geometric derivatives: a comparison with density functional results.

Abstract: Second harmonic generation hyperpolarizabilities and their geometric derivatives have been calculated for HF, HCl, CO, and LiF, using the coupled cluster hierarchies, CCS, CC2, CCSD, CC3, and large correlation-consistent basis sets The full configuration interaction results have been used to test the accuracy of the coupled cluster results The CCS and CC2 methods do not improve on the Hartree-Fock results while CCSD is robust and gives significant improvements compared to CCS and CC2 The effects of triples in CC3 are in some cases substantial Higher order correlation effects are significant for LiF Including core-valence correlation effects is required only if high accuracy is desired The coupled cluster results have been used as benchmarks for the results obtained by means of density functional theory using various exchange-correlation functionals For the hyperpolarizability B3LYP was found to perform best, ie, to give the results closest to the CC3 ones, while for the geometric derivatives none of the considered functionals was able to give a consistent description for all the considered molecules

Extensions of r12 corrections to CC2-R12 for excited states.

Abstract: As known since about two decades, R12 methods, which include terms linear in the interelectronic distance r12 in the wave function, improve substantially the basis set convergence of the ground state correlation energy. In a previous study, however, it was found that the same approach does not give a similar systematic improvement if applied to excited states in the framework of coupled cluster response theory. In the present work, we examine the reason for this behavior and show that the inclusion of additional orbitals in the construction of the r12 pair functions leads to an enhanced basis set convergence (and thus a balanced description) also for the excited states.

Femtosecond pump/probe photoelectron spectroscopy of isolated C60 negative ions

Abstract: We have measured pump/probe photoelectron spectra of mass-selected, near room temperature C60− in the gas phase. The lifetime of the vibrationally excited B−(Eg2) state at a calculated energy of 1.26eV was found to be τ=2.2±0.2ps. The dominant decay process corresponds to intramolecular radiationless transitions into ground state C60−. This is in contrast to C60 for which pumping at the absorption onset (1.95eV) leads to predominantly intersystem crossing.

Accurate Nonlinear Optical Properties for Small Molecules

Abstract: During the last decade it became possible to calculate by quantum chemical ab initio methods not only static but also frequency-dependent properties with high accuracy. Today, the most important tools for such calculations are coupled cluster response methods in combination with systematic hierarchies of correlation consistent basis sets. Coupled cluster response methods combine a computationally efficient treatment of electron correlation with a qualitatively correct pole structure and frequency dispersion of the response functions. Both are improved systematically within a hierarchy of coupled cluster models. The present contribution reviews recent advances in the highly accurate calculation of frequency-dependent properties of atoms and small molecules, electronic structure methods, basis set convergence and extrapolation techniques. Reported applications include first and second hyperpolarizabilities, Faraday, Buckingham and Cotton–Mouton effects as well as Jones and magneto-electric birefringence

Microwave and theoretical investigation of the internal rotation in m-cresol

Abstract: The microwave spectrum of m-cresol (3-methylphenol) has been investigated using a molecular beam Fourier transform microwave spectrometer in the frequency range from 3 to 26.5 GHz. The rotation of the hydroxy group into two different unequal energetic minima leads to different spectra for the syn- and anticonformers. Because of a high potential barrier both conformers can be analyzed independently. The methyl group is undergoing an almost free internal rotation which is only hindered by small barriers and splits the vibrational ground state in two states of internal rotation denoted as A and E species. The spacing between the species is found to be up to 10 GHz. The potential for the internal rotation can be determined from the spectra and analyzed in terms of the Fourier components V3 and V6. For syn-m-cresol these parameters were determined as V3=673(3) GHz and V6=-335(24) GHz and for anti-m-cresol V3=95(5) GHz and V6=-416(46) GHz. The barriers to internal rotation were furthermore calculated with second-order Moller-Plesset perturbation theory and second-order coupled-cluster singles- and-doubles model (CC2) in the electronic ground state and with CC2 in the first excited state. The CC2 method is found to be an appropriate method to calculate potential barriers in electronic excited states of such compounds.

Optical properties of N-succinimidyl bithiophene and the effects of the binding to biomolecules: comparison between coupled-cluster and time-dependent density functional theory calculations and experiments.

Abstract: We report a joint theoretical−experimental study on the optical properties of 5-N-succinimidyl-2,2‘-bithiophene (NS-2T), a prototype system for a new class of biomarkers. Time-dependent density functional theory (TD-DFT) and approximate coupled-cluster single and doubles (CC2) calculations are performed in the ground and excited states. Theoretical results are compared with absorption, photoluminescence (PL), time-resolved PL, and PL quantum efficiency measurements. The excited state of NS-2T has a larger dipole moment as compared to that of the ground state, explaining the experimental shift of the PL peak in solvents of different polarity, and a smaller intersystem crossing (ISC) rate as compared to that of isolated bithiophene (2T), explaining the increased PL quantum efficiency. We also studied two model systems to describe the effects of the covalent binding of NS-2T to biomolecules and proteins with the e-NH2 lysine groups. These model systems show optical properties closer to 2T, as the PL quantum ...

Density dependence of electric properties of binary mixtures of inert gases

Abstract: A computational study of the density dependence of the refractivity and dielectric constants, the electric-field induced second harmonic generation (ESHG) hyperpolarizabilities and of the Kerr constants of binary mixtures of helium, neon and argon is presented. Potentials and interaction properties of the homonuclear A2 and heteronuclear AB dimers (A, B = He, Ne, Ar) are taken from a previous study [J. Lopez Cacheiro, B. Fernandez, D. Marchesan, S. Coriani, C. Hattig, A. Rizzo. Molec. Phys., 102, 101 (2004)]. Dispersion coefficients for the second virial coefficients allow for the determination of the density dependence at any frequency far from the lowest resonance. Fully quantum mechanical results are presented and a comparison with the corresponding classical estimates is discussed. Deep minima are predicted to occur in the ESHG second virial coefficient curve drawn as a function of the molar fraction of one of the components in binary mixtures of He/Ar and Ne/Ar. This phenomenon, observed over a wide range of temperatures, should be easily verifiable experimentally.