using-xdiag

Exact diagonalization in Julia is performed using the XDiag.jl stack to construct Hilbert-space blocks and compute ground states via matrix-free eigensolvers. It handles symmetry-adapted basis constraints and operator summation for quantum lattice systems.

Lanczos and LOBPCG solvers are supported for quantum many-body ground state calculations. The environment also provides Krylov methods for computing excited states and performing real or imaginary time evolution dynamical calculations.

Symmetry-adapted exact diagonalization handles U(1) charge sectors by constructing memory-efficient basis encodings and applying spatial permutation symmetries. This ensures scientifically rigorous Hilbert-space block construction for quantum many-body lattice systems.

A valid Julia environment configuration is required to run exact diagonalization computations. You must adhere to the symmetry-adapted basis constraints defined by the method card to ensure correct operator summation and matrix-free eigensolver execution.

Computing dynamical spectral functions for a Heisenberg chain is best done using exact diagonalization with advanced solvers like Lanczos. This approach ensures memory-efficient basis encoding and correct symmetry handling for accurate spectral calculations.

Matrix-free eigensolvers are used for quantum many-body lattice systems to compute ground states and dynamics efficiently without explicitly storing large Hamiltonian matrices. This approach enables exact diagonalization of larger system sizes within memory constraints.