◼ 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. Q-Chem input files
To investigate bond delocalization effects in the C7H7FeCl metallacycle at the CAM-B3LYP/def2-SVP theory level, use the following input file:
$molecule
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
$end
$rem
JOBTYP SP
METHOD CAM-B3LYP
BASIS def2-SVP
GUI 2
$end
Run the calculation using Q-Chem (tested version: 6.4.0-1) to get the formatted checkpoint file c7h7fecl.fchk.
terminal
$ qchem -nt 4 c7h7fecl.in >& c7h7fecl.out
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.fchk --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: 49.79370 per atom: 3.11211
Total EDLB_G population: 41.50637 per atom: 2.59415
Total EDDB_G population: 8.28733 per atom: 0.51796
> 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: 122x71x116 (~100^3)
Progress: ######################################## 100%
Integrated density: 8.28793 (100.0% of 8.28733)
Time: 3s
The integrated density (~8.288e) almost perfectly matches the analytical EDDBG population (~8.287e) — this confirms the cube accurately represents the density.
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.fchk --verbose
First, the program reports what was parsed from the wavefunction file:
output
> Loading and parsing input file... 0s
Source file: c7h7fecl.fchk
File format: Gaussian Formatted Checkpoint File (FCHK)
Number of atoms: 16
Number of basis functions: 182 (spherical)
Number of contracted shells: 82
Number of primitives: 166
Number of molecular orbitals: 182
Effective core potential: No
Charge / multiplicity: 0 / 1
Total electrons: 92
Unrestricted: No
Reference type: HF / DFT
Electron population (Tr[D*S]): 92.00000
Next, it summarises the internal AO→NAO localization procedure:
output
> Building natural atomic orbitals (NAO)... 0s
Total number of NAOs: 182 = 21 (Cor) + 45 (Val) + 116 (Ryd)
Total electron population: 92.00000
- core population: 41.99889 (~45.7 %)
- valence population: 49.79370 (~54.1 %)
- rydberg population: 0.20741 (~0.2 %)
Target NAO representation: Natural valence basis (NVB = Val)
- NVB population: 49.79370 (~54.1 %)
Out of the 182 NAOs constructed by the internal driver, 21 are classified as core (Cor), 45 as valence (Val), and 116 as Rydberg (Ryd). The default NAO representation (NVB) keeps only the 45 valence NAOs that capture ~49.8e — the analysis below operates on this subspace.
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):
output
> Printing results of electron population analysis in orbital resolution...
+------------+-------------+---------------------------+-----------------------------+
| Atom | Orbital | Electron delocalization | Electron population |
+------------+-------------+---------------------------+-----------------------------+
Index Sym Index NVB Total Alpha Beta Total Alpha Beta
----- --- ----- ---- ======= ------- ------- ------- ------- -------
1 Fe 1 4s 0.04506 - - 0.41546 - -
1 Fe 2 3d1 0.11904 - - 1.61962 - -
1 Fe 3 3d2 0.23922 - - 0.99373 - -
1 Fe 4 3d3 0.80025 - - 1.40782 - -
1 Fe 5 3d4 0.06577 - - 1.93022 - -
1 Fe 6 3d5 0.07334 - - 0.85220 - -
2 C 7 2s 0.06723 - - 1.15492 - -
2 C 8 2p1 0.05403 - - 1.09689 - -
2 C 9 2p2 0.62691 - - 0.84263 - -
2 C 10 2p3 0.13906 - - 1.20120 - -
3 C 11 2s 0.06723 - - 1.15492 - -
3 C 12 2p1 0.05403 - - 1.09689 - -
3 C 13 2p2 0.62691 - - 0.84263 - -
3 C 14 2p3 0.13906 - - 1.20120 - -
4 C 15 2s 0.03171 - - 0.96152 - -
4 C 16 2p1 0.05678 - - 1.10321 - -
4 C 17 2p2 0.84733 - - 1.01893 - -
4 C 18 2p3 0.04574 - - 1.19031 - -
5 C 19 2s 0.03171 - - 0.96152 - -
5 C 20 2p1 0.05678 - - 1.10321 - -
5 C 21 2p2 0.84733 - - 1.01893 - -
5 C 22 2p3 0.04574 - - 1.19031 - -
6 C 23 2s 0.03189 - - 0.94738 - -
6 C 24 2p1 0.05783 - - 1.13633 - -
6 C 25 2p2 0.77213 - - 0.90451 - -
6 C 26 2p3 0.03621 - - 1.12834 - -
7 C 27 2s 0.03189 - - 0.94738 - -
7 C 28 2p1 0.05783 - - 1.13633 - -
7 C 29 2p2 0.77213 - - 0.90451 - -
7 C 30 2p3 0.03621 - - 1.12834 - -
8 C 31 2s 0.03028 - - 0.93835 - -
8 C 32 2p1 0.06255 - - 1.18218 - -
8 C 33 2p2 0.90585 - - 1.02758 - -
8 C 34 2p3 0.03867 - - 1.08611 - -
9 H 35 1s 0.02102 - - 0.77156 - -
10 H 36 1s 0.02303 - - 0.79899 - -
11 H 37 1s 0.02303 - - 0.79899 - -
12 H 38 1s 0.02652 - - 0.77642 - -
13 H 39 1s 0.02652 - - 0.77642 - -
14 H 40 1s 0.01932 - - 0.77325 - -
15 H 41 1s 0.01932 - - 0.77325 - -
16 Cl 42 3s 0.01555 - - 1.95185 - -
16 Cl 43 3p1 0.06848 - - 1.70454 - -
16 Cl 44 3p2 0.03668 - - 1.90569 - -
16 Cl 45 3p3 0.09412 - - 1.93713 - -
----- --- ----- ---- ======= ------- ------- ------- ------- -------
The data reveal that the main delocalization channel involves the 3d3 orbital on Fe, occupied by ~1.4e and contributing ~0.8e (~60%) to delocalization, together with the 2p2 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.34267 - - 7.21904 - -
2 C 0.88724 - - 4.29563 - -
3 C 0.88724 - - 4.29563 - -
4 C 0.98156 - - 4.27397 - -
5 C 0.98156 - - 4.27397 - -
6 C 0.89805 - - 4.11658 - -
7 C 0.89805 - - 4.11658 - -
8 C 1.03735 - - 4.23421 - -
9 H 0.02102 - - 0.77156 - -
10 H 0.02303 - - 0.79899 - -
11 H 0.02303 - - 0.79899 - -
12 H 0.02652 - - 0.77642 - -
13 H 0.02652 - - 0.77642 - -
14 H 0.01932 - - 0.77325 - -
15 H 0.01932 - - 0.77325 - -
16 Cl 0.21483 - - 7.49921 - -
----- --- ========= --------- --------- --------- --------- ---------
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.68018e. 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 3d3 orbital on Fe look like? To find out, simply run:
and open the generated c7h7fecl.NVB.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 NVB Total Alpha Beta Total Alpha Beta
----- --- ----- ---- ======= ------- ------- ------- ------- -------
1 Fe 1 4s 0.00826 - - 0.01134 - -
1 Fe 2 3d1 0.00277 - - 0.00379 - -
1 Fe 3 3d2 0.09884 - - 0.21490 - -
1 Fe 4 3d3 0.57797 - - 1.25698 - -
1 Fe 5 3d4 0.00007 - - 0.00009 - -
1 Fe 6 3d5 0.01233 - - 0.01693 - -
2 C 7 2s 0.00030 - - 0.00054 - -
2 C 8 2p1 0.01445 - - 0.02168 - -
2 C 9 2p2 0.52823 - - 0.79709 - -
2 C 10 2p3 0.00004 - - 0.00008 - -
3 C 11 2s 0.00030 - - 0.00054 - -
3 C 12 2p1 0.01445 - - 0.02168 - -
3 C 13 2p2 0.52823 - - 0.79709 - -
3 C 14 2p3 0.00004 - - 0.00008 - -
4 C 15 2s 0.00001 - - 0.00001 - -
4 C 16 2p1 0.01723 - - 0.02322 - -
4 C 17 2p2 0.73294 - - 0.98696 - -
4 C 18 2p3 0.00002 - - 0.00003 - -
5 C 19 2s 0.00001 - - 0.00001 - -
5 C 20 2p1 0.01723 - - 0.02322 - -
5 C 21 2p2 0.73294 - - 0.98696 - -
5 C 22 2p3 0.00002 - - 0.00003 - -
6 C 23 2s 0.00000 - - 0.00000 - -
6 C 24 2p1 0.01877 - - 0.02485 - -
6 C 25 2p2 0.65902 - - 0.87155 - -
6 C 26 2p3 0.00008 - - 0.00011 - -
7 C 27 2s 0.00000 - - 0.00000 - -
7 C 28 2p1 0.01877 - - 0.02485 - -
7 C 29 2p2 0.65902 - - 0.87155 - -
7 C 30 2p3 0.00008 - - 0.00011 - -
8 C 31 2s 0.00000 - - 0.00000 - -
8 C 32 2p1 0.02362 - - 0.03010 - -
8 C 33 2p2 0.77542 - - 0.98751 - -
8 C 34 2p3 0.00000 - - 0.00000 - -
----- --- ----- ---- ======= ------- ------- ------- ------- -------
the 3d3 (Fe) orbital has index 4. A visual inspection of Molecular Orbital no. 4 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.
