1
|
Hunter KE, Mao Y, Chin AW, Zuehlsdorff TJ. Environmentally Driven Symmetry Breaking Quenches Dual Fluorescence in Proflavine. J Phys Chem Lett 2024; 15:4623-4632. [PMID: 38647005 DOI: 10.1021/acs.jpclett.4c00677] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
Nonadiabatic couplings between several electronic excited states are ubiquitous in many organic chromophores and can significantly influence optical properties. A recent experimental study demonstrated that the proflavine molecule exhibits surprising dual fluorescence in the gas phase, which is suppressed in polar solvent environments. Here, we uncover the origin of this phenomenon by parametrizing a linear-vibronic coupling Hamiltonian from spectral densities of system-bath coupling constructed along molecular dynamics trajectories, fully accounting for interactions with the condensed-phase environment. The finite-temperature absorption, steady-state emission, and time-resolved emission spectra are then computed using powerful, numerically exact tensor network approaches. We find that the dual fluorescence in vacuum is driven by a single well-defined coupling mode but is quenched in solution due to dynamic solvent-driven symmetry breaking that mixes the two low-lying electronic states. We expect the computational framework developed here to be widely applicable to the study of non-Condon effects in complex condensed-phase environments.
Collapse
Affiliation(s)
- Kye E Hunter
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, United States
| | - Yuezhi Mao
- Department of Chemistry and Biochemistry, San Diego State University, San Diego, California 92182, United States
| | - Alex W Chin
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 place Jussieu, Paris 75005, France
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, United States
| |
Collapse
|
2
|
Allan L, Zuehlsdorff TJ. Taming the third order cumulant approximation to linear optical spectroscopy. J Chem Phys 2024; 160:074108. [PMID: 38380749 DOI: 10.1063/5.0182745] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Accepted: 01/23/2024] [Indexed: 02/22/2024] Open
Abstract
The second order cumulant method offers a promising pathway to predicting optical properties in condensed phase systems. It allows for the computation of linear absorption spectra from excitation energy fluctuations sampled along molecular dynamics (MD) trajectories, fully accounting for vibronic effects, direct solute-solvent interactions, and environmental polarization effects. However, the second order cumulant approximation only guarantees accurate line shapes for energy gap fluctuations obeying Gaussian statistics. A third order correction has recently been derived but often yields unphysical spectra or divergent line shapes for moderately non-Gaussian fluctuations due to the neglect of higher order terms in the cumulant expansion. In this work, we develop a corrected cumulant approach, where the collective effect of neglected higher order contributions is approximately accounted for through a dampening factor applied to the third order cumulant term. We show that this dampening factor can be expressed as a function of the skewness and kurtosis of energy gap fluctuations and can be parameterized from a large set of randomly sampled model Hamiltonians for which exact spectral line shapes are known. This approach is shown to systematically remove unphysical contributions in the form of negative absorbances from cumulant spectra in both model Hamiltonians and condensed phase systems sampled from MD and dramatically improves over the second order cumulant method in describing systems exhibiting Duschinsky mode mixing effects. We successfully apply the approach to the coumarin-153 dye in toluene, obtaining excellent agreement with experiment.
Collapse
Affiliation(s)
- Lucas Allan
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| |
Collapse
|
3
|
Ribó EG, Mao Z, Hirschi JS, Linsday T, Bach K, Walter ED, Simons CR, Zuehlsdorff TJ, Nyman M. Implementing vanadium peroxides as direct air carbon capture materials. Chem Sci 2024; 15:1700-1713. [PMID: 38303956 PMCID: PMC10829016 DOI: 10.1039/d3sc05381d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Accepted: 11/22/2023] [Indexed: 02/03/2024] Open
Abstract
Direct air capture (DAC) removal of anthropogenic CO2 from the atmosphere is imperative to slow the catastrophic effects of global climate change. Numerous materials are being investigated, including various alkaline inorganic metal oxides that form carbonates via DAC. Here we explore metastable early d0 transition metal peroxide molecules that undergo stabilization via multiple routes, including DAC. Specifically here, we describe via experiment and computation the mechanistic conversion of A3V(O2)4 (tetraperoxovanadate, A = K, Rb, Cs) to first a monocarbonate VO(O2)2(CO3)3-, and ultimately HKCO3 plus KVO4. Single crystal X-ray structures of rubidium and cesium tetraperoxovanadate are reported here for the first time, likely prior-challenged by instability. Infrared spectroscopy (FTIR), powder X-ray diffraction (PXRD), 51V solid state NMR (nuclear magnetic resonance), tandem thermogravimetry-mass spectrometry (TGA-MS) along with calculations (DFT, density functional theory) all converge on mechanisms of CO2 capture and release that involve the vanadium centre, despite the end product of a 300 days study being bicarbonate and metavanadate. Electron Paramagnetic Resonance (EPR) Spectroscopy along with a wet chemical assay and computational studies evidence the presense of ∼5% adventitous superoxide, likely formed by peroxide reduction of vanadium, which also stabilizes via the reaction with CO2. The alkalis have a profound effect on the stability of the peroxovanadate compounds, stability trending K > Rb > Cs. While this translates to more rapid CO2 capture with heavier alkalis, it does not necessarily lead to capture of more CO2. All compounds capture approximately two equivalents CO2 per vanadium centre. We cannot yet explain the reactivity trend of the alkali peroxovanadates, because any change in speciation of the alkalis from reactions to product is not quantifiable. This study sets the stage for understanding and implementing transition metal peroxide species, including peroxide-functionalized metal oxides, for DAC.
Collapse
Affiliation(s)
| | - Zhiwei Mao
- Department of Chemistry, Oregon State University Corvallis OR 97331 USA
| | - Jacob S Hirschi
- Department of Chemistry, Oregon State University Corvallis OR 97331 USA
| | - Taylor Linsday
- Department of Chemistry, Oregon State University Corvallis OR 97331 USA
| | - Karlie Bach
- Department of Chemistry, Oregon State University Corvallis OR 97331 USA
| | - Eric D Walter
- Pacific Northwest National Laboratory, Environmental Molecular Sciences Laboratory Richland WA 99352 USA
| | | | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University Corvallis OR 97331 USA
| | - May Nyman
- Department of Chemistry, Oregon State University Corvallis OR 97331 USA
| |
Collapse
|
4
|
Wiethorn ZR, Hunter KE, Zuehlsdorff TJ, Montoya-Castillo A. Beyond the Condon limit: Condensed phase optical spectra from atomistic simulations. J Chem Phys 2023; 159:244114. [PMID: 38153146 DOI: 10.1063/5.0180405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Accepted: 12/06/2023] [Indexed: 12/29/2023] Open
Abstract
While dark transitions made bright by molecular motions determine the optoelectronic properties of many materials, simulating such non-Condon effects in condensed phase spectroscopy remains a fundamental challenge. We derive a Gaussian theory to predict and analyze condensed phase optical spectra beyond the Condon limit. Our theory introduces novel quantities that encode how nuclear motions modulate the energy gap and transition dipole of electronic transitions in the form of spectral densities. By formulating the theory through a statistical framework of thermal averages and fluctuations, we circumvent the limitations of widely used microscopically harmonic theories, allowing us to tackle systems with generally anharmonic atomistic interactions and non-Condon fluctuations of arbitrary strength. We show how to calculate these spectral densities using first-principles simulations, capturing realistic molecular interactions and incorporating finite-temperature, disorder, and dynamical effects. Our theory accurately predicts the spectra of systems known to exhibit strong non-Condon effects (phenolate in various solvents) and reveals distinct mechanisms for electronic peak splitting: timescale separation of modes that tune non-Condon effects and spectral interference from correlated energy gap and transition dipole fluctuations. We further introduce analysis tools to identify how intramolecular vibrations, solute-solvent interactions, and environmental polarization effects impact dark transitions. Moreover, we prove an upper bound on the strength of cross correlated energy gap and transition dipole fluctuations, thereby elucidating a simple condition that a system must follow for our theory to accurately predict its spectrum.
Collapse
Affiliation(s)
- Zachary R Wiethorn
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Kye E Hunter
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | | |
Collapse
|
5
|
Chen MS, Mao Y, Snider A, Gupta P, Montoya-Castillo A, Zuehlsdorff TJ, Isborn CM, Markland TE. Elucidating the Role of Hydrogen Bonding in the Optical Spectroscopy of the Solvated Green Fluorescent Protein Chromophore: Using Machine Learning to Establish the Importance of High-Level Electronic Structure. J Phys Chem Lett 2023; 14:6610-6619. [PMID: 37459252 DOI: 10.1021/acs.jpclett.3c01444] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
Hydrogen bonding interactions with chromophores in chemical and biological environments play a key role in determining their electronic absorption and relaxation processes, which are manifested in their linear and multidimensional optical spectra. For chromophores in the condensed phase, the large number of atoms needed to simulate the environment has traditionally prohibited the use of high-level excited-state electronic structure methods. By leveraging transfer learning, we show how to construct machine-learned models to accurately predict the high-level excitation energies of a chromophore in solution from only 400 high-level calculations. We show that when the electronic excitations of the green fluorescent protein chromophore in water are treated using EOM-CCSD embedded in a DFT description of the solvent the optical spectrum is correctly captured and that this improvement arises from correctly treating the coupling of the electronic transition to electric fields, which leads to a larger response upon hydrogen bonding between the chromophore and water.
Collapse
Affiliation(s)
- Michael S Chen
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Andrew Snider
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Prachi Gupta
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Andrés Montoya-Castillo
- Department of Chemistry, University of Colorado, Boulder, Boulder, Colorado 80309, United States
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, United States
| | - Christine M Isborn
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
6
|
Smith KT, Hunter K, Chiu NC, Zhuang H, Jumrusprasert P, Stickle WF, Reimer JA, Zuehlsdorff TJ, Stylianou K. Hypsochromically‐shifted Emission of Metal‐organic Frameworks Generated through Post‐synthetic Ligand Reduction. Angew Chem Int Ed Engl 2023. [DOI: 10.1002/ange.202302123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/18/2023]
Affiliation(s)
| | - Kye Hunter
- Oregon State University Chemistry UNITED STATES
| | | | - Hao Zhuang
- UC Berkeley: University of California Berkeley Chemical and Biomolecular Engineering UNITED STATES
| | | | | | - Jeffrey A. Reimer
- UC Berkeley: University of California Berkeley Chemical and Biomolecular Engineering UNITED STATES
| | | | - Kyriakos Stylianou
- Oregon State University Chemistry 153 Gilbert Hall 97330 Corvallis UNITED STATES
| |
Collapse
|
7
|
Smith KT, Hunter K, Chiu NC, Zhuang H, Jumrusprasert P, Stickle WF, Reimer JA, Zuehlsdorff TJ, Stylianou K. Hypsochromically-shifted Emission of Metal-organic Frameworks Generated through Post-synthetic Ligand Reduction. Angew Chem Int Ed Engl 2023:e202302123. [PMID: 36929127 DOI: 10.1002/anie.202302123] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2023] [Revised: 03/14/2023] [Accepted: 03/16/2023] [Indexed: 03/18/2023]
Abstract
Luminescent materials with tunable emission are becoming increasingly desirable as we move towards needing efficient Light Emitting Diodes (LEDs) for displays. Key to developing better displays is the advancement of strategies for rationally designing emissive materials that are tunable and efficient. We report a series of emissive metal-organic frameworks (MOFs) generated using BUT-10 (BUT: Beijing University of Technology) that emits green light with λmax at 525 nm. Post-synthetic reduction of the ketone on the fluorenone ligand in BUT-10 generates new materials, BUT-10-M and BUT-10-R. The emission for BUT-10-R is hypsochromically-shifted by 113 nm. Multivariate BUT-10-M structures demonstrate emission with two maxima corresponding to the emission of both fluorenol and fluorenone moieties present in their structures. Our study represents a novel post-synthetic ligand reduction strategy for producing emissive MOFs with tunable emission ranging from green, white-blue to deep blue.
Collapse
Affiliation(s)
- Kyle T Smith
- Oregon State University, chemistry, UNITED STATES
| | - Kye Hunter
- Oregon State University, Chemistry, UNITED STATES
| | | | - Hao Zhuang
- UC Berkeley: University of California Berkeley, Chemical and Biomolecular Engineering, UNITED STATES
| | | | | | - Jeffrey A Reimer
- UC Berkeley: University of California Berkeley, Chemical and Biomolecular Engineering, UNITED STATES
| | | | - Kyriakos Stylianou
- Oregon State University, Chemistry, 153 Gilbert Hall, 97330, Corvallis, UNITED STATES
| |
Collapse
|
8
|
Abou Taka A, Lu SY, Gowland D, Zuehlsdorff TJ, Corzo HH, Pribram-Jones A, Shi L, Hratchian HP, Isborn CM. Comparison of Linear Response Theory, Projected Initial Maximum Overlap Method, and Molecular Dynamics-Based Vibronic Spectra: The Case of Methylene Blue. J Chem Theory Comput 2022; 18:3039-3051. [PMID: 35472264 DOI: 10.1021/acs.jctc.1c01127] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The simulation of optical spectra is essential to molecular characterization and, in many cases, critical for interpreting experimental spectra. The most common method for simulating vibronic absorption spectra relies on the geometry optimization and computation of normal modes for ground and excited electronic states. In this report, we show that the utilization of such a procedure within an adiabatic linear response (LR) theory framework may lead to state mixings and a breakdown of the Born-Oppenheimer approximation, resulting in a poor description of absorption spectra. In contrast, computing excited states via a self-consistent field method in conjunction with a maximum overlap model produces states that are not subject to such mixings. We show that this latter method produces vibronic spectra much more aligned with vertical gradient and molecular dynamics (MD) trajectory-based approaches. For the methylene blue chromophore, we compare vibronic absorption spectra computed with the following: an adiabatic Hessian approach with LR theory-optimized structures and normal modes, a vertical gradient procedure, the Hessian and normal modes of maximum overlap method-optimized structures, and excitation energy time-correlation functions generated from an MD trajectory. Because of mixing between the bright S1 and dark S2 surfaces near the S1 minimum, computing the adiabatic Hessian with LR theory and time-dependent density functional theory with the B3LYP density functional predicts a large vibronic shoulder for the absorption spectrum that is not present for any of the other methods. Spectral densities are analyzed and we compare the behavior of the key normal mode that in LR theory strongly couples to the optical excitation while showing S1/S2 state mixings. Overall, our study provides a note of caution in computing vibronic spectra using the excited-state adiabatic Hessian of LR theory-optimized structures and also showcases three alternatives that are less sensitive to adiabatic state mixing effects.
Collapse
Affiliation(s)
- Ali Abou Taka
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Shao-Yu Lu
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Duncan Gowland
- Department of Physics, King's College London, London WC2R 2LS, United Kingdom
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, United States
| | - Hector H Corzo
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Aurora Pribram-Jones
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Hrant P Hratchian
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Christine M Isborn
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
9
|
Abstract
A water-stable, porphyrin-based metal-organic framework (MOF) produces a distinct colour change in response to acids' pKa and concentrations. This colour change is associated with the protonation of the N-atoms within the porphyrin ligand present in the MOF structure. As a proof-of-concept, we demonstrate the use of this MOF for detecting traces of different acidic pesticides present in water samples spontaneously.
Collapse
Affiliation(s)
- Kyle T Smith
- Department of Chemistry, Oregon State University, 153 Gilbert Hall, Corvallis, OR, 97331-4003, USA.
| | - Chloe A Ramsperger
- Department of Chemistry, Oregon State University, 153 Gilbert Hall, Corvallis, OR, 97331-4003, USA.
| | - Kye E Hunter
- Department of Chemistry, Oregon State University, 153 Gilbert Hall, Corvallis, OR, 97331-4003, USA.
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, 153 Gilbert Hall, Corvallis, OR, 97331-4003, USA.
| | - Kyriakos C Stylianou
- Department of Chemistry, Oregon State University, 153 Gilbert Hall, Corvallis, OR, 97331-4003, USA.
| |
Collapse
|
10
|
Lu SY, Zuehlsdorff TJ, Hong H, Aguirre VP, Isborn CM, Shi L. The Influence of Electronic Polarization on Nonlinear Optical Spectroscopy. J Phys Chem B 2021; 125:12214-12227. [PMID: 34726915 DOI: 10.1021/acs.jpcb.1c05914] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The environment surrounding a chromophore can dramatically affect the energy absorption and relaxation process, as manifested in optical spectra. Simulations of nonlinear optical spectroscopy, such as two-dimensional electronic spectroscopy (2DES) and transient absorption (TA), will be influenced by the computational model of the environment. We here compare a fixed point charge molecular mechanics model and a quantum mechanical (QM) model of the environment in computed 2DES and TA spectra of Nile red in water and the chromophore of photoactive yellow protein (PYP) in water and protein environments. In addition to simulating these nonlinear optical spectra, we directly juxtapose the computed excitation energy correlation function to the dynamic Stokes shift function often used to analyze environment dynamics. Overall, we find that for the three systems studied here the mutual electronic polarization provided by the QM environment manifests in broader 2DES signals, as well as a larger reorganization energy and a larger static Stokes shift due to stronger coupling between the chromophore and the environment.
Collapse
Affiliation(s)
- Shao-Yu Lu
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, United States
| | - Hanbo Hong
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Vincent P Aguirre
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Christine M Isborn
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
11
|
Dunnett AJ, Gowland D, Isborn CM, Chin AW, Zuehlsdorff TJ. Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue. J Chem Phys 2021; 155:144112. [PMID: 34654312 DOI: 10.1063/5.0062950] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute-solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute-solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.
Collapse
Affiliation(s)
- Angus J Dunnett
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 place Jussieu, 75005 Paris, France
| | - Duncan Gowland
- Department of Physics, King's College London, London WC2R 2LS, United Kingdom
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Alex W Chin
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 place Jussieu, 75005 Paris, France
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| |
Collapse
|
12
|
Zuehlsdorff TJ, Shedge SV, Lu SY, Hong H, Aguirre VP, Shi L, Isborn CM. Vibronic and Environmental Effects in Simulations of Optical Spectroscopy. Annu Rev Phys Chem 2021; 72:165-188. [DOI: 10.1146/annurev-physchem-090419-051350] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.
Collapse
Affiliation(s)
- Tim J. Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Sapana V. Shedge
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Shao-Yu Lu
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Hanbo Hong
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Vincent P. Aguirre
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Liang Shi
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Christine M. Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| |
Collapse
|
13
|
Shedge SV, Zuehlsdorff TJ, Khanna A, Conley S, Isborn CM. Explicit environmental and vibronic effects in simulations of linear and nonlinear optical spectroscopy. J Chem Phys 2021; 154:084116. [PMID: 33639769 DOI: 10.1063/5.0038196] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurately simulating the linear and nonlinear electronic spectra of condensed phase systems and accounting for all physical phenomena contributing to spectral line shapes presents a significant challenge. Vibronic transitions can be captured through a harmonic model generated from the normal modes of a chromophore, but it is challenging to also include the effects of specific chromophore-environment interactions within such a model. We work to overcome this limitation by combining approaches to account for both explicit environment interactions and vibronic couplings for simulating both linear and nonlinear optical spectra. We present and show results for three approaches of varying computational cost for combining ensemble sampling of chromophore-environment configurations with Franck-Condon line shapes for simulating linear spectra. We present two analogous approaches for nonlinear spectra. Simulated absorption spectra and two-dimensional electronic spectra (2DES) are presented for the Nile red chromophore in different solvent environments. Employing an average Franck-Condon or 2DES line shape appears to be a promising method for simulating linear and nonlinear spectroscopy for a chromophore in the condensed phase.
Collapse
Affiliation(s)
- Sapana V Shedge
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Ajay Khanna
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Stacey Conley
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
14
|
Chen MS, Zuehlsdorff TJ, Morawietz T, Isborn CM, Markland TE. Exploiting Machine Learning to Efficiently Predict Multidimensional Optical Spectra in Complex Environments. J Phys Chem Lett 2020; 11:7559-7568. [PMID: 32808797 DOI: 10.1021/acs.jpclett.0c02168] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
The excited-state dynamics of chromophores in complex environments determine a range of vital biological and energy capture processes. Time-resolved, multidimensional optical spectroscopies provide a key tool to investigate these processes. Although theory has the potential to decode these spectra in terms of the electronic and atomistic dynamics, the need for large numbers of excited-state electronic structure calculations severely limits first-principles predictions of multidimensional optical spectra for chromophores in the condensed phase. Here, we leverage the locality of chromophore excitations to develop machine learning models to predict the excited-state energy gap of chromophores in complex environments for efficiently constructing linear and multidimensional optical spectra. By analyzing the performance of these models, which span a hierarchy of physical approximations, across a range of chromophore-environment interaction strengths, we provide strategies for the construction of machine learning models that greatly accelerate the calculation of multidimensional optical spectra from first principles.
Collapse
Affiliation(s)
- Michael S Chen
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Tobias Morawietz
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
15
|
Zuehlsdorff TJ, Hong H, Shi L, Isborn CM. Nonlinear spectroscopy in the condensed phase: The role of Duschinsky rotations and third order cumulant contributions. J Chem Phys 2020; 153:044127. [PMID: 32752702 DOI: 10.1063/5.0013739] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
First-principles modeling of nonlinear optical spectra in the condensed phase is highly challenging because both environment and vibronic interactions can play a large role in determining spectral shapes and excited state dynamics. Here, we compute two dimensional electronic spectroscopy (2DES) signals based on a cumulant expansion of the energy gap fluctuation operator, with specific focus on analyzing mode mixing effects introduced by the Duschinsky rotation and the role of the third order term in the cumulant expansion for both model and realistic condensed phase systems. We show that for a harmonic model system, the third order cumulant correction captures effects introduced by a mismatch in curvatures of ground and excited state potential energy surfaces, as well as effects of mode mixing. We also demonstrate that 2DES signals can be accurately reconstructed from purely classical correlation functions using quantum correction factors. We then compute nonlinear optical spectra for the Nile red and methylene blue chromophores in solution, assessing the third order cumulant contribution for realistic systems. We show that the third order cumulant correction is strongly dependent on the treatment of the solvent environment, revealing the interplay between environmental polarization and the electronic-vibrational coupling.
Collapse
Affiliation(s)
- Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Hanbo Hong
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Liang Shi
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
16
|
Prentice JCA, Aarons J, Womack JC, Allen AEA, Andrinopoulos L, Anton L, Bell RA, Bhandari A, Bramley GA, Charlton RJ, Clements RJ, Cole DJ, Constantinescu G, Corsetti F, Dubois SMM, Duff KKB, Escartín JM, Greco A, Hill Q, Lee LP, Linscott E, O'Regan DD, Phipps MJS, Ratcliff LE, Serrano ÁR, Tait EW, Teobaldi G, Vitale V, Yeung N, Zuehlsdorff TJ, Dziedzic J, Haynes PD, Hine NDM, Mostofi AA, Payne MC, Skylaris CK. The ONETEP linear-scaling density functional theory program. J Chem Phys 2020; 152:174111. [PMID: 32384832 DOI: 10.1063/5.0004445] [Citation(s) in RCA: 51] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
We present an overview of the onetep program for linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the density matrix and the NGWFs are optimized with localization constraints. By taking advantage of localization, onetep is able to perform calculations including thousands of atoms with computational effort, which scales linearly with the number or atoms. The code has a large and diverse range of capabilities, explored in this paper, including different boundary conditions, various exchange-correlation functionals (with and without exact exchange), finite electronic temperature methods for metallic systems, methods for strongly correlated systems, molecular dynamics, vibrational calculations, time-dependent DFT, electronic transport, core loss spectroscopy, implicit solvation, quantum mechanical (QM)/molecular mechanical and QM-in-QM embedding, density of states calculations, distributed multipole analysis, and methods for partitioning charges and interactions between fragments. Calculations with onetep provide unique insights into large and complex systems that require an accurate atomic-level description, ranging from biomolecular to chemical, to materials, and to physical problems, as we show with a small selection of illustrative examples. onetep has always aimed to be at the cutting edge of method and software developments, and it serves as a platform for developing new methods of electronic structure simulation. We therefore conclude by describing some of the challenges and directions for its future developments and applications.
Collapse
Affiliation(s)
- Joseph C A Prentice
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Jolyon Aarons
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - James C Womack
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Alice E A Allen
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Lampros Andrinopoulos
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Lucian Anton
- UKAEA, Culham Science Centre, Abingdon OX14 3DB, United Kingdom
| | - Robert A Bell
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Arihant Bhandari
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Gabriel A Bramley
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Robert J Charlton
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Rebecca J Clements
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Daniel J Cole
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne NE1 7RU, United Kingdom
| | - Gabriel Constantinescu
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Fabiano Corsetti
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Simon M-M Dubois
- Institute of Condensed Matter and Nanosciences, Université Catholique de Louvain, Louvain-la-Neuve, Belgium
| | - Kevin K B Duff
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - José María Escartín
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Andrea Greco
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Quintin Hill
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Louis P Lee
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Edward Linscott
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - David D O'Regan
- School of Physics, AMBER, and CRANN Institute, Trinity College Dublin, The University of Dublin, Dublin 2, Ireland
| | - Maximillian J S Phipps
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Laura E Ratcliff
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Álvaro Ruiz Serrano
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Edward W Tait
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Gilberto Teobaldi
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Valerio Vitale
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Nelson Yeung
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Peter D Haynes
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Nicholas D M Hine
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Arash A Mostofi
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Mike C Payne
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
17
|
Affiliation(s)
- Tim J. Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Hanbo Hong
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Christine M. Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
18
|
Zuehlsdorff TJ, Montoya-Castillo A, Napoli JA, Markland TE, Isborn CM. Optical spectra in the condensed phase: Capturing anharmonic and vibronic features using dynamic and static approaches. J Chem Phys 2019; 151:074111. [PMID: 31438704 DOI: 10.1063/1.5114818] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Abstract
Simulating optical spectra in the condensed phase remains a challenge for theory due to the need to capture spectral signatures arising from anharmonicity and dynamical effects, such as vibronic progressions and asymmetry. As such, numerous simulation methods have been developed that invoke different approximations and vary in their ability to capture different physical regimes. Here, we use several models of chromophores in the condensed phase and ab initio molecular dynamics simulations to rigorously assess the applicability of methods to simulate optical absorption spectra. Specifically, we focus on the ensemble scheme, which can address anharmonic potential energy surfaces but relies on the applicability of extreme nuclear-electronic time scale separation; the Franck-Condon method, which includes dynamical effects but generally only at the harmonic level; and the recently introduced ensemble zero-temperature Franck-Condon approach, which straddles these limits. We also devote particular attention to the performance of methods derived from a cumulant expansion of the energy gap fluctuations and test the ability to approximate the requisite time correlation functions using classical dynamics with quantum correction factors. These results provide insights as to when these methods are applicable and able to capture the features of condensed phase spectra qualitatively and, in some cases, quantitatively across a range of regimes.
Collapse
Affiliation(s)
- Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | | | - Joseph A Napoli
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
19
|
Shedge SV, Zuehlsdorff TJ, Servis MJ, Clark AE, Isborn CM. Effect of Ions on the Optical Absorption Spectra of Aqueously Solvated Chromophores. J Phys Chem A 2019; 123:6175-6184. [DOI: 10.1021/acs.jpca.9b03163] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Sapana V. Shedge
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Tim J. Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Michael J. Servis
- Department of Chemistry and the Material Science and Engineering Program, Washington State University, Pullman, Washington 99164, United States
| | - Aurora E. Clark
- Department of Chemistry and the Material Science and Engineering Program, Washington State University, Pullman, Washington 99164, United States
- Voiland School of Chemical Engineering and Bioengineering, Washington State University, Pullman, Washington 99164, United States
| | - Christine M. Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
20
|
Zuehlsdorff TJ, Napoli JA, Milanese JM, Markland TE, Isborn CM. Unraveling electronic absorption spectra using nuclear quantum effects: Photoactive yellow protein and green fluorescent protein chromophores in water. J Chem Phys 2018; 149:024107. [PMID: 30007372 DOI: 10.1063/1.5025517] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Many physical phenomena must be accounted for to accurately model solution-phase optical spectral line shapes, from the sampling of chromophore-solvent configurations to the electronic-vibrational transitions leading to vibronic fine structure. Here we thoroughly explore the role of nuclear quantum effects, direct and indirect solvent effects, and vibronic effects in the computation of the optical spectrum of the aqueously solvated anionic chromophores of green fluorescent protein and photoactive yellow protein. By analyzing the chromophore and solvent configurations, the distributions of vertical excitation energies, the absorption spectra computed within the ensemble approach, and the absorption spectra computed within the ensemble plus zero-temperature Franck-Condon approach, we show how solvent, nuclear quantum effects, and vibronic transitions alter the optical absorption spectra. We find that including nuclear quantum effects in the sampling of chromophore-solvent configurations using ab initio path integral molecular dynamics simulations leads to improved spectral shapes through three mechanisms. The three mechanisms that lead to line shape broadening and a better description of the high-energy tail are softening of heavy atom bonds in the chromophore that couple to the optically bright state, widening the distribution of vertical excitation energies from more diverse solvation environments, and redistributing spectral weight from the 0-0 vibronic transition to higher energy vibronic transitions when computing the Franck-Condon spectrum in a frozen solvent pocket. The absorption spectra computed using the combined ensemble plus zero-temperature Franck-Condon approach yield significant improvements in spectral shape and width compared to the spectra computed with the ensemble approach. Using the combined approach with configurations sampled from path integral molecular dynamics trajectories presents a significant step forward in accurately modeling the absorption spectra of aqueously solvated chromophores.
Collapse
Affiliation(s)
- Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Joseph A Napoli
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Joel M Milanese
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
21
|
Zuehlsdorff TJ, Isborn CM. Combining the ensemble and Franck-Condon approaches for calculating spectral shapes of molecules in solution. J Chem Phys 2018; 148:024110. [PMID: 29331131 DOI: 10.1063/1.5006043] [Citation(s) in RCA: 52] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The correct treatment of vibronic effects is vital for the modeling of absorption spectra of many solvated dyes. Vibronic spectra for small dyes in solution can be easily computed within the Franck-Condon approximation using an implicit solvent model. However, implicit solvent models neglect specific solute-solvent interactions on the electronic excited state. On the other hand, a straightforward way to account for solute-solvent interactions and temperature-dependent broadening is by computing vertical excitation energies obtained from an ensemble of solute-solvent conformations. Ensemble approaches usually do not account for vibronic transitions and thus often produce spectral shapes in poor agreement with experiment. We address these shortcomings by combining zero-temperature vibronic fine structure with vertical excitations computed for a room-temperature ensemble of solute-solvent configurations. In this combined approach, all temperature-dependent broadening is treated classically through the sampling of configurations and quantum mechanical vibronic contributions are included as a zero-temperature correction to each vertical transition. In our calculation of the vertical excitations, significant regions of the solvent environment are treated fully quantum mechanically to account for solute-solvent polarization and charge-transfer. For the Franck-Condon calculations, a small amount of frozen explicit solvent is considered in order to capture solvent effects on the vibronic shape function. We test the proposed method by comparing calculated and experimental absorption spectra of Nile red and the green fluorescent protein chromophore in polar and non-polar solvents. For systems with strong solute-solvent interactions, the combined approach yields significant improvements over the ensemble approach. For systems with weak to moderate solute-solvent interactions, both the high-energy vibronic tail and the width of the spectra are in excellent agreement with experiments.
Collapse
Affiliation(s)
- T J Zuehlsdorff
- School of Natural Sciences, University of California Merced, N. Lake Road, Merced, California 95344, USA
| | - C M Isborn
- School of Natural Sciences, University of California Merced, N. Lake Road, Merced, California 95344, USA
| |
Collapse
|
22
|
Zuehlsdorff TJ, Haynes PD, Payne MC, Hine NDM. Predicting solvatochromic shifts and colours of a solvated organic dye: The example of nile red. J Chem Phys 2018; 146:124504. [PMID: 28388154 DOI: 10.1063/1.4979196] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The solvatochromic shift, as well as the change in colour of the simple organic dye nile red, is studied in two polar and two non-polar solvents in the context of large-scale time-dependent density-functional theory (TDDFT) calculations treating large parts of the solvent environment from first principles. We show that an explicit solvent representation is vital to resolve absorption peak shifts between nile red in n-hexane and toluene, as well as acetone and ethanol. The origin of the failure of implicit solvent models for these solvents is identified as being due to the strong solute-solvent interactions in form of π-stacking and hydrogen bonding in the case of toluene and ethanol. We furthermore demonstrate that the failures of the computationally inexpensive Perdew-Burke-Ernzerhof (PBE) functional in describing some features of the excited state potential energy surface of the S1 state of nile red can be corrected for in a straightforward fashion, relying only on a small number of calculations making use of more sophisticated range-separated hybrid functionals. The resulting solvatochromic shifts and predicted colours are in excellent agreement with experiment, showing the computational approach outlined in this work to yield very robust predictions of optical properties of dyes in solution.
Collapse
Affiliation(s)
- T J Zuehlsdorff
- School of Natural Sciences, University of California Merced, 5200 N. Lake Road, Merced, California 95343, USA
| | - P D Haynes
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - M C Payne
- Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - N D M Hine
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
23
|
Joseph J, Baumann KN, Koehler P, Zuehlsdorff TJ, Cole DJ, Weber J, Bohndiek SE, Hernández-Ainsa S. Distance dependent photoacoustics revealed through DNA nanostructures. Nanoscale 2017; 9:16193-16199. [PMID: 29043366 DOI: 10.1039/c7nr05353c] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Molecular rulers that rely on the Förster resonance energy transfer (FRET) mechanism are widely used to investigate dynamic molecular processes that occur on the nanometer scale. However, the capabilities of these fluorescence molecular rulers are fundamentally limited to shallow imaging depths by light scattering in biological samples. Photoacoustic tomography (PAT) has recently emerged as a high resolution modality for in vivo imaging, coupling optical excitation with ultrasound detection. In this paper, we report the capability of PAT to probe distance-dependent FRET at centimeter depths. Using DNA nanotechnology we created several nanostructures with precisely positioned fluorophore-quencher pairs over a range of nanoscale separation distances. PAT of the DNA nanostructures showed distance-dependent photoacoustic signal enhancement and demonstrated the ability of PAT to reveal the FRET process deep within tissue mimicking phantoms. Further, we experimentally validated these DNA nanostructures as a novel and biocompatible strategy to augment the intrinsic photoacoustic signal generation capabilities of small molecule fluorescent dyes.
Collapse
Affiliation(s)
- James Joseph
- Cavendish Laboratory, Department of Physics, University of Cambridge, Cambridge, UK.
| | | | | | | | | | | | | | | |
Collapse
|
24
|
Zuehlsdorff TJ, Hine NDM, Payne MC, Haynes PD. Linear-scaling time-dependent density-functional theory beyond the Tamm-Dancoff approximation: Obtaining efficiency and accuracy with in situ optimised local orbitals. J Chem Phys 2015; 143:204107. [PMID: 26627950 DOI: 10.1063/1.4936280] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present a solution of the full time-dependent density-functional theory (TDDFT) eigenvalue equation in the linear response formalism exhibiting a linear-scaling computational complexity with system size, without relying on the simplifying Tamm-Dancoff approximation (TDA). The implementation relies on representing the occupied and unoccupied subspaces with two different sets of in situ optimised localised functions, yielding a very compact and efficient representation of the transition density matrix of the excitation with the accuracy associated with a systematic basis set. The TDDFT eigenvalue equation is solved using a preconditioned conjugate gradient algorithm that is very memory-efficient. The algorithm is validated on a small test molecule and a good agreement with results obtained from standard quantum chemistry packages is found, with the preconditioner yielding a significant improvement in convergence rates. The method developed in this work is then used to reproduce experimental results of the absorption spectrum of bacteriochlorophyll in an organic solvent, where it is demonstrated that the TDA fails to reproduce the main features of the low energy spectrum, while the full TDDFT equation yields results in good qualitative agreement with experimental data. Furthermore, the need for explicitly including parts of the solvent into the TDDFT calculations is highlighted, making the treatment of large system sizes necessary that are well within reach of the capabilities of the algorithm introduced here. Finally, the linear-scaling properties of the algorithm are demonstrated by computing the lowest excitation energy of bacteriochlorophyll in solution. The largest systems considered in this work are of the same order of magnitude as a variety of widely studied pigment-protein complexes, opening up the possibility of studying their properties without having to resort to any semiclassical approximations to parts of the protein environment.
Collapse
Affiliation(s)
- T J Zuehlsdorff
- Cavendish Laboratory, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - N D M Hine
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - M C Payne
- Cavendish Laboratory, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - P D Haynes
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| |
Collapse
|