First-principles calculations of anharmonic and deuteration effects on the photophysical properties of polyacenes and porphyrinoids

Molecular phosphorescence enhancement by the plasmon field of metal nanoparticles

A theoretical model is proposed that allows the estimation of the quantum yield of phosphorescence of dye molecules in the vicinity of plasmonic nanoparticles. For this purpose, the rate constants of the radiative and nonradiative intramolecular transitions for rhodamine 123 (Rh123) and brominated rhodamine (Rh123-2Br) dyes have been calculated. The plasmon effect of Ag nanoparticles on various types of luminescence processes has been studied both theoretically and experimentally. We show that in the presence of a plasmonic nanoparticle, the efficiency of the immediate and delayed fluorescence increases significantly. The phosphorescence rate of the rhodamine dyes also increases near plasmonic nanoparticles. The long-lived luminescence i.e., delayed fluorescence and phosphorescence is more enhanced for Rh123-2Br than for Rh123. The largest phosphorescence quantum yield is obtained when the dye molecule is at a distance of 4-6 nm from the nanoparticle surface. Our results can be used in the design of plasmon-enhancing nanostructures for light-emitting media, organic light-emitting diodes, photovoltaic devices, and catalysts for activation of molecular oxygen.

Vibronic Coupling Drives the Ultrafast Internal Conversion in a Functionalized Free-Base Porphyrin

Internal conversion (IC) is a common radiationless transition in polyatomic molecules. Theory predicts that molecular vibrations assist IC between excited states, and ultrafast experiments can provide insight into their structure-function relationship. Here we elucidate the dynamics of the vibrational modes driving the IC process within the Q band of a functionalized porphyrin molecule. Through a combination of ultrafast multidimensional spectroscopies and theoretical modeling, we observe a 60 fs Qy-Qx IC and demonstrate that it is driven by the interplay among multiple high-frequency modes. Notably, we identify 1510 cm-1 as the leading tuning mode that brings the porphyrin to an optimal geometry for energy surface crossing. By employing coherent wave packet analysis, we highlight a set of short-lived vibrations (1200-1400 cm-1), promoting the IC within ≈60 fs. Furthermore, we identify one coupling mode (1350 cm-1) that is responsible for vibronic mixing within the Q states. Our findings indicate that porphyrin-core functionalization modulates IC effectively, offering new opportunities in photocatalysis and optoelectronics.

Intramolecular rate-constant calculations based on the correlation function using temperature dependent quantum Green's functions

A theoretical method for calculating rate constants for internal conversion (IC), intersystem crossing (ISC) and radiative (R) electronic transitions is presented. The employed method uses temperature-dependent quantum Green's functions, which give the opportunity to consider almost any nth-order polynomial perturbation operator and the influence of external electromagnetic fields on the rate constants. The rate constants of the IC, ISC and R processes are calculated for two important indocyanine molecules namely indocyanine green (ICG) and heptamethine cyanine (IR808) at the Franck-Condon level using the temperature-dependent quantum Green's function approach. Calculations at the time-dependent density functional theory level with the MN15 functional show that ICG and IR808 have only one triplet state below the S1 state. The main deactivation channel of the S1 state is the IC process with a large (kIC(S1 → S0)) rate constant of ∼109-1011 s-1. The estimated quantum yield of fluorescence (φfl) is ∼0.001-0.24 for the two studied molecules, which agrees rather well with experimental values. Thus, the present approach enables calculations of the three kinds of rate constants and the quantum yield of fluorescence using the same computational methodology.

Internal conversion induced by external electric and magnetic fields

We have developed a new methodology for calculating contributions to the rate constants (kIC) of internal conversion that are induced by external electric (kIC-E) or magnetic (kIC-M) fields. The influence of the external electric and magnetic fields on the kIC was estimated for seven representative molecules. We show that the kIC-E contribution calculated at a field strength of 1011 V m-1 is generally as large as the kIC rate constant in the absence of the external field. For indocyanine green, azaoxa[8]circulene, and pyromitene 567, the kIC-E contribution is as large as kIC already at a field strength of 109 V m-1. Such electric-field strengths occur for example in plasmonic studies and in strong laser-field experiments. The induced effect on the kIC rate constant should be accounted for in calculations of photophysical properties of molecules involved in such experiments. The induced effect of an external magnetic field on kIC can be neglected in experiments on Earth because the magnetic contribution becomes significant only at very strong magnetic fields of 104-105 T that cannot be achieved on Earth. However, the magnetic effect on the rate constant of internal conversion can be important in astrophysical studies, where extremely strong magnetic fields occur near neutron stars and white dwarfs.

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.

Theoretical studies on thermally activated delayed fluorescence of "carbene-metal-amide" Cu and Au complexes: geometric structures, excitation characters, and mechanisms

"Carbene-metal(I)-amide" (CMA) complexes have garnered significant attention due to their remarkable properties and potential TADF applications in organic electronics. However, the atomistic working mechanism is still elusive. Herein, we chose two CMA complexes, i.e., cyclic (alkyl)(amino) carbene-copper[gold](I)-carbazole (CAAC-Cu[Au]-Cz), and employed both DFT and TD-DFT methods, in combination with radiative and nonradiative rate calculations, to investigate geometric and electronic structures of these two complexes in the ground and excited states, including orbital compositions, electronic transitions, absorption and emission spectra, and the luminescence mechanism. It is found that the coplanar or perpendicular conformations are coexistent in the ground state (S0), the lowest excited singlet state (S1), and the triplet state (T1). Both the coplanar and perpendicular S1 and T1 states have similar ligand-to-ligand charge transfer (LLCT) character between CAAC and Cz, and some charge-transfer character between metal atoms and ligands, which is beneficial to minimize the singlet-triplet energy gaps (ΔEST) and increase the spin-orbit coupling (SOC). An interesting three-state (S0, S1, T1) model involving two regions (coplanar and perpendicular) is proposed to rationalize the experimental TADF phenomena in the CMA complexes. In addition to the coplanar ones, the perpendicular S1 and T1 states also play a role in promoting the repopulation of the coplanar S1 exciton, which is a primary source for the delayed fluorescence.

Influence of plasmons on the luminescence properties of solvatochromic merocyanine dyes with different solvatochromism

The effect of localized surface plasmon resonance (LSPR) of a system consisting of a highly dipolar merocyanine dye and a silver nanoparticle (NP) was studied experimentally and theoretically. A theoretical model for estimating the fluorescence quantum yield (φfl) using quantum chemical calculations of intramolecular and intermolecular electronic transition rate constants was developed. Calculations show that the main deactivation channels of the lowest excited singlet state of the studied merocyanines are internal conversion (kIC(S1 → S0)) and fluorescence (kr(S1 → S0)). The intersystem-crossing transition has a low probability due to the large energy difference between the singlet and triplet levels. In the presence of plasmonic NPs, the fluorescence quantum yield is increased by a factor of two according to both experiment and computations. The calculated values of φfl, when considering changes in kr(S1 → S0) and the energy-transfer rate constant (ktransfer) from the dye to the NP was also twice as large at distances of 6-8 nm between the NP and the dye molecule. We also found that the LSPR effect can be increased or decreased depending on the value of the dielectric constant (εm) of the environment.

Internal conversion rate constant calculations considering Duschinsky, anharmonic and Herzberg-Teller effects

A novel method for calculating rate constants for internal conversion (k_IC) that simultaneously accounts for Duschinsky, anharmonic and Herzberg-Teller effects has been developed and implemented. This method has been applied to robust planar molecules like tetraoxa[8]circulene (4B), free-base porphyrin (H2P) and pyrometene (PM567) with small Duschinsky rotation (i.e. with almost identical normal coordinates in the ground and excited states) and to poly[n]fluorenes (P[n]F) (n = 2-14) with a substantial Duschinsky rotation. The obtained results show that the Duschinsky effect is large in the harmonic approximation, whereas it is in general much smaller in the anharmonic approximation. The Duschinsky effect is found to be large for high frequency vibrational modes with energies of ∼3300 cm^-1 such as the X-H (X = C, N and O) stretching modes that mix in the S₁ → S₀ electronic transition. However, even in this case, the increase in k_IC due to the Duschinsky effect does not exceed one order of magnitude. The calculations show that anharmonic contributions to k_IC are larger than Herzberg-Teller contributions which in turn are larger than contributions from the Duschinsky effect ANH > HT > Du. We also show that an approximation, where only X-H bonds are considered in the k_IC calculation, is accurate even for P[n]F (n = 2-14).

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.

Interexcited State Photophysics I: Benchmarking Density Functionals for Computing Nonadiabatic Couplings and Internal Conversion Rate Constants

We present the first benchmarking study of nonadiabatic matrix coupling elements (NACMEs) calculated using different density functionals. Using the S₁ → S₀ transition in perylene solvated in toluene as a case study, we calculate the photophysical properties and corresponding rate constants for a variety of density functionals from each rung of Jacob's ladder. The singlet photoluminescence quantum yield (sPLQY) is taken as a measure of accuracy, measured experimentally here as 0.955. Important quantum chemical parameters such as geometries, absorption, emission, and adiabatic energies, NACMEs, Hessians, and transition dipole moments were calculated for each density functional basis set combination (data set) using density functional theory based multireference configuration interaction (DFT/MRCI) and compared to experiment where possible. We were able to derive simple relations between the TDDFT and DFT/MRCI photophysical properties; with semiempirical damping factors of ∼0.843 ± 0.017 and ∼0.954 ± 0.064 for TDDFT transition dipole moments and energies to DFT/MRCI level approximations, respectively. NACMEs were dominated by out-of-plane derivative components belonging to the center-most ring atoms with weaker contributions from perturbations along the transverse and longitudinal axes. Calculated theoretical spectra compared well to both experiment and literature, with fluorescence lifetimes between 7.1 and 12.5 ns, agreeing within a factor of 2 with experiment. Internal conversion (IC) rates were then calculated and were found to vary wildly between 10⁶-10¹⁶ s^-1 compared with an experimental rate of the order 10⁷ s^-1. Following further testing by mixing data sets, we found a strong dependence on the method used to obtain the Hessian. The 5 characterized data sets ranked in order of most promising are PBE0/def2-TZVP, ωB97XD/def2-TZVP, HCTH407/TZVP, PBE/TZVP, and PBE/def2-TZVP.

The dominant nature of Herzberg-Teller terms in the photophysical description of naphthalene compared to anthracene and tetracene

The first order and second order corrected photoluminescence quantum yields are computed and compared to experiment for naphthalene in this manuscript discussing negative results. Results for anthracene and tetracene are recalled from previous work (Manian et al. in J Chem Phys 155:054108, 2021), and the results for all three polyacenes are juxtaposed to each other. While at the Franck-Condon point, each of the three noted polyacenes were found to possess a quantum yield near unity. Following the consideration of Herzberg-Teller effects, quantum yields stabilised for anthracene and tetracene to 0.19 and 0.08, respectively. Conversely, the second order corrected quantum yield for naphthalene was found to be 0.91. Analysis of this result showed that while the predicted non-radiative pathways correlate well with what should be expected, the approximation used to calculate second order corrected fluorescence, which yielded very positive results for many other molecular systems, here is unable to account for strong second order contributions, resulting in a grossly overestimated rate of fluorescence. However, substitution of an experimental radiative rate results in a quantum yield of 0.33. This work extols the importance of Herzberg-Teller terms in photophysical descriptions of chromophores, and highlights those cases in which a treatment beyond the above approximation is required.

Spectral Signatures of Hydrogen Thioperoxide (HOSH) and Hydrogen Persulfide (HSSH): Possible Molecular Sulfur Sinks in the Dense ISM

For decades, sulfur has remained underdetected in molecular form within the dense interstellar medium (ISM), and somewhere a molecular sulfur sink exists where it may be hiding. With the discovery of hydrogen peroxide (HOOH) in the ISM in 2011, a natural starting point may be found in sulfur-bearing analogs that are chemically similar to HOOH: hydrogen thioperoxide (HOSH) and hydrogen persulfide (HSSH). The present theoretical study couples the accuracy in the anharmonic fundamental vibrational frequencies from the explicitly correlated coupled cluster theory with the accurate rotational constants provided by canonical high-level coupled cluster theory to produce rovibrational spectra for use in the potential observation of HOSH and HSSH. The ν6 mode for HSSH at 886.1 cm−1 is within 0.2 cm−1 of the gas-phase experiment, and the B0 rotational constant for HSSH of 6979.5 MHz is within 9.0 MHz of the experimental benchmarks, implying that the unknown spectral features (such as the first overtones and combination bands) provided herein are similarly accurate. Notably, a previous experimentally-attributed 2ν1 mode, at 7041.8 cm−1, has been reassigned to the ν1+ν5 combination band based on the present work’s ν1+ν5 value at 7034.3 cm−1. The most intense vibrational transitions for each molecule are the torsions, with HOSH having a more intense transition of 72 km/mol compared to HSSH’s intensity of 14 km/mol. Furthermore, HOSH has a larger net dipole moment of 1.60 D compared to HSSH’s 1.15 D. While HOSH may be the more likely candidate of the two for possible astronomical observation via vibrational spectroscopy due to the notable difference in their intensities, both HSSH and HOSH have large enough net dipole moments to be detectable by rotational spectroscopy to discover the role these molecules may have as possible molecular sulfur sinks in the dense ISM.

Exciton Dynamics of a Diketo-Pyrrolopyrrole Core for All Low-Lying Electronic Excited States Using Density Functional Theory-Based Methods

Ab initio treatments of interexcited state internal conversion (IC) are more often than not missing from exciton dynamic descriptions, because of their inherent complexity. Here, we define "interexcited state IC" as a same-spin nonradiative transition between states i and j, where i ≠ j ≠ 0. Competing directly with multiexciton processes such as singlet fission or triplet photoupconversion, inclusion of this mechanism in the narrative of molecular photophysics would allow for strategic synthesis of chromophores for more efficient photon-harvesting applications. Herein, we present a robust formalism which can model these rates using density functional theory (DFT)-based methods within the Franck-Condon and Herzberg-Teller regime. Using an unsubstituted diketo-pyrrolopyrrole (DPP) core as a case study, we illustrate the exciton dynamics along the first four excited states for both singlet and triplet manifolds, showing ultrafast same-spin transfer mechanisms due to all excited states, excluding the first triplet level, being in close energetic proximity (within 0.8 eV of each other). The resulting electron same-spin rates outcompete the electron spin-flipping intersystem crossing (ISC) rates, with excitons firmly obeying Kasha's rule as they cascade down from the high-lying excited states toward the lower states. Furthermore, we calculated that only the first singlet excited state displayed a reasonable probability of triplet exciton generation, of ∼40%, with a near-zero chance of the exciton reverting to the singlet manifold once the electron-hole pair are of parallel spin.

The combination of skeleton-engineering and periphery-engineering: a design strategy for organic doublet emitters

A series of radicals based on tris(2,4,6-trichlorophenyl)methyl (TTM) were theoretically designed and evaluated by combining skeleton-engineering and periphery-engineering strategies.

Anharmonic Vibrational Frequencies of Water Borane and Associated Molecules

Water borane (BH3OH2) and borinic acid (BH2OH) have been proposed as intermediates along the pathway of hydrogen generation from simple reactants: water and borane. However, the vibrational spectra for neither water borane nor borinic acid has been investigaged experimentally due to the difficulty of isolating them in the gas phase, making accurate quantum chemical predictions for such properties the most viable means of their determination. This work presents theoretical predictions of the full rotational and fundamental vibrational spectra of these two potentially application-rich molecules using quartic force fields at the CCSD(T)-F12b/cc-pCVTZ-F12 level with additional corrections included for the effects of scalar relativity. This computational scheme is further benchmarked against the available gas-phase experimental data for the related borane and HBO molecules. The differences are found to be within 3 cm−1 for the fundamental vibrational frequencies and as close as 15 MHz in the B0 and C0 principal rotational constants. Both BH2OH and BH3OH2 have multiple vibrational modes with intensities greater than 100 km mol−1, namely ν2 and ν4 in BH2OH, and ν1, ν3, ν4, ν9, and ν13 in BH3OH2. Finally, BH3OH2 has a large dipole moment of 4.24 D, which should enable it to be observable by rotational spectroscopy, as well.

Unified Framework for Photophysical Rate Calculations in TADF Molecules

One of the challenges in organic light-emitting diodes research is finding ways to increase device efficiency by making use of the triplet excitons that are inevitably generated in the process of electroluminescence. One way to do so is by thermally activated delayed fluorescence (TADF), a process in which triplet excitons undergo upconversion to singlet states, allowing them to relax radiatively. The discovery of this phenomenon has ensued a quest for new materials that are able to effectively take advantage of this mechanism. From a theoretical standpoint, this requires the capacity to estimate the rates of the various processes involved in the photophysics of candidate molecules, such as intersystem crossing, reverse intersystem crossing, fluorescence, and phosphorescence. Here, we present a method that is able to, within a single framework, compute all of these rates and predict the photophysics of new molecules. We apply the method to two TADF molecules and show that results compare favorably with other theoretical approaches and experimental results. Finally, we use a kinetic model to show how the calculated rates act in concert to produce different photophysical behavior.

Modeling radiative and non-radiative pathways at both the Franck-Condon and Herzberg-Teller approximation level

Here, we present a concise model that can predict the photoluminescent properties of a given compound from first principles, both within and beyond the Franck-Condon approximation. The formalism required to compute fluorescence, Internal Conversion (IC), and Inter-System Crossing (ISC) is discussed. The IC mechanism, in particular, is a difficult pathway to compute due to difficulties associated with the computation of required bosonic configurations and non-adiabatic coupling elements. Here, we offer a discussion and breakdown on how to model these pathways at the Density Functional Theory (DFT) level with respect to its computational implementation, strengths, and current limitations. The model is then used to compute the photoluminescent quantum yield (PLQY) of a number of small but important compounds: anthracene, tetracene, pentacene, diketo-pyrrolo-pyrrole (DPP), and Perylene Diimide (PDI) within a polarizable continuum model. Rate constants for fluorescence, IC, and ISC compare well for the most part with respect to experiment, despite triplet energies being overestimated to a degree. The resulting PLQYs are promising with respect to the level of theory being DFT. While we obtained a positive result for PDI within the Franck-Condon limit, the other systems require a second order correction. Recomputing quantum yields with Herzberg-Teller terms yields PLQYs of 0.19, 0.08, 0.04, 0.70, and 0.99 for anthracene, tetracene, pentacene, DPP, and PDI, respectively. Based on these results, we are confident that the presented methodology is sound with respect to the level of quantum chemistry and presents an important stepping stone in the search for a tool to predict the properties of larger coupled systems.

Fast estimation of the internal conversion rate constant in photophysical applications

An efficient method for estimating non-adiabatic coupling matrix elements (NACME) and rate constants for internal conversion (k_IC) is presented. The method, based on Plotnikov's theory, requires only calculations of the electronic wave functions and the corresponding electronic excitation energies. Computationally expensive calculations of the derivatives of the electronic wave function with respect to the nuclear coordinates are avoided. When the main accepting modes of the electronic excitation energy are X-H vibrations, the present method can be used for estimating the efficiency of the energy transfer between donor and acceptor molecules. It can also be used in studies of the influence of hydrogen bonding or solvent effect on fluorescence quenching, in studies of vibronic effects of TADF (thermally activated delayed fluorescence) emitters, and for calculating k_IC. Here, k_IC and NACME are calculated for free-base porhyrin, magnesium porphyrin, azulene, naphthalene, pyrene and fluorenone interacting with a solvent molecule. Reverse k_IC and NACME are further calculated for the T₁→ T₂ transition of dibenzothiophene-S,S-dioxide (PTZ-DBTO2), which is used in TADF applications. Finally, we estimate the efficiency of the energy transfer between two large porphyrinoid dimers.

Efficient enumeration of bosonic configurations with applications to the calculation of non-radiative rates

This work presents algorithms for the efficient enumeration of configuration spaces following Boltzmann-like statistics, with example applications to the calculation of non-radiative rates, and an open-source implementation. Configuration spaces are found in several areas of physics, particularly wherever there are energy levels that possess variable occupations. In bosonic systems, where there are no upper limits on the occupation of each level, enumeration of all possible configurations is an exceptionally hard problem. We look at the case where the levels need to be filled to satisfy an energy criterion, for example, a target excitation energy, which is a type of knapsack problem as found in combinatorics. We present analyses of the density of configuration spaces in arbitrary dimensions and how particular forms of kernel can be used to envelope the important regions. In this way, we arrive at three new algorithms for enumeration of such spaces that are several orders of magnitude more efficient than the naive brute force approach. Finally, we show how these can be applied to the particular case of internal conversion rates in a selection of molecules and discuss how a stochastic approach can, in principle, reduce the computational complexity to polynomial time.

Computational study of IR, Raman, and NMR spectra of 4-methylmethcathinone drug

Molecular electronic structure, IR, UV, and NMR spectra of the most popular cathinone, known as mephedrone or 4-methylmethcathinone (4-MMC), is studied thoroughly by quantum chemical calculation in terms of the density functional theory (DFT). Geometry optimization of 4-MMC and its hydrochloride complex is performed with the B3LYP functional, and all vibrational frequencies are analyzed in all details. On this background, the IR and Raman spectra are interpreted. The importance of low-frequency terahertz and Raman spectra is stressed for distinguishing of various MMC isomers. The UV spectrum is calculated by time-dependent DFT method which allows complete interpretation of intense absorption bands at 270 and 210 nm as combinations of various ππ*, nπ*, and charge transfer excitations in amino-phenyl moieties. Very informative analysis of UV absorption and NMR spectra provides useful details on the structure-activity relationship for mephedrone molecule.

A computational study of the vibronic effects on the electronic spectra and the photophysics of aza[7]helicene

We report a computational study on vibronic effects in the spectroscopy, photoinduced processes and decay back to the ground state of aza[7]helicene, a helicene with an unusually high fluorescence quantum yield (QY = 0.39). In a first step, we compute and assign the absorption and electronic circular dichroism (ECD) spectra in its full frequency range from 2.7 to 5.0 eV, accounting for nonadiabatic effects. Then we compute the quantum dynamics of the cascade of ultrafast internal conversions of the highly-excited singlet states to the lowest-energy one S1. Finally we adopt Fermi golden rule rates to compute the QY of the dye, taking into account the competition between the radiative decay and the nonradiative decays to the ground state and to the energy-accessible triplet states. We use time-dependent density functional theory (TD-DFT), including solvent (dichloromethane) effects within the polarizable continuum model, to parameterize a linear vibronic coupling (LVC) model involving the first lowest 12 singlet states and all the normal coordinates. Nonadiabatic spectra and internal conversions dynamics are then computed through wavepacket propagations with the Multilayer (ML) extension of the Multiconfigurational Time Dependent Hartree method (ML-MCTDH). We highlight the molecular vibrations playing a major role in determining the shape of the spectra and analyse the effect of inter-state couplings. At the same time we report a breakdown of perturbative Herzberg-Teller approach. The computed QY is in perfect agreement with experiment and allows us to ascertain that intersystem crossings are the processes limiting the fluorescence from S1. They involve the three lowest triplet states and are made effective by spin-orbit coupling and vibronic effects.