UV/vis absorption/emission

UV/vis absorption/emission#

Generalized eigenvalue equation#

The standard method to calculate UV/vis absorption and emission spectra is to solve the generalized eigenvalue equation [NRS18].

In the case of SCF theory, it is commonly referred to as the time-dependent density functional theory or Hartree–Fock (TDDFT or TDHF) approach. TDHF is also known as the Random Phase Approximation (RPA).

If electron de-excitations are ignored in the formation of the electronic Hessian, then one arrives at the Tamm–Dancoff approximation and which can be invoked with a keyword in the input file.

VeloxChem implements a reduced-space Davidson algorithm to solve the equation for the N lowest eigenvalues (bottom-up). Based on these eigenvalues, or transition frequencies, and the associated transition moments, the dimensionless oscillator strengths are calculated according to

fn0=2meωn03ℏe2∑α=x,y,z|⟨0|μ^α|n⟩|2

With oscillator strengths and transition frequencies, the linear absorption cross section can be determined from the expression [NRS18]

σ(ω)=2π2e2ω(4πε0)mec∑n>0f(ω;ωn0,γ)fn0ωn0

where f is the Cauchy distribution.

@jobs
task: response
@end

@method settings
xcfun: b3lyp
basis: def2-svp
@end

@response
property: absorption
! tamm_dancoff: yes
nstates: 3
@end

@molecule
charge: 0
multiplicity: 1
xyz:  
...
@end

Complex polarization propagator approach#

The linear absorption cross section can be determined directly from the imaginary part of the polarizability [NRS18]

σ(ω)=ωϵ0cIm{α―(−ω;ω)}

where

α―=13(αxx+αyy+αzz)

and

ααβ(−ω;ω)=−⟨⟨μ^α;μ^β⟩⟩ωγ

The polarizability is complex and calculated with a damping term, ℏγ, associated with the inverse finite lifetime of the excited states. The default program setting for this parameter is 0.124 eV (or 0.004556 a.u.).

The resulting values for σ(ω) are presented in atomic units and can be converted to the SI unit of m2 by multiplying with a factor of a02.

The arbitrary frequency region is specified in the input file together with a requested frequency resolution.

@jobs
task: response
@end

@method settings
xcfun: b3lyp
basis: def2-svp
@end

@response
property: absorption (cpp)
! frequency region (and resolution)
frequencies: 0.0-0.15 (0.0025)
damping: 0.0045563  ! this is the default value
@end

@molecule
charge: 0
multiplicity: 1
xyz:  
...
@end
Contents