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.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

For an ab initio calculation of the magnetic excitations: RelaxSE!

In this paper, we present a novel efficient and parallel implementation, RelaxSE, for the calculation of the low-lying excited states and energies of strongly correlated systems. RelaxSE is based on the fully uncontracted multi-reference method of Selected Active Space + Single excitations. This method has been specifically designed to be able to tackle systems with numerous open shells per atoms. It is, however, computationally challenging due to the rapid scaling of the number of determinants and their non-trivial ordering induced by the selection process. We propose a combined determinant-driven and integral-driven approach designed for hybrid OpenMP/MPI parallelization. The performances of RelaxSE are evaluated on a controlled test set and show linear scaling with respect to the number of determinants and a small overhead due to the parallelization. Systems with up to 1 × 10⁹ determinants are successfully computed.

Pressure-dependent X-ray diffraction of the multiferroics RMn2O5

The room-temperature structural properties of the RMn₂O₅ multiferroics have been investigated under pressure, using powder X-ray scattering and density functional theory (DFT) calculations. It was possible to determine the lattice parameters and the main atomic positions as a function of pressure. Good agreement was observed between the X-ray and DFT results for most of the determined crystallographic data. From the DFT calculations, it was possible to infer the pressure evolution of the exchange interactions, and this analysis led to the conclusion that the onset of the q = (½, 0, ½) magnetic structure under pressure is related to the increase in the J₁ super-exchange terms (due to the reduction in the Mn-O distances) compared with the Mn-R exchange interactions. In addition, the 1D antiferromagnetic character of the compounds should be reinforced under pressure.

Divide-Expand-Consolidate Second-Order Møller-Plesset Theory with Periodic Boundary Conditions

We present a generalization of the divide-expand-consolidate (DEC) framework for local coupled-cluster calculations to periodic systems and test it at the second-order Møller-Plesset (MP2) level of theory. For simple model systems with periodicity in one, two, and three dimensions, comparisons with extrapolated molecular calculations and the local MP2 implementation in the Cryscor program show that the correlation energy errors of the extended DEC (X-DEC) algorithm can be controlled through a single parameter, the fragment optimization threshold. Two computational bottlenecks are identified: the size of the virtual orbital spaces and the number of pair fragments required to achieve a given accuracy of the correlation energy. For the latter, we propose an affordable algorithm based on cubic splines interpolation of a limited number of pair-fragment interaction energies to determine a pair cutoff distance in accordance with the specified fragment optimization threshold.

Excitation energies along a range-separated adiabatic connection

We present a study of the variation of total energies and excitation energies along a range-separated adiabatic connection. This connection links the non-interacting Kohn-Sham electronic system to the physical interacting system by progressively switching on the electron-electron interactions whilst simultaneously adjusting a one-electron effective potential so as to keep the ground-state density constant. The interactions are introduced in a range-dependent manner, first introducing predominantly long-range, and then all-range, interactions as the physical system is approached, as opposed to the conventional adiabatic connection where the interactions are introduced by globally scaling the standard Coulomb interaction. Reference data are reported for the He and Be atoms and the H2 molecule, obtained by calculating the short-range effective potential at the full configuration-interaction level using Lieb's Legendre-transform approach. As the strength of the electron-electron interactions increases, the excitation energies, calculated for the partially interacting systems along the adiabatic connection, offer increasingly accurate approximations to the exact excitation energies. Importantly, the excitation energies calculated at an intermediate point of the adiabatic connection are much better approximations to the exact excitation energies than are the corresponding Kohn-Sham excitation energies. This is particularly evident in situations involving strong static correlation effects and states with multiple excitation character, such as the dissociating H2 molecule. These results highlight the utility of long-range interacting reference systems as a starting point for the calculation of excitation energies and are of interest for developing and analyzing practical approximate range-separated density-functional methodologies.

Solid-phase synthesis of peptide-viologen conjugates

This paper presents a robust method for the conjugation of viologens to peptides using an amide coupling strategy that is compatible with standard Fmoc solid-phase peptide synthesis. Methodology is presented for monitoring the milligram scale process quantitatively by UV spectroscopy. This chemistry enables the synthesis of a broad range of asymmetric viologens in high yield at room temperature and is compatible with a wide range of functional groups, including amine, guanidinyl, thiol, carboxylic acid, phenol, and indole.