Theoretical S1 ← S0 Absorption Energies of the Anionic Forms of Oxyluciferin by Variational Monte Carlo and Many-Body Green’s Function Theory

Time-Resolved Excited-State Analysis of Molecular Electron Dynamics by TDDFT and Bethe-Salpeter Equation Formalisms

In this work, a theoretical and computational set of tools to study and analyze time-resolved electron dynamics in molecules, under the influence of one or more external pulses, is presented. By coupling electronic-structure methods with the resolution of the time-dependent Schrödinger equation, we developed and implemented the time-resolved induced density of the electronic wavepacket, the time-resolved formulation of the differential projection density of states (ΔPDOS), and of transition contribution map (TCM) to look at the single-electron orbital occupation and localization change in time. Moreover, to further quantify the possible charge transfer, we also defined the energy-integrated ΔPDOS and the fragment-projected TCM. We have used time-dependent density-functional theory (TDDFT), as implemented in ADF software, and the Bethe-Salpeter equation, as provided by MolGW package, for the description of the electronic excited states. This suite of postprocessing tools also provides the time evolution of the electronic states of the system of interest. To illustrate the usefulness of these postprocessing tools, excited-state populations have been computed for HBDI (the chromophore of GFP) and DNQDI molecules interacting with a sequence of two pulses. Time-resolved descriptors have been applied to study the time-resolved electron dynamics of HBDI, DNQDI, LiCN (being a model system for dipole switching upon highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) electronic excitation), and Ag22. The computational analysis tools presented in this article can be employed to help the interpretation of fast and ultrafast spectroscopies on molecular, supramolecular, and composite systems.

H2O-Bridged Proton-Transfer Channel in Emitter Species Formation in Obelin Bioluminescence

Bioluminescence of a number of marine organisms is conditioned by Ca²⁺-regulated photoprotein (CaRP) with coelenterazine as the reaction substrate. The reaction product, coelenteramide, at the first singlet excited state (S₁) is the emitter of CaRP. The S₁-state coelenteramide is produced via the decomposition of coelenterazine dioxetanone. Experiments suggested that the neutral S₁-coelenteramide is the primary emitter species. This supposition contradicts with theoretical calculations showing that the anionic S₁-coelenteramide is a primary product of the decomposition of coelenterazine dioxetanone. In this study, applying molecular dynamic (MD) simulations and the hybrid quantum mechanics/molecular mechanics (QM/MM) method, we investigated a proton-transfer (PT) process taking place in CaRP obelin from Obelia longissima for emitter formation. Our calculations demonstrate a concerted PT process with a water molecule as a bridge between anionic S₁-coelenteramide and the nearest histidine residue. The low activation barrier as well as the strong hydrogen-bond network between the proton donor and the proton acceptor suggests a fast PT process comparable with that of the lifetime of excited anionic S₁-coelenteramide. The existence of the PT process eliminates the discrepancy between experimental and theoretical studies. The fast PT process at emitter formation can also take place in other CaRPs.

Quantum-mechanical hydration plays critical role in the stability of firefly oxyluciferin isomers: State-of-the-art calculations of the excited states

Stabilizing mechanisms of three possible isomers (phenolate-keto, phenolate-enol, and phenol-enolate) of the oxyluciferin anion hydrated with quantum explicit water molecules in the first singlet excited state were investigated using first-principles Born-Oppenheimer molecular dynamics simulations for up to 1.8 ns (or 3.7 × 10⁶ MD steps), revealing that the surrounding water molecules were distributed to form clear single-layered structures for phenolate-keto and multi-layered structures for phenolate-enol and phenol-enolate isomers. The isomers employed different stabilizing mechanisms compared to the ground state. Only the phenolate-keto isomer became attracted to the water molecules in its excited state and was stabilized by increasing the number of hydrogen bonds with nearby water molecules. The most stable isomer in the excited state was the phenolate-keto, and the phenolate-enol and phenol-enolate isomers were higher in energy by ∼0.38 eV and 0.57 eV, respectively, than the phenolate-keto. This was in contrast to the case of ground state in which the phenolate-enol was the most stable isomer.

Many-body perturbation theory calculations using the yambo code

yambo is an open source project aimed at studying excited state properties of condensed matter systems from first principles using many-body methods. As input, yambo requires ground state electronic structure data as computed by density functional theory codes such as Quantum ESPRESSO and Abinit. yambo's capabilities include the calculation of linear response quantities (both independent-particle and including electron-hole interactions), quasi-particle corrections based on the GW formalism, optical absorption, and other spectroscopic quantities. Here we describe recent developments ranging from the inclusion of important but oft-neglected physical effects such as electron-phonon interactions to the implementation of a real-time propagation scheme for simulating linear and non-linear optical properties. Improvements to numerical algorithms and the user interface are outlined. Particular emphasis is given to the new and efficient parallel structure that makes it possible to exploit modern high performance computing architectures. Finally, we demonstrate the possibility to automate workflows by interfacing with the yambopy and AiiDA software tools.

QMCPACK: an open source ab initio quantum Monte Carlo package for the electronic structure of atoms, molecules and solids

QMCPACK is an open source quantum Monte Carlo package for ab initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wavefunctions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary-field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit and graphical processing unit systems. We detail the program's capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://qmcpack.org.

The Bethe–Salpeter equation in chemistry: relations with TD-DFT, applications and challenges

We review the Bethe–Salpeter formalism and analyze its performances for the calculation of the excited state properties of molecular systems.