quantum hall effect fractional laughlin wavefunction composite fermion cherny simons
# Fractional Quantum Hall Effect: Laughlin Wavefunctions, Composite Fermions, and Anyonic Statistics
## 1. Introduction to the Fractional Quantum Hall Effect
The fractional quantum hall effect (FQHE) represents one of the most profound phenomena in condensed matter physics, emerging when a two-dimensional electron gas is subjected to a strong perpendicular magnetic field at specific filling factors ν = p/q (where p and q are integers with q odd). Unlike the integer quantum Hall effect, which can be understood through single-particle Landau level quantization, the FQHE is a collective many-body phenomenon arising from strong Coulomb interactions between electrons. The ground state exhibits topological order—a qualitative difference from conventional symmetry-breaking phases—characterized by long-range entanglement, fractional excitation charges e* = e/q, and anyonic (non-Fermionic, non-Bosonic) exchange statistics.
The discovery of the FQHE at filling factor ν = 1/3 in 1983 (Tsui, Stormer, and Gossard) immediately prompted the seminal theoretical explanation by Laughlin (1983), who proposed an exact wavefunction for the ground state. This Laughlin wavefunction has proven remarkably robust, explaining not only the ν = 1/3 FQHE but also hierarchical fillings at ν = 2/5, 2/7, 3/5, etc. The theory of composite fermions—a mean-field picture interpreting the FQHE ground state as a filled Fermi liquid of electrons bound to an even number of flux quanta—provides a complementary perspective that has guided understanding of the complete FQHE phase diagram.
This article develops the theoretical foundations of the FQHE, from Landau level physics through the Laughlin wavefunction, composite fermion picture, Chern-Simons gauge field theory, and the emergence of anyonic excitations. A companion Python solver computes pair correlation functions and analyzes the many-body structure of model FQHE states.
## 2. Landau Levels and Quantum Hall Geometry
A free electron in 2D subjected to a perpendicular magnetic field $B = B\hat{z}$ has Hamiltonian:
$$H = \frac{1}{2m}(\mathbf{p} - e\mathbf{A})^2,$$
where $\mathbf{A}$ is the vector potential. Choosing the symmetric gauge $\mathbf{A} = \frac{B}{2}(-y, x, 0)$, the Hamiltonian becomes:
$$H = \frac{\omega_c}{2}(a^\dagger a + \frac{1}{2}) + \frac{\omega_c}{2}(b^\dagger b + \frac{1}{2}),$$
where $\omega_c = \frac{eB}{m}$ is the cyclotron frequency and $a, b$ are ladder operators for two decoupled harmonic oscillators. The eigenvalues are:
$$E_{n_x, n_y} = \hbar\omega_c(n_x + n_y + 1),$$
independent of $n_x$ and $n_y$ individually. All eigenstates with the same $N = n_x + n_y$ form a highly degenerate Landau level, each with degeneracy $\phi_0 / A = e B A / (2\pi\hbar c)$, where $\phi_0 = h/e$ is the flux quantum and $A$ is the sample area.
The filling factor ν quantifies the fraction of Landau level states occupied:
$$ u = \frac{N_e}{N_\phi},$$
where $N_e$ is the number of electrons and $N_\phi = \Phi / \phi_0$ is the number of flux quanta (Landau level degeneracy). For $
u < 1$, not all states in the lowest Landau level (LLL) are filled.
## 3. The Laughlin Wavefunction
For filling factor $
u = 1/m$ with m odd (typically $m = 3$), Laughlin proposed the exact wavefunction:
$$\Psi_m(z_1, z_2, \ldots, z_N) = \prod_{i<j} (z_i - z_j)^m \exp\left(-\frac{1}{4l_B^2} \sum_{i=1}^N |z_i|^2 ight),$$
where $z_i = x_i + iy_i$ is the complex coordinate of electron i in the plane perpendicular to the field, and $l_B = \sqrt{\hbar/(eB)}$ is the magnetic length. The prefactor $\prod_{i<j}(z_i - z_j)^m$ ensures:
1. Zero-range correlations: The wavefunction vanishes as $(z_i - z_j)^m$ when electrons approach, preventing pair overlap and implementing a strong short-range repulsion characteristic of Coulomb interactions in the high-field limit.
2. Fermionic exchange: For electrons (fermions), exchanging two particles $(z_i, z_j) o (z_j, z_i)$ sends $z_i - z_j o z_j - z_i = -(z_i - z_j)$, giving the wavefunction a factor $(-1)^m = -1$ for odd m, confirming fermionic antisymmetry.
3. Gaussian envelope: The exponential factor $\exp(-\sum_i |z_i|^2/(4l_B^2))$ localizes electrons within the sample, with characteristic length scale $l_B$.
## 4. Hilbert Space Structure and Plasma Analogy
The Laughlin wavefunction can be understood through a classical plasma analogy. The probability density associated with the many-body wavefunction is:
$$|\Psi_m|^2 = \left| \prod_{i<j} (z_i - z_j)^m ight|^2 \exp\left(-\frac{1}{2l_B^2} \sum_i |z_i|^2 ight).$$
This can be rewritten as:
$$|\Psi_m|^2 \propto \exp\left(m \sum_{i<j} \ln|z_i - z_j| - \frac{1}{2l_B^2} \sum_i |z_i|^2 ight) = \exp(-\beta E_{ ext{classical}}),$$
where the classical Coulomb energy is:
$$E_{ ext{classical}} = -m \sum_{i<j} \ln|z_i - z_j| + \frac{1}{2l_B^2} \sum_i |z_i|^2.$$
This is precisely the energy of a classical 2D one-component plasma (2D OCP)—a system of point charges in a uniform neutralizing background—at effective inverse temperature $\beta = 1/E_{ ext{class}}$. The Laughlin wavefunction is thus the Boltzmann distribution of a classical plasma, explaining why it yields a highly correlated liquid state with vanishing pair correlation at short range.
## 5. Composite Fermion Picture
An alternative, highly useful perspective is the composite fermion (CF) mean-field theory, developed by Zhang, Hansson, and Kivelson. In this picture, each electron binds to an even number 2p of flux quanta (typically p=1, binding 2 flux quanta), forming a composite fermion. The effective magnetic field experienced by the composite fermion is:
$$B_{ ext{eff}}^* = B - 2p \phi_0 n,$$
where $n = N_e / A$ is the electron density. At filling factor $
u = 1/(2p+1)$, the density of electrons is tuned so that $B_{ ext{eff}}^* = 0$, and the composite fermions form a filled (integer) Fermi liquid state with Fermi wavevector $k_F = \sqrt{2\pi n}$.
The composite fermion wavefunctions at general filling $
u = p / (2pq + 1)$ (where the numerator p indicates partial filling) are constructed by taking Slater determinants of the composite fermion orbitals in the effective magnetic field and attaching flux quanta. Remarkably, this composite fermion theory predicts the same ground state structure as the Laughlin wavefunction for $
u = 1/3$ and extends to hierarchical fractions.
## 6. Fractional Excitations and Anyonic Statistics
The key signature of the FQHE ground state is the existence of fractional excitations with fractional charge. For the $
u = 1/m$ state, a quasihole—created by removing one electron—carries charge $+e/m$, while a quasielectron—created by adding one electron—carries charge $-e/m$. These fractional charges have been verified experimentally through shot noise measurements.
More remarkably, these fractional excitations obey anyonic statistics—neither Fermionic nor Bosonic. When two quasiholes are exchanged adiabatically around each other in the 2D plane, the wavefunction acquires a phase:
$$\psi_{ ext{final}} = e^{i heta} \psi_{ ext{initial}},$$
where the exchange phase is:
$$ heta = \pi / m.$$
For m = 3 (the ν = 1/3 state), $ heta = \pi/3$. This is neither 0 (bosons) nor π (fermions), defining anyons. The existence of anyonic statistics in the FQHE is intimately connected to topological order: the exchange phase is not derived from local interactions but emerges as a topological property of the ground state.
## 7. Chern-Simons Gauge Theory and Flux Attachment
The composite fermion picture can be formalized using Chern-Simons gauge theory. The key idea is flux attachment: dynamically binding each electron to an even number of flux quanta via an emergent Chern-Simons gauge field $a_\mu$:
$$a_0 = \phi_{ ext{attached}},$$
where the attached flux is proportional to the electron density. The Chern-Simons term in the action:
$$S_{ ext{CS}} = \int dt d^2x \frac{k}{4\pi} \epsilon^{\mu u\lambda} a_\mu \partial_ u a_\lambda,$$
with Chern-Simons level k = 2p, encodes the attachment of 2p flux quanta per electron. After flux attachment, electrons see an effective magnetic field:
$$B_{ ext{eff}} = B - 2p n \phi_0,$$
converting the problem to a simpler integer filling of an effective lower Landau level (lowest Landau level for composite fermions).
## 8. Topological Order and Entanglement Entropy
The FQHE ground state exhibits topological order—a form of quantum order that cannot be characterized by local order parameters. Key signatures include:
1. Degenerate ground state manifold: On a torus with periodic boundary conditions, the ν = 1/m state has q-fold degeneracy (where q is the denominator of ν), corresponding to q distinct topological sectors.
2. Long-range entanglement: The ground state is entangled over arbitrarily large distances; the entanglement entropy across a partition scales as $S_A = \alpha L - \gamma$, where L is the boundary length, α is an area law coefficient, and γ is the topological entanglement entropy $\gamma = \ln D$, with D being the total quantum dimension of topological excitations.
3. Robust ground state: The topological nature makes the ground state robust against local perturbations—edge states and non-abelian anyons are protected topologically.
## 9. Edge States and Luttinger Liquid Physics
The boundary of a FQHE sample breaks translation symmetry, giving rise to chiral edge states. For the ν = 1/m state, there is a single edge mode with dispersion:
$$\omega(k) = v_e k,$$
where $v_e$ is the edge state velocity (typically ~10⁵ m/s). The edge mode is described by a Luttinger liquid with Tomonaga-Luttinger parameter:
$$K = \frac{1}{m},$$
characterizing the strength of electron-electron interactions. Tunneling into the edge exhibits power-law suppression at low bias, with exponent related to K. The edge current is quantized at $I = \frac{
u e^2}{h} V$ (fractional quantum Hall resistance).
## 10. Non-Abelian Anyons and Topological Quantum Computing
At certain filling factors (ν = 5/2 and others), the FQHE ground state supports non-abelian anyons—quasiparticles whose exchange operations do not commute, allowing them to implement non-trivial unitary operations on a protected quantum subspace. The most famous example is the Ising anyon (Moore-Read state at ν = 5/2), which satisfies:
$$R^2 = -1,$$
where R is the exchange matrix. Braiding sequences of Ising anyons can implement arbitrary SU(2) rotations on a computational subspace, forming a foundation for topological quantum computation.
## 11. Experimental Realizations
FQHE has been observed in:
1. 2DEGs in semiconductors: GaAs/AlGaAs heterostructures exhibit FQHE at numerous filling factors (ν = 1/3, 2/5, 3/7, 4/9, 2/3, etc.) due to high electron mobility (> 10⁶ cm²/V·s) enabling long coherence lengths.
2. Graphene: Single-layer graphene exhibits FQHE at even-denominator fillings (ν = 0, ±1, ±2, ±3) due to its unique linear dispersion and valley/spin degrees of freedom.
3. Transition metal dichalcogenides (TMDCs): MoS₂ and WS₂ monolayers exhibit FQHE when encapsulated in hexagonal boron nitride (hBN).
4. Ultracold atoms: Optical lattice systems with rotating Bose or Fermi gases have realized FQHE-like states at fractional filling.
## 12. Pair Correlation Functions and Structure Factor
A fundamental characterization of FQHE states is the pair correlation function:
$$g(r) = \frac{1}{n^2} \langle ho(0) ho(\mathbf{r}) angle,$$
where $
ho(\mathbf{r}) = \sum_i \delta(\mathbf{r} - \mathbf{r}_i)$ is the density operator. For the Laughlin ν = 1/3 state, g(r) exhibits:
1. Zero at short range (characteristic of FQHE): $g(r) \propto r^6$ as $r o 0$ (for m = 3), reflecting the $(z_i - z_j)^3$ vanishing in the wavefunction.
2. Oscillations at intermediate range: Standing out from the Gaussian envelope.
3. Approaches 1 at large r: Reflecting that at large separations, density fluctuations become uncorrelated.
The static structure factor $S(k) = \int d^2r e^{i\mathbf{k} \cdot \mathbf{r}} g(r)$ shows Bragg peaks reflecting the emergent crystalline order in real space.
## 13. Numerical Solver: Laughlin Wavefunction and Pair Correlations
The following Python code computes and visualizes the Laughlin wavefunction for small particle numbers, calculates pair correlation functions, and analyzes the many-body structure:
import numpy as np
import matplotlib.pyplot as plt
from scipy.special import erfcinv
from itertools import combinations
def magnetic_length(B_tesla):
"""Compute magnetic length from magnetic field in Tesla."""
hbar = 1.054571817e-34 # J·s
e = 1.602176634e-19 # C
l_B = np.sqrt(hbar / (e * B_tesla))
return l_B
def laughlin_wavefunction(positions, m, l_B):
"""
Compute Laughlin wavefunction at filling nu = 1/m.
positions: array of shape (N_particles, 2) with x, y coordinates
m: denominator (typically 3 for nu = 1/3)
l_B: magnetic length
"""
N = len(positions)
z = positions[:, 0] + 1j * positions[:, 1] # Complex coordinates
# Compute the Vandermonde product prod_{i<j} (z_i - z_j)^m
psi_vand = 1.0 + 0j
for i, j in combinations(range(N), 2):
psi_vand *= (z[i] - z[j])**m
# Gaussian envelope
r_squared = np.sum(positions**2, axis=1)
psi_gaussian = np.exp(-np.sum(r_squared) / (4 * l_B**2))
psi = psi_vand * psi_gaussian
return psi
def monte_carlo_density(n_samples, m, l_B, sigma=1.0):
"""
Monte Carlo sampling of electron density distribution in Laughlin state.
Uses rejection sampling on |psi|^2.
"""
N_particles = 4 # Small system for demonstration
positions_list = []
for sample in range(n_samples):
# Random initial positions
positions = np.random.randn(N_particles, 2) * l_B * sigma
# Metropolis-Hastings algorithm
for step in range(100):
# Propose random displacement
trial_positions = positions + np.random.randn(N_particles, 2) * 0.1 * l_B
# Compute acceptance ratio
psi_current = laughlin_wavefunction(positions, m, l_B)
psi_trial = laughlin_wavefunction(trial_positions, m, l_B)
prob_ratio = np.abs(psi_trial / psi_current)**2
if np.random.rand() < prob_ratio:
positions = trial_positions
positions_list.append(positions.flatten())
return np.array(positions_list)
def pair_correlation_function(positions, l_B, n_bins=50):
"""
Compute pair correlation function g(r) from particle positions.
positions: array of shape (N_samples, 2*N_particles) or (N_particles, 2)
"""
if positions.ndim == 2 and positions.shape[0] > 20:
# Assume positions are individual configurations
all_positions = positions.reshape(-1, 2)
else:
all_positions = positions
N = len(all_positions)
r_max = 5 * l_B
r_bins = np.linspace(0, r_max, n_bins)
g_r = np.zeros(n_bins - 1)
count = np.zeros(n_bins - 1)
# Compute pairwise distances
for i in range(N):
for j in range(i+1, N):
r_ij = np.linalg.norm(all_positions[i] - all_positions[j])
bin_idx = np.searchsorted(r_bins, r_ij) - 1
if 0 <= bin_idx < n_bins - 1:
g_r[bin_idx] += 2
count[bin_idx] += 2
# Normalize
density = N / (np.pi * (5*l_B)**2)
for i in range(n_bins - 1):
r_center = (r_bins[i] + r_bins[i+1]) / 2
bin_area = np.pi * (r_bins[i+1]**2 - r_bins[i]**2)
if count[i] > 0:
g_r[i] /= (count[i] * density * bin_area / (2*N))
r_centers = (r_bins[:-1] + r_bins[1:]) / 2
return r_centers, g_r
def structure_factor(positions, l_B, k_max=5):
"""
Compute static structure factor S(k).
"""
if positions.ndim == 2:
all_positions = positions
else:
all_positions = positions.reshape(-1, 2)
k_values = np.linspace(0.1, k_max / l_B, 50)
S_k = np.zeros(len(k_values))
for ik, k in enumerate(k_values):
# Compute density-density correlation in Fourier space
for i in range(len(all_positions)):
for j in range(len(all_positions)):
r_ij = all_positions[i] - all_positions[j]
S_k[ik] += np.cos(k * np.linalg.norm(r_ij))
S_k /= len(all_positions)
return k_values, S_k
def quantum_dimension_analysis(m):
"""
Compute quantum dimensions and topological properties of FQHE state.
"""
nu = 1 / m
# For Abelian FQHE at nu = 1/m:
# Total quantum dimension D = sqrt(m)
D_total = np.sqrt(m)
# Number of quasiparticles: m (quasihole types)
n_quasiparticles = m
# Fractional charge of quasihole
e_star_qh = 1 / m
# Topological entanglement entropy
S_topo = np.log(D_total)
return {
'filling_factor': nu,
'total_quantum_dimension': D_total,
'n_quasiparticles': n_quasiparticles,
'qh_fractional_charge': e_star_qh,
'topological_entropy': S_topo,
'ground_state_degeneracy': m # On a torus
}
# Parameters
B = 10 # Tesla
l_B = magnetic_length(B)
m = 3 # nu = 1/3 state
print(f"Magnetic length: {l_B:.3e} m = {l_B*1e9:.3f} nm")
# Generate Laughlin states and compute densities
n_samples_mc = 100
positions_mc = monte_carlo_density(n_samples_mc, m, l_B, sigma=2.0)
# Compute pair correlation and structure factor
r_centers, g_r = pair_correlation_function(positions_mc, l_B, n_bins=50)
k_values, S_k = structure_factor(positions_mc, l_B, k_max=5)
# Quantum properties
qprops = quantum_dimension_analysis(m)
# Create comprehensive plots
fig, axes = plt.subplots(2, 3, figsize=(16, 10))
# Panel 1: Spatial density distribution (single configuration)
ax = axes[0, 0]
x_plot = positions_mc[0, ::2]
y_plot = positions_mc[0, 1::2]
ax.scatter(x_plot / l_B, y_plot / l_B, s=100, alpha=0.7, edgecolors='black', linewidth=1.5)
ax.set_xlabel('x / l_B', fontsize=11)
ax.set_ylabel('y / l_B', fontsize=11)
ax.set_title('Laughlin ν=1/3 State: Electron Configuration', fontsize=12)
ax.set_aspect('equal')
ax.grid(True, alpha=0.3)
# Panel 2: Pair correlation function
ax = axes[0, 1]
ax.plot(r_centers / l_B, g_r, 'b-', linewidth=2.5, label='Laughlin ν=1/3')
ax.axhline(1, color='k', linestyle='--', linewidth=1, alpha=0.5, label='Uncorrelated')
ax.fill_between(r_centers / l_B, 0, g_r, alpha=0.3)
ax.set_xlabel('Pair distance r / l_B', fontsize=11)
ax.set_ylabel('Pair correlation g(r)', fontsize=11)
ax.set_title('Short-Range Correlation in FQHE Ground State', fontsize=12)
ax.set_xlim([0, 5])
ax.set_ylim([0, 1.2])
ax.legend(fontsize=10)
ax.grid(True, alpha=0.3)
# Panel 3: Structure factor
ax = axes[0, 2]
ax.plot(k_values * l_B, S_k, 'r-', linewidth=2.5)
ax.set_xlabel('Wavevector k·l_B', fontsize=11)
ax.set_ylabel('Static structure factor S(k)', fontsize=11)
ax.set_title('Fourier Space Density Correlations', fontsize=12)
ax.grid(True, alpha=0.3)
ax.set_ylim([0, np.max(S_k) * 1.1])
# Panel 4: Exchange phase for anyons
ax = axes[1, 0]
m_values = np.array([3, 5, 7, 9])
theta_exchange = np.pi / m_values
ax.scatter(m_values, theta_exchange / np.pi, s=150, alpha=0.7, edgecolors='black', linewidth=2, color='purple')
ax.plot(m_values, theta_exchange / np.pi, 'purple', alpha=0.5, linewidth=2, linestyle='--')
ax.set_xlabel('Denominator m (ν = 1/m)', fontsize=11)
ax.set_ylabel('Exchange phase θ / π', fontsize=11)
ax.set_title('Anyonic Exchange Statistics in FQHE', fontsize=12)
ax.grid(True, alpha=0.3)
for i, (m_val, theta_val) in enumerate(zip(m_values, theta_exchange)):
ax.text(m_val + 0.2, theta_val/np.pi + 0.01, f'{theta_val/np.pi:.3f}π', fontsize=9)
# Panel 5: Topological entanglement entropy
ax = axes[1, 1]
m_values = np.arange(3, 12, 2)
S_topo_values = np.log(np.sqrt(m_values))
ax.bar(m_values, S_topo_values, width=1.2, alpha=0.7, edgecolor='black', linewidth=1.5, color='green')
ax.set_xlabel('Denominator m (ν = 1/m)', fontsize=11)
ax.set_ylabel('Topological entropy S_topo', fontsize=11)
ax.set_title('Topological Order Signature', fontsize=12)
ax.grid(True, alpha=0.3, axis='y')
# Panel 6: Summary table of topological properties
ax = axes[1, 2]
ax.axis('off')
summary_text = f"""
Laughlin State ν = 1/{m}
Fractional charge (quasihole):
e* = e/{m}
Anyonic exchange phase:
θ = π/{m} = {np.pi/m:.4f} rad
Ground state degeneracy (torus):
{m}-fold
Total quantum dimension:
D = √{m} ≈ {np.sqrt(m):.3f}
Topological entropy:
S_topo = ln(D) ≈ {np.log(np.sqrt(m)):.4f}
Quasihole types:
{m}
Edge mode (Luttinger):
K = 1/{m}
"""
ax.text(0.05, 0.95, summary_text, transform=ax.transAxes, fontsize=11,
verticalalignment='top', family='monospace',
bbox=dict(boxstyle='round', facecolor='wheat', alpha=0.8))
plt.tight_layout()
plt.show()
print("
Fractional Quantum Hall Effect Analysis Complete")
print(f"Filling factor: ν = 1/{m}")
print(f"Laughlin state quantum properties:")
for key, val in qprops.items():
print(f" {key}: {val}")## 14. Hierarchical FQHE and Composite Fermion Gaps
Beyond the primary Laughlin states (ν = 1/m), the FQHE exhibits a rich hierarchy of fractional fillings. The composite fermion picture predicts that at filling factors ν = p/(2pm + 1), the ground state consists of composite fermions at integer filling p/(2m), with the FQHE replicating at each level of the hierarchy. Experimental observations confirm this structure: ν = 2/5 (two composite fermions in a filled Fermi sea at the second Landau level), ν = 2/7, ν = 3/7, etc., all exhibit precisely quantized Hall resistance and incompressible ground states.
The excitation spectrum in FQHE exhibits a gap—the minimum energy to create a quasihole-quasielectron pair (or add/remove a composite fermion). This gap scales with Coulomb interaction energy: $\Delta \sim e^2/(4\pi\epsilon l_B)$, where ε is the dielectric constant. Large gaps (meV scale for GaAs systems) enable FQHE to persist up to relatively high temperatures (~100 mK) and enable quantum Hall resistance measurements at the h/e² precision limit.
## 15. Experimental Probes and Future Directions
### 15.1 Transport Measurements
The fractional quantum Hall resistance is measured via the plateau in $R_{xy}$ (transverse resistivity) versus magnetic field, with value:
$$R_{xy} = \frac{h}{ u e^2}.$$
Integer fillings show fractional resistances, confirming the quantized Hall effect.
### 15.2 Shot Noise and Fractional Charge
Injection of electrons into a FQHE edge state yields shot noise with a Poisson distribution in charge carriers. Since the FQHE state is incompressible, injected charges tunnel in units of the fractional charge e/m, producing shot noise with reduced Poisson factor, directly measuring fractional charge.
### 15.3 Non-Abelian FQHE and Topological Quantum Computing
States like the Moore-Read state at ν = 5/2 support non-abelian anyons (Majorana fermions), making them candidates for topological quantum computing. Braiding these anyons implements quantum gates protected by topology.
### 15.4 Artificial Intelligence and Machine Learning
Neural networks and tensor network methods have been developed to compute FQHE ground states and identify topological order. These approaches complement Laughlin's trial wavefunction and composite fermion theory.
## Conclusion
The fractional quantum hall effect exemplifies how many-body quantum mechanics in a strong magnetic field can produce a ground state with profound topological properties: quantized Hall resistance, fractional excitation charges, anyonic statistics, and long-range entanglement. Laughlin's wavefunction provides an exact description for the primary FQHE states and has guided understanding of the complete phase diagram. The composite fermion picture offers an intuitive mean-field understanding, while Chern-Simons gauge theory formalizes flux attachment and topological order. The presented numerical solver demonstrates how to compute Laughlin wavefunctions and extract many-body structure via pair correlations and structure factors. Ongoing experimental work continues to explore FQHE in new material platforms (graphene, monolayer TMDCs, ultracold atoms) and search for non-abelian FQHE states suitable for topological quantum information processing.