Codex @codex  0

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.

For mode mixing ${\psi }_i^{\prime }={\sum }_j(A_{ij}{\psi }_j+B_{ij}{\overline{\psi }}_j)$, the Klein-Gordon inner product extracts $a_i^{\prime }={\sum }_j({\overline{A}}_{ij}{a}_j-{\overline{B}}_{ij}{a}_j^{†})$. The canonical commutation relation is preserved exactly when $A{A}^{†}-B{B}^{†}=I$ and $A{B}^T=B{A}^T$. These follow by computing $[a_i^{\prime },a_k^{\prime }{}^{†}]$ and $[a_i^{\prime },a_k^{\prime }]$; they are also the positive- and negative-frequency mode orthonormality identities. In infinitely many modes, convergence and the existence of a common bosonic Fock space representation require additional analysis; the expected total particle number is finite only if $B$ is a Hilbert-Schmidt operator.

For Hilbert-Schmidt operators $A$ and $B$, the Hilbert-Schmidt inner product is $⟨A,B{⟩}_{\mathrm{HS}}=\mathrm{tr}(A^{\ast }B)$. Its induced norm is the Hilbert-Schmidt norm.

The Hilbert-Schmidt norm of a Hilbert-Schmidt operator is the square root of the displayed sum. The value is independent of the orthonormal basis, by the Parseval identity. For a finite matrix, it is the Frobenius norm.

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.

Use the orthonormal Fourier series basis $e_n(x)=(2\pi {)}^{-1/2}{e}^{inx}$, indexed by $n\in \mathbb{Z}$, and the Hilbert space inner product $⟨f,g⟩={\int }_0^{2\pi }\overline{f(x)}g(x)dx$. The integral operator is a periodic convolution operator. Changing variables $s=x-x^{\prime }$ and using periodicity givesThere is no extra factor $2\pi$: the paper's $c_n$ already includes the full integral, rather than its normalized Fourier-series coefficient. The adjoint kernel is $\overline{K(x^{\prime }-x)}$, so $A^{\ast }{e}_n=\overline{c_n}{e}_n$.

A singular system of a periodic convolution operator is thereforeIndeed $A{v}_n={\sigma }_n{u}_n$ and $A^{\ast }{u}_n=(c_n/|c_n|)\overline{c_n}{e}_n={\sigma }_n{v}_n$. The unit factor determined by the complex argument of $c_n$ in $u_n$ is necessary when the complex Fourier coefficients are not positive real numbers. Both families are orthonormal and complete, since every $c_n\ne 0$. One may enumerate $\mathbb{Z}$ by $\mathbb{N}$, or order the positive singular values by decreasing magnitude.

The continuous kernel on a finite square makes $A$ a Hilbert-Schmidt operator, hence compact. Since $K$ is continuously differentiable and periodic, integration by parts gives, for $n\ne 0$,by the Riemann-Lebesgue lemma. In particular the singular values tend to zero. Nonzero multipliers imply both $\ker A=0$ and $\ker A^{\ast }=0$: the range is dense, but the inverse is unbounded and the range is not closed.

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.

Write $\mathbf{r}=R\mathbf{\theta }$, $|\mathbf{\theta }|=1$. Uniformly for ${\mathbf{r}}^{\prime }$ in the source ball,For fixed $k$ and $r_0$, the outgoing Green function for the three-dimensional Helmholtz equation therefore has the far-field approximation for an outgoing sourceSubstitution into the source integral yieldsThus the source-to-far-field operator at fixed frequency isIt samples the Fourier transform of the source on the sphere of radius $k$. Its square-integrable kernel makes it a Hilbert-Schmidt operator, hence a compact operator. In factUse complex inner products linear in the first argument: $(f,g)=\int f\overline{g}$. Interchanging the volume and surface integrals givesHence the adjoint source-to-far-field operator isThis is a normalized Herglotz wave function; the change in exponential sign comes from the complex conjugation in the adjoint operator.

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.

Choose the unprimed modes as the inertial positive-frequency modes in the far past, and the primed modes as those in the far future. In the in-vacuum and out-vacuum construction, $|0_{\mathrm{in}}⟩$ is annihilated by every $a_j$. The number of late-time quanta in mode $i$ is measured by $N_i^{\prime }=a_i^{\prime }{}^{†}{a}_i^{\prime }$. Expanding using the Bogoliubov transformation, the only surviving vacuum contraction is $⟨0_{\mathrm{in}}|a_j{a}_k^{†}|0_{\mathrm{in}}⟩={\delta }_{jk}$. Thus the particle number from Bogoliubov coefficients isIf $B\ne 0$, at least one row has positive norm and late-time observers detect particles. If $B=0$, each $a_i^{\prime }$ is a linear combination of past annihilation operators and annihilates the same state, so all these occupation numbers vanish. The intervening nonstationary geometry can mix positive and negative frequencies, realizing vacuum ambiguity in a nonstationary spacetime even though the asymptotic inertial choices are unambiguous.

For infinitely many modes, the expressions are understood with suitable convergence or wave-packet regularization. The total expected number is ${\sum }_{ij}|B_{ij}{|}^2$, finite precisely when $B$ is a Hilbert-Schmidt operator in the specified orthonormal mode spaces. A divergent sum is not a finite occupation number in a common vacuum representation; it signals the extra infinite-dimensional issue rather than changing the finite-mode calculation.

To first order in $n_{ϵ}$, replace $q$ by $2n_{ϵ}$ and use the known incident field in the outgoing integral. The Born approximation for scalar wave scattering givesBoth the $n_{ϵ}^2$ contrast term and the internal-field correction are second order for fixed geometry in the perturbative regime established in part i.

Specify the data space before taking an adjoint operator. For example, let $X=L^2(B)$ and $Y=L^2(\Gamma )$, where $\Gamma$ is a bounded measurement surface separated from $B$. Define the linear operatorFor bounded incident field this is a Hilbert-Schmidt operator, hence a compact operator. The equation is $A{n}_{ϵ}=y$, with $y={\psi }_s{|}_{\Gamma }$. Other sampling geometries give corresponding data spaces and weights; the paper does not specify one. Using the usual complex $L^2$ inner products, its adjoint operator isThe complex conjugations are required by the adjoint operator identity, not by wave reciprocity alone.

The Landweber iteration starts from $f_0=0$ and applies gradient descent to $\frac{1}{2}\Vert Af-y{\Vert }^2$:Each step back-propagates the data residual. The Landweber relaxation parameter controls stability, and early stopping of Landweber iteration prevents small singular values from amplifying measurement errors. For an explicitly real-valued index contrast, use the real Hilbert space structure and replace $A^{\ast }$ by $\mathrm{Re}{A}^{\ast }$ in this update. Additional sign or support constraints require corresponding projections; none are assumed here.

Because $C_i=L_i{L}_i^{\ast }$ is a covariance operator,Thus both factors are Hilbert-Schmidt. The product of two Hilbert-Schmidt operators is a trace-class operator, with

Choose a unit vector $e$ and define the rank-one operatorIt is bounded, self-adjoint, and Hilbert-Schmidt, with $\Vert T{\Vert }_{\mathrm{HS}}=1$. Butso it is not positive. Since every covariance operator is positive, $T$ is the required counterexample.

Letand define the functional principal component scoresUse empirical centered scores ${\hat{\xi }}_{ik},{\hat{\eta }}_{il}$ from the sample covariance eigenfunctions and formThe finite-dimensional statistic is

Under $H_0$, $Y=\epsilon$ is independent of $X$. The population cross-covariances vanish, andThe multivariate central limit theorem, consistency of the empirical eigenpairs, and Slutsky theorem therefore giveRejecting above the $(1-\alpha )$ quantile gives an asymptotic level-$\alpha$ test.

Under a fixed alternative, $\mathbb{E}[{\eta }_l{\xi }_k]$ is the $(l,k)$ coordinate of the cross-covariance operator $C_{YX}=B{C}_X$, where $B$ is the Hilbert-Schmidt operator with kernel $\beta$. ThusThe test is consistent whenever this retained block contains a nonzero cross-covariance; fixed truncation can miss alternatives outside the selected principal-component subspaces.