PubMed Health⌕ Search

Biomedical subjects

Hans-Joachim Werner

Publications and source records attributed to Hans-Joachim Werner.

15 recordsLinked to original sources
PubMed2006

Calculation of smooth potential energy surfaces using local electron correlation methods.

The geometry dependence of excitation domains in local correlation methods can lead to noncontinuous potential energy surfaces. We propose a simple domain merging procedure which eliminates this problem in many situations. The method is applied to heterolytic bond dissociations of ketene and propadienone, to SN2 reactions of Cl(-) with alkylchlorides, and in a quantum mechanical/molecular mechanical study of the chorismate mutase enzyme. It is demonstrated that smooth potentials are obtained in all cases. Furthermore, basis set superposition error effects are reduced in local calculations, and it is found that this leads to better basis set convergence when computing barrier heights or weak interactions. When the electronic structure strongly changes between reactants or products and the transition state, the domain merging procedure leads to a balanced description of all structures and accurate barrier heights.

Journal Article↗
PubMed2006

Calculation of intermolecular interactions in the benzene dimer using coupled-cluster and local electron correlation methods.

Potential energy curves for the parallel-displaced, T-shaped and sandwich structures of the benzene dimer are computed with density fitted local second-order Møller-Plesset perturbation theory (DF-LMP2) as well as with the spin-component scaled (SCS) variant of DF-LMP2. While DF-LMP2 strongly overestimates the dispersion interaction, in common with canonical MP2, the DF-SCS-LMP2 interaction energies are in excellent agreement with the best available literature values along the entire potential energy curves. The DF-SCS-LMP2 dissociation energies for the three structures are also compared with new complete basis set estimates of the interaction energies obtained from accurate coupled cluster (CCSD(T)) and DF-SCS-MP2 calculations. Since LMP2 is essentially free of basis set superposition errors, counterpoise corrections are not required. As a result, DF-SCS-LMP2 is computationally inexpensive and represents an attractive method for the study of larger pi-stacked systems such as truncated sections of DNA.

Journal Article↗
PubMed2006

Application of Gaussian-type geminals in local second-order Møller-Plesset perturbation theory.

In this work Gaussian-type Geminals (GTGs) are applied in local second-order Moller-Plesset perturbation theory to improve the basis set convergence. Our implementation is based on the weak orthogonality functional of Szalewicz et al., [Chem. Phys. Lett. 91, 169 (1982); J. Chem. Phys. 78, 1420 (1983)] and a newly developed program for calculating the necessary many-electron integrals. The local approximations together with GTGs in the treatment of the correlation energy are introduced and tested. First results for correlation energies of H(2)O, CH(4), CO, C(2)H(2), C(2)H(4), H(2)CO, and N(2)H(4) as well as some reaction and activation energies are presented. More than 97% of the valence-shell correlation energy is recovered using aug-cc-pVDZ basis sets and six GTGs per electron pair. The results are compared with conventional calculations using correlation-consistent basis sets as well as with MP2-R12 results.

Journal Article↗
PubMed2006

Accurate potential energy surface and quantum reaction rate calculations for the H+CH4-->H2+CH3 reaction.

Calculations for the cumulative reaction probability N(E) (for J=0) and the thermal rate constant k(T) of the H+CH(4)-->H(2)+CH(3) reaction are presented. Accurate electronic structure calculations and a converged Shepard-interpolation approach are used to construct a potential energy surface which is specifically designed to allow the precise calculation of k(T) and N(E). Accurate quantum dynamics calculations employing flux correlation functions and multiconfigurational time-dependent Hartree wave packet propagation compute N(E) and k(T) based on this potential energy surface. The present work describes in detail the various convergence test performed to investigate the accuracy of the calculations at each step. These tests demonstrate the predictive power of the present calculations. In addition, approximate approaches for reaction rate calculations are discussed. A quite accurate approximation can be obtained from a potential energy surface which includes only interpolation points on the minimum energy path.

Journal Article↗
PubMed2006

Explicitly correlated local second-order perturbation theory with a frozen geminal correlation factor.

The recently introduced MP2-R122*A(loc) and LMP2-R122*A(loc) methods are modified to use a short-range correlation factor expanded as a fixed linear combination of Gaussian geminals. Density fitting is used to reduce the effort for integral evaluation, and local approximations are introduced to improve the scaling of the computational resources with molecular size. The MP2-F122*A(loc) correlation energies converge very rapidly with respect to the atomic orbital basis set size. Already with the aug-cc-pVTZ basis the correlation energies computed for a set of 21 small molecules are found to be within 0.5% of the MP2 basis set limit. Furthermore the short-range correlation factor leads to an improved convergence of the resolution of the identity, and eliminates problems with long-range errors in density fitting caused by the linear r12 factor. The DF-LMP2-F122*A(loc) method is applied to compute second-order correlation energies for molecules with up to 49 atoms and more than 1600 basis functions.

Journal Article↗
PubMed2006

Impact of local and density fitting approximations on harmonic vibrational frequencies.

Harmonic vibrational frequencies are computed using second-order Møller-Plesset perturbation theory (MP2) with and without local (LMP2) and density fitting (DF) approximations. Results for a test set of 17 small and medium size molecules (366 normal modes) are presented, and frequency scaling factors for LMP2 in combination with two different basis sets are determined. Comparison of the MP2 and LMP2 frequencies with experimental data reveals that the introduction of local approximations leads to a slightly better agreement with experiment. This is attributed to the reduction of basis set superposition errors in local calculations. Introduction of DF approximations within the LMP2 formalism leads to negligible deviations but significantly reduces the computational cost. These facts extend the applicability of the method to larger systems with large basis sets. As an example, the method is applied to a full DF-LMP2/cc-pVTZ frequency calculation for testosterone (49 atoms).

Journal Article↗
PubMed2006

Explicitly correlated second-order perturbation theory using density fitting and local approximations.

Three major obstacles in electronic structure theory are the steep scalings of computer time with respect to system size and basis size and the slow convergence of correlation energies in orbital basis sets. Three solutions to these are, respectively, local methods, density fitting, and explicit correlation; in this work, we combine all three to produce a low-order scaling method that can achieve accurate MP2 energies for large systems. The errors introduced by the local approximations into the R12 treatment are analyzed for 16 chemical reactions involving 21 molecules. Weak pair approximations, as well as local resolution of the identity approximations, are tested for molecules with up to 49 atoms, over 100 correlated electrons, and over 1000 basis functions.

Journal Article↗
PubMed2005

Towards accurate ab initio calculations on the vibrational modes of the alkaline earth metal hydrides.

The fundamental modes of the alkaline earth metal hydrides (BeH2, MgH2, CaH2) and their dimers, HX(H)2XH, have been studied by vibrational configuration-interaction calculations based on very accurate potential energy surfaces. Comparison with experimental data obtained from matrix isolation and gas phase measurements is provided and the agreement was found to be excellent for the monomers but poor for the dimers. In addition, many fundamental bands are predicted which have not yet been detected experimentally.

Journal Article↗
PubMed2004

First-principles theory for the H + CH4 --> H2 + CH3 reaction.

A full-dimensional quantum dynamics simulation of a hydrogen atom reacting with methane on an accurate ab initio potential energy surface is reported. Based on first-principles theory, thermal rate constants are predicted with an accuracy comparable to (or even exceeding) experimental precision. The theoretical prediction is within the range of the significantly varied experimental rate constants reported by different groups. This level of accuracy has previously been achieved only for smaller, three-or four-atom reactive systems. Comparison with classical transition state theory confirms the importance of quantum mechanical tunneling for the rate constant below 400 kelvin.

Journal Article↗
PubMed2004

Analytical energy gradients for local second-order Møller-Plesset perturbation theory using density fitting approximations.

An efficient method to compute analytical energy derivatives for local second-order Møller-Plesset perturbation energy is presented. Density fitting approximations are employed for all 4-index integrals and their derivatives. Using local fitting approximations, quadratic scaling with molecular size and cubic scaling with basis set size for a given molecule is achieved. The density fitting approximations have a negligible effect on the accuracy of optimized equilibrium structures or computed energy differences. The method can be applied to much larger molecules and basis sets than any previous second-order Møller-Plesset gradient program. The efficiency and accuracy of the method is demonstrated for a number of organic molecules as well as for molecular clusters. Examples of geometry optimizations for molecules with 100 atoms and over 2000 basis functions without symmetry are presented.

Journal Article↗
PubMed2004

Dynamically weighted multiconfiguration self-consistent field: multistate calculations for F+H2O-->HF+OH reaction paths.

A novel method of dynamically adjusted weighting factors in state-averaged multiconfiguration self-consistent-field calculations (SA-MCSCF) is described that is applicable to systems of arbitrary dimensionality. The proposed dynamically weighted approach automatically weights the relevant electronic states in each region of the potential energy surface, smoothly adjusting between these regions with an energy dependent functional. This method is tested on the F(2P)+H2O-->HF+OH(2Pi) reaction, which otherwise proves challenging to describe with traditional SA-MCSCF methods due to (i) different asymptotic degeneracies of reactant (threefold) and product (twofold) channels, and (ii) presence of low-lying charge transfer configurations near the transition state region. The smoothly varying wave functions obtained by dynamically weighted multiconfigurational self-consistent field represent excellent reference states for high-level multireference configuration interaction calculations and offer an ideal starting point for construction of multiple state potential energy surfaces.

Journal Article↗
PubMed2004

Details and consequences of the nonadiabatic coupling in the Cl(2P) + H2 reaction.

In an investigation of nonadiabaticity in the Cl + H2 reaction we examine the various coupling terms responsible. By neglecting various of these terms, we show that the spin-orbit coupling in the Cl atom is primarily responsible for non Born-Oppenheimer effects in this reaction, and that the anisotropies in the contributing electronic potential energy surfaces as well as Coriolis terms, are considerably weaker in their effect. Neglect of all coupling except spin-orbit leads to little change in the calculated reaction probabilities both for the Born-Oppenheimer allowed [Cl(2P3/2) + H2] and Born-Oppenheimer forbidden [Cl(2P1/2) + H2] reactions. This implies that one can safely neglect the coupling between different values of the electronic projection quantum number omega. Consequently, the computational effort involved in an accurate quantum scattering treatment of this (and, hopefully, other similar reactions) can be substantially reduced.

Journal Article↗
PubMed2003

Ab initio excited-state dynamics of the photoactive yellow protein chromophore.

The photoisomerization mechanism of the neutral form of the photoactive yellow protein (PYP) chromophore is investigated using ab initio quantum chemistry and first-principles nonadiabatic molecular dynamics (ab initio multiple spawning or AIMS). We identify the nature of the two lowest-lying excited states, characterize the short-time behavior of molecules excited directly to S2, and explain the origin of the experimentally observed wavelength-dependent photoisomerization quantum yield.

Bacterial Proteins↗
PubMed2002

Theoretical study of the validity of the Born-Oppenheimer approximation in the Cl + H2 --> HCl + H reaction.

Reactivity of the excited spin-orbit state of Cl with H2 to yield ground-state HCl products is forbidden by the Born-Oppenheimer (BO) approximation. We used new ab initio potential energy surfaces and exact quantum scattering calculations to explore the extent of electronic nonadiabaticity in this reaction. In direct contrast to recent experiments, we predict that the BO-allowed reaction of the ground spin-orbit state will be much more efficient than the BO-forbidden reaction of the excited spin-orbit state. Also, Coriolis coupling opens up an electronically nonadiabatic inelastic channel, which competes substantially with reaction.

Journal Article↗