Density Matrices¶

After completing the CASSCF optimization, we have obtained a multiconfigurational wavefunction,

\[ \Psi=\sum_I c_I\Phi_I. \]

Although this wavefunction completely describes the electronic structure of the molecule, it is not convenient for calculating molecular properties directly.

Instead, quantum chemistry programs construct density matrices.

The density matrix provides a compact description of the electronic distribution and serves as the starting point for calculating almost every molecular property.

Why Do We Need Density Matrices?¶

Imagine trying to calculate the electron density directly from the complete CASSCF wavefunction.

For a calculation containing hundreds or thousands of Configuration State Functions,

Ψ

↓

Thousands of CSFs

↓

Very Complicated

Performing property calculations directly from the wavefunction would be extremely inefficient.

Instead,

the wavefunction is transformed into a much smaller mathematical object,

Wavefunction

↓

Density Matrix

↓

Properties

Almost all post-processing calculations use the density matrix rather than the wavefunction itself.

What Does the Density Matrix Describe?¶

The density matrix tells us

where electrons are located,
how electrons are distributed among molecular orbitals,
how strongly orbitals interact,
how electron correlation changes orbital occupations.

It is therefore one of the most information-rich quantities produced by a CASSCF calculation.

One-Particle Density Matrix¶

The most commonly used density matrix is the one-particle density matrix,

\[ \boxed{ D_{pq} = \langle \Psi | a_p^\dagger a_q | \Psi \rangle } \]

where

\(a_p^\dagger\) creates an electron in orbital \(p\),
\(a_q\) removes an electron from orbital \(q\),
\(D_{pq}\) measures the probability of finding an electron in orbital \(p\) while removing one from orbital \(q\).

Although this expression appears abstract, its interpretation is straightforward.

Physical Meaning¶

Each element of the density matrix describes the relationship between two molecular orbitals.

For example,

Orbital p

↓

Orbital q

↓

Electron Distribution

The matrix therefore summarizes how electrons occupy and interact between orbitals.

Matrix Representation¶

For four orbitals,

the density matrix has the form

\[ \mathbf D= \begin{pmatrix} D_{11}&D_{12}&D_{13}&D_{14}\\ D_{21}&D_{22}&D_{23}&D_{24}\\ D_{31}&D_{32}&D_{33}&D_{34}\\ D_{41}&D_{42}&D_{43}&D_{44} \end{pmatrix} \]

Each element has a different physical meaning.

Diagonal Elements¶

The diagonal elements,

\[ D_{11}, D_{22}, D_{33}, ... \]

represent the electron population of each orbital.

For example,

Orbital 1

Occupation ≈2

↓

Nearly Full

or

Orbital 8

Occupation ≈0

↓

Nearly Empty

These values provide direct information about electron occupation.

Off-Diagonal Elements¶

The off-diagonal elements,

\[ D_{12}, D_{13}, D_{24}, ... \]

measure the interaction between different orbitals.

Large values indicate

orbital mixing,
electron delocalization,
strong correlation between orbitals.

Small values indicate weak interaction.

Density Matrix in RHF¶

For a closed-shell RHF calculation,

the density matrix is very simple.

Occupied orbitals have

Occupation = 2

Virtual orbitals have

Occupation = 0

The density matrix therefore resembles

\[ \mathbf D= \begin{pmatrix} 2&0&0&0\\ 0&2&0&0\\ 0&0&0&0\\ 0&0&0&0 \end{pmatrix} \]

The occupations are fixed because RHF assumes a single electronic configuration.

Density Matrix in CASSCF¶

In CASSCF,

electrons are shared among many configurations.

As a result,

orbital occupations become fractional.

For example,

\[ \mathbf D= \begin{pmatrix} 1.98&0&0&0\\ 0&1.75&0&0\\ 0&0&0.31&0\\ 0&0&0&0.02 \end{pmatrix} \]

Notice that

the occupations are no longer limited to

0

or

2

Instead,

they reflect the multiconfigurational nature of the wavefunction.

Why Fractional Occupations Matter¶

Fractional occupations indicate that electrons are distributed over several configurations.

For example,

Occupation

2.00

↓

Fully Occupied

1.50

↓

Strongly Correlated

0.80

↓

Partially Occupied

0.02

↓

Nearly Empty

This behaviour is one of the clearest signatures of static electron correlation.

Properties Calculated from the Density Matrix¶

Once the density matrix has been constructed,

many molecular properties can be calculated directly.

These include

electron density,
Mulliken populations,
dipole moments,
atomic charges,
expectation values,
natural orbitals,
natural occupation numbers.

Thus,

Wavefunction

↓

Density Matrix

↓

Observable Properties

One-Particle vs Two-Particle Density Matrices¶

CASSCF also constructs the two-particle density matrix.

One-Particle Density

↓

Single Electron

Two-Particle Density

↓

Pairs of Electrons

The two-particle density matrix is required for

electron correlation energy,
orbital optimization,
analytical gradients,
post-CASSCF methods.

Although more computationally demanding,

it contains additional information about electron-electron interactions.

Connection with Orbital Optimization¶

During every CASSCF iteration,

the density matrices are recomputed.

CI Calculation

↓

Density Matrix

↓

Orbital Gradient

↓

Orbital Rotation

↓

Repeat

Thus,

the density matrix is an essential ingredient in the self-consistent optimization procedure.

Looking Ahead¶

Although the density matrix itself is useful,

its true power appears after diagonalization.

Diagonalizing the one-particle density matrix produces

natural orbitals,
natural occupation numbers,

which provide one of the most powerful tools for analysing multiconfigurational wavefunctions and assessing the quality of the chosen active space.

This is the subject of the next chapter.

Key Takeaways¶

The density matrix provides a compact description of the electronic structure.
Most molecular properties are calculated from the density matrix rather than directly from the wavefunction.
The diagonal elements correspond to orbital occupations.
The off-diagonal elements describe orbital interactions and electron delocalization.
RHF produces occupations of approximately 2 or 0, whereas CASSCF naturally yields fractional occupations.
Density matrices are fundamental for orbital optimization, property calculations, and post-CASSCF methods.