Engineering topics

Van Neck, Dimitri

Publications and source records attributed to Van Neck, Dimitri.

5 records

DOE OSTI2016

DMRG-CASPT2 study of the longitudinal static second hyperpolarizability of all-trans polyenes

We have implemented internally contracted complete active space second order perturbation theory (CASPT2) with the density matrix renormalization group (DMRG) as active space solver [Y. Kurashige and T. Yanai, J. Chem. Phys. 135, 094104 (2011)]. Internally contracted CASPT2 requires to contract the generalized Fock matrix with the 4-particle reduced density matrix (4-RDM) of the reference wavefunction. The required 4-RDM elements can be obtained from 3-particle reduced density matrices (3-RDM) of different wavefunctions, formed by symmetry-conserving single-particle excitations op top of the reference wavefunction. In our spin-adapted DMRG code CHEMPS2https://github.com/sebwouters/chemps2, we decompose these excited wavefunctions as spin-adapted matrix product states and calculate their 3-RDM in order to obtain the required contraction of the generalized Fock matrix with the 4-RDM of the reference wavefunction. In this work, we study the longitudinal static second hyperpolarizability of all-trans polyenes C{sub 2n}H{sub 2n+2} [n = 4–12] in the cc-pVDZ basis set. DMRG-SCF and DMRG-CASPT2 yield substantially lower values and scaling with system size compared to RHF and MP2, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Normal mode analysis of macromolecular systems with the mobile block Hessian method

Until recently, normal mode analysis (NMA) was limited to small proteins, not only because the required energy minimization is a computationally exhausting task, but also because NMA requires the expensive diagonalization of a 3N{sub a}×3N{sub a} matrix with N{sub a} the number of atoms. A series of simplified models has been proposed, in particular the Rotation-Translation Blocks (RTB) method by Tama et al. for the simulation of proteins. It makes use of the concept that a peptide chain or protein can be seen as a subsequent set of rigid components, i.e. the peptide units. A peptide chain is thus divided into rigid blocks with six degrees of freedom each. Recently we developed the Mobile Block Hessian (MBH) method, which in a sense has similar features as the RTB method. The main difference is that MBH was developed to deal with partially optimized systems. The position/orientation of each block is optimized while the internal geometry is kept fixed at a plausible - but not necessarily optimized - geometry. This reduces the computational cost of the energy minimization. Applying the standard NMA on a partially optimized structure however results in spurious imaginary frequencies and unwanted coordinate dependence. The MBH avoids these unphysical effects by taking into account energy gradient corrections. Moreover the number of variables is reduced, which facilitates the diagonalization of the Hessian. In the original implementation of MBH, atoms could only be part of one rigid block. The MBH is now extended to the case where atoms can be part of two or more blocks. Two basic linkages can be realized: (1) blocks connected by one link atom, or (2) by two link atoms, where the latter is referred to as the hinge type connection. In this work we present the MBH concept and illustrate its performance with the crambin protein as an example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

A hybrid configuration interaction treatment based on seniority number and excitation schemes

We present a configuration interaction method in which the Hamiltonian of an N-electron system is projected on Slater determinants selected according to the seniority-number criterion along with the traditional excitation-based procedure. This proposed method is especially useful to describe systems which exhibit dynamic (weak) correlation at determined geometric arrangements (where the excitation-based procedure is more suitable) but show static (strong) correlation at other arrangements (where the seniority-number technique is preferred). The hybrid method amends the shortcomings of both individual determinant selection procedures, yielding correct shapes of potential energy curves with results closer to those provided by the full configuration interaction method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

Communication: DMRG-SCF study of the singlet, triplet, and quintet states of oxo-Mn(Salen)

We use CHEMPS2, our free open-source spin-adapted implementation of the density matrix renormalization group (DMRG) [S. Wouters, W. Poelmans, P. W. Ayers, and D. Van Neck, Comput. Phys. Commun. 185, 1501 (2014)], to study the lowest singlet, triplet, and quintet states of the oxo-Mn(Salen) complex. We describe how an initial approximate DMRG calculation in a large active space around the Fermi level can be used to obtain a good set of starting orbitals for subsequent complete-active-space or DMRG self-consistent field calculations. This procedure mitigates the need for a localization procedure, followed by a manual selection of the active space. Per multiplicity, the same active space of 28 electrons in 22 orbitals (28e, 22o) is obtained with the 6-31G{sup *}, cc-pVDZ, and ANO-RCC-VDZP basis sets (the latter with DKH2 scalar relativistic corrections). Our calculations provide new insight into the electronic structure of the quintet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2011

Faddeev random-phase approximation for molecules

The Faddeev random-phase approximation is a Green's function technique that makes use of Faddeev equations to couple the motion of a single electron to the two-particle-one-hole and two-hole-one-particle excitations. This method goes beyond the frequently used third-order algebraic diagrammatic construction method: all diagrams involving the exchange of phonons in the particle-hole and particle-particle channel are retained, but the phonons are now described at the level of the random-phase approximation, which includes ground-state correlations, rather than at the Tamm-Dancoff approximation level, where ground-state correlations are excluded. Previously applied to atoms, this paper presents results for small molecules at equilibrium geometry.

74 ATOMIC AND MOLECULAR PHYSICS↗