RESUMEN
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.
RESUMEN
Ultrafast non-adiabatic dynamics of the small heteroaromatic compound 2-thiopyridone has been studied with surface hopping simulations based on multi-configurational quantum chemistry. Initial excitation of the bright S 2 ( π , π * ) state is found to promptly relax to S 1 ( n , π * ) through in-plane motion. The subsequent dynamics are oppositely driven by out-of-plane motion, which results in both complex population transfers among all the available states and intersystem crossing predominantly through the "El-Sayed forbidden" S 1 ( n , π * ) to T 2 ( n , π * ) channel, through significant mixing of electronic excitation characters. Despite this complexity, the femto- to picosecond triplet population, expected from several spectroscopic measurements, is well described as a simple exponential decay of the singlet state manifold. No proton transfer is found in the reported trajectories, but two mechanisms for its possible mediation in previously reported experiments are proposed based on the observed structural dynamics: (i) ultrafast intra-molecular transfer driven by the initially coherent in-plane motion and (ii) inter-molecular solvent-mediated transfer driven by the out-of-plane modes that dominate the later motion.
RESUMEN
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.
RESUMEN
2-Thiopyridone (2-TP), a common model system for excited-state proton transfer, has been simulated in aqueous solution with ab initio molecular dynamics. The interplay of electronic structure, protonation, and solvatization is investigated by comparison of three differently protonated molecular forms and between the lowest singlet and triplet electronic states. An interdependence clearly manifests in the mixed-character T1 state for the 2-TP form, systematic structural distortions of the 2-mercaptopyridine (2-MP) form, and photobase protolysis of the 2-TP- form, in the aqueous phase. In comparison, simplified continuum models for the solvatization are found to be significantly inaccurate for several of the species. To facilitate future computational studies, we therefore present a minimal representative solvatization complex for each stable form and electronic state. Our findings demonstrate the importance of explicit solvatization of the compound and sets the stage for including it also in future studies.
RESUMEN
Excited-state proton transfer (ESPT) is a fundamental process in biomolecular photochemistry, but its underlying mediators often evade direct observation. We identify a distinct pathway for ESPT in aqueous 2-thiopyridone, by employing transient N 1s X-ray absorption spectroscopy and multi-configurational spectrum simulations. Photoexcitations to the singlet S2 and S4 states both relax promptly through intersystem crossing to the triplet T1 state. The T1 state, through its rapid population and near nanosecond lifetime, mediates nitrogen site deprotonation by ESPT in a secondary intersystem crossing to the S0 potential energy surface. This conclusively establishes a dominant ESPT pathway for the system in aqueous solution, which is also compatible with previous measurements in acetonitrile. Thereby, the hitherto open questions of the pathway for ESPT in the compound, including its possible dependence on excitation wavelength and choice of solvent, are resolved.
RESUMEN
Understanding and controlling properties of transition metal complexes is a crucial step towards tailoring materials for sustainable energy applications. In a systematic approach, we use resonant inelastic X-ray scattering to study the influence of ligand substitution on the valence electronic structure around an aqueous iron(ii) center. Exchanging cyanide with 2-2'-bipyridine ligands reshapes frontier orbitals in a way that reduces metal 3d charge delocalization onto the ligands. This net decrease of metal-ligand covalency results in lower metal-centered excited state energies in agreement with previously reported excited state dynamics. Furthermore, traces of solvent-effects were found indicating a varying interaction strength of the solvent with ligands of different character. Our results demonstrate how ligand exchange can be exploited to shape frontier orbitals of transition metal complexes in solution-phase chemistry; insights upon which future efforts can built when tailoring the functionality of photoactive systems for light-harvesting applications.
RESUMEN
A combination of multi-configurational restricted active space calculations with a Dyson orbital formalism has been applied for accurate simulations of 4d photo-electron spectra of the I3- molecular ion. The analysis based on the occupation numbers of natural orbitals allowed to predict and rationalize the spectral fingerprints of solvent-induced nuclear asymmetry. In particular, it demonstrates how the nuclear asymmetry directly causes an increase of shake-up intensity. The relative intensity of shake-up and main features of the I 4d XPS spectrum could therefore serve as a simplified experimental observable of structural asymmetry, complementary to changes in the shape of the main spectral features.
RESUMEN
Soft X-ray spectroscopies are ideal probes of the local valence electronic structure of photocatalytically active metal sites. Here, we apply the selectivity of time-resolved resonant inelastic X-ray scattering at the iron L-edge to the transient charge distribution of an optically excited charge-transfer state in aqueous ferricyanide. Through comparison to steady-state spectra and quantum chemical calculations, the coupled effects of valence-shell closing and ligand-hole creation are experimentally and theoretically disentangled and described in terms of orbital occupancy, metal-ligand covalency, and ligand field splitting, thereby extending established steady-state concepts to the excited-state domain. π-Back-donation is found to be mainly determined by the metal site occupation, whereas the ligand hole instead influences σ-donation. Our results demonstrate how ultrafast resonant inelastic X-ray scattering can help characterize local charge distributions around catalytic metal centers in short-lived charge-transfer excited states, as a step toward future rationalization and tailoring of photocatalytic capabilities of transition-metal complexes.
RESUMEN
We describe how inversion symmetry separation of electronic state manifolds in resonant inelastic soft X-ray scattering (RIXS) can be applied to probe excited-state dynamics with compelling selectivity. In a case study of Fe L3-edge RIXS in the ferricyanide complex Fe(CN)63-, we demonstrate with multi-configurational restricted active space spectrum simulations how the information content of RIXS spectral fingerprints can be used to unambiguously separate species of different electronic configurations, spin multiplicities, and structures, with possible involvement in the decay dynamics of photo-excited ligand-to-metal charge-transfer. Specifically, we propose that this could be applied to confirm or reject the presence of a hitherto elusive transient Quartet species. Thus, RIXS offers a particular possibility to settle a recent controversy regarding the decay pathway, and we expect the technique to be similarly applicable in other model systems of photo-induced dynamics.
RESUMEN
The femtosecond excited-state dynamics following resonant photoexcitation enable the selective deformation of N-H and N-C chemical bonds in 2-thiopyridone in aqueous solution with optical or X-ray pulses. In combination with multiconfigurational quantum-chemical calculations, the orbital-specific electronic structure and its ultrafast dynamics accessed with resonant inelastic X-ray scattering at the Nâ 1s level using synchrotron radiation and the soft X-ray free-electron laser LCLS provide direct evidence for this controlled photoinduced molecular deformation and its ultrashort timescale.
RESUMEN
Simple models for friction are typically one-dimensional, but real interfaces are two-dimensional. We investigate the effects of the second dimension on static and dynamic friction by using the Frenkel-Kontorova (FK) model. We study the two most straightforward extensions of the FK model to two dimensions and simulate both the static and dynamic properties. We show that the behavior of the static friction is robust and remains similar in two dimensions for physically reasonable parameter values. The dynamic friction, however, is strongly influenced by the second dimension and the accompanying additional dynamics and parameters introduced into the models. We discuss our results in terms of the thermal equilibration and phonon dispersion relations of the lattices, establishing a physically realistic and suitable two-dimensional extension of the FK model. We find that the presence of additional dissipation channels can increase the friction and produces significantly different temperature dependence when compared to the one-dimensional case. We also briefly study the anisotropy of the dynamic friction and show highly nontrivial effects, including that the friction anisotropy can lead to motion in different directions depending on the value of the initial velocity.
