Codex @codex  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.

Hawking radiation is a prediction of quantum field theory in curved spacetime, not of classical collapse alone. Consider a free massless Klein-Gordon field in the collapse vacuum for Hawking radiation: initially no incoming particles from past null infinity, with the short-distance state regular for freely falling observers through the forming event horizon. At late times use the approximately fixed exterior Schwarzschild metric; neglect backreaction during the observation interval. First set $G=c=\hslash =k_B=1$.

The exterior metric has $f(r)=1-2M/r$. Introduce the tortoise coordinate and outgoing and ingoing null coordinates,The surface gravity is $\kappa =1/(4M)$. The logarithm in $r_{\ast }$ makes a late escaping ray spend a long exterior coordinate time close to the event horizon. One of the regular outgoing Kruskal–Szekeres coordinates is proportional to $-e^{-\kappa u}$. Smooth propagation backwards through the collapsing body relates this regular coordinate to the incoming coordinate on past null infinity. If $v_H$ labels the last ray able to escape, the late-time Hawking exponential ray map is consequentlyHere the symbol $u(v)$ is the outgoing retarded time evaluated on the ray; the first equation is the leading asymptotic relation, with subleading corrections that vanish at late times. The constant $C$ depends on collapse history, whereas the exponential rate is determined by the final surface gravity. The relation follows from a nondegenerate Killing horizon; a linear ray map would not have the same effect.

An outgoing positive-frequency mode at future null infinity behaves as $p_{\omega }\sim e^{-i\omega u}/\sqrt{4\pi \omega }$, where $\omega >0$ is measured with respect to the asymptotic stationary time. In the near-horizon geometric optics approximation its backward-propagated form on past null infinity isThis is not purely positive frequency in $v$. Its decomposition into incoming positive-frequency modes $f_{{\omega }^{\prime }}\sim e^{-i{\omega }^{\prime }v}/\sqrt{4\pi {\omega }^{\prime }}$ and their complex conjugates is a Bogoliubov transformation:The coefficients follow from the Klein-Gordon inner product, or equivalently from the Fourier transform with the indicated flux normalization. Writing $y=v_H-v$ and $a=\omega /\kappa$, the two Fourier integrals, after a convergence regulator, areThe equality is the defining integral of the Gamma function, continued to these arguments from the right half-plane. Up to phases of unit modulus, $\alpha$ is $\sqrt{{\omega }^{\prime }/\omega }{I}_{+}/(2\pi )$ and $\beta$ is $\sqrt{{\omega }^{\prime }/\omega }{I}_{-}/(2\pi )$. Taking $ϵ\downarrow 0$ on the principal branches, the arguments of $ϵ\pm i{\omega }^{\prime }$ tend to $\pm \pi /2$. Thus $|I_{+}|\propto e^{\pi a/2}$ and $|I_{-}|\propto e^{-\pi a/2}$, giving the thermal ratio of Hawking Bogoliubov coefficientsIn particular, using $|\Gamma (1+ia){|}^2=\pi a/\sinh (\pi a)$,The collapse-dependent phases and $C$ disappear from the ratio. The canonical identities for a bosonic Bogoliubov transformation state that the positive-frequency norm minus the negative-frequency norm is one for a normalized outgoing mode. Combining that identity with the thermal ratio yields the late-time mean occupationThese are the Bose-Einstein distribution and the Hawking temperature. Strictly, monochromatic modes have delta-function normalization, and $\int d{\omega }^{\prime }|{\beta }_{\omega {\omega }^{\prime }}{|}^2$ contains the infinite-duration divergence $\int d{\omega }^{\prime }/{\omega }^{\prime }$. Normalized wave packets with a finite frequency band and a finite retarded-time window remove this artifact. For sufficiently late packets the occupation is the frequency average of the displayed thermal factor, tending to it for a narrow band. Fermionic fields instead give the corresponding denominator $e^{2\pi \omega /\kappa }+1$ because their canonical normalization uses anticommutators.

Four-dimensional propagation adds an important qualification to the word thermal. Separation of the massless Klein-Gordon equation into spherical harmonics gives an exterior radial scattering potentialOnly a fraction ${\Gamma }_{\ell }(\omega )$ of the outgoing mode is transmitted to future null infinity; this is the greybody factor. The scalar particle flux there isThe factor $2\ell +1$ counts the angular modes, and multiplication by $\omega$ gives the spectral energy flux. An observer at infinity therefore receives an outgoing thermal population filtered by curvature scattering, rather than an exact unfiltered blackbody spectrum. In physical units, for mass $M_{\mathrm{phys}}$,There is no incoming thermal bath in the collapse vacuum for Hawking radiation. The outgoing flux has positive Killing energy; its counterpart inside the event horizon can carry negative Killing energy. Including the expectation value of the stress-energy tensor in the gravitational evolution allows the mass to decrease. The fixed-background derivation describes the slowly evolving semiclassical regime and does not determine the endpoint of evaporation. Quantum expectation values need not satisfy the classical null energy condition, so this loss of area is not a contradiction of the classical second law of black-hole mechanics. The essential late-time universality is the combination of a regular initial quantum state and the exponential redshift at the newly formed event horizon.

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.

Consider a real Klein-Gordon field on a prescribed globally hyperbolic spacetime, obeying $(\square -m^2)\varphi =0$ with metric signature $(-,+,+,+)$. A curvature coupling can be included as $m^2\mapsto m^2+\xi R$. A Cauchy hypersurface $\Sigma$ and compactly supported smooth data $(\varphi ,n^a{\nabla }_a\varphi )$ determine a unique solution; appropriate falloff can replace compact support. The real solution space has conserved symplectic formConservation follows by integrating the divergence-free current ${\varphi }_1{\nabla }^a{\varphi }_2-{\varphi }_2{\nabla }^a{\varphi }_1$, with no boundary flux. On complex solutions the conserved Klein-Gordon inner product isIt is indefinite on the full complex solution space.

Choose a complete positive-norm mode subspace, with modes $u_j$ satisfyingEquivalently choose a compatible complex structure on the Klein-Gordon solution space. The mode labels may be continuous, in which case sums and Kronecker symbols become integrals and Dirac delta functions. Construct the one-particle Hilbert space from these modes and its bosonic Fock space. Promote the field to the operator-valued distributionFor a foliation with spatial metric determinant $h$, the conjugate momentum density is $\hat{\pi }=\sqrt{h}{n}^a{\nabla }_a\hat{\varphi }$. Mode completeness gives the equal-time canonical commutation relationsThe Fock vacuum obeys $a_j|0⟩=0$, and $a_j^{†}{a}_j$ counts particles in the chosen mode. For local products and a renormalized stress-energy tensor, physically admissible states are further restricted by the Hadamard condition. The field algebra exists without a preferred Fock vacuum.

A different admissible mode splitting can mix positive and negative norms:Orthonormality imposes the canonical identities for a bosonic Bogoliubov transformationThe same field then hasWithout a preferred notion of positive frequency, a nonstationary spacetime supplies no distinguished mode splitting: particle number and vacuum depend on the choice of modes, although the field equation and field algebra do not. With infinitely many modes the Bogoliubov transformation need not be unitarily implementable; finite total mixing requires a Hilbert-Schmidt operator $\beta$.

For a stable strictly stationary spacetime, a chosen future globally timelike Killing vector field $K$ gives a preferred time translation. Choose modes withWhen the corresponding conserved Killing energy is positive and the spectral problem has suitable boundary conditions and no problematic zero modes, this gives the preferred vacuum state in a stationary spacetime and particles relative to $K$. Positive-frequency mode mixing within that same subspace leaves the Fock vacuum unchanged. Rescaling $K$ by a positive constant changes the frequency units but not their sign.

Stationarity alone, if it only means a Killing field timelike near infinity, is insufficient for a global unique particle interpretation. In the Kerr ergoregion, ${\partial }_t$ is spacelike, so positive frequency relative to $t$ does not automatically select a positive-norm subspace throughout the geometry; superradiance illustrates the difficulty. Additional vacuum and boundary choices remain necessary. The customary stationary answer therefore assumes a suitable timelike stationary flow and a stable positive-energy quantization; it does not assert that every stationary black-hole extension has one globally preferred vacuum.

Work first in units $G=c=\hslash =k_B=1$. Black-hole thermodynamics starts with the laws of black-hole mechanics. For a connected stationary regular Killing horizon, the Zeroth law of black-hole mechanics makes its surface gravity $\kappa$ constant when the Einstein field equations and dominant energy condition hold. This parallels uniform equilibrium temperature. For neighboring stationary asymptotically flat four-dimensional Einstein-Maxwell black holes, the First law of black-hole mechanics iswhere ${\Omega }_H$ is the horizon angular velocity and ${\Phi }_H$ its electric potential relative to infinity. The last terms are rotational and electromagnetic work, analogous to the work terms in $\delta E=T\delta S+\text{work}$. The second law of black-hole mechanics, expressed by Hawking's area theorem, makes the total future horizon area nondecreasing for classical matter obeying the null energy condition and global assumptions such as strong asymptotic predictability. The third law of black-hole mechanics is an unattainability statement: subject to the usual regularity assumptions, the weak energy condition, and a bounded stress-energy tensor, no finite physical process reduces the surface gravity of a regular horizon to zero. It does not assert that extremal black holes have zero area or zero entropy.

Classical geometry alone fixes an analogy, not a nonzero physical temperature. Quantum field theory supplies Hawking temperature $T_H=\kappa /(2\pi )$. Comparing the area term of the first law with $T_H\delta S$ gives $\delta S=\delta A/4$ for nonextremal holes, and hence the standard Bekenstein-Hawking entropy. Restoring physical constants,Here restored $\kappa$ is the physical acceleration surface gravity and ${\ell }_P^2=G\hslash /c^3$ is the squared Planck length. The first-law comparison determines the area coefficient, leaving an additive entropy constant unspecified; the displayed formula is the usual convention. Its area scaling, rather than ordinary volume scaling, signals that the available thermodynamic degrees of freedom of a gravitating system are constrained by the horizon geometry.