In this Article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density matrix renormalization group (DMRG) methods, and hybrid multiconfigurational wave function and density functional theory models. Some of these implementations include an array of additional options and functionalities. The paper proceeds and describes developments related to explorations of potential energy surfaces. Here we present methods for the optimization of conical intersections, the simulation of adiabatic and nonadiabatic molecular dynamics, and interfaces to tools for semiclassical and quantum mechanical nuclear dynamics. Furthermore, the Article describes features unique to simulations of spectroscopic and magnetic phenomena such as the exact semiclassical description of the interaction between light and matter, various X-ray processes, magnetic circular dichroism, and properties. Finally, the paper describes a number of built-in and add-on features to support the OpenMolcas platform with postcalculation analysis and visualization, a multiscale simulation option using frozen-density embedding theory, and new electronic and muonic basis sets.

/pdf/openmolcas-from-source-code-to-insight-1v2zixllur.pdf

OpenMolcas : From Source Code to Insight

A series of two-coordinate complexes of iron(II) were prepared and studied for single-molecule magnet behavior. Five of the compounds, Fe[N(SiMe3)(Dipp)]2 (1), Fe[C(SiMe3)3]2 (2), Fe[N(H)Ar′]2 (3), Fe[N(H)Ar*]2 (4), and Fe(OAr′)2 (5) feature a linear geometry at the FeII center, while the sixth compound, Fe[N(H)Ar#]2 (6), is bent with an N–Fe–N angle of 140.9(2)° (Dipp = C6H3-2,6-Pri2; Ar′ = C6H3-2,6-(C6H3-2,6-Pri2)2; Ar* = C6H3-2,6-(C6H2-2,4,6-Pri2)2; Ar# = C6H3-2,6-(C6H2-2,4,6-Me3)2). Ac magnetic susceptibility data for all compounds revealed slow magnetic relaxation under an applied dc field, with the magnetic relaxation times following a general trend of 1 > 2 > 3 > 4 > 5 ≫ 6. Arrhenius plots created for the linear complexes were fit by employing a sum of tunneling, direct, Raman, and Orbach relaxation processes, resulting in spin reversal barriers of Ueff = 181, 146, 109, 104, and 43 cm−1 for 1–5, respectively. CASSCF/NEVPT2 calculations on the crystal structures were performed to explore the influence of deviations from rigorous D∞h geometry on the d-orbital splittings and the electronic state energies. Asymmetry in the ligand fields quenches the orbital angular momentum of 1–6, but ultimately spin–orbit coupling is strong enough to compensate and regenerate the orbital moment. The lack of simple Arrhenius behavior in 1–5 can be attributed to a combination of the asymmetric ligand field and the influence of vibronic coupling, with the latter possibility being suggested by thermal ellipsoid models to the diffraction data.

Slow magnetization dynamics in a series of two-coordinate iron(II) complexes

Physical Content, Analytical Derivation, and Rigorous Extraction of Magnetic Hamiltonians Jean Paul Malrieu,† Rosa Caballol,‡ Carmen J. Calzado, Coen de Graaf,‡,∥ and Nathalie Guiheŕy*,† †Laboratoire de Chimie et Physique Quantiques, Universite ́ de Toulouse 3, 118 route de Narbonne, 31062 Toulouse, France ‡Departament de Química Física i Inorgaǹica, Universitat Rovira i Virgili, Marcel·lí Domingo s/n, 43007 Tarragona, Spain Departamento de Química Física, Universidad de Sevilla, Profesor Garcia Gonzalez s/n, 41012 Sevilla, Spain Institucio ́ Catalana de Recerca i Estudis Avanca̧ts (ICREA), Passeig Lluis Companys 23, 08010 Barcelona, Spain

Magnetic interactions in molecules and highly correlated materials: physical content, analytical derivation, and rigorous extraction of magnetic Hamiltonians.

We present a second-order perturbation theory based on a density matrix renormalization group self-consistent field (DMRG-SCF) reference function. The method reproduces the solution of the complete active space with second-order perturbation theory (CASPT2) when the DMRG reference function is represented by a sufficiently large number of renormalized many-body basis, thereby being named DMRG-CASPT2 method. The DMRG-SCF is able to describe non-dynamical correlation with large active space that is insurmountable to the conventional CASSCF method, while the second-order perturbation theory provides an efficient description of dynamical correlation effects. The capability of our implementation is demonstrated for an application to the potential energy curve of the chromium dimer, which is one of the most demanding multireference systems that require best electronic structure treatment for non-dynamical and dynamical correlation as well as large basis sets. The DMRG-CASPT2/cc-pwCV5Z calculations were performed with a large (3d double-shell) active space consisting of 28 orbitals. Our approach using large-size DMRG reference addressed the problems of why the dissociation energy is largely overestimated by CASPT2 with the small active space consisting of 12 orbitals (3d4s), and also is oversensitive to the choice of the zeroth-order Hamiltonian.

Second-order perturbation theory with a density matrix renormalization group self-consistent field reference function: theory and application to the study of chromium dimer

First principles approach to the electronic structure, magnetic anisotropy and spin relaxation in mononuclear 3d-transition metal single molecule magnets

A formulation of the n-electron valence state perturbation theory (NEVPT) at the third order of perturbation is presented. The present implementation concerns the so-called strongly contracted variant of NEVPT, where only a subspace of the first-order interacting space is taken into account. The resulting strongly contracted NEVPT3 approach is discussed in three test cases: (a) the energy difference between the 3B1 and 1A1 states of the methylene molecule, (b) the potential-energy curve of the N2 molecule ground state, and (c) the chromium dimer (Cr2) ground-state potential-energy profile. Particular attention is devoted to the last case where large basis sets comprising also h orbitals are adopted and where remarkable differences between the second- and third-order results show up.

Third-order multireference perturbation theory: the n-electron valence state perturbation-theory approach.

The potential energy curves of the molecules Mo2, W2, and CrMo have been studied ab initio using large basis sets and the “n-electron valence state perturbation theory” up to the third order in the energy. The third order results for Mo2 and W2 reproduce the equilibrium distances re and the harmonic frequencies ωe in fairly good accordance with the experimental values but tend to underestimate the dissociation energy. The CrMo molecule, for which experimental dissociation energy data do not exist yet, is predicted to have a value for De of ≈2.5eV.

Ground states of the Mo2, W2, and CrMo molecules: a second and third order multireference perturbation theory study.

The potential energy curves of the ground state and of some excited states of the manganese dimer have been calculated over a wide range of internuclear distances using the second order n-electron valence state perturbation theory applied to a complete active space self-consistent field reference wave function. The ground state of Mn(2), for which also the third order NEVPT has been used, is calculated to be a singlet belonging to the Sigma(g) (+) symmetry, characterized by a large equilibrium internuclear distance R(e) of 3.7-3.8 A, by a low dissociation energy D(e) of 0.07-0.08 eV, and by a small harmonic frequency omega(e) of 43 cm(-1). The experimental evidence that Mn(2) is a van der Waals molecule is thus confirmed. Among the excited states, (11)Pi(u), which is usually indicated as the ground state by density functional theory studies, appears as a low-lying state with R(e)=2.50 A, D(e)=1.35 eV, and omega(e)=246 cm(-1).

An ab initio multireference perturbation theory study on the manganese dimer

A short summary of the n-electron valence state perturbation theory (NEVPT) is provided, and the defects inherent in cases of quasi-degeneration are underscored. A new version of the theory, called QD-NEVPT, based on quasi-degenerate perturbation theory, is presented. The new theory introduces an effective interaction among the zero-order wave functions and overcomes the difficulties of state-specific NEVPT. An application concerning the interaction among the π π* and Rydberg singlets of formaldehyde is presented. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006

Alex Cavallini

Papers

Third-order multireference perturbation theory: the n-electron valence state perturbation-theory approach.

Ground states of the Mo2, W2, and CrMo molecules: a second and third order multireference perturbation theory study.

An ab initio multireference perturbation theory study on the manganese dimer

Developments in the n‐electron valence state perturbation theory