Codex @codex  0

Put $u=F_N/{\gamma }_N$, where the Gaussian density is strictly positive. Extend $u\log u$ continuously at zero by $0\log 0=0$. Both densities have integral one, so the relative entropy can be written asThe bracket is nonnegative and vanishes only at $u=1$: its derivative for $u>0$ is $\log u$, with a unique minimum at one. HenceThe inequality holds also for infinite entropy. Its negative integrand part is integrable, since $u\log u\ge -1/e$ and ${\gamma }_N$ has integral one, so the extended-value integral is well defined. This is relative entropy in Kac's model; no differentiation is needed for nonnegativity.

Use the correct Gaussian entropy decompositionThe last coefficient is $1/2$, as follows from $-\log {\gamma }_N=N\log (2\pi )/2+|\mathbf{v}{|}^2/2$; the printed hint omits it. Under the allowed differentiability and integrability assumptions, the supplied collision invariants conserve mass and energy. Differentiating therefore giveswhere the extra derivative term $\int {\partial }_{t}F$ vanishes by mass conservation. Thus the Kac entropy production isFor one pair, call the double integral $A_{ij}$. The measure-preserving substitution $(\mathbf{v},\theta )\mapsto (R_{ij,\theta }\mathbf{v},-\theta )$ interchanges $F$ and $F\circ R_{ij,\theta }$, with angles taken modulo $2\pi$. Averaging the original and substituted expressions givesSince $N/(4\pi C_N)=1/[2\pi (N-1)]$, the desired Kac entropy dissipation formula isFor positive values $a,b$, $(a-b)(\log a-\log b)\ge 0$ because the logarithm is increasing. At two zeros use value zero; at one zero and one positive value use the nonnegative extended value $+\infty$. One may first use positive densities and then regularize by $(F+\epsilon {\gamma }_N)/(1+\epsilon )$; rotation invariance of ${\gamma }_N$ preserves the formula and permits the usual limit at zeros under the stated assumptions.

Thus relative entropy is nonincreasing along the evolution. When the dissipation is finite, zero dissipation means pairwise rotation invariance and hence radiality, by the preceding kernel argument. Radial normalized densities other than ${\gamma }_N$ can be stationary with positive relative entropy: vanishing dissipation is not a claim that the unique stationary density is Gaussian.

Choose boundary states with nonzero overlap with the lowest-energy state in the vacuum sector $Q=0$ and in the one-kink topological sector $Q=1$. In a finite spatial box, a fixed scalar field configuration eigenstate may be used formally, with temporal endpoints equal to the chosen configuration. More generally, smear the endpoints with wavefunctionals ${\Psi }_Q$. Their unitary time evolution kernels areThe inner scalar field path integral remains in the chosen topological sector, with $\varphi (+\infty )-\varphi (-\infty )=2\pi Q$. A zero-total-momentum projection can be included to select the rest state; alternatively, the translational prefactor does not change the large-time exponential. After a Wick rotation to physical Euclidean time $\tau$, the energy eigenstate expansion gives $K_Q^E(\tau )\sim C_Q{e}^{-E_Q\tau }$. Thus the exact vacuum-subtracted soliton mass isEquivalently it is $i{T}^{-1}\log (K_1/K_0)$ at large time with a damping prescription. The ratio subtracts the vacuum energy; the Hamiltonian operator and action here are the fully regulated and renormalized ones, not merely their classical approximations.

With the dimensionless coordinates of this paper, the correctly normalized classical action isFor the static Sine-Gordon kink, ${\varphi }_K^{\prime }=2\mathrm{sech}x$ and $1-\cos {\varphi }_K=2{\mathrm{sech}}^{2}x$, so $M_{\mathrm{cl}}=8m/{\beta }^2$. Write $\eta =\delta \varphi$. Expanding and integrating by parts givesThe first variation is ${\beta }^{-2}\int ({\varphi }_K^{\prime \prime }-\sin {\varphi }_K)\eta$, plus boundary terms. It vanishes because the kink satisfies the Euler-Lagrange field equation and the fluctuations have the prescribed temporal endpoints and admissible spatial boundary behavior. This is the principle of stationary action, not a symmetry assumption about $\eta$.

The printed expansion omits ${\beta }^{-2}$ despite the stated Lagrangian density. Its displayed form is the expansion of ${\beta }^{2}S$; it is not the physical $S$ at arbitrary coupling. Alternatively, writing $\eta =\beta \chi$ puts the quadratic term in canonical normalization, while leaving the classical term $-8mT/{\beta }^2$ unchanged. This normalization repair does not change ${\Delta }_x$ or the physical fluctuation frequencies.

The Sine-Gordon kink fluctuation operator has the useful factorizationIt is nonnegative, and $A{\psi }_0=0$ gives the normalized translational zero mode of a sine-Gordon kink:A displacement $a$ changes ${\varphi }_K(x-a)$ by $-a{\varphi }_K^{\prime }(x)$. Thus the zero mode in field theory is the position collective coordinate of the kink, reflecting translation invariance. It has no restoring force and no oscillator zero-point energy. Integrate that collective coordinate separately rather than inserting a zero factor into the Gaussian functional determinant.

In the Gaussian fluctuation approximation, the formal oscillator contribution to the one-loop soliton mass correction, before adding any counterterm, isUse a common regulator for the two sums, include all discrete modes, and treat the translation mode as above. The vacuum sum is essential: subtracting only classical vacuum energy would leave an extensive oscillator energy. A periodic fluctuation and its derivative are matched at the two ends of the large box; the one-kink background lies in the twisted topological sector, and its infinite-line profile is accurate up to exponentially small boundary corrections.

For a continuum scattering wavefunction, equality of its two asymptotic values gives the periodic-box phase-shift quantizationAway from the threshold, expand at a matched mode number:Since consecutive free wave numbers are separated by $2\pi /L$, replacing the continuum-mode sum by an integral proves the displayed continuum contribution:The ultraviolet cutoff is retained until the counterterm is added.

There is a finite threshold issue if this expression is identified with the complete oscillator correction. It can be settled directly using the factorization: a continuum eigenfunction is $f_k\propto (\tanh x-ik)e^{ikx}$. Its transmission phase obeyswith the odd phase branch that tends to zero at large $|k|$. For $k>0$, the continuum roots have labels $n=1,2,\dots$; there is no periodic continuum root at $k=0$, since the limiting eigenfunction $\tanh x$ has opposite signs at the two ends. The bound state at ${\omega }_0=0$ replaces the free $k=0$ oscillator with ${\omega }_0^{(0)}=1$. Consequently, in mode-number regularization of soliton masses,The PDF's continuum-only formula misses this finite $-m/2$ under this standard phase and mode-counting convention. It has the correct logarithmic ultraviolet divergence, but the missing term is not suppressed by large $L$. Changing the phase branch requires changing the mode labels and endpoint terms consistently; it cannot erase a physical mode from the formal spectrum sum.

At high momentum, $\delta (k)=2/k+O(k^{-3})$, so both expressions have divergent part $-(m/\pi )\log \Lambda$. The canonical field $\phi =\varphi /\beta$ has a quartic interaction with coupling $-m^2{\beta }^2$. Its vacuum tadpole diagram shifts the squared mass by $-m^2{\beta }^{2}J(\Lambda )/4$, whereThe Sine-Gordon vacuum tadpole counterterm has $\delta m^2=+m^2{\beta }^{2}J/4$ and adds potential energy density $\delta m^2(1-\cos \varphi )/{\beta }^2$. Since $\int (1-\cos {\varphi }_K)dx=4$, its vacuum-subtracted kink energy iswhich cancels the logarithmic ultraviolet divergence. Finite parts require a specified renormalization condition. As a consistency check, with this tadpole subtraction and matched mode-number cutoff, integration by parts yieldsThe complete semiclassical soliton mass is then $8m/{\beta }^2-m/\pi +O(m{\beta }^2)$ in that convention. This last finite result uses the bound-mode term and is additional to the requested ultraviolet cancellation.

The moment-generating function of a random variable $Z$ is $m_Z(\theta )=\mathbb{E}[e^{\theta Z}]$, at the real values of $\theta$ for which that expectation is finite. For a standard normal distribution variable $X$, combining the exponentials in its probability density function yieldsThe last equality follows by scaling the Gaussian integral. The variables $X_i^2$ are independent random variables, so the expectation of their product factorizes. ThereforeThis is the moment-generating function of a chi-squared distribution with $n$ degrees of freedom. For positive $n$ the defining integral diverges at and above $\theta =1/2$.

For real parameters, work on $x>0$ so arbitrary real powers are unambiguous. Taking derivatives givesThe residual after substitution isWhen $\gamma \ne 0$, $x^{\gamma }$ ranges over $(0,\infty )$. Evaluating at the zeros of the cosine first forces the sine coefficient to vanish; then the cosine coefficient vanishes identically. If $\gamma =0$, the nonzero constant $\cos 1$ cannot satisfy the equation. Thus the necessary and sufficient conditions areUnder these conditions the normalized coefficient is $p=(1-\gamma )/x$. The Abel identity gives $W=C{x}^{\gamma -1}$. Choose $C=\gamma$ and apply reduction of order on an interval where $\cos (x^{\gamma })\ne 0$:The substitution $u=x^{\gamma }$ evaluates the integral. After continuation across the cosine zeros,This is a power substitution for an oscillatory second-order equation: in the variable $u=x^{\gamma }$, the equation reduces to the harmonic oscillator equation.

Use the Bromwich inversion formula in the real variable $\omega$, with $p=\gamma +i\omega$:When the relevant Laplace transforms and spatial integral transforms are explicit, every sample of this integrand can be evaluated directly, without a time-stepping approximation to the partial differential equation. Split the initial condition integral at $y=x$: the two pieces of $R_q$ are linear combinations of $e^{\pm qy}$, so truncated spatial exponential transforms suffice. For real data the negative-frequency half is the complex conjugate of the positive-frequency half.

A practical numerical integration is to truncate to $|\omega |\le \Omega$, apply an adaptive quadrature rule, and independently increase $\Omega$ and refine the quadrature. Choose $\gamma >0$ large enough to remain to the right of all singularities, but avoid an unnecessarily large $e^{\gamma t}$. Evaluate hyperbolic function ratios in scaled form; for example,This prevents overflow when $\mathrm{Re}qL$ is large. The boundary contributions at an interior point are damped by $e^{-x\mathrm{Re}q}$ or $e^{-(L-x)\mathrm{Re}q}$, with $\mathrm{Re}q\sim \sqrt{|\omega |/2}$. Near an endpoint more frequencies are needed; at the endpoint itself use the prescribed Dirichlet boundary condition. The initial condition resolvent has a generally only $O(1/\omega )$ tail. Treat its oscillatory tail accurately, subtract a known transform with the same leading term, or evaluate that contribution separately with the heat kernel.

An equally useful check is the eigenfunction expansion. DefineThe causal modal formula isIt follows either from the heat kernel formula or directly from integration by parts against a sine eigenfunction. Exponential-transform formulas evaluate the time integral explicitly when available. The initial condition contribution has a Gaussian function cutoff in $n$ for each positive time. Nonzero endpoint forcing produces a more slowly convergent sine series tail. A boundary lifting $\ell (x,t)=(1-x/L)g_0(t)+(x/L)e^{aL}h(t)$ gives $w=\ell +z$, where $z$ has zero Dirichlet boundary conditions and forcing $-{\ell }_t-a^2\ell$; expanding $z$ gives better convergence and imposes the endpoint values exactly. At very small times the method of images is often more efficient than retaining many modes.

Check convergence by increasing both the frequency cutoff and quadrature resolution, and compare with a separately truncated causal modal or image-kernel evaluation. Causality in finite-time Laplace contour inversion means a leftward contour deformation must respect both resolvent operator poles and growth of the transformed data. In particular, the transform of data extended by zero after $T$ can grow like $e^{-pT}$ in the left half-plane: for $t<T$ one cannot close that contour indiscriminately and discard its large arc. The modal forcing integral with upper limit $t$ avoids this causality error.

The expansion is singular near the moving lower limit, so expanding $e^{-\epsilon /x}$ throughout the original integral would lose its constant term. Introduce the symmetric exponential integralUnder $x=\epsilon /t$, the discarded piece becomesThis is a uniform Taylor expansion because $t\ge 1$. To derive the logarithmic expansion of a symmetric exponential integral without needing another special function, first exploit the same substitution to writeThe estimate follows from $0\le 1-e^{-\epsilon /x}\le \epsilon /x$. Here $\gamma$ is Euler's constant and the small-argument expansion of the exponential integral fixes the constant as well as the logarithm.

Differentiate under the integral and integrate the derivative of $e^{-x-\epsilon /x}/x$ over $(0,\infty )$; the exponentially small endpoint terms vanish. This gives the ordinary differential equationFor the Frobenius method at its regular singular point, write $J=A\log \epsilon +B+\epsilon (C_1\log \epsilon +D)+\cdots$. Substitution gives $C_1=A$ and $D=B-2A$. The preceding leading calculation sets $A=-1$, $B=-2\gamma$, henceSubtracting the moving-cutoff correction to a symmetric exponential integral yields the requested answer, including the logarithmically enhanced first-order term:Thus the $\epsilon \log \epsilon$ switchback term must be retained when asking for all contributions through order $\epsilon$.

Define the nested integralThe order-$s^2$ equation in the inner expansion is ${\mathcal{L}}_0{Y}_2=2z{Y}_1$. The preceding result decomposes its forcing into $2z^{2}E+2bzE-2bz+2bzg$. The particular solutions combine toFor example, ${\mathcal{L}}_{0}H=E$ follows directly by differentiating the nested integral, so the equation can also be verified without relying on the supplied particular-solution formulas. At zero $Y_{2,p}=b^2$; add $-b^2+DE$ to impose the left boundary condition. As $z\to \infty$ this givesThe outer expansion in the overlap, derived in the previous section, requires the remaining constant to be $-\frac{1}{2}\log s$. Thus $D=C+b^2-\frac{1}{2}\log s$. The logarithmic matching constant is an error-function logarithmic switchback.

The complete answer through order $\epsilon$ in the inner expansion isHere $E=\mathrm{erf}z$, $g=e^{-z^2}$ and $b=2/\sqrt{\pi }$. This expression is zero at $z=0$, and its large-$z$ expansion agrees with $1+x+x^2/2-(\epsilon /2)\log x$. Together with $y_{\mathrm{out}}=e^x[1-(\epsilon /2)\log x]$, it covers both requested regions. The overlap uses $1\ll z\ll {\epsilon }^{-1/2}$, with the usual termwise ordering restrictions when matching higher terms.

If a single continuous approximation is desired, an additive composite expansion isThe apparent logarithmic divergence at $x=0$ cancels: its residual is $-(\epsilon /2)(e^x-1)\log x\to 0$. The composite therefore extends continuously to the left endpoint and satisfies the right boundary condition to the requested order.

For smooth rapidly decaying initial data, apply the Fourier transform in $x$ and differentiate under the integral. The Fourier transform of a derivative turns the wave equation intoSolving this constant-coefficient ordinary differential equation givesAt $k=0$ the quotient means its continuous limit $t$, so $\tilde{v}(0,t)=\tilde{f}(0)+t\tilde{g}(0)$.

Since $\cos kt=(e^{ikt}+e^{-ikt})/2$, the Translation property of the Fourier transform inverts the first term as $(f(x+t)+f(x-t))/2$. Also $K_t(x)=\frac{1}{2}{1}_{[-t,t]}(x)$ has Fourier transform $\sin kt/k$, including its limiting value at zero. The convolution theorem therefore inverts the second term asThis proves the D'Alembert formulaThe formula itself extends beyond the rapidly decaying class: $f\in C^2$ and $g\in C^1$ suffice for a classical solution, verified by differentiation and the two initial conditions. The decay assumption above justifies the ordinary transform derivation without distributional qualifications.

Take a smooth variation $u_{\epsilon }=u+\epsilon \varphi$ with $\varphi =0$ on the boundary. Assume enough differentiability for the following operations, for example a smooth bounded open set, $u\in C^2$, and a twice continuously differentiable integrand. Differentiation under the integral givesBy integration by parts, the boundary term vanishes and the first variation isStationarity for every such variation, including every compactly supported smooth one, and the fundamental lemma of the calculus of variations give the Euler-Lagrange equationFor the quadratic integrand, $L_u=2k^{2}u$, $L_{u_x}=2u_x$, and $L_{u_y}=2u_y$. Thuswith the prescribed Dirichlet boundary condition. For real $k$, a stationary function is also the unique minimizer, if it exists: for any nonzero admissible $\varphi$, the linear term vanishes and $I[u+\varphi ]-I[u]={\int }_{\mathcal{D}}(|\nabla \varphi {|}^2+k^2{\varphi }^2)>0$. This last positivity uses the zero boundary values; it also holds at $k=0$.

Work on $([0,1],\mathcal{B},dt)$ as a probability space, with Lebesgue measure of total mass one. Let ${\mathcal{F}}_n$ be the filtration generated by the level-$n$ half-open dyadic cells together with the separate null cell $\{1\}$. Define the dyadic slope martingale byThe endpoint $t=1$ may be assigned any value, since it is a null set for Lebesgue measure. The two child slopes average to their parent slope, by telescoping the two increments. Therefore $\mathbb{E}[G_{n+1}\mid {\mathcal{F}}_n]=G_n$ almost everywhere, so this is a martingale. The Lipschitz condition gives $|G_n|\le K$ everywhere except possibly at the freely chosen endpoint, where we take zero.

Apply the Lp martingale convergence theorem with $p=2$. Its limit $g$ has $|g|\le K$ almost everywhere and $G_n\to g$ in L1 norm as well, by Cauchy-Schwarz inequality. Choose a measurable representative of $g$ and set it to zero on any exceptional null set; it is then a bounded measurable function on the entire interval.

Set $f_n(x)=f(0)+{\int }_0^x{G}_n(t)dt$. Telescoping at the grid points shows that $f_n$ is the linear interpolation of $f$ on the dyadic grid. The Lipschitz condition gives $\Vert f_n-f{\Vert }_{\infty }\le K{2}^{-n}\to 0$. AlsoThe two uniform limits coincide, giving the absolutely continuous function representationThe chosen bound holds for every $t$ after the null-set modification; the integral identity holds for every $x$ simultaneously.

Use the same dyadic slope martingale $G_n$ and dyadic filtration as in (c). Each $G_n$ is integrable, since it takes finitely many finite values. On each dyadic cell, the absolute slope is $2^n$ times the absolute endpoint increment. Consequently the hypothesis in the PDF is exactlyThus $(G_n)$ is uniformly integrable. It is also bounded in L1 norm: choose a finite ${\lambda }_0>0$ at which the supremum of the tails is finite, and use $\Vert G_n{\Vert }_1\le {\lambda }_0+{\sup }_j{V}_j(f,{\lambda }_0)$. The uniformly integrable martingale convergence theorem supplies $g\in L^1([0,1])$ with $G_n\to g$ in L1 norm.

The functions $f_n(x)=f(0)+{\int }_0^x{G}_n(t)dt$ are again the dyadic linear interpolations of $f$. Since $f$ is continuous on a compact interval, it is uniformly continuous, and $\Vert f_n-f{\Vert }_{\infty }\le {\omega }_f(2^{-n})\to 0$, where ${\omega }_f$ is its modulus of continuity. On the other hand, the integral of $G_n-g$ is uniformly bounded in absolute value by $\Vert G_n-g{\Vert }_1\to 0$. Hence the dyadic slope-tail criterion for absolute continuity givesNo boundedness of $g$ is asserted here; the tail condition permits integrable densities that are unbounded.

Adaptive rejection sampling is exact rejection sampling for a differentiable log-concave probability density. Given interior evaluation points $s_j$, construct the upper envelopewhere $h$ is the unnormalized target probability density function. The concave function $\ell$ lies below all its tangent lines, so $e^u\ge h$. Interpolation between adjacent evaluation points gives a lower envelope $l\le \ell$ there, allowing a squeeze test. The normalized upper envelope $g=e^u/\int e^u$ is a piecewise exponential proposal distribution.

Draw $X\sim g$ and an independent $W$ from the uniform distribution, and accept when $W\le e^{\ell (X)-u(X)}$. A squeeze test $W\le e^{l(X)-u(X)}$ can establish acceptance without evaluating $\ell (X)$. Otherwise evaluate $\ell (X)$ and, whenever it is evaluated, add $X$ to the envelope's points. Repeat until acceptance. On each linear piece $u(x)=rx+s$, its integral is $e^s(e^{rd}-e^{rc})/r$ on $[c,d]$, with value $e^s(d-c)$ when $r=0$. These integrals select a piece; inverse transform sampling samples within it.

Here, with $a=y_4$ and $b=y_2+y_3$, the log-posterior on $(0,1)$ hasChoose evaluation points strictly inside $(0,1)$, avoiding any endpoint singularities of the logarithm. Every finite set of tangent lines gives a finite upper-envelope integral on this bounded interval; slope conditions needed for unbounded supports are unnecessary. When $N=0$ the posterior distribution is simply the uniform distribution and may be sampled directly.

For any fixed envelope, proposal times acceptance is proportional to $h(x)$. Conditional on all past evaluations and rejections, the next accepted value therefore has the same normalized posterior density. To see that adaptation does not spoil this, at every attempted draw its accepted mass is a scalar times that fixed posterior density; summing over possible rejection histories preserves it. The acceptance probability stays bounded below by its positive initial-envelope value because insertion only lowers the upper envelope. Thus acceptance occurs almost surely. With fresh independent randomness, successive accepted values are exact independent samples from the posterior distribution, even when the envelope is retained and improved between them.

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.

Particle production explains the quantum input. For a real scalar field obeying the massless covariant wave equation, the Klein-Gordon inner product on solutions isThe wave equation makes its current divergence-free, so the product is independent of a Cauchy hypersurface when the boundary flux vanishes. Use suitably normalized wave packets or the appropriate continuum distributions. Asymptotic past and future Minkowski spacetimes give preferred positive-frequency solution mode bases $u_j^{\mathrm{in}}$ and $u_i^{\mathrm{out}}$, normalized by $(u_i,u_j)={\delta }_{ij}$, $(u_i^{\ast },u_j^{\ast })=-{\delta }_{ij}$ and $(u_i,u_j^{\ast })=0$. Between those regions, a nonstationary geometry generally has no preferred positive-frequency splitting.

The mode bases are related by a Bogoliubov transformation,The last identities follow from conservation of the Klein-Gordon inner product and are the canonical identities for a bosonic Bogoliubov transformation. Expanding the field in either complete basis and extracting its positive-frequency coefficient givesThus the state annihilated by every $a_j^{\mathrm{in}}$ hasThis particle number from Bogoliubov coefficients is nonzero when time evolution mixes positive and negative frequencies. It describes one state as vacuum in the early particle basis and populated in the late basis; it does not require an arbitrary choice of a vacuum at each intermediate time. For a finite implementable transformation the state remains a squeezed pure quantum state. In infinitely many modes, a common unitary bosonic Fock space implementation requires additional conditions, such as the beta map being a Hilbert-Schmidt operator; finite wave packet observables need not share every global divergence of an ideal continuum calculation.

For gravitational collapse, the future is not globally Minkowski: it contains a black hole. The relevant late out-basis includes modes reaching future null infinity together with modes entering the future event horizon. Modes on infinity alone are not a complete Cauchy basis. Nevertheless the same mode-mixing calculation determines the outgoing particle flux. Near a nonextremal horizon, the logarithm in the tortoise coordinate converts regular early null coordinates into the Hawking exponential ray mapwhere $U$ is the early affine null coordinate and $u$ is late retarded time at infinity. Backward propagation of a late outgoing mode $e^{-i\omega u}$ gives a factor $(U_H-U{)}^{i\omega /\kappa }$ for $U<U_H$. Its positive- and negative-frequency Fourier integrals differ by analytic continuation around the logarithmic branch. With a convergence regulator they reduce toTaking $ϵ\downarrow 0$ gives the thermal ratio of Hawking Bogoliubov coefficients, $|\beta {|}^2/|\alpha {|}^2=e^{-2\pi \omega /\kappa }$. Combining this with the canonical normalization yields the bosonic occupation $[e^{2\pi \omega /\kappa }-1{]}^{-1}$. Late-time wave packets turn formal continuum coefficients into a finite number flux.

A Schwarzschild black hole has $\kappa =1/(4M)$, so this spectrum has $T_H=1/(8\pi M)$. The exterior curvature potential partly reflects the outgoing modes. Its transmission probabilities are the greybody factors, giving, for one massless scalar species in the stationary late-time approximation,Thus the horizon temperature is universal, while the spectrum received at infinity is not a perfect featureless blackbody radiation spectrum. The derivation assumes the near-horizon quantum state inherited from a regular collapse vacuum and ignores rapid backreaction over the timescale of the packets. Exponentially blueshifted precursor frequencies indicate the usual short-distance assumption in the semiclassical calculation, not an independently demonstrated quantum-gravity description of the endpoint.

The outgoing positive energy flux drives black-hole evaporation. At leading order for a large isolated Schwarzschild black hole, area scales as $M^2$ and temperature as $M^{-1}$, giving $dM/du\sim -\gamma /M^2$ and a lifetime of order $M^3$, with the coefficient depending on species and transmission factors. Its negative heat capacity of a Schwarzschild black hole, $dM/d{T}_H=-8\pi M^2$, means that it heats up as it loses mass and cannot be in stable canonical equilibrium with an unlimited thermal reservoir. The approximation fails when quantum-gravitational scales are reached and does not fix whether the endpoint is complete evaporation, a remnant, or something else.

Area loss does not contradict the classical area theorem, because the quantum stress-energy tensor need not obey the classical null energy condition; negative horizon energy accompanies the positive outgoing flux. The appropriate thermodynamic statement is the generalized second law, involving $S_{\mathrm{BH}}+S_{\mathrm{outside}}$. Finally, Hawking radiation raises the black hole information paradox: if a pure quantum state completely evaporates into an exactly thermal mixed quantum state and its partners disappear, the result conflicts with ordinary unitary time evolution. Early outgoing radiation can be mixed simply because of its entanglement with interior modes while the complete state remains a pure quantum state; that fact alone is not information loss. Resolving the fate of all correlations through the evaporation endpoint requires more than the leading semiclassical flux calculation.

Put $u=\lambda v$ and let $n$ now mean number per area. During time $T$, the two strips of impact parameters between $b$ and $b+db$, on opposite sides of the trajectory, have combined area $2uTdb$. ThusThe same straight-line impulse approximation gives a kick $2Gm/(bu)$. Summing its squared size over these encounters givesassuming $b_{\max }\gg b_{\min }$. The relevant change is now the small random speed $u$, not the full circular speed $v$. Setting $D_v{t}_{\mathrm{rel}}\sim u^2$ yieldsThe razor-thin disk relaxation integral is dominated by nearby encounters, rather than equally by logarithmic intervals. Because kicks cease to be weak near the cutoff, its numerical coefficient is an estimate.