The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas 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 OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations

Applications to metal K pre-edges of transitionmetal dimers illustrate the approximate origin independence for the intensities in the length representation

X-ray absorption spectroscopy (XAS) in the metal K pre-edge is a standard probe of electronic and geometric structure of transition metal complexes. Simulating the K pre-edge spectra requires contributions beyond the electric dipole, but if that term is non-zero, the second-order terms, e. g. electric quadrupoles, are no longer origin-independent. In the velocity representation, complete origin independence can be achieved by including all terms to the same order in the oscillator strength. Here, we implement that approach in the length representation and use it for restricted active space (RAS) simulations of metal K pre-edges of iron monomers and dimers. Complete origin independence is not achieved and the size of the remaining errors depends on the electric dipole oscillator strength and its ratio in length and velocity representations. The error in the origin independence is in the ANO basis sets two orders of magnitude smaller than the value of the individual contributions. For systemswith strong electric dipole contributions, the errors are not significant within 3 angstrom from a metal centre, far enough to handlemany multi-metal systems. Furthermore, we discuss the convergence of the multipole expansion, the possibility to assign spectral contributions, and the origin of negative absorption intensities. [GRAPHICS

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 alpha (2p -> 1s) or K beta (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. Compared to experiment, relative energies of the two [Fe(CN)(6)](3-) resonances deviate by 0.2 eV in both incident energy and energy transfer directions, and multiplet splittings in [Fe(CN)(6)](4-) are reproduced within 0.1 eV. These values are similar to what can be expected for valence excitations. The development opens the modeling of hard X-ray scattering processes for both solution catalysts and enzymatic systems

OpenMolcas : From Source Code to Insight

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