Chapter 10

Rydberg and Continuum States; Fano Resonances

10.1Introduction

Every chapter so far has been preparation for this one.

We have a molecule (Chapter 7) with several ionization thresholds (Chapter 9), supporting Rydberg series that converge to the higher ones and therefore lie above the lower ones. Such a state is bound to one ionic core and degenerate with the continuum of another, and the electron–electron repulsion connects them. The tools are in place: continuum normalization (Chapter 3), Fermi's golden rule (Chapter 6), and the \((n^{*})^{-3}\) scaling of amplitude near the core (Chapter 4).

This chapter produces three things:

  1. the width \(\Gamma\) of an autoionizing state and hence its lifetime;

  2. the Fano profile, the asymmetric line shape produced when a resonant and a direct route lead to the same final state, together with the asymmetry parameter \(q\);

  3. the scaling laws that parameterise a whole Rydberg series from measured level positions and a single anchor.

By the end you will reproduce numbers from a published research table by hand. Chapter 11 then makes the whole thing time-dependent.

10.2Continuum Normalization

A bound state is normalizable, \(\braket{\psi}{\psi}=1\); a free electron is not. Something must replace normalization, and the choice matters because factors of \(2\pi\) and densities of states hide inside it.

Definition 10.2.1 - Energy normalization

Continuum states labelled by energy are normalized to a Dirac delta,

\begin{equation} \label{eq:FANO-energy-norm} \braket{\varepsilon}{\varepsilon'}=\delta(\varepsilon-\varepsilon'). \end{equation}

Two consequences are used constantly.

Theorem 10.2.1 - Mixed completeness

For a system with both bound and free states,

\begin{equation} \label{eq:FANO-completeness} \ket{\Psi}=\sum_n c_n\ket{n} +\int d\varepsilon\;c(\varepsilon)\ket{\varepsilon}, \qquad \sum_n|c_n|^{2}+\int d\varepsilon\,|c(\varepsilon)|^{2}=1 . \end{equation}

Equation (10.2) is the structure of the state vector in Chapter 11, and Def. 10.2.1 is what makes the two terms dimensionally compatible: \(|c(\varepsilon)|^{2}\) is a probability density per unit energy, while \(|c_n|^{2}\) is a probability.

The second consequence is Remark 6.5.1 of Chapter 6: the density of states is absorbed into the states themselves, so in atomic units Fermi's golden rule reduces to

\begin{equation} \label{eq:FANO-golden-normalized} \Gamma=2\pi\bigl|\bra{\varepsilon}\hat V\ket{i}\bigr|^{2} . \end{equation}

That is the reason for the convention: one fewer factor to get wrong.

10.3A Discrete State Coupled to a Continuum

Here is the model, stripped to essentials. Two ingredients: a discrete state \(\ket{a}\) at energy \(E_a\) — an autoionizing Rydberg state — and a continuum \(\ket{\varepsilon}\) spanning \(E_a\), normalized as in Def. 10.2.1. Neither is an eigenstate of the full Hamiltonian, because a coupling connects them:

\begin{equation} \label{eq:FANO-coupling} V_\varepsilon=\bra{\varepsilon}\hat H\ket{a} . \end{equation}

Physically \(\hat H\) here is the electron–electron repulsion of Chapter 5, called configuration interaction in this context.

The model Hamiltonian is therefore

\begin{equation} \label{eq:FANO-hamiltonian} \hat H=E_a\ket{a}\bra{a} +\int d\varepsilon\;\varepsilon\ket{\varepsilon}\bra{\varepsilon} +\int d\varepsilon\Bigl[V_\varepsilon\ket{\varepsilon}\bra{a} +V_\varepsilon^{*}\ket{a}\bra{\varepsilon}\Bigr]. \end{equation}
Theorem 10.3.1 - Autoionization width

Because the coupling is to a continuum rather than to a single state, \(\ket{a}\) does not shift and split — it decays, with

\begin{equation} \label{eq:FANO-width} \boxed{\;\Gamma=2\pi|V_\varepsilon|^{2}\Big|_{\varepsilon=E_a}\;} \end{equation}

and lifetime

\begin{equation} \label{eq:FANO-lifetime} \tau=\frac{\hbar}{\Gamma}, \qquad\text{i.e.}\qquad \Gamma\tau=\hbar . \end{equation}

Apply Fermi's golden rule in the form Eq. (10.3), with \(\hat V\) the coupling term of Eq. (10.5) and the final states the continuum at the same energy. Remark 6.5.2 of Chapter 6 established that the golden rule holds for any coupling, not only a radiative one.

In practical units, Eq. (10.7) reads

\begin{equation} \label{eq:FANO-units} \boxed{\;\Gamma\,[\mathrm{meV}]\times\tau\,[\mathrm{fs}]=658\;} \end{equation}

so a state of width \(37\) meV lives \(18\) fs (helium's \(2s2p\) resonance from Chapter 5), and a state living a picosecond has a width under a microelectronvolt. Widths and lifetimes are the same information.

It is often convenient to read Eq. (10.6) backwards: the width is what experiments measure, the coupling is what a calculation needs. Taking the coupling real and constant across the narrow resonance,

\begin{equation} \label{eq:FANO-effective-coupling} \boxed{\;\bra{\varepsilon}\hat H\ket{a}=\sqrt{\frac{\Gamma}{2\pi}}\;} \end{equation}

This expression appears verbatim in the research equations of Chapter 11.

10.4Fano's Diagonalisation

Now solve Eq. (10.5) exactly. The eigenstates of the full Hamiltonian are not \(\ket{a}\) and \(\ket{\varepsilon}\) but mixtures of them. Write the eigenstate at energy \(E\) as

\begin{equation} \label{eq:FANO-ansatz} \ket{\Psi_E}=\alpha(E)\ket{a} +\int d\varepsilon\;\beta_E(\varepsilon)\ket{\varepsilon}, \end{equation}

and impose \(\hat H\ket{\Psi_E}=E\ket{\Psi_E}\).

Projecting onto \(\bra{a}\) and onto \(\bra{\varepsilon}\) gives two equations:

\begin{align} \label{eq:FANO-eq-a} E_a\alpha+\int d\varepsilon\;V_\varepsilon^{*}\beta_E(\varepsilon)&=E\alpha,\\ \label{eq:FANO-eq-eps} \varepsilon\,\beta_E(\varepsilon)+V_\varepsilon\alpha &=E\,\beta_E(\varepsilon). \end{align}

Equation (10.12) rearranges to \((E-\varepsilon)\beta_E(\varepsilon)=V_\varepsilon\alpha\), which cannot simply be divided through because \(E-\varepsilon\) vanishes somewhere in the integration range. The correct general solution of such an equation adds a homogeneous piece proportional to \(\delta(E-\varepsilon)\) and interprets the quotient as a Cauchy principal value:

\begin{equation} \label{eq:FANO-beta} \beta_E(\varepsilon) =\left[\mathcal{P}\frac{1}{E-\varepsilon} +z(E)\,\delta(E-\varepsilon)\right]V_\varepsilon\,\alpha , \end{equation}

with \(z(E)\) fixed by substituting back into Eq. (10.11):

\begin{equation} \label{eq:FANO-z} z(E)=\frac{E-E_a-F(E)}{|V_E|^{2}}, \qquad F(E)=\mathcal{P}\!\int d\varepsilon\, \frac{|V_\varepsilon|^{2}}{E-\varepsilon} . \end{equation}

The function \(F(E)\) is a level shift: it displaces the resonance from \(E_a\) to \(E_a+F\), and in practice it is absorbed into the measured resonance position. What remains is the interesting part.

Theorem 10.4.1 - The modified discrete state

The admixture of the discrete state in the exact eigenstate at energy \(E\) is

\begin{equation} \label{eq:FANO-alpha} |\alpha(E)|^{2} =\frac{|V_E|^{2}} {\bigl[E-E_a-F(E)\bigr]^{2}+\pi^{2}|V_E|^{4}} =\frac{\Gamma/2\pi}{(E-E_r)^{2}+\Gamma^{2}/4}, \end{equation}

a Lorentzian of full width \(\Gamma=2\pi|V_E|^{2}\) centred on the shifted position \(E_r=E_a+F\).

So the discrete state is not destroyed — it is spread over a range of exact eigenstates of width \(\Gamma\). That is the spectral statement of a finite lifetime, and Eq. (10.15) is where Theorem 10.3.1 comes from in an exact treatment rather than a perturbative one.

10.5The Fano Profile

Now add the laser. Light can reach the continuum from the ground state \(\ket{g}\) in two ways:

  1. directly, \(\ket{g}\to\ket{\varepsilon}\), with amplitude \(\bra{\varepsilon}\hat D\ket{g}\);

  2. resonantly, \(\ket{g}\to\ket{a}\to\ket{\varepsilon}\), via the discrete state, which then autoionizes.

Both end in the same final state at the same energy, so they are indistinguishable and their amplitudes add before squaring. The cross term is the entire Fano effect.

The transition rate is \(|\bra{\Psi_E}\hat D\ket{g}|^{2}\), and from Eq. (10.10) the matrix element has one piece from each route. Carrying the algebra through with Eqs. (10.13) and (10.14) gives Fano's 1961 result.

Theorem 10.5.1 - Fano profile
\begin{equation} \label{eq:FANO-profile} \boxed{\; \sigma(E)=\sigma_b\,\frac{(q+\epsilon)^{2}}{1+\epsilon^{2}}\;} \end{equation}

where \(\sigma_b\) is the background cross section from the direct route alone,

\begin{equation} \label{eq:FANO-reduced-energy} \epsilon=\frac{E-E_r}{\Gamma/2}=\frac{2(E-E_r)}{\Gamma} \end{equation}

is the detuning in units of the half-width, and

\begin{equation} \label{eq:FANO-q-def} q=\frac{\bra{\tilde a}\hat D\ket{g}} {\pi V_E^{*}\bra{\varepsilon}\hat D\ket{g}} \end{equation}

is the asymmetry parameter: the ratio of the amplitude for exciting the (shifted, “dressed”) discrete state \(\ket{\tilde a}\) to that for direct ionization.

Read the formula by taking limits.

Two limits of \(q\) deserve names:

So \(q\) interpolates between a peak and a hole. Nothing about the state changes; only the balance of the two routes to it.

Remark 10.5.1 - Asymmetry means interference

A symmetric line — Lorentzian or Gaussian — can always be read as a rate: something decays, and the spectrum is the Fourier transform of an exponential. An asymmetric line cannot. Asymmetry requires an amplitude whose phase varies across the resonance, interfering with something else. When you see an asymmetric profile in a spectrum, you are looking at interference, and the asymmetry measures a relative phase.

Example 10.5.1 - Extracting \(\Gamma\) and \(q\) from a measured line

A photoabsorption line shows a minimum at \(59.9\) eV and a maximum at \(60.2\) eV, with the peak four times the far-off background. Find \(q\), \(E_r\) and \(\Gamma\).

Solution.

The peak-to-background ratio is \(1+q^{2}=4\), so \(q^{2}=3\) and \(q=\pm\sqrt{3}=\pm1.73\). Take \(q>0\) for now.

The zero is at \(\epsilon=-q\) and the maximum at \(\epsilon=1/q\), so in energy

\[ E_{\min}=E_r-\frac{q\Gamma}{2}, \qquad E_{\max}=E_r+\frac{\Gamma}{2q} . \]

Subtracting,

\[ E_{\max}-E_{\min}=\frac{\Gamma}{2}\left(q+\frac{1}{q}\right) =\frac{\Gamma}{2}\cdot\frac{q^{2}+1}{q} =\frac{\Gamma}{2}\cdot\frac{4}{1.73}=1.156\,\Gamma . \]

With \(E_{\max}-E_{\min}=0.3\) eV,

\[ \Gamma=\frac{0.3}{1.156}=0.26~\mathrm{eV}, \]

and by Eq. (10.8) the lifetime is \(\tau=658/260=2.5\) fs. Then

\[ E_r=E_{\min}+\frac{q\Gamma}{2}=59.9+\frac{1.73\times0.26}{2}=60.1~\mathrm{eV}. \]

Note that the resonance position is at neither the peak nor the dip: it lies between them. Quoting the peak as the resonance energy is a common and avoidable error, and for small \(q\) it can be badly wrong.

10.6A Rydberg Series of Autoionizing States

Put the resonance into a Rydberg series and everything becomes predictable, because Chapter 4 already established how Rydberg quantities scale.

Autoionization happens at the core: the Rydberg electron must be close in for \(\hat V_{ee}\) to transfer energy. Chapter 4's central result, Eq. (4.62), was that the probability of being near the core falls as \((n^{*})^{-3}\). Combining with Theorem 10.3.1:

Theorem 10.6.1 - Rydberg scaling of widths and dipoles
\begin{equation} \label{eq:FANO-gamma-scaling} \Gamma_n\propto|V_\varepsilon|^{2}\propto (n^{*})^{-3}, \qquad \tau_n\propto (n^{*})^{3}, \end{equation}

and, since a dipole matrix element is linear in the wave function where the width is quadratic,

\begin{equation} \label{eq:FANO-dipole-scaling} \bra{\mathcal R_n}\hat D\ket{\varepsilon}\propto (n^{*})^{-3/2}. \end{equation}

But \(n^{*}\) is not what an experiment measures — level positions are. Eliminate it using \(n^{*}=[\mathcal{E}_{\mathrm{Ryd}}/(E_{\mathrm{th}}-E_n)]^{1/2}\) from Chapter 5, Eq. (5.18). Since \((n^{*})^{3}\propto(E_{\mathrm{th}}-E_n)^{-3/2}\),

\begin{equation} \label{eq:FANO-tau-working} \boxed{\; \tau_n=\tau_{\mathrm{ref}} \left(\frac{E_{\mathrm{th}}-E_{\mathrm{ref}}} {E_{\mathrm{th}}-E_n}\right)^{3/2}\;} \end{equation}

and since \((n^{*})^{-3/2}\propto(E_{\mathrm{th}}-E_n)^{3/4}\),

\begin{equation} \label{eq:FANO-dipole-working} \boxed{\; \bra{\mathcal R_n}\hat D\ket{\varepsilon} =\bra{\mathcal R_{\mathrm{ref}}}\hat D\ket{\varepsilon} \left(\frac{E_{\mathrm{th}}-E_n} {E_{\mathrm{th}}-E_{\mathrm{ref}}}\right)^{3/4}\;} \end{equation}

Notice what has dropped out: the quantum defect. It cancels between the two applications of the \(n^{*}\) formula. Measured level positions plus one anchor give every other member's lifetime and coupling, with no structure calculation at all.

Example 10.6.1 - Reproducing a published parameter table

The CO\(_2\) Henning sharp series, \((3\sigma_u)^{-1}nd\sigma_g\), converges to \(E_B=18.08\) eV. The \(5d\) member lies at \(17.5395\) eV with a measured lifetime \(372.6\) fs. Predict the lifetime and relative infrared dipole of the \(4d\) member at \(17.2831\) eV.

Solution.

Binding energies below threshold:

\[ E_B-E_{5d}=18.08-17.5395=0.5405~\mathrm{eV}, \qquad E_B-E_{4d}=18.08-17.2831=0.7969~\mathrm{eV}. \]

Before computing, predict the direction. The \(4d\) state is more deeply bound, hence more compact, hence closer to the core — so it should autoionize faster and couple to the continuum more strongly.

Lifetime, from Eq. (10.21):

\[ \tau_{4d}=372.6\left(\frac{0.5405}{0.7969}\right)^{3/2} =372.6\times(0.6782)^{3/2} =372.6\times0.5585=208.1~\mathrm{fs}. \]

Dipole, from Eq. (10.22), taking the \(5d\) value as \(1\):

\[ \bra{\mathcal R_{4d}}\hat D_{\mathrm{IR}}\ket{\varepsilon} =\left(\frac{0.7969}{0.5405}\right)^{3/4}=(1.4744)^{3/4}=1.338 . \]

Both agree with the published values to every digit quoted. The lifetime fell and the dipole rose, as predicted.

Two arithmetic warnings. The exponents \(\tfrac32\) and \(\tfrac34\) come with the ratio the other way up — invert one and you get a plausible but wrong answer, the commonest error here. And as a sanity check, neighbouring members near \(n^{*}\approx5\) should differ by factors of order \((5/4)^{3}\approx2\), which \(0.56\) and \(1.34\) bracket.

10.6.1How \(q\) Varies Along a Series

From Eq. (10.18), \(q\) is the ratio of resonant to direct excitation amplitude. The denominator involves the continuum coupling, fixed by the width through Eq. (10.9), so

\begin{equation} \label{eq:FANO-q-scaling} q_n\simeq \frac{\bra{g}\hat D_{\mathrm{XUV}}\ket{\mathcal R_n}} {\sqrt{\pi\Gamma_n/2}\;\bra{g}\hat D_{\mathrm{XUV}}\ket{\varepsilon}} \;\propto\; \bra{g}\hat D_{\mathrm{XUV}}\ket{\mathcal R_n}\sqrt{\tau_n} . \end{equation}

A longer-lived state — a narrower resonance — therefore has a larger \(q\), hence a more peak-like, less asymmetric line. That is intuitive in hindsight: a narrow resonance concentrates the resonant amplitude into a small energy range where it dominates the smooth background.

Notice which quantity obeys no scaling law: the bound–bound dipole \(\bra{g}\hat D_{\mathrm{XUV}}\ket{\mathcal R_n}\) in the numerator. It depends on overlap with the ground state and on core details that no \(n^{*}\) power captures, and in practice it is read off measured photoabsorption peak heights state by state. Knowing which quantities scale and which must be measured is most of the skill in using a model like this.

10.7Two Series, and Why They Differ

CO\(_2\) has two Rydberg series converging to \(B^{2}\Sigma_u^{+}\): the Henning sharp series \((3\sigma_u)^{-1}nd\sigma_g\) and the Henning diffuse series \((3\sigma_u)^{-1}ns\sigma_g\). They interleave in the same energy range, autoionize into the same coupled continuum (\(\mathcal C\), here \(X^{2}\Pi_g\)), and have lifetimes differing by an order of magnitude — hundreds to thousands of femtoseconds for the sharp states, tens for the diffuse.

Everything needed is already in hand. From Chapter 5, the quantum defects: \(\delta_d\approx0\) and \(\delta_s\approx1.3\), because the centrifugal barrier of Chapter 3 keeps a \(d\) electron out of the core while an \(s\) electron dives through it. From this chapter, the width is proportional to the probability of being at the core. An \(s\) Rydberg electron is there far more often, so it autoionizes far faster.

One mechanism — penetration — fixes the quantum defects, the autoionization widths and the asymmetry parameters of both series.

The names are historical and inverted: the sharp series is the \(d\) series, because a long-lived state gives a narrow line by \(\Gamma\tau=\hbar\). Sharp means slow.

10.8Summary

10.9Exercises

  1. Widths and lifetimes. Verify Eq. (10.8) from \(\Gamma\tau=\hbar\) with \(\hbar=0.658\) eV fs. Then: a state of width \(1.8\) meV lives how long? A state living \(372.6\) fs has what width, in meV and \(\mu\)eV?

  2. Anatomy of the profile. For Eq. (10.16):

    1. Confirm the zero at \(\epsilon=-q\) and the value \(q^{2}\) at \(\epsilon=0\).

    2. Carry out the differentiation in Sec. 10.5 and confirm the maximum is at \(\epsilon=1/q\) with value \(1+q^{2}\).

    3. Sketch the profile for \(q=0\), \(1\), \(3\) and \(q\to\infty\) on one set of axes.

  3. Extracting parameters. Repeat Example 10.5.1 for a line whose minimum is at \(60.05\) eV, maximum at \(60.20\) eV, with peak-to-background \(8\). Find \(q\), \(\Gamma\), \(\tau\) and \(E_r\).

  4. Building a parameter table. Using \(E_B=18.08\) eV, the \(5d\) anchor (\(E=17.5395\) eV, \(\tau=372.6\) fs, dipole \(\equiv1.000\)) and Eqs. (10.21)(10.22), compute \(n^{*}\), \(\tau_n\) and the relative dipole for the sharp-series members at \(17.8650\) eV (\(8d\)) and \(17.9407\) eV (\(10d\)). Then, with \(q_{5d}=2.00\) and XUV dipoles \(2.0196\times10^{-2}\) (\(5d\)) and \(1.3687\times10^{-2}\) (\(8d\)), use Eq. (10.23) to predict \(q_{8d}\). The published value is \(2.71\); if you disagree, find the arithmetic error rather than adjusting anything.

  5. The diffuse series. Repeat for the diffuse series anchored at \(6s\) (\(E=17.4579\) eV, \(\tau=33.0\) fs) for the \(9s\) state at \(17.8457\) eV. Then answer: \(9s\) and \(9d\) lie only \(0.06\) eV apart yet their lifetimes differ by more than tenfold. Which quantity in Eq. (10.21) carries the difference, and which does not? Explain the physical origin in two sentences.

10.10Project: Fano Profiles from a Measured Series

The problem. Build the absorption spectrum of a whole Rydberg series of autoionizing states from measured level positions and one anchor, and watch the interference structure emerge.

The data are the CO\(_2\) Henning sharp series converging to \(E_B=18.08\) eV:

State\(4d\)\(5d\)\(6d\)\(7d\)\(8d\)\(10d\)
\(E_n\) (eV)17.283117.539517.701817.800017.865017.9407

with \(\tau_{5d}=372.6\) fs measured and \(q_{5d}=2.00\) as anchor.

  1. On paper. Compute \(n^{*}\) for each state from Eq. (5.18), hence the quantum defect. Confirm \(\delta_d\) is small and roughly constant — that constancy is Chapter 5's claim, and this is the test of it.

  2. On paper. Using Eq. (10.21), compute \(\tau_n\) for all six, then convert to widths \(\Gamma_n\) in meV using Eq. (10.8).

  3. On the computer. Plot \(\Gamma_n\) against \(n^{*}\) on log–log axes and fit a power law. Report the fitted exponent; it should be \(-3\).

  4. On the computer. Take \(q_n\) from Eq. (10.23), assuming for simplicity a common XUV bound dipole — an approximation the real calculation avoids, and one you should comment on. Plot the total spectrum

    \[ \sigma(E)=\sigma_b\sum_n\frac{(q_n+\epsilon_n)^{2}}{1+\epsilon_n^{2}}, \qquad \epsilon_n=\frac{2(E-E_n)}{\Gamma_n}, \]

    over \(17.2\) to \(18.0\) eV. Work out how many energy points are needed to resolve the narrowest resonance before plotting.

  5. On the computer. Zoom in on one resonance and mark the zero at \(\epsilon_n=-q_n\), the peak at \(\epsilon_n=1/q_n\) and the centre \(E_n\). Then vary \(q\) for that resonance from \(0\) to \(10\) and show the shape morphing from a window to a peak.

  6. What has to move. One animation: sweep \(q\) upward from \(0\) for the whole series at once, so every line turns from a hole into a peak together. Caption it with what to watch — the zero crossing sliding through the resonance as \(q\) grows.

The check. Part (c) is the real test, and it is a consistency check rather than a measurement: since (b) was derived from a \((n^{*})^{-3}\) law, recovering \(-3\) confirms the arithmetic, not the physics. Say so plainly in your write-up.

Be ready to answer. Your spectrum shows the higher members crowding together as they approach \(18.08\) eV while their widths shrink. Both trends are powers of \(n^{*}\). Which power governs each — and what therefore happens to the ratio of width to spacing as \(n\) grows? What does that tell you about whether neighbouring resonances can be treated as independent, an assumption the research model makes explicitly?