Excited States with TDDFT

Introduction

Color, photochemistry and fluorescence involve electronic states above the ground state. Time-dependent density-functional theory, usually abbreviated TDDFT, offers a practical route to excitation energies and transition intensities for many molecules. Its output may list dozens of states with energies, wavelengths and orbital transitions. The useful interpretation starts by asking which molecular geometry is held fixed, how strongly a transition couples to light, and whether the chosen functional can describe the state's electronic character. A TDDFT stick spectrum is a prediction under assumptions, not a direct replacement for a measured solution absorption spectrum.

Core explanation

In linear-response TDDFT, a small time-dependent perturbation probes how the ground-state electron density responds. The response equations yield excitation energies and transition properties within the chosen approximate exchange–correlation functional and kernel. A vertical absorption excitation is the difference between the ground state and an excited state at the ground-state nuclear geometry. Nuclei are effectively fixed during the electronic transition, consistent with the Franck–Condon picture. The energy may be reported in eV or converted approximately to a photon wavelength by λ(nm) ≈ 1240/E(eV). This wavelength conversion is arithmetic, not a guarantee that a visible band will be observed.

An oscillator strength f is a dimensionless measure linked to the transition dipole and excitation energy for an electric-dipole transition. A state with f near zero can be dark under that approximation even if its excitation energy lies in the visible range. A spectrum can contain vibronic intensity borrowing, solvent effects and broadening not represented by a simple vertical stick. To simulate a plotted absorption envelope, broaden and weight the computed transitions using an explicitly chosen line shape and width. Original benchmark work on TDDFT oscillator strengths shows that transition properties merit separate validation from energies.

An excited-state configuration is often described as a mixture of one-electron orbital promotions, such as HOMO → LUMO. Such labels can help identify π→π or lone-pair→π character, but a TDDFT state can mix several promotions; orbital order and character depend on the method. Analyze transition density or natural transition orbitals when a single pair of canonical orbitals does not explain the state clearly. The ground-state HOMO–LUMO gap is not the same as the response excitation energy, even within one calculation.

Standard local and semilocal approximations, and some global hybrids, can substantially underestimate long-range charge-transfer excitation energies. The underlying approximate functional and its long-range electron–hole interaction matter. Range-separated hybrids can improve such cases, though parameter choice and validation remain necessary. Primary long-range-corrected functional research explicitly evaluated charge-transfer as well as ordinary valence excitations. Rydberg states require diffuse basis functions and a suitable long-range potential. States with strong double-excitation character are poorly handled by common adiabatic linear-response TDDFT because its usual frequency-independent kernel does not generate the required excitation structure reliably. These limitations mean that a method successful for a local π→π transition may fail on another state in the same molecule.

Absorption and emission refer to different nuclear geometries. After excitation, a molecule may relax on an excited-state surface and then emit from that relaxed geometry. A vertical emission energy calculated there differs from a vertical absorption energy at the ground-state geometry. Solvent polarization can respond on different timescales to electronic and nuclear rearrangement, so nonequilibrium solvation treatment may matter. A state may also cross or mix with another during relaxation; following the state label by ordinal number alone can be unsafe when states reorder. Experimental band maxima include vibronic and environmental structure, so a vertical energy should be compared with an appropriate experimental feature, not automatically with any convenient peak.

Step-by-step reasoning

1. Optimize and verify the ground-state geometry, charge, spin and relevant conformers. 2. Select a functional and basis suited to the expected local, Rydberg or charge-transfer character; include diffuse functions where needed. 3. Compute enough excited states to cover the spectral region and inspect energies, oscillator strengths and transition character. 4. Convert energies to wavelength only with correct units and state that the values are vertical at the chosen geometry. 5. For plotted spectra, disclose broadening and, for solution comparisons, solvent and conformer treatment. 6. Investigate suspicious states or strong method sensitivity with benchmarks or a suitable wavefunction method.

Visual explanation

Draw two potential-energy curves against a nuclear coordinate: ground state below and excited state above. A vertical arrow from the ground-state minimum to the upper curve represents absorption. A second downward vertical arrow from the excited-state minimum represents emission. Beside the diagram, show three TDDFT sticks at different wavelengths: one tall for large f, one short for small f and one absent for f ≈ 0. Their heights portray electric-dipole intensity, not state populations.

Real-world analogy

Taking a fast photograph of a person mid-jump freezes their body position even though the jump continues afterward. Vertical excitation similarly changes the electronic state before nuclei appreciably rearrange. An emission event after relaxation starts from a different position on the excited-state surface. The analogy does not describe quantum transitions directly, but it clarifies why absorption and emission energies can differ without any arithmetic inconsistency.

Real-world example

A team designs an organic dye with a donor and acceptor connected by a conjugated bridge. A conventional TDDFT calculation predicts a low-energy visible transition with appreciable oscillator strength. Inspection of the transition density shows electron density moving from donor toward acceptor, suggesting charge-transfer character. The team checks a range-separated functional and solvent treatment rather than trusting the first wavelength. It also computes a relaxed excited-state geometry before discussing fluorescence. A successful prediction must be judged against a relevant experimental spectrum or high-level benchmark, not solely an appealing orbital picture.

Why?

Why can an optically dark state still matter chemically? Oscillator strength measures coupling to electric-dipole absorption from a particular initial state, not whether the excited state exists. A dark state can be populated through internal conversion, intersystem crossing, energy transfer or symmetry-breaking interactions. It may control photochemical pathways even if it barely appears in a direct absorption stick spectrum. Thus a low f should change how a state is sought experimentally, not erase it from mechanistic consideration.

Common misconception

“TDDFT wavelength equals the solution color exactly.” Real spectra reflect vibronic structure, solvent, conformers and method error. “The lowest excited state must have the strongest absorption” is false when its transition dipole is weak. “A HOMO–LUMO gap is a TDDFT result” confuses ground-state eigenvalues with response excitations. “A single functional works equally well for local, charge-transfer, Rydberg and double-excitation states” ignores distinct known failure modes.

Worked example

Suppose an illustrative TDDFT calculation gives a vertical state at 2.50 eV with f = 0.40 and another at 2.20 eV with f = 0.002. Their approximate wavelengths are 1240/2.50 = 496 nm and 1240/2.20 ≈ 564 nm. The lower-energy state does not necessarily dominate the visible absorption because its electric-dipole oscillator strength is much smaller. If a relaxed excited-state calculation gives an emission energy of 2.10 eV, its approximate photon wavelength is 590 nm, but comparing that number with an observed emission maximum still requires solvent and vibronic context. The example does not imply a particular dye or that peak height is exactly proportional to f after all experimental effects.

Quick check

1. What stays fixed in a vertical electronic excitation calculation? Answer: Nuclear geometry, at the selected ground-state configuration for vertical absorption. 2. Does f ≈ 0 mean an excited state cannot exist or participate in photochemistry? Answer: No. It means the modeled electric-dipole transition from the specified initial state is very weak.

Exam focus

Use λ(nm) ≈ 1240/E(eV) correctly and label a state as vertical at a specified geometry. Distinguish energy from oscillator strength, absorption from relaxed-state emission, and TDDFT response energy from an orbital gap. Identify at least one reason to test a charge-transfer or Rydberg assignment with a different method or larger basis. Interpret a stick spectrum with its broadening and experimental conditions stated.

Advanced insight

Excitation energies and oscillator strengths are different observables and can have different error patterns. Natural transition orbitals compress a many-orbital transition into dominant hole and particle patterns, but their interpretation still depends on the underlying TDDFT state. In a chromophore embedded in a protein or solid, polarization and electronic coupling to neighbors can shift or mix states beyond an isolated-molecule calculation. Nonadiabatic pathways require information about multiple surfaces and couplings, not merely a list of vertical energies.

Summary

TDDFT predicts electronic response at specified structures, yielding vertical excitation energies and transition strengths. A wavelength converted from energy locates a possible transition; oscillator strength helps indicate whether it may appear strongly in electric-dipole absorption. Functional, basis, solvent, geometry and state character determine reliability, with charge-transfer, Rydberg and double-excitation states requiring special scrutiny. Absorption and emission involve different geometries, so a useful spectral comparison must name the calculated observable.

Practice questions

1. What wavelength corresponds approximately to a 3.10 eV vertical excitation? Answer: 1240/3.10 ≈ 400 nm. 2. A 2.0-eV state has f = 0.0001 while a 2.4-eV state has f = 0.5. Which is likely stronger in a simple electric-dipole stick spectrum? Answer: The 2.4-eV state, despite its higher energy. 3. Why might a Rydberg excitation need diffuse basis functions? Answer: Its excited electron density extends far from the nuclei and requires spatially flexible functions. 4. Why is a vertical absorption energy generally different from a relaxed fluorescence emission energy? Answer: They are evaluated at different nuclear geometries and may have different solvent polarization states.