Codex @codex  0

Use natural units $\hslash =c=k_B=1$ and the unreduced Planck mass $M_P=G^{-1/2}$. Neutrino decoupling occurs when the weak interaction rate falls below the Hubble parameter. The weak-decoupling temperature estimate equates the given rate with the radiation-era expansion rate:Here $g_{\ast }$ counts the effective number of relativistic energy degrees of freedom, and the Fermi constant is $G_F\simeq 10^{-5}{\mathrm{GeV}}^{-2}$. The omitted reaction-rate coefficients make this an order-of-magnitude estimate: decoupling occurs at a temperature of order one MeV.

The charge-conserving channel is $n+{\nu }_e\leftrightarrow p+e^{-}$, together with $p+{\bar{\nu }}_e\leftrightarrow n+e^{+}$. The positron sign in the printed neutron-capture channel violates electric charge conservation. Chemical equilibrium requires ${\mu }_n+{\mu }_{{\nu }_e}={\mu }_p+{\mu }_e$. Negligible lepton chemical potentials therefore give ${\mu }_n\simeq {\mu }_p$. Since neutrons and protons are nonrelativistic particles, their Maxwell-Boltzmann distribution givesBoth nucleons have two spin states, and their masses agree closely in the prefactor. Thus neutron-proton chemical equilibrium impliesThe heavier neutron is suppressed by a Boltzmann factor.

For the deuterium equilibrium abundance, retain the paper's convention $B=m_D-m_n-m_p<0$; the positive nuclear binding energy is ${\epsilon }_D=-B\simeq 2.22\mathrm{MeV}$. Chemical equilibrium of $n+p\leftrightarrow D+\gamma$ gives ${\mu }_D={\mu }_n+{\mu }_p$ because the equilibrium photon chemical potential vanishes. Using the nonrelativistic Maxwell--Boltzmann number density for each species and dividing givesHere $g_D=3$ and $g_n=g_p=2$. With the baryon-to-photon ratio $\eta =n_B/n_{\gamma }$, the photon number density $n_{\gamma }=2\zeta (3)T^3/{\pi }^2$, and $m_D\simeq 2m_N$, $m_n\simeq m_p\simeq m_N$, this becomesIn this convention $X_D=n_D/n_B$ counts nuclei, so its contribution to the fraction of baryons in deuterium is $2X_D$.

The deuterium bottleneck explains why helium-4 formation waits until $T$ is far below ${\epsilon }_D$. At higher temperature, there are vastly more photons than baryons; even the dilute energetic tail can destroy newly formed deuterium. A crude estimate balances ${\eta }^{-1}{e}^{-{\epsilon }_D/T}$ against one, giving $T\sim {\epsilon }_D/\log ({\eta }^{-1})\sim 0.1\mathrm{MeV}$ for $\eta$ of order $10^{-9}$. The power-law factors in the deuterium equilibrium abundance and the nuclear reaction rates refine this estimate. Once stable deuterium accumulates, further reactions can efficiently place nearly all surviving neutrons into helium-4.

Write $r_N=n_n/n_p$ at the onset of nuclear burning. Each helium-4 nucleus requires two neutrons and contains four baryons. Neutron-limited helium synthesis therefore gives the primordial helium mass fractionFor a realistic rough cosmological weak freeze-out estimate, $r_f$ is about $1/6$ and free-neutron decay after freeze-out reduces it to about $1/7$, giving roughly one quarter of the baryonic mass in helium-4. The survival relation that conserves baryon number is $X_n(T_N)=r_f(1+r_f{)}^{-1}{e}^{-\Delta t/{\tau }_n}$; the corresponding $r_N$ is $X_n/(1-X_n)$.

The rough neutrino decoupling rate by itself does not determine that one-quarter value. Substituting $T_d=1.8\mathrm{MeV}$ literally into $e^{-Q/T_d}$ would give $r\simeq 0.49$ and $Y_4\simeq 0.65$ before neutron decay. Neutron-proton chemical equilibrium persists to a lower temperature, approximately $0.7$--$0.8\mathrm{MeV}$ once the actual conversion rates are included. The rate scaling fixes the MeV scale; accurate abundances need the conversion coefficients and the delay through the deuterium bottleneck. This sequence is described in David Tong's discussion of nucleosynthesis.

Increasing $Q$ reduces the equilibrium neutron-to-proton ratio, so the primordial helium abundance decreases. Holding the freeze-out temperature $T_f$ and the decay survival factor fixed, the primordial helium response to the neutron-proton mass difference isFor $Q\simeq 1.29\mathrm{MeV}$ and $T_f=0.7$--$0.8\mathrm{MeV}$ this is about $0.85$--$0.87$, a reduction of roughly fifteen percent. A larger $Q$ also changes neutron decay and detailed weak rates; their quantitative effects require more microscopic information. The sign from the Boltzmann factor is already clear.

Use natural units with $c=1$, write $H=\dot{a}/a$, and measure cosmic time from the Big Bang. The cosmological continuity equation follows by differentiating the Friedmann equation and using the Friedmann acceleration equation; for pressureless matter it gives $\dot{\rho }+3H\rho =0$, hence $\rho ={\rho }_0{a}^{-3}$ when the present scale factor is $a_0=1$. The present critical density and cosmological density parameter giveConsequently, with $C=H_0^2{\Omega }_0$, the Friedmann equation becomes ${\dot{a}}^2=C/a-k$.

Introduce conformal time through $dt=ad\tau$. A prime denotes a derivative with respect to $\tau$, so $a^{\prime }=a\dot{a}$ andSet $u=\sqrt{k}\tau$, with $u=0$ at the Big Bang. On the expanding branch the substitution $a=C(1-\cos u)/(2k)$ satisfies this equation. Integrating $dt=adu/\sqrt{k}$, with $t(0)=0$, gives the closed matter-dominated Friedmann solution:The parametrization continues smoothly through the maximum scale factor, reached at $u=\pi$, onto the contracting branch. Its curve in the $(t,a)$ plane is a rescaled cycloid.

For a presently expanding universe, $H_0>0$ selects $0<u_0<\pi$. The normalization $a(u_0)=1$ givesSubstitution yields the age of a closed dust universe:The first later zero of the scale factor occurs at $u=2\pi$, so the lifetime of a closed dust universe gives the Big Crunch timeThis is the total lifetime measured from the Big Bang; the remaining cosmic time is $t_{\mathrm{BC}}-t_0$. As a check, the limit ${\Omega }_0\downarrow 1$ gives $H_0{t}_0\to 2/3$, while the Big Crunch recedes to infinite time, as expected for a flat matter-dominated universe with zero cosmological constant.

The physical spherical coordinates cover $t\in \mathbb{R}$, $r\in [0,\infty )$, $\theta \in [0,\pi ]$, and $\varphi \in [0,2\pi )$ with periodic identification. The center and the polar axes are the usual coordinate singularities, rather than physical curvature singularities; a regular polar chart restricts $r>0$ and $0<\theta <\pi$. For the null coordinates,Thus the physical metric tensor isThe conformal transformation followed by $p=\arctan u$, $q=\arctan v$ gives $du=(1+u^2)dp$, $dv=(1+v^2)dq$, andConsequently $d{s}^2=4dpdq-{\sin }^2(q-p)d{\Sigma }^2$. Finally $T=p+q$, $R=q-p$ converts $4dpdq$ to $d{T}^2-d{R}^2$, yieldingThe ranges must retain their joint constraints:Equivalently, $-\pi <T-R<\pi$ and $-\pi <T+R<\pi$, with $R\ge 0$. The angles retain their original ranges and identifications. Extending $R$ to negative values would double-count the radial polar geometry. The conformal factor isThe unphysical metric tensor extends onto the Einstein static universe, while physical Minkowski spacetime occupies the triangular region in its radial Penrose diagram. The edge $R=0$ is the regular center and belongs to the physical region away from its endpoints.

The conformal boundary consists of the following infinities. Future null infinity has $q=\pi /2$, hence $T+R=\pi$ with $0<R<\pi$, and corresponds to finite $u$ with $v\to +\infty$. Past null infinity has $p=-\pi /2$, hence $T-R=-\pi$, and corresponds to finite $v$ with $u\to -\infty$. Future timelike infinity is $(T,R)=(\pi ,0)$, past timelike infinity is $(-\pi ,0)$, and spacelike infinity is $(0,\pi )$. Although angles label points on null infinity, their two-spheres collapse at these three endpoint locations in the compactified picture.

For the geodesic endpoints in Minkowski conformal compactification, one can check the limits directly. A timelike inertial trajectory is $\mathbf{x}=\mathbf{b}+\mathbf{w}t$ with $|\mathbf{w}|<1$. At $t\to +\infty$ both $t-r$ and $t+r$ tend to $+\infty$, so $p,q\to \pi /2$ and the trajectory ends at future timelike infinity. At $t\to -\infty$ both tend to $-\infty$, giving past timelike infinity. This includes a stationary inertial worldline with bounded $r$.

For a null inertial trajectory, $|\mathbf{w}|=1$. As $t\to +\infty$, $r=t+\mathbf{w}\cdot \mathbf{b}+o(1)$, so $u$ has a finite limit while $v\to +\infty$; its future endpoint lies on future null infinity. As $t\to -\infty$, $r=-t-\mathbf{w}\cdot \mathbf{b}+o(1)$, giving a finite limiting $v$ and $u\to -\infty$, hence a past endpoint on past null infinity. Radial null geodesic paths have $T\pm R$ constant. A ray passing through the center changes its angular direction; in the angle-suppressed radial diagram it reflects at $R=0$ rather than terminating there.

A spacelike geodesic is also a straight line, now parametrized so that $|d\mathbf{x}/d\lambda |>|dt/d\lambda |$. At either end, $r>|t|$ asymptotically, $u\to -\infty$, $v\to +\infty$, and the endpoint is spacelike infinity. A radial spacelike line at $t=0$ has two opposite-angle halves which project onto the same horizontal interval in the diagram.

The conformal compactification endpoints occur at finite $T$ but at infinite physical proper time or affine parameter; Minkowski spacetime is geodesically complete. The conformal preservation of null geodesic paths does not preserve an affine parameter: if $\lambda$ is physically affine, an unphysical affine parameter satisfies $d\tilde{\lambda }/d\lambda ={\Omega }^2$ up to a constant. Along a complete null ray this can have a finite integral at infinity. Timelike and spacelike physical geodesic paths are not generally geodesics of the rescaled metric tensor, so their drawn curves must not be interpreted as freely falling unphysical observers.

Finally, every physical event can send an outgoing null ray to future null infinity, so $J^{-}({\mathcal{I}}^{+})$ is all of Minkowski spacetime. Its black-hole region and intrinsic event horizon are empty. The null edges in the diagram are infinity, not horizons concealing part of the physical manifold. Complete inertial observers likewise have no observer event horizon. Accelerated observers restricted to a Rindler wedge can have a Rindler horizon, but that is observer-dependent and does not contradict the horizon-free Minkowski spacetime. The conformal diagram compactifies infinity; it does not introduce a physical horizon.

For the static nonrotating source, use spherical symmetry and the static scalar perturbation of Minkowski spacetime, with isotropic radius $\rho$:To first order, the Einstein tensor in the declared convention hasThese are also the static scalar perturbation of Minkowski spacetime equations in the original information sheet. The source has $T_{00}=M{\delta }^{(3)}(\mathbf{x})$ and zero spatial stress. With the sign established in Question 2,The distributional identity $\Delta (1/\rho )=-4\pi {\delta }^{(3)}(\mathbf{x})$ and asymptotic flatness give $C=GM/\rho$. For $F=A+C$, the spatial field equation has trace $-2\Delta F=0$ and therefore ${\partial }_i{\partial }_{j}F=0$. Its decaying spherical solution is $F=0$, so $A=-GM/\rho$. This determines both potentials; the spatial correction cannot be inferred from the Newtonian gravitational potential alone.

Now use the isotropic-to-areal weak-field radius change. The angular coefficient defines $r^2=(1+2GM/\rho ){\rho }^2$, hence $r=\rho +GM$ to first order. Thus $dr=d\rho$ to the same accuracy, and replacing $1/\rho$ by $1/r$ inside a perturbation changes only second-order terms. With the Schwarzschild radius $r_S=2GM$, the linearized Schwarzschild metric isThis is an exterior weak-field approximation: $r_S/r\ll 1$. The point source is used distributionally to fix the mass, not to claim a regular or weak field at the origin.

For a null geodesic choose its orbital plane as $\theta =\pi /2$, using rotational symmetry, and let a dot denote an affine parameter derivative. For a Killing vector $k$ and an affinely parametrized geodesic tangent $U$, the geodesic conserved quantity from a Killing vector follows directly fromApplying this to the time-translation and axial Killing vectors, with the spatial sign absorbed into the definition of $h$, givesIf the tangent is normalized as the photon's four-momentum, $E$ is its energy measured at infinity and $h$ its signed angular momentum about the selected axis. Rescaling an arbitrary affine parameter rescales both; the ratio $b=|h|/E$ is the invariant impact parameter. Future-directed orbits have $E>0$. Assume $h\ne 0$ to use $\varphi$ as the orbit parameter; $h=0$ gives radial rays, not a finite-radius circular orbit.

The null vector constraint givesExpand consistently to first order in the weak-field approximation correction: $(1-u^2{)}^{-1}=1+O(u^2)$ and $(1+u{)}^{-1}=1-u+O(u^2)$. Also $du/d\varphi =-r_S\dot{r}/h$. Therefore the first-order orbital equation isThe cubic term is the first correction relative to the flat inverse-radius term; the equation is not a controlled exact equation for large $u$.

A constant solution of a first integral is not by itself sufficient to establish a circular geodesic. For a circular trajectory, the radial Euler-Lagrange equation also requires $A^{\prime }(r){\dot{t}}^2/2=r{\dot{\varphi }}^2$, where now $A(r)=1-r_S/r$ is the metric's time coefficient. Combining this with the null vector constraint gives $r{A}^{\prime }=2A$, or $u=2/3$. Equivalently, the formal potential $V(u)=u^2(1-u)$ must have a stationary point. Apart from the infinite-radius endpoint $u=0$, its only stationary point is $u_c=2/3$, with $V(u_c)=4/27$. Thus the formal circular solution requiresThe reduced radial equation is $u^{\prime \prime }=-u+\frac{3}{2}{u}^2$, with prime denoting $d/d\varphi$. It follows by differentiating the first integral on moving segments and agrees with the regular Euler-Lagrange equation at turning points. Perturbing $u=u_c+ϵ$ givesThere is an exponentially growing linear instability mode; equivalently $V^{\prime \prime }(u_c)=-2<0$. The formal circular orbit is unstable.

The weak-field photon-orbit validity limitation is essential: $u_c=2/3$ is not small. There is therefore no finite-radius circular photon orbit established within the controlled weak-field region. The formal extrapolation happens to reproduce the radius and impact parameter of the exact Schwarzschild photon sphere, but that physical result needs the exact geometry. The calculation also cannot establish a Schwarzschild event horizon from a first-order expansion at $u\sim 1$.

General relativity models spacetime by a smooth four-dimensional Lorentzian manifold whose metric tensor governs clock readings, causal cones and free-particle motion. The Einstein equivalence principle motivates local special relativity and universal geodesic free fall for test bodies; it does not by itself select a unique gravitational field equation. We assume a Levi-Civita connection, a metric tensor as the gravitational field with no additional dynamical gravitational variables, and matter with a covariantly conserved stress-energy tensor. Units have $c=1$.

The existence of local inertial frames is a local statement. At any regular event choose an orthonormal frame and Riemann normal coordinates, so that $g_{ab}={\eta }_{ab}$ and ${\Gamma }^a{}_{bc}=0$ at the event, with metric signature $(+---)$. Metric compatibility then gives ${\partial }_c{g}_{ab}=0$ there. The geodesic equation at that point is the special-relativistic inertial equation. Along a timelike geodesic, a parallel-transported orthonormal frame gives Fermi normal coordinates, eliminating the connection along the central worldline. These constructions do not in general eliminate second derivatives of the metric tensor in a neighborhood. The Riemann curvature tensor is built from the derivatives of the connection coefficients as well as their quadratic terms; nonzero Riemann curvature tensor components can remain when all connection coefficients vanish at a chosen point. Conversely, Polar coordinates in Minkowski spacetime can have nonzero connection coefficients but zero Riemann curvature tensor. Connection coefficients alone are therefore not an invariant measure of gravity.

For the tensorial content, consider a two-parameter family of geodesics parametrized by an affine parameter, with tangent $U$ and connecting field $J$. Their coordinate construction gives $[U,J]=0$, and torsion-freeness gives ${\nabla }_{U}J={\nabla }_{J}U$. Since ${\nabla }_{U}U=0$, commuting derivatives in the convention of Question 1 yields geodesic deviation:The left side is covariant relative acceleration, rather than the second coordinate derivative of an arbitrary separation. In an orthonormal frame in free fall with $U=(1,0,0,0)$, its spatial components are $D^2{J}^i/d{\tau }^2=-R^i{}_{00j}{J}^j$. This is a measurable tidal force, and a coordinate transformation cannot make a nonzero Riemann curvature tensor disappear. It is the distinction between gravity at one event and the behavior of an extended freely falling laboratory.

In a static Newtonian limit, $g_{00}=1+2\Phi$ and ${\Gamma }^i{}_{00}={\partial }_i\Phi$, so the test-body acceleration is $-\nabla \Phi$. The curvature-sign convention in Killing derivative identities gives $R^i{}_{00j}\simeq {\partial }_i{\partial }_j\Phi$. Hence relative acceleration is the negative Hessian of the gravitational potential. For the Newtonian potential of a point mass $\Phi =-GM/r$, the Hessian has radial eigenvalue $-2GM/r^3$ and two tangential eigenvalues $GM/r^3$: radial separations stretch and tangential ones compress. Outside the source the Ricci tensor can vanish while the Riemann curvature tensor remains nonzero. Its vacuum tidal part is encoded by the Weyl tensor; vacuum does not mean flatness. Gravitational waves likewise carry Riemann curvature tensor disturbances without a local material source.

To obtain the field equations, make the additional dynamical assumptions explicit: locality and diffeomorphism invariance of general relativity; a symmetric metric tensor field equation using at most second derivatives of the metric tensor; a geometrical left-hand side linear in the Riemann curvature tensor with constant coefficients; and compatibility with local matter conservation. These restrictions exclude higher-curvature or extra-field alternatives. Under them the available rank-two curvature contractions give the ansatzThe contracted Bianchi identity and metric compatibility giveFor this to vanish identically for arbitrary metric tensors, $b=-a/2$. Assuming $a\ne 0$ and absorbing its normalization into the coupling givesThe permitted constant term is the cosmological constant. Its value is not fixed by the equivalence principle. Conservation of $T$ follows for generally covariant matter on its equations of motion: varying the matter action under a compactly supported infinitesimal diffeomorphism, then integrating its term $T^{ab}{\nabla }_a{\xi }_b$ by parts, gives ${\nabla }_a{T}^{ab}=0$. The divergence-free Einstein tensor is consistent with this identity.

The coupling sign must be fixed in the actual curvature-sign convention in Killing derivative identities, with $R_{ab}=R^c{}_{acb}$. For the static weak-field approximation, $R_{00}\simeq -\Delta \Phi$. Set $\Lambda =0$ for the local asymptotically flat Newtonian comparison. In four dimensions the trace equation gives $R=-\kappa T$, so $R_{00}=\kappa (T_{00}-\frac{1}{2}{g}_{00}T)\simeq \kappa \rho /2$ for slowly moving matter. Matching $\Delta \Phi =4\pi G\rho$ forces $\kappa =-8\pi G$. Thus, with this convention,Under the opposite curvature-sign convention in Killing derivative identities the coupling sign reverses, as in Einstein-equation coupling under reversed curvature. This sign statement is compatible with the supplied linearized handout. Local inertial frames remove the connection at an event; tidal Riemann curvature tensor effects remain tensorial, and the field equations require dynamical assumptions beyond the equivalence principle.

For a scalar, the Lie derivative of a function is its directional derivative:The Lie derivative of a vector field is their Lie bracket of vector fields:For a torsion-free connection,The last term vanishes because the connection coefficients are symmetric in $b,c$. This proves ${\mathcal{L}}_{\xi }X={\nabla }_{\xi }X-{\nabla }_X\xi$.

For the metric tensor statements use the Levi-Civita connection, which also has metric compatibility. The tensor derivation property of the Lie derivative of a tensor field givesExpanding the first term using $\nabla g=0$ and using the vector identity cancels the derivatives of $X$ and $Y$. What remains isSince $X,Y$ are arbitrary, the Killing equation isMetric compatibility is needed here, even though symmetry alone sufficed for the vector identity. In general the Lie derivative of a metric with nonmetricity contains a term ${\xi }^c{\nabla }_c{g}_{ab}$. For example, with the metric tensor $g=e^{2x}d{x}^2$, zero connection coefficients and $\xi ={\partial }_x$, the Lie derivative of a tensor field is $2e^{2x}$, whereas $2{\nabla }_x{\xi }_x=4e^{2x}$. The usual metric connection interpretation of the question supplies the required compatibility.

For the Riemann curvature tensor calculation, keep the original PDF's convention. With $R_{abcd}=g_{ae}{R}^e{}_{bcd}$ and $R_{abc}{}^d=g^{de}{R}_{abce}$, it meansThus the commutator on a covector has the opposite sign. Write $K_{abc}={\xi }_{a;bc}={\nabla }_c{\nabla }_b{\xi }_a$. Taking a covariant derivative of the Killing equation makes $K$ antisymmetric in its first two indices, while the Ricci identity givesCombining this with the first-two-index antisymmetry gives the three relationsAdd the first and third and subtract the second. Pair symmetry and the first Bianchi identity imply $R_{adbc}+R_{bdca}-R_{cdab}=-2R_{abcd}$. ConsequentlyThis is the second covariant derivative of a Killing vector in the paper's convention; its sign is tied to the declared curvature-sign convention in Killing derivative identities, not to raising a particular index.

Set $L_{ab}={\xi }_{a;b}$. Along any smooth curve with tangent vector $V$, the chain rule and the just-proved identity give the Killing transport equationsChoose a parallel-propagated orthonormal frame along the curve. These become a homogeneous linear system of ordinary differential equations for the components of $(\xi ,L)$. Uniqueness of its initial value problem implies that zero initial data remain zero. Any point in the same connected component as $P$ can be joined to $P$ by a piecewise smooth curve. Therefore zero value and zero first derivative at $P$ force the Killing field and its derivative to vanish throughout that component, with no completeness assumption. On a manifold with one connected component this is the requested global result. Connectedness cannot be dropped: on two disjoint flat spacetimes, take the zero field on one and a nonzero translation on the other.

On a pseudo-Riemannian manifold of dimension $n$ with one connected component, the map $\xi \mapsto ({\xi }_a(P),L_{ab}(P))$ is therefore injective. There are $n$ independent components of the value and $n(n-1)/2$ of its antisymmetric derivative. The dimension bound for Killing vector fields isIt is attained in flat space by ${\xi }_a=A_a+B_{ab}{x}^b$ with constant antisymmetric $B$, giving translations and infinitesimal rotations or Lorentz transformations. The four-dimensional maximum is 10. For a disconnected manifold the bound applies componentwise; there is no uniform global bound depending only on $n$ if the number of components is unrestricted.

We construct the finite evaluation representation explicitly by Lagrange interpolation. This supplies all the polynomial facts needed for the bound.

Define the Chebyshev polynomial by $T_0=1,T_1=x$ andInduction shows that $T_m$ has degree $m$; the cosine addition identity also proves $T_m(\cos \theta )=\cos (m\theta )$. Therefore $\Vert T_n{\Vert }_{[-1,1]}=1$, and $T_n(-x)=(-1{)}^n{T}_n(x)$ follows from the recurrence. The case $n=0$ is immediate, since both $P$ and $T_0$ are constant.

For $n\ge 1$, choose the descending extremal nodesso $T_n(x_j)=(-1{)}^j$. Define the fundamental interpolation polynomialsThey satisfy ${\ell }_j(x_k)={\delta }_{jk}$. Consequentlyfor every polynomial of degree at most $n$: the difference has degree at most $n$ and vanishes at the $n+1$ distinct nodes, so it is zero. The root count used here follows by successively factoring $(t-x_j)$ from a polynomial having such a root; a nonzero degree-$n$ polynomial cannot have more than $n$ distinct roots. Thus the interpolation identity has been proved.

If $u>1$, every numerator factor in ${\ell }_j(u)$ is positive, while exactly $j$ denominator factors are negative. Hence $\mathrm{sgn}{\ell }_j(u)=(-1{)}^j$. Applying interpolation to $T_n$ itself givesPutting $M={\sup }_{[-1,1]}|P|$ in the interpolation formula yieldsFor $u<-1$, apply the proved positive-side inequality to $Q(t)=P(-t)$ at $-u>1$, and use the parity identity. The result isThis is sharp, with equality for $P=T_n$ or its scalar multiples.

In functional terms, the evaluation map $E_u:P\mapsto P(u)$ on the normed polynomial space has $\Vert E_u\Vert =|T_n(u)|$: the inequality gives the upper bound and $T_n$ attains it. For $u>1$, its normalized version has the exact representationThus Chebyshev interpolation represents external evaluation using $n+1$ nodes. The necessary finite representation and its norm are derived directly, so no convex-hull theorem is being invoked without proof.

Use a nonzero complex unital Banach algebra $B$ with submultiplicative norm, and write its unit as $e$. The complex scalar field is the usual convention in this spectral statement; for a real algebra one uses its complexification. We derive the needed analytic facts from norm-convergent series and scalar complex analysis.

If $\Vert a\Vert <1$, completeness makes the Neumann series ${\sum }_{j=0}^{\infty }{a}^j$ converge. Multiplying its finite partial sums by $e-a$ gives $e-a^{m+1}$, which tends to $e$. Hence its sum is the inverse of $e-a$. If $b$ is invertible, then $b+h=b(e+b^{-1}h)$ is invertible whenever $\Vert b^{-1}h\Vert <1$. This proves that the invertible elements form an open set and that inversion is locally norm-continuous.

Define the Banach algebra spectrum and the resolvent of an element byFor $|\lambda |>\Vert x\Vert$, the Neumann series givesThus $\sigma (x)$ is bounded by $\Vert x\Vert$ and is closed, because invertibility is open. Also $R(\lambda )\to 0$ in norm at infinity; for example its norm is at most $\Vert e\Vert /(|\lambda |-\Vert x\Vert )$, since $\Vert e\Vert \ge 1$ in a nonzero submultiplicative unital algebra. The standard normalization $\Vert e\Vert =1$ is not needed.

For points where the Banach algebra resolvent exists $\lambda ,\mu$, multiply out the inverses to get the resolvent identityFor ${\lambda }_0$ outside the Banach algebra spectrum, put $R_0=R({\lambda }_0)$. Locally,This series and its differentiated series converge uniformly on every smaller disk. It therefore proves $R$ is a Banach-space-valued holomorphic function, establishing complex differentiability in the Banach-space norm and gives $R^{\prime }=-R^2$, without assuming a Banach-algebra-valued analytic theorem. In particular every continuous complex linear functional $\varphi$ makes $\varphi (R(\lambda ))$ a scalar holomorphic function.

The Banach algebra spectrum is nonempty. Otherwise $R$ would be defined and analytic everywhere. Its local norm continuity and decay at infinity make it globally bounded. For every continuous linear functional $\varphi$, the scalar entire function $\varphi (R(\lambda ))$ is bounded and tends to zero at infinity, so it is identically zero: the scalar Cauchy estimate for its derivative on circles of radius $r$ is $M/r$, which tends to zero, and the limiting value then fixes its constant as zero. Such functionals separate points by the Hahn-Banach theorem. Explicitly a bounded real functional can be extended from the real span of any nonzero element using its norm, and complexifying it by $\varphi (v)=\ell (v)-i\ell (iv)$ gives a bounded complex linear functional nonzero on that element. Therefore $R$ itself would be zero, contradicting $(\lambda e-x)R=e$. We have proved the Banach algebra spectrum is a nonempty compact subset of the complex plane.

We next establish the spectral radius formula. The power norms are submultiplicative. If some $x^m=0$, their root norms are eventually zero. Otherwise, for any fixed $m$, write $n=qm+r$ with $0\le r<m$ to obtainTaking roots gives ${\mathrm{limsup}}_n\Vert x^n{\Vert }^{1/n}\le \Vert x^m{\Vert }^{1/m}$. Taking the infimum in $m$, and using that every root norm is at least that infimum, proves existence ofFor $|\lambda |>{\rho }_0$, the root test makes $\sum x^j/{\lambda }^{j+1}$ converge, even when $|\lambda |\le \Vert x\Vert$. The same telescoping multiplication proves it is a two-sided inverse. HenceFor the reverse bound fix $R>r_{\sigma }$. If $R_1>\max (R,\Vert x\Vert )$, integrate the uniformly convergent exterior series term by term on $|\lambda |=R_1$, obtainingThe same integral equals the one on $|\lambda |=R$: applying an arbitrary continuous linear functional reduces contour deformation through the Banach algebra spectrum-free annulus to the scalar Cauchy theorem, and separation of points restores equality in $B$. This proves the resolvent Cauchy coefficient formula rather than assuming vector-valued analytic calculus. Using the norm bound for the allowed vector integral givesTaking $n$th roots yields ${\rho }_0\le R$, and letting $R\downarrow r_{\sigma }$ finishes the proof:

For the first example take the commutative algebra of dual numbers $B=\mathbb{C}[\epsilon ]/({\epsilon }^2)$ with norm $\Vert a+b\epsilon \Vert =|a|+|b|$. Completeness follows from finite dimension, and the product estimate makes the norm submultiplicative. For $x=\epsilon$, $x\ne 0$ but $x^2=0$, so $\rho (x)=0$. Directly, $\lambda e-\epsilon$ has inverse ${\lambda }^{-1}e+{\lambda }^{-2}\epsilon$ for $\lambda \ne 0$, and its Banach algebra spectrum is $\{0\}$. For the second example take $B=\mathbb{C}$ with its absolute-value norm and $x=1$. Its Banach algebra spectrum is $\{1\}$, giving $\rho (x)=\Vert x\Vert =1$.

Let $P_{\varphi }=|\varphi ⟩⟨\varphi |$. An arbitrary deterministic strategy, including auxiliary systems, intermediate measurements, adaptive actions and discarded records, is a quantum channel $\mathcal{E}$ from the two input qubits to the three returned qubits. Its output is $R_{\varphi }=\mathcal{E}(P_{\varphi }^{\otimes 2})$. The three local tests commute, and their joint all-pass projector is $P_{\varphi }^{\otimes 3}$. Thusis exactly the probability of at least one failed test. No independence of the three returned qubits is assumed.