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

We present a wave-function based method to solve the time-dependent many-electron Schrödinger equation (TDSE) with special emphasis on strong-field ionization phenomena. The theory builds on the configuration-interaction (CI) approach supplemented by the generalized-active-space (GAS) concept from quantum chemistry. The latter allows for a controllable reduction in the number of configurations in the CI expansion by imposing restrictions on the active orbital space. The method is similar to the recently formulated time-dependent restricted-active-space (TD-RAS) CI method [D. Hochstuhl, and M. Bonitz, Phys. Rev. A 86, 053424 (2012)]. We present details of our implementation and address convergence properties with respect to the active spaces and the associated account of electron correlation in both ground state and excitation scenarios. We apply the TD-GASCI theory to strong-field ionization of polar diatomic molecules and illustrate how the method allows us to uncover a strong correlation-induced shift of the preferred direction of emission of photoelectrons.

Simulations of iron K pre-edge X-ray absorption spectra using the core restricted active space method.Physical Chemistry, http://dx.doi.org/10.1039/c5cp07487hAccess to the published version may require subscription. N.B. When citing this work, cite the original published paper. Permanent link to this version:http://urn.kb.se/resolve?urn=urn:nbn:se:uu:diva-243571Simulations of iron K pre-edge X-ray absorption spectra using the restricted active space method † Meiyuan Guo, a Lasse Kragh Sørensen, a Mickaël G. Delcey, a,b Rahul V. Pinjari, a,c and Marcus Lundberg * aThe intensities and relative energies of metal K pre-edge features are sensitive to both geometric and electronic structure. With the possibility to collect high-resolution spectral data it important to find theoretical methods that include all important spectral effects: ligand-field splitting, multiplet structures, 3d-4p orbital hybridization, and charge-transfer excitations. Here the restricted active space (RAS) method is used for the first time to calculate metal K pre-edge spectra of open-shell systems, and its performance is tested against six iron complexes: [FeCl 6 ] n-, [FeCl 4 ] n-, and [Fe(CN) 6 ] n-in ferrous and ferric oxidation states. The method gives good descriptions of the spectral shapes for all six systems. The mean absolute deviation for the relative energies of different peaks is only 0.1 eV. For the two systems that lack centrosymmetry [FeCl 4 ] 2-/1-, the ratios between dipole and quadrupole intensity contributions are reproduced with an error of 10%, which leads to good descriptions of the integrated pre-edge intensities. To gain further chemical insight, the origins of the pre-edge features have been analyzed with a chemically intuitive molecular orbital picture that serves as a bridge between the spectra and the electronic structures. The RAS method can thus be used to predict and rationalize the effects of changes in both oxidation state and ligand environment in a number of hard X-ray studies of small and medium-sized molecular systems.

The developments of the open-source chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes can address, while showing that is an attractive platform for state-of-the-art atomistic computer simulations.

Electronically excited states play important roles in many chemical reactions and spectroscopic techniques. In quantum chemistry, a common technique to solve excited states is the multiroot Davidson algorithm, but it is not designed for processes like X‐ray spectroscopy that involves hundreds of highly excited states. We show how the use of a restricted active space wavefunction together with a projection operator to remove low‐lying electronic states offers an efficient way to reach single and double‐core‐hole states. Additionally, several improvements to the stability and efficiency of the configuration interaction (CI) algorithm for a large number of states are suggested. When applied to a series of transition metal complexes the new CI algorithm does not only resolve divergence issues but also leads to typical reduction in computational time by 70%, with the largest savings for small molecules and large active spaces. Together, the projection operator and the improved CI algorithm now make it possible to simulate a wide range of single‐ and two‐photon spectroscopies. © 2019 Wiley Periodicals, Inc.

Molecular orbital simulations of metal 1s2p resonant inelastic X-ray scattering. For first-row transition metals, high-resolution 3d electronic structure information can be obtained using resonant inelastic X-ray scattering (RIXS). In the hard X-ray region, a K pre-edge (1s → 3d) excitation can be followed by monitoring the dipole-allowed Kα (2p → 1s) or Kβ (3p → 1s) emission, processes labeled 1s2p or 1s3p RIXS. Here the restricted active space (RAS) approach, which is a molecular orbital method, is used for the first time to study hard X-ray RIXS processes. This is achieved by including the two sets of core orbitals in different partitions of the active space. Transition intensities are calculated using both first-and second-order expansions of the wave vector, including, but not limited to, electric dipoles and quadrupoles. The accuracy of the approach is tested for 1s2p RIXS of iron hexacyanides [Fe(CN) 6 ] n-in ferrous and ferric oxidation states. RAS simulations accurately describe the multiplet structures and the role of 2p and 3d spin-orbit coupling on energies and selection rules. Journal of Physical Chemistry

PostprintThis is the accepted version of a paper published in Molecular Physics. This paper has been peer-reviewed but does not include the final publisher proof-corrections or journal pagination. Citation for the original published paper (version of record):Sørensen, L K., Guo, M., Lindh, R., Lundberg, M. (2017) Applications to metal K pre-edges of transitionmetal dimers illustrate the approximate origin independence for the intensities in the length representation Molecular Physics, 115(1-2): 174-189

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 post calculation analysis and visualization, a multiscale simulation option using frozen-density embedding theory and new electronic and muonic basis sets.