Quantum optimal control of molecular isomerization in the presence of a competing dissociation channel

The quantum optimal control of isomerization in the presence of a competing dissociation channel is simulated on a two-dimensional model. The control of isomerization of a hydrogen atom is achieved through vibrational transitions on the ground-state surface as well as with the aid of an excited-state surface. The effects of different competing dissociation channel configurations on the isomerization control are explored. Suppression of the competing dissociation dynamics during the isomerization control on the ground-state surface becomes easier with an increase in the spatial separation between the isomerization and dissociation regions and with a decrease in the dissociation channel width. Isomerization control first involving transfer of amplitude to an excited-state surface is less influenced by the dissociation channel configuration on the ground-state surface, even in cases where the excited-state surface allows for a moderate spreading of the excited wave packet.     

Femtosecond study on the isomerization dynamics of NK88. II. Excited-state dynamics

The molecule 3,3(')-diethyl-2,2(')-thiacyanine isomerizes after irradiation with light of the proper wavelength. After excitation, it undergoes a transition, in which one or more conical intersections are involved, back to the ground state to form different product photoisomers. The dynamics before and directly after the transition back to the ground state is investigated by transient absorption spectroscopy in a wavelength region of 360-950 nm, as well as by fluorescence upconversion. It is shown that the excited-state dynamics are governed by two time scales: a short one with a decay time of less than 2 ps and a long one with about 9 ps. A thorough comparison of the experimental results with those of configuration interaction singles and time-dependent density functional theory calculations suggests that these dynamics are related to two competing pathways differing in the molecular twisting on the excited surface after photoexcitation. From the experimental point of view this picture arises taking into account the time scales for ground-state bleach, excited-state absorption, stimulated emission, fluorescence, and assumed hot ground-state absorption both in the solvent methanol and ethylene glycol.     

On the Formation of Cyclobutane Pyrimidine Dimers in UV‐irradiated DNA: Why are Thymines More Reactive?¶

The reaction pathways for thermal and photochemical formation of cyclobutane pyrimidine dimers in DNA are explored using density functional theory techniques. Although it is found that the thermal [2 + 2] cycloadditions of thymine + thymine (T + T --> T x T), cytosine + cytosine (C + C --> C x C) and cytosine + thymine (C + T --> C x T) all are similarly unfavorable in terms of energy barriers and reaction energies, the excited-state energy curves associated with the corresponding photochemical cycloadditions display differences that--in line with experimental findings--unanimously point to the predominance of T x T in UV-irradiated DNA. It is shown that the photocycloaddition of thymines is facilitated by the fact that the S1 state of the corresponding reactant complex lies comparatively high in energy. Moreover, at a nuclear configuration coinciding with the ground-state transition structure, the excited-state energy curve displays an absolute minimum only for the T + T system. Finally, the T + T system is also associated with the most favorable excited-state energy barriers and has the smallest S2-S0 energy gap at the ground-state transition structure.     

Time-dependent Density Functional Theory Study on the Electronically Excited States of N-Methylformamide in Aqueous Solution

Excited-state hydrogen-bonding dynamics of N-methylformamide (NMF) in water has been investigated by time-dependent density functional theory (TDDFT) method. The ground-state geometry optimizations were calculated by density functional theory (DFT) method, while the electronic transition energies and corresponding oscillation strengths of the low-lying electronically excited states of isolated NMF, water monomers and the hydrogen-bonded NMF-H 2 O were calculated by TDDFT method. According to Zhao's rule on the excited-state hydrogen bonding dynamics, our results demonstrate that the intermolecular hydrogen bond C=O···O-H is strengthened and weakened in different electronically excited states. The hydrogen bond strengthening and weakening in the electronically excited state plays an important role in the photophysics of NMF in solutions.     

Solvent effect on the excited-state dynamics of analogues of the photoactive yellow protein chromophore

We previously reported that two analogues of the Photoactive Yellow Protein chromophore, trans-p-hydroxycinnamic acid (pCA(2-)) and its amide derivative (pCM-) in their deprotonated forms, undergo a trans-cis photoisomerization whereas the thioester derivative, trans-p-hydroxythiophenyl cinnamate (pCT-), does not. pCT- is also the only one to exhibit a short-lived intermediate on its excited-state deactivation pathway. We here further stress the existence of two different relaxation mechanisms for these molecules and examine the reaction coordinates involved. We looked at the effect of the solvent properties (viscosity, polarity, solvation dynamics) on their excited-state relaxation dynamics, probed by ultrafast transient absorption spectroscopy. Sensitivity to the solvent properties is found to be larger for pCT- than for pCA(2-) and pCM-. This difference is considered to reveal that either the relaxation pathway or the reaction coordinate is different for these two classes of analogues. It is also found to be correlated to the electron donor-acceptor character of the molecule. We attribute the excited-state deactivation of analogues bearing a weaker acceptor group, pCA(2-) and pCM-, to a stilbene-like photoisomerization mechanism with the concerted rotation of the ethylenic bond and one adjacent single bond. For pCT-, which contains a stronger acceptor group, we consider a photoisomerization mechanism mainly involving the single torsion of the ethylenic bond. The excited-state deactivation of pCT- would lead to the formation of a ground-state intermediate at the "perp" geometry, which would return to the initial trans conformation without net isomerization.     

Excited-state intermolecular proton transfer reactions of 7-azaindole(MeOH)(n) (n = 1-3) clusters in the gas phase: on-the-fly dynamics simulation

Ultrafast excited-state intermolecular proton transfer (PT) reactions in 7-azaindole(methanol)(n) (n = 1-3) [7AI(MeOH)(n=1-3)] complexes were performed using dynamics simulations. These complexes were first optimized at the RI-ADC(2)/SVP-SV(P) level in the gas phase. The ground-state structures with the lowest energy were also investigated and presented. On-the-fly dynamics simulations for the first-excited state were employed to investigate reaction mechanisms and time evolution of PT processes. The PT characteristics of the reactions were confirmed by the nonexistence of crossings between S(ππ*) and S(πσ*) states. Excited-state dynamics results for all complexes exhibit excited-state multiple-proton transfer (ESmultiPT) reactions via methanol molecules along an intermolecular hydrogen-bonded network. In particular, the two methanol molecules of a 7AI(MeOH)(2) cluster assist the excited-state triple-proton transfer (ESTPT) reaction effectively with highest probability of PT.     

Electronic excited-state energies from a linear response theory based on the ground-state two-electron reduced density matrix

Ground-state two-particle reduced density matrices (2-RDMs) are used to calculate excited-state energy spectra. Solving the Schrodinger equation for excited states dominated by single excitations from the ground-state wavefunction requires the ground-state 2- and 3-RDMs. The excited states, however, can be obtained without a knowledge of the ground-state 3-RDM by two methods: (i) cumulant expansion methods which build the 3-RDM from the 2-RDM, and (ii) double commutator methods which eliminate the 3-RDM. Previous work [Mazziotti, Phys. Rev. A 68, 052501 (2003)] examined the accuracy of excited states extracted from ground-state 2-RDMs, which were calculated by full configuration interaction or the variational 2-RDM method. In this work we employ (i) advances in semidefinite programming to treat the excited states of water and hydrogen fluoride and chains of hydrogen atoms, and (ii) the addition of partial three-particle N-representability conditions to compute more accurate ground-state 2-RDMs. With the hydrogen chains we examine the metal-to-insulator transition as measured by the band gap (the difference between the ground-state and the first excited-state energies), which is difficult for excited-state methods to capture.     

Ultrafast pump-probe study of excited-state charge-transfer dynamics in umecyanin from horseradish root

We have applied femtosecond pump-probe spectroscopy to investigate the excited-state dynamics of umecyanin from horseradish roots, by exciting its 600-nm ligand-to-metal charge-transfer band with a 15-fs pulse and probing over a broad range in the visible region. The decay of the pump-induced ground-state bleaching is modulated by clearly visible oscillations and occurs exponentially with a time constant depending on the observed spectral component of the transmission difference signal, ranging from 270 fs up to 700 fs. The slower decaying process characterizes the spectral component corresponding to the metal-to-ligand charge-transfer transition. The excited-state decay rate is significantly lower than in other blue copper proteins, probably because of the larger energy gap between ligand- and metal-based orbitals in umecyanin. Wavelength dependence of the recovery times could be due to either the excitation of several transitions or the occurrence of intramolecular vibrational relaxation within the excited state. We also find evidence of a hot ground-state absorption, at 700 nm, persisting for several picoseconds. The vibrational coherence induced by the ultrashort pump pulse allows vibrational activity to be observed, mainly in the ground state, as expected in a system with fast excited-state decay. However, we find evidence of a rapidly damped oscillation, which we assign to the excited state. Finally, the Fourier transform of the oscillatory component of the signal presents additional bands in the low-frequency region which are assigned to collective motions of the protein.     

Watching ultrafast barrierless excited-state isomerization of pseudocyanine in real time

The photoinduced excited-state processes in 1,1'-diethyl-2,2'-cyanine iodine are investigated using femtosecond time-resolved pump-probe spectroscopy. Using a broad range of probe wavelengths, the relaxation of the initially prepared excited-state wavepacket can be followed down to the sink region. The data directly visualize the directed downhill motion along the torsional reaction coordinate and suggest a barrierless excited-state isomerization in the short chain cyanine dye. Additionally, ultrafast ground-state hole and excited-state hole replica broadening is observed. While the narrow excited-state wavepacket broadens during pump-probe overlap, the ground-state hole burning dynamics takes place on a significantly longer time-scale. The experiment reported can be considered as a direct monitoring of the shape and the position of the photoprepared wavepacket on the excited-state potential energy surface.     

Coherent vibrational dynamics of Au144(SR)60 nanoclusters

The coherent vibrational dynamics of gold nanoclusters (NCs) provides important information on the coupling between vibrations and electrons as well as their mechanical properties, which is critical for understanding the evolution from a metallic state to a molecular state with diminishing size. Coherent vibrations have been widely explored in small-sized atomically precise gold NCs, while it remains a challenge to observe them in large-sized gold NCs. In this work, we report the coherent vibrational dynamics of atomically precise Au₁₄₄(SR)₆₀ NCs via temperature-dependent femtosecond transient absorption (TA) spectroscopy. The population dynamics of Au₁₄₄(SR)₆₀ consists of three relaxation processes: internal conversion, core–shell charge transfer and relaxation to the ground state. After removing the population dynamics from the TA kinetics, fast Fourier transform analysis on the residual oscillation reveals distinct vibrational modes at 1.5 THz (50 cm⁻¹) and 2 THz (67 cm⁻¹), which arise from the wavepacket motions along the ground-state and excited-state potential energy surfaces (PES), respectively. These results are helpful for understanding the physical properties of gold nanostructures with a threshold size that lies in between those of molecular-like NCs and metallic-state nanoparticles.

The coherent vibrational dynamics of Au₁₄₄(SR)₆₀ nanoclusters was revealed by temperature-dependent ultrafast transient absorption spectroscopy. Both excited-state and ground-state wavepacket motions contribute to the vibrational coherence.     

Second-chance forward isomerization dynamics of the red/green cyanobacteriochrome NpR6012g4 from Nostoc punctiforme

The primary ultrafast Z-to-E isomerization photodynamics of the phytochrome-related cyanobacteriochrome NpR6012g4 from Nostoc punctiforme was studied by transient absorption pump-dump-probe spectroscopy. A 2 ps dump pulse resonant with the stimulated emission band depleted 21% of the excited-state population, while the initial photoproduct Lumi-R was depleted by only 11%. We observed a red-shifted ground-state intermediate (GSI) that we assign to a metastable state that failed to isomerize fully. Multicomponent global analysis implicates the generation of additional Lumi-R from the GSI via crossing over the ground-state thermal barrier for full isomerization, explaining the discrepancy between excited-state and Lumi-R depletion by the dump pulse. This second-chance ground-state dynamics provides a plausible explanation for the unusually high quantum yield of 40% for the primary isomerization step in the forward reaction of NpR6012g4.     

Experimental studies of light-induced charge transfer and charge redistribution in (X(2)-bipyridine)Re(I)(CO)(3)Cl complexes

Stark emission spectroscopy, transient DC photoconductivity (TDCP), and ground-state dipole moment measurements have been used to evaluate charge transfer (CT) within various (X(2)-bipyridine)Re(I)(CO)(3)Cl complexes following (3)MLCT excited-state formation. The Stark technique reports on vector differences between ground-state (mu(g)) and excited-state (mu(e)) dipole moments, while TDCP, when combined with independently obtained mu(g) information, reports on scalar differences. For systems featuring collinear, same-signed ground- and excited-state dipole moments, the scalar and vector differences are equivalent. However, for the low symmetry systems studied here, they are distinctly different. The vector difference yields the effective adiabatic one-electron-transfer distance (R(12)), while the combined vector and scalar data yield information about dipole rotation upon ground-state/excited-state interconversion. For the systems examined, charge transfer distances are substantially smaller than geometric electron-donor/electron-acceptor site separation distances. The measured distances are significantly affected by changes in acceptor ligand substituent composition. Electron-donating substituents decrease CT distances, while electron-withdrawing substituents increase CT distances, as do aromatic substituents that are capable of expanding the bipyridyl ligand (acceptor ligand) pi system. The Stark measurements additionally indicate that the CT vector and the transition dipole moment are significantly orthogonal, a consequence of strong polarization of the Re-Cl bond (orthogonal to the metal/acceptor-ligand plane) in the ground electronic state and relaxation of the polarization in the upper state. The ground-state Re-Cl bond polarization is sufficiently large that the overall ground-state scalar dipole moment exceeds the overall excited-state scalar dipole moment, despite transfer of an electron from the metal center to the diimine ligand. This finding provides an explanation for the otherwise puzzling negative solvatochromism exhibited in this family of compounds. Combining TDCP and Stark results, we find that the dipole moment can be rotated in some instances by more than 90 degrees upon (3)MLCT excited-state formation. The degree of rotation or reorientation can be modulated by changing the identity of the acceptor ligand substituents. Reorientational effects are smallest when the compounds feature aromatic substituents capable of spatially extending the pi system of the acceptor ligand.     

Photodissociation of HBr. 1. Electronic structure, photodissociation dynamics, and vector correlation coefficients

Ab initio potential energy curves, transition dipole moments, and spin-orbit coupling matrix elements are computed for HBr. These are then used, within the framework of time-dependent quantum-mechanical wave-packet calculations, to study the photodissociation dynamics of the molecule. Total and partial integral cross sections, the branching fraction for the formation of excited-state bromine atoms Br(2P(1/2)), and the lowest order anisotropy parameters, beta, for both ground and excited-state bromine are calculated as a function of photolysis energy and compared to experimental and theoretical data determined previously. Higher order anisotropy parameters are computed for the first time for HBr and compared to recent experimental measurements. A new expression for the Re[a1(3) (parallel, perpendicular)] parameter describing coherent parallel and perpendicular production of ground-state bromine in terms of the dynamical functions is given. Although good agreement is obtained between the theoretical predictions and the experimental measurements, the discrepancies are analyzed to establish how improvements might be achieved. Insight is obtained into the nonadiabatic dynamics by comparing the results of diabatic and fully adiabatic calculations.     

Generalized two-dimensional heterocorrelation analysis of spectrally resolved and temporally resolved fluorescence of the 8-anilino-1-naphthalenesulfonate-apomyoglobin complex with pH perturbation

We demonstrate two-dimensional heterocorrelation analysis between spectrally resolved and temporally resolved fluorescence to investigate the decay dynamics of the 8-anilino-1-naphthalenesulfonate- (ANS-) apomyoglobin complex. The dynamic changes of the lifetime components are disclosed across the emission spectrum with an external pH-perturbation. Two different fluorescence lifetime schemes of the ANS-apomyoglobin complex are revealed. From pH 8.5 to 4.5, the transition of protein conformation from the native state to the folding intermediate, a short lifetime component is found to correlate with a short-wavelength emission whose population diminishes with decreasing pH. The lifetime components reflect the excited-state populations of the nascent and the charge-transfer species. From pH 4.2 to 1.0, the transition from the folding intermediate to the acid-unfolded state, the short lifetime is responsible for a long-wavelength emission and the fraction of this component increases when the solution becomes more acidic. In this pH range, the decay components reflect the ground-state populations of microenvironments. The relative decay dynamics across the emission spectrum are revealed without collecting decays at each wavelength. More importantly, these conclusions are reached without the necessity of statistical fitting of the decay data with an a priori decay model.     

A comparison between different semiclassical approximations for optical response functions in nonpolar liquid solutions

The temporal behavior of optical response functions (ORFs) reflects the quantum dynamics of an electronic superposition state, and as such lacks a well-defined classical limit. In this paper, we consider the importance of accounting for the quantum nature of the dynamics when calculating ORFs of different types. To this end, we calculated the ORFs associated with the linear absorption spectrum and the nonlinear two-pulse photon-echo experiment, via the following approaches: (1) the semiclassical forward-backward approach; (2) an approach based on linearizing the path-integral forward-backward action in terms of the difference between the forward and backward paths; (3) an approach based on ground state nuclear dynamics. The calculations were performed on a model that consists of a two-state chromophore solvated in a nonpolar liquid. The different methods were found to yield very similar results for the absorption spectrum and "diagonal" two-pulse photon echo (i.e., the homodyne-detected signal at time t=t(0) after the second pulse, where t(0) is the time interval between the two pulses). The different approximations yielded somewhat different results in the case of the time-integrated photon-echo signal. The reasons for the similarity between the predictions of different approximations are also discussed.     

Excited-state dynamics of carotenoids in light-harvesting complexes. 1. Exploring the relationship between the S1 and S* states

Dispersed transient absorption spectra collected at variable excitation intensities in combination with time-resolved signals were used to explore the underlying connectivity of the electronic excited-state manifold of the carotenoid rhodopin glucoside in the light-harvesting 2 complex isolated from Rhodopseudomonas acidophila. We find that the S state, which was recently identified as an excited state in carotenoids bound in bacterial light-harvesting complexes, exhibits a different response to the increase of excitation intensity than the S(1) state, which suggests that the models used so far to describe the excited states of carotenoids are incomplete. We propose two new models that can describe both the time-resolved and the intensity-dependent data; the first postulates that S(1) and S* are not populated in parallel after the decay of the initially excited S(2) state but instead result from the excitation of distinct ground-state subpopulations. The second model introduces a resonantly enhanced light-induced transition during excitation, which promotes population to higher-lying excited states that favors the formation of S* over S(1). Multiwavelength target analysis of the time-resolved and excitation-intensity dependence measurements were used to characterize the involved states and their responses. We show that both proposed models adequately fit the measured data, although it is not possible to determine which model is most apt. The physical origins and implications of both models are explored.     

Photodissociation of HI and DI: testing models for electronic structure via polarization of atomic photofragments

The photodissociation dynamics of HI and DI are examined using time-dependent wave-packet techniques. The orientation and alignment parameters aQ(K) (p) are determined as a function of photolysis energy for the resulting ground-state I(2P(3/2)) and excited-state I(2P(1/2)) atoms. The aQ(K) (p) parameters describe the coherent and incoherent contributions to the angular momentum distributions from the A 1pi(1), a 3pi(1), and t 3sigma(1) electronic states accessed by perpendicular excitation and the a 3pi(0+) state accessed by a parallel transition. The outcomes of the dynamics based on both shifted ab initio results and three empirical models for the potential-energy curves and transition dipole moments are compared and contrasted. It is demonstrated that experimental measurement of the aQ(K) (p) parameters for the excitation from the vibrational ground state (upsilon=0) would be able to distinguish between the available models for the HI potential-energy curves and transition dipole moments. The differences between the aQ(K) (p) parameters for the excitation from upsilon=0 stand in sharp contrast to the scalar properties, i.e., total cross section and I* branching fraction, which require experimental measurement of photodissociation from excited vibrational states (upsilon>0) to distinguish between the models.     

The effect of hydrogen bonding on the excited-state proton transfer in 2-(2'-hydroxyphenyl)benzothiazole: a TDDFT molecular dynamics study

The dynamics of the excited-state proton transfer (ESPT) in a cluster of 2-(2'-hydroxyphenyl)benzothiazole (HBT) and hydrogen-bonded water molecules was investigated by means of quantum chemical simulations. Two different enol ground-state structures of HBT interacting with the water cluster were chosen as initial structures for the excited-state dynamics: (i) an intramolecular hydrogen-bonded structure of HBT and (ii) a cluster where the intramolecular hydrogen bond in HBT is broken by intermolecular interactions with water molecules. On-the-fly dynamics simulations using time-dependent density functional theory show that after photoexcitation to the S(1) state the ESPT pathway leading to the keto form strongly depends on the initial ground state structure of the HBT-water cluster. In the intramolecular hydrogen-bonded structures direct excited-state proton transfer is observed within 18 fs, which is a factor two faster than proton transfer in HBT computed for the gas phase. Intermolecular bonded HBT complexes show a complex pattern of excited-state proton transfer involving several distinct mechanisms. In the main process the tautomerization proceeds via a triple proton transfer through the water network with an average proton transfer time of approximately 120 fs. Due to the lack of the stabilizing hydrogen bond, intermolecular hydrogen-bonded structures have a significant degree of interring twisting already in the ground state. During the excited state dynamics, the twist tends to quickly increase indicating that internal conversion to the electronic ground state should take place at the sub-picosecond scale.     

Ultrafast dynamics of isolated model photoactive yellow protein chromophores: "Chemical perturbation theory" in the laboratory

Pump-probe and pump-dump probe experiments have been performed on several isolated model chromophores of the photoactive yellow protein (PYP). The observed transient absorption spectra are discussed in terms of the spectral signatures ascribed to solvation, excited-state twisting, and vibrational relaxation. It is observed that the protonation state has a profound effect on the excited-state lifetime of p-coumaric acid. Pigments with ester groups on the coumaryl tail end and charged phenolic moieties show dynamics that are significantly different from those of other pigments. Here, an unrelaxed ground-state intermediate could be observed in pump-probe signals. A similar intermediate could be identified in the sinapinic acid and in isomerization-locked chromophores by means of pump-dump probe spectroscopy; however, in these compounds it is less pronounced and could be due to ground-state solvation and/or vibrational relaxation. Because of strong protonation-state dependencies and the effect of electron donor groups, it is argued that charge redistribution upon excitation determines the twisting reaction pathway, possibly through interaction with the environment. It is suggested that the same pathway may be responsible for the initiation of the photocycle in native PYP.     

Analytical approach for the excited-state Hessian in time-dependent density functional theory: formalism, implementation, and performance

The paper presents the formalism, implementation, and performance of the analytical approach for the excited-state Hessian in the time-dependent density functional theory (TDDFT) that extends our previous work [J. Liu and W. Z. Liang, J. Chem. Phys. 135, 014113 (2011)] on the analytical Hessian in TDDFT within Tamm-Dancoff approximation (TDA) to full TDDFT. In contrast to TDA-TDDFT, an appreciable advantage of full TDDFT is that it maintains the oscillator strength sum rule, and therefore yields more precise results for the oscillator strength and other related physical quantities. For the excited-state harmonic vibrational frequency calculation, however, full TDDFT does not seem to be advantageous since the numerical tests demonstrate that the accuracy of TDDFT with and without TDA are comparable to each other. As a common practice, the computed harmonic vibrational frequencies are scaled by a suitable scale factor to yield good agreement with the experimental fundamental frequencies. Here we apply both the optimized ground-state and excited-state scale factors to scale the calculated excited-state harmonic frequencies and find that the scaling decreases the root-mean-square errors. The optimized scale factors derived from the excited-state calculations are slightly smaller than those from the ground-state calculations.