Computing UV/vis spectra from the adiabatic and vertical Franck-Condon schemes with the use of Cartesian and internal coordinates

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.

First-Principles Calculations of Excited-State Decay Rate Constants in Organic Fluorophores

In this Perspective, we discuss recent advances made to evaluate from first-principles the excited-state decay rate constants of organic fluorophores, focusing on the so-called static strategy. In this strategy, one essentially takes advantage of Fermi's golden rule (FGR) to evaluate rate constants at key points of the potential energy surfaces, a procedure that can be refined in a variety of ways. In this way, the radiative rate constant can be straightforwardly obtained by integrating the fluorescence line shape, itself determined from vibronic calculations. Likewise, FGR allows for a consistent calculation of the internal conversion (related to the non-adiabatic couplings) in the weak-coupling regime and intersystem crossing rates, therefore giving access to estimates of the emission yields when no complex photophysical phenomenon is at play. Beyond outlining the underlying theories, we summarize here the results of benchmarks performed for various types of rates, highlighting that both the quality of the vibronic calculations and the accuracy of the relative energies are crucial to reaching semiquantitative estimates. Finally, we illustrate the successes and challenges in determining the fluorescence quantum yields using a series of organic fluorophores.

Pitfall in simulations of vibronic TD-DFT spectra: diagnosis and assessment

In this Communication, we study the effect of spurious oscillations in the profiles of energy derivatives with respect to nuclear coordinates calculated with density functional approximations (DFAs) for formaldehyde, pyridine, and furan in their ground and electronic excited states. These spurious oscillations, which can only be removed using extensive integration grids that increase enormously the CPU cost of DFA calculations, are significant in the case of third- and fourth-order energy derivatives of the ground and excited states computed by M06-2X and ωB97X functionals. The errors in question propagate to anharmonic vibronic spectra computed under the Franck-Condon approximation, i.e., positions and intensities of vibronic transitions are affected to a large extent (shifts as significant as hundreds of cm-1 were observed). On the other hand, the LC-BLYP and CAM-B3LYP functionals show a much less pronounced effect due to spurious oscillations. Based on the results presented herein, we recommend either LC-BLYP or CAM-B3LYP with integration grids (250, 974) (or larger) for numerically stable simulations of vibronic spectra including anharmonic effects.

Chemical Transformations Can Occur during DMS Separations: Lessons Learned from Beer's Bittering Compounds

While developing a DMS-based separation method for beer's bittering compounds, we observed that the argentinated forms of humulone tautomers (i.e., [Hum + Ag]⁺) were partially resolvable in a N₂ environment seeded with 1.5 mol % of isopropyl alcohol (IPA). Attempting to improve the separation by introducing resolving gas unexpectedly caused the peaks for the cis-keto and trans-keto tautomers of [Hum + Ag]⁺ to coalesce. To understand why resolution loss occurred, we first confirmed that each of the tautomeric forms (i.e., dienol, cis-keto, and trans-keto) responsible for the three peaks in the [Hum + Ag]⁺ ionogram were assigned to the correct species by employing collision-induced dissociation, UV photodissociation spectroscopy, and hydrogen-deuterium exchange (HDX). The observation of HDX indicated that proton transfer was stimulated by dynamic clustering processes between IPA and [Hum + Ag]⁺ during DMS transit. Because IPA accretion preferentially occurs at Ag⁺, which can form pseudocovalent bonds with a suitable electron donor, solvent clustering also facilitated the formation of exceptionally stable microsolvated ions. The exceptional stability of these microsolvated configurations disproportionately impacted the compensation voltage (CV) required to elute each tautomer when the temperature within the DMS cell was varied. The disparity in CV response caused the peaks for the cis- and trans-keto species to merge when a temperature gradient was induced by the resolving gas. Moreover, simulations showed that microsolvation with IPA mediates dienol to trans-keto tautomerization during DMS transit, which, to the best of our knowledge, is the first observation of keto/enol tautomerization occurring within an ion-mobility device.

A (TD-)DFT study on photo-NHC catalysis: photoenolization/Diels-Alder reaction of acid fluorides catalyzed by N-heterocyclic carbenes

A comprehensive mechanistic study on the N-heterocyclic carbene (NHC) catalyzed photoenolization/Diels-Alder (PEDA) reaction of acid fluorides was performed in the framework of (time-dependent) density functional theory ((TD)-DFT). The 1,5-hydrogen atom transfer (1,5-HAT) during photoenolization of an ortho-toluoyl azolium salt was found to be feasible via, first, singlet excitation and photoenolization, and then, after crossing to the triplet manifold, populating a biradical dienol which allows for the formation of two ortho-quinodimethane (o-QDM) isomers due to a low rotational barrier. The (Z)-isomer is mostly unproductive through sigmatropic rearrangement back to the starting material while the (E)-isomer reacts in a subsequent concerted Diels-Alder reaction likely as the deprotonated dienolate. The experimentally observed diastereoselectivity is correctly predicted by theory and is determined by a more favorable endo trajectory in the cycloaddition step. These findings demonstrate that ortho-toluoyl azolium species exhibit similar photophysical properties as structurally related benzophenones, highlighting the unique ability of the NHC organocatalyst to transiently alter the excited state properties of an otherwise photoinactive carboxylic acid derivative, thereby expanding the scope of classical carbonyl photochemistry.

Cost-Effective Simulations of Vibrationally-Resolved Absorption Spectra of Fluorophores with Machine-Learning-Based Inhomogeneous Broadening

The results of electronic and vibrational structure simulations are an invaluable support for interpreting experimental absorption/emission spectra, which stimulates the development of reliable and cost-effective computational protocols. In this work, we contribute to these efforts and propose an efficient first-principle protocol for simulating vibrationally-resolved absorption spectra, including nonempirical estimations of the inhomogeneous broadening. To this end, we analyze three key aspects: (i) a metric-based selection of density functional approximation (DFA) so to benefit from the computational efficiency of time-dependent density function theory (TD-DFT) while safeguarding the accuracy of the vibrationally-resolved spectra, (ii) an assessment of two vibrational structure schemes (vertical gradient and adiabatic Hessian) to compute the Franck-Condon factors, and (iii) the use of machine learning to speed up nonempirical estimations of the inhomogeneous broadening. In more detail, we predict the absorption band shapes for a set of 20 medium-sized fluorescent dyes, focusing on the bright ππ^★ S₀ → S₁ transition and using experimental results as references. We demonstrate that, for the studied 20-dye set which includes structures with large structural variability, the preselection of DFAs based on an easily accessible metric ensures accurate band shapes with respect to the reference approach and that range-separated functionals show the best performance when combined with the vertical gradient model. As far as band widths are concerned, we propose a new machine-learning-based approach for determining the inhomogeneous broadening induced by the solvent microenvironment. This approach is shown to be very robust offering inhomogeneous broadenings with errors as small as 2 cm^-1 with respect to genuine electronic-structure calculations, with a total CPU time reduced by 98%.

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/.

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.

Effect of protein matrix on CP29 spectra and energy transfer pathways

We investigate energy transfer pathways between strongly coupled chlorophylls (Chls) in the CP29 (LHCII B4.1) antenna complex of Pisum sativum, including the possibility of higher energy states. We test for the environmental effects caused by the protein, membrane and solvent using a hybrid QM/MM approach. Classical molecular dynamics simulations of the full CP29 complex embedded in a DOPC membrane have been performed, followed by calculations of the time dependent DFT spectra of all Chls at several timesteps. The relative orientations of transition dipole moments (TDMs) were specifically analyzed, including and excluding the point charge field (PCF) of the surrounding environment. The PCF is found to drastically shift the spectra of specific Chls, while the majority of Chls is mostly unaffected. The net effect on the sum spectrum is however found to be negligible: The few strong changes in Chl spectra cancel each other due to being opposite in sign. We further find that the spectra of the Chls coordinating to water show a blue shift upon introduction of the environment. Conversely, the spectra of the Chls coordinating to glutamine show a red shift upon activation of the PCF. As the main influence of the PCF for tuning the couplings, we identify the energetic position of the individual chromophores. The fine-tuning, especially for states energetically above the Qy state, is however controlled by the changes in the TDM orientations. We also find an indication for the PCF to steer potentially harmful high energy excitations away from the PSII core complex.

UVPD Spectroscopy of Differential Mobility-Selected Prototropic Isomers of Rivaroxaban

Two ion populations of protonated Rivaroxaban, [C₁₉H₁₈ClN₃O₅S + H]⁺, are separated under pure N₂ conditions using differential mobility spectrometry prior to characterization in a hybrid triple quadrupole linear ion trap mass spectrometer. These populations are attributed to bare protonated Rivaroxaban and to a proton-bound Rivaroxaban-ammonia complex, which dissociates prior to mass-selecting the parent ion. Ultraviolet photodissociation (UVPD) and collision-induced dissociation (CID) studies indicate that both protonated Rivaroxaban ion populations are comprised of the computed global minimum prototropic isomer. Two ion populations are also observed when the collision environment is modified with 1.5% (v/v) acetonitrile. In this case, the protonated Rivaroxaban ion populations are produced by the dissociation of the ammonium complex and by the dissociation of a proton-bound Rivaroxaban-acetonitrile complex prior to mass selection. Again, both populations exhibit a similar CID behavior; however, UVPD spectra indicate that the two ion populations are associated with different prototropic isomers. The experimentally acquired spectra are compared with computed spectra and are assigned to two prototropic isomers that exhibit proton sharing between distal oxygen centers.

Investigating the biological actions of some Schiff bases using density functional theory study

Schiff base ligands have wide varieties of application in several fields. One of which is the biological actions they possess such as anti-fungal, anti-bacterial, anti-malarial, and anti-viral characteristics. In this study, some synthesized phenylimino-based Schiff bases were investigated using density functional theory (DFT) to unravel their biological descriptors. The gas-phase quantum chemical calculation was done on the Schiff base 3-((E)-(phenylimino)methyl)benzene-1,2-diol and other synthesized analogues to evaluate their reactivity and stability properties including the substituent effect on the basic molecule. The Coulomb-attenuating method (CAM-B3LYP) functional was employed for the theoretical calculations. The Nuclear Magnetic Resonance (NMR), Fourier Transform-Infrared (FT-IR), Ultraviolet/visible spectroscopies calculated agrees with the experimental values. The obtained charge transfer and electronic features provide useful information regarding the active sites for biological application in the compounds.

A computational study of the vibrationally-resolved electronic circular dichroism spectra of single-chain transoid and cisoid oligothiophenes in chiral conformations

We simulate the vibronic profile of the electronic circular dichroism (ECD) spectra of oligothiophenes in cisoid and transoid chiral arrangements. We consider oligomers of different lengths, from two to fifteen units, and investigate extensively how the ECD spectral shapes depend on the inter-ring torsions. In general, the molecular structures we consider are not stationary points of the ground state potential energy surface. Therefore, in order to perform vibronic calculations, we present a new computational protocol able to define reduced-dimensionality models where the effect of the off-equilibrium modes is removed. This is done adopting a description of the vibrational motions in curvilinear internal coordinates, and vertical harmonic models coupled with an iterative application of projectors to define energy Hessians, and therefore effective normal modes, in the space complementary to the one of the off-equilibrium coordinates. Although we consider both Franck-Condon and Herzberg-Teller contributions, the results show that transoid twisted ribbons always give rise to monosignated ECD spectra, while bi-signated and multi-signated spectra are expected for cisoid helices. These findings are explained on the basis of the different transition strengths of the lowest electronic states imparted by the different spatial arrangement, that is almost linear for transoid structures and more globular for cisoid ones. We predicted the chiroptical response of a large number of possible molecular arrangements. These data are employed to critically discuss the experimental ECD of polythiophenes in different experimental conditions, forming either aggregates or host-guest complexes. The method here proposed to perform vibronic calculations in reduced-dimensionality models is of general applicability and its potential interest goes beyond the practical application presented here.

Absorption and luminescence spectroscopy of mass-selected flavin adenine dinucleotide mono-anions

We report the absorption profile of isolated Flavin Adenine Dinucleotide (FAD) mono-anions recorded using photo-induced dissociation action spectroscopy. In this charge state, one of the phosphoric acid groups is deprotonated and the chromophore itself is in its neutral oxidized state. These measurements cover the first four optical transitions of FAD with excitation energies from 2.3 to 6.0 eV (210-550 nm). The S₀ → S₂ transition is strongly blue shifted relative to aqueous solution, supporting the view that this transition has a significant charge-transfer character. The remaining bands are close to their solution-phase positions. This confirms that the large discrepancy between quantum chemical calculations of vertical transition energies and solution-phase band maxima cannot be explained by solvent effects. We also report the luminescence spectrum of FAD mono-anions in vacuo. The gas-phase Stokes shift for S₁ is 3000 cm^-1, which is considerably larger than any previously reported for other molecular ions and consistent with a significant displacement of the ground and excited state potential energy surfaces. Consideration of the vibronic structure is thus essential for simulating the absorption and luminescence spectra of flavins.

Elucidation of the intrinsic optical properties of hydrogen-bonded and protonated flavin chromophores by photodissociation action spectroscopy

The intrinsic optical properties of the flavin chromophore when engaged in hydrogen bonding or being protonated were elucidated by photo-induced action spectroscopy and computations.

Advances in Quantum Mechanochemistry: Electronic Structure Methods and Force Analysis

In quantum mechanochemistry, quantum chemical methods are used to describe molecules under the influence of an external force. The calculation of geometries, energies, transition states, reaction rates, and spectroscopic properties of molecules on the force-modified potential energy surfaces is the key to gain an in-depth understanding of mechanochemical processes at the molecular level. In this review, we present recent advances in the field of quantum mechanochemistry and introduce the quantum chemical methods used to calculate the properties of molecules under an external force. We place special emphasis on quantum chemical force analysis tools, which can be used to identify the mechanochemically relevant degrees of freedom in a deformed molecule, and spotlight selected applications of quantum mechanochemical methods to point out their synergistic relationship with experiments.

Revisiting Vertical Models To Simulate the Line Shape of Electronic Spectra Adopting Cartesian and Internal Coordinates

Vertical models for the simulation of spectroscopic line shapes expand the potential energy surface (PES) of the final state around the equilibrium geometry of the initial state. These models provide, in principle, a better approximation of the region of the band maximum. At variance, adiabatic models expand each PES around its own minimum. In the harmonic approximation, when the minimum energy structures of the two electronic states are connected by large structural displacements, adiabatic models can breakdown and are outperformed by vertical models. However, the practical application of vertical models faces the issues related to the necessity to perform a frequency analysis at a nonstationary point. In this contribution we revisit vertical models in harmonic approximation adopting both Cartesian (x) and valence internal curvilinear coordinates (s). We show that when x coordinates are used, the vibrational analysis at nonstationary points leads to a deficient description of low-frequency modes, for which spurious imaginary frequencies may even appear. This issue is solved when s coordinates are adopted. It is however necessary to account for the second derivative of s with respect to x, which here we compute analytically. We compare the performance of the vertical model in the s-frame with respect to adiabatic models and previously proposed vertical models in x- or Q₁-frame, where Q₁ are the normal coordinates of the initial state computed as combination of Cartesian coordinates. We show that for rigid molecules the vertical approach in the s-frame provides a description of the final state very close to the adiabatic picture. For sizable displacements it is a solid alternative to adiabatic models, and it is not affected by the issues of vertical models in x- and Q₁-frames, which mainly arise when temperature effects are included. In principle the G matrix depends on s, and this creates nonorthogonality problems of the Duschinsky matrix connecting the normal modes of initial and final states in adiabatic approaches. We highlight that such a dependence of G on s is also an issue in vertical models, due to the necessity to approximate the kinetic term in the Hamiltonian when setting up the so-called GF problem. When large structural differences exist between the initial and the final-state minima, the changes in the G matrix can become too large to be disregarded.

General formulation of vibronic spectroscopy in internal coordinates

Our general platform integrating time-independent and time-dependent evaluations of vibronic effects at the harmonic level for different kinds of absorption and emission one-photon, conventional and chiral spectroscopies has been extended to support various sets of internal coordinates. Thanks to the implementation of analytical first and second derivatives of different internal coordinates with respect to cartesian ones, both vertical and adiabatic models are available, with the inclusion of mode mixing and, possibly, Herzberg-Teller contributions. Furthermore, all supported non-redundant sets of coordinates are built from a fully automatized algorithm using only a primitive redundant set derived from a bond order-based molecular topology. Together with conventional stretching, bending, and torsion coordinates, the availability of additional coordinates (including linear and out-of-plane bendings) allows a proper treatment of specific systems, including, for instance, inter-molecular hydrogen bridges. A number of case studies are analysed, showing that cartesian and internal coordinates are nearly equivalent for semi-rigid systems not experiencing significant geometry distortions between initial and final electronic states. At variance, delocalized (possibly weighted) internal coordinates become much more effective than their cartesian counterparts for flexible systems and/or in the presence of significant geometry distortions accompanying electronic transitions.

Insights into the Photoelectron Spectroscopy of Chlorofluoroethenes Studied by Density-Functional and Coupled-Cluster Theories

The first two ionic states of chlorofluoroethenes were studied by using both time-independent and time-dependent density-functional theories. We calculated the equilibrium geometries and harmonic vibrational frequencies of 1,1-, cis-, and trans-C2H2FCl and their cations by using the B3LYP and B3PW91 functionals together with the cc-pVTZ and aug-cc-pVTZ basis sets. Franck-Condon factors were computed by the method developed in our group, in which the Duschinsky effect was treated explicitly. A new technique, named alignment transformation, followed by Euler transformations was developed to achieve the Eckart conditions. The adiabatic ionization energies were calculated by the CCSD(T) method extrapolated to the complete basis set limit. Insights into the simulated photoelectron spectra of C2H2FCl indicate that the resolutions of recent threshold photoelectron experiments are not high enough to detect individual transitions. The high-resolution photoelectron spectra of C2H2FCl are predicted for future reference. The computed adiabatic ionization energies of the three isomers of C2H2FCl are in accord with the experiments with the absolute deviations ranging from 0.004 to 0.021 eV. We suggest that the agreement between experimental and theoretical spectra should be a key criterion to judge whether a spectral assignment is reasonable.

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.

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.

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.