1
|
Pieri E, Walker AR, Zhu M, Martínez TJ. Conical Intersection Accessibility Dictates Brightness in Red Fluorescent Proteins. J Am Chem Soc 2024; 146:17646-17658. [PMID: 38885641 DOI: 10.1021/jacs.4c00458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/20/2024]
Abstract
Red fluorescent protein (RFP) variants are highly sought after for in vivo imaging since longer wavelengths improve depth and contrast in fluorescence imaging. However, the lower energy emission wavelength usually correlates with a lower fluorescent quantum yield compared to their green emitting counterparts. To guide the rational design of bright variants, we have theoretically assessed two variants (mScarlet and mRouge) which are reported to have very different brightness. Using an α-CASSCF QM/MM framework (chromophore and all protein residues within 6 Å of it in the QM region, for a total of more than 450 QM atoms), we identify key points on the ground and first excited state potential energy surfaces. The brighter variant mScarlet has a rigid scaffold, and the chromophore stays largely planar on the ground state. The dimmer variant mRouge shows more flexibility and can accommodate a pretwisted chromophore conformation which provides easier access to conical intersections. The main difference between the variants lies in the intersection seam regions, which appear largely inaccessible in mScarlet but partially accessible in mRouge. This observation is mainly related with changes in the cavity charge distribution, the hydrogen-bonding network involving the chromophore and a key ARG/THR mutation (which changes both charge and steric hindrance).
Collapse
Affiliation(s)
- Elisa Pieri
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599, United States
| | - Alice R Walker
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - Mingning Zhu
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
2
|
Shu Y, Truhlar DG. Generalized Semiclassical Ehrenfest Method: A Route to Wave Function-Free Photochemistry and Nonadiabatic Dynamics with Only Potential Energies and Gradients. J Chem Theory Comput 2024; 20:4396-4426. [PMID: 38819014 DOI: 10.1021/acs.jctc.4c00424] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
3
|
Kim SS, Rhee YM. Potential energy interpolation with target-customized weighting coordinates: application to excited-state dynamics of photoactive yellow protein chromophore in water. Phys Chem Chem Phys 2024; 26:9021-9036. [PMID: 38440829 DOI: 10.1039/d3cp05643k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2024]
Abstract
Interpolation of potential energy surfaces (PESs) can provide a practical route to performing molecular dynamics simulations with a reliability matching a high-level quantum chemical calculation. An obstacle to its widespread use is perhaps the lack of general and optimal interpolation settings that can be applied in a black-box manner for any given molecular system. How to set up the weights for interpolation is one such task, and we still need to diversify the approaches in order to treat various systems. Here, we develop a new interpolation weighting scheme, which allows us to choose the weighting coordinates in a system-specific manner, by amplifying the contribution from specific internal coordinates. The new weighting scheme with an appropriate selection of coordinates is proved to be effective in reducing the interpolation error along the reaction pathway. As a demonstration, we consider the photoactive yellow protein chromophore system, as it constitutes itself as an interesting target that bears long-standing questions related to excited-state dynamics inside protein environments. We build its two-state diabatic interpolated PES with the new weighting scheme. We indeed see the utility of our scheme by conducting nonadiabatic molecular dynamics simulations with the required semi-global PES based on a limited number of data points.
Collapse
Affiliation(s)
- Seung Soo Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| |
Collapse
|
4
|
List NH, Jones CM, Martínez TJ. Chemical control of excited-state reactivity of the anionic green fluorescent protein chromophore. Commun Chem 2024; 7:25. [PMID: 38316834 PMCID: PMC10844232 DOI: 10.1038/s42004-024-01099-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2023] [Accepted: 01/05/2024] [Indexed: 02/07/2024] Open
Abstract
Controlling excited-state reactivity is a long-standing challenge in photochemistry, as a desired pathway may be inaccessible or compete with other unwanted channels. An important example is internal conversion of the anionic green fluorescent protein (GFP) chromophore where non-selective progress along two competing torsional modes (P: phenolate and I: imidazolinone) impairs and enables Z-to-E photoisomerization, respectively. Developing strategies to promote photoisomerization could drive new areas of applications of GFP-like proteins. Motivated by the charge-transfer dichotomy of the torsional modes, we explore chemical substitution on the P-ring of the chromophore as a way to control excited-state pathways and improve photoisomerization. As demonstrated by methoxylation, selective P-twisting appears difficult to achieve because the electron-donating potential effects of the substituents are counteracted by inertial effects that directly retard the motion. Conversely, these effects act in concert to promote I-twisting when introducing electron-withdrawing groups. Specifically, 2,3,5-trifluorination leads to both pathway selectivity and a more direct approach to the I-twisted intersection which, in turn, doubles the photoisomerization quantum yield. Our results suggest P-ring engineering as an effective approach to boost photoisomerization of the anionic GFP chromophore.
Collapse
Affiliation(s)
- Nanna H List
- Department of Chemistry, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-10044, Stockholm, Sweden.
| | - Chey M Jones
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, CA, 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, CA, 94305, USA.
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA.
| |
Collapse
|
5
|
Liu Y, Sanchez DM, Ware MR, Champenois EG, Yang J, Nunes JPF, Attar A, Centurion M, Cryan JP, Forbes R, Hegazy K, Hoffmann MC, Ji F, Lin MF, Luo D, Saha SK, Shen X, Wang XJ, Martínez TJ, Wolf TJA. Rehybridization dynamics into the pericyclic minimum of an electrocyclic reaction imaged in real-time. Nat Commun 2023; 14:2795. [PMID: 37202402 DOI: 10.1038/s41467-023-38513-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2022] [Accepted: 04/28/2023] [Indexed: 05/20/2023] Open
Abstract
Electrocyclic reactions are characterized by the concerted formation and cleavage of both σ and π bonds through a cyclic structure. This structure is known as a pericyclic transition state for thermal reactions and a pericyclic minimum in the excited state for photochemical reactions. However, the structure of the pericyclic geometry has yet to be observed experimentally. We use a combination of ultrafast electron diffraction and excited state wavepacket simulations to image structural dynamics through the pericyclic minimum of a photochemical electrocyclic ring-opening reaction in the molecule α-terpinene. The structural motion into the pericyclic minimum is dominated by rehybridization of two carbon atoms, which is required for the transformation from two to three conjugated π bonds. The σ bond dissociation largely happens after internal conversion from the pericyclic minimum to the electronic ground state. These findings may be transferrable to electrocyclic reactions in general.
Collapse
Affiliation(s)
- Y Liu
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY, 11790, USA
| | - D M Sanchez
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
- Department of Chemistry, Stanford University, 333 Campus Drive, Stanford, CA, 94305, USA
- Design Physics Division, Lawrence Livermore National Laboratory, Livermore, CA, USA
| | - M R Ware
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - E G Champenois
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - J Yang
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
- Center of Basic Molecular Science, Department of Chemistry, Mong Man Wai Building of Science and Technology, S-1027 Tsinghua University, Beijing, China
| | - J P F Nunes
- Department of Physics and Astronomy, University of Nebraska-Lincoln, Theodore Jorgensen Hall 208, 855 N 16th Street, Lincoln, NE, 68588, USA
- Diamond Light Source, Harwell Science Campus, Fermi Ave, Didcot, OX11 0DE, UK
| | - A Attar
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - M Centurion
- Department of Physics and Astronomy, University of Nebraska-Lincoln, Theodore Jorgensen Hall 208, 855 N 16th Street, Lincoln, NE, 68588, USA
| | - J P Cryan
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - R Forbes
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - K Hegazy
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - M C Hoffmann
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - F Ji
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - M-F Lin
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - D Luo
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - S K Saha
- Department of Physics and Astronomy, University of Nebraska-Lincoln, Theodore Jorgensen Hall 208, 855 N 16th Street, Lincoln, NE, 68588, USA
| | - X Shen
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - X J Wang
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - T J Martínez
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA.
- Department of Chemistry, Stanford University, 333 Campus Drive, Stanford, CA, 94305, USA.
| | - T J A Wolf
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA.
| |
Collapse
|
6
|
Substitution Effects on the Mechanism of Light-Induced 2,5-Diaryltetrazole-Naphthoquinone 1,3-Dipolar Cycloaddition: A Theoretical Study. COMPUT THEOR CHEM 2023. [DOI: 10.1016/j.comptc.2023.114060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
|
7
|
Jones CM, List NH, Martínez TJ. Steric and Electronic Origins of Fluorescence in GFP and GFP-like Proteins. J Am Chem Soc 2022; 144:12732-12746. [PMID: 35786916 DOI: 10.1021/jacs.2c02946] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Fluorescent proteins have become routine tools for biological imaging. However, their nanosecond lifetimes on the excited state present computational hurdles to a full understanding of these photoactive proteins. In this work, we simulate approximately 0.5 nanoseconds of ab initio molecular dynamics to elucidate steric and electronic features responsible for fluorescent protein behavior. Using green fluorescent protein (GFP) and Dronpa2─widely used fluorescent proteins with contrasting functionality─as case studies, we leverage previous findings in the gas phase and solution to explore the deactivation mechanisms available to these proteins. Starting with ground-state analyses, we identify steric (the distribution of empty pockets near the chromophore) and electronic (electric fields exerted on chromophore moieties) factors that offer potential avenues for rational design. Picosecond timescale simulations on the excited state reveal that the chromophore can access twisted structures in Dronpa2, while the chromophore is largely confined to planarity in GFP. We couple ab initio multiple spawning (AIMS) and enhanced sampling simulations to discover and characterize conical intersection seams that facilitate internal conversion, which is a rare event in both systems. Our AIMS simulations correctly capture the relative fluorescence profiles of GFP and Dronpa2 within the first few picoseconds, and we attribute the diminished fluorescence intensity of Dronpa2, relative to GFP, to flexible chromophore intermediates on the excited state. Furthermore, we predict that twisted chromophore intermediates produce red-shifted intensities in the Dronpa2 fluorescence spectrum. If confirmed experimentally, this spectroscopic signature would provide valuable insights when screening and developing novel fluorescent proteins.
Collapse
Affiliation(s)
- Chey M Jones
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Nanna H List
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
8
|
Raucci U, Weir H, Bannwarth C, Sanchez DM, Martínez TJ. Chiral photochemistry of achiral molecules. Nat Commun 2022; 13:2091. [PMID: 35440559 PMCID: PMC9019051 DOI: 10.1038/s41467-022-29662-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2021] [Accepted: 03/25/2022] [Indexed: 11/23/2022] Open
Abstract
Chirality is a molecular property governed by the topography of the potential energy surface (PES). Thermally achiral molecules interconvert rapidly when the interconversion barrier between the two enantiomers is comparable to or lower than the thermal energy, in contrast to thermally stable chiral configurations. In principle, a change in the PES topography on the excited electronic state may diminish interconversion, leading to electronically prochiral molecules that can be converted from achiral to chiral by electronic excitation. Here we report that this is the case for two prototypical examples – cis-stilbene and cis-stiff stilbene. Both systems exhibit unidirectional photoisomerization for each enantiomer as a result of their electronic prochirality. We simulate an experiment to demonstrate this effect in cis-stilbene based on its interaction with circularly polarized light. Our results highlight the drastic change in chiral behavior upon electronic excitation, opening up the possibility for asymmetric photochemistry from an effectively nonchiral starting point. The authors report non-adiabatic first principles molecular dynamics to show how an achiral molecule can be converted to a chiral one upon photoexcitation. These results demonstrate the possibility of asymmetric photochemistry starting from achiral reactants.
Collapse
Affiliation(s)
- Umberto Raucci
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, CA, 94305, USA.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA.,Italian Institute of Technology, Genova, GE, Italy
| | - Hayley Weir
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, CA, 94305, USA.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - Christoph Bannwarth
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, CA, 94305, USA.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA.,Institute of Physical Chemistry, RWTH Aachen University, Aachen, Germany
| | - David M Sanchez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, CA, 94305, USA.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA.,Design Physics Division, Lawrence Livermore National Laboratory, Livermore, CA, USA
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, CA, 94305, USA. .,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA.
| |
Collapse
|
9
|
Song C. State-averaged CASSCF with polarizable continuum model for studying photoreactions in solvents: Energies, analytical nuclear gradients, and non-adiabatic couplings. J Chem Phys 2022; 156:104102. [DOI: 10.1063/5.0085855] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
This paper presents state-averaged complete active space self-consistent field in polarizable continuum model (PCM) for studies of photoreactions in solvents. The wavefunctions of the solute and the PCM surface charges of the solvent are optimized simultaneously such that the state-averaged free energy is variationally minimized. The method supports both fixed weights and dynamic weights where the weights are automatically adjusted based on the energy gaps. The corresponding analytical nuclear gradients and non-adiabatic couplings are also derived. Furthermore, we show how the new method can be entirely formulated in terms of seven basic operations, which allows the implementation to benefit from existing high-performance libraries on graphical processing units. Results demonstrating the accuracy and performance of the implementation are presented and discussed. We also apply the new method to the study of minimal conical intersection search and photoreaction energy pathways in solvents. Effects from the polarity of the solvents and different formulas of dynamic weights are compared and discussed.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Chemistry, University of California Davis, Davis, California 95616, USA
| |
Collapse
|
10
|
Champenois EG, Sanchez DM, Yang J, Figueira Nunes JP, Attar A, Centurion M, Forbes R, Gühr M, Hegazy K, Ji F, Saha SK, Liu Y, Lin MF, Luo D, Moore B, Shen X, Ware MR, Wang XJ, Martínez TJ, Wolf TJA. Conformer-specific photochemistry imaged in real space and time. Science 2021; 374:178-182. [PMID: 34618569 DOI: 10.1126/science.abk3132] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Abstract
[Figure: see text].
Collapse
Affiliation(s)
- E G Champenois
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - D M Sanchez
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA.,Department of Chemistry, Stanford University, Stanford, CA, USA
| | - J Yang
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA.,SLAC National Accelerator Laboratory, Menlo Park, CA, USA.,Center of Basic Molecular Science, Department of Chemistry, Tsinghua University, Beijing, China
| | - J P Figueira Nunes
- Department of Physics and Astronomy, University of Nebraska, Lincoln, NE, USA
| | - A Attar
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - M Centurion
- Department of Physics and Astronomy, University of Nebraska, Lincoln, NE, USA
| | - R Forbes
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - M Gühr
- Institut für Physik und Astronomie, Universität Potsdam, Potsdam, Germany
| | - K Hegazy
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA.,Department of Physics, Stanford University, Stanford, CA, USA
| | - F Ji
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - S K Saha
- Department of Physics and Astronomy, University of Nebraska, Lincoln, NE, USA
| | - Y Liu
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY, USA
| | - M-F Lin
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - D Luo
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - B Moore
- Department of Physics and Astronomy, University of Nebraska, Lincoln, NE, USA
| | - X Shen
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - M R Ware
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - X J Wang
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - T J Martínez
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA.,Department of Chemistry, Stanford University, Stanford, CA, USA
| | - T J A Wolf
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| |
Collapse
|
11
|
Jones CM, List NH, Martínez TJ. Resolving the ultrafast dynamics of the anionic green fluorescent protein chromophore in water. Chem Sci 2021; 12:11347-11363. [PMID: 34667545 PMCID: PMC8447926 DOI: 10.1039/d1sc02508b] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Accepted: 07/13/2021] [Indexed: 11/21/2022] Open
Abstract
The chromophore of the green fluorescent protein (GFP) is critical for probing environmental influences on fluorescent protein behavior. Using the aqueous system as a bridge between the unconfined vacuum system and a constricting protein scaffold, we investigate the steric and electronic effects of the environment on the photodynamical behavior of the chromophore. Specifically, we apply ab initio multiple spawning to simulate five picoseconds of nonadiabatic dynamics after photoexcitation, resolving the excited-state pathways responsible for internal conversion in the aqueous chromophore. We identify an ultrafast pathway that proceeds through a short-lived (sub-picosecond) imidazolinone-twisted (I-twisted) species and a slower (several picoseconds) channel that proceeds through a long-lived phenolate-twisted (P-twisted) intermediate. The molecule navigates the non-equilibrium energy landscape via an aborted hula-twist-like motion toward the one-bond-flip dominated conical intersection seams, as opposed to following the pure one-bond-flip paths proposed by the excited-state equilibrium picture. We interpret our simulations in the context of time-resolved fluorescence experiments, which use short- and long-time components to describe the fluorescence decay of the aqueous GFP chromophore. Our results suggest that the longer time component is caused by an energetically uphill approach to the P-twisted intersection seam rather than an excited-state barrier to reach the twisted intramolecular charge-transfer species. Irrespective of the location of the nonadiabatic population events, the twisted intersection seams are inefficient at facilitating isomerization in aqueous solution. The disordered and homogeneous nature of the aqueous solvent environment facilitates non-selective stabilization with respect to I- and P-twisted species, offering an important foundation for understanding the consequences of selective stabilization in heterogeneous and rigid protein environments.
Collapse
Affiliation(s)
- Chey M Jones
- Department of Chemistry and the PULSE Institute, Stanford University Stanford CA 94305 USA
- SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Nanna H List
- Department of Chemistry and the PULSE Institute, Stanford University Stanford CA 94305 USA
- SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University Stanford CA 94305 USA
- SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| |
Collapse
|
12
|
Abstract
In this article, we review nonadiabatic molecular dynamics (NAMD) methods for modeling spin-crossover transitions. First, we discuss different representations of electronic states employed in the grid-based and direct NAMD simulations. The nature of interstate couplings in different representations is highlighted, with the main focus on nonadiabatic and spin-orbit couplings. Second, we describe three NAMD methods that have been used to simulate spin-crossover dynamics, including trajectory surface hopping, ab initio multiple spawning, and multiconfiguration time-dependent Hartree. Some aspects of employing different electronic structure methods to obtain information about potential energy surfaces and interstate couplings for NAMD simulations are also discussed. Third, representative applications of NAMD to spin crossovers in molecular systems of different sizes and complexities are highlighted. Finally, we pose several fundamental questions related to spin-dependent processes. These questions should be possible to address with future methodological developments in NAMD.
Collapse
Affiliation(s)
- Saikat Mukherjee
- Institut de Chimie Radicalaire, CNRS 7273, Aix-Marseille University, 13013 Marseille, France;
| | - Dmitry A Fedorov
- Oak Ridge Associated Universities, Oak Ridge, Tennessee 37830, USA;
| | - Sergey A Varganov
- Department of Chemistry, University of Nevada, Reno, Nevada 89557-0216, USA;
| |
Collapse
|
13
|
Li X, Yu H, Lou F, Feng J, Whangbo MH, Xiang H. Spin Hamiltonians in Magnets: Theories and Computations. Molecules 2021; 26:803. [PMID: 33557181 PMCID: PMC7913993 DOI: 10.3390/molecules26040803] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2020] [Revised: 02/01/2021] [Accepted: 02/02/2021] [Indexed: 11/16/2022] Open
Abstract
The effective spin Hamiltonian method has drawn considerable attention for its power to explain and predict magnetic properties in various intriguing materials. In this review, we summarize different types of interactions between spins (hereafter, spin interactions, for short) that may be used in effective spin Hamiltonians as well as the various methods of computing the interaction parameters. A detailed discussion about the merits and possible pitfalls of each technique of computing interaction parameters is provided.
Collapse
Affiliation(s)
- Xueyang Li
- Key Laboratory of Computational Physical Sciences (Ministry of Education), State Key Laboratory of Surface Physics, Department of Physics, Fudan University, Shanghai 200433, China; (X.L.); (H.Y.); (F.L.); (J.F.)
- Shanghai Qi Zhi Institute, Shanghai 200232, China
| | - Hongyu Yu
- Key Laboratory of Computational Physical Sciences (Ministry of Education), State Key Laboratory of Surface Physics, Department of Physics, Fudan University, Shanghai 200433, China; (X.L.); (H.Y.); (F.L.); (J.F.)
- Shanghai Qi Zhi Institute, Shanghai 200232, China
| | - Feng Lou
- Key Laboratory of Computational Physical Sciences (Ministry of Education), State Key Laboratory of Surface Physics, Department of Physics, Fudan University, Shanghai 200433, China; (X.L.); (H.Y.); (F.L.); (J.F.)
- Shanghai Qi Zhi Institute, Shanghai 200232, China
| | - Junsheng Feng
- Key Laboratory of Computational Physical Sciences (Ministry of Education), State Key Laboratory of Surface Physics, Department of Physics, Fudan University, Shanghai 200433, China; (X.L.); (H.Y.); (F.L.); (J.F.)
- School of Physics and Materials Engineering, Hefei Normal University, Hefei 230601, China
| | - Myung-Hwan Whangbo
- Department of Chemistry, North Carolina State University, Raleigh, NC 27695-8204, USA;
| | - Hongjun Xiang
- Key Laboratory of Computational Physical Sciences (Ministry of Education), State Key Laboratory of Surface Physics, Department of Physics, Fudan University, Shanghai 200433, China; (X.L.); (H.Y.); (F.L.); (J.F.)
- Shanghai Qi Zhi Institute, Shanghai 200232, China
| |
Collapse
|
14
|
Seritan S, Bannwarth C, Fales BS, Hohenstein EG, Isborn CM, Kokkila‐Schumacher SIL, Li X, Liu F, Luehr N, Snyder JW, Song C, Titov AV, Ufimtsev IS, Wang L, Martínez TJ. TeraChem
: A graphical processing unit
‐accelerated
electronic structure package for
large‐scale
ab initio molecular dynamics. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1494] [Citation(s) in RCA: 53] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Stefan Seritan
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Christoph Bannwarth
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Bryan S. Fales
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Edward G. Hohenstein
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Christine M. Isborn
- Department of Chemistry University of California Merced Merced California USA
| | | | - Xin Li
- Division of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health KTH Royal Institute of Technology Stockholm Sweden
| | - Fang Liu
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge Massachusetts USA
| | | | | | - Chenchen Song
- Department of Physics University of California Berkeley Berkeley California USA
- Molecular Foundry Lawrence Berkeley National Laboratory Berkeley California USA
| | | | - Ivan S. Ufimtsev
- Department of Structural Biology Stanford University School of Medicine Stanford California USA
| | - Lee‐Ping Wang
- Department of Chemistry University of California Davis Davis California USA
| | - Todd J. Martínez
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| |
Collapse
|
15
|
Bannwarth C, Yu JK, Hohenstein EG, Martínez TJ. Hole-hole Tamm-Dancoff-approximated density functional theory: A highly efficient electronic structure method incorporating dynamic and static correlation. J Chem Phys 2020; 153:024110. [PMID: 32668944 DOI: 10.1063/5.0003985] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The study of photochemical reaction dynamics requires accurate as well as computationally efficient electronic structure methods for the ground and excited states. While time-dependent density functional theory (TDDFT) is not able to capture static correlation, complete active space self-consistent field methods neglect much of the dynamic correlation. Hence, inexpensive methods that encompass both static and dynamic electron correlation effects are of high interest. Here, we revisit hole-hole Tamm-Dancoff approximated (hh-TDA) density functional theory for this purpose. The hh-TDA method is the hole-hole counterpart to the more established particle-particle TDA (pp-TDA) method, both of which are derived from the particle-particle random phase approximation (pp-RPA). In hh-TDA, the N-electron electronic states are obtained through double annihilations starting from a doubly anionic (N+2 electron) reference state. In this way, hh-TDA treats ground and excited states on equal footing, thus allowing for conical intersections to be correctly described. The treatment of dynamic correlation is introduced through the use of commonly employed density functional approximations to the exchange-correlation potential. We show that hh-TDA is a promising candidate to efficiently treat the photochemistry of organic and biochemical systems that involve several low-lying excited states-particularly those with both low-lying ππ* and nπ* states where inclusion of dynamic correlation is essential to describe the relative energetics. In contrast to the existing literature on pp-TDA and pp-RPA, we employ a functional-dependent choice for the response kernel in pp- and hh-TDA, which closely resembles the response kernels occurring in linear response and collinear spin-flip TDDFT.
Collapse
Affiliation(s)
- Christoph Bannwarth
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Jimmy K Yu
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Edward G Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
16
|
Song C, Martínez TJ. Reduced scaling extended multi-state CASPT2 (XMS-CASPT2) using supporting subspaces and tensor hyper-contraction. J Chem Phys 2020; 152:234113. [PMID: 32571032 DOI: 10.1063/5.0007417] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
We present a reduced scaling formulation of the extended multi-state CASPT2 (XMS-CASPT2) method, which is based on our recently developed state-specific CASPT2 (SS-CASPT2) formulation using supporting subspaces and tensor hyper-contraction. By using these two techniques, the off-diagonal elements of the effective Hamiltonian can be computed with only O(N3) operations and O(N2) memory, where N is the number of basis functions. This limits the overall computational scaling to O(N4) operations and O(N2) memory. Thus, excited states can now be obtained at the same reduced (relative to previous algorithms) scaling we achieved for SS-CASPT2. In addition, we also investigate how the energy denominators can be factorized with the Laplace quadrature when some of the denominators are negative, which is critical for excited state calculations. An efficient implementation of the method has been developed using graphical processing units while also exploiting spatial sparsity in tensor operations. We benchmark the accuracy of the new method by comparison to non-THC formulated XMS-CASPT2 for the excited states of various molecules. In our tests, the THC approximation introduces negligible errors (≈0.01 eV) compared to the non-THC reference method. Scaling behavior and computational timings are presented to demonstrate performance. The new method is also interfaced with quantum mechanics/molecular mechanics (QM/MM). In an example study of green fluorescent protein, we show how the XMS-CASPT2 potential energy surfaces and excitation energies are affected by increasing the size of the QM region up to 278 QM atoms with more than 2300 basis functions.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
17
|
Tran LN, Shea JAR, Neuscamman E. Tracking Excited States in Wave Function Optimization Using Density Matrices and Variational Principles. J Chem Theory Comput 2019; 15:4790-4803. [DOI: 10.1021/acs.jctc.9b00351] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Lan Nguyen Tran
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Ho Chi Minh City Institute of Physics, VAST, Ho Chi Minh City 700000, Vietnam
| | - Jacqueline A. R. Shea
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
18
|
Levine BG, Esch MP, Fales BS, Hardwick DT, Peng WT, Shu Y. Conical Intersections at the Nanoscale: Molecular Ideas for Materials. Annu Rev Phys Chem 2019; 70:21-43. [DOI: 10.1146/annurev-physchem-042018-052425] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The ability to predict and describe nonradiative processes in molecules via the identification and characterization of conical intersections is one of the greatest recent successes of theoretical chemistry. Only recently, however, has this concept been extended to materials science, where nonradiative recombination limits the efficiencies of materials for various optoelectronic applications. In this review, we present recent advances in the theoretical study of conical intersections in semiconductor nanomaterials. After briefly introducing conical intersections, we argue that specific defects in materials can induce conical intersections between the ground and first excited electronic states, thus introducing pathways for nonradiative recombination. We present recent developments in theoretical methods, computational tools, and chemical intuition for the prediction of such defect-induced conical intersections. Through examples in various nanomaterials, we illustrate the significance of conical intersections for nanoscience. We also discuss challenges facing research in this area and opportunities for progress.
Collapse
Affiliation(s)
- Benjamin G. Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Michael P. Esch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - B. Scott Fales
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Dylan T. Hardwick
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Wei-Tao Peng
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Yinan Shu
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| |
Collapse
|
19
|
Levine BG, Peng WT, Esch MP. Locality of conical intersections in semiconductor nanomaterials. Phys Chem Chem Phys 2019; 21:10870-10878. [PMID: 31106323 DOI: 10.1039/c9cp01584a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A predictive theory connecting atomic structure to the rate of recombination would enable the rational design of semiconductor nanomaterials for optoelectronic applications. Recently our group has demonstrated that the theoretical study of conical intersections can serve this purpose. Here we review recent work in this area, focusing on the thesis that low-energy conical intersections in nanomaterials share a common feature: locality. We define a conical intersection as local if (a) the intersecting states differ by the excitation of an electron between spatially local orbitals, and (b) the intersection is accessed when the energies of these orbitals are tuned by local distortions of the geometry. After illustrating the locality of the conical intersection responsible for recombination at dangling bond defects in silicon, we demonstrate the locality of low-energy conical intersections in cases where locality may be a surprise. First, we demonstrate the locality of low-energy self-trapped conical intersections in a pristine silicon nanocrystal, which has no defects that one would expect to serve as the center of a local intersection. Second, we demonstrate that the lowest energy intersection in a silicon system with two neighboring dangling bond defects localizes to a single defect site. We discuss the profound implications of locality for predicting the rate of recombination and suggest that the locality of intersections could be exploited in the experimental study of recombination, where spectroscopic studies of molecular models of defects could provide new insights.
Collapse
Affiliation(s)
- Benjamin G Levine
- Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA.
| | - Wei-Tao Peng
- Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA.
| | - Michael P Esch
- Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA.
| |
Collapse
|
20
|
The photochemical ring-opening of 1,3-cyclohexadiene imaged by ultrafast electron diffraction. Nat Chem 2019; 11:504-509. [DOI: 10.1038/s41557-019-0252-7] [Citation(s) in RCA: 111] [Impact Index Per Article: 22.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2018] [Accepted: 03/07/2019] [Indexed: 11/09/2022]
|
21
|
|
22
|
Sand AM, Hoyer CE, Truhlar DG, Gagliardi L. State-interaction pair-density functional theory. J Chem Phys 2018; 149:024106. [DOI: 10.1063/1.5036727] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Affiliation(s)
- Andrew M. Sand
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute, The University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Chad E. Hoyer
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute, The University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute, The University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute, The University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
23
|
Crespo-Otero R, Barbatti M. Recent Advances and Perspectives on Nonadiabatic Mixed Quantum–Classical Dynamics. Chem Rev 2018; 118:7026-7068. [DOI: 10.1021/acs.chemrev.7b00577] [Citation(s) in RCA: 301] [Impact Index Per Article: 50.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rachel Crespo-Otero
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| | | |
Collapse
|
24
|
Simine L, Lammert H, Sun L, Onuchic JN, Rossky PJ. Fluorescent Proteins Detect Host Structural Rearrangements via Electrostatic Mechanism. J Am Chem Soc 2018; 140:1203-1206. [PMID: 29328673 DOI: 10.1021/jacs.7b10851] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The rational design of genetically encoded fluorescent biosensors, which can detect rearrangements of target proteins via interdomain allostery, is hindered by the absence of mechanistic understanding of the underlying photophysics. Here, we focus on the modulation of fluorescence by mechanical perturbation in a popular biological probe: enhanced Green Fluorescent Protein (eGFP). Using a combination of molecular dynamics (MD) simulations and quantum chemistry, and a set of physically motivated assumptions, we construct a map of fluorescence quantum yield as a function of a 2D electric field imposed by the protein environment on the fluorophore. This map is transferable between Tsien's Class 2 GFP's, and it allows one to estimate the shifts in fluorescence intensity due to mechanical perturbations directly from MD simulations. We use it in combination with steered MD simulations to put forward a hypothesis for the mechanism of a genetically encoded voltage probe (ArcLight) whose mechanism is currently under debate.
Collapse
Affiliation(s)
- Lena Simine
- Departments of Chemistry, and ‡Physics, and §the Center for Theoretical Biological Physics, Rice University , Houston, Texas 77030, United States
| | - Heiko Lammert
- Departments of Chemistry, and ‡Physics, and §the Center for Theoretical Biological Physics, Rice University , Houston, Texas 77030, United States
| | - Li Sun
- Departments of Chemistry, and ‡Physics, and §the Center for Theoretical Biological Physics, Rice University , Houston, Texas 77030, United States
| | - José N Onuchic
- Departments of Chemistry, and ‡Physics, and §the Center for Theoretical Biological Physics, Rice University , Houston, Texas 77030, United States
| | - Peter J Rossky
- Departments of Chemistry, and ‡Physics, and §the Center for Theoretical Biological Physics, Rice University , Houston, Texas 77030, United States
| |
Collapse
|
25
|
Suchan J, Hollas D, Curchod BFE, Slavíček P. On the importance of initial conditions for excited-state dynamics. Faraday Discuss 2018; 212:307-330. [DOI: 10.1039/c8fd00088c] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The vast majority of ab initio excited-state simulations are performed within semiclassical, trajectory-based approaches. Apart from the underlying electronic-structure theory, the reliability of the simulations is controlled by a selection of initial conditions for the classical trajectories. We discuss appropriate choices of initial conditions for simulations of different experimental arrangements: dynamics initiated by continuum-wave (CW) laser fields or triggered by ultrashort laser pulses.
Collapse
Affiliation(s)
- Jiří Suchan
- Department of Physical Chemistry
- University of Chemistry and Technology, Prague
- 16628 Prague
- Czech Republic
| | - Daniel Hollas
- Department of Physical Chemistry
- University of Chemistry and Technology, Prague
- 16628 Prague
- Czech Republic
| | | | - Petr Slavíček
- Department of Physical Chemistry
- University of Chemistry and Technology, Prague
- 16628 Prague
- Czech Republic
| |
Collapse
|
26
|
Hollas D, Šištík L, Hohenstein EG, Martínez TJ, Slavíček P. Nonadiabatic Ab Initio Molecular Dynamics with the Floating Occupation Molecular Orbital-Complete Active Space Configuration Interaction Method. J Chem Theory Comput 2017; 14:339-350. [DOI: 10.1021/acs.jctc.7b00958] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Daniel Hollas
- Department
of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague 6, Czech Republic
| | - Lukáš Šištík
- Department
of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague 6, Czech Republic
| | - Edward G. Hohenstein
- Department
of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
- PhD
Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
| | - Todd J. Martínez
- Department
of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Petr Slavíček
- Department
of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague 6, Czech Republic
- J. Heyrovský Institute of Physical Chemistry, Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
27
|
Gozem S, Luk HL, Schapiro I, Olivucci M. Theory and Simulation of the Ultrafast Double-Bond Isomerization of Biological Chromophores. Chem Rev 2017; 117:13502-13565. [DOI: 10.1021/acs.chemrev.7b00177] [Citation(s) in RCA: 175] [Impact Index Per Article: 25.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Samer Gozem
- Department
of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| | - Hoi Ling Luk
- Chemistry
Department, Bowling Green State University, Overman Hall, Bowling Green, Ohio 43403, United States
| | - Igor Schapiro
- Fritz
Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Massimo Olivucci
- Chemistry
Department, Bowling Green State University, Overman Hall, Bowling Green, Ohio 43403, United States
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro
2, 53100 Siena, Italy
| |
Collapse
|
28
|
Shu Y, Fales BS, Peng WT, Levine BG. Understanding Nonradiative Recombination through Defect-Induced Conical Intersections. J Phys Chem Lett 2017; 8:4091-4099. [PMID: 28799771 DOI: 10.1021/acs.jpclett.7b01707] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Defects are known to introduce pathways for the nonradiative recombination of electronic excitations in semiconductors, but implicating a specific defect as a nonradiative center remains challenging for both experiment and theory. In this Perspective, we present recent progress toward this goal involving the identification and characterization of defect-induced conical intersections (DICIs), points of degeneracy between the ground and first excited electronic states of semiconductor materials that arise from the deformation of specific defects. Analysis of DICIs does not require the assumption of weak correlation between the electron and hole nor of stationary nuclei. It is demonstrated that in some cases an energetically accessible DICI is present even when no midgap state is predicted by single-particle theories (e.g., density functional theory). We review recent theoretical and computational developments that enable the location of DICIs in semiconductor nanomaterials and present insights into the photoluminescence of silicon nanocrystals gleaned from DICIs.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, University of Minnesota , Minneapolis, Minnesota 55455, United States
| | - B Scott Fales
- Department of Chemistry and the PULSE Institute, Stanford University , Stanford, California 94305, United States
- SLAC National Accelerator Laboratory , Menlo Park, California 94025, United States
| | - Wei-Tao Peng
- Department of Chemistry, Michigan State University , East Lansing, Michigan 48824, United States
| | - Benjamin G Levine
- Department of Chemistry, Michigan State University , East Lansing, Michigan 48824, United States
| |
Collapse
|