Gauge-Covariant Husimi Q-function

• The gauge-covariant Husimi Q-function is a positive quasi-probability distribution defined on phase space that incorporates gauge invariance using Wilson lines and covariant kernels.
• It is constructed using gauge-covariant coherent states or by smoothing the Wigner function, enabling strict invariance under local gauge transformations in charged quantum systems.
• This formulation facilitates the analysis of physical observables such as flux, thermodynamic quantities, and phase-space localization in systems like Landau levels and curved geometries.

The gauge-covariant Husimi Q-function is a positive quasi-probability distribution formulated on phase space that generalizes the standard Husimi Q-function to charged quantum systems in the presence of gauge fields, ensuring strict invariance under local gauge transformations. It achieves this by systematically modifying either the coherent-state construction or the Wigner-to-Husimi transformation to incorporate the gauge structure, notably through the explicit inclusion of Wilson lines, magnetic translation operators, or covariant kernels. This construction is central for phase-space analysis of quantum dynamics in electromagnetic backgrounds, non-Abelian gauge theories, and systems of Landau levels in both Euclidean and curved (hyperbolic) spaces.

1. Formal Construction and Gauge Invariance

The gauge-covariant Husimi Q-function $Q(\mathbf{r},\mathbf{p})$ is constructed so as to remain strictly invariant under gauge transformations $\mathbf{A} \to \mathbf{A} + \nabla\chi$, with the quantum state (density matrix) transforming as $$\hat\rho \to e^{i q \chi(\hat{\mathbf{r}})} \hat\rho e^{-i q \chi(\hat{\mathbf{r}})}$$ for charge $q$. Two dominant approaches exist:

a) Gauge-Covariant Coherent States Approach:

For a charged particle in a static vector potential $A(x)$, the minimal-uncertainty coherent state centered at $(x_0, p_0)$ incorporates a Wilson-line phase: $$\langle x|x_0, p_0;\sigma\rangle_A = N(\sigma)\exp\left\{-\frac{(x - x_0)^2}{4\sigma^2} + i p_0 \cdot (x - x_0) + i q \int_{x_0}^x A(\xi) \cdot d\xi\right\}$$ This state transforms covariantly, picking up only a local (unphysical) phase under a gauge transformation. The Q-function is defined as

$$Q_A(x_0,p_0) = |\langle x_0,p_0;\sigma|\psi\rangle|^2$$

which is strictly gauge-invariant owing to cancellation of the gauge phase factors in the overlap and its modulus squared (Mason et al., 2012, Datseris et al., 2019, Korennoy, 2018).

b) Gauge-Covariant Smoothing of the Wigner Function:

The Stratonovich-Wigner function $W_g(\mathbf{r},\mathbf{p},t)$ incorporates a straight-line Wilson phase: $$W_g(\mathbf{r},\mathbf{p},t) = \frac{1}{(2\pi\hbar)^3} \int d^3u \exp\left\{\frac{i}{\hbar} \mathbf{u}\cdot\left[\mathbf{p} + \frac{e}{c}\int_{-1/2}^{1/2}\!d\tau\,\mathbf{A}(\mathbf{r} + \tau \mathbf{u}, t)\right]\right\} \rho(\mathbf{r} - \tfrac{\mathbf{u}}{2}, \mathbf{r} + \tfrac{\mathbf{u}}{2}, t)$$ The Husimi Q-function $\mathbf{A} \to \mathbf{A} + \nabla\chi$0 is then obtained by Gaussian smoothing: $\mathbf{A} \to \mathbf{A} + \nabla\chi$1 Guaranteeing positivity and gauge invariance (Korennoy, 2018).

2. Fundamental Properties and Transformation Rules

The gauge-covariant Husimi Q-function satisfies:

• Positivity: $\mathbf{A} \to \mathbf{A} + \nabla\chi$2 everywhere due to the smoothing of the (generally non-positive) Wigner function by a Gaussian kernel.
• Normalization: $\mathbf{A} \to \mathbf{A} + \nabla\chi$3 (for pure state) or $\mathbf{A} \to \mathbf{A} + \nabla\chi$4 (general density matrix).
• Gauge Invariance: Under $\mathbf{A} \to \mathbf{A} + \nabla\chi$5, and $\mathbf{A} \to \mathbf{A} + \nabla\chi$6, the Q-function remains unchanged due to explicit cancellation of all induced local gauge phases.
• Classical Limit: As $\mathbf{A} \to \mathbf{A} + \nabla\chi$7, $\mathbf{A} \to \mathbf{A} + \nabla\chi$8 reduces to a classical probability distribution; the evolution equation for $\mathbf{A} \to \mathbf{A} + \nabla\chi$9 recovers the Liouville equation with Lorentz force contributions (Korennoy, 2018).
• Phase-Space Localization: For Landau levels and other eigenstates, the Q-function is sharply peaked in phase space at the classical guiding center orbits.

For non-Abelian gauge fields (Yang-Mills), the gauge-covariant Husimi functional for field variables $$\hat\rho \to e^{i q \chi(\hat{\mathbf{r}})} \hat\rho e^{-i q \chi(\hat{\mathbf{r}})}$$0 is obtained by Gaussian coarse-graining the Wigner functional, utilizing only gauge-invariant combinations such as $$\hat\rho \to e^{i q \chi(\hat{\mathbf{r}})} \hat\rho e^{-i q \chi(\hat{\mathbf{r}})}$$1, ensuring invariance under local gauge rotations (Tsukiji et al., 2016).

3. Relation to Physical Observables: Flux and Thermodynamic Quantities

The gauge-covariant Q-function encodes local phase-space information, and its moments are directly related to physical observables:

• Covariant Probability Current: In the limit $$\hat\rho \to e^{i q \chi(\hat{\mathbf{r}})} \hat\rho e^{-i q \chi(\hat{\mathbf{r}})}$$2, the Husimi Q-function recovers the standard covariant flux operator:

$$\hat\rho \to e^{i q \chi(\hat{\mathbf{r}})} \hat\rho e^{-i q \chi(\hat{\mathbf{r}})}$$3

where $$\hat\rho \to e^{i q \chi(\hat{\mathbf{r}})} \hat\rho e^{-i q \chi(\hat{\mathbf{r}})}$$4 is the gauge-invariant probability current at $$\hat\rho \to e^{i q \chi(\hat{\mathbf{r}})} \hat\rho e^{-i q \chi(\hat{\mathbf{r}})}$$5 (Mason et al., 2012).

• Thermodynamic Quantities: For thermal (Gibbs) states, the Husimi Q-function yields explicit expressions for mean phase-space distributions, variances, and provides lower bounds for the grand canonical thermodynamic potential via the Berezin–Lieb inequality (Mouayn et al., 2021).

4. Explicit Realizations: Electrons, Landau Levels, and Curved Geometries

a) Ballistic Electrons in Magnetic Fields:

In the analysis of electronic transport under magnetic fields, Husimi Q-functions defined using magnetic translation operators are essential to maintain gauge invariance. The peaks of the Q-function follow classical cyclotron orbits, caustics, and edge states, enabling visualization of quantum transport features such as Klein tunneling or skipping orbits in graphene nanodevices (Datseris et al., 2019).

b) Landau and Hyperbolic Landau Levels:

For planar and hyperbolic (Poincaré disk) Landau problems, gauge-covariant coherent states are constructed either via Wilson lines or explicit gauge phases $$\hat\rho \to e^{i q \chi(\hat{\mathbf{r}})} \hat\rho e^{-i q \chi(\hat{\mathbf{r}})}$$6. The resulting Q-functions possess closed forms involving Jacobi or Laguerre polynomials, Kampé de Fériet functions, and yield characteristic functions and moments analytically. In the flat ($$\hat\rho \to e^{i q \chi(\hat{\mathbf{r}})} \hat\rho e^{-i q \chi(\hat{\mathbf{r}})}$$7) limit, these results reduce to the standard Landau-level Q-functions on the complex plane (Mouayn et al., 2021).

c) Non-Abelian Gauge Fields:

For Yang-Mills theory on the lattice, the Wigner functional is defined in terms of $$\hat\rho \to e^{i q \chi(\hat{\mathbf{r}})} \hat\rho e^{-i q \chi(\hat{\mathbf{r}})}$$8 and $$\hat\rho \to e^{i q \chi(\hat{\mathbf{r}})} \hat\rho e^{-i q \chi(\hat{\mathbf{r}})}$$9, and the gauge-covariant Husimi functional is constructed by coarse-graining with a product of local gauge-invariant Gaussians. This enables the computation of Husimi-Wehrl entropy and its time evolution in semiclassical approximations (Tsukiji et al., 2016).

5. Evolution Equations and Operator Approach

The time evolution of the gauge-covariant Husimi Q-function is governed by an exact Moyal-type equation, involving “smeared” electromagnetic (or chromodynamic) fields:

$q$0

Here $q$1 are quantum-corrected (smeared) fields. In the classical limit, this reduces to the Liouville equation (Korennoy, 2018).

All phase-space quasiprobabilities of this type are operator traces: $q$2 with explicit dequantizer and quantizer operators furnished in both coordinate and Weyl-operator forms (Korennoy, 2018).

6. Computational and Practical Aspects

Numerical evaluation of gauge-covariant Husimi Q-functions in high-dimensional field theories employs efficient techniques:

• Product Ansatz: Factorizing the multi-field Husimi functional over degrees of freedom, providing a controlled overestimate of entropy measures (10–20% for moderate dimensions) (Tsukiji et al., 2016).
• Test-Particle and Parallel Test-Particle Methods: Monte Carlo ensembles of classical field configurations (test particles) effectively sample the Wigner and hence Husimi functionals, enabling dynamical studies in the semiclassical regime.
• Analytic and Statistical Parameters: Closed forms for mean, variance, and characteristic functions in Landau and hyperbolic Landau models enable explicit statistical analysis, including thermodynamic limits and semi-classical approximations (Mouayn et al., 2021).

7. Alternative Formulations and Generalizations

Alternative gauge-invariant Husimi functions arise by modifying the manner in which the Wilson phase is distributed (e.g., “non-Stratonovich” options). While these approaches remain positive and gauge-invariant, they differ in operator-ordering ambiguities or the distribution of quantum corrections, merging in the semi-classical or classical limit (Korennoy, 2018).

The formalism is fully compatible with generalizations to time-dependent fields, higher dimensions, and non-Abelian structures, as long as gauge-invariant smearing and projection is respected.

---

References:

• "Extending the Concept of Probability Flux" (Mason et al., 2012)
• "Gauge-independent Husimi functions of charged quantum particles in the electro-magnetic field" (Korennoy, 2018)
• "Husimi function for electrons moving in magnetic fields" (Datseris et al., 2019)
• "Husumi Q-functions attached to hyperbolic Landau levels" (Mouayn et al., 2021)
• "Entropy production from chaoticity in Yang-Mills field theory with use of the Husimi function" (Tsukiji et al., 2016)