In this paper, we examine decay and fragmentation of core-excited and core-ionized water molecules combining quantum chemical calculations and electron-energy-resolved electron–ion coincidence spectroscopy. The experimental technique allows us to connect electronic decay from core-excited states, electronic transitions between ionic states, and dissociation of the molecular ion. To this end, we calculate the minimum energy dissociation path of the core-excited molecule and the potential energy surfaces of the molecular ion. Our measurements highlight the role of ultra-fast nuclear motion in the 1a1−14a1 core-excited molecule in the production of fragment ions. OH+ fragments dominate for spectator Auger decay. Complete atomization after sequential fragmentation is also evident through detection of slow H+ fragments. Additional measurements of the non-resonant Auger decay of the core-ionized molecule (1a1−1) to the lower-energy dication states show that the formation of the OH+ + H+ ion pair dominates, whereas sequential fragmentation OH+ + H+ → O + H+ + H+ is observed for transitions to higher dication states, supporting previous theoretical investigations.

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

In this paper, we examine decay and fragmentation of core-excited and core-ionized water molecules combining quantum chemical calculations and electron-energy-resolved electron–ion coincidence spectroscopy. The experimental technique allows us to connect electronic decay from core-excited states, electronic transitions between ionic states, and dissociation of the molecular ion. To this end, we calculate the minimum energy dissociation path of the core-excited molecule and the potential energy surfaces of the molecular ion. Our measurements highlight the role of ultra-fast nuclear motion in the 1a 1 4a1 core-excited molecule in the production of fragment ions. OH + fragments dominate for spectator Auger decay. Complete atomization after sequential fragmentation is also evident through detection of slow H fragments. Additional measurements of the non-resonant Auger decay of the core-ionized molecule (1a 1 ) to the lower-energy dication states show that the formation of the OH + H ion pair dominates, whereas sequential fragmentation OH + H → O + H + H is observed for transitions to higher dication states, supporting previous theoretical investigations. Published under license by AIP Publishing. https://doi.org/10.1063/1.5141414., s

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.

The Complete Active Space Second-Order Perturbation Theory (CASPT2) is well-established as a high-accuracy electronic structure method. It was originally implemented in the early 1990s to an efficient computer code in the molcas program suite, and this implementation has been extensively used as a standard tool. Here, we report a comparison of it against two independent computer-aided implementations of the CASPT2 method, revealing that the CASPT2 energies provided by the original code of molcas (version 8 or earlier) are inconsistent with the predictions of the lately developed computer-aided implementations. It is shown that this error is associated with the projections of the first-order equation onto the fully internally contracted multireference bases which are partially inconsistent between the left- and right-hand sides. The degree of the errors is assessed by performing illustrative CASPT2 calculations. The errors in total CASPT2 energies are demonstrated to be negligible relative to chemical accuracy in many cases, while there is a difficult case where they may substantially alter chemical description. The incorporation of the consistent projections into molcas has been carried out, which is available in the version 8 sp1.

The chromium dimer has long been a benchmark molecule to evaluate the performance of different computational methods ranging from density functional theory to wave function methods. Among the latter, multiconfigurational perturbation theory was shown to be able to reproduce the potential energy surface of the chromium dimer accurately. However, for modest active space sizes, it was later shown that different definitions of the zeroth-order Hamiltonian have a large impact on the results. In this work, we revisit the system for the third time with multiconfigurational perturbation theory, now in order to increase the active space of the reference wave function. This reduces the impact of the choice of zeroth-order Hamiltonian and improves the shape of the potential energy surface significantly. We conclude by comparing our results of the dissocation energy and vibrational spectrum to those obtained from several highly accurate multiconfigurational methods and experiment. For a meaningful comparison, we used the extrapolation to the complete basis set for all methods involved.

This chapter introduces the multiconfigurational SCF (MCSCF) theory. It discusses the multiconfigurational (MC) approach in terms of a partitioning of the exact wave function (in a finite one-particle basis), as represented by a linear combination of all possible electron configurations-the full Configuration Interaction (CI) wave function. It presents a background on the need to develop an MCSCF theory. The chapter devotes to the technique of determining the MCSCF wave function. It discusses some aspects of different representations and higher order derivatives with respect to the wave function parameters. The chapter presents two special types of MCSCF wave functions-the complete active space SCF (CASSCF) and restrictive active space SCF (RASSCF)-that have had tremendous impact on the use of the MCSCF method. The chapter concludes with some general rules to be applied in the selection of active orbitals for CASSCF and RASSCF calculations.

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. © 2015 Wiley Periodicals, Inc.

Human pigmentation is a complex phenomenon commonly believed to serve a photoprotective function through the generation and strategic localization of black insoluble eumelanin biopolymers in sun exposed areas of the body. Despite compelling biomedical relevance to skin cancer and melanoma, eumelanin photoprotection is still an enigma: What makes this pigment so efficient in dissipating the excess energy brought by harmful UV-light as heat? Why has Nature selected 5,6-dihydroxyindole-2-carboxylic acid (DHICA) as the major building block of the pigment instead of the decarboxylated derivative (DHI)? By using pico- and femtosecond fluorescence spectroscopy we demonstrate herein that the excited state deactivation in DHICA oligomers is 3 orders of magnitude faster compared to DHI oligomers. This drastic effect is attributed to their specific structural patterns enabling multiple pathways of intra- and interunit proton transfer. The discovery that DHICA-based scaffolds specifically confer uniquely robust photoprotective properties to natural eumelanins settles a fundamental gap in the biology of human pigmentation and opens the doorway to attractive advances and applications.

In this work, we present a parallel approach to complete and restricted active space second‐order perturbation theory, (CASPT2/RASPT2). We also make an assessment of the performance characteristics of its particular implementation in the Molcas quantum chemistry programming package. Parallel scaling is limited by memory and I/O bandwidth instead of available cores. Significant time savings for calculations on large and complex systems can be achieved by increasing the number of processes on a single machine, as long as memory bandwidth allows, or by using multiple nodes with a fast, low‐latency interconnect. We found that parallel efficiency drops below 50% when using 8–16 cores on the shared‐memory architecture, or 16–32 nodes on the distributed‐memory architecture, depending on the calculation. This limits the scalability of the implementation to a moderate amount of processes. Nonetheless, calculations that took more than 3 days on a serial machine could be performed in less than 5 h on an InfiniBand cluster, where the individual nodes were not even capable of running the calculation because of memory and I/O requirements. This ensures the continuing study of larger molecular systems by means of CASPT2/RASPT2 through the use of the aggregated computational resources offered by distributed computing systems. © 2013 Wiley Periodicals, Inc.

The well-known natural orbitals are defined as eigenfunctions of a one-particle reduced density operator, and can be obtained from a computed density matrix by diagonalization. Similarly, in this article we define the binatural orbitals, which are obtained for a pair of wave functions by a singular value decomposition of a reduced transition density matrix. The pair of states would usually be eigenstates of the electronic Hamiltonian, and the binatural orbitals then serve as a useful tool for the analysis of the transition between these states. More generally, application to any two state functions gives important information as to how the two states differ. Some examples are shown.

Bjorn Roos is one of the pioneers in the development and usage of multiconfigurational methods, in particular, the complete active space self-consistent field method and the perturbational complete active space perturbation theory through second order. To perform multiconfigurational calculations using these methods, a set of active orbitals must be selected, and the success of the methods depends on the choice of this set. This is not only sometimes easy but also sometimes difficult, especially for use of the more recent RASSCF and RASPT2 methods (which use a "restricted active space" rather than the complete one). Although an automated procedure for selecting the active orbitals would be a preferable solution, this does not seem feasible yet. An account of the problem is given, with examples and some approaches that usually work. (C) 2011 Wiley Periodicals, Inc. Int J Quantum Chem 111: 3329-3338, 2011

Laser evaporation of carbon-rich uranium/carbon alloys followed by atom reactions in a solid argon matrix and trapping at 8 K gives weak infrared absorptions for CUO at 852 and 804 cm(-1). A new band at 827 cm(-1) becomes a doublet with mixed carbon 12 and 13 isotopes and exhibits the 1.0381 isotopic frequency ratio, which is appropriate for the UC diatomic molecule, and another new band at 891 cm(-1) gives a three-band mixed isotopic spectrum with the 1.0366 isotopic frequency ratio, which is characteristic of the linear CUC molecule. CASPT2 calculations with dynamical correlation find the C[triple bond]U[triple bond]C ground state as linear 3Sigma(u)+ with 1.840 A bond length and molecular orbital occupancies for an effective bond order of 2.83. Similar calculations with spin-orbit coupling show that the U[triple bond]C diatomic molecule has a quintet (Lambda = 5, Omega = 3) ground state, a similar 1.855 A bond length, and a fully developed triple bond of 2.82 effective bond order.

On the right are shown results of parallelization over nodes up to 24 nodes. The tests were made on LUNARC cluster with relatively slow Gigabit connection. Numerical gradient shows a perfect scaling until the number of nodes is bigger than the number of displacements. Computing and usage of integrals needs small data transfer and thus shows an almost perfect scaling. The result of RASPT2 for a polyyne system is also shown, to be representative of small molecules with a big active space. The scaling is clearly better for this system since only the many-particle density matrices generation in the RASPT2/CASPT2 code is parallelized yet. Leading CPU manufacturers claim that in a near future, CPUs with 96 and more cores will be available. Quantum chemical codes should be adapted in order to use efficiently these technologies.

Some of the new unique features of the MOLCAS quantum chemistry package version 7 are presented in this report. In particular, the Cholesky decomposition method applied to some quantum chemical methods is described. This approach is used both in the context of a straight forward approximation of the two‐electron integrals and in the generation of so‐called auxiliary basis sets. The article describes how the method is implemented for most known wave functions models: self‐consistent field, density functional theory, 2nd order perturbation theory, complete‐active space self‐consistent field multiconfigurational reference 2nd order perturbation theory, and coupled‐cluster methods. The report further elaborates on the implementation of a restricted‐active space self‐consistent field reference function in conjunction with 2nd order perturbation theory. The average atomic natural orbital basis for relativistic calculations, covering the whole periodic table, are described and associated unique properties are demonstrated. Furthermore, the use of the arbitrary order Douglas‐Kroll‐Hess transformation for one‐component relativistic calculations and its implementation are discussed. This section especially focuses on the implementation of the so‐called picture‐change‐free atomic orbital property integrals. Moreover, the ElectroStatic Potential Fitted scheme, a version of a quantum mechanics/molecular mechanics hybrid method implemented in MOLCAS, is described and discussed. Finally, the report discusses the use of the MOLCAS package for advanced studies of photo chemical phenomena and the usefulness of the algorithms for constrained geometry optimization in MOLCAS in association with such studies. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010

Some of the new unique features of the MOLCAS quantum chemistry package version 7 are presented inthis report. In particular, the Cholesky decomposition method applied to some quantum chemical methods is described. This approach is used both in the context of a straight forward approximation of the two-electron integrals and in the generation of so-called auxiliary basis sets. The article describes how the method is implemented for most known wave functions models: self-consistent field, density functional theory, 2nd order perturbation theory, complete-active space self-consistent field multiconfigurational reference 2nd order perturbation theory, and coupled-cluster methods. The report further elaborates on the implementation of a restricted-active space self-consistent field reference function in conjunction with 2nd order perturbation theory. The average atomic natural orbital basis for relativistic calculations, covering the whole periodic table, are described and associated unique properties are demonstrated. Furthermore, the use of the arbitrary order Douglas-Kroll-Hess transformation for one-component relativistic calculations and its implementation are discussed. This section especially focuses on the implementation of the so-called picture-change-free atomic orbital property integrals. Moreover, the ElectroStatic Potential Fitted scheme, a version of a quantum mechanics/molecular mechanics hybrid method implemented in MOLCAS, is described and discussed. Finally, the report discusses the use of the MOLCAS package for advanced studies of photo chemical phenomena and the usefulness of the algorithms for constrained geometry optimization in MOLCAS in association with such studies. (c) 2009 Wiley Periodicals, Inc. J Corn put Chem 31: 224-247, 2010

Laser-ablated group 6 metal atoms react with NF3 and PF3 to form the simple lowest energy N[triple bond]MF3 and P[triple bond]MX3 products, and this investigation has been extended to AsF3. Mo and W atoms react with AsF3 upon excitation by laser ablation or UV irradiation to form stable trigonal As[triple bond]MF3 terminal arsenides. These molecules are identified by comparison of the closely related infrared spectra of the analogous phosphide species and with frequencies calculated by density functional theory and multiconfigurational second order perturbation theory (CASSCF/CASPT2). Computed CASSCF/CASPT2 triple bond lengths for the As[triple bond]MoF3 and As[triple bond]WF3 molecules are 2.240 A and 2.250 A, respectively. The natural bond orders calculated by CASSCF/CASPT2 decrease from 2.67 to 2.60 for P[triple bond]MoF3 to As[triple bond]MoF3 and from 2.74 to 2.70 for P[triple bond]WF3 to As[triple bond]WF3 as the arsenic valence orbitals are less effective than those of phosphorus in bonding to each metal atom and the larger metal orbital size becomes more compatible with the arsenic valence orbitals. The Cr atom reaction gives the arsinidene AsF=CrF2 product instead of the higher energy As[triple bond]CrF3 molecule as the Cr (VI) state is not supported by the softer pnictides.