Homodyne Nonclassical Area Study

• Homodyne nonclassical area is defined as an integrated scalar quantity extracted from phase-sensitive homodyne data that signals nonclassical behavior.
• Multiple formulations, including tomographic, click-based, and phase-resolved approaches, use quadrature variances and negativities to certify nonclassicality.
• Experimental protocols apply numerical integration with efficiency and thermal corrections to robustly measure nonclassical area in single and multimode optical setups.

Searching arXiv for the cited papers and related homodyne nonclassical area literature. Homodyne nonclassical area denotes a family of integrated quantities extracted from homodyne data that quantify nonclassical behavior either in optical phase space, in the local-oscillator phase domain, or directly from optical tomograms. The term is not used in a single universal sense across the literature. In one line of work it refers to the integrated negative part of a nonclassicality quasiprobability reconstructed from balanced homodyne detection; in another it denotes an area obtained from phase-resolved click-homodyne witnesses; and in a tomographic formulation it is the excess area traced by the quadrature standard deviation over the local-oscillator phase relative to the coherent-state reference. These constructions are united by a common principle: homodyne measurements provide phase-sensitive marginals or directly sampled phase-space functionals from which a scalar, integrated witness of nonclassicality can be formed without requiring a unique canonical quasiprobability (Kühn et al., 2014).

1. Terminological scope and conceptual setting

The most widely cited explicit use of the term is the tomographic one. For a single optical mode with rotated quadrature

$$\hat X_\theta \equiv \frac{\hat a\,e^{-i\theta}+\hat a^\dagger e^{i\theta}}{\sqrt{2}},$$

and optical tomogram

$$w(X,\theta)=\langle X,\theta|\hat\rho|X,\theta\rangle,$$

the quadrature moments follow from

$$\langle \hat X_\theta^n\rangle=\int_{-\infty}^{\infty}dX\,X^n\,w(X,\theta),$$

with variance

$$\mathrm{Var}(\hat X_\theta)=\langle \hat X_\theta^2\rangle-\langle \hat X_\theta\rangle^2,$$

and standard deviation

$$\Delta X_\theta=\sqrt{\mathrm{Var}(\hat X_\theta)}.$$

The homodyne nonclassical area is then defined as

$$\sigma(\rho):=\int_0^{2\pi}d\theta\,\Delta X_\theta,$$

with coherent-state reference $\sqrt{2}\pi$, so that $\mathrm{NA}\equiv \sigma(\rho)-\sqrt{2}\pi$ (Rohith et al., 2018).

A broader usage emerges in related homodyne literatures. In balanced-homodyne reconstructions of regularized Glauber–Sudarshan distributions, a natural area measure is the integrated negativity

$$A_{\mathrm{nc}}(w)=\int_{P_\Omega(\alpha;w)<0}\!\! |P_\Omega(\alpha;w)|\,d^2\alpha,$$

although the source paper certifies nonclassicality directly by negativities of the filtered distribution rather than by an explicitly named “nonclassical area” (Kühn et al., 2014). In click-based homodyne schemes, analogous area-like measures are formed by integrating the negative part of sampled phase-space functions or the negative part of phase-resolved witness functions over the local-oscillator phase (Luis et al., 2014). The common object is therefore not a single invariant, but an integrated homodyne-accessible witness.

This plurality of definitions reflects the diversity of homodyne architectures. Standard balanced homodyne with linear photodiodes naturally yields optical tomograms and quadrature variances. Balanced or unbalanced homodyne with on-off arrays yields click moments and click-based phase-space functions. Generalized homodyne with nonclassical local oscillators yields non-Gaussian positive-operator-valued measures whose interference structure can itself be assigned an integrated nonclassical area (Combes et al., 2022).

2. Tomogram-based homodyne nonclassical area

In the tomographic formulation, the lower bound

$\sigma(\rho)\ge \sqrt{2}\,\pi$

follows from the uncertainty relation

$$w(X,\theta)=\langle X,\theta|\hat\rho|X,\theta\rangle,$$0

together with the symmetry

$$w(X,\theta)=\langle X,\theta|\hat\rho|X,\theta\rangle,$$1

Equality holds for pure coherent states, so for pure single-mode states $$w(X,\theta)=\langle X,\theta|\hat\rho|X,\theta\rangle,$$2 is a sufficient indicator of nonclassicality, whereas $$w(X,\theta)=\langle X,\theta|\hat\rho|X,\theta\rangle,$$3 characterizes displaced coherent states (Rohith et al., 2018).

For mixed states, the same quantity is not a strict witness, because classical mixing can increase quadrature variances. To address this, the mixed-state correction uses the Wigner–Yanase skew information

$$w(X,\theta)=\langle X,\theta|\hat\rho|X,\theta\rangle,$$4

and defines the “quantum” variance

$$w(X,\theta)=\langle X,\theta|\hat\rho|X,\theta\rangle,$$5

with corrected area

$$w(X,\theta)=\langle X,\theta|\hat\rho|X,\theta\rangle,$$6

For a single-mode thermal state of mean photon number $$w(X,\theta)=\langle X,\theta|\hat\rho|X,\theta\rangle,$$7, this construction yields $$w(X,\theta)=\langle X,\theta|\hat\rho|X,\theta\rangle,$$8 and hence $$w(X,\theta)=\langle X,\theta|\hat\rho|X,\theta\rangle,$$9, eliminating the spurious positive value produced by classical thermal broadening (Rohith et al., 2018).

The tomographic area is directly computable from homodyne samples. For each phase $$\langle \hat X_\theta^n\rangle=\int_{-\infty}^{\infty}dX\,X^n\,w(X,\theta),$$0, one acquires quadrature samples, estimates $$\langle \hat X_\theta^n\rangle=\int_{-\infty}^{\infty}dX\,X^n\,w(X,\theta),$$1, calibrates the vacuum shot-noise level so that coherent states satisfy $$\langle \hat X_\theta^n\rangle=\int_{-\infty}^{\infty}dX\,X^n\,w(X,\theta),$$2, and approximates the integral by numerical quadrature over $$\langle \hat X_\theta^n\rangle=\int_{-\infty}^{\infty}dX\,X^n\,w(X,\theta),$$3. Efficiency corrections use

$$\langle \hat X_\theta^n\rangle=\int_{-\infty}^{\infty}dX\,X^n\,w(X,\theta),$$4

or, with thermal contamination,

$$\langle \hat X_\theta^n\rangle=\int_{-\infty}^{\infty}dX\,X^n\,w(X,\theta),$$5

This makes the quantity experimentally tractable without density-matrix reconstruction or quasiprobability tomography (Rohith et al., 2018).

The same formalism extends to multimode systems. For $$\langle \hat X_\theta^n\rangle=\int_{-\infty}^{\infty}dX\,X^n\,w(X,\theta),$$6 modes, one defines

$$\langle \hat X_\theta^n\rangle=\int_{-\infty}^{\infty}dX\,X^n\,w(X,\theta),$$7

and

$$\langle \hat X_\theta^n\rangle=\int_{-\infty}^{\infty}dX\,X^n\,w(X,\theta),$$8

with coherent-state bound

$$\langle \hat X_\theta^n\rangle=\int_{-\infty}^{\infty}dX\,X^n\,w(X,\theta),$$9

The corresponding multimode homodyne nonclassical area is $$\mathrm{Var}(\hat X_\theta)=\langle \hat X_\theta^2\rangle-\langle \hat X_\theta\rangle^2,$$0 (Rohith et al., 2018).

3. State dependence, invariances, and dynamical behavior

The tomographic homodyne nonclassical area captures both isotropic and anisotropic departures from coherent-state behavior. For number states $$\mathrm{Var}(\hat X_\theta)=\langle \hat X_\theta^2\rangle-\langle \hat X_\theta\rangle^2,$$1,

$$\mathrm{Var}(\hat X_\theta)=\langle \hat X_\theta^2\rangle-\langle \hat X_\theta\rangle^2,$$2

independent of $$\mathrm{Var}(\hat X_\theta)=\langle \hat X_\theta^2\rangle-\langle \hat X_\theta\rangle^2,$$3, so

$$\mathrm{Var}(\hat X_\theta)=\langle \hat X_\theta^2\rangle-\langle \hat X_\theta\rangle^2,$$4

which is strictly positive for $$\mathrm{Var}(\hat X_\theta)=\langle \hat X_\theta^2\rangle-\langle \hat X_\theta\rangle^2,$$5 and increases monotonically with $$\mathrm{Var}(\hat X_\theta)=\langle \hat X_\theta^2\rangle-\langle \hat X_\theta\rangle^2,$$6 (Rohith et al., 2018).

For squeezed vacuum $$\mathrm{Var}(\hat X_\theta)=\langle \hat X_\theta^2\rangle-\langle \hat X_\theta\rangle^2,$$7 with $$\mathrm{Var}(\hat X_\theta)=\langle \hat X_\theta^2\rangle-\langle \hat X_\theta\rangle^2,$$8,

$$\mathrm{Var}(\hat X_\theta)=\langle \hat X_\theta^2\rangle-\langle \hat X_\theta\rangle^2,$$9

and

$$\Delta X_\theta=\sqrt{\mathrm{Var}(\hat X_\theta)}.$$0

where $$\Delta X_\theta=\sqrt{\mathrm{Var}(\hat X_\theta)}.$$1 and $$\Delta X_\theta=\sqrt{\mathrm{Var}(\hat X_\theta)}.$$2 is the complete elliptic integral of the second kind. The area is independent of the squeezing angle $$\Delta X_\theta=\sqrt{\mathrm{Var}(\hat X_\theta)}.$$3 and increases monotonically with $$\Delta X_\theta=\sqrt{\mathrm{Var}(\hat X_\theta)}.$$4 (Rohith et al., 2018). More generally, for squeezed Fock states $$\Delta X_\theta=\sqrt{\mathrm{Var}(\hat X_\theta)}.$$5 the variance factorizes and the area increases with both photon number and squeezing strength.

Displacement invariance is a persistent property. Coherent displacements change tomographic means but not variances, so the area is unchanged under phase-space translations. This is explicit for coherent states, squeezed coherent states, and displaced squeezed Fock states in the tomographic treatment (Rohith et al., 2018). Likewise, phase-space rotations merely shift the $$\Delta X_\theta=\sqrt{\mathrm{Var}(\hat X_\theta)}.$$6 origin and do not alter the integrated area.

A later tomogram-based study used the related quantity

$$\Delta X_\theta=\sqrt{\mathrm{Var}(\hat X_\theta)}.$$7

to track nonclassicality dynamics in Kerr and cubic media. There, $$\Delta X_\theta=\sqrt{\mathrm{Var}(\hat X_\theta)}.$$8 exhibits recurrences at the revival time, local minima at fractional revivals, monotonic decay to zero under amplitude damping, and saturation to a positive constant under phase damping because classical mixing broadens quadrature distributions (Athira et al., 4 May 2026). The same work emphasizes that for mixed states this area is not a strict witness: phase damping can leave a positive asymptotic value even when the surviving structure is classical in origin.

4. Phase-space negativity areas from balanced homodyne reconstruction

A distinct lineage associates homodyne nonclassical area with negativity regions of a regularized phase-space distribution. In the nonclassicality-filter approach, one starts from the normally ordered characteristic function $$\Delta X_\theta=\sqrt{\mathrm{Var}(\hat X_\theta)}.$$9 and defines a filtered quasiprobability

$$\sigma(\rho):=\int_0^{2\pi}d\theta\,\Delta X_\theta,$$0

Negativities of $$\sigma(\rho):=\int_0^{2\pi}d\theta\,\Delta X_\theta,$$1 certify nonclassicality. The filters are chosen so that $$\sigma(\rho):=\int_0^{2\pi}d\theta\,\Delta X_\theta,$$2 is regular, its negativities arise only from the state, and in the invertible case the full state information is preserved (Kühn et al., 2014).

For balanced homodyne detection, the filtered distribution can be sampled directly from quadrature-phase pairs $$\sigma(\rho):=\int_0^{2\pi}d\theta\,\Delta X_\theta,$$3 via

$$\sigma(\rho):=\int_0^{2\pi}d\theta\,\Delta X_\theta,$$4

with pattern function

$$\sigma(\rho):=\int_0^{2\pi}d\theta\,\Delta X_\theta,$$5

and

$$\sigma(\rho):=\int_0^{2\pi}d\theta\,\Delta X_\theta,$$6

The variance estimate is

$$\sigma(\rho):=\int_0^{2\pi}d\theta\,\Delta X_\theta,$$7

Statistical significance is then quantified by

$$\sigma(\rho):=\int_0^{2\pi}d\theta\,\Delta X_\theta,$$8

A target such as $$\sigma(\rho):=\int_0^{2\pi}d\theta\,\Delta X_\theta,$$9 determines the sample size needed for high-significance certification (Kühn et al., 2014).

Within this framework, a natural negativity-based area is

$\sqrt{2}\pi$0

Numerically, one reconstructs $\sqrt{2}\pi$1 on an $\sqrt{2}\pi$2 grid, identifies negative cells, and approximates the integral by a weighted sum. Error bars can be estimated either by quadrature of the pointwise variances or by bootstrap resampling of the homodyne data (Kühn et al., 2014).

This construction extends to multimode radiation fields. For two modes,

$\sqrt{2}\pi$3

and the corresponding negative-volume measure becomes

$\sqrt{2}\pi$4

In the reported examples, negativities were visualized for phase-randomized squeezed vacuum, single-photon-added thermal states, and two-mode heralded states, with the band-limited filter $\sqrt{2}\pi$5 minimizing the required data volume in many regimes (Kühn et al., 2014).

5. Click-counting, phase-resolved witnesses, and area constructions

Homodyne nonclassical area also appears in click-based homodyne detection. In unbalanced homodyne with an array of on-off detectors, the displaced signal is measured through click statistics

$\sqrt{2}\pi$6

with

$\sqrt{2}\pi$7

The directly sampled click-based phase-space function is

$\sqrt{2}\pi$8

For even $\sqrt{2}\pi$9, any negativity of $\mathrm{NA}\equiv \sigma(\rho)-\sqrt{2}\pi$0 with $\mathrm{NA}\equiv \sigma(\rho)-\sqrt{2}\pi$1 certifies nonclassicality (Luis et al., 2014).

A corresponding area is the integrated negativity

$\mathrm{NA}\equiv \sigma(\rho)-\sqrt{2}\pi$2

estimated on a displacement grid by

$\mathrm{NA}\equiv \sigma(\rho)-\sqrt{2}\pi$3

This construction is especially useful in the weak-light regime, where small even detector numbers can reveal negativities efficiently. The same formalism underlies later experimental work showing significant phase-space negativities in a high-loss regime with click-counting detectors, where a related quantity was written as

$\mathrm{NA}\equiv \sigma(\rho)-\sqrt{2}\pi$4

for the click-based function

$\mathrm{NA}\equiv \sigma(\rho)-\sqrt{2}\pi$5

with pointwise statistical error

$\mathrm{NA}\equiv \sigma(\rho)-\sqrt{2}\pi$6

In that experiment, clear negativities were observed at $\mathrm{NA}\equiv \sigma(\rho)-\sqrt{2}\pi$7, and for $\mathrm{NA}\equiv \sigma(\rho)-\sqrt{2}\pi$8 with $\mathrm{NA}\equiv \sigma(\rho)-\sqrt{2}\pi$9 the maximal signed significance reached $$A_{\mathrm{nc}}(w)=\int_{P_\Omega(\alpha;w)<0}\!\! |P_\Omega(\alpha;w)|\,d^2\alpha,$$0 at $$A_{\mathrm{nc}}(w)=\int_{P_\Omega(\alpha;w)<0}\!\! |P_\Omega(\alpha;w)|\,d^2\alpha,$$1 (Bohmann et al., 2017).

A different click-homodyne route defines phase-domain areas from witness functions. Balanced homodyne with on-off detector arrays introduces a nonlinear quadrature operator

$$A_{\mathrm{nc}}(w)=\int_{P_\Omega(\alpha;w)<0}\!\! |P_\Omega(\alpha;w)|\,d^2\alpha,$$2

with nonclassicality certified by

$$A_{\mathrm{nc}}(w)=\int_{P_\Omega(\alpha;w)<0}\!\! |P_\Omega(\alpha;w)|\,d^2\alpha,$$3

or, more generally, by negative principal minors of the click-moment matrix. An area-like measure is then

$$A_{\mathrm{nc}}(w)=\int_{P_\Omega(\alpha;w)<0}\!\! |P_\Omega(\alpha;w)|\,d^2\alpha,$$4

where $$A_{\mathrm{nc}}(w)=\int_{P_\Omega(\alpha;w)<0}\!\! |P_\Omega(\alpha;w)|\,d^2\alpha,$$5 may be the normally ordered variance or a higher-order determinant. This turns a phase-resolved witness into an integrated measure of the strength and phase extent of nonclassicality (Sperling et al., 2014).

6. Generalized homodyne measurements, multimode bright beams, and operational interpretations

Not all homodyne nonclassical areas are attached to reconstructed states; some are attached to the measurement itself. When the local oscillator is a Schrödinger cat state rather than a coherent state, balanced homodyne no longer implements a standard quadrature projector. The resulting POVM density contains interference terms

$$A_{\mathrm{nc}}(w)=\int_{P_\Omega(\alpha;w)<0}\!\! |P_\Omega(\alpha;w)|\,d^2\alpha,$$6

with $$A_{\mathrm{nc}}(w)=\int_{P_\Omega(\alpha;w)<0}\!\! |P_\Omega(\alpha;w)|\,d^2\alpha,$$7, and after integrating over the summed photocount variable $$A_{\mathrm{nc}}(w)=\int_{P_\Omega(\alpha;w)<0}\!\! |P_\Omega(\alpha;w)|\,d^2\alpha,$$8 the measurement reduces to the reflection-symmetric projector

$$A_{\mathrm{nc}}(w)=\int_{P_\Omega(\alpha;w)<0}\!\! |P_\Omega(\alpha;w)|\,d^2\alpha,$$9

In this setting a nonclassical area can be defined operationally as the integrated excess probability density generated by the interference terms,

$\sigma(\rho)\ge \sqrt{2}\,\pi$0

Retaining $\sigma(\rho)\ge \sqrt{2}\,\pi$1, or at least the parity of $\sigma(\rho)\ge \sqrt{2}\,\pi$2, preserves the interference and maximizes this area; discarding $\sigma(\rho)\ge \sqrt{2}\,\pi$3 drives the interference envelope to zero in the large-$\sigma(\rho)\ge \sqrt{2}\,\pi$4 limit (Combes et al., 2022).

A more engineering-oriented interpretation appears in homodyne detection of bright multimode nonclassical states. There the measured normalized variance is

$\sigma(\rho)\ge \sqrt{2}\,\pi$5

with

$\sigma(\rho)\ge \sqrt{2}\,\pi$6

where unmatched bright modes contribute excess noise. The “nonclassical area” is then the set of modes or frequencies for which the measured quadrature noise remains below shot noise. In that picture the area shrinks as $\sigma(\rho)\ge \sqrt{2}\,\pi$7 grows at fixed mode mismatch, and squeezing is observable only if

$\sigma(\rho)\ge \sqrt{2}\,\pi$8

For Gaussian entanglement, a practical rule extracted in that analysis is $\sigma(\rho)\ge \sqrt{2}\,\pi$9 if entanglement is to remain detectable at macroscopic brightness (Usenko et al., 2015).

This suggests a useful taxonomy. One may distinguish a phase-space area, defined through negativity regions of a reconstructed quasiprobability; a phase-domain area, defined by integrating a phase-resolved witness; and a tomographic area, defined directly from the $$w(X,\theta)=\langle X,\theta|\hat\rho|X,\theta\rangle,$$00-dependent quadrature standard deviation. The literature uses all three, depending on the measurement model and the aspect of nonclassicality under study.