Hans Lischka

MOLCAS 8: New Capabilities for Multiconfigurational Quantum Chemical Calculations across the Periodic Table

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas–Kroll–Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC‐PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large‐scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

PNO–CI (pair natural orbital configuration interaction) and CEPA–PNO (coupled electron pair approximation with pair natural orbitals) calculations of molecular systems. II. The molecules BeH2, BH, BH3, CH4, CH−3, NH3 (planar and pyramidal), H2O, OH+3, HF and the Ne atom

PNO–CI and CEPA–PNO calculations are performed for the molecules MgH2, AlH, AlH3, SiH4, PH3, H2S, HCl, and the Ar atom. Two types of Gaussian basis sets are used; both sets contain one p‐set on H. The ’’small’’ set includes one d‐set on the heavy atom, the ’’standard’’ basis two d‐sets and one f‐set. Both for MgH2 and Ar, a ’’large’’ and a ’’very large’’ basis are used as well, which contain additional polarization functions. The energy improvement due to the different polarization functions is analyzed. Hartree–Fock limits for the molecular energies are estimated. The computed valence shell correlation energies are analyzed in terms of quantities defined in part I, in particular in terms of the IEPA (independent electron pair) correlation energies eμIEPA and the error ΔEIEPA of the IEPA scheme. Both the valence shell interorbital pair correlation energies and the IEPA error are smaller in absolute value than those of the corresponding first row hydrides, provided that one uses the localized representatio...

Relaxation mechanisms of UV-photoexcited DNA and RNA nucleobases

A comprehensive effort in photodynamical ab initio simulations of the ultrafast deactivation pathways for all five nucleobases adenine, guanine, cytosine, thymine, and uracil is reported. These simulations are based on a complete nonadiabatic surface-hopping approach using extended multiconfigurational wave functions. Even though all five nucleobases share the basic internal conversion mechanisms, the calculations show a distinct grouping into purine and pyrimidine bases as concerns the complexity of the photodynamics. The purine bases adenine and guanine represent the most simple photodeactivation mechanism with the dynamics leading along a diabatic ππ* path directly and without barrier to the conical intersection seam with the ground state. In the case of the pyrimidine bases, the dynamics starts off in much flatter regions of the ππ* energy surface due to coupling of several states. This fact prohibits a clear formation of a single reaction path. Thus, the photodynamics of the pyrimidine bases is much richer and includes also nπ* states with varying importance, depending on the actual nucleobase considered. Trapping in local minima may occur and, therefore, the deactivation time to the ground state is also much longer in these cases. Implications of these findings are discussed (i) for identifying structural possibilities where singlet/triplet transitions can occur because of sufficient retention time during the singlet dynamics and (ii) concerning the flexibility of finding other deactivation pathways in substituted pyrimidines serving as candidates for alternative nucleobases.

Analytic evaluation of nonadiabatic coupling terms at the MR-CI level. I. Formalism

An efficient and general method for the analytic computation of the nonandiabatic coupling vector at the multireference configuration interaction (MR-CI) level is presented. This method is based on a previously developed formalism for analytic MR-CI gradients adapted to the use for the computation of nonadiabatic coupling terms. As was the case for the analytic energy gradients, very general, separate choices of invariant orbital subspaces at the multiconfiguration self-consistent field and MR-CI levels are possible, allowing flexible selections of MR-CI wave functions. The computational cost for the calculation of the nonadiabatic coupling vector at the MR-CI level is far below the cost for the energy calculation. In this paper the formalism of the method is presented and in the following paper [Dallos et al., J. Chem. Phys. 120, 7330 (2004)] applications concerning the optimization of minima on the crossing seam are described.

Analytic evaluation of nonadiabatic coupling terms at the MR-CI level. II. Minima on the crossing seam: Formaldehyde and the photodimerization of ethylene

The method for the analytic calculation of the nonadiabatic coupling vector at the multireference configuration-interaction (MR-CI) level and its program implementation into the COLUMBUS program system described in the preceding paper [Lischka et al., J. Chem. Phys. 120, 7322 (2004)] has been combined with automatic searches for minima on the crossing seam (MXS). Based on a perturbative description of the vicinity of a conical intersection, a Lagrange formalism for the determination of MXS has been derived. Geometry optimization by direct inversion in the iterative subspace extrapolation is used to improve the convergence properties of the corresponding Newton-Raphson procedure. Three examples have been investigated: the crossing between the 1(1)B1/2(1)A1 valence states in formaldehyde, the crossing between the 2(1)A1/3(1)A1 pi-pi* valence and ny-3py Rydberg states in formaldehyde, and three crossings in the case of the photodimerization of ethylene. The methods developed allow MXS searches of significantly larger systems at the MR-CI level than have been possible before and significantly more accurate calculations as compared to previous complete-active space self-consistent field approaches.

A general multireference configuration interaction gradient program

An efficient and general method for the computation of analytic energy gradients and energy response properties for general MRCI (multireference configuration interaction) and ACPF (averaged coupled pair functional) wave functions is presented. This methodology includes a general approach, based on successive orbital transformations, for the inclusion of the effects of various orbital resolution (canonicalization) constraints. Initial implementation in the columbus Program System demonstrates, particularly for large‐scale multireference wave functions, that the additional computational effort required for the energy gradient is a small fraction of that required for the energy. For polyatomic molecules, the computational resources required for the energy gradient do not depend explicitly on the number of constituent atoms. This combination of features represents a major step forward in the computation and characterization of molecular potential energy surfaces.

Nonadiabatic Deactivation of 9H-Adenine: A Comprehensive Picture Based on Mixed Quantum−Classical Dynamics

Mixed quantum-classical dynamics simulations at the multireference configuration interaction (MR-CIS) level were performed for 9 H-adenine in order to understand its ultrafast nonradiative decay process. Dynamics simulations were also performed for the model system 6-aminopyrimidine. MR-CIS and complete active space perturbation theory (CASPT2) have been employed to characterize a large variety of qualitatively different conical intersections, the branches of the crossing seam connecting them, and the reaction paths from the Franck-Condon region for 9 H-adenine. The results show a two-step process consisting of ultrashort deactivation from S 3 to S 1 and a longer exponential decay step corresponding to the conversion from S 1 to S 0.

Critical appraisal of excited state nonadiabatic dynamics simulations of 9H-adenine

In spite of the importance of nonadiabatic dynamics simulations for the understanding of ultrafast photo-induced phenomena, simulations based on different methodologies have often led to contradictory results. In this work, we proceed through a comprehensive investigation of on-the-fly surface-hopping simulations of 9H-adenine in the gas phase using different electronic structure theories (ab initio, semi-empirical, and density functional methods). Simulations that employ ab initio and semi-empirical multireference configuration interaction methods predict the experimentally observed ultrafast deactivation of 9H-adenine with similar time scales, however, through different internal conversion channels. Simulations based on time-dependent density functional theory with six different hybrid and range-corrected functionals fail to predict the ultrafast deactivation. The origin of these differences is analyzed by systematic calculations of the relevant reaction pathways, which show that these discrepancies can always be traced back to topographical features of the underlying potential energy surfaces.

Analysis of Excitonic and Charge Transfer Interactions from Quantum Chemical Calculations

A procedure for a detailed analysis of excited states in systems of interacting chromophores is proposed. By considering the one-electron transition density matrix, a wealth of information is recovered that may be missed by manually analyzing the wave function. Not only are the position and spatial extent given, but insight into the intrinsic structure of the exciton is readily obtained as well. For example, the method can differentiate between excitonic and charge resonance interactions even in completely symmetric systems. Four examples are considered to highlight the utility of the approach: interactions between the nπ* states in a formaldehyde dimer, excimer formation in the naphthalene dimer, stacking interaction in an adenine dimer, and the excitonic band structure in a conjugated phenylenevinylene oligomer.

The Multiradical Character of One- and Two-Dimensional Graphene Nanoribbons

When is an acene stable? The pronounced multiradical character of graphene nanoribbons of different size and shape was investigated with high-level multireference methods. Quantitative information based on the number of effectively unpaired electrons leads to specific estimates of the chemical stability of graphene nanostructures.