Projective spin adaptation for the exact diagonalization of isotropic spin clusters

Shadan Ghassemi Tabrizi,Thomas D. Kühne

2024-08-09

Abstract:Spin Hamiltonians, like the Heisenberg model, are used to describe magnetic properties of exchange-coupled molecules and solids. For finite clusters, physical quantities such as heat capacities, magnetic susceptibilities or neutron-scattering spectra, can be calculated based on energies and eigenstates obtained by exact diagonalization (ED). Utilizing spinrotational symmetry SU(2) to factor the Hamiltonian with respect to total spin S facilitates ED, but the conventional approach to spin-adapting the basis is more intricate than selecting states with a given magnetic quantum number M (the spin z-component), as it relies on irreducible tensor-operator techniques and spin-coupling coefficients. Here, we present a simpler technique based on applying a spin projector to uncoupled basis states. As an alternative to Löwdin's projection operator, we consider a group-theoretical formulation of the projector, which can be evaluated either exactly or approximately using an integration grid. An important aspect is the choice of uncoupled basis states. We present an extension of Löwdin's theorem for s = 1/2 to arbitrary local spin quantum numbers s, which allows for the direct selection of configurations that span a complete, linearly independent basis in an S sector upon spin projection. We illustrate the procedure with a few examples.

Strongly Correlated Electrons

What problem does this paper attempt to address?

The paper is primarily dedicated to addressing the issue of spin adaptation when dealing with isotropic spin clusters in the exact diagonalization method. Specifically, it proposes a new method based on spin projection to simplify the traditionally complex spin adaptation techniques. In detail, the authors aim to:

1. **Simplify the spin adaptation process**: Traditional spin adaptation techniques rely on Irreducible Tensor Operator (ITO) techniques and spin coupling coefficients, which make the calculations complex. This paper proposes a new method based on projecting uncoupled basis states through spin projectors, avoiding the complex ITO techniques.

2. **Select appropriate uncoupled basis states**: To ensure that the set of states generated after spin projection is complete and linearly independent, it is necessary to correctly select the initial uncoupled basis states. The paper extends Löwdin's theorem to apply to cases with arbitrary local spin quantum numbers and proposes an intuitive method to select these states.

3. **Implement spin projection**: The paper discusses two techniques for implementing spin projection: one based on the Löwdin projection operator method and the other based on the group theory projection operator method, the latter of which can be solved directly without discretizing the relevant integrals. Additionally, it explores how to improve numerical stability by constructing an optimally conditioned basis.

4. **Verify the effectiveness and accuracy of the method**: Specific examples are provided to demonstrate how to select basis states according to the extended Löwdin theorem and compare the numerical accuracy of different projection methods. The results show that the method based on group theory projection operators provides higher accuracy.

In summary, this paper aims to simplify the spin adaptation process in exact diagonalization calculations and improve computational efficiency and numerical stability by introducing a new spin projection method. This is of significant importance for understanding the properties of magnetic materials such as multinuclear transition metal complexes.