Time-Dependent Approach to Resonance Raman Spectra Including Duschinsky Rotation and Herzberg–Teller Effects: Formalism and Its Realistic Applications

Time-dependent surface-enhanced Raman scattering: A theoretical approach

A new procedure for computing the time-dependent Raman scattering of molecules in the proximity of plasmonic nanoparticles (NPs) is proposed, drawing inspiration from the pioneering Lee and Heller's theory. This strategy is based on a preliminary simulation of the molecular vibronic wavefunction in the presence of a plasmonic nanostructure and an incident light pulse. Subsequently, the Raman signal is evaluated through an inverse Fourier Transform of the coefficients' dynamics. Employing a multiscale approach, the system is treated by coupling the quantum mechanical description of the molecule with the polarizable continuum model for the NP. This method offers a unique advantage by providing insights into the time evolution of the plasmon-enhanced Raman signal, tracking the dynamics of the incident electric field. It not only provides for the total Raman signal at the process's conclusion but also gives transient information. Importantly, the flexibility of this approach allows for the simulation of various incident electric field profiles, enabling a closer alignment with experimental setups. This adaptability ensures that the method is relevant and applicable to diverse real-world scenarios.

A theoretical investigation of benzothiadiazole derivatives for high efficiency OLEDs

It is of great importance and worthy of efforts to give a clear structure-property relationship and microscopic mechanism of fluorescence emitters with high quantum yield. In this work, we perform a detailed computational investigation to give an explanation to the high efficiency of a fluorescence emitter XBTD-NPh based TADF sensitized fluorescence (TSF) OLEDs, and construct a symmetry structure DSBNA-BTD. Theoretical calculations show that XBTD-NPh is a long-time phosphorescent material at 77 K and TADF is attributed to the RISC of T1 to S1 state. For DSBNA-BTD, excitons arrived at T1 state comes to a large rate of nonradiatively path to the ground state, meaning it is may not be an efficient TADF molecule. For both molecules, the fast IC between T2 and T1 state results in that the hot exciton channel T1-Tn-S1 makes no contribution to the TADF.

Time-Dependent Resonant Inelastic X-ray Scattering of Pyrazine at the Nitrogen K-Edge: A Quantum Dynamics Approach

We calculate resonant inelastic X-ray scattering spectra of pyrazine at the nitrogen K-edge in the time domain including wavepacket dynamics in both the valence and core-excited state manifolds. Upon resonant excitation, we observe ultrafast non-adiabatic population transfer between core-excited states within the core-hole lifetime, leading to molecular symmetry distortions. Importantly, our time-domain approach inherently contains the ability to manipulate the dynamics of this process by detuning the excitation energy, which effectively shortens the scattering duration. We also explore the impact of pulsed incident X-ray radiation, which provides a foundation for state-of-the-art time-resolved experiments with coherent pulsed light sources.

The Resonance Raman Spectrum of Cytosine in Water: Analysis of the Effect of Specific Solute-Solvent Interactions and Non-Adiabatic Couplings

In this contribution, we report a computational study of the vibrational Resonance Raman (vRR) spectra of cytosine in water, on the grounds of potential energy surfaces (PES) computed by time-dependent density functional theory (TD-DFT) and CAM-B3LYP and PBE0 functionals. Cytosine is interesting because it is characterized by several close-lying and coupled electronic states, challenging the approach commonly used to compute the vRR for systems where the excitation frequency is in quasi-resonance with a single state. We adopt two recently developed time-dependent approaches, based either on quantum dynamical numerical propagations of vibronic wavepackets on coupled PES or on analytical correlation functions for cases in which inter-state couplings were neglected. In this way, we compute the vRR spectra, considering the quasi-resonance with the eight lowest-energy excited states, disentangling the role of their inter-state couplings from the mere interference of their different contributions to the transition polarizability. We show that these effects are only moderate in the excitation energy range explored by experiments, where the spectral patterns can be rationalized from the simple analysis of displacements of the equilibrium positions along the different states. Conversely, at higher energies, interference and inter-state couplings play a major role, and the adoption of a fully non-adiabatic approach is strongly recommended. We also investigate the effect of specific solute-solvent interactions on the vRR spectra, by considering a cluster of cytosine, hydrogen-bonded by six water molecules, and embedded in a polarizable continuum. We show that their inclusion remarkably improves the agreement with the experiments, mainly altering the composition of the normal modes, in terms of internal valence coordinates. We also document cases, mostly for low-frequency modes, in which a cluster model is not sufficient, and more elaborate mixed quantum classical approaches, in explicit solvent models, need to be applied.

FCclasses3: Vibrationally-resolved spectra simulated at the edge of the harmonic approximation

We introduce FCclasses3, a code to carry out vibronic simulations of electronic spectra and nonradiative rates, based on the harmonic approximation. Key new features are: implementation of the full family of vertical and adiabatic harmonic models, vibrational analysis in curvilinear coordinates, extension to several electronic spectroscopies and implementation of time-dependent approaches. The use of curvilinear valence internal coordinates allows the adoption of quadratic model potential energy surfaces (PES) of the initial and final states expanded at arbitrary configurations. Moreover, the implementation of suitable projectors provides a robust framework for defining reduced-dimensionality models by sorting flexible coordinates out of the harmonic subset, so that they can then be treated at anharmonic level, or with mixed quantum classical approaches. A set of tools to facilitate input preparation and output analysis is also provided. We show the program at work in the simulation of different spectra (one and two-photon absorption, emission and resonance Raman) and internal conversion rate of a typical rigid molecule, anthracene. Then, we focus on absorption and emission spectra of a series of flexible polyphenyl molecules, highlighting the relevance of some of the newly implemented features. The code is freely available at http://www.iccom.cnr.it/en/fcclasses/.

Time Resolved Raman Scattering of Molecules: A Quantum Mechanics Approach with Stochastic Schroedinger Equation

Raman scattering is a very powerful tool employed to characterize molecular systems. Here we propose a novel theoretical strategy to calculate the Raman cross-section in time domain, by computing the cumulative Raman signal emitted during the molecular evolution in time. Our model is based on a numerical propagation of the vibronic wave function under the effect of a light pulse of arbitrary shape. This approach can therefore tackle a variety of experimental setups. Both resonance and nonresonance Raman scattering can be retrieved, and also the time-dependent fluorescence emission is computed. The model has been applied to porphyrin considering both resonance and nonresonance conditions and varying the incident field duration. Moreover the effect of the vibrational relaxation, which should be taken into account when its time scale is similar to that of the Raman emission, has been included through the stochastic Schroedinger equation approach.

Nonadiabatic Vibrational Resonance Raman Spectra from Quantum Dynamics Propagations with LVC Models. Application to Thymine

We present a viable protocol to compute vibrational resonance Raman (vRR) spectra for systems with several close-lying and potentially coupled electronic states. It is based on the parametrization of linear vibronic coupling (LVC) models from time-dependent density functional theory (TD-DFT) calculations and quantum dynamics propagations of vibronic wavepackets with the multilayer version of the multiconfiguration time-dependent Hartree (ML-MCTDH) method. Our approach is applied to thymine considering seven coupled electronic states, comprising the three lowest bright states, and all vibrational coordinates. Computed vRR at different excitation wavelengths are in good agreement with the available experimental data. Up to 250 nm the signal is dominated by the lowest HOMO → LUMO transition, whereas at 233 nm, in the valley between the two lowest energy absorption bands, the contributions of all the three bright states, and their interferences and couplings, are important. Inclusion of solvent (water) effects improves the agreement with experiment, reproducing the coalescence of vibrational bands due to CC and C═O stretchings. With our approach we disentangle and assess the effect of interferences between the contribution of different quasi-resonant states to the transition polarizability and the effect of interstate couplings. Our findings strongly suggest that in cases of close-lying and potentially coupled states a simple inclusion of interference effects is not sufficient, and a fully nonadiabatic computation should instead be performed. We also document that for systems with strong couplings and quasi-degenerate states, the use of HT perturbative approach, not designed for these cases, may lead to large artifacts.

Efficient simulation of resonance Raman spectra with tight-binding approximations to Density Functional Theory

Resonance Raman spectroscopy has long been established as one of the most sensitive techniques for detection, structure characterization and probing the excited-state dynamics of biochemical systems. However, the analysis of resonance Raman spectra is much facilitated when measurements are accompanied by Density Functional Theory (DFT) calculations which are expensive for large biomolecules. In this work, resonance Raman spectra are therefore computed with the Density Functional Tight-Binding (DFTB) method in the time-dependent excited-state gradient approximation. To test the accuracy of the tight-binding approximations, this method is first applied to typical resonance Raman benchmark molecules like β-carotene and compared to results obtained with pure and range-separated exchange-correlation (xc) functionals. We then demonstrate the efficiency of the approach by considering a computationally challenging heme variation. Overall, we find that the vibrational frequencies and excited-state properties (energies and gradients) which are needed to simulate the spectra are reasonably accurate and suitable for interpretation of experiments. We can therefore recommend DFTB as a fast computational method to interpret resonance Raman spectra.

Vibronic Coupling Effect on the Vibrationally Resolved Electronic Spectra and Intersystem Crossing Rates of a TADF Emitter: 7-PhQAD

Assessing and improving the performance of organic light-emitting diode (OLED) materials require quantitative prediction of rate coefficients for the intersystem crossing (ISC) and reverse ISC (RISC) processes, which are determined not only by the energy gap and the direct spin-orbit coupling (SOC) between the first singlet and triplet excited-states at a thermal equilibrium position of the initial electronic state but also by the non-Condon effects such as the Herzberg-Teller-like vibronic coupling (HTVC) and the spin-vibronic coupling (SVC). Here we apply the time-dependent correlation function approaches to quantitatively calculate the vibrationally resolved absorption and fluorescence spectra and ISC/RISC rates of a newly synthesized multiple-resonance-type (MR-type) thermally activated delayed fluorescence (TADF) emitter, 7-phenylquinolino[3,2,1-de]acridine-5,9-dione (7-PhQAD), with the inclusion of the Franck-Condon (FC), HTVC, and Duschinsky rotation (DR) effects. The SVC effect on the rates has also been approximately evaluated. We find that the experimentally measured ISC rates of 7-PhQAD originate predominantly from the vibronic coupling, consistent with the previous reports on other MR-type TADF emitters. The SVC effect on ISC rates is about 10 times larger than the HTVC effect, and the latter increases the ISC rates by more than 1 order of magnitude while it slightly affects the vibrationally resolved absorption and fluorescence spectra. The discrepancy between the theoretical and experimental results is attributed to inaccurately describing excited-states calculated by the time-dependent density functional theory as well as to not fully accounting for the complex experimental conditions. This work provides a demonstration of what proportion of ISC and RISC rate coefficients of a MR-type TADF emitter can be covered by the HTVC effect, and it opens design routes that go beyond the FC approximation for the future development of high-performance OLED devices.

Ultrafast stimulated resonance Raman signatures of lithium polysulfides for shuttling effect characterization: An ab initio study

The shuttling effect is a crucial obstacle to the practical deployment of lithium sulfur batteries (LSBs). This can be ascribed to the generation of lithium polysulfide (LiPS) redox intermediates that are soluble in the electrolyte. The detailed mechanism of the shuttling, including the chemical structures responsible for the loss of effective mass and the dynamics/kinetics of the redox reactions, are not clear so far. To obtain this microscopic information, characterization techniques with high spatial and temporal resolutions are required. Here, we propose that resonance Raman spectroscopy combined with ultrafast broadband pulses is a powerful tool to reveal the mechanism of the shuttling effect. By combining the chemical bond level spatial resolution of resonance Raman and the femtosecond scale temporal resolution of the ultrafast pulses, this novel technique holds the potential of capturing the spectroscopic fingerprints of the LiPS intermediates during the working stages of LSBs. Using ab initio simulations, we show that, in addition to the excitation energy selective enhancement, resonance Raman signals of different LiPS intermediates are also characteristic and distinguishable. These results will facilitate the real-time in situ monitoring of LiPS species and reveal the underlying mechanism of the shuttling effect.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Sum-over-state expressions including second-order Herzberg-Teller effects for the calculation of absorption and resonance Raman intensities

The sum-over-state expressions are derived to calculate the second-order Herzberg-Teller (HT) effects in absorption and resonance Raman spectroscopies. These effects depend on the second derivatives of the transition dipole moment with respect to the vibrational coordinates. The method is applied to the molecule of 1,3-butadiene using density functional theory calculations. It is found that the second-order HT effects are significant for both absorption and resonance Raman intensities, and that the calculated spectra are in good agreement with the experimental data. The second-order HT effects originate from diagonal elements of the second derivatives matrix, whereas non-diagonal elements have a negligible impact on the intensities of 1,3-butadiene.

Herzberg-Teller Effect Predominates in Sum-Frequency Vibrational Spectroscopy of Limonene Chiral Liquids

We theoretically study the bulk sum-frequency vibrational spectroscopy of chiral liquids under the influence of the Franck-Condon, Herzberg-Teller, and nonadiabatic effects. With quantum chemistry computations we calculate the chiral spectra for the R-limonene molecule. When we compare the theoretical and experimental spectra, we find that the Herzberg-Teller effect under the Born-Oppenheimer approximation, instead of the nonadiabatic effect, predominates in the chiral spectra.

Relative contributions of Franck-Condon to Herzberg-Teller terms in charge transfer surface-enhanced Raman scattering spectroscopy

We have theoretically modeled charge transfer (CT) surface enhanced raman scattering (SERS) spectroscopy using pyridine bound to a planar Ag₆ metal nanocluster. CT states were determined by natural transition orbital hole-particle plots and CT distance D_CT and the amount of charge transferred q_CT indices. We first consider a resonance Raman (RR) model based on the Albrecht approach and calculate the ratio of the Herzberg-Teller (HT) B or C term to the Franck-Condon (FC) A term for a totally symmetric a₁ vibrational mode exciting in the lowest energy CT state. Using a dimensionless upper limit to the displacement factor ∆ = 0.05 in the FC term based on the examination of overtones in experimental spectra and a calculated HT coupling constant h_CT = 0.439 eV/Å(amu)^1/2 in the HT term, we calculated the scattering ratio of the HT to FC intensities as 147. This example indicated that for totally symmetric modes, the scattering intensity would all come from HT scattering. To further verify this result, we used the general time-dependent-RR formulation of Baiardi, Bloino, and Barone with the adiabatic Hessian model to calculate the FC, the Frank-Condon and Herzberg-Teller (FCHT), and the HT terms for pyridine in the C_2v Ag₆-pyridine complexes. For all cases we studied with pyridine in two orientations either parallel or perpendicular to the planar Ag₆ cluster, the HT terms, FCHT + HT, dominate the FC term in the CT RR spectrum. These results indicate that for CT SERS, the intensity of all the totally and non-totally symmetric vibrational modes should come from the HT effect.

Calculation of vibrationally resolved absorption spectra of acenes and pyrene

Absorption spectra of polycyclic aromatic hydrocarbons have been simulated by using a real-time generating function method that combines adiabatic electronic excitation energies with vibrational energies of the excited states.

Time-Dependent Formulation of Resonance Raman Optical Activity Spectroscopy

In this work, we extend the theoretical framework recently developed for the simulation of resonance Raman (RR) spectra of medium-to-large sized systems to its chiral counterpart, namely, resonance Raman optical activity (RROA). The theory is based on a time-dependent (TD) formulation, with the transition tensors obtained as half-Fourier transforms of the appropriate cross-correlation functions. The implementation has been kept as general as possible, supporting adiabatic and vertical models for the PES representation, both in Cartesian and internal coordinates, with the possible inclusion of Herzberg-Teller (HT) effects. Thanks to the integration of this TD-RROA procedure within a general-purpose quantum-chemistry program, both solvation and leading anharmonicity effects can be included in an effective way. The implementation is validated on one of the smallest chiral molecule (methyloxirane). Practical applications are illustrated with three medium-size organic molecules (naproxen-OCD₃, quinidine and 2-Br-hexahelicene), whose simulated spectra are compared to the corresponding experimental data.

Probing Dynamics in Higher-Lying Electronic States with Resonance-Enhanced Femtosecond Stimulated Raman Spectroscopy

Femtosecond stimulated Raman scattering (FSRS) measurements typically probe the structural dynamics of a molecule in the first electronically excited state, S₁. While these measurements often rely on an electronic resonance condition to increase signal strength or enhance species selectivity, the effects of the resonance condition are usually neglected. However, mode-specific enhancements of the vibrational transitions in an FSRS spectrum contain detailed information about the resonant (upper) electronic state. Analogous to ground-state resonance Raman spectroscopy, the relative intensities of the Raman bands reveal displacements of the upper potential energy surface due to changes in the bonding pattern upon S _n ← S₁ electronic excitation, and therefore provide a sensitive probe of the ultrafast dynamics in the higher-lying state, S _n. Raman gain profiles from the wavelength-dependent FSRS spectrum of the model compound 2,5-diphenylthiophene (DPT) reveal several modes with large displacement in the upper potential energy surface, including strong enhancement of a delocalized C-S-C stretching and ring deformation mode. The experimental results provide a benchmark for comparison with calculated spectra using time-dependent density functional theory (TD-DFT) and equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD), where the calculations are based on the time-dependent formalism for resonance Raman spectroscopy. The simulated spectra are obtained from S₁-S _n transition strengths and the energy gradients of the upper (S _n) potential energy surfaces along the S₁ normal mode coordinates. The experimental results provide a stringent test of the computational approach, and indicate important limitations based on the level of theory and basis set. This work provides a foundation for making more accurate assignments of resonance-enhanced excited-state Raman spectra, as well as extracting novel information about higher-lying excited states in the transient absorption spectrum of a molecule.

Theories of quantum dissipation and nonlinear coupling bath descriptors

The quest of an exact and nonperturbative treatment of quantum dissipation in nonlinear coupling environments remains in general an intractable task. In this work, we address the key issues toward the solutions to the lowest nonlinear environment, a harmonic bath coupled both linearly and quadratically with an arbitrary system. To determine the bath coupling descriptors, we propose a physical mapping scheme, together with the prescription reference invariance requirement. We then adopt a recently developed dissipaton equation of motion theory [R. X. Xu et al., Chin. J. Chem. Phys. 30, 395 (2017)], with the underlying statistical quasi-particle ("dissipaton") algebra being extended to the quadratic bath coupling. We report the numerical results on a two-level system dynamics and absorption and emission line shapes.

The role of Herzberg-Teller effects on the resonance Raman spectrum of trans-porphycene investigated by time dependent density functional theory

The S₁ excited state properties as well as the associated absorption and resonance Raman (RR) spectra of trans-porphycene are investigated by means of time dependent density functional theory calculations. The relative magnitude of the Franck-Condon (FC) contribution and of the Herzberg-Teller (HT) effects is evaluated for both the absorption and RR intensities. The accuracy of the calculated spectra is assessed by employing different theoretical approximations and by comparing with experimental data. The obtained results show that Duschinsky effects lead to noticeable modifications in the absorption intensities but are nearly negligible in the RR spectrum. By contrast, the HT effects are stronger for the RR intensities compared to the absorption intensities, and these effects significantly improve the agreement with the experimental RR spectrum. Moreover, the HT effects produce different values of the RR depolarization ratios, which can be used to quantify the relative importance of the FC and HT contributions. Generally, it is found that the HT effects have a significant role on the RR spectrum of trans-porphycene and that their inclusion in the computational scheme is mandatory to accurately predict the RR intensities.

Identifying Cu(ii)-amyloid peptide binding intermediates in the early stages of aggregation by resonance Raman spectroscopy: a simulation study

The aggregation of amyloid beta (Aβ) peptides plays a crucial role in the pathology and etiology of Alzheimer's disease. Experimental evidence shows that copper ion is an aggregation-prone species with the ability to coordinately bind to Aβ and further induce the formation of neurotoxic Aβ oligomers. However, the detailed structures of Cu(ii)-Aβ complexes have not been illustrated, and the kinetics and dynamics of the Cu(ii) binding are not well understood. Two Cu(ii)-Aβ complexes have been proposed to exist under physiological conditions, and another two might exist at higher pH values. By using ab initio simulations for the spontaneous resonance Raman and time domain stimulated resonance Raman spectroscopy signals, we obtained the characteristic Raman vibronic features of each complex. These signals contain rich structural information with high temporal resolution, enabling the characterization of transient states during the fast Cu-Aβ binding and interconversion processes.

Theoretical Assessment of Excited State Gradients and Resonance Raman Intensities for the Azobenzene Molecule

The ground state geometries and vibrational frequencies as well as the excitation energies and excited state gradients of the S₁(nπ*) and S₂(ππ*) states of trans- and cis-azobenzene are investigated by several DFT methods, namely B3LYP, PBE, M06-2X, CAM-B3LYP, and ωB97X. Excited state properties and in particular gradients are also assessed using the wave function based methods EOM-CCSD and RASPT2/RASSCF. Comparison with experimental data shows that the B3LYP functional gives the most accurate results for the ground state geometry and vibrational frequencies. The analysis of the vertical excitation energies reveals that the RASPT2 approach provides the most accurate excitation energies with deviations of the order of 0.1 eV. Among the TDDFT methods, the CAM-B3LYP functional shows the best performance on the excitation energies. By assessing the excited state gradients with respect to the reference RASPT2 data, the most accurate gradients are obtained with B3LYP, whereas other functionals as well as the EOM-CCSD and RASSCF calculations give less consistent results. Overall, despite the tendency of B3LYP to underestimate the excitation energies, this functional provides the most balanced description of both ground and excited state properties for both isomers of azobenzene in the Franck-Condon region.

Organic light-emitting diodes: theoretical understanding of highly efficient materials and development of computational methodology

Theoretical understanding of organic light-emitting diodes started from the quest to the nature of the primary excitation in organic molecular and polymeric materials. We found the electron correlation strength, bond-length alternation as well as the conjugation extent have strong influences on the orderings of the lowest lying excited states through the first application of density matrix renormalization group theory to quantum chemistry. The electro-injected free carriers (with spin 1/2) can form both singlet and triplet bound states. We found that the singlet exciton formation ratio can exceed the conventional 25% spin statistics limit. We proposed a vibration correlation function formalism to evaluate the excited-state decay rates, which is shown to not only give reasonable estimations for the quantum efficiency but also a quantitative account for the aggregation-induced emission (AIE). It is suggested to unravel the AIE mechanism through resonance Raman spectroscopy.

Assessment of Franck-Condon Methods for Computing Vibrationally Broadened UV-vis Absorption Spectra of Flavin Derivatives: Riboflavin, Roseoflavin, and 5-Thioflavin

We address the performance of the vertical and adiabatic Franck-Condon (VFC/AFC) approaches combined with time-independent or time-dependent (TI/TD) formalisms in simulating the one-photon absorption spectra of three flavin compounds with distinct structural features. Calculations were done in the gas phase and in two solvents (water, benzene) for which experimental reference measurements are available. We utilized the independent mode displaced harmonic oscillator model without or with frequency alteration (IMDHO/IMDHO-FA) and also accounted for Duschinsky mixing effects. In the initial validation on the first excited singlet state of riboflavin, the range-separated functionals, CAM-B3LYP and ωB97xD, showed the best performance, but B3LYP also gave a good compromise between peak positions and spectral topology. Large basis sets were not mandatory to obtain high-quality spectra for the selected systems. The presence of a symmetry plane facilitated the computation of vibrationally broadened spectra, since different FC variants yield similar results and the harmonic approximation holds rather well. Compared with the AFC approach, the VFC approach performed equally well or even better for all three flavins while offering several advantages, such as avoiding error-prone geometry optimization procedures on excited-state surfaces. We also explored the advantages of curvilinear displacements and of a Duschinsky treatment for the AFC spectra in cases when a rotatable group is present on the chromophore. Taken together, our findings indicate that the combination of the VFC approach with the TD formalism and the IMDHO-FA model offers the best overall performance.

Vibrationally resolved optical spectra of modified diamondoids obtained from time-dependent correlation function methods

Optical properties of modified diamondoids have been studied theoretically using vibrationally resolved electronic absorption, emission and resonance Raman spectra. A time-dependent correlation function approach has been used for electronic two-state models, comprising a ground state (g) and a bright, excited state (e), the latter determined from linear-response, time-dependent density functional theory (TD-DFT). The harmonic and Condon approximations were adopted. In most cases origin shifts, frequency alteration and Duschinsky rotation in excited states were considered. For other cases where no excited state geometry optimization and normal mode analysis were possible or desired, a short-time approximation was used. The optical properties and spectra have been computed for (i) a set of recently synthesized sp(2)/sp(3) hybrid species with C[double bond, length as m-dash]C double-bond connected saturated diamondoid subunits, (ii) functionalized (mostly by thiol or thione groups) diamondoids and (iii) urotropine and other C-substituted diamondoids. The ultimate goal is to tailor optical and electronic features of diamondoids by electronic blending, functionalization and substitution, based on a molecular-level understanding of the ongoing photophysics.

Effect of TiO2particles on normal and resonance Raman spectra of coumarin 343: a theoretical investigation

Both the normal Raman spectra (NRS) and resonance Raman spectra (RRS) can be used to figure out the isomers and their interfacial structures. Furthermore, the differences in RRS between the locally excited state and the charge transfer state of C343–TiO₂are helpful to understand and control the electron transfer at the interface.

Assessment of mode-mixing and Herzberg-Teller effects on two-photon absorption and resonance hyper-Raman spectra from a time-dependent approach

A time-dependent approach is presented to simulate the two-photon absorption (TPA) and resonance hyper-Raman scattering (RHRS) spectra including Duschinsky rotation (mode-mixing) and Herzberg-Teller (HT) vibronic coupling effects. The computational obstacles for the excited-state geometries, vibrational frequencies, and nuclear derivatives of transition dipole moments, which enter the expressions of TPA and RHRS cross sections, are further overcome by the recently developed analytical excited-state energy derivative approaches in the framework of time-dependent density functional theory. The excited-state potential curvatures are evaluated at different levels of approximation to inspect the effects of frequency differences, mode-mixing and HT on TPA and RHRS spectra. Two types of molecules, one with high symmetry (formaldehyde, p-difluorobenzene, and benzotrifluoride) and the other with non-centrosymmetry (cis-hydroxybenzylidene-2,3-dimethylimidazolinone in the deprotonated anion state (HDBI(-))), are used as test systems. The calculated results reveal that it is crucial to adopt the exact excited-state potential curvatures in the calculations of TPA and RHRS spectra even for the high-symmetric molecules, and that the vertical gradient approximation leads to a large deviation. Furthermore, it is found that the HT contribution is evident in the TPA and RHRS spectra of HDBI(-) although its one- and two-photon transitions are strongly allowed, and its effect results in an obvious blueshift of the TPA maximum with respect to the one-photon absorption maximum. With the HT and solvent effects getting involved, the simulated blueshift of 1291 cm(-1) agrees well with the experimental measurement.

Ultrafast charge transfer in nickel phthalocyanine probed by femtosecond Raman-induced Kerr effect spectroscopy

The recently developed technique of femtosecond stimulated Raman spectroscopy, and its variant, femtosecond Raman-induced Kerr effect spectroscopy (FRIKES), offer access to ultrafast excited-state dynamics via structurally specific vibrational spectra. We have used FRIKES to study the photoexcitation dynamics of nickel(II) phthalocyanine with eight butoxy substituents, NiPc(OBu)8. NiPc(OBu)8 is reported to have a relatively long-lived ligand-to-metal charge-transfer (LMCT) state, an essential characteristic for efficient electron transfer in photocatalysis. Following photoexcitation, vibrational transitions in the FRIKES spectra, assignable to phthalocyanine ring modes, evolve on the femtosecond to picosecond time scales. Correlation of ring core size with the frequency of the ν10 (asymmetric C-N stretching) mode confirms the identity of the LMCT state, which has a ∼500 ps lifetime, as well as that of a precursor d-d excited state. An even earlier (∼0.2 ps) transient is observed and tentatively assigned to a higher-lying Jahn-Teller-active LMCT state. This study illustrates the power of FRIKES spectroscopy in elucidating ultrafast molecular dynamics.

Molecular properties of excited electronic state: formalism, implementation, and applications of analytical second energy derivatives within the framework of the time-dependent density functional theory/molecular mechanics

This work extends our previous works [J. Liu and W. Z. Liang, J. Chem. Phys. 135, 014113 (2011); J. Liu and W. Z. Liang, J. Chem. Phys. 135, 184111 (2011)] on analytical excited-state energy Hessian within the framework of time-dependent density functional theory (TDDFT) to couple with molecular mechanics (MM). The formalism, implementation, and applications of analytical first and second energy derivatives of TDDFT/MM excited state with respect to the nuclear and electric perturbations are presented. Their performances are demonstrated by the calculations of adiabatic excitation energies, and excited-state geometries, harmonic vibrational frequencies, and infrared intensities for a number of benchmark systems. The consistent results with the full quantum mechanical method and other hybrid theoretical methods indicate the reliability of the current numerical implementation of developed algorithms. The computational accuracy and efficiency of the current analytical approach are also checked and the computational efficient strategies are suggested to speed up the calculations of complex systems with many MM degrees of freedom. Finally, we apply the current analytical approach in TDDFT/MM to a realistic system, a red fluorescent protein chromophore together with part of its nearby protein matrix. The calculated results indicate that the rearrangement of the hydrogen bond interactions between the chromophore and the protein matrix is responsible for the large Stokes shift.

A robust and effective time-independent route to the calculation of Resonance Raman spectra of large molecules in condensed phases with the inclusion of Duschinsky, Herzberg-Teller, anharmonic, and environmental effects

We present an effective time-independent implementation to model vibrational resonance Raman (RR) spectra of medium-large molecular systems with the inclusion of Franck-Condon (FC) and Herzberg-Teller (HT) effects and a full account of the possible differences between the harmonic potential energy surfaces of the ground and resonant electronic states. Thanks to a number of algorithmic improvements and very effective parallelization, the full computations of fundamentals, overtones, and combination bands can be routinely performed for large systems possibly involving more than two electronic states. In order to improve the accuracy of the results, an effective inclusion of the leading anharmonic effects is also possible, together with environmental contributions under different solvation regimes. Reduced-dimensionality approaches can further enlarge the range of applications of this new tool. Applications to imidazole, pyrene, and chlorophyll a1 in solution are reported, as well as comparisons with available experimental data.