A Genetic Algorithm Approach for Compact Wave Function Representations in Spin-Adapted Bases

The accurate treatment of many-unpaired-electron systems remains a central challenge in quantum chemistry, due to the exponential growth of the many-electron wave function with the number of correlated electrons. Quantum Anamorphosis addresses this challenge through physically motivated localization of molecular orbitals and site reordering, which yield unique block-diagonal Hamiltonian matrices and compact spin-adapted many-body wave functions. In this work, we introduce a genetic algorithm to identify optimal orbital/site orderings that enhance wave function compactness, thereby enabling the study of larger systems than previously possible. Crucially, we propose fitness functions based on approximate measures of the wave function compactness, which enable inexpensive genetic algorithm searches. We benchmark the strategy against one- and two-dimensional nearest-neighbor Heisenberg models, the one-dimensional next-nearest-neighbor Heisenberg model, and selected ground and excited states of the nitrogenase P-cluster, employing intermediate CAS(48,40) and large CAS(114,73) active space ab initio Hamiltonians. These results highlight the applicability and scalability of the genetic-algorithm-driven approach for tackling systems with many unpaired electrons. The P-cluster test case is particularly relevant, as it demonstrates that wave function compression applies to both ground and excited states, and allows selective targeting of specific low-lying states within the same spin multiplicity.