Many-electron systems with fractional electron number and spin: exact properties above and below the equilibrium total spin value

This paper rigorously derives the exact properties of the ensemble ground state for many-electron systems with fractional electron numbers and spin projections, resolving ground-state ambiguities in low-spin cases via entropy maximization, characterizing high-spin dependencies, and establishing generalized ionization potential theorems and new derivative discontinuities to advance density functional theory approximations.

The Gist

Imagine you are trying to describe a crowd of people (electrons) in a room. In the world of quantum physics, things get weird: sometimes the crowd doesn't have a whole number of people (maybe it's 5.8 people), and sometimes the people are spinning in different directions (spin).

This paper is like a rulebook for figuring out exactly what this "fractional crowd" looks like when it's in its most relaxed, lowest-energy state. The authors, Yuli Goshen and Eli Kraisler, are solving a puzzle that has been tricky for scientists trying to predict how materials behave.

Here is the breakdown using simple analogies:

1. The Problem: The "Fuzzy" Crowd

In standard physics, we usually deal with whole numbers of electrons (1, 2, 3...). But in real-world chemistry and materials science, we often need to describe systems that are "in between" states, like an atom that has lost a tiny bit of an electron, or a magnet that is partially flipped.

Think of it like a smoothie. If you have a whole apple (1 electron) and a whole orange (1 electron), you have 2 fruits. But if you blend them, you have a mix. The paper asks: If I have 5.8 electrons, what is the exact recipe of the "smoothie" that makes the system happiest (lowest energy)?

2. The "Low Spin" Zone: The Traffic Light

The authors divide the problem into two zones based on how much the electrons are "spinning" (magnetic alignment).

Zone 1: The Calm Zone (Low Spin)
Imagine a traffic light that can be Red, Yellow, or Green.

• The Old Rule: Scientists knew that if you have a fractional electron, the system is a mix of the state with the whole number below (Red) and the whole number above (Green).
• The New Discovery: The authors found that in this calm zone, there is a mystery. There isn't just one way to mix the Red and Green states to get the Yellow state. There are many ways to do it, and they all result in the same total energy. It's like having many different recipes for a smoothie that all taste exactly the same.
• The Solution: To pick the one correct recipe, the authors suggest using a concept called Entropy. Think of entropy as "disorder" or "freedom." Nature loves to be as free and disordered as possible. By choosing the mix that has the maximum freedom (maximum entropy), they found a unique, perfect answer. It's like saying, "If all these smoothie recipes taste the same, let's pick the one that uses the most variety of fruits."

3. The "High Spin" Zone: The Wild West

Zone 2: The Chaotic Zone (High Spin)
Now, imagine the electrons are spinning wildly, like a mosh pit.

• The Rule: In this chaotic zone, the "smoothie" recipe changes depending on the specific type of atom (the system). There is no single universal rule like in the calm zone.
• The Discovery: The authors proved three things:
  1. The mix only uses specific "pure" states that are on the edge of the chaos.
  2. The mix is made of at most three specific ingredients (pure states).
  3. As you slowly add more spin, the third ingredient stays the same for a while, even if you change the total number of electrons slightly.
• The Analogy: It's like building a tower with blocks. In the calm zone, you can use any blocks. In the wild zone, you can only use three specific blocks, and the shape of the tower depends entirely on the specific type of blocks you have (the specific atom).

4. The "Jump" in the Energy (The Bumpy Road)

One of the most exciting parts of the paper is about Derivative Discontinuities.

Imagine you are driving a car on a road that represents the energy of the system.

• The Old View: The road was a smooth ramp. As you added electrons, the energy went up or down smoothly.
• The New View: The authors show that the road has sudden jumps or cliffs. When you cross a specific line (like when the spin reaches a certain limit), the "slope" of the road changes instantly.
• Why it matters: In the language of the computer programs scientists use (called DFT), these jumps mean the "virtual landscape" the electrons see changes abruptly. It's like driving over a speed bump that you didn't see coming. The authors provide a map of exactly where these bumps are and how high they are.

5. Why Should You Care?

This might sound very technical, but it's the foundation for better technology.

• Better Batteries and Solar Cells: To design better materials, scientists need to predict how electrons move and interact. Current computer models often get these "fractional" and "spinning" situations wrong, leading to inaccurate predictions.
• The Fix: By giving these computer models the new "rules of the road" (the exact properties derived in this paper), scientists can build much more accurate simulations. This means we can design new materials for electronics, medicine, and energy without having to build them in a lab first.

Summary

• The Goal: Figure out what a system of electrons looks like when it has a "fractional" amount of charge and spin.
• The Calm Zone: There are many answers, but the "most free" (maximum entropy) one is the right one.
• The Chaotic Zone: The answer depends on the specific atom, but it's limited to a mix of three specific states.
• The Result: A new map of "energy jumps" that helps computer programs predict how materials behave much more accurately.

In short, the authors have written a more precise instruction manual for the quantum world, helping us understand the "in-between" states of matter that are crucial for the future of technology.

Technical Summary

Here is a detailed technical summary of the paper "Many-electron systems with fractional electron number and spin: exact properties above and below the equilibrium total spin value" by Yuli Goshen and Eli Kraisler.

1. Problem Statement

The paper addresses the fundamental theoretical challenge of describing the ensemble ground state of general, finite, many-electron systems at zero temperature when both the total electron number ($N_{tot}$) and the total $z$-projection of the spin ($M_{tot}$) are fractional.

While the behavior of systems with fractional electron numbers is well-understood (piecewise linearity of energy), the simultaneous variation of fractional spin introduces complex ambiguities and structural dependencies. Specifically:

• Low-Spin Regime: Within a specific region defined by a boundary spin $M_B$ (Region 1), the ground state form was previously derived, but issues regarding the non-uniqueness of the ensemble weights (and consequently the spin density) remained.
• High-Spin Regime: Outside this region ($|M_{tot}| > M_B$), the structure of the ground state ensemble was not rigorously characterized, and its dependence on specific system properties was unclear.
• DFT Implications: The lack of exact conditions for these regimes limits the accuracy of Density Functional Theory (DFT), particularly regarding derivative discontinuities in Kohn-Sham (KS) potentials and the prediction of ionization potentials (IP) and fundamental gaps for spin-polarized systems.

2. Methodology

The authors employ a rigorous theoretical framework based on ensemble density functional theory and variational principles, supported by extensive numerical analysis.

• Variational Derivation: The authors minimize the total energy expectation value $E_{tot} = \text{Tr}\{\hat{\Lambda}\hat{H}\}$ subject to constraints on $N_{tot}$ and $M_{tot}$. They utilize the convexity conjecture for energy sequences (ionization potentials and electron affinities) to prove which pure states contribute to the ground-state ensemble.
• Entropy Maximization: To resolve the ambiguity in the low-spin regime, the authors propose maximizing the von Neumann entropy ($S = -\text{Tr}\{\hat{\Lambda} \ln \hat{\Lambda}\}$) as a physical criterion to uniquely determine the ensemble weights.
• Geometric Analysis: The problem is mapped onto the $N_\uparrow - N_\downarrow$ plane. The authors analyze the geometry of "tiles" (triangular regions) formed by pure states to determine the ensemble composition in the high-spin regime.
• Kohn-Sham DFT Connection: They relate the derived ensemble properties to KS eigenvalues ($\varepsilon^\sigma_{ho}$) and potentials ($v^\sigma_{KS}$), deriving expressions for derivative discontinuities at boundaries where the energy slope changes.
• Numerical Validation: The theoretical predictions are validated against the National Institute of Standards and Technology (NIST) Atomic Spectra Database for atoms with atomic numbers $1 \le Z \le 54$.

3. Key Contributions and Results

A. Low-Spin Regime ($|M_{tot}| \le M_B$)

• Alternative Proof: The authors provide an alternative, rigorous proof for the form of the ensemble ground state, confirming it is a linear combination of ground states with $N_0$ and $N_0+1$ electrons.
• Resolution of Non-Uniqueness: They demonstrate that for fractional $N_{tot}$ and $M_{tot}$, the spin densities $n_\sigma(\mathbf{r})$ are not uniquely determined by $N_{tot}$ and $M_{tot}$ alone.
  • They show that maximizing the system's entropy uniquely determines the statistical weights $\lambda_{N,M}$.
  • This results in a unique ensemble where weights follow a geometric progression ($\lambda_{N,M} \propto z^M$), distinct from other criteria like minimizing spin variance ($\Delta S_z$).
• Piecewise Linearity: Confirmed that total energy and any operator commuting with the spin ladder operator are piecewise-linear in $N_{tot}$ and independent of $M_{tot}$ within this region.

B. High-Spin Regime ($|M_{tot}| > M_B$)

• Three General Properties: The authors prove three exact properties characterizing the high-spin ground state:
  1. State Selection: Only pure states with $M \ge S_{min}(N)$ contribute to the ensemble (states must lie outside or on the edge of Region 1).
  2. Triplet Composition: The ensemble consists of at most three pure states (forming a triangle in the $N_\uparrow - N_\downarrow$ plane), provided no accidental degeneracy exists.
  3. Boundary Behavior: As $M_{tot}$ approaches $M_B$ from above, two of the three states are fixed ($|N_0, S_0\rangle$ and $|N_0+1, S_1\rangle$), while the third state is determined by the specific system's energy landscape but is independent of the fractional part $\alpha$.
• System Dependence: Unlike the low-spin case, the high-spin ensemble structure is system-dependent. The specific third state is chosen to minimize the energy, often involving spin-flip energies and fundamental gaps.

C. Kohn-Sham DFT Implications

• Generalized IP Theorem: The authors extend the IP theorem to fractional $N_{tot}$ and high-spin cases. They derive exact relations between KS frontier orbital energies and total energy differences (IP, spin-flip energies).
• New Derivative Discontinuities: They identify that crossing boundaries between different "tiles" in the energy profile causes abrupt jumps in KS eigenvalues.
  • This necessitates derivative discontinuities ($\Delta^\sigma$) in the KS potentials $v^\sigma_{KS}(\mathbf{r})$.
  • These jumps can occur even when $N_{tot}$ is fractional, provided the boundary is not parallel to the axes in the $N_\uparrow - N_\downarrow$ plane.

D. Numerical Analysis

• Using NIST data, the authors constructed ensemble maps for various atoms (He, C, Fe, etc.).
• Non-Convexity: They identified rare cases (e.g., $Fe^{7+}$) where the pure state energy is not convex with respect to $N_\sigma$, leading to ground states that are ensembles of integer states rather than a single pure state.
• Tile Structures: The analysis confirmed that ensemble regions are not always simple right triangles; complex "tiles" appear depending on the relative magnitudes of spin-flip energies and fundamental gaps.

4. Significance

• Theoretical Foundation: The work provides a complete, exact description of the ensemble ground state for fractional electron numbers and spins, filling a critical gap in many-body quantum theory.
• DFT Development: The derived exact conditions (piecewise linearity, flat-plane conditions, and new derivative discontinuities) serve as rigorous constraints for developing next-generation exchange-correlation functionals. Current approximations often fail to capture these spin-dependent discontinuities.
• Predictive Power: By generalizing the IP theorem and defining the behavior of KS potentials at high spins, the paper offers a pathway to more accurate predictions of ionization potentials, fundamental gaps, and spin-related properties in materials science and chemistry.
• Entropy Principle: The proposal to use entropy maximization to resolve ground-state ambiguity offers a physically motivated and unique solution, distinguishing it from ad-hoc constraints.

In summary, this paper establishes a comprehensive theoretical framework for fractional electron-spin systems, proving that while low-spin ensembles are universal and uniquely defined by entropy, high-spin ensembles are system-specific and governed by a strict set of geometric and energetic rules, with profound implications for the accuracy of electronic structure calculations.