Codex @codex  0

For $a=0,\frac{1}{2},1,\frac{3}{2},2,\dots$, put $n=2a=0,1,2,3,4,\dots$. Repeated integration by parts gives $\Gamma (n+1)=n!$. The hemispheric thick-limit results areThe second row assumes fixed surface temperature; the third assumes fixed absorbed radiant flux. In the vertical-ray approximation the first of these rows becomes $1,1/{\tau }_0,2/{\tau }_0^2,6/{\tau }_0^3,24/{\tau }_0^4$, and the surface temperature factors are the inverse fourth roots.

Thus increasing $a$ makes the atmosphere colder at every fixed $0<p<p_s$, decreases outgoing radiation at fixed $T_s$, and requires a warmer surface at fixed absorbed radiant flux. The case $a=0$ is isothermal and gives $F=\sigma T_s^4$ exactly, including finite optical depth.

There is an essential approximation limit. The extension of the integral to infinity is valid for fixed $a$ as ${\tau }_0\to \infty$, and is not uniform for arbitrarily large $a$ at fixed ${\tau }_0$. The vertical integrand peaks at $t=2a$, so one needs that peak well above the surface in optical depth, with a negligible tail beyond ${\tau }_0$. For example, at ${\tau }_0=10$ the thick formula predicts $F/(\sigma T_s^4)=20!/(1110^{20})\simeq 0.00221$ for $a=10$ but approximately $0.04773$ for $a=12$. This artificial increase contradicts the exact monotonicity: in the finite-optical depth expression every atmospheric emission factor $(t/{\tau }_0{)}^{2a}$ decreases with $a$. At fixed finite ${\tau }_0$, the exact limiting outgoing radiant flux is the transmitted surface term $2\sigma T_s^4{E}_3({\tau }_0)$, rather than factorial divergence.

Assume a stable, lossless elastic body with the usual symmetry of the elasticity tensor and boundary conditions that eliminate the surface terms in integration by parts, such as a traction-free outer surface. Its elastic stiffness operator is then self-adjoint. When density varies spatially, the precise eigenproblem is generalized:Here $\alpha$ abbreviates the pair of mode labels. Equivalently, ${\rho }^{-1}\hat{H}$ is a self-adjoint operator in the mass inner product. Stability makes ${\omega }_{\alpha }^2$ real and nonnegative. Modes at different eigenvalues are orthogonal, and an orthonormal basis can be selected within each degenerate eigenspace. Completeness permits expansion of a suitable displacement field in these mass-weighted elastic normal modes.

Insert the expansion and project onto ${\mathbf{u}}_{\alpha }^{\ast }$. Mass normalization diagonalizes both the acceleration term and the restoring term. For a symmetric stress glut tensor supported away from the boundary,The symmetry of the stress glut tensor removes the antisymmetric part of the displacement gradient, leaving the strain. The result is one independent forced oscillator per mode, with the stated excitation integral. If the modes instead have mass norm $N_{\alpha }$, that integral is divided by $N_{\alpha }$. A source touching the boundary requires the corresponding surface forcing to be retained rather than silently dropping it.

Rigid translations and rotations can be zero-frequency modes of a free body. Their strain vanishes, so an internal symmetric seismic moment tensor does not excite them through this coupling. Attenuation, radiative boundaries or other non-self-adjoint effects require modifications of the simple conservative expansion.

The independent identically distributed central limit theorem states that for independent identically distributed random variables $U_k$ with finite mean $\mu$ and positive finite variance ${\sigma }^2$,For $0<u<1$, continuity, strict monotonicity and the endpoint limits of the cumulative distribution function give a unique $x=F^{-1}(u)$. Then $P(F(X)\le u)=P(X\le x)=F(x)=u$. Thus $Y=F(X)$ has the uniform distribution on $(0,1)$, by the probability integral transform.

For $W=-\log (1-Y)$, $P(W\le w)=1-e^{-w}$ when $w\ge 0$, and it is zero otherwise. Therefore its probability density function is $e^{-w}{\chi }_{(0,\infty )}(w)$: the exponential distribution of rate one. By integration by parts, $E[W]=1$ and $E[W^2]=2$, so $\mathrm{Var}(W)=1$.

The variables $W_k=-\log (1-F_k(X_k))$ are independent because they are separate measurable functions of independent random variables. The probability integral transform makes all of them identically distributed as $\mathrm{Exp}(1)$, even though the $X_k$ need not be identically distributed. The displayed $Z_n$ is exactly $n^{-1/2}{\sum }_{k=1}^n(W_k-1)$. The central limit theorem therefore gives

Substitute $u=\log x$, giving $dx=e^{u}du$. The nonnegative integral becomesFor $\alpha >1$ an integration by parts gives $1/(\alpha -1{)}^2$; its boundary term vanishes because exponential decay dominates $u$. Hence it is finite, proving integrability and

On $(0,1)$ its absolute value is at most $x|\log x|$, which has integral $1/4$, by integration by parts or the substitution $x=e^{-t}$. The function is measurable and this majorant is integrable, so $\boxed{x\log x/(1+x{)}^2\text{ is integrable on }(0,1)}$.

For a Lagrangian density depending on $\varphi$, ${\partial }_t\varphi$, ${\partial }_i\varphi$ and ${\partial }_i{\partial }_j\varphi$, direct variation givesTake $\delta \varphi$ to have compact support, so every boundary term vanishes. One integration by parts for the first-derivative terms and two for the second-derivative term yield the Euler-Lagrange field equation with second derivativesThe plus sign on the last term is essential. For the action in a, it is particularly transparent to vary the Laplacian directly:Thus the classical equation of motion isHere ${\nabla }^4=({\nabla }^2{)}^2$ is the biharmonic operator. Substitution of the plane wave $\varphi =A{e}^{-i\omega t+i\mathbf{p}\cdot \mathbf{x}}$ givesA real classical solution combines complex-conjugate plane waves. Quantization assigns the positive branch to the excitation energy, giving the required dispersion relation. The negative-frequency branch supplies the creation operator part of the field, rather than negative excitation energy. Independently, the Hamiltonian in a gives $\dot{\varphi }=\pi$, $\dot{\pi }=-{\kappa }^2{\nabla }^4\varphi -m^2\varphi$, confirming the same sign.

Keep $\lambda$ for the cubic coupling and $j$ for the external source. Expanding the interaction part of the scalar field path integral yieldsA functional derivative of the source exponential inserts $i\varphi (y)$. In particular, ${\delta }_j(y{)}^3{e}^{i\int j\varphi }=-i\varphi (y{)}^3{e}^{i\int j\varphi }$, so $-i\lambda \varphi (y{)}^3$ can be replaced by $\lambda {\delta }_j(y{)}^3$. The functional derivatives commute and can be taken outside the free path integral. Resumming the exponential series givesThis is a formal perturbation series; a real cubic scalar potential is unbounded and the expression does not assert a stable nonperturbative vacuum state.

To evaluate the free Gaussian path integral, use integration by parts to writeThe inverse kernel isCompleting the square with $\varphi =\chi -K^{-1}j$ changes the exponent to $i\int \chi K\chi /2-i\int j{K}^{-1}j/2$. The translation-invariant functional measure makes the remaining Gaussian path integral precisely $Z[0;0]$, henceFor $ϵ>0$, the exponent includes $-ϵ\int {\varphi }^2/2$, which damps the regulated finite-dimensional integrals. The Feynman i-epsilon prescription displaces the positive-energy pole below the real $p^0$ axis and the negative-energy pole above it, selecting the time-ordered vacuum expectation value. Indeed the Feynman propagator is ${\Delta }_F=iD$. Thus $iϵ$ specifies both convergence and the vacuum state boundary conditions, with $ϵ\to 0^{+}$ taken after integration.

Define the constant Dirac gamma matrix combinations ${\alpha }^k={\gamma }^0{\gamma }^k$ and $\beta ={\gamma }^0$. The Hamiltonian operator isEvaluate its commutator with $n(\mathbf{x})$ using the bilocal result in a. The mass term cancels immediately when $\mathbf{z}=\mathbf{y}$, because its two Dirac delta distributions coincide. For the derivative term first keep the two points distinct:The terms in which the derivative acts on ${\psi }_c$ cancel against each other; the derivative of the second Dirac delta distribution remains. ThusThe last equality is integration by parts, or the defining derivative rule for a distribution. Consequently the literal conserved current whose divergence equals $C$ isBecause a commutator of two Hermitian operators is anti-Hermitian, this $\mathbf{J}$ includes $-i$. The physical Hermitian Dirac current is instead $j^k=\bar{\psi }{\gamma }^k\psi$, so $C=-i\nabla \cdot \mathbf{j}$. The Heisenberg equation of motion now givesThis establishes the continuity equation directly from the required commutator, without assuming the Dirac equation or invoking the Noether theorem as a substitute for the computation.

In the scalar two-particle model, the fermionic exchange symmetry and bosonic exchange symmetry conditions areFor particles with spin, exchange acts on all their labels; the stated coordinate-only model fixes the spatial symmetry sector. Define $S\Psi (x_1,x_2)=\Psi (x_2,x_1)$. Equal masses make the kinetic operator invariant under exchange, and the assumed potential is invariant too, so $HS=SH$. Therefore $H\Psi =E\Psi$ implies $H(S\Psi )=E(S\Psi )$. The componentsare respectively symmetric and antisymmetric eigenfunctions whenever nonzero. Their sum is the original nonzero quantum state, so at least one is nonzero. An eigenspace may contain both sectors if the energy is degenerate.

The normal modes of two harmonically coupled identical particles provide the separation. The coordinate transformation to $y_{+},y_{-}$ is orthogonal, with absolute Jacobian one. Exchanging the particles leaves $y_{+}$ fixed and changes $y_{-}$ to $-y_{-}$. Hence bosonic quantum states are even and fermionic quantum states odd in the relative coordinate. The chain rule givesAlso $x_1^2+x_2^2=y_{+}^2+y_{-}^2$ and $(x_1-x_2{)}^2=2y_{-}^2$. Combining the two trapping potentials with the interaction therefore gives the separated Hamiltonian operatorThese are independent oscillators of frequencies $\omega$ and $\omega \cos \alpha$. Both frequencies are positive in the stated range. Their product eigenstates have energiesThe parity of an oscillator level is $(-1{)}^n$. Consequently bosons allow even $n_{-}$ and fermions odd $n_{-}$, while $n_{+}$ is unrestricted.

The bosonic ground state uses $n_{+}=n_{-}=0$, and the fermionic ground state uses $n_{+}=0,n_{-}=1$. The oscillator length in the relative mode isUsing the normalized Gaussian for each mode givesIndeed ${\varphi }_{\tau }(y)=\sqrt{2}y{\psi }_{\tau }(y)/\tau$ by differentiating the Gaussian, so it is the normalized odd first excited state. The ground energies are $E_{\mathrm{B}}=\hslash \omega (1+\cos \alpha )/2$ and $E_{\mathrm{F}}=\hslash \omega (1+3\cos \alpha )/2$.

Choose normalized factors for a normalized product quantum state; this can always be done by reciprocal rescaling of its factors. Since the coordinate Jacobian has absolute value one,If the supplied factors are not individually normalized, the general expression has an additional factor $\Vert \psi {\Vert }^2$; choosing normalized factors is the convention in the displayed formula.

For the ground states, put $\rho (y)=|{\psi }_{\tau }(y){|}^2=e^{-y^2/{\tau }^2}/(\sqrt{\pi }\tau )$. Then $|{\varphi }_{\tau }(y){|}^2=2y^2\rho (y)/{\tau }^2$ and ${\rho }^{\prime }=-2y\rho /{\tau }^2$. Thus, for differentiable $f$ with suitable Gaussian integrability and vanishing boundary terms,This is integration by parts, and every polynomial satisfies the required endpoint conditions. The relative-coordinate moments in even and odd oscillator ground states follow by substituting $f(y)=y^n$:For odd $n$, both expectations vanish by parity, so this is an equality of zeros rather than an increase in signed mean separation. The even moments show the physical difference. In particular ${\mathbb{E}}_{\mathrm{B}}[y_{-}^2]={\tau }^2/2$ and ${\mathbb{E}}_{\mathrm{F}}[y_{-}^2]=3{\tau }^2/2$, so the fermionic root-mean-square particle separation is $\sqrt{3}$ times the bosonic value. The fermionic probability density vanishes at coincidence $y_{-}=0$, expressing Pauli exclusion principle in this spatial symmetry sector. This is an exchange constraint on the wave function, not an additional classical force. As $\alpha$ approaches $\pi /2$, the relative oscillator softens and its width diverges; the excluded endpoint would no longer have this normalizable oscillator ground state.

Choose an orthonormal eigenbasis $H_0|n⟩=E_n|n⟩$ and the phase convention $⟨0|{\psi }^{\prime }⟩=0$. Comparing coefficients of the eigenvalue equation at first order yieldsProjecting onto $|0⟩$ and the other eigenvectors givesProjection of the next-order equation onto $|0⟩$ givesThese are coefficients of the printed series, so $E^{\prime \prime }$ denotes the quadratic coefficient rather than the second derivative itself. Without a phase convention an arbitrary purely imaginary multiple of ${\psi }_0$ can be added to the first vector correction. Continuous spectrum, if present, requires the corresponding spectral integral; the oscillator has a discrete complete eigenbasis.

Put $K=m{\omega }^2/2$ and let $g(x)=(2\pi {\sigma }^2{)}^{-1/2}{e}^{-x^2/(2{\sigma }^2)}$ denote the ground-state density. Then $H^{\prime }=K(a^2-x^2)1_{[-a,a]}$ and $p={\int }_{-a}^{a}g(x)dx$. Because $g^{\prime }=-xg/{\sigma }^2$, integration by parts yieldsThus the two first-order shifts areThe first is the supplied ground-state expression; the second follows because the excited-state density is $x^{2}g/{\sigma }^2$, so $E_1^{\prime }=K(a^2{I}_2-I_4)/{\sigma }^2$.

Subtracting, and using $K{\sigma }^2=\hslash \omega /4$, givesThe bracket is zero at $a=0$, tends to one at infinity, and has derivative $2a^{2}g(a)/{\sigma }^2\ge 0$. Therefore the first-order predicted level gap isFor every $\lambda >2$ it becomes negative for all sufficiently large $a$, crossing when $p-2ag=2/\lambda$. It is still positive at $a=0$: the stated conclusion must be read with the large-$a$ regime, rather than for every interval width.

For $\lambda >1$ the inner potential curves downward, placing a central barrier between two symmetric minima at $\pm a$. The even unperturbed ground state has its greatest density at the central barrier, while the odd first excited state has a node there and is less penalized. This explains the relative first-order shifts. A negative predicted gap signals failure of that perturbation approximation, not inversion of the exact lowest two states: the true one-dimensional ground state remains nodeless and even, with the odd state above it. At a high barrier the actual low levels instead form a close double-well pair. The apparent level crossing from a finite-interval oscillator perturbation is useful precisely as a diagnostic of the approximation's range.

Let $u=rF$. Multiplication of the radial equation by $-2mr/{\hslash }^2$ gives $u^{\prime \prime }+[A/r-\lambda (\lambda +1)/r^2-{\alpha }^2]u=0$. With $u=e^{-\alpha r}f$ the constant terms cancel, leavingThe indicial equation is $c(c-1)=\lambda (\lambda +1)$, with roots $c=\lambda +1$ and $c=-\lambda$. The usual physical Coulomb radial origin condition selects $c=\lambda +1$: it imposes the regular reduced radial boundary condition $u(0)=0$. Normalizability alone at the origin is insufficient to make this selection when an irregular branch is locally square-integrable, notably at $\lambda =0$.

Write $f=r^{\lambda +1}{\sum }_{j\ge 0}{a}_j{r}^j$. Coefficient comparison givesA nonterminating regular series has the exponentially growing asymptotic branch, suggested by $a_{j+1}/a_j\sim 2\alpha /j$: it supplies an $e^{2\alpha r}$ component, which the prefactor $e^{-\alpha r}$ cannot normalize. Equivalently the regular confluent-hypergeometric solution grows unless its first parameter is a nonpositive integer. Thus a bound state requires the numerator to vanish at some $j=k\ge 0$, givingThis is series-termination quantization of the Coulomb radial equation. It requires $A>0$ for the assumed $\alpha >0$; at $A=0$ there are no such negative-energy Coulomb bound states.

For $k=0$, $F=N{r}^{\lambda }{e}^{-\alpha r}$. The angular integral cancels when computing $⟨1/r⟩$, while the radial volume element contributes $r^{2}dr$. ConsequentlyThe ratio follows by one integration by parts and remains valid for real $\lambda \ge 0$; no integer-only factorial identity is needed. Since now $\alpha =A/[2(\lambda +1)]$,For a classical particle with that negative total energy, nonnegative kinetic energy requires $E-V={\hslash }^2(A/r-{\alpha }^2)/(2m)\ge 0$, hence $\boxed{r\le A/{\alpha }^2}$. This is a necessary radial bound; extra angular momentum can impose a stronger restriction.

At $k=\lambda =0$, $\alpha =A/2$ and the classical boundary is $r=2/\alpha$. The normalized radial probability density is $4{\alpha }^3{r}^2{e}^{-2\alpha r}$. With $z=2\alpha r$, the exterior quantum probability isThe final integration gives $(4^2+2\cdot 4+2)e^{-4}/2$. This is the nonzero quantum tail beyond the classical turning region.

Choose an orthonormal eigenbasis $H_0|n⟩=E_n|n⟩$ and the phase convention $⟨0|{\psi }^{\prime }⟩=0$. Comparing coefficients of the eigenvalue equation at first order yieldsProjecting onto $|0⟩$ and the other eigenvectors givesProjection of the next-order equation onto $|0⟩$ givesThese are coefficients of the printed series, so $E^{\prime \prime }$ denotes the quadratic coefficient rather than the second derivative itself. Without a phase convention an arbitrary purely imaginary multiple of ${\psi }_0$ can be added to the first vector correction. Continuous spectrum, if present, requires the corresponding spectral integral; the oscillator has a discrete complete eigenbasis.

Put $K=m{\omega }^2/2$ and let $g(x)=(2\pi {\sigma }^2{)}^{-1/2}{e}^{-x^2/(2{\sigma }^2)}$ denote the ground-state density. Then $H^{\prime }=K(a^2-x^2)1_{[-a,a]}$ and $p={\int }_{-a}^{a}g(x)dx$. Because $g^{\prime }=-xg/{\sigma }^2$, integration by parts yieldsThus the two first-order shifts areThe first is the supplied ground-state expression; the second follows because the excited-state density is $x^{2}g/{\sigma }^2$, so $E_1^{\prime }=K(a^2{I}_2-I_4)/{\sigma }^2$.

Subtracting, and using $K{\sigma }^2=\hslash \omega /4$, givesThe bracket is zero at $a=0$, tends to one at infinity, and has derivative $2a^{2}g(a)/{\sigma }^2\ge 0$. Therefore the first-order predicted level gap isFor every $\lambda >2$ it becomes negative for all sufficiently large $a$, crossing when $p-2ag=2/\lambda$. It is still positive at $a=0$: the stated conclusion must be read with the large-$a$ regime, rather than for every interval width.

For $\lambda >1$ the inner potential curves downward, placing a central barrier between two symmetric minima at $\pm a$. The even unperturbed ground state has its greatest density at the central barrier, while the odd first excited state has a node there and is less penalized. This explains the relative first-order shifts. A negative predicted gap signals failure of that perturbation approximation, not inversion of the exact lowest two states: the true one-dimensional ground state remains nodeless and even, with the odd state above it. At a high barrier the actual low levels instead form a close double-well pair. The apparent level crossing from a finite-interval oscillator perturbation is useful precisely as a diagnostic of the approximation's range.

In the scalar two-particle model, the fermionic exchange symmetry and bosonic exchange symmetry conditions areFor particles with spin, exchange acts on all their labels; the stated coordinate-only model fixes the spatial symmetry sector. Define $S\Psi (x_1,x_2)=\Psi (x_2,x_1)$. Equal masses make the kinetic operator invariant under exchange, and the assumed potential is invariant too, so $HS=SH$. Therefore $H\Psi =E\Psi$ implies $H(S\Psi )=E(S\Psi )$. The componentsare respectively symmetric and antisymmetric eigenfunctions whenever nonzero. Their sum is the original nonzero quantum state, so at least one is nonzero. An eigenspace may contain both sectors if the energy is degenerate.

The normal modes of two harmonically coupled identical particles provide the separation. The coordinate transformation to $y_{+},y_{-}$ is orthogonal, with absolute Jacobian one. Exchanging the particles leaves $y_{+}$ fixed and changes $y_{-}$ to $-y_{-}$. Hence bosonic quantum states are even and fermionic quantum states odd in the relative coordinate. The chain rule givesAlso $x_1^2+x_2^2=y_{+}^2+y_{-}^2$ and $(x_1-x_2{)}^2=2y_{-}^2$. Combining the two trapping potentials with the interaction therefore gives the separated Hamiltonian operatorThese are independent oscillators of frequencies $\omega$ and $\omega \cos \alpha$. Both frequencies are positive in the stated range. Their product eigenstates have energiesThe parity of an oscillator level is $(-1{)}^n$. Consequently bosons allow even $n_{-}$ and fermions odd $n_{-}$, while $n_{+}$ is unrestricted.

The bosonic ground state uses $n_{+}=n_{-}=0$, and the fermionic ground state uses $n_{+}=0,n_{-}=1$. The oscillator length in the relative mode isUsing the normalized Gaussian for each mode givesIndeed ${\varphi }_{\tau }(y)=\sqrt{2}y{\psi }_{\tau }(y)/\tau$ by differentiating the Gaussian, so it is the normalized odd first excited state. The ground energies are $E_{\mathrm{B}}=\hslash \omega (1+\cos \alpha )/2$ and $E_{\mathrm{F}}=\hslash \omega (1+3\cos \alpha )/2$.

Choose normalized factors for a normalized product quantum state; this can always be done by reciprocal rescaling of its factors. Since the coordinate Jacobian has absolute value one,If the supplied factors are not individually normalized, the general expression has an additional factor $\Vert \psi {\Vert }^2$; choosing normalized factors is the convention in the displayed formula.

For the ground states, put $\rho (y)=|{\psi }_{\tau }(y){|}^2=e^{-y^2/{\tau }^2}/(\sqrt{\pi }\tau )$. Then $|{\varphi }_{\tau }(y){|}^2=2y^2\rho (y)/{\tau }^2$ and ${\rho }^{\prime }=-2y\rho /{\tau }^2$. Thus, for differentiable $f$ with suitable Gaussian integrability and vanishing boundary terms,This is integration by parts, and every polynomial satisfies the required endpoint conditions. The relative-coordinate moments in even and odd oscillator ground states follow by substituting $f(y)=y^n$:For odd $n$, both expectations vanish by parity, so this is an equality of zeros rather than an increase in signed mean separation. The even moments show the physical difference. In particular ${\mathbb{E}}_{\mathrm{B}}[y_{-}^2]={\tau }^2/2$ and ${\mathbb{E}}_{\mathrm{F}}[y_{-}^2]=3{\tau }^2/2$, so the fermionic root-mean-square particle separation is $\sqrt{3}$ times the bosonic value. The fermionic probability density vanishes at coincidence $y_{-}=0$, expressing Pauli exclusion principle in this spatial symmetry sector. This is an exchange constraint on the wave function, not an additional classical force. As $\alpha$ approaches $\pi /2$, the relative oscillator softens and its width diverges; the excluded endpoint would no longer have this normalizable oscillator ground state.

Inside the infinite square well, the time-independent Schrödinger equation and Dirichlet boundary conditions areand the wavefunction vanishes outside the interval. For a nonzero state, integration by parts givesThe strict inequality follows because a function with zero derivative and zero endpoint values is zero. Put $k=\sqrt{2mE}/\hslash$. The general solution is $A\sin kx+B\cos kx$; the first boundary condition gives $B=0$, and the second gives $ka=n\pi$, $n=1,2,\dots$. Since ${\int }_0^a{\sin }^2(n\pi x/a)dx=a/2$, the normalized energy eigenstates and eigenvalues areThe normalization fixes the magnitude of the coefficient; its arbitrary global phase has been chosen real and positive.