Real-Space Second Chern Number

• Real-space second Chern number is a 4D topological invariant defined through position-space operators and spectral projectors, quantizing the quantum Hall response.
• It replaces momentum derivatives with spatial commutators and employs finite-size estimators and kernel polynomial methods for robust computations in disordered systems.
• The framework generalizes lower-dimensional Chern markers and underpins quantized responses in systems including photonic fibers and hyperbolic lattices.

Searching arXiv for papers on real-space second Chern number and closely related formulations. The real-space second Chern number is a higher-dimensional topological invariant formulated directly in terms of position-space operators and the occupied-state projector, rather than exclusively through Berry curvature on a translationally invariant Brillouin torus. In four-dimensional gapped systems, it is the natural counterpart of the first Chern number in two dimensions, and it quantizes the topology underlying the 4D quantum Hall response. Recent work has established exact real-space formulas, finite-system estimators, and kernel-polynomial implementations for disordered 4D Chern insulators, while related constructions have shown that a second Chern number can also emerge in mixed momentum–real-space parameter spaces and in hyperbolic lattices with an effective 4D Brillouin zone (Shiina et al., 20 Feb 2025, Chen et al., 25 Jul 2025, Lu et al., 2016, Zhang et al., 2023).

1. Definition and topological setting

In momentum space, the second Chern number for occupied bands over a four-dimensional Brillouin zone is written as

$$C_2 = \frac{1}{32\pi^2} \int_{\mathrm{BZ}} d^4k \,\epsilon^{ijkl}\, \mathrm{Tr}\!\left(F_{ij}(k)\,F_{kl}(k)\right),$$

which is the standard 4D Berry-curvature expression in translationally invariant systems (Shiina et al., 20 Feb 2025). In wedge-product form, one also encounters

$$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$

with $\Omega=dA-iA\wedge A$ and $$A_\mu^{mn}(k)=i\langle u_m(k)|\partial_{k_\mu}u_n(k)\rangle$$ for the occupied-band subbundle (Zhang et al., 2023). These normalizations are paper-specific conventions, but both encode the same integer invariant.

The real-space program replaces momentum derivatives by position-space derivations and expresses the invariant using the spectral projector $P$ onto occupied states. The exact marker form derived for translation-invariant systems is

$$C_2 = \frac{1}{2!}\left(\frac{2\pi}{i}\right)^2 \epsilon^{\mu\nu\rho\sigma} \, \mathrm{tr}\!\left[P\,[X_\mu,P]\,[X_\nu,P]\,[X_\rho,P]\,[X_\sigma,P]\right]_{ii},$$

where the trace is taken over internal degrees of freedom at a fixed site $i$ (Shiina et al., 20 Feb 2025). For translationally invariant systems, the expression is independent of $i$; for disordered systems, the appropriate replacement is a trace per unit volume or a windowed finite-size estimator (Shiina et al., 20 Feb 2025).

A closely related projected-position form used for large-scale kernel-polynomial computations is

$$C_2 = - 2 \pi^2 \epsilon^{j_1 j_2 j_3 j_4} \mathrm{Tr}\!\left[ P X_{j_1} P X_{j_2} P X_{j_3} P X_{j_4} P \right],$$

evaluated as a bulk trace per volume in open systems (Chen et al., 25 Jul 2025). The same work notes that the commutator form

$$C_2 = - 2 \pi^2 \epsilon^{abcd} \mathrm{Tr}\!\left[ P [X_a, P] [X_b, P] [X_c, P] [X_d, P] \right]$$

is equivalent for gapped systems with a well-defined trace per volume, although the calculations are performed with the projected-position product form (Chen et al., 25 Jul 2025).

This framework places the real-space second Chern number within noncommutative geometry. Defining $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$0, one has

$$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$1

which matches the marker form through $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$2 (Shiina et al., 20 Feb 2025).

2. Exact real-space formulas and finite-system estimators

An important development is the extension of Kitaev’s real-space formulation of the first Chern number to the second Chern number (Shiina et al., 20 Feb 2025). In addition to the site-local marker, that work derives a projector-partition form using half-space lattice projectors $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$3:

$$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$4

and equivalently

$$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$5

This is presented as the recommended finite-system estimator in practice (Shiina et al., 20 Feb 2025).

The same paper also proves a five-region deformation-invariant partition formula. Partitioning the 4D lattice into five disjoint regions $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$6 with projectors $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$7, one obtains

$$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$8

This form is described as strictly topological in the infinite system, although on finite systems with truncation it can be slightly less accurate than the half-space projector formula (Shiina et al., 20 Feb 2025).

For disordered systems, quantization requires a spectral gap at the Fermi energy $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$9, or more generally a mobility gap, together with locality of the Fermi projector $\Omega=dA-iA\wedge A$0 and decay of the commutators $\Omega=dA-iA\wedge A$1 (Shiina et al., 20 Feb 2025). Under these conditions, the real-space formulas yield an integer stable under weak disorder and smooth perturbations (Shiina et al., 20 Feb 2025). A practical estimator uses a window projector $\Omega=dA-iA\wedge A$2 selecting a hypercube of linear size $\Omega=dA-iA\wedge A$3:

$\Omega=dA-iA\wedge A$4

with $\Omega=dA-iA\wedge A$5 (Shiina et al., 20 Feb 2025). In numerics, choosing $\Omega=dA-iA\wedge A$6 for a system of linear size $\Omega=dA-iA\wedge A$7 is reported to effectively cancel boundary terms (Shiina et al., 20 Feb 2025).

A recurrent point in this literature is the necessity of bulk windowing. Taking the full finite-system trace can produce cancellations between artificial boundaries and drive the estimator to zero, whereas intermediate windows suppress boundary artifacts (Shiina et al., 20 Feb 2025). This suggests that the real-space second Chern number is best regarded not as a naive finite-volume trace, but as a bulk trace-per-volume quantity regularized to exclude nonuniversal edge contributions.

3. Computational frameworks

Two complementary computational strategies appear in the cited work. The first is a direct finite-size evaluation based on the half-space projector formula or the equivalent commutator marker (Shiina et al., 20 Feb 2025). The second is a scalable kernel polynomial method (KPM) implementation that avoids diagonalization and allows system sizes far beyond exact projector construction (Chen et al., 25 Jul 2025).

For the projector-based approach, the central object is the Fermi projector $\Omega=dA-iA\wedge A$8. One route is spectral flattening: define $\Omega=dA-iA\wedge A$9 and $$A_\mu^{mn}(k)=i\langle u_m(k)|\partial_{k_\mu}u_n(k)\rangle$$0 when the spectrum is gapped around $$A_\mu^{mn}(k)=i\langle u_m(k)|\partial_{k_\mu}u_n(k)\rangle$$1 (Shiina et al., 20 Feb 2025). The same work also describes KPM/Chebyshev approximations to the Heaviside step function and rational approximations via shifted inverses $$A_\mu^{mn}(k)=i\langle u_m(k)|\partial_{k_\mu}u_n(k)\rangle$$2 (Shiina et al., 20 Feb 2025). The paper emphasizes that for disordered models it is often more efficient to access $$A_\mu^{mn}(k)=i\langle u_m(k)|\partial_{k_\mu}u_n(k)\rangle$$3 implicitly as a linear map on vectors, avoiding storage of dense matrices (Shiina et al., 20 Feb 2025).

The KPM study makes this implicit strategy explicit. After rescaling the Hamiltonian to $$A_\mu^{mn}(k)=i\langle u_m(k)|\partial_{k_\mu}u_n(k)\rangle$$4 with spectrum in $$A_\mu^{mn}(k)=i\langle u_m(k)|\partial_{k_\mu}u_n(k)\rangle$$5,

$$A_\mu^{mn}(k)=i\langle u_m(k)|\partial_{k_\mu}u_n(k)\rangle$$6

the projector is approximated by a truncated Chebyshev series with Jackson kernel,

$$A_\mu^{mn}(k)=i\langle u_m(k)|\partial_{k_\mu}u_n(k)\rangle$$7

with moments

$$A_\mu^{mn}(k)=i\langle u_m(k)|\partial_{k_\mu}u_n(k)\rangle$$8

and Jackson coefficients

$$A_\mu^{mn}(k)=i\langle u_m(k)|\partial_{k_\mu}u_n(k)\rangle$$9

The recursive action on a vector uses

$P$0

(Chen et al., 25 Jul 2025).

The trace per volume is then estimated stochastically on a bulk subregion $P$1:

$P$2

where $P$3 are random phase vectors supported in $P$4 (Chen et al., 25 Jul 2025). For the reported 4D calculations, the paper uses $P$5 and $P$6 (Chen et al., 25 Jul 2025). Each $P$7 evaluation proceeds by repeated applications of $P$8 and coordinate operators $P$9, summed over the $$C_2 = \frac{1}{2!}\left(\frac{2\pi}{i}\right)^2 \epsilon^{\mu\nu\rho\sigma} \, \mathrm{tr}\!\left[P\,[X_\mu,P]\,[X_\nu,P]\,[X_\rho,P]\,[X_\sigma,P]\right]_{ii},$$0 nonzero Levi-Civita terms (Chen et al., 25 Jul 2025).

The two approaches differ mainly in scale and implementation details. The exact projector formulas establish the invariant and its finite-size regularizations (Shiina et al., 20 Feb 2025). The KPM formulation supplies the numerical route for large sparse Hamiltonians, including systems with $$C_2 = \frac{1}{2!}\left(\frac{2\pi}{i}\right)^2 \epsilon^{\mu\nu\rho\sigma} \, \mathrm{tr}\!\left[P\,[X_\mu,P]\,[X_\nu,P]\,[X_\rho,P]\,[X_\sigma,P]\right]_{ii},$$1 sites, and extends the same logic to an exploratory six-dimensional calculation of the third Chern number (Chen et al., 25 Jul 2025).

4. Wilson–Dirac benchmark and disorder

The canonical test bed for real-space evaluations of $$C_2 = \frac{1}{2!}\left(\frac{2\pi}{i}\right)^2 \epsilon^{\mu\nu\rho\sigma} \, \mathrm{tr}\!\left[P\,[X_\mu,P]\,[X_\nu,P]\,[X_\rho,P]\,[X_\sigma,P]\right]_{ii},$$2 is the four-dimensional Wilson–Dirac model. One formulation used in the projector-based study is

$$C_2 = \frac{1}{2!}\left(\frac{2\pi}{i}\right)^2 \epsilon^{\mu\nu\rho\sigma} \, \mathrm{tr}\!\left[P\,[X_\mu,P]\,[X_\nu,P]\,[X_\rho,P]\,[X_\sigma,P]\right]_{ii},$$3

with periodic boundary conditions on sites $$C_2 = \frac{1}{2!}\left(\frac{2\pi}{i}\right)^2 \epsilon^{\mu\nu\rho\sigma} \, \mathrm{tr}\!\left[P\,[X_\mu,P]\,[X_\nu,P]\,[X_\rho,P]\,[X_\sigma,P]\right]_{ii},$$4 and truncation window $$C_2 = \frac{1}{2!}\left(\frac{2\pi}{i}\right)^2 \epsilon^{\mu\nu\rho\sigma} \, \mathrm{tr}\!\left[P\,[X_\mu,P]\,[X_\nu,P]\,[X_\rho,P]\,[X_\sigma,P]\right]_{ii},$$5 selecting $$C_2 = \frac{1}{2!}\left(\frac{2\pi}{i}\right)^2 \epsilon^{\mu\nu\rho\sigma} \, \mathrm{tr}\!\left[P\,[X_\mu,P]\,[X_\nu,P]\,[X_\rho,P]\,[X_\sigma,P]\right]_{ii},$$6 (Shiina et al., 20 Feb 2025). The paper reports results for $$C_2 = \frac{1}{2!}\left(\frac{2\pi}{i}\right)^2 \epsilon^{\mu\nu\rho\sigma} \, \mathrm{tr}\!\left[P\,[X_\mu,P]\,[X_\nu,P]\,[X_\rho,P]\,[X_\sigma,P]\right]_{ii},$$7, $$C_2 = \frac{1}{2!}\left(\frac{2\pi}{i}\right)^2 \epsilon^{\mu\nu\rho\sigma} \, \mathrm{tr}\!\left[P\,[X_\mu,P]\,[X_\nu,P]\,[X_\rho,P]\,[X_\sigma,P]\right]_{ii},$$8, and $$C_2 = \frac{1}{2!}\left(\frac{2\pi}{i}\right)^2 \epsilon^{\mu\nu\rho\sigma} \, \mathrm{tr}\!\left[P\,[X_\mu,P]\,[X_\nu,P]\,[X_\rho,P]\,[X_\sigma,P]\right]_{ii},$$9 (Shiina et al., 20 Feb 2025). In the clean phase diagram, the real-space estimator reproduces quantized plateaus as $i$0 varies; for example, at $i$1 with $i$2, the exact value is $i$3, while the finite-size estimate gives $i$4 for $i$5 and approaches the integer value as the system size increases (Shiina et al., 20 Feb 2025).

The KPM study uses the $i$6-dimensional Wilson–Dirac family

$i$7

with

$i$8

on a hypercubic lattice of linear size $i$9, volume $i$0, fixing $i$1 and $i$2 in all numerics (Chen et al., 25 Jul 2025). For the 4D case, the explicit gamma matrices are

$i$3

(Chen et al., 25 Jul 2025). The paper computes $i$4 versus mass for $i$5 and reports rapid convergence to the momentum-space phase diagram. At $i$6, where one expects $i$7, the deviation $i$8 decays exponentially with $i$9 (Chen et al., 25 Jul 2025).

Both studies treat disorder directly in real space. In the KPM work, Anderson-type on-site disorder is added as

$$C_2 = - 2 \pi^2 \epsilon^{j_1 j_2 j_3 j_4} \mathrm{Tr}\!\left[ P X_{j_1} P X_{j_2} P X_{j_3} P X_{j_4} P \right],$$0

with independent $$C_2 = - 2 \pi^2 \epsilon^{j_1 j_2 j_3 j_4} \mathrm{Tr}\!\left[ P X_{j_1} P X_{j_2} P X_{j_3} P X_{j_4} P \right],$$1 uniformly distributed in $$C_2 = - 2 \pi^2 \epsilon^{j_1 j_2 j_3 j_4} \mathrm{Tr}\!\left[ P X_{j_1} P X_{j_2} P X_{j_3} P X_{j_4} P \right],$$2 (Chen et al., 25 Jul 2025). The disorder-averaged phase diagram $$C_2 = - 2 \pi^2 \epsilon^{j_1 j_2 j_3 j_4} \mathrm{Tr}\!\left[ P X_{j_1} P X_{j_2} P X_{j_3} P X_{j_4} P \right],$$3 agrees well with boundaries predicted by the self-consistent Born approximation at weak-to-moderate disorder (Chen et al., 25 Jul 2025). The SCBA self-energy is written as

$$C_2 = - 2 \pi^2 \epsilon^{j_1 j_2 j_3 j_4} \mathrm{Tr}\!\left[ P X_{j_1} P X_{j_2} P X_{j_3} P X_{j_4} P \right],$$4

and the practical comparison uses the lowest-order Born approximation with $$C_2 = - 2 \pi^2 \epsilon^{j_1 j_2 j_3 j_4} \mathrm{Tr}\!\left[ P X_{j_1} P X_{j_2} P X_{j_3} P X_{j_4} P \right],$$5 on the right-hand side to determine where the renormalized spectrum closes its gap (Chen et al., 25 Jul 2025).

The projector-based study distinguishes two disorder mechanisms in the 4D Wilson–Dirac model: random hopping and random mass (Shiina et al., 20 Feb 2025). With $$C_2 = - 2 \pi^2 \epsilon^{j_1 j_2 j_3 j_4} \mathrm{Tr}\!\left[ P X_{j_1} P X_{j_2} P X_{j_3} P X_{j_4} P \right],$$6, $$C_2 = - 2 \pi^2 \epsilon^{j_1 j_2 j_3 j_4} \mathrm{Tr}\!\left[ P X_{j_1} P X_{j_2} P X_{j_3} P X_{j_4} P \right],$$7, and averages over 10 disorder realizations on the $$C_2 = - 2 \pi^2 \epsilon^{j_1 j_2 j_3 j_4} \mathrm{Tr}\!\left[ P X_{j_1} P X_{j_2} P X_{j_3} P X_{j_4} P \right],$$8 system, random hopping drives a direct transition from the $$C_2 = - 2 \pi^2 \epsilon^{j_1 j_2 j_3 j_4} \mathrm{Tr}\!\left[ P X_{j_1} P X_{j_2} P X_{j_3} P X_{j_4} P \right],$$9 phase to a trivial phase as the half-filled gap closes and reopens, while random mass produces a broad nearly gapless regime in which $$C_2 = - 2 \pi^2 \epsilon^{abcd} \mathrm{Tr}\!\left[ P [X_a, P] [X_b, P] [X_c, P] [X_d, P] \right]$$0 is ill-defined rather than a trivial gapped phase (Shiina et al., 20 Feb 2025). These results place the real-space second Chern number in the same operational category as lower-dimensional Chern markers: it remains informative in the presence of disorder provided that the occupied projector remains local and the relevant gap structure persists.

5. Mixed-space and non-Euclidean realizations

Although the phrase “real-space second Chern number” often refers to projector formulas in physical space, related work shows that the same invariant can arise when one of the four coordinates is a real-space geometric parameter. In a magnetized double-gyroid photonic crystal, a helical modulation

$$C_2 = - 2 \pi^2 \epsilon^{abcd} \mathrm{Tr}\!\left[ P [X_a, P] [X_b, P] [X_c, P] [X_d, P] \right]$$1

produces localized one-way modes along the helix axis (Lu et al., 2016). Far from the helix core, the Bloch Hamiltonian depends smoothly on $$C_2 = - 2 \pi^2 \epsilon^{abcd} \mathrm{Tr}\!\left[ P [X_a, P] [X_b, P] [X_c, P] [X_d, P] \right]$$2, so the eigenstates are periodic in all four variables and define a 4D torus $$C_2 = - 2 \pi^2 \epsilon^{abcd} \mathrm{Tr}\!\left[ P [X_a, P] [X_b, P] [X_c, P] [X_d, P] \right]$$3 (Lu et al., 2016). The second Chern number is then defined on this mixed momentum–real-space parameter space as

$$C_2 = - 2 \pi^2 \epsilon^{abcd} \mathrm{Tr}\!\left[ P [X_a, P] [X_b, P] [X_c, P] [X_d, P] \right]$$4

using the paper’s convention and normalization (Lu et al., 2016).

In that construction, time-reversal symmetry must be broken for $$C_2 = - 2 \pi^2 \epsilon^{abcd} \mathrm{Tr}\!\left[ P [X_a, P] [X_b, P] [X_c, P] [X_d, P] \right]$$5 because $$C_2 = - 2 \pi^2 \epsilon^{abcd} \mathrm{Tr}\!\left[ P [X_a, P] [X_b, P] [X_c, P] [X_d, P] \right]$$6 with $$C_2 = - 2 \pi^2 \epsilon^{abcd} \mathrm{Tr}\!\left[ P [X_a, P] [X_b, P] [X_c, P] [X_d, P] \right]$$7 is odd under time reversal while $$C_2 = - 2 \pi^2 \epsilon^{abcd} \mathrm{Tr}\!\left[ P [X_a, P] [X_b, P] [X_c, P] [X_d, P] \right]$$8 is even (Lu et al., 2016). The simplest helical modulation yields $$C_2 = - 2 \pi^2 \epsilon^{abcd} \mathrm{Tr}\!\left[ P [X_a, P] [X_b, P] [X_c, P] [X_d, P] \right]$$9 for $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$00, and the observed number and chirality of one-way fiber modes satisfy $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$01 (Lu et al., 2016). This is not a projector-based real-space marker, but it is a real-space manifestation of the second Chern number in which the extra coordinate is the physical helix angle rather than a synthetic momentum dimension (Lu et al., 2016).

A different route appears in hyperbolic band topology. In the $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$02 hyperbolic tiling, four independent translation directions generated by the Fuchsian group $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$03 are equipped with $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$04 twist phases, producing a 4D Brillouin zone $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$05 with angles $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$06 (Zhang et al., 2023). The corresponding 4D Bloch Hamiltonian is

$$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$07

with

$$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$08

and the second Chern number again given by

$$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$09

(Zhang et al., 2023).

For $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$10, $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$11, the central gap has $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$12 while the first Chern numbers vanish; for $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$13, $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$14, the central 4D gap closes and the side gaps have $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$15 but nontrivial first Chern numbers on certain 2D slices, such as $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$16 near $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$17 (Zhang et al., 2023). The authors do not implement a projector-based real-space marker; instead they construct abelian periodic-boundary-condition clusters whose discrete translation operators can be simultaneously diagonalized, enabling a discrete momentum-space evaluation of $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$18 even in finite hyperbolic lattices (Zhang et al., 2023). This suggests a broader taxonomy: projector markers, mixed-space invariants, and twist-space invariants are distinct constructions that all realize second-Chern topology outside the standard 4D crystal setting.

6. Physical interpretation, responses, and limitations

The second Chern number is the topological invariant of the 4D quantum Hall effect and controls a quantized nonlinear response coefficient in four dimensions (Chen et al., 25 Jul 2025). Dimensional reduction connects 4D Chern insulators to 3D time-reversal invariant topological insulators, where pumping along the extra dimension yields the magnetoelectric response $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$19, and $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$20 controls the topological term proportional to $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$21 in the reduced description (Shiina et al., 20 Feb 2025). In this sense, the real-space second Chern number is not merely a numerical substitute for Berry-curvature integration; it is the disorder-compatible form of the bulk invariant governing higher-dimensional response.

The photonic mixed-space realization gives an explicit transport interpretation. A nonzero $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$22 produces one-way fiber modes localized on the helix axis, with immunity to Rayleigh and Mie scattering, strong inhibition of stimulated Brillouin scattering and backward Raman scattering, elimination of Fresnel reflection, and single-polarization one-way propagation (Lu et al., 2016). In that system, the coupling between momentum-space Berry curvature and the $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$23-winding of the Dirac mass yields chiral one-dimensional defect transport (Lu et al., 2016).

The hyperbolic circuit realization provides a different boundary signature. Midgap boundary states appear only in bandgaps with nonzero Chern invariants; for the central gap with $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$24 at $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$25, $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$26, a significant boundary-node impedance peak occurs in the frequency range $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$27 MHz, and the spatial impedance map at $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$28 MHz is localized around partially open boundary units (Zhang et al., 2023). Because the mismatch between the 4D momentum torus and 2D real-space hyperbolic geometry complicates fixed-direction chiral transport, the observed correspondence is formulated in terms of bulk–boundary spectral signatures rather than a direct transport coefficient (Zhang et al., 2023).

Several limitations are common across formulations. Exact quantization in real space requires a spectral gap or mobility gap and locality of $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$29; near gap closures, finite-size estimators deviate from integers and become sensitive to boundary conditions (Shiina et al., 20 Feb 2025). In KPM, kernel broadening is of order $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$30, so narrow gaps require larger $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$31; otherwise the projector is smeared and $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$32 can lose quantization (Chen et al., 25 Jul 2025). Higher dimensions amplify the computational burden: in 6D, the third-Chern calculation involves $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$33 nonzero Levi-Civita terms and much larger Hilbert spaces, so the reported results remain qualitative rather than fully quantized (Chen et al., 25 Jul 2025).

A common misconception is that any “real-space” manifestation of $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$34 must be a local projector marker. The literature considered here shows that this is too narrow. The exact real-space formulas of the Kitaev type provide one rigorous route (Shiina et al., 20 Feb 2025), but mixed momentum–real-space constructions in photonics (Lu et al., 2016) and twist-generated 4D parameter spaces in hyperbolic lattices (Zhang et al., 2023) also realize second-Chern topology in systems lacking an ordinary 4D crystal momentum. A plausible implication is that “real-space second Chern number” now names a family of related higher-dimensional invariants whose common feature is independence from a conventional 4D translationally invariant Bloch setting.

7. Outlook and generalization

The recent literature establishes a computational and conceptual progression. First, exact real-space formulas now place the second Chern number on the same footing as established lower-dimensional Chern markers, with explicit equivalence to the momentum-space invariant and practical finite-system estimators (Shiina et al., 20 Feb 2025). Second, large-scale KPM calculations demonstrate that these formulas are usable at the scale of four-dimensional lattices with $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$35 sites and can capture disorder-driven topology in agreement with self-consistent Born approximation predictions (Chen et al., 25 Jul 2025). Third, the invariant has been realized in nonstandard geometries and parameter spaces, including mixed-space photonic fibers and hyperbolic circuit networks (Lu et al., 2016, Zhang et al., 2023).

The generalization to higher Chern classes is already explicit. The KPM work extends the projected-position formula to the third Chern number in six dimensions,

$$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$36

and the exact real-space framework states the broader index-theoretic structure for $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$37 dimensions (Chen et al., 25 Jul 2025, Shiina et al., 20 Feb 2025). Although full quantization in 6D remains limited by finite-size effects, the reported qualitative agreement suggests that higher-dimensional real-space topological invariants are computationally accessible (Chen et al., 25 Jul 2025).

Open directions are identified most explicitly in the hyperbolic setting: developing a robust real-space local marker for $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$38 on hyperbolic lattices, establishing chiral response measurements linked to $$C_2 = \frac{1}{8\pi^2}\int_{\mathrm{BZ}_{4\mathrm{D}}} \mathrm{d}^4k\,\mathrm{tr}\bigl(\Omega\wedge\Omega\bigr),$$39 in hyperbolic geometries, scaling to larger clusters while mitigating losses and parasitics, exploring disorder and interactions, and extending to other tessellations and symmetries (Zhang et al., 2023). In Euclidean lattices, the main frontier is algorithmic: combining KPM with tensor-network techniques to push accessible Hilbert spaces further, especially for higher Chern numbers (Chen et al., 25 Jul 2025). Taken together, these results indicate that the second Chern number has moved from a largely momentum-space invariant of ideal 4D band theory to a real-space computable quantity relevant to disorder, finite samples, curved lattices, and synthetic higher-dimensional platforms.