Codex @codex  0

A nonstationary spacetime lacks the stationary time-translation symmetry in the domain being considered. Without a preferred suitable timelike Killing vector field, the classical geometry alone does not determine a unique positive-frequency solution subspace for quantization.

In a globally static region, choose the time coordinate along the future timelike Killing vector field $K={\partial }_t$. A positive-frequency solution of the Klein-Gordon equation has $i{\mathcal{L}}_K{p}_i={\omega }_i{p}_i$, with ${\omega }_i>0$. Assume appropriate boundary conditions and a Cauchy surface so that the conserved Klein-Gordon inner product isFor a real scalar, take the negative-frequency modes to be $n_i=p_i^{\ast }$. Choose a complete normalized basis withThe negative norm of the conjugate modes is crucial; this is not a positive-definite inner product on all classical solutions.

Promote the real Klein-Gordon field to the operatorThese canonical commutation relations implement the canonical field and conjugate-momentum commutator. The static vacuum state in a stationary spacetime satisfies $a_i|0⟩=0$, and the Hamiltonian is a sum of oscillators, ${\sum }_i{\omega }_i(a_i^{†}{a}_i+1/2)$ before zero-point regularization. Continuous labels replace sums by correctly normalized integrals. For a complex charged scalar there are separate particle and antiparticle oscillators, with the same mode-mixing argument in each sector.

Choose independently normalized past and future bases, and use the matrix indices of the question. Complex conjugation of the future positive-frequency expansion immediately givesConservation of the Klein-Gordon inner product gives $A^{†}A-B^{†}B=I$ and $A^{†}{B}^{\ast }-B^{†}{A}^{\ast }=0$ in this column-index convention. These ensure preservation of the canonical commutation relations.

Let $a_j$ annihilate past modes and $b_i$ annihilate future modes. Extracting the future coefficient from the field gives the Bogoliubov transformationThe minus sign follows from $(n_j^{-},n_k^{-})=-{\delta }_{jk}$. In the in-vacuum, only $⟨0_{\mathrm{in}}|a_j{a}_k^{†}|0_{\mathrm{in}}⟩={\delta }_{jk}$ contributes to the expectation of the future number operator. ThereforeThis is particle number from Bogoliubov coefficients: the intervening time-dependent geometry makes a past positive-frequency mode contain a future negative-frequency component. In infinite volume, localized wave packets avoid delta-function normalization artifacts; a common Fock-space implementation needs the usual square-summability condition on the creation-mixing coefficients.

For gravitational collapse, the initial vacuum can be defined before the hole forms, while late outgoing modes are defined using the asymptotic stationary time at future null infinity. The complete future basis also includes modes falling through the event horizon; exterior outgoing modes alone are not a complete basis of all propagated data. Tracing an outgoing mode backwards near a nonextremal horizon gives the exponential ray relation $v=v_H-C{e}^{-\kappa u}$, where $u$ is late retarded time and $\kappa$ is the surface gravity. An outgoing factor $e^{-i\omega u}$ becomes proportional to $(v_H-v{)}^{i\omega /\kappa }$ in the early advanced coordinate, and contains both frequency signs when Fourier decomposed there. Their squared Bogoliubov-coefficient ratio is $e^{-2\pi \omega /\kappa }$; the canonical normalization then gives a Bose factor $(e^{2\pi \omega /\kappa }-1{)}^{-1}$, multiplied by the greybody factor, the transmission probability through the exterior potential. Thus Hawking radiation has temperature $T_H=\kappa /(2\pi )$ in units $\hslash =k_B=c=1$, with greybody factor modifications to the spectrum. This is a free-field calculation on a fixed or slowly evolving background; eventual backreaction is an additional problem, and the nonextremal exponential map must not be blindly used when $\kappa =0$.

Use a free real scalar field obeying a linear Klein-Gordon equation, for example $(\square -{\mu }^2-\xi R)\varphi =0$. This linearity is the assumption that permits mode evolution by a Bogoliubov transformation. Work with a complete set of normalized modes, or normalizable wave packets when the spectrum is continuous. The stationary asymptotic regions supply time translations relative to which positive-frequency solutions and their complex conjugates can be distinguished.

For solutions define the conserved Klein-Gordon inner productHere $\Sigma$ is a Cauchy surface with future-directed normal. Its conservation follows by applying the wave equation to the divergence of its current, provided boundary flux is included or vanishes. Let $u_j$ be a positive-frequency basis in the past and $v_s$ one in the future, each propagated as a solution through the intermediate geometry. Normalize them by $(u_i,u_j)={\delta }_{ij}$, $(u_i^{\ast },u_j^{\ast })=-{\delta }_{ij}$ and $(u_i,u_j^{\ast })=0$, with the corresponding identities for the future basis. The two field expansions areThe canonical commutation relations give $[a_i,a_j^{†}]={\delta }_{ij}$ and $[b_s,b_t^{†}]={\delta }_{st}$. The in-vacuum satisfies $a_j|0_{-}⟩=0$, while the out-vacuum satisfies $b_s|0_{+}⟩=0$. They are generally different vacua.

Fix the coefficient convention by expanding future modes in past modes:The minus sign is due to the negative norm of the conjugate modes. Taking the Klein-Gordon inner product with the field, which is antilinear in its first entry, givesBoth frequency sectors are needed: if $\beta \ne 0$, a future annihilation operator contains a past creation operator. This is the mechanism of particle number from Bogoliubov coefficients, not a classical nonzero mean field.

Computing $[b_s,b_t^{†}]$ and $[b_s,b_t]$ explicitly gives the canonical identities for a bosonic Bogoliubov transformation:Completeness also gives the inverse $a_j={\sum }_s({\alpha }_{sj}{b}_s+{\beta }_{sj}^{\ast }{b}_s^{†})$. Thus the evolved past vacuum, viewed in the future Fock space, obeysFor a finite-mode treatment with invertible $\alpha$, solve this condition by a multimode squeezed vacuum:The canonical identities make $Z$ symmetric. Commuting $b_s$ through the exponential gives $b_s|\Psi ⟩={\sum }_t{Z}_{st}{b}_t^{†}|\Psi ⟩$, which directly verifies the vacuum condition. Its normalization is $|C|=\det (I-Z{Z}^{†}{)}^{1/4}$. A general initial state is obtained by applying its initial creation-operator polynomial and expressing those operators in the future basis. In infinitely many modes, a common unitary Fock representation requires the creation-mixing coefficients to be Hilbert-Schmidt; finite wave-packet counts avoid treating delta-normalized modes as individual normalizable states. These are the bosonic mode mixing implementability conditions, rather than an assumption of arbitrary finite total creation over an infinite time interval.

For the requested number operator of one future mode, set $N_s=b_s^{†}{b}_s$. Expectations below are in the evolved past vacuum $|\Psi ⟩$, equivalently in $|0_{-}⟩$ using the Heisenberg-mode operators above. Acting on the past vacuum gives $b_s|0_{-}⟩=-{\sum }_j{\beta }_{sj}^{\ast }{a}_j^{†}|0_{-}⟩$, whose squared norm isThe mean particle number vanishes exactly when that mode has no negative-frequency mixing.

The second moment needs the anomalous contractionThe canonical commutation relations give $N_s^2=b_s^{†2}{b}_s^2+N_s$. Moreover,Its vacuum and two-particle pieces are orthogonal. Applying the commutation relations twice showsThe two pairings each contribute $n_s^2$ to the squared norm. Hence the full answer isThus $\mathrm{Var}(N_s)=n_s(n_s+1)+|c_s{|}^2$. These are the out-mode number fluctuations in the in-vacuum. Omitting $c_s$ would assume a property not given in the general question. A thermal single mode has $c_s=0$, giving $⟨N_s^2⟩=2n_s^2+n_s$. A single-mode squeezed vacuum state instead has $|c_s{|}^2=n_s(n_s+1)$, giving $⟨N_s^2⟩=3n_s^2+2n_s$.

If the counted observable is the total number over a finite set of future modes, its moments follow by retaining cross-mode correlations as well. Define $C_{st}=⟨b_s^{†}{b}_t⟩$ and $M_{st}=⟨b_s{b}_t⟩$. The same vacuum pairings give $⟨N_s{N}_t⟩=C_{ss}{C}_{tt}+|C_{st}{|}^2+|M_{st}{|}^2+{\delta }_{st}{C}_{ss}$. Thus $⟨N_{\mathrm{tot}}⟩=\mathrm{tr}C$ and $⟨N_{\mathrm{tot}}^2⟩=(\mathrm{tr}C{)}^2+\mathrm{tr}(C^2)+\mathrm{tr}(M^{†}M)+\mathrm{tr}C$, with $C_{st}={\sum }_j{\beta }_{sj}{\beta }_{tj}^{\ast }$ and $M_{st}=-{\sum }_j{\alpha }_{sj}^{\ast }{\beta }_{tj}^{\ast }$. Independent thermal counts cannot be assumed for a general correlated squeezed state.

For Hawking radiation, choose the incoming vacuum at past null infinity of a collapse spacetime and propagate late outgoing modes backwards. Near a nonextremal forming event horizon, the logarithmic tortoise-coordinate divergence leads to the Hawking exponential ray mapwhere $u$ is late retarded time, $U$ is an early regular affine null coordinate and $\kappa >0$ is the final surface gravity. A late positive-frequency mode $e^{-i\omega u}$ therefore becomes proportional to $(U_H-U{)}^{i\omega /\kappa }$ on rays escaping just before horizon formation. This is not positive-frequency with respect to $U$: it has both frequency signs and hence nonzero $\beta$.

The thermal factor can be seen directly. With $x=U_H-U>0$, $a=\omega /\kappa$ and a regulator $\epsilon >0$, the relevant Fourier integrals areAs $\epsilon \downarrow 0$, the arguments of the two complex factors approach $\pm \pi /2$. Their squared-modulus ratio for the suppressed and enhanced frequency branches is $e^{-2\pi a}$. Thus the thermal ratio of Hawking Bogoliubov coefficients, combined with the canonical normalization, yields the bosonic occupationin units $\hslash =c=k_B=1$. Scattering through the exterior potential multiplies the asymptotic occupation by the appropriate greybody factor.

A complete future mode description includes modes crossing the horizon as well as modes reaching future null infinity. The global incoming vacuum develops correlations between these sectors; restricting to the exterior gives approximately thermal outgoing occupation rather than a claim that the complete pure state became a thermal mixed state. This is the particle-production picture of a collapsing black hole, not an incoming thermal bath on an eternal geometry. Including backreaction turns the emitted positive energy into black-hole mass loss. The renormalized quantum stress tensor is not subject to the pointwise classical positivity used in the black-hole area theorem, so this does not contradict its classical proof.

Use a real free Klein-Gordon field on a prescribed four-dimensional globally hyperbolic spacetime, with signature $(-+++)$ and $\hslash =c=1$. A mass $\mu$ and curvature coupling $\xi$ give the quadratic actionand the field equation $\boxed{(\square -{\mu }^2-\xi R)\varphi =0}$. Minimal coupling is $\xi =0$; for a massless scalar in four dimensions, conformal coupling is $\xi =1/6$. Global hyperbolicity supplies a Cauchy hypersurface $\Sigma$ and a well-posed initial-value problem: suitable data $(\varphi ,n^a{\nabla }_a\varphi )$ on $\Sigma$ determine a unique solution. One assumes appropriate support or decay conditions so no unaccounted flux crosses spatial infinity. The geometry is classical here; quantum backreaction is a further problem.

The real solution space has the conserved symplectic form on scalar-field solutionsIts complexification defines the Klein-Gordon inner productThe divergence of its current is $i(u^{\ast }\square v-v\square u^{\ast })=0$ by the field equation, proving independence of the chosen Cauchy hypersurface. This Hermitian form is not positive definite on the whole complex solution space: conjugating a solution reverses its norm.

Choose a complete positive-norm mode space and basis $u_i$ satisfyingHere completeness means that every real solution expands in these modes and their conjugates. Equivalently, one chooses a compatible complex structure on the Klein-Gordon solution space; this is extra data, not fixed by the field equation. Quantization replaces the mode coefficients by annihilation operators and creation operators:In a foliation with induced metric determinant $h$, the canonical momentum is $\pi =\sqrt{h}{n}^a{\nabla }_a\varphi$. The normalized complete mode expansion implements the equal-time canonical commutation relations $[\hat{\varphi }(x),\hat{\pi }(y)]=i{\delta }^3(x-y)$ and vanishing field-field and momentum-momentum commutators on $\Sigma$. Propagation by the field equation gives microcausality: field observables commute at spacelike separation. Thus there is a local quantum field theory before one has selected a particle interpretation.

The chosen positive-mode space completes to a one-particle Hilbert space ${\mathcal{H}}_1$. Its bosonic Fock space is ${⨁}_{n=0}^{\infty }{\mathrm{Sym}}^n{\mathcal{H}}_1$. The Fock vacuum obeys $a_i|0⟩=0$; acting with $a_i^{†}$ produces multiparticle states, and $N_i=a_i^{†}{a}_i$ counts particles in mode $i$. The vacuum two-point function is $W(x,y)={\sum }_i{u}_i(x)u_i(y{)}^{\ast }$. The Hadamard condition restricts physically admissible states by prescribing their local short-distance singularity; it permits renormalization of the stress-energy tensor but does not select one unique vacuum.

The particle ambiguity is that there is generally no distinguished choice of positive-mode space. A second complete normalized basis can beThe Bogoliubov transformation preserves the mode normalization precisely whenIts operators obey $b_i=(v_i,\hat{\varphi }{)}_{KG}={\sum }_j({\alpha }_{ij}^{\ast }{a}_j-{\beta }_{ij}^{\ast }{a}_j^{†})$. Consequently a state with no $a$-particles generally has $b$-particles if $\beta \ne 0$. A unitary change of basis within the same positive-mode space has $\beta =0$ and merely relabels modes; a change mixing the positive and negative subspaces changes the vacuum and particle notion. Arbitrary coordinates called “time” do not resolve this choice. This is the vacuum ambiguity in a nonstationary spacetime, while the field equation and local commutator remain unchanged.

In the standard stationary setting, a complete future-directed timelike Killing vector field $K$ generates a time-translation symmetry. With a stable positive energy operator and suitable boundary conditions, choose positive-frequency solutions byThe positive stationary Hamiltonian selects a particle splitting: positive spectral modes define annihilation operators, and the normal-ordered Hamiltonian is $H_K={\sum }_i{\omega }_i{a}_i^{†}{a}_i$. Its ground state is the vacuum state in a stationary spacetime. The Killing symmetry transports this splitting unchanged, so no positive-negative mixing is induced by stationary evolution. Degeneracies allow unitary basis changes without changing the vacuum. Stationarity fixes the ground-state particle interpretation relative to this time flow; it does not forbid excited or thermal states.

The qualification is that asymptotic stationarity alone, in the sense used for a Kerr exterior, is insufficient for a global positive-energy construction: $K$ becomes spacelike in an ergoregion of a stationary spacetime. The usual unambiguous construction applies to a strictly stationary spacetime with the stated stability and spectral assumptions, or to appropriate stationary asymptotic regions for scattering. Zero modes or instability need separate treatment. Likewise different observer time flows on different regions need not define the same particles.

For particle creation by a nonstationary spacetime, choose normalized positive-frequency modes $u_j^{\mathrm{in}}$ in the stationary early region and $u_i^{\mathrm{out}}$ in the stationary late region. Propagate both sets by the Klein-Gordon equation through the intermediate time-dependent geometry. Compute their overlaps on any common Cauchy hypersurface:The initial in-vacuum has $a_j^{\mathrm{in}}|0_{\mathrm{in}}⟩=0$. Since $a_i^{\mathrm{out}}={\sum }_j({\alpha }_{ij}^{\ast }{a}_j^{\mathrm{in}}-{\beta }_{ij}^{\ast }{a}_j^{\mathrm{in}†})$, contracting the in operators yields the particle number from Bogoliubov coefficients:Thus the procedure is to fix asymptotic mode choices, solve the classical wave equation, project using the conserved inner product, and take squared negative-frequency overlaps. Particle energy comes from the time-dependent background, so particle creation does not contradict the linear deterministic field equation. For an initial state diagonal in in-mode occupation numbers $n_j$, with no anomalous correlations, the more general expectation is ${\sum }_j[|{\alpha }_{ij}{|}^2{n}_j+|{\beta }_{ij}{|}^2(n_j+1)]$.

For continuous spectra, sums become integrals and normalized wave packets avoid meaningless squares of delta functions. A finite total count requires the appropriate convergence of the double sum or integral. The bosonic mode mixing implementability condition is that $\beta$ be Hilbert-Schmidt; otherwise the two Fock representations need not be unitarily equivalent, even though the local field algebra and finite-mode particle calculations still make sense. This distinction between deterministic field evolution and a chosen representation is essential to quantum field theory in curved spacetime.

Work initially in units $G=c=\hslash =k_B=1$ and consider gravitational collapse settling to a nonextremal Schwarzschild black hole. Classically, outgoing rays sufficiently close to the forming event horizon suffer an arbitrarily large redshift. Quantum mechanically, the meaning of positive frequency at past null infinity differs from that at future null infinity; a Bogoliubov transformation between these two mode descriptions can therefore mix creation operators and annihilation operators. Hawking's argument computes that mixing and obtains a thermal occupation, even when the incoming state has no particles.

To see why the frequency mixing is universal, use final Schwarzschild time $t$, the Schwarzschild tortoise coordinate $r_{\ast }$, and retarded time $u=t-r_{\ast }$. Near the final event horizon, $r_{\ast }\sim (2\kappa {)}^{-1}\log (r-2M)$, with Schwarzschild surface gravity $\kappa =1/(4M)$. On a fixed advanced-time section, $u=v-2r_{\ast }$, so $r-2M$ is proportional to $e^{-\kappa u}$. Tracing a ray back through the regular collapsing region maps this small separation smoothly to its separation from the last escaping incoming ray. If $U$ is an affine incoming null coordinate and $U_H$ labels that limiting ray, the Hawking exponential ray map isThe logarithm at a simple horizon root fixes the exponent. Smooth propagation through the earlier collapse changes $A$ and phases but not the leading exponential redshift.

For a massless bosonic test scalar field, an outgoing mode of frequency $\omega >0$ at future null infinity behaves as $e^{-i\omega u}$. Trace it backward in the high-frequency geometric optics approximation. Its relevant incoming dependence is, up to normalization and a phase,Here $\Theta$ is the Heaviside step function. Incoming rays with $U>U_H$ do not escape. This truncated logarithmic phase is not purely positive frequency with respect to $U$. Decomposing it into incoming positive-frequency solutions $e^{-i{\omega }^{\prime }U}$ and negative modes $e^{+i{\omega }^{\prime }U}$ gives the Bogoliubov transformation. With $x=U_H-U$, the two relevant Fourier transform integrals areThe former is the positive-frequency coefficient and the latter is the negative-frequency coefficient; the Klein-Gordon inner product introduces the same frequency normalization into their modulus ratio. By the defining integral of the Gamma function and analytic continuation in its Laplace parameter,For a complex number $z$ in the right half-plane, $|z^{-1-ia}|=|z{|}^{-1}{e}^{a\arg z}$. As $ϵ\downarrow 0$, the two arguments tend to $+\pi /2$ and $-\pi /2$. This gives the thermal ratio of Hawking Bogoliubov coefficients:The regulator fixes the branches and is essential to the sign of the thermal exponent.

Expand the outgoing annihilator as $b=\int d{\omega }^{\prime }({\alpha }_{\omega {\omega }^{\prime }}^{\ast }{a}_{{\omega }^{\prime }}-{\beta }_{\omega {\omega }^{\prime }}^{\ast }{a}_{{\omega }^{\prime }}^{†})$, including a complete set of incoming channels. The canonical identities for a bosonic Bogoliubov transformation express $[b,b^{†}]=1$. For normalized outgoing wave packets, they give $\int (|\alpha {|}^2-|\beta {|}^2)d{\omega }^{\prime }=1$. In the incoming Fock vacuum, the outgoing mean of the particle number operator is $N_{\omega }=\int |\beta {|}^{2}d{\omega }^{\prime }$. Combining the two identities with the thermal ratio, in the late-time narrow-frequency packet limit, yieldsWave packets avoid interpreting the divergent normalization of infinitely long continuum modes as a finite particle count. The temperature is the Hawking temperature, and the occupation is the zero-chemical-potential Bose-Einstein distribution. Restoring constants, the Schwarzschild black hole temperature is $T_H=\hslash c^3/(8\pi GM{k}_B)$ for physical mass $M$.

The preceding thermal occupation describes the horizon-originating channel before exterior scattering. The exterior angular and curvature potential partly reflects it. For a scalar field, the number flux at future null infinity is consequentlywhere the greybody factor ${\Gamma }_{\ell }$ is the transmission probability of a partial wave. Multiplication by $\omega$ gives its energy flux. Thus the asymptotic spectrum has the thermal denominator but is not an exact blackbody spectrum. For fermionic fields the canonical anticommutation relations replace the minus sign in the occupation denominator by a plus sign, at the same Hawking temperature.

The relevant state is the collapse vacuum for Hawking radiation: no incoming thermal bath at past null infinity, and regular short-distance behavior in freely falling coordinates at the forming event horizon. It predicts an outgoing flux, unlike equilibrium with a thermal bath on an eternal black hole. The derivation uses quantum field theory in curved spacetime on a prescribed classical geometry; it does not require outgoing particles to follow forbidden classical paths out of the interior. Positive energy carried to infinity is accompanied, in the semiclassical description, by negative Killing-energy flux into the black hole. Including that flux in slow backreaction decreases its mass, giving black-hole evaporation.

The approximation applies to the late radiation of a large nonextremal black hole while its surface gravity changes slowly on a time scale ${\kappa }^{-1}$. Backward propagation involves very high locally measured frequencies, so the calculation assumes the usual regular short-distance quantum-field state. It determines the leading radiation and temperature, not the Planck-scale endpoint of black-hole evaporation. The exponential horizon redshift forces positive/negative frequency mixing with a thermal ratio, giving outgoing Hawking radiation at $T_H=\kappa /(2\pi )$.

Canonical mode mixing and its inverse. Treat annihilation operators as column vectors and keep complex conjugation, transpose and Hermitian adjoint distinct. Computing the canonical commutation relations for $a^{\prime }=\alpha a+\beta a^{†}$ givesConsequently the required conditions areFor complete invertible mode mixing, putThe canonical identities for a bosonic Bogoliubov transformation are $\mathcal{T}J{\mathcal{T}}^{†}=J$, giving ${\mathcal{T}}^{-1}=J{\mathcal{T}}^{†}J$. Its upper blocks yieldEquivalently, the column identities are ${\alpha }^{†}\alpha -{\beta }^T\overline{\beta }=I$ and ${\alpha }^{†}\beta ={\beta }^T\overline{\alpha }$. In finite dimension invertibility follows from the canonical matrix identity. In infinitely many modes, the row identities alone need not give a complete inverse: $a_i^{\prime }=a_{i+1}$ preserves them but discards one mode. Here completeness of the two field-mode expansions supplies the needed invertible transformation.

Obtaining coefficients from the modes. Use the Klein-Gordon inner product, antilinear in its first argument,The normalized positive-frequency solutions satisfy $(p_i,p_j)={\delta }_{ij}$, $({\overline{p}}_i,{\overline{p}}_j)=-{\delta }_{ij}$ and $(p_i,{\overline{p}}_j)=0$. Current conservation makes this product independent of the Cauchy hypersurface when boundary flux vanishes. Taking the product of the field with a mode gives $a_i=(p_i,\hat{\Phi })$ and henceWith this operator convention, the corresponding mode expansion isThe minus sign comes from the negative norm of conjugate modes. Using a different convention for mode coefficients can move this sign and the complex conjugates; the operator formula fixes them unambiguously here.

Particle count. In the in-vacuum, only the contraction $⟨0|a_j{a}_k^{†}|0⟩={\delta }_{jk}$ survives. ThereforeThis is particle number from Bogoliubov coefficients: the old vacuum contains that expected number of particles in the new mode $i$. Nonzero negative-frequency mixing is the source of particle production. For continuum modes one uses normalized wave packets and replaces the sum by the corresponding integral.

The squeezed-state relation. First use finitely many modes, or an implementable infinite-mode limit. SeekThe canonical commutation relations giveAll further nested commutators vanish since $\hat{F}$ contains only creation operators. ThusThe old annihilation conditions hold exactly whenThe inverse transformation's canonical identity $A{B}^T=B{A}^T$ shows $A^{-1}B=(A^{-1}B{)}^T$, so the required symmetry is automatic. Also $A{A}^{†}-B{B}^{†}=I$ givesHence the singular values of $M$ are less than one. The Autonne-Takagi factorization $M=U\mathrm{diag}(s_j)U^T$ changes to independent canonical oscillators. For each factor, the squared norm of $\exp (s_j{b}_j^{†2}/2)|0⟩$ is the even-occupation seriesMultiplying these norms gives the finite-mode normalizationThe phase of $C$ is arbitrary. The resulting multimode squeezed vacuum is annihilated by all old annihilation operators, so uniqueness of the normalized Fock vacuum identifies it with $|0⟩$.

The PDF's squeezed-vacuum claim needs qualification for infinitely many modes. Bosonic mode mixing implementability requires $\beta$ to be a Hilbert-Schmidt operator, ${\sum }_{ij}|{\beta }_{ij}{|}^2<\infty$, for a common ordinary bosonic Fock space. For a counterexample, take $\alpha =\cosh sI$ and $\beta =\sinh sI$ on countably many modes with fixed $s>0$. The algebraic canonical conditions hold, but $M=\tanh sI$ and the finite-$N$ normalization is $C_N=(\cosh s{)}^{-N/2}\to 0$. Every mode has a fixed positive probability of nonzero occupation in the required product state; the probability that all but finitely many modes are empty is zero. Ordinary Fock space vectors instead have total occupation finite with probability one, even if their expected occupation is infinite. Thus there is no nonzero common-Fock-space vector of the prescribed form. The squeezed expression is a normalizable vacuum relation with a finite-mode regulator or the implementability condition, not solely from the commutator identities.

Black-hole radiation. In a collapse spacetime, choose early positive-frequency solutions with respect to affine incoming time $v$ on past null infinity, and late outgoing modes proportional to $e^{-i\omega u}$ on future null infinity. The early state is their in-vacuum. The Hawking exponential ray map hasfor late outgoing rays, with $\kappa >0$. Tracing a late mode backwards therefore gives a profile proportional toIt is not a pure positive-frequency incoming wave, so the Bogoliubov transformation has nonzero $\beta$.

To see the thermal factor, put $x=v_H-v>0$ and $a=\omega /\kappa$. The positive- and negative-frequency Fourier pieces have, apart from common normalization and phases, the regulated integralsAs $\epsilon \downarrow 0$, the positive-frequency piece has squared-modulus factor $e^{+\pi a}$ and the negative-frequency piece has $e^{-\pi a}$. Consequently the thermal ratio of Hawking Bogoliubov coefficients isFor a normalized narrow-frequency outgoing wave packet, combine this ratio with ${\sum }_j(|{\alpha }_{ij}{|}^2-|{\beta }_{ij}{|}^2)=1$. The occupation is the Bose-Einstein distributionThe Hawking temperature is thus the same one found by Euclidean regularity. For Schwarzschild, $\kappa =1/(4M)$ gives $T_H=1/(8\pi M)$ in natural units. Propagation through the exterior potential multiplies the asymptotic occupation by the greybody factor ${\Gamma }_{\ell }(\omega )$.

The outgoing radiation has partner-mode correlations across the horizon: the total state can be pure while the reduced outgoing state is thermal. A complete mode basis must include modes entering the horizon as well as those reaching infinity. The particle statement is made for normalized packets and finite observation intervals; idealized infinite-duration emission need not define the two vacua in one global Fock space. This explains both the squeezed-state structure and the physically measurable Hawking radiation.

Fix a classical globally hyperbolic spacetime with metric signature $(-+++)$; in this solution take $\hslash =1$. A free real Klein-Gordon field can be specified bywhere $\mu$ is its mass, $\xi$ its curvature coupling and $\mathcal{R}$ the scalar curvature. Global hyperbolicity ensures a well-posed initial-value problem on a Cauchy hypersurface $\Sigma$ and the existence of retarded and advanced propagators. Thus compactly supported field and normal-derivative data determine a classical solution. This fixes the dynamics, but not a Fock vacuum.

For real solutions with suitable support or falloff, the symplectic form on scalar-field solutions iswith $n$ the future unit normal. The field equation makes the current conserved, so this symplectic form is independent of $\Sigma$ when boundary flux vanishes. Quantize the initial data by the canonical commutation relation: with $\pi =n^a{\nabla }_a\varphi$ and the delta function defined relative to $d\Sigma$, $[\hat{\varphi }(x),\hat{\pi }(y)]=i{\delta }_{\Sigma }(x,y)$ and the two equal-field commutators vanish. Equivalently, construct the field algebra using the causal propagator. A state on that algebra is additional input.

To construct a particle representation, complexify the classical solution space. Its conserved Klein-Gordon inner product isIt is indefinite on all complex solutions. Choose a complete positive-norm subspace and an orthonormal mode basis $u_i$ with $(u_i,u_j)={\delta }_{ij}$, $(u_i^{\ast },u_j^{\ast })=-{\delta }_{ij}$ and $(u_i,u_j^{\ast })=0$. Such a choice is encoded by a compatible complex structure on the Klein-Gordon solution space. Its positive subspace gives the one-particle Hilbert space, and the associated bosonic Fock space contains symmetrized many-particle states. The field expansion isThe creation operator $a_i^{†}$ adds a particle in mode $i$, the annihilation operator removes one, and the number operator is $N_i=a_i^{†}{a}_i$. For continuous mode labels the sums and Kronecker deltas become integrals and delta functions, or one can work with normalized wave packets.

The ambiguity is precisely that the field equation and global hyperbolicity do not select that positive subspace. A different normalized basis may mix $u_j$ and $u_j^{\ast }$ by a Bogoliubov transformation, and then its annihilation operators mix $a_j$ and $a_j^{†}$. Its Fock vacuum and number operators differ. The Hadamard condition constrains physically acceptable short-distance singularities and allows local renormalization, but it still leaves many states. Hence there is generally no observer-independent particle count on an arbitrary dynamical geometry.

In a stable strictly stationary spacetime, a globally future timelike Killing vector field $K$ gives a preferred time translation. Fix its normalization and suitable boundary conditions, and choose positive-frequency solutions satisfying $i{\mathcal{L}}_{K}u=\omega u$ with $\omega >0$. The corresponding positive spectral subspace gives the preferred vacuum state in a stationary spacetime. Unitary changes of basis within it leave the vacuum and the particle notion unchanged. This construction assumes a well-defined positive stationary generator; stationarity by itself is insufficient if $K$ becomes spacelike, as in a Kerr ergoregion, or if unstable or zero modes obstruct the ground-state construction. It selects a preferred ground state under the stated assumptions, not a unique state among all thermal and excited states.

If the geometry is suitably stationary in the asymptotic past and future, choose those preferred mode spaces separately, giving the in-vacuum and out-vacuum. Propagate the past modes through the intervening region using the field equation and compare them to the future modes using the conserved Klein-Gordon inner product. Adopt the conventionThe canonical identities for a bosonic Bogoliubov transformation read $\alpha {\alpha }^{†}-\beta {\beta }^{†}=I$ and $\alpha {\beta }^T=\beta {\alpha }^T$. Extracting the future annihilation operator with the same inner product givesIn the in-vacuum, only $⟨a_j{a}_k^{†}⟩={\delta }_{jk}$ contributes to $⟨b_i^{†}{b}_i⟩$. Therefore the particle number from Bogoliubov coefficients isNonzero $\beta$ is the production of future particles from the past vacuum. Summing over future modes gives the total expected particle number when that sum is finite. For infinitely many modes, a Hilbert-Schmidt operator $\beta$ is the condition for unitary implementability between these pure bosonic Fock space representations; finite-volume or wave-packet calculations must respect the relevant measures and convergence. Particle production is determined by the negative-frequency mixing, rather than by identifying a single instantaneous vacuum throughout the time-dependent region.