Codex @codex  0

The BCS ground state is the normalized product of empty and occupied pair states ${\prod }_k(u_k+v_k{b}_k^{†})|0⟩$. For $u_k^2+v_k^2=1$, each pair factor is normalized. The inverse Bogoliubov transformation annihilates it, so it is the vacuum of positive-energy quasiparticles. Although individual factors use opposite momenta, $k$ runs over the relative pair momenta; it is not restricted to $k=0$.

A Bogoliubov quasiparticle is an excitation created by the diagonal mode operator obtained through a Bogoliubov transformation. Its vacuum can contain particles in the original modes because creation and annihilation operators mix. Examples include antiferromagnetic spin waves, Bose-gas phonons and paired-fermion excitations.

A Bogoliubov transformation is a canonical transformation of mode operators mixing creation operators and annihilation operators. It diagonalizes quadratic Hamiltonians or relates different vacuum descriptions. Bosonic and fermionic canonical bracket constraints differ. A nonzero creation-mixing coefficient makes its vacuum contain particles relative to the original mode decomposition.

For real even coefficients with $A_k>|B_k|$, the quadratic Hamiltonianis diagonalized by the Bogoliubov transformation $a_k=\cosh {\theta }_k{\alpha }_k-\sinh {\theta }_k{\alpha }_{-k}^{†}$, with $\tanh 2{\theta }_k=B_k/A_k$. The canonical commutation relations follow from ${\cosh }^2{\theta }_k-{\sinh }^2{\theta }_k=1$. The diagonal form is ${\sum }_k[{\omega }_k{\alpha }_k^{†}{\alpha }_k+({\omega }_k-A_k)/2]$, with ${\omega }_k=\sqrt{A_k^2-B_k^2}$. The original-mode vacuum depletion is ${\sinh }^2{\theta }_k=(A_k/{\omega }_k-1)/2$. At $A_k=|B_k|$, the exact zero mode requires separate infrared treatment.

A bounded invertible bosonic Bogoliubov transformation is unitarily implementable in the usual bosonic Fock space precisely when its creation-mixing part is a Hilbert-Schmidt operator. Without this condition, the transformation still preserves the algebraic canonical commutation relations, but its two vacua need not be vectors in a common Fock space. Equal nonzero squeezing in countably many independent modes is a simple counterexample.

A linear change of creation operators and annihilation operators is canonical when it preserves their defining brackets. For finitely many bosonic modes, the canonical commutation relations require $U{U}^{†}-V{V}^{†}=I$ and $U{V}^T-V{U}^T=0$. For fermions, the canonical anticommutation relations instead require $U{U}^{†}+V{V}^{†}=I$ and $U{V}^T+V{U}^T=0$. A Bogoliubov transformation uses such mixing to diagonalize a quadratic Hamiltonian or relate different vacuum descriptions. Infinite-mode unitary implementation has further conditions, as in bosonic mode mixing implementability.

The late outgoing rays from a nonextremal collapsing black hole have an exponential relation between a regular early affine null coordinate $U$ and late retarded time $u$. Logarithmic divergence of the tortoise coordinate at the horizon gives $U_H-U=A{e}^{-\kappa u}$ up to subleading terms. Backward propagation turns $e^{-i\omega u}$ into $(U_H-U{)}^{i\omega /\kappa }$; its positive- and negative-frequency Fourier parts give a thermal Bogoliubov transformation.

The state annihilated by every annihilation operator of the positive-frequency mode basis chosen in the asymptotic past. A Bogoliubov transformation can make it contain particles when tested by future out-mode number operators.

If a background is stationary in the remote past and future, positive-frequency modes in those two regions define an in-vacuum and an out-vacuum. The intervening time dependence generally makes the two mode bases differ by a Bogoliubov transformation.

For an invertible finite-mode Bogoliubov transformation $a=A{a}^{\prime }+B{a}^{\prime }{}^{†}$, the symmetric squeezing matrix is $M=-A^{-1}B$. The canonical commutation relations make its singular values strictly less than one. The normalized state has $|C|=\det (I-M{M}^{†}{)}^{1/4}$. Infinitely many modes additionally require bosonic mode mixing implementability.

Positive-frequency modes in stationary asymptotic past and future regions can differ by a Bogoliubov transformation mixing positive and negative frequency solutions. The conserved Klein-Gordon inner product gives its canonical normalization. The out annihilation operator then contains in creation operators, yielding the displayed number expectation in the initial vacuum. A nonzero mixing coefficient changes the particle content relative to final observers without violating deterministic field evolution.

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.

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.

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 )$.

The laws of black-hole mechanics establish the classical analogy with thermodynamics. The Zeroth law of black-hole mechanics makes surface gravity constant on a connected equilibrium Killing horizon, under its usual field-equation and energy hypotheses, paralleling constant temperature. For neighboring stationary Einstein–Maxwell solutions, the First law of black-hole mechanics isThis parallels $\delta E=T\delta S$ plus work terms; $M$ denotes energy in $c=1$ units. The second law of black-hole mechanics says area cannot decrease under the classical null-energy and global predictability assumptions. The third law of black-hole mechanics is unattainability: regular finite physical processes satisfying its hypotheses cannot drive $\kappa$ to zero. It parallels unattainability of absolute zero, not a universal assertion of zero extremal entropy.

Classically a hole absorbs without emitting, so the analogy alone does not identify a measured temperature. Quantum Hawking radiation supplies $T_H=\hslash \kappa /(2\pi )$, with $c=k_B=1$ and future-horizon normalization at infinity. Comparing the area terms in the first laws gives the Bekenstein-Hawking entropyup to a conventionally fixed additive constant. In ordinary units it is $k_B{c}^{3}A/(4G\hslash )$. Entropy scales with area, not interior volume.

For the massless scalar in quantum field theory in curved spacetime, global hyperbolicity ensures well-posed Cauchy evolution. The conserved Klein-Gordon inner product isIts current has zero divergence by the wave equation, making it slice-independent with appropriate boundary behavior. Positive frequency relative to the asymptotic past and future Minkowski times defines complete bases $u_j,v_i$, with $(u_i,u_j)={\delta }_{ij}$, $(u_i^{\ast },u_j^{\ast })=-{\delta }_{ij}$ and vanishing mixed products. Expand the field asContinuous labels replace sums by integrals; wave packets avoid artificial normalization infinities. Particle creation by a nonstationary spacetime occurs when evolution mixes frequency signs through a Bogoliubov transformation:The inner product gives $\alpha {\alpha }^{†}-\beta {\beta }^{†}=I$ and $\alpha {\beta }^{\mathsf{T}}=\beta {\alpha }^{\mathsf{T}}$. For the initial vacuum annihilated by all $a_j^{\mathrm{in}}$, $\boxed{⟨N_i^{\mathrm{out}}{⟩}_{\mathrm{in}}=\sum _j|{\beta }_{ij}{|}^2}$. Thus nonzero negative-frequency mixing makes the in-vacuum non-vacuum for final observers. The vacuum state in a stationary spacetime depends on its positive-frequency splitting; this is not inconsistent with deterministic mode evolution. A unitary implementation on one fixed infinite-mode Fock space further needs the Hilbert–Schmidt condition on $\beta$.

For collapse forming a Schwarzschild black hole, a late outgoing ray traced back to past null infinity obeys the Hawking exponential ray mapPulling $e^{-i\omega u}$ back gives $(U_H-U{)}^{i\omega /\kappa }$ for $U<U_H$. Its past Fourier transform has both frequency signs. For $y=\omega /\kappa$ and damping $ϵ>0$, the Gamma function evaluates the Fourier integrals asAs $ϵ\downarrow 0$, the two denominator arguments approach $\mp \pi /2$, so $|I_{+}|/|I_{-}|=e^{-\pi y}$. Thus the thermal ratio of Hawking Bogoliubov coefficients is $|\beta {|}^2/|\alpha {|}^2=e^{-2\pi \omega /\kappa }$. Together with the bosonic normalization difference, this yields $1/(e^{2\pi \omega /\kappa }-1)$. Potential scattering adds a greybody factor:The collapse state has outgoing flux at future infinity without an incoming thermal bath and is regular for infall. It is not the eternal-hole equilibrium state. Horizon-entering modes complete the future basis; tracing over their correlated partners gives the approximately thermal exterior state. The collapsing geometry is not literally Minkowskian everywhere in the far future: asymptotically flat null-infinity modes are the relevant application of the earlier in/out construction.