First principles simulation of the UV absorption spectrum of ethylene using the vertical Franck-Condon approach

Toward the improvement of vibronic spectra and non-radiative rate constants using the vertical Hessian method

For the computation of vibrationally resolved electronic spectra, various approaches can be employed. Adiabatic approaches simulate vibronic transitions using harmonic potentials of the initial and final states, while vertical approaches extrapolate the final state potential from the gradients and Hessian at the Franck-Condon point, avoiding a full exploration of the potential energy surface of the final state. Our implementation of the vertical Hessian (VH) method has been validated with a benchmark set of four small molecules, each presenting unique challenges, such as complex topologies, problematic low-frequency vibrations, or significant geometrical changes upon electronic excitation. We assess the quality of both adiabatic and vertical approaches for simulating vibronic transitions. For two types of donor-acceptor compounds with promising thermally activated delayed fluorescence properties, our computations confirm that the vertical approaches outperform the adiabatic ones. The VH method significantly reduces computational costs and yields meaningful emission spectra, where adiabatic models fail. More importantly, we pioneer the use of the VH method for the computation of rate constants for non-radiative processes, such as intersystem crossing and reverse intersystem crossing along a relaxed interpolated pathway of a donor-acceptor compound. This study highlights the potential of the VH method to advance computational vibronic spectroscopy by providing meaningful simulations of intricate decay pathway mechanisms in complex molecular systems.

Calculating absorption and fluorescence spectra for chromophores in solution with ensemble Franck-Condon methods

Accurately modeling absorption and fluorescence spectra for molecules in solution poses a challenge due to the need to incorporate both vibronic and environmental effects, as well as the necessity of accurate excited state electronic structure calculations. Nuclear ensemble approaches capture explicit environmental effects, Franck-Condon methods capture vibronic effects, and recently introduced ensemble-Franck-Condon approaches combine the advantages of both methods. In this study, we present and analyze simulated absorption and fluorescence spectra generated with combined ensemble-Franck-Condon approaches for three chromophore-solvent systems and compare them to standard ensemble and Franck-Condon spectra, as well as to the experiment. Employing configurations obtained from ground and excited state ab initio molecular dynamics, three combined ensemble-Franck-Condon approaches are directly compared to each other to assess the accuracy and relative computational time. We find that the approach employing an average finite-temperature Franck-Condon line shape generates spectra nearly identical to the direct summation of an ensemble of Franck-Condon spectra at one-fourth of the computational cost. We analyze how the spectral simulation method, as well as the level of electronic structure theory, affects spectral line shapes and associated Stokes shifts for 7-nitrobenz-2-oxa-1,3-diazol-4-yl and Nile red in dimethyl sulfoxide and 7-methoxy coumarin-4-acetic acid in methanol. For the first time, our studies show the capability of combined ensemble-Franck-Condon methods for both absorption and fluorescence spectroscopy and provide a powerful tool for simulating linear optical spectra.

Time dependent vibrational electronic coupled cluster (VECC) theory for non-adiabatic nuclear dynamics

A time-dependent vibrational electronic coupled-cluster (VECC) approach is proposed to simulate photo-electron/UV-VIS absorption spectra as well as time-dependent properties for non-adiabatic vibronic models, going beyond the Born-Oppenheimer approximation. A detailed derivation of the equations of motion and a motivation for the ansatz are presented. The VECC method employs second-quantized bosonic construction operators and a mixed linear and exponential ansatz to form a compact representation of the time-dependent wave-function. Importantly, the method does not require a basis set, has only a few user-defined inputs, and has a classical (polynomial) scaling with respect to the number of degrees of freedom (of the vibronic model), resulting in a favorable computational cost. In benchmark applications to small models and molecules, the VECC method provides accurate results compared to multi-configurational time-dependent Hartree calculations when predicting short-time dynamical properties (i.e., photo-electron/UV-VIS absorption spectra) for non-adiabatic vibronic models. To illustrate the capabilities, the VECC method is also successfully applied to a large vibronic model for hexahelicene with 14 electronic states and 63 normal modes, developed in the group by Aranda and Santoro [J. Chem. Theory Comput. 17, 1691, (2021)].

Prediction of fluorescence quantum yields using the extended thawed Gaussian approximation

Spontaneous emission and internal conversion rates are calculated within harmonic approximations and compared to the results obtained within the semi-classical extended thawed Gaussian approximation (ETGA). This is the first application of the ETGA in the calculation of internal conversion and emission rates for real molecular systems, namely, formaldehyde, fluorobenzene, azulene, and a dicyano-squaraine dye. The viability of the models as black-box tools for prediction of spontaneous emission and internal conversion rates is assessed. All calculations were done using a consistent protocol in order to investigate how different methods perform without previous experimental knowledge using density functional theory (DFT) and time-dependent DFT (TD-DFT) with B3LYP, PBE0, ωB97XD, and CAM-B3LYP functionals. Contrasting the results with experimental data shows that there are further improvements required before theoretical predictions of emission and internal conversion rates can be used as reliable indicators for the photo-luminescence properties of molecules. We find that the ETGA performs rather similar to the vertical harmonical model. Including anharmonicities in the calculation of internal conversion rates has a moderate effect on the quantitative results in the studied systems. The emission rates are fairly stable with respect to computational parameters, but the internal conversion rate reveals itself to be highly dependent on the choice of the spectral line shape function, particularly the width of the Lorentzian function, associated with homogeneous broadening.

Unravelling the Origin of the Vibronic Spectral Signatures in an Excitonically Coupled Indocarbocyanine Cy3 Dimer

The indocarbocyanine Cy3 dye is widely used to probe the dynamics of proteins and DNA. Excitonically coupled Cy3 dimers exhibit very unique spectral signatures that depend on the interchromophoric geometrical orientation induced by the environment, making them powerful tools to infer the dynamics of their surroundings. Understanding the origin of the dimeric spectral signatures is a necessity for an accurate interpretation of the experimental results. In this work, we simulate the vibronic spectrum of an experimentally well-studied Cy3 dimer, and we explain the origin of the experimental signatures present in its linear absorption spectrum. The Franck-Condon harmonic approximations, among other tests, are used to probe the factors contributing to the spectrum. It is found that the first peak in the absorption spectrum originates from the lower energy excitonic state, while the next two peaks are vibrational progressions of the higher energy excitonic state. The polar solvent plays a crucial role in the appearance of the spectrum, being responsible for the localized S1 minimum, which results in an increased intensity of the first peak.

Extensive Analysis of the Parameters Influencing Radiative Rates Obtained through Vibronic Calculations

Defining a theoretical model systematically delivering accurate ab initio predictions of the fluorescence quantum yields of organic dyes is highly desirable for designing improved fluorophores in a systematic rather than trial-and-error way. To this end, the first required step is to obtain reliable radiative rates (kr), as low kr typically precludes effective emission. In the present contribution, using a series of 10 substituted phenyls with known experimental kr, we analyze the impact of the computational protocol on the kr determined through the thermal vibration correlation function (TVCF) approach on the basis of time-dependent density functional theory (TD-DFT) calculations of the energies, structures, and vibrational parameters. Both the electronic structure (selected exchange-correlation functional, application or not of the Tamm-Dancoff approximation) and the vibronic parameters (line-shape formalism, coordinate system, potential energy surface model, and dipole expansion) are tackled. Considering all possible combinations yields more than 3500 cases, allowing to extract statistically-relevant information regarding the impact of each computational parameter on the magnitude of the estimated kr. It turns out that the selected vibronic model can have a significant impact on the computed kr, especially the potential energy surface model. This effect is of the same order of magnitude as the difference noted between B3LYP and CAM-B3LYP estimates. For the treated compounds, all evaluated functionals do deliver reasonable trends, fitting the experimental values.

Identification via Virtual Screening of Emissive Molecules with a Small Exciton-Vibration Coupling for High Color Purity and Potential Large Exciton Delocalization

A sequence of quantum chemical computations of increasing accuracy was used in this work to identify molecules with small exciton reorganization energy (exciton-vibration coupling), of interest for light emitting devices and coherent exciton transport, starting from a set of ∼4500 known molecules. We validated an approximate computational approach based on single-point calculations of the force in the excited state, which was shown to be very efficient in identifying the most promising candidates. We showed that a simple descriptor based on the bond order could be used to find molecules with potentially small exciton reorganization energies without performing excited state calculations. A small set of chemically diverse molecules with a small exciton reorganization energy was analyzed in greater detail to identify common features leading to this property. Many such molecules display an A-B-A structure where the bonding/antibonding patterns in the fragments A are similar in HOMO and LUMO. Another group of molecules with small reorganization energy displays instead HOMO and LUMO with a strong nonbonding character.

Internal conversion rates from the extended thawed Gaussian approximation: Theory and validation

The theoretical prediction of the rates of nonradiative processes in molecules is fundamental in assessing their emissive properties. In this context, global harmonic models have been widely used to simulate vibronic spectra as well as internal conversion rates and to predict photoluminescence quantum yields. However, these simplified models suffer from the limitations that are inherent to the harmonic approximation and can have a severe effect on the calculated internal conversion rates. Therefore, the development of more accurate semiclassical methods is highly desirable. Here, we introduce a procedure for the calculation of nonradiative rates in the framework of the time-dependent semi-classical Extended Thawed Gaussian Approximation (ETGA). We systematically investigate the performance of the ETGA method by comparing it to the adiabatic and vertical harmonic methods, which belong to the class of widely used global harmonic models. Its performance is tested in potentials that cannot be treated adequately by global harmonic models, beginning with Morse potentials of varying anharmonicity followed by a double well potential. The calculated radiative and nonradiative internal conversion rates are compared to reference values based on exact quantum dynamics. We find that the ETGA has the capability to predict internal conversion rates in anharmonic systems with an appreciable energy gap, whereas the global harmonic models prove to be insufficient.

Franck-Condon spectra of unbound and imaginary-frequency vibrations via correlation functions: a branch-cut free, numerically stable derivation

Molecular electronic spectra can be represented in the time domain as auto-correlation functions of the initial vibrational wavepacket. We present a derivation of the harmonic vibrational auto-correlation function that is valid for both real and imaginary harmonic frequencies. The derivation rests on Lie algebra techniques that map otherwise complicated exponential operator arithmetic to simpler matrix formulae. The expressions for the zero- and finite-temperature harmonic auto-correlation functions have been carefully structured both to be free of branch-cut discontinuities and to remain numerically stable with finite-precision arithmetic. Simple extensions correct the harmonic Franck-Condon approximation for the lowest-order anharmonic and Herzberg-Teller effects. Quantitative simulations are shown for several examples, including the electronic absorption spectra of F2, HOCl, CH2NH, and NO2.

Modeling the Electronic Absorption Spectra of the Indocarbocyanine Cy3

Accurate modeling of optical spectra requires careful treatment of the molecular structures and vibronic, environmental, and thermal contributions. The accuracy of the computational methods used to simulate absorption spectra is limited by their ability to account for all the factors that affect the spectral shapes and energetics. The ensemble-based approaches are widely used to model the absorption spectra of molecules in the condensed-phase, and their performance is system dependent. The Franck-Condon approach is suitable for simulating high resolution spectra of rigid systems, and its accuracy is limited mainly by the harmonic approximation. In this work, the absorption spectrum of the widely used cyanine Cy3 is simulated using the ensemble approach via classical and quantum sampling, as well as, the Franck-Condon approach. The factors limiting the ensemble approaches, including the sampling and force field effects, are tested, while the vertical and adiabatic harmonic approximations of the Franck-Condon approach are also systematically examined. Our results show that all the vertical methods, including the ensemble approach, are not suitable to model the absorption spectrum of Cy3, and recommend the adiabatic methods as suitable approaches for the modeling of spectra with strong vibronic contributions. We find that the thermal effects, the low frequency modes, and the simultaneous vibrational excitations have prominent contributions to the Cy3 spectrum. The inclusion of the solvent stabilizes the energetics significantly, while its negligible effect on the spectral shapes aligns well with the experimental observations.

On the Importance of Well-Defined Thermal Correlation Functions in Simulating Vibronic Spectra

Two difficulties associated with the computations of thermal vibrational correlation functions are discussed. The first one is the lack of a well-behaved expression that is valid at both high-temperature and T → 0 K limits. Specifically, if the partition function and the propagator are considered separately, then thermal vibrational correlation functions may have an indeterminate form 0/0 in the limit T → 0 K. This difficulty is resolved when the partition function and the propagator are jointly considered in the harmonic approximation, which allows a problematic term that emanates from the zero-point energy to be canceled out, thereby producing a thermal correlation function with a determinate form in the T → 0 K limit. The second difficulty is related to the multivaluedness of the vibrational correlation function. We show numerically that an improper selection of branch leads to discontinuities in the computed correlation function and incorrect vibronic spectra. We propose a phase tracking procedure that ensures continuity of both real and imaginary parts of the correlation function to recover the correct spectra. We support our findings by simulating the UV-vis absorption spectra of pentacene at 4 K and benzene at 298 K. Both are found to be in good agreement with their experimental counterparts.

Comparison of Linear Response Theory, Projected Initial Maximum Overlap Method, and Molecular Dynamics-Based Vibronic Spectra: The Case of Methylene Blue

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 S₁ and dark S₂ surfaces near the S₁ 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 S₁/S₂ 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.

Vibronic Spectra of π-Conjugated Systems with a Multitude of Coupled States: A Protocol Based on Linear Vibronic Coupling Models and Quantum Dynamics Tested on Hexahelicene

Hexahelicene is a prototype of an extended π-conjugated system with axial chirality. Its absorption (ABS) and electronic circular dichroism (ECD) spectra show vibronic features and strong nonadiabatic effects, challenging currently available computational methods. Here, we compute the nonadiabatic ABS and ECD vibronic spectra of hexahelicene in the full energy range, covering ∼2 eV and 14-18 coupled electronic states, including all of the relevant nuclear coordinates. To this end, we exploit a recently proposed protocol that uses time-dependent density functional theory to parameterize linear vibronic coupling models comprising several electronic states. Spectra are computed through quantum dynamical propagations with multiconfigurational time-dependent Hartree methods. Our results nicely reproduce the experimental spectra providing an assignment of the main observed bands. On the contrary, we document that the application of the Herzberg-Teller intensity-borrowing theory leads to large artifacts. The proposed approach is of general applicability for rigid systems and represents a viable tool for studying the photophysical properties of π-conjugated systems characterized by a dense manifold of interacting electronic states.

Insights for vibronic effects on spectral shapes of electronic circular dichroism and circularly polarized luminescence of aza[7]helicene

We present a systematic study of vibrationally resolved absorption (ABS), electronic circular dichroism (ECD), emission (EMI), and circularly polarized luminescence (CPL) of aza[7]helicene. Because of the rare experience of theoretical CPL calculation, a variety of harmonic models have been employed to compute the vibronic structures. To fully understand the vibronic effects on the spectral shapes, Franck-Condon (FC) and Herzberg-Teller (HT) contributions, Duschinsky mixings and temperature effect have all been taken into consideration. The performance of different alternative approximate methods has been carefully compared and discussed in detail. The results show that Vertical Hessian (VH) model has a slight better performance on the spectral shapes than Adiabatic Hessian (AH), especially for CPL spectrum. The thermal excitation effect has led to a reduced resolution and a broader spectral width. The moderate HT effects on the different spectral shapes have been addressed. The dissymmetry factors have been correctly reproduced and the main vibronic features of the four different spectral shapes have been successfully captured. A good estimation of the overall spectral width, relative position and relative height of different spectral bands has been presented. The nice agreement with the experiment allows us to present a detailed interpretation of the spectral shapes.

Vibrationally resolved electronic circular dichroism and circularly polarized luminescence spectra of a boron-fused double helicene: A theoretical study

In this work, we present the theoretical study of the vibrationally resolved absorption (ABS), electronic circular dichroism (ECD), emission (EMI), and circularly polarized luminescence (CPL) spectra of a boron-fused double helicene, with a detailed and complete discussion of the alternative possible approximate methods. Given the fact that few examples of CPL calculations exist, the potential energy surfaces (PESs) have been constructed and compared with Adiabatic (AH) and Vertical Hessian (VH) models. All the vibronic calculations have accounted for Duschinsky mixings, Franck-Condon (FC) effect and Herzberg-Teller (HT) contribution. Moreover, different HT expansions have been checked and compared, by computing the derivatives of the electric and magnetic dipole transition moments around the equilibrium geometries of the initial and final states. Our results show that both AH and VH models have well reproduced the experimental vibronic structures and VH model shows a better performance in the simulation of spectral lineshapes. They also show that HT effects dominate the shapes of EMI and CPL, tuning the relative heights of the different vibronic peaks, improving the agreement with the experiment for EMI. Moreover, HT effects are the main reason for the mirror-symmetry breaking between ECD and CPL spectra. Furthermore, interesting interference effects between FC and HT contributions have also been addressed.

Extremely solvent-enhanced absorbance and fluorescence of carbazole interpreted using a damped Franck-Condon simulation

Extremely solvent-enhanced absorption and fluorescence spectra of carbazole were investigated by performing a generalized multi-set damped Franck-Condon spectral simulation. Experimental absorption and fluorescence spectra of carbazole in the gas phase were first well reproduced by performing an un-damped Franck-Condon simulation, but a one-set scaling damped Franck-Condon simulation severely underestimated the intensities of the peaks of experimental absorption and fluorescence spectra of carbazole in n-hexane. Then, a multi-set scaling damped Franck-Condon simulation was proposed and carried out for simulating the extremely solvent-enhanced absorbance and fluorescence, and here, the simulated spectra agreed well with the experimental ones. Five (four) representative solvent-enhanced normal modes corresponding to the combination of ring stretching and ring breathing vibrational motions were determined to be responsible for enhanced absorbance (fluorescence) in n-hexane solution. Furthermore, different scalings were applied to the ground and first-excited states, resulting in different enhancement of absorbance and fluorescence, and this analysis revealed atoms in the carbazole interacting with n-hexane solvent molecules and, hence, leading to different normal-mode vibrational vector patterns in the ground and first-excited states, respectively. Basically, the same conclusion was drawn from a simulation with HF-CIS and the three functionals (TD)B3LYP, (TD)B3LYP-35, and (TD)BHandHLYP. The present multi-set scaling damped Franck-Condon simulation scheme was demonstrated to successfully interpret extremely solvent-enhanced absorbance and fluorescence of carbazole in n-hexane-solvent.

Unspecified verticality of Franck–Condon transitions, absorption and emission spectra of cyanine dyes, and a classically inspired approximation

The computed vertical energy, E_v,a/f, from the equilibrium geometry of the initial electronic state is frequently considered as representative of the experimental excitation/emission energy, E_abs/fl = hc/λ_max. Application of the quantum mechanical version of the Franck–Condon principle does not involve precise specification of nuclear positions before, after, or during an electronic transition. Moreover, the duration of an electronic transition is not experimentally accessible in spectra with resolved vibrational structure. It is shown that computed vibronic spectra based on TDDFT methods and application of quantum mechanical FC analysis predict E_abs = hc/λ_max with a 10-fold improvement in accuracy compared to E_v,a for nine cyanine dyes. It is argued that part of the reason for accuracy when this FC analysis is compared to experiment as opposed to E_v,a/f is the unspecified verticality of transitions in the context of the quantum version of the FC principle. Classical FC transitions that preserve nuclear kinetic energy before and after an electronic transition were previously found to occur at a weighted average of final and initial electronic state molecular geometries known as the r-centroid. Inspired by this approach a qualitative method using computed vertical and adiabatic energies and the harmonic approximation is developed and applied yielding a 5-fold improvement in accuracy compared to E_v,a. This improvement results from the dominance of low frequency vibronic transitions in the cyanine dye major band. The model gives insight into the nature of the redshift when qPCR dye EvaGreen is complexed to λDNA and is applicable to the low frequency band of similar non cyanine dyes such as curcumin. It is found that the computed vibronic cyanine dye spectra from time-dependent FC analysis at 0 K and 298 K show decreased intensity at higher temperature suggestive of increased intensity with restricted motion shown when cyanine dyes are used in biomedical imaging. A 2-layer ONIOM model of the DNA minor groove indicates restricted motion of the TC-1 dye excited state in this setting indicative of enhanced fluorescence.

Insight into cyanine dye λ_max from quantum and classical FC principle; high accuracy with classically inspired approach.

On the simulation of vibrationally resolved electronic spectra of medium-size molecules: the case of styryl substituted BODIPYs

BODIPY dyes are used in a variety of applications because of their peculiar spectroscopic and photo-physical properties that vary depending on the stereochemistry of the functional groups attached to the boron-dipyrromethene core structure. In this work, we have applied several computational methods, adapted for semi-rigid molecules based on the Franck-Condon principle, for the study of the optical properties of BODIPY systems and for the understanding of the influence of functional groups on their spectroscopic features. We have analyzed the electronic spectra of two styryl substituted BODIPY molecules of technological interest, properly taking into account the vibronic contribution. For comparison with recently recorded experimental data in methanol, the vibrationally resolved electronic spectra of these systems were computed using both Time-Independent (TI) and Time-Dependent (TD) formalisms. The first step toward the analysis of optical properties of the styryl modified BODIPYs was a benchmark of several density functionals, to select the most appropriate one. We have found that all benchmarked functionals provide good results in terms of band shape but some of them show strong discrepancies in terms of band position. Beyond the issue of the electronic structure calculation method, different levels of sophistication can be adopted for the calculation of vibronic transitions. In this study, the effect of mode couplings and the influence of the Herzberg-Teller terms on the theoretical spectra has been investigated. It has been found that all levels of theory considered give reproducible results for the investigated systems: band positions and shapes are similar at all levels and little improvements have been found in terms of band shape with the inclusion of Herzberg-Teller effect. Inclusion of temperature effects proved to be challenging due to the important impact of large amplitude motions. Better agreement can be achieved by adopting a suitable set of coordinates coupled with a reduced-dimensionality scheme.

Assessment of Electron Propagator Methods for the Simulation of Vibrationally Resolved Valence and Core Photoionization Spectra

The analysis of photoelectron spectra is usually facilitated by quantum mechanical simulations. Because of the recent improvement of experimental techniques, the resolution of experimental spectra is rapidly increasing, and the inclusion of vibrational effects is usually mandatory to obtain a reliable reproduction of the spectra. With the aim of defining a robust computational protocol, a general time-independent formulation to compute different kinds of vibrationally resolved electronic spectra has been generalized to also support photoelectron spectroscopy. The electronic structure data underlying the simulation are computed using different electron propagator approaches. In addition to the more standard approaches, a new and robust implementation of the second-order self-energy approximation of the electron propagator based on a transition operator reference (TOEP2) is presented. To validate our implementation, a series of molecules has been used as test cases. The result of the simulations shows that, for ultraviolet photoionization spectra, the more accurate nondiagonal approaches are needed to obtain a reliable reproduction of vertical ionization energies but that diagonal approaches are sufficient for energy gradients and pole strengths. For X-ray photoelectron spectroscopy, the TOEP2 approach, besides being more efficient, is also the most accurate in the reproduction of both vertical ionization energies and vibrationally resolved bandshapes.

High-level Ab Initio Absorption Spectra Simulations of Neutral, Anionic and Neutral+ Chromophore of Green Fluorescence Protein Chromophore Models in Gas Phase and Solution

Semiclassical ab initio simulations of the absorption spectra of neutral and anionic p-hydroxybenzylidene-2,3-dimethylimidazolinone (p-HBDI), a model chromophore of green fluorescent protein (GFP) and of a positively charged neutral (N+)-HBDI chromophore model, were performed in gas phase with the resolution-of-identity algebraic diagrammatic construction through second-order (RI-ADC(2)) method. The calculated absorption spectra in gas phase are composed of one band centered at 3.51 eV (HBDI), 2.50 eV (HBDI^- ) and 3.02 eV ((N+)-HBDI) owing to the absorption of the first ¹ ππ* transition. Band maxima are redshifted by ~0.1 eV with respect to the corresponding vertical energies. The COSMO-RI-ADC(2) calculations of the first vertical excitation energy of HBDI, HBDI^- and (N+)-HBDI forms in polar solution including microsolvation simulate the observed solvent redshift for neutral HBDI and the solvent blueshift of the HBDI^- and (N+)-HBDI forms. The state-specific solvation approach applied to TDDFT calculations reproduced the experimental solvent shifts for the three HBDI forms, demonstrating a more accurate theoretical description as compared to the linear-response TDDFT approach.

A proposed new scheme for vibronically resolved time-dependent photoelectron spectroscopy: pump-repump-continuous wave-photoelectron spectroscopy (prp-cw-pes)

We propose a new scheme for time-resolved photoelectron spectroscopy denoted as pump-repump-continuous wave-photoelectron spectroscopy (prp-cw-pes). This scheme is comprised of two femtosecond laser (pump) pulses under cw illumination (probe). By changing the time-delay between pump and repump lasers one can manipulate the populations of vibronic levels in electronic excited states. The cw laser acts on for a long time and establishes resonance between excited states and the continuum photo-ionized states. Sharp spectra can be obtained from the resonance condition [formula: see text]. The intensities in the spectra are sensitive to the time-delay between the pump-repump pulses, but only depend on the populations of excited states, not the phase relations (coherence). As a result, each time-delayed snapshot spectrum is a weighted sum of so-called fingerprints, where a fingerprint is the vibrationally resolved photoelectron spectrum for a single vibronic excited state. The latter information can potentially be simulated reliably using vibronic models and wave packet propagation methods. In the easiest application of the experiment, different time-delays produce different spectra, for a single molecular system. This wealth of experimental data can be fitted to an, ideally small, set of theoretical fingerprints by adjusting the populations as fitting parameters. This technique might be able to distinguish between closely related molecular species. Adopting a different viewpoint, the proposed scheme can also be employed to monitor the time-dependent dynamics by changing the phase relationship between the pump and repump laser that can be viewed as a "control mechanism" employed in wave packet interferometry. Simplifications arise as the change in the spectra is due to the changing populations, not because of the coherence. In this paper, we outline the ideas behind the scheme and illustrate the ideas theoretically using simple model systems.

A Gaussian Wave Packet Propagation Approach to Vibrationally Resolved Optical Spectra at Non-Zero Temperatures

An effective time dependent approach based on a method that is similar to the Gaussian wave packet propagation (GWP) technique of Heller is developed for the computation of vibrationally resolved electronic spectra at finite temperatures in the harmonic, Franck-Condon/Hertzberg-Teller approximations. Since the vibrational thermal density matrix of the ground electronic surface and the time evolution operator on that surface commute, it is possible to write the spectrum generating correlation function as a trace of the time evolved doorway state. In the stated approximations, the doorway state is a superposition of the harmonic oscillator zero and one quantum eigenfunctions and thus can be propagated by the GWP. The algorithm has an O(N(3)) dependence on the number of vibrational modes. An application to pyrene absorption spectrum at two temperatures is presented as a proof of the concept.

Theoretical investigation of the broad one-photon absorption line-shape of a flexible symmetric carbazole derivative

The one-photon absorption spectrum of a carbazole derivative has been studied by employing density functional response theory combined with a mixed quantum/classical approach to simulate the spectral shape.

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.

Vibronic-structure tracking: a shortcut for vibrationally resolved UV/Vis-spectra calculations

The vibrational coarse structure and the band shapes of electronic absorption spectra are often dominated by just a few molecular vibrations. By contrast, the simulation of the vibronic structure even in the simplest theoretical models usually requires the calculation of the entire set of normal modes of vibration. Here, we exploit the idea of the mode-tracking protocol [M. Reiher and J. Neugebauer, J. Chem. Phys. 118, 1634 (2003)] in order to directly target and selectively calculate those normal modes which have the largest effect on the vibronic band shape for a certain electronic excitation. This is achieved by defining a criterion for the importance of a normal mode to the vibrational progressions in the absorption band within the so-called "independent mode, displaced harmonic oscillator" (IMDHO) model. We use this approach for a vibronic-structure investigation for several small test molecules as well as for a comparison of the vibronic absorption spectra of a truncated chlorophyll a model and the full chlorophyll a molecule. We show that the method allows to go beyond the often-used strategy to simulate absorption spectra based on broadened vertical excitation peaks with just a minimum of computational effort, which in case of chlorophyll a corresponds to about 10% of the cost for a full simulation within the IMDHO approach.

Environmental and complexation effects on the structures and spectroscopic signatures of organic pigments relevant to cultural heritage: the case of alizarin and alizarin-Mg(II)/Al(III) complexes

An integrated computational approach allowed an unbiased analysis of optical and structural properties of alizarin-based pigments, which can be directly compared with experimental results. Madder lake pigments have been modeled by Mg(II)- and Al(III)-coordinated alizarin taking into account solvation and metal-linkage effects, responsible for colour modifications. Moreover, different environmental conditions have been analyzed for free alizarin, showing in all cases semi-quantitative agreement with experimental spectroscopic data (UV-VIS). Our results point out the ability of in silico approaches to unravel the subtle interplay of stereo-electronic, dynamic, and environmental effects in tuning the physico-chemical properties of pigments relevant to cultural heritage.

Franck–Condon factors—Computational approaches and recent developments

Algorithms and methodologies for the calculation of Franck–Condon integrals are reviewed. Starting from the standard approach based on the Cartesian representation of the normal modes and the use of Duschinsky's transformation and recursion formulas, methods for treating large displacements of the equilibrium positions and anharmonic and non-Condon effects are considered. Finally, focusing attention on problems arising from the use of recurrence relations, some of the proposed solutions are critically reviewed along with new alternative approaches.

Comparison of vertical and adiabatic harmonic approaches for the calculation of the vibrational structure of electronic spectra

The calculation of the vibrational structure associated to electronic spectra in large molecules requires a Taylor expansion of the initial and final state potential energy surface (PES) around some reference nuclear structure. Vertical (V) and adiabatic (A) approaches expand the final state PES around the initial-state (V) or final-state (A) equilibrium structure. Simplest models only take into account displacements of initial- and final-state minima, intermediate ones also allow for difference in frequencies and more accurate models introduce the Dushinsky effect through the computation of the Hessians of both the initial and final state. In this contribution we summarize and compare the mathematical expressions of the complete hierarchy of V and A harmonic models and we implement them in a numerical code, presenting a detailed comparison of their performance on a number of prototypical systems. We also address non-Condon effects through linear expansions of the transition dipole as a function of nuclear coordinates (Herzberg-Teller effect) and compare the results of expansions around initial and final state equilibrium geometries. By a throughout analysis of our results we highlight a number of general trends in the relative performance of the models that can provide hints for their proper choice. Moreover we show that A and V models including final state PES Hessian outperform the simpler ones and that discrepancies in their predictions are diagnostic for failure of harmonic approximation and/or of Born-Oppenheimer approximation (existence of remarkable geometry-dependent mixing of electronic states).