Q-Chem features a new module for extended excited-state analysis, which is interfaced to the ADC, CC/EOM-CC, CIS, and TDDFT/SF-TDDFT methods.760, 762, 759, 39, 761, 654 These analyses are based on the state, transition and difference density matrices of the excited states; the theoretical background for such analysis is given in Chapter 7.15.
||Largest NTO occupation numbers|
||, sum of NTO occupation numbers|
||Energy of hole NTO,|
||Energy of particle NTO,|
||NTO participation ratio|
||Replace in the above two formulas|
||Expec. value of the particle-hole permutation operator,|
||Mean position of hole|
||Mean position of electron|
||Linear e/h distance|
||RMS hole size:|
||RMS elec. size:|
The transition-density (1TDM) based analyses include the construction and export of natural transition orbitals636 (NTOs) and electron and hole densities,762 the evaluation of charge transfer numbers,760 an analysis of exciton multipole moments,39, 761, 654 and quantification of electron-hole entanglement.763 NTOs are obtained by singular value decomposition (SVD) of the 1TDM:
where is diagonal matrix containing singular values and unitary matrices and contain the respective particle and hole NTOs. Note that:
Furthermore, the formation and export of state-averaged NTOs, and the decomposition of the excited states into transitions of state-averaged NTOs are implemented.762 The difference and/or state densities can be exported themselves, as well as employed to construct and export natural orbitals, natural difference orbitals, and attachment and detachment densities.354 Furthermore, two measures of unpaired electrons are computed.362 In addition, a Mulliken or Löwdin population analysis and an exciton analysis can be performed based on the difference/state densities. The main descriptors of the various analyses that are printed for each excited state are given in Tables 10.1 and 10.2. For a detailed description with illustrative examples, see Refs. 762 and 759.
||Number of unpaired electrons|
||Number of unpaired electrons|
||NO participation ratio|
||Promotion number and|
||D/A participation ratio and|
||Mean position of detachment density|
||Mean position of attachment density|
||Linear D/A distance|
||RMS size of detachment density|
||RMS size of attachment density|
To activate any excited-state analysis STATE_ANALYSIS has to be set
to TRUE. For individual analyses there is currently only a limited
amount of fine grained control. The construction and export of any type of
orbitals is controlled by MOLDEN_FORMAT to export the orbitals as
MolDen files, and NTO_PAIRS which specifies the number of
important orbitals to print (note that the same keyword controls the number of
natural orbitals, the number of natural difference orbitals, and the number of
NTOs to be printed). Setting MAKE_CUBE_FILES to TRUE
triggers the construction and export of densities in “cube file”
format378 (see Section 10.5.4 for details).
Activating transition densities in $plots will generate cube files for the
transition density, the electron density, and the hole density of the
respective excited states, while activating state densities or
attachment/detachment densities will generate cube files for the state density,
the difference density, the attachment density and the detachment density.
Setting IQMOL_FCHK = TRUE (equivalently, GUI = 2)
will export data to the “
.fchk” formatted checkpoint file, and
switches off the generation of cube files. The population analyses are
controlled by POP_MULLIKEN and LOWDIN_POPULATION. Setting
the latter to TRUE will enforce Löwdin population analysis to be
employed for regular populations as well as CT numbers,
while by default the Mulliken population analysis is used.
Any MolDen or cube files generated by the excited state analyses can be
found in the directory
plots in the job’s scratch directory. Their names
always start with a unique identifier of the excited state (the exact form of
this human readable identifier varies with the excited state method). The
names of MolDen files are then followed by either
_nto.mo depending on the type of orbitals they
contain. In case of cube files the state identifier is followed by
_hole for state, difference, transition, attachment,
detachment, electron, or hole densities, respectively. All cube files have the
.cube. In unrestricted calculations an additional part is added
to the file name before
.cube which indicates (
_b) spin. The only exception is the state density for which
_sd are added indicating the total or spin-density parts
of the state density.
Analysis of relaxed CIS/TDDFT densities can be triggered by
The corresponding output files are marked by
Computation of ESPs for state, transition, and electron/hole densities (see Ref. 475) can be
triggered by setting ESP_GRID=-3.
These are indicated by
_esp as part of the file name.
ctnum_*.om file created in the main directory serves as input
for a charge transfer number analysis, as explained, e.g., in
Refs. 760, 653. Use the
TheoDORE program (theodore-qc.sourceforge.net) to
create electron/hole correlation plots and to compute fragment based descriptors.
When doing excited-state calculations from an open-shell reference, libwfa will perform the analysis for both and transition densities. Make sure you look at the correct one. The way to figure it out is to remember that in open-shell references , e.g., in doublet references, the unpaired electron is and the hole is . Thus, for transitions of the unpaired electron into the unoccupied orbitals you need block, whereas for the transitions from doubly occupied orbitals into the singly un-occupied orbital (the hole) you need the block.
In Hermitian formalisms, is a Hermitian conjugate
of , but in non-Hermitian approaches, such as
coupled-cluster theory, the two are slightly different. While for quantitative
interstate properties both and are
computed, the qualitative trends in exciton properties derived from
and are very similar. Only
one 1TDM is analyzed for EOM-CC.
In spin-restricted calculations, the libwfa module computes NTOs for
the block of transition density. Thus, when computing NTOs for
the transitions between open-shell EOM-IP/EA states make sure to specify
correct spin states. For example, use EOM_EA_ALPHA to visualize
transitions involving the extra electron.
$rem basis = def2-sv(p) n_frozen_core = fc method = adc(2) ee_singlets = [0,1,1,0] state_analysis = true $end $molecule 0 1 6 0 0 0.523383 8 -0 0 -0.671856 1 0.931138 0 1.11728 1 -0.931138 0 1.11728 $end
$rem method = pbe0 basis = def2-sv(p) cis_n_roots = 4 cis_singlets = true cis_triplets = true rpa = false state_analysis = true molden_format = true nto_pairs = 3 make_cube_files = true esp_grid = -3 $end $molecule 0 1 6 1.19438 1.10251 -0 6 -0.00836561 1.69243 -0 7 -1.1696 0.978035 -0 6 -1.21206 -0.402293 -0 7 0.0346914 -0.97914 0 6 1.28159 -0.348737 0 8 -2.24342 -1.02375 -0 8 2.29918 -0.995854 0 1 -0.12316 2.76714 -0 1 -2.06144 1.4441 -0 1 0.0448178 -1.98999 0 1 2.10472 1.67984 -0 $end $plots Write cube files for all 4 states 70 -3.5 3.5 70 -3.5 3.5 30 -1.5 1.5 0 4 0 0 1 2 3 4 $end
Other examples of libwfa uses:
Example 7.3.7 illustrates wave-function analysis of the SF-DFT states in para-benzyne;
Example 126.96.36.199 illustrates wave-function analysis of XAS transitions within CVS-EOM-EE;
Example 7.10.26 illustrates wave-function analysis for transitions between EOM-IP states;
Example 7.10.16 illustrates wave-function analysis of complex-valued densities within CAP-EOM-CCSD.