Fermionic Star-Exponential

• Fermionic star-exponential is a central construct in deformation quantization for fermionic systems, enabling closed-form expressions for quantum evolution and ground-state energies.
• It employs Grassmann coherent-state path integrals and the Stratonovich–Weyl quantization map to derive explicit formulas for time evolution operators and propagators.
• The approach offers practical computational techniques for simple and driven fermionic oscillators, effectively bypassing traditional operator-based methods.

The fermionic star-exponential is a central object in the deformation-quantization (Weyl–Wigner–Moyal, WWM) approach for fermionic quantum systems, representing the symbol of the quantum evolution operator in a phase space parameterized by Grassmann (anticommuting) variables. It generalizes the bosonic star-exponential, allowing direct computation of quantum dynamics entirely within the formalism of functions on fermionic phase space via the star-product. The rigorous construction of the fermionic star-exponential, enabling closed-form expressions for quantum time evolution and ground-state energies, relies on fermionic coherent-state path integrals and the inversion of the Stratonovich–Weyl quantization map. This framework notably yields a phase-space version of the Feynman–Kac formula for fermions, bypassing operator-based techniques and proving effective in explicit computations for simple and driven fermionic oscillators (Kafuri, 30 Jan 2026).

1. Algebraic Structure: Star-Product and Star-Exponential

On a $2n$-dimensional fermionic phase space with coordinates $(\eta_j, v_j)$, satisfying anticommutation relations $\{\eta_j,\eta_k\}_+=\{v_j,v_k\}_+=0$, $\{\eta_j,v_k\}_+=\hbar\,\delta_{jk}$, the deformation-quantization star-product for smooth functions $f(\eta,v)$ and $g(\eta,v)$ is defined as

$$(f * g)(\eta,v) = f(\eta,v) \exp\left[\frac{i}{\hbar} \sum_{j=1}^n \big( \overleftarrow{\partial}_{\eta_j} \overrightarrow{\partial}_{v_j} + \overleftarrow{\partial}_{v_j} \overrightarrow{\partial}_{\eta_j} \big) \right] g(\eta,v).$$

This differential form is equivalent to an integral representation over Grassmann variables: $$(f * g)(\eta,v) = 2^{-n} \int d\eta' dv' d\xi du\, f(\eta+\eta',v+v')\,g(\eta+\xi, v+u) \exp\left\{ -\frac{2}{\hbar} \sum_j \left[\eta_j' u_j + \xi_j v_j' \right] \right\}.$$ The fermionic star-exponential of a Hamiltonian symbol $H(\eta, v)$ is the formal series

$$\exp_{*}\left(-\frac{i}{\hbar} t H\right) = \sum_{k=0}^{\infty} \frac{1}{k!} \left(-\frac{it}{\hbar}\right)^k H^{*k}$$

with $H^{*k}$ indicating the $k$-fold star-product.

2. Grassmann Coherent States and Path Integral Symbol Calculus

Fermionic coherent states are defined as

$$|v\rangle = \exp\left\{\frac{1}{\hbar} \sum_j v_j \hat\psi_j^\dagger \right\} |0\rangle, \qquad \langle \eta| = \langle 0| \exp\left\{\frac{1}{\hbar}\sum_j \eta_j \hat\psi_j\right\}$$

and satisfy

$$\int d\eta dv \exp\left\{-\frac{1}{\hbar} \eta \cdot v\right\} |v\rangle \langle \eta| = \mathbb{I}, \quad \eta\cdot v = \sum_j \eta_j v_j.$$

The time evolution kernel (propagator) in the $v$-basis is

$$K(v_f,t;v_i,0) = \langle v_f|\mathcal{T}\exp\left\{-\frac{i}{\hbar} \int_0^t dt'\, \hat H\right\}|v_i\rangle$$

with a path-integral representation,

$$K(v_f, t; v_i, 0) = \int \mathcal{D}\eta\, \mathcal{D}v\, \exp\left\{ \frac{i}{\hbar} \int_0^t dt' [ i\eta \cdot \dot v - H(\eta, v) ] \right\}$$

where $\mathcal{D}\eta\, \mathcal{D}v$ denotes integration over antiperiodic Grassmann trajectories.

By inverting the Stratonovich–Weyl quantization map, the star-exponential is expressed as

$$\exp_{*}\left(-\frac{i}{\hbar} t H(\eta, v)\right) = 2^n \int d\Delta\eta\, d\Delta v\; \exp\left\{ -\frac{2}{\hbar} \eta \cdot \Delta v \right\} K(v+\Delta v, t; v-\Delta v, 0).$$

3. Closed-Form Solutions: Meticulous Formula and Practical Computation

For the single-mode fermionic oscillator with $H(\eta, v)=\hbar\omega\,\eta v$, explicit integration yields the closed-form meticulous formula: $$\boxed{ \exp_{*}\left(-\frac{i}{\hbar} t H(\eta, v)\right) = \exp\left( -\frac{i}{2} \omega t \right) \exp\left[ \frac{2}{\hbar} e^{-i\omega t} \eta v \right] \exp\left[ \frac{2}{\hbar} (1 - e^{-i\omega t}) v\eta \right] }$$ The derivation involves evaluating the path-integral for the propagator, variable changes, reduction to a double Grassmann Gaussian integral, and completing the square, concluding with application of the general formula for fermionic Gaussian integrals. For $n$ modes, a fully analogous result is obtained with the exterior algebra structure $\eta v \in \wedge^{2n}$ considered.

4. Fermionic Feynman–Kac Formula in Phase Space

Expanding in the basis of star-eigenfunctions and performing a Wick rotation ($t\rightarrow -iT$), the ground-state energy $E_0$ of a system with nondegenerate ground state is found by

$$E_0 = -\lim_{T\to+\infty} \frac{1}{T} \ln\left[ \int d\eta\, dv\, \exp_{*}\big( -T H(\eta, v) \big) \right]$$

For the oscillator $H = \hbar\omega \eta v$, the integration yields

$$\int d\eta\,dv\,\exp_{*}( - T \hbar\omega \eta v ) = \exp\left( -\frac{\omega T}{2} \right) (1 + e^{-\omega T} )$$

and hence

$E_0 = \frac{\omega}{2},$

reproducing the expected zero-point energy purely within the phase-space framework.

5. Exemplary Applications: Simple and Driven Oscillators

For the simple fermionic oscillator ($n=1$, $H=\hbar\omega[\eta v - \tfrac{1}{2}]$), the meticulous star-exponential and Feynman–Kac formula reproduce the correct $E_0 = \omega/2$. For a driven oscillator with Hamiltonian $$H_{\rm drv}(\eta,v) = \hbar\omega \eta v + a\,\eta + a^* v$$ (with Grassmann-valued force $a$), the propagator is computed using the Heisenberg equation or generating functionals, inserted into the integral star-exponential formula, and yields a closed-form result: $$\exp_{*}\left(-\frac{i}{\hbar}t\, H_{\rm drv}\right) = \exp\left(-\frac{i}{2} \omega t\right) \exp\left[ \frac{2}{\hbar} e^{-i\omega t} ( \eta v + F_1(a,t) ) \right] \exp\left[ \frac{2}{\hbar} (1 - e^{-i\omega t}) ( v\eta + F_2(a,t) ) \right],$$ where $F_{1,2}(a,t)$ are certain two-point source-dependent functions. The Feynman–Kac formula yields the eigenvalues $E_\pm = \pm \tfrac{1}{2} \sqrt{\omega^2 + 4|a|^2}$, and in the weak-coupling limit $|a|\ll \omega$, the “naive” star-exponential recovered by ad-hoc prescriptions emerges as the leading-order approximation.

6. Methodological Considerations and Convergence

Several methodological points warrant emphasis:

• The convergence of the formal star-exponential series is an unresolved question in general; the propagator integral representation circumvents practical issues.
• Weyl-symmetric (Stratonovich–Weyl) ordering ensures consistency between operator multiplication and the phase-space star-product.
• Grassmann integrals are treated according to strict rules of parity and the Berezin calculus: $\int d\eta\,\eta=1$, $\int d\eta\,1=0$, with Gaussian integration yielding determinants in the numerator.
• Boundary condition “remediation” via Grassmann delta functions $$\delta(\eta-\eta') := \int d\xi\, \exp[i \xi (\eta-\eta')]$$ is only necessary in naive approaches and is entirely subsumed within the full coherent-state path integral formalism.

7. Significance and Computational Outlook

The amalgamation of Grassmann coherent-state path integrals and the Stratonovich–Weyl quantization map provides a robust method for deriving explicit, closed-form star-exponentials in fermionic quantum systems. This enables direct computation of propagators, spectral properties, and ground-state energies—exemplified by the fermionic Feynman–Kac formula—within deformation quantization, without recourse to operator techniques. The methodology is immediately applicable to not only canonical quadratic (oscillator-type) Hamiltonians but also more general situations where the propagator can be constructed or approximated, with the rigorous (meticulous) scheme encompassing, as a limiting case, naive ad-hoc approaches for weak coupling (Kafuri, 30 Jan 2026).