X

Search Results

Searching....

7.15 Visualization of Excited States

7.15.2 Attachment/Detachment Density Analysis

(April 13, 2024)

Consider the one-particle density matrices of the initial and final states of interest, 𝐏1 and 𝐏2 respectively. Assuming that each state is represented in a finite basis of spin-orbitals, such as the molecular orbital basis, and each state is at the same geometry. Subtracting these matrices yields the difference density

𝚫=𝐏1-𝐏2 (7.131)

Now, the eigenvectors of the one-particle density matrix 𝐏 describing a single state are termed the natural orbitals, and provide the best orbital description that is possible for the state, in that a CI expansion using the natural orbitals as the single-particle basis is the most compact. The basis of the attachment/detachment analysis is to consider what could be termed natural orbitals of the electronic transition and their occupation numbers (associated eigenvalues). These are defined as the eigenvectors 𝐔 defined by

𝐔𝚫𝐔=δ (7.132)

The sum of the occupation numbers δp of these orbitals is then

tr(𝚫)=p=1Nδp=n (7.133)

where n is the net gain or loss of electrons in the transition. The net gain in an electronic transition which does not involve ionization or electron attachment will obviously be zero.

The detachment density

𝐃=𝐔𝐝𝐔 (7.134)

is defined as the sum of all natural orbitals of the difference density with negative occupation numbers, weighted by the absolute value of their occupations where 𝐝 is a diagonal matrix with elements

dp=-min(δp,0) (7.135)

The detachment density corresponds to the electron density associated with single particle levels vacated in an electronic transition or hole density.

The attachment density

𝐀=𝐔𝐚𝐔 (7.136)

is defined as the sum of all natural orbitals of the difference density with positive occupation numbers where 𝐚 is a diagonal matrix with elements

ap=max(δp,0). (7.137)

The attachment density corresponds to the electron density associated with the single particle levels occupied in the transition or particle density. The difference between the attachment and detachment densities yields the original difference density matrix

𝚫=𝐀-𝐃. (7.138)

Within a CIS or TDDFT calculation, where the transitions are strictly one-electron in nature, the matrices 𝐀 and 𝐃 are the particle (virtual-virtual) and hole (occupied-occupied) components of the unrelaxed difference density matrix.