◼ Cyclic π-delocalization in metallacycle C7H7FeCl
This example comes from the original study
Electron delocalization in planar metallacycles: Hückel or Möbius aromatic?
D.W. Szczepanik (), M. Solà () ChemistryOpen 8 (2019) 219−227. DOI: 10.1002/open.201900014. URL
Abstract:
In this work the relationship between the formal number of π-electrons, d-orbital conjugation topology, π-electron delocalization and aromaticity in d-block metallacycles is investigated in the context of recent findings concerning the correlation of π-HOMO topology and the magnetic aromaticity indices in these species. It is demonstrated that for π-electron rich d-metallacycles the direct link between aromaticity, the number of π-electrons and the frontier π-orbital topology does not strictly hold and for such systems it is very difficult to unambiguously associate their aromaticity with the '4n+2' (Hückel) and '4n' (Möbius) rules. It is also shown that the recently proposed electron density of delocalized bonds (EDDB) method can successfully be used not only to quantify and visualize aromaticity in such difficult cases, but also – in contrast to magnetic aromaticity descriptors – to provide a great deal of information on the real role of d-orbitals in metallacycles without the ambiguity of bookkeeping of electrons in the π-subsystem of the molecular ring. Interestingly, some of the metallacycles studied cannot be classified exclusively as Hückel or Möbius because they have a hybrid Hückel-Möbius or even quasi-aromatic nature.
1. NBO7 archive workflow (legacy)
Note — the .49 lane is the legacy path. In earlier releases of runEDDB (and the original runEDDB.R script), the NBO7 archive .49 was the only supported input format. Since runEDDB-nvbao the in-house NAO driver reconstructs the NAO basis directly from .fchk / .json / .molden / .tmol / .gms inputs — no NBO7 installation is needed. The .49 lane is kept for backward compatibility and for users who already have NBO7 archives.
The .49 file is an NBO7 archive produced by post-processing a wavefunction from any QC program interfaced with NBO7 (Gaussian, ORCA, Q-Chem, GAMESS, Turbomole, …). In this tutorial we use Gaussian G16 as the wavefunction source. The flow is a single step:
Prepare a Gaussian input file with Pop(NBO7Read) (this invokes the externally-installed NBO7 rather than the older NBO 3.1 shipped inside Gaussian) and append the NBO-directive line that asks NBO7 to write the runEDDB-ready .49 archive directly:
# CAM-B3LYP/def2SVP Pop(NBO7Read)
C7H7FeCl
0 1
Fe -0.064251938266 0.000000000000 0.383381287583
C 0.102879051008 1.467437846806 -0.581782848285
C 0.102879051008 -1.467437846806 -0.581782848285
C 0.042058436784 1.945591795572 -1.877580721536
C 0.042058436784 -1.945591795572 -1.877580721535
C -0.012197738912 1.304915633429 -3.124178523721
C -0.012197738912 -1.304915633430 -3.124178523721
C -0.022766917881 0.000000000000 -3.622055704846
H -0.061079076601 0.000000000000 -4.715595953142
H 0.167516111578 2.281514394810 0.169634195203
H 0.167516111578 -2.281514394810 0.169634195204
H 0.029959638034 3.038668532273 -1.974142731372
H 0.029959638034 -3.038668532273 -1.974142731372
H -0.055473924074 2.032073158396 -3.942519373400
H -0.055473924074 -2.032073158396 -3.942519373400
Cl 1.502523783826 0.000000000000 1.822684376427
$NBO SKIPBO FILE=c7h7fecl AOINFO=49 AONAO=W49 DMNAO=W49 $END
The Pop(NBO7Read) keyword is essential: Pop(NBORead) (without the 7) invokes the much older NBO 3.1 module bundled inside Gaussian, whose archive precision and AOINFO conventions are not fully compatible with the runEDDB cube generator. The trailing $NBO … $END block instructs NBO7 to SKIPBO (skip the standard bond-orbital analysis — runEDDB does its own), write the AOINFO, AONAO (AO ↔ NAO transformation), and DMNAO (density in the NAO basis) into a single .49 archive named after FILE=c7h7fecl.
Run Gaussian (tested version: ES64L-G16RevC.01) with an NBO7-interfaced executable; this writes the runEDDB-ready c7h7fecl.49 directly:
terminal
$ g16 c7h7fecl.gjf
2. Global delocalization pattern
First, we are going to investigate the distribution of electron density of all delocalized bonds in the entire molecule using the EDDBG function (the subscript G stands for global: all atoms and bonds in the molecule are included). To do so, run the following command in the terminal:
terminal
$ runEDDB --input c7h7fecl.49 --quiet
The --quiet (or -q) option enables quiet mode: only the EDDBG summary populations are printed; useful for batch jobs where only the numerical result is needed:
output
> Printing summary of EDDB calculations:
Total ED_G population: 91.79267 per atom: 5.73704
Total EDLB_G population: 83.48845 per atom: 5.21803
Total EDDB_G population: 8.30422 per atom: 0.51901
> Calculation completed on 2026-05-11. Total time: 0s.
The default analysis basis is the Natural Valence Basis (NVB), which retains only the valence-shell Natural Atomic Orbitals (NAO). Of the ~49.8 valence electrons in the NVB, about 41.5e are localized (lone pairs and two-center bonds) and ~8.3e are shared through the system of delocalized bonds.
To quickly visualize the EDDBG-layer of the one-electron density, add the --output-eddb option:
This generates a Gaussian cube file (c7h7fecl.EDDB.cube) for visualization in, e.g., Avogadro2, IQmol, GaussView, VMD, or Chemcraft; a typical isosurface value is 0.015-0.020. The --ncores (or -n) option parallelizes the grid evaluation, which is the most compute-intensive step. It is recommended to use as many cores as available.
The grid quality is controlled by --cube-size (or -c) option (default: 100). Running without --quiet shows the cube export details:
output
> Exporting EDDB_G to c7h7fecl.EDDB.cube...
Grid scope: Global
Grid size: 73x120x115 (~100^3)
Progress: ######################################## 100%
Integrated density: 8.30501 (100.0% of 8.30422)
Time: 2s
The integrated density (~8.305e) almost perfectly matches the analytical EDDBG population (~8.304e) — this confirms the cube accurately represents the density. The .49 lane uses a parallel NBO-convention cube generator that consumes the AOINFO / AONAO transformation matrices supplied by NBO7 directly, so the grid integration agrees with the analytical trace to roughly the same precision as the native lanes.
To visualize the generated cube file (c7h7fecl.EDDB.cube) in Avogadro2, select in program menu Analyze → Create Surfaces and set parameters as in the figure below:
3. Detailed population analysis
Now, switch to the verbose mode using --verbose (or -v) option:
terminal
$ runEDDB --input c7h7fecl.49 --verbose
First, the program reports what was parsed from the wavefunction file:
output
> Loading and parsing input file... 0s
Source file: c7h7fecl.49
File format: NBO7 archive (.49)
Number of atoms: 16
Number of basis functions: 182 (spherical)
Number of contracted shells: 82
Number of primitives: 166
Number of molecular orbitals: 0 (.49 carries D in NAO basis directly)
Effective core potential: No
Charge / multiplicity: 0 / 1
Total electrons: 92
Unrestricted: No
Reference type: HF / DFT
Electron population (Tr[DMNAO]): 92.00000
Next, it summarises the internal AO→NAO localization procedure:
The .49 archive carries a flat NMB flag per NAO (0 = Rydberg, 1 = NMB) with no Cor / Val split, so the parser cannot separate the 21 core NAOs from the 45 valence NAOs — all 66 non-Rydberg NAOs are routed through NMB (Cor + Val) which carries ~91.8e. The 116 Rydberg NAOs are excluded as usual. To avoid silent confusion, runEDDB makes -b nmb the effective default for .49 input (the -b nvb default of the native lanes does not apply); a one-line clarifier appears in the Stage-2 detail block under Target NAO representation: in verbose mode.
In the next stage program prints the BOP configuration and performs calculations:
output
> Performing multicenter bond-orbital projections (BOP)...
Type of EDDB function: EDDB_G (entire molecular system)
BOP algorithm: Standard BOP (HF/DFT 1-densities)
Atoms included in BOP: 16
Pairs included in BOP: 69
Bond-order threshold: 0.00100
Progress: ######################################## 100%
Number of BOP cycles: 1
BOP time: 0s
In the entire system, 69 atom pairs exceed the Wiberg bond-order threshold (10−3) and enter the BOP cycle.
After the BOP step, program prints three electron-population tables. The first lists Natural Orbitals for Bond Delocalization (NOBD):
From the table above it is clear that the first four NOBDs make the dominant contributions to electron delocalization. Their orbital decomposition shows that the first three NOBDs have almost purely π-symmetry (~1.7-1.8e, 98-100% p-character), while the fourth involves the Fe atom through its 3d orbitals (~26% p- and ~74% d-character). The remaining NOBDs have much smaller occupancies associated with σ-delocalization tails.
The next table presents the population analysis in the atomic-orbital resolution — one row per NAO in the target representation (here NVB, valence-only):
The data reveal that the main delocalization channel involves the 3d2 orbital on Fe (NMB index 12), occupied by ~1.5e and contributing ~0.9e (~60%) to delocalization, together with the 2p1 orbitals on carbon atoms, occupied by about 0.8–1.0e and contributing about 0.6–0.9e to delocalization. Electrons on chlorine and hydrogen atoms are predominantly localized.
Finally, the last table collapses the atomic-orbital data to atomic resolution:
output
> Printing results of electron population analysis in atomic resolution...
+------------+-----------------------------------+-----------------------------------+
| Atom | Electron delocalization | Electron population |
+------------+-----------------------------------+-----------------------------------+
Index Sym Total Alpha Beta Total Alpha Beta
----- --- ========= --------- --------- --------- --------- ---------
1 Fe 1.34453 - - 25.21789 - -
2 C 0.88945 - - 6.29494 - -
3 C 0.88945 - - 6.29494 - -
4 C 0.98428 - - 6.27454 - -
5 C 0.98428 - - 6.27454 - -
6 C 0.89853 - - 6.11607 - -
7 C 0.89853 - - 6.11607 - -
8 C 1.03824 - - 6.23459 - -
9 H 0.02104 - - 0.77158 - -
10 H 0.02336 - - 0.79915 - -
11 H 0.02336 - - 0.79915 - -
12 H 0.02679 - - 0.77646 - -
13 H 0.02679 - - 0.77646 - -
14 H 0.01945 - - 0.77325 - -
15 H 0.01945 - - 0.77325 - -
16 Cl 0.21667 - - 17.49978 - -
----- --- ========= --------- --------- --------- --------- ---------
The numbers confirm that only Fe and the C atoms in the ring contribute significantly to the global delocalization. This atomic-resolution table is the only one printed by default (i.e., without --verbose option); it is the starting point for most discussions and is sufficient for many quick assessments.
4. Cyclic delocalization of π-electrons
Aromatic stabilization in molecular rings with π-conjugated bonds is associated with cyclic delocalization of π-electrons, whereas cross-ring orbital interactions are often locally destabilizing (cf. the antibonding interaction between para-related carbon atoms in benzene). To assess the effect of cyclic π-electron delocalization in the 8-membered ring (8MR) of the C7H7FeCl system, we perform an EDDBP calculation (the subscript P stands for pathway):
The --pathway (or -p) option specifies the bond connectivity along the ring. The --output-nobd option generates a new Gaussian FCHK file (c7h7fecl.NOBD.fchk) in which the Alpha MO coefficients are replaced by the spinless NOBDs. The --auto-pi (or -a) option automatically selects only π-type NOBDs based on their orbital composition; also, this flag triggers printing of the NOBD-resolution table by default, so one can immediately see which orbitals were classified as π-type.
The --auto-pi option correctly identifies the first four NOBDs as π-orbitals (marked with *); together they account for ~8.0e π-electrons of which ~5.4e is cyclically delocalized. To visually inspect the NOBDs and verify the automatic π-NOBD selection, open the generated c7h7fecl.NOBD.fchk file in Avogadro2 and select particular NOBD in Molecular Orbitals window (in generated c7h7fecl.NOBD.fchk file α-MOs ordered by increasing energy are replaced by NOBDs ordered by decreasing occupation number):
As follows from the summary of EDDB calculations, the average per-atom population of cyclically delocalized π-electrons in the 8MR is 0.67954e. However, the EDDBP atomic populations reveal that the distribution of electrons is not uniform: a noticeable π-deficiency on atoms 2 and 3 is a fingerprint of high-energy ionic resonance structures contributing to the wavefunction and reducing the aromatic stabilization energy. In this sense, the capacity of the π-channel to connect all eight atomic orbitals in a coherent cycle is limited by those two carbon atoms.
The Cyclic Delocalization Index (CDI), printed at the end of the summary, equals the smallest per-atom EDDBP value — the maximum population of π-electrons (per atom) that can be uniformly delocalized around the cycle. Here, the cyclic delocalization index is about ~0.54; for comparison, the CDI for benzene is ~0.89, so the relative π-aromaticity in the 8MR of C7H7FeCl can be estimated as ~61%.
5. 3d-orbital conjugation topology
The remaining open question is: what does the 3d2 orbital on Fe look like? To find out, simply run:
and open the generated c7h7fecl.NMB.fchk file in Avogadro2 (select Analyze → Create Surfaces). According to the resulting NVB-resolution table,
output
> Printing results of electron population analysis in orbital resolution...
+------------+-------------+---------------------------+-----------------------------+
| Atom | Orbital | Electron delocalization | Electron population |
+------------+-------------+---------------------------+-----------------------------+
Index Sym Index NMB Total Alpha Beta Total Alpha Beta
----- --- ----- ---- ======= ------- ------- ------- ------- -------
1 Fe 1 1s 0.00000 - - 0.00000 - -
1 Fe 2 2s 0.00000 - - 0.00000 - -
1 Fe 3 3s 0.00000 - - 0.00000 - -
1 Fe 4 4s 0.00827 - - 0.01135 - -
1 Fe 5 2p1 0.00000 - - 0.00000 - -
1 Fe 6 3p1 0.00000 - - 0.00000 - -
1 Fe 7 2p2 0.00000 - - 0.00000 - -
1 Fe 8 3p2 0.00000 - - 0.00000 - -
1 Fe 9 2p3 0.00000 - - 0.00000 - -
1 Fe 10 3p3 0.00000 - - 0.00000 - -
1 Fe 11 3d1 0.01020 - - 0.01400 - -
1 Fe 12 3d2 0.64008 - - 1.39783 - -
1 Fe 13 3d3 0.03492 - - 0.07622 - -
1 Fe 14 3d4 0.00192 - - 0.00264 - -
1 Fe 15 3d5 0.00275 - - 0.00377 - -
2 C 16 1s 0.00000 - - 0.00000 - -
2 C 17 2s 0.00029 - - 0.00052 - -
2 C 18 2p1 0.54147 - - 0.81764 - -
2 C 19 2p2 0.00001 - - 0.00001 - -
2 C 20 2p3 0.00004 - - 0.00007 - -
3 C 21 1s 0.00000 - - 0.00000 - -
3 C 22 2s 0.00029 - - 0.00052 - -
3 C 23 2p1 0.54147 - - 0.81764 - -
3 C 24 2p2 0.00001 - - 0.00001 - -
3 C 25 2p3 0.00004 - - 0.00007 - -
4 C 26 1s 0.00000 - - 0.00000 - -
4 C 27 2s 0.00001 - - 0.00001 - -
4 C 28 2p1 0.75021 - - 1.01068 - -
4 C 29 2p2 0.00007 - - 0.00010 - -
4 C 30 2p3 0.00002 - - 0.00003 - -
5 C 31 1s 0.00000 - - 0.00000 - -
5 C 32 2s 0.00001 - - 0.00001 - -
5 C 33 2p1 0.75021 - - 1.01068 - -
5 C 34 2p2 0.00007 - - 0.00010 - -
5 C 35 2p3 0.00002 - - 0.00003 - -
6 C 36 1s 0.00000 - - 0.00000 - -
6 C 37 2s 0.00000 - - 0.00000 - -
6 C 38 2p1 0.67720 - - 0.89585 - -
6 C 39 2p2 0.00003 - - 0.00004 - -
6 C 40 2p3 0.00008 - - 0.00011 - -
7 C 41 1s 0.00000 - - 0.00000 - -
7 C 42 2s 0.00000 - - 0.00000 - -
7 C 43 2p1 0.67720 - - 0.89585 - -
7 C 44 2p2 0.00003 - - 0.00004 - -
7 C 45 2p3 0.00008 - - 0.00011 - -
8 C 46 1s 0.00000 - - 0.00000 - -
8 C 47 2s 0.00000 - - 0.00000 - -
8 C 48 2p1 0.79919 - - 1.01793 - -
8 C 49 2p2 0.00010 - - 0.00013 - -
8 C 50 2p3 0.00000 - - 0.00000 - -
----- --- ----- ---- ======= ------- ------- ------- ------- -------
the 3d2 (Fe) orbital has index 12 in the NMB layout (the first 11 NMB orbitals on Fe are core 1s/2s/3s/4s/2p/3p, then come the five 3d). A visual inspection of Molecular Orbital no. 12 reveals this is a 3dXY orbital, which means the metallacycle is a Möbius-Craig aromatic with a predominantly δ-type conjugation topology (see the reference paper for more examples of metallacycles with π- and δ-conjugation topology).
6. Final conclusion
The EDDB analysis at the CAM-B3LYP/def2-SVP theory level reveals that the model C7H7FeCl metallacycle (in its ground-state configuration) can be classified as an 8π-electron Möbius-Craig aromatic species with predominating δ-conjugation topology involving 3dXY (Fe) orbital. The effectiveness of cyclic delocalization of π-bonds in this 8-membered ring is relatively high - about 61% of the one in archetypical π-aromatic benzene.
