---
title: State-Subtraction Methods
url: https://www.emergentmind.com/topics/state-subtraction-method
type: topic
---

# State-Subtraction Methods

State-subtraction method denotes a family of subtraction-based constructions whose meaning depends on disciplinary context. In microwave and continuous-variable quantum optics, it refers to conditional photon-subtraction operations used for non-Gaussian state engineering; in SU(1,1) interferometry, to internal multi-photon subtraction of the probe state; in perturbative QCD, to local subtraction of unresolved radiation or infrared counterterms associated with specific phase-space configurations; and in nonequilibrium kinetic theory, to real intermediate state subtraction that removes on-shell resonant contributions already represented by decay and inverse-decay channels [1210.4042, 2311.14612, 1806.09570, 2309.16615]. This suggests that the expression is field-dependent, but the common structural motif is the controlled removal of a designated state component in order to expose a finite, non-Gaussian, or non-double-counted remainder.

## 1. Quantum-optical operator content

In the microwave setting of dichotomic single-photon detection, the central object is the subtraction operator
\[
\hat{B}\equiv\sum_{n=1}^\infty |n-1\rangle\langle n|,
\]
which acts as \(\hat{B}|n\rangle=|n-1\rangle\) for \(n\ge 1\) and \(\hat{B}|0\rangle=0\). The corresponding click-conditioned state update is
\[
\rho \longrightarrow \rho'=\frac{\hat{B}\rho\hat{B}^\dagger}{\mathrm{Tr}[\hat{B}\rho\hat{B}^\dagger]}.
\]
The same paper emphasizes that \(\hat{B}\) is not the annihilation operator \(\hat a\), but is related to it by \(\hat a=\hat B\sqrt{\hat n}\); the absence of the \(\sqrt n\) weighting is the essential nonlinearity of the measurement back-action in the Josephson Photomultiplier regime [1210.4042].

For bipartite continuous-variable systems, coherent subtraction generalizes fixed-order subtraction to coherent superpositions of many subtraction patterns. The subtraction operator is written as
\[
C_S=\sum_{m_1,m_2=0}^{\infty} A_{m_1m_2}\,a_1^{m_1}a_2^{m_2},
\qquad
\rho_S=C_S\rho_G C_S^\dagger,
\]
with \(\rho_G\) a Gaussian kernel. Jiang et al. introduce the still more general operator
\[
C_{SA}=\sum_{m_1,m_2,m_3,m_4=0}^{\infty}
A_{m_1m_2m_3m_4}\,
a_1^{\dagger m_1}a_2^{\dagger m_2}a_1^{m_3}a_2^{m_4},
\]
and note that any two-mode CV state can be written in this form [1211.5826].

The resulting distinction is technically important. Single-mode subtraction, fixed local subtraction \(a_1^m a_2^{m'}\), coherent subtraction \(\sum A_{m_1m_2}a_1^{m_1}a_2^{m_2}\), and annihilation-based subtraction \(a\) are different operations with different number-basis nonlinearities, separability properties, and implementation assumptions. In the cited literature, “state-subtraction method” therefore does not mean a unique operator identity, but a class of conditional or coherent subtraction maps.

## 2. Measurement-based state engineering in microwave cavities

The measurement-based version starts from a coherent state
\[
|\alpha\rangle=e^{-|\alpha|^2/2}\sum_{n=0}^{\infty}\frac{\alpha^n}{\sqrt{n!}}|n\rangle,
\]
and applies one or more click-conditioned subtraction events. After one subtraction,
\[
\hat B|\alpha\rangle
=
e^{-|\alpha|^2/2}\,\alpha\,
\sum_{m=0}^{\infty}\frac{\alpha^m}{\sqrt{(m+1)!}}|m\rangle,
\]
which is not a coherent state; the photon-number distribution is no longer Poissonian, and the resulting state is non-Gaussian and nonclassical. After \(N\) clicks, the unnormalized state is \(\hat B^N|\alpha\rangle\), and the success probability is
\[
P_N
=
1-e^{-|\alpha|^2}\sum_{n=0}^{N-1}\frac{|\alpha|^{2n}}{n!}
=
1-\frac{\Gamma(N,|\alpha|^2)}{\Gamma(N)}.
\]
The protocol then exploits the noncommutativity of subtraction and coherent displacement,
\[
\hat B\hat D(\gamma)\neq \hat D(\gamma)\hat B,
\]
to engineer families of target states by sequences of displacements and \(N\)-photon subtraction steps [1210.4042].

The paper identifies three principal output families: ordinary squeezed vacuum states (\(k=2\)), generalized squeezed states with \(k\)-fold rotational symmetry (\(k\ge 3\)), and squeezed multi-component Schrödinger cat states. Their phase-space diagnostics are the Wigner and Husimi–\(Q\) functions; the cited simulations show \(k\)-fold rotational symmetry, strong squeezing along specific quadratures, multiple coherent peaks for cat-like states, and interference fringes with negativity. The same work reports generalized squeezed states (\(k=2,3,4\)) with fidelities \(>99\%\) with respect to
\[
|\Psi_k\rangle=S^{(k)}(z)|0\rangle,
\qquad
S^{(k)}(z)=\exp\!\Big[-\tfrac12\left(z(\hat a^\dagger)^k-z^*\hat a^k\right)\Big],
\]
and success probabilities \(P\) between \(0.995\) and \(1.0\) for modest \(|\alpha|\sim 4\text{--}7\). For squeezed multi-component cat states, the paper gives fidelities \(\gtrsim 94\%\) and success probabilities typically between \(\sim 0.25\) and \(\sim 10^{-3}\text{--}10^{-2}\) in the reported examples [1210.4042].

Implementation is discussed in the setting of a high-\(Q\) superconducting microwave resonator strongly coupled to a Josephson Photomultiplier. In the operating regime of interest, the click Kraus operator is effectively \(\hat B\), while relaxation and dark counts are treated as deviations from the ideal model. Superconducting resonators with \(Q\)-factors approaching \(10^7\) and photon lifetimes \(\sim 10^5\) ns, together with JPM measurement times \(O(10\text{--}100)\) ns or less, imply that on the order of \(10^2\text{--}10^3\) detection operations can be performed within a single cavity lifetime [1210.4042].

## 3. Coherent subtraction and entanglement in bipartite continuous variables

For bipartite Gaussian kernels, coherent subtraction is analyzed primarily as an entanglement-engineering and entanglement-detection procedure. If \(\rho_G^{II}\) is a separable Gaussian kernel in second standard form, then its coherent-state decomposition yields
\[
\rho_G^{II}
=
\int P(\alpha_1,\alpha_2)\,
|\alpha_1,\alpha_2\rangle\langle\alpha_1,\alpha_2|\,
d^2\alpha_1\,d^2\alpha_2,
\]
with \(P(\alpha_1,\alpha_2)\) Gaussian and positive-definite. Acting with coherent subtraction gives
\[
\rho_S
=
\int P(\alpha_1,\alpha_2)\,
\left|
\sum_{m,m'}A_{mm'}\alpha_1^m\alpha_2^{m'}
\right|^2
|\alpha_1,\alpha_2\rangle\langle\alpha_1,\alpha_2|\,
d^2\alpha_1\,d^2\alpha_2,
\]
which remains separable. Accordingly, coherent subtraction cannot create entanglement from a separable second-standard-form Gaussian kernel [1211.5826].

The symmetric coherent subtraction of a symmetric two-mode squeezed thermal state,
\[
\rho_{SSst}
=
\sum_{m=0}^{\infty}A_m a_1^m a_2^m\,\rho_{st}\,
\sum_{n=0}^{\infty}A_n^* a_1^{\dagger n}a_2^{\dagger n},
\]
has a sharper result. For the Gaussian kernel, the necessary and sufficient separability condition is
\[
2b_0+2c_1\ge 1.
\]
Jiang et al. show that the separability condition of \(\rho_{SSst}\) is exactly the same as that of its Gaussian kernel. In this case, symmetric coherent subtraction preserves the separability boundary rather than shifting it [1211.5826].

A different behavior occurs for a generic first-standard-form Gaussian kernel at the boundary \(\tau=1\), with
\[
\tau=\frac{1}{4(b_0+c_1)^2-c_2^2}.
\]
For
\[
\rho_{SS}
=
\sum_m A_m a_1^m a_2^m\,\rho_G\,
\sum_n A_n^* a_1^{\dagger n}a_2^{\dagger n},
\]
the realignment criterion leads to a sufficient entanglement condition in terms of
\[
B_{mn}=m!n!(A_mA_n^*+A_nA_m^*).
\]
If \(B_{mn}\ge 0\) for all \(m,n\) and \(B_{mn}>0\) for at least one pair \(m\neq n\), then the symmetric coherent subtraction of the boundary-separable Gaussian kernel is entangled. In phase language, constructive interference among subtraction amplitudes is the relevant design rule [1211.5826].

This body of results places coherent subtraction between two limits. In one limit it is a non-entangling dressing of a separable Gaussian mixture; in the other it is an entanglement-generating non-Gaussian operation whose efficacy depends on the coefficient structure of the coherent superposition.

## 4. Internal multi-photon subtraction in SU(1,1) interferometry

In SU(1,1) interferometry, the subtraction step is placed inside the interferometer, after the first two-mode parametric amplifier and before the phase shifter and second amplifier. With input
\[
|\Psi_{\mathrm{in}}\rangle=|\alpha\rangle_a\otimes|0\rangle_b,
\]
the subtraction operator is
\[
U_P=a^m\otimes b^n,
\]
so that the normalized post-selected state is
\[
|\psi_{m,n}\rangle
=
\frac{a^m b^n|\psi\rangle}
{\sqrt{\langle\psi|(a^\dagger)^m a^m (b^\dagger)^n b^n|\psi\rangle}}.
\]
The phase is encoded only on mode \(a\) through \(U_\phi=e^{i\phi a^\dagger a}\), and homodyne detection measures
\[
X=\frac{a+a^\dagger}{\sqrt2}.
\]
Phase sensitivity is evaluated by
\[
\Delta\phi
=
\frac{\sqrt{\langle X^2\rangle-\langle X\rangle^2}}
{\left|\frac{\partial}{\partial\phi}\langle X\rangle\right|}.
\]
The paper analyzes phase sensitivity, quantum Fisher information, and quantum Cramér–Rao bound for this internally photon-subtracted probe [2311.14612].

The principal findings are systematic. Internal photon subtraction always improves phase sensitivity relative to the standard SU(1,1) interferometer without subtraction, and increasing the number of subtracted photons improves sensitivity further. The QFI,
\[
F=4\langle\Delta^2 n_a\rangle,
\]
also increases with gain \(g\), coherent amplitude \(\alpha\), and subtraction number. For fixed \(N=4\), the standard SU(1,1) interferometer cannot surpass SQL, whereas symmetric multi-photon subtraction on both modes can surpass SQL over a broad parameter range, even under significant internal losses such as \(T_k=0.5\) [2311.14612].

Mode asymmetry is intrinsic because the phase shift and homodyne measurement both act on mode \(a\). As a result, single-mode subtraction on \(a\) and on \(b\) has different metrological consequences. For homodyne-based phase sensitivity, subtraction on mode \(b\) tends to be better at small \(\phi\) or small \(g\), while subtraction on mode \(a\) becomes superior at larger \(\phi\) or larger \(g\). For QFI and QCRB, subtraction on mode \(b\) is generally slightly better overall. The same analysis under loss yields
\[
F_L
=
\frac{4F\eta\langle n_a\rangle}{(1-\eta)F+4\eta\langle n_a\rangle},
\]
and the photon-subtracted schemes retain improved robustness against internal photon losses [2311.14612].

In this metrological usage, state subtraction is not primarily a state-generation primitive for standalone nonclassicality, but a controlled non-Gaussian transformation of the intracavity probe state.

## 5. Infrared subtraction in perturbative QCD

In perturbative QCD, subtraction methods are introduced because real and virtual corrections are separately infrared divergent. The subtraction term locally reproduces the soft and collinear behavior of the real matrix element in phase space, is simple enough to be integrated analytically over unresolved degrees of freedom, and is then added back in integrated form so that pole cancellation with virtual terms becomes explicit. The cited literature explicitly interprets this logic as a state-subtraction method in the sense that one subtracts local counterterms associated with a given color-ordered or sector-defined phase-space configuration [1001.2397, 1806.09570].

In NNLO antenna subtraction with one hadronic initial state, the required ingredients are initial–final three-parton tree antennae \(X^0_{i,jk}\), four-parton tree antennae \(X^0_{i,jkl}\), and one-loop three-parton antennae \(X^1_{i,jk}\). The framework organizes double-real, real–virtual, and double-virtual pieces by subtracting antenna counterterms built from color-connected radiators and reduced matrix elements with mapped momenta. For the integration of initial–final antennae, the phase-space integrals are rewritten as cut loop integrals and reduced by IBP and Lorentz-invariance identities to master integrals; the paper reports \(9\) master integrals for NNLO double-real initial–final antennae and \(6\) master integrals for NNLO one-loop initial–final antennae. The integrated antennae are distributions in the Bjorken-like variable \(z\), involving \(\delta(1-z)\) and plus-distributions, and were cross-checked against known NNLO DIS coefficient functions [1001.2397].

Local analytic sector subtraction reorganizes the same infrared problem with sector functions and Catani–Seymour-type mappings. At NNLO for final-state radiation, the method is defined as local, sectorized, and analytic: counterterms match singular behavior point-by-point in phase space, singular regions are isolated by sector partitions, and the counterterms are integrated analytically over unresolved radiation. The construction uses NLO sector functions \(W_{ij}\), NNLO sector functions \(W_{ijkl}\), projectors such as \(\mathbf S_i\), \(\mathbf C_{ij}\), \(\mathbf S_{ij}\), \(\mathbf C_{ijk}\), and \(\mathbf C_{ijkl}\), and produces finite combinations of the schematic form
\[
RR-\mathbf 1_{ij}RR-\mathbf 2_{ijkl}RR+\mathbf 1_{ij}\mathbf 2_{ijkl}RR.
\]
The corresponding integrated counterterms cancel the poles of real–virtual and double-virtual contributions [1806.09570, 1912.09368].

The NLO extension to initial- and final-state radiation preserves the same philosophy. The real phase space is partitioned by positive sector functions \(W_{ij}\) satisfying \(\sum_i\sum_{j\neq i}W_{ij}=1\), with each sector containing at most one soft leg and one collinear pair. The local counterterm is
\[
K=\sum_i\sum_{j\neq i}\big[{\cal S}_i+{\cal C}_{ij}-{\cal S}_i{\cal C}_{ij}\big]R\,W_{ij},
\]
and its analytical integral yields \(I\) and \(J\) terms that cancel, respectively, the poles of virtual corrections and PDF counterterms. The same work introduces smooth damping factors with tunable exponents \(\alpha,\beta,\gamma\) to improve numerical stability without changing the pole structure, and numerically validates the method for processes including \(e^+e^-\to jj\), \(e^+e^-\to jjj\), \(pp\to Z\), \(pp\to Zj\), and \(pp\to W^+W^-j\) [2209.09123].

Within QCD, therefore, state subtraction is neither measurement-based nor heralded. It is an analytic device for removing unresolved radiation contributions at the integrand level and restoring them in a form where infrared cancellation is manifest.

## 6. Real intermediate state subtraction in kinetic theory

In nonequilibrium kinetic theory, the subtraction target is an on-shell resonance already included as an explicit species in the Boltzmann network. For a process
\[
aa\to b\to cc
\]
with unstable \(b\), the propagator is
\[
D_b(s)=\frac{1}{s-m_b^2+i\,m_b\Gamma_b}.
\]
If the kinetic equations also contain \(b\leftrightarrow aa\) and \(b\leftrightarrow cc\), then the resonant contribution to \(aa\leftrightarrow cc\) is double-counted unless the real intermediate state is subtracted from the \(2\to2\) channel [2309.16615].

The cited analysis traces the ambiguity of RIS subtraction to the on-shell approximation that reduces the Schwinger–Dyson and Kadanoff–Baym description to the Boltzmann limit. Standard RIS subtraction and principal-value subtraction share the same on-shell limit,
\[
D^{\rm on}(a)=-i\pi\delta(a-m^2),
\qquad
|D^{\rm on}(a)|^2=\frac{\pi}{m\Gamma}\delta(a-m^2),
\]
but differ in the off-shell Hermitian part of the propagator. The paper argues that both SRS and PVS can yield negative effective scattering rates because they overestimate the on-shell weight and force the off-shell remainder to become negative in parts of the resonance neighborhood [2309.16615].

To avoid this, the paper introduces a cut-subtraction scheme. With a cut function \(\Theta(a-m^2,\Delta)\), the off-shell propagator is defined as
\[
D^\Delta(a)=\big[1-\Theta(a-m^2,\Delta)\big]D(a),
\]
and the on-shell weight is measured by
\[
R(\Delta)
=
\frac{1}{\pi}\int ds\,\Theta(s-m^2,\Delta)\,m\Gamma\,|D(s)|^2.
\]
Here \(R(\Delta)\) is an effective one-particle weight function: \(R(\Delta)\approx 1\) indicates that the near-resonant region is well approximated by an on-shell particle description, whereas \(R(\Delta)\ll 1\) signals that off-shell contributions dominate and the Boltzmann treatment of the resonance as an independent thermal species is not reliable [2309.16615].

This version of state subtraction has a different objective from both quantum-optical subtraction and QCD infrared subtraction. It is a consistency prescription for kinetic equations, ensuring that decay/inverse-decay terms and resonant \(2\to2\) channels do not represent the same physical intermediate state twice.

Source: https://www.emergentmind.com/topics/state-subtraction-method