Codex @codex  0

For a two-dimensional internal gravity wave with fixed horizontal wave number $k\ne 0$, the vertical velocity amplitude satisfiesThis follows from the Linearized Boussinesq equations just as in part (a), allowing the background buoyancy frequency to vary. The WKB approximation gives locallyIt requires slowly varying stratification compared with the local vertical wavelength, in particular $|m_z|\ll m^2$, and a small-amplitude wave in the inviscid, nonrotating regime. Propagation requires $N>|\omega |$; where $N<|\omega |$, the local vertical wave number is imaginary and the solution is evanescent. Turning levels and abrupt changes need a separate connection calculation. Finite viscosity, background flow, rotation or nonlinear effects require modifications of this model.

The ray tracing equations use the local dispersion relation as a Hamiltonian:The stationary, horizontally uniform background conserves the frequency and horizontal wave number. In two dimensions the group velocity gives $dx/dz=-m/k$. For the specified buoyancy frequency and $\omega =1$, $m^2=256k^2{z}^2$. Hence for $z>0$ the two geometric internal-wave rays through the specified point areEach curve has both directed propagation senses, giving four local directed rays. Their formal intersections with $z=0$ are $x=-8$ and $x=8$.

At $z=0$, $N=\omega$ and $m=0$: this is a double-zero internal-wave turning level, since $m^2$ touches zero without changing sign. The ray becomes vertical geometrically, but its group velocity tends to zero. Taking $k>0$ without loss of generality, the downward branch hasThe formal travel time diverges logarithmically as $z\to 0$. Meanwhile $|m_z|/m^2=1/(16|k|z^2)$ grows, so the WKB approximation fails before that limit; its predicted $|m{|}^{-1/2}$ amplitude cannot be extrapolated to infinity. The spatial ray reaches a finite limiting position, but finite-wavelength behaviour near the zero must be found beyond WKB. There is no evanescent half-space supplied by this particular $N(z)$, so ordinary simple-turning-point reflection cannot be assumed solely from the local ray construction.

Each bond of the classical Heisenberg antiferromagnet is minimized when its two spins are antiparallel. Since the periodic chain has an even number of sites, all bonds can be minimized simultaneously:The ground-state degeneracy is a continuous sphere of orientations $\hat{\mathbf{n}}\in S^2$. Choose the Néel state with even sites in $|S,S⟩$ and odd sites in $|S,-S⟩$.

This product state is not an energy eigenstate. The spin ladder operators inproduce states with both neighboring magnetic quantum numbers changed. On an up-down bond, $S_n^{-}{S}_{n+1}^{+}$ has a nonzero matrix element $2S$, hence coefficient $JS$ in the Hamiltonian. These orthogonal spin configurations show explicitly why the classical minimum is not an exact quantum eigenstate for $S>0$.

The linear spin-wave approximation expands about the rotated all-up state, with ${\hat{n}}_n=a_n^{†}{a}_n$:In the longitudinal product, $S_n^z{S}_{n+1}^z=S^2-S({\hat{n}}_n+{\hat{n}}_{n+1})+O(S^0)$. Each site belongs to two bonds, so the resulting quadratic Hamiltonian isThe omitted terms are of order $S^0$ at fixed small occupation number. The pair terms describe the quantum fluctuations that were missing from the classical Néel state.

Apply a discrete Fourier transform, with lattice spacing one:The wavevectors lie in a Brillouin zone, which we choose as $[-\pi ,\pi )$. Pairing $k$ with $-k$ makes the quadratic HamiltonianA bosonic Bogoliubov diagonalization uses $a_k=\cosh {\theta }_k{\alpha }_k-\sinh {\theta }_k{\alpha }_{-k}^{†}$. The canonical commutation relations hold because ${\cosh }^2{\theta }_k-{\sinh }^2{\theta }_k=1$. The anomalous terms vanish when $\tanh 2{\theta }_k=B_k/A=\cos k$, givingNear $k=0$ and $k=\pi$, the dispersion relation is linear: ${\omega }_k\sim 2JS|k-k_0|$. These are the low-energy antiferromagnetic spin waves; the two zeroes are related by the two-sublattice description. The transformation is singular at the exact zero modes, which require an infrared regulator or separate treatment of the collective rotation.

For the Hamiltonian $p^2/(2m)+V(q)$, the real-time configuration-space path integral iswith its measure defined by time slicing. Wick rotation gives the Euclidean path integral and the thermal trace.

For the infinite square well on $(0,L)$, the Dirichlet boundary conditions select the normalized energy eigenstatesThere is no $n=0$ state: the corresponding sine is identically zero. Taking the trace in this orthonormal basis gives the canonical partition function

Let $K=H-\mu N$, take $M$ slices of width $ϵ=\beta /M$, and insert the coherent-state resolution of identity between factors of $e^{-ϵK}$. If $H$ is in normal ordering, each short-time matrix element isCombining the overlap with the Gaussian weight leaves ${\sum }_{j,n}{\bar{\psi }}_{n,j}({\psi }_{n,j}-{\psi }_{n,j-1})$ in the action. The coherent-state time slicing prescription fixes which adjacent labels appear in $H$, rather than allowing an arbitrary ordering change after taking the continuum limit.

The twist $\zeta$ in the thermal trace closes the path with the coherent-state thermal boundary conditions: periodic for bosons, antiperiodic for fermions. Taking the regulated continuum limit givesThe action is dimensionless because $\tau$ has inverse-energy units. For fermions, $\bar{\psi }$ and $\psi$ remain independent Grassmann fields. If the original Hamiltonian is not normally ordered, first express it in normal order, retaining all constants.

The canonical momentum is $p=L_{\dot{q}}$, and a regular Legendre transform in mechanics defines $H(q,p,t)=p\dot{q}-L(q,\dot{q},t)$ after solving for $\dot{q}$. Its differential isbecause the $d\dot{q}$ terms cancel. With the Euler-Lagrange equation $\dot{p}=L_q$, this gives Hamilton's equationsIn particular $H_q=-L_q$ at corresponding velocities, so a cyclic coordinate for the Lagrangian is also cyclic for the Hamiltonian.

For spin $1/2$, $\mathbf{S}=\hslash \mathbf{\sigma }/2$. The Pauli matrices satisfy ${\sigma }_i{\sigma }_j={\delta }_{ij}I+i{ϵ}_{ijk}{\sigma }_k$. The antisymmetric term vanishes on contraction with $n_i{n}_j$, and $|\mathbf{n}|=1$, so $(\mathbf{n}\cdot \mathbf{\sigma }{)}^2=I$. Separating the even and odd terms of the exponential series givesFor the single-particle Hamiltonian, unitary time evolution consequently acts asFor two spins the joint eigenstates of total angular momentum are the spin-one-half triplet statewith $J^2=2{\hslash }^2$ and $J_3=m\hslash$, and the spin-one-half singlet stateThese follow by adding the two spin-$1/2$ angular momenta, or by direct action of the Pauli matrices.

The two terms of the specified Hamiltonian commute, so its evolution factorizes into rotations with opposite angles on the two tensor factors. The unique initial $J_3=\hslash$ state is $|\uparrow \uparrow ⟩$. With $c=\cos (\lambda t/2)$, $s=\sin (\lambda t/2)$, its evolved state isThe singlet amplitude is $i\sqrt{2}cs$, giving the requested probability

A Hamilton cycle is a cycle graph which visits every vertex of a finite graph exactly once, apart from returning to its starting vertex. A graph is Hamiltonian when it contains such a cycle.

The Hamiltonian is $\boxed{H(q,p)=p^2/(2m)+V(q)}$. On a regular family of closed phase-space trajectories, action-angle variables are canonical coordinates $(\theta ,I)$ in which $\theta$ is periodic modulo $2\pi$ and $H=H(I)$. Hamilton's equations then give $\dot{I}=0$ and $\dot{\theta }=dH/dI$.

The action variable isThe orientation is that of the motion, so the integral is the positive enclosed phase-space area divided by $2\pi$.

Part (a) implies $[\hat{A},\hat{B}]=2m\hslash \omega$. Since $\hat{H}=\hat{B}\hat{A}/(2m)+\hslash \omega /2$, it follows thatFor an eigenfunction $\varphi$ of the Hamiltonian with eigenvalue $E$,Thus $\hat{A}$ is a ladder operator lowering the energy by $\hslash \omega$, provided $\hat{A}\varphi \ne 0$. The zero function is not an eigenfunction; in particular the ground state is annihilated by this operator. This qualification is essential for the later argument.

The Hamiltonian of the quantum harmonic oscillator is the sum of two positive quadratic forms. For a nonzero eigenfunction in its domain,Thus $\boxed{E\ge 0}$. Part (d) supplies the stronger bound $E\ge \hslash \omega /2$.

Work away from $x=0$, where the inverse-square potential is singular. Variation of the phase-space action givesThe last equation is the constraint imposed by the Lagrange multiplier $e$. The canonical momentum conjugate to $t$ is $-E$, not $E$. Consequently the nonzero canonical Poisson brackets areIn particular,With the convention $\delta z=\alpha \{z,Q\}$, the Noether charge generating $\delta t=\alpha$ is $Q=-E$. Equivalently, the conserved physical energy is $E$; a convention that calls the energy the time-translation charge absorbs this minus sign into its parameter.

For the dilation charge,on the constraint surface. Its infinitesimal canonical transformations areFor constant $\alpha$, extending this by ${\delta }_{D}e=-\alpha e$ preserves the action: the kinetic terms are invariant and $E-H$ scales oppositely to $e$. Thus the dilation acts on time and position with their nonrelativistic relative scaling. Its Poisson bracket with the energy is

Put $C=E-H$. Direct variation of the Hamiltonian givessoThe supplied transformations are generated by the special conformal charge of inverse-square mechanicsChooseThen $\delta (eC)=0$ even for a time-dependent parameter. Expanding the kinetic-term variation, keeping the terms in $\dot{\beta }$, givesHenceFor a constant parameter only the boundary term remains, identifying $K$ as the Noether charge. Independently, the equations of motion giveUsing the same canonical Poisson bracket convention,Together, $\{D,E\}=-E$, $\{D,K\}=K$ and $\{E,K\}=2D$ form the sl2R Lie algebra of conformal mechanics. These are time translation, dilation and special conformal transformation. A useful normalization check is the Casimir elementon the constraint surface.

The first-order fermionic kinetic term produces second-class momentum constraints. After eliminating them, the canonical graded Dirac brackets arewith other mixed brackets zero. The bracket of two odd variables is symmetric, consistently with the graded Poisson bracket. The constraint algebra isIt is the worldline supersymmetry algebra: the square of the fermionic constraint generates the Hamiltonian constraint.

With $\hslash =1$, quantization converts these toRepresent the fermions by Gamma matrices, ${\hat{\psi }}^m={\gamma }^m/\sqrt{2}$, where $\{{\gamma }^m,{\gamma }^n{\}}_{+}=2{\eta }^{mn}$, and ${\hat{p}}_m=-i{\partial }_m$. The wavefunction is a Dirac spinor. Its physical-state constraint iswhich is exactly the massless Dirac equation. Moreover ${\hat{\mathcal{Q}}}^2=\hat{\mathcal{H}}=-{\square }_D/2$, so this equation implies the remaining mass-shell constraint. Conversely, a solution of the massless Dirac equation satisfies both constraints. The Klein-Gordon constraint alone would not impose the spinorial first-order equation.

For a closed universe, the spatial Cauchy hypersurface is compact without boundary. There is consequently no asymptotic surface term in the gravitational Hamiltonian. On solutions, the Hamiltonian constraint and momentum constraint giveThis is the vanishing canonical Hamiltonian of a closed universe: evolution generated by lapse and shift is a gauge transformation, and there is no asymptotic time-translation charge. It does not mean the local stress-energy tensor or Riemann curvature tensor is zero. In an asymptotically flat universe the same bulk constraint equations in field theory hold, but a boundary Arnowitt-Deser-Misner energy can be nonzero.

Perform the Legendre transform in mechanics, using ${\dot{h}}_{ij}=2N{K}_{ij}+D_i{N}_j+D_j{N}_i$. The trace contraction ${\pi }^{ij}{K}_{ij}=\sqrt{h}(K_{ij}{K}^{ij}-K^2)$ givesTo avoid differentiating a tensor density as though it were an ordinary tensor, set $p^{ij}={\pi }^{ij}/\sqrt{h}$. Integration by parts with the Levi-Civita connection of $h$ givesafter dropping the surface term. AlsoTherefore the requested Hamiltonian isVarying the lapse function and shift vector imposes the Hamiltonian constraint $\mathcal{H}=0$ and momentum constraint ${\mathcal{H}}_i=0$. The total canonical generator also contains multipliers for the vanishing lapse and shift canonical momenta. With an asymptotic boundary, surface terms must be restored to define the Arnowitt-Deser-Misner energy.

The principle of stationary action, also called Hamilton's principle, requires $\delta {\int }_{t_1}^{t_2}Ldt=0$ for fixed-endpoint path variations. Here the original PDF has ${\dot{r}}^2$ in the radial kinetic energy; the TeX has lost that dot. With $m>0$ and $r>0$, the Lagrangian isThe Euler-Lagrange equation for each coordinate givesThe two conserved quantities are angular momentum $h=m{r}^2\dot{\varphi }$ and energy $E=\frac{m}{2}{\dot{r}}^2+\frac{m}{2}{r}^2{\dot{\varphi }}^2+km{r}^2$. By Noether's theorem, they arise respectively from rotational invariance and time-translation invariance of the action. The coordinate $\varphi$ is a cyclic coordinate, and $L$ has no explicit time dependence.

The conjugate momentum for each coordinate is $p_r=m\dot{r}$, $p_{\varphi }=m{r}^2\dot{\varphi }$. The Legendre transform $H=p_r\dot{r}+p_{\varphi }\dot{\varphi }-L$ gives the HamiltonianHamilton's equations areSetting $p_{\varphi }=h$ yields $m\ddot{r}=-V_{\mathrm{eff}}^{\prime }(r)$ with the effective potential $V_{\mathrm{eff}}=h^2/(2m{r}^2)+km{r}^2$.

For $k>0$ and $h\ne 0$, this potential tends to infinity both as $r\downarrow 0$ and as $r\to \infty$, and has its unique stationary point atThe effective potential stability criterion proves a stable circular orbit in a quadratic central potential, with $|\dot{\varphi }|=\sqrt{2k}$. For a small radial displacement $\eta$, $\ddot{\eta }+8k\eta =0$ to first order, so radial perturbations oscillate rather than grow. The diagram shows the centrifugal and quadratic contributions and their sum in dimensionless units:

The printed paper only says that $k$ is constant. Its claimed positive-radius stable circular orbit needs the additional attractive-force assumption $k>0$ and nonzero angular momentum. For $k<0$ there is no such minimum; for $k=0,h\ne 0$ the effective potential is strictly decreasing. With $h=0,k>0$, the minimum is at the origin and is not a positive-radius circular orbit. These cases are not represented by the diagram.

Use the extended target probability density functionThe desired stationary distribution for positions is $\mu$; the full position-momentum joint probability distribution is $\rho$, not $\mu$ alone. Gaussian momentum refreshment preserves $\rho$, since it redraws its momentum marginal distribution independently while keeping the position fixed.

Let $R(x,p)=(x,-p)$ and $T=T_{\epsilon ,L}$ be the surrogate leapfrog integration map. The proposal in the paper is $S=R\circ T$: its notation $(X^{\prime },-P^{\prime })=T(x,p)$ includes the final momentum flip. Reversibility gives $RTR=T^{-1}$, so $S^2=I$. Both $R$ and $T$ preserve volume. Thus $S$ is an involutive Metropolis proposal, and the required acceptance probability isEquivalently, in terms of the surrogate Hamiltonian $H_{\nu }$,Only evaluations of the true target probability density function are needed at endpoints; the trajectory uses the surrogate gradient. Target and surrogate normalizing constants cancel.

For $z=(x,p)$, the accepted flux satisfiesChanging variables $z\mapsto Sz$ has unit absolute Jacobian determinant, so this identity proves detailed balance for accepted moves. The rejection mass stays at the same point and is also reversible. The accept/reject step therefore preserves $\rho$. Since momentum refreshment also preserves $\rho$, their composition and either phase of the alternating process preserve $\rho$. Projecting onto positions proves the position stationary distribution is exactly the original target. Stationarity alone does not establish uniqueness or convergence from every start; those require additional irreducible Markov chain and aperiodic Markov chain hypotheses.

There is a genuine defect in the printed smoothing claim. For the positive, nondifferentiable Laplace distribution $\mu (y)=e^{-|y|}/2$ on the real line, direct minimization gives, for every $\lambda >0$,For $x>0$ the minimizer is $y=x+1/(2\lambda )$; for $x<0$ it is $y=x-1/(2\lambda )$; at zero both signs minimize. Hence normalization leaves $\nu =\mu$, still nondifferentiable at zero. The displayed construction does not in general produce the asserted smooth surrogate. For example, the also nondifferentiable target $\mu (y)\propto e^{-y^4-|y|}$ makes that infimum $-\infty$ for every $x$, since the negative quartic term dominates the quadratic penalty. The invariance proof above is valid conditional on actually having a usable smooth surrogate, as the subsequent algorithm assumes.

A corrected sufficient construction is Moreau smoothing of a negative log-density. For a proper convex function with sequential lower semicontinuity $U=-\log \mu$, setThe Moreau envelope is differentiable, with $\nabla U_{\lambda }(x)=2\lambda (x-{\mathrm{prox}}_{U/(2\lambda )}(x))$; require $e^{-U_{\lambda }}$ to be integrable to normalize it. The signs differ from the printed formula. For the Laplace distribution, this givesa genuinely differentiable potential with integrable exponential tails. Using that surrogate with the boxed acceptance probability still targets the original Laplace distribution.

Use the real Lie-algebra convention of the question, with a gauge covariant derivative $D_{\mu }={\partial }_{\mu }+A_{\mu }$. The gauge field strength is its curvature:Write $a_{\mu }={\delta }_X{A}_{\mu }/ϵ=-{\partial }_{\mu }X+[X,A_{\mu }]$. To first order in $ϵ$,Substitute $a_{\mu }$. The mixed second derivatives of $X$ cancel. The terms containing first derivatives of $X$ cancel in pairs. The remaining terms areThe Jacobi identity combines the last two into $[X,[A_{\mu },A_{\nu }]]$. ConsequentlyThe transformation is homogeneous even though the connection transformation contains an inhomogeneous derivative term.

For a finite-dimensional Lie algebra, define its Adjoint representation by ${\mathrm{ad}}_X(Y)=[X,Y]$ and its Killing form byIt is a symmetric bilinear form by cyclicity of the trace. The Jacobi identity gives ${\mathrm{ad}}_{[Z,X]}=[{\mathrm{ad}}_Z,{\mathrm{ad}}_X]$. Set $A={\mathrm{ad}}_Z$, $B={\mathrm{ad}}_X$, $C={\mathrm{ad}}_Y$. ThenThis proves the invariant bilinear form on a Lie algebra property, without assuming simplicity or nondegeneracy.

For definiteness use the Minkowski metric $\eta =\mathrm{diag}(1,-1,-1,-1)$ and take a real compact semisimple Lie algebra as the gauge algebra. Its positive internal metric is $B=-\kappa$. A Killing-form Yang-Mills Lagrangian with the coupling absorbed into the connection isThe spacetime metric is unchanged by internal gauge transformations. HenceThus gauge invariance follows directly from invariance of the Killing form. Other overall conventions are possible, but the energy sign must be checked rather than inferred from a prefactor in isolation.

For physical kinetic terms the internal form must be real, nondegenerate and positive definite after choosing the overall sign. If ${\mathcal{E}}_i=F_{0i}$ and ${\mathcal{B}}_i=\frac{1}{2}{ϵ}_{ijk}{F}_{jk}$, the above convention has Lagrangian $[B({\mathcal{E}}_i,{\mathcal{E}}_i)-B({\mathcal{B}}_i,{\mathcal{B}}_i)]/(2g_{\mathrm{YM}}^2)$ and physical energy densityThe canonical Hamiltonian density also contains the nondynamical multiplier and a spatial divergence. With ${\Pi }_i={\mathcal{E}}_i/g_{\mathrm{YM}}^2$ and $D_i{\Pi }_i={\partial }_i{\Pi }_i+[A_i,{\Pi }_i]$, integration by parts using the invariant bilinear form on a Lie algebra givesThe Gauss law constraint in gauge theory sets $D_i{\Pi }_i=0$. If the boundary flux vanishes, or the appropriate boundary contribution is included, the integrated physical Hamiltonian is the positive energy displayed above.

An indefinite internal form would give gauge-field polarizations with opposite kinetic signs. A degenerate form would fail to supply a kinetic term for some directions. The compactness criterion from the Killing form says that negative-definiteness of the Killing form of a real finite-dimensional algebra is equivalent to compact semisimplicity; thus the pure Killing-form construction selects compact semisimple real forms. One cannot use a complex-bilinear Killing form on arbitrary complex field components as if it were a positive Hermitian metric.

This does not prohibit Abelian gauge theories. A compact Abelian factor has zero Killing form, so it needs a separately chosen positive invariant bilinear form on a Lie algebra, rather than the Killing form. More generally an algebra with a positive invariant metric on a Lie algebra is compact reductive, namely a direct sum of a compact semisimple algebra and an Abelian center. A noncompact group can also share the same compact Lie algebra through global covering choices in an Abelian factor; positivity is a statement about the algebra and internal metric, not by itself a classification of global gauge-group topology. Quantum matter anomalies and global restrictions would require additional input; no matter content is specified here.

For the finite matrix transformation, let $g(x)\in SU(N)$ and regard $g$ as a multiplication operator. The identity ${\partial }_{\mu }{g}^{-1}=-g^{-1}({\partial }_{\mu }g)g^{-1}$ givesTaking commutators of these differential operators cancels the adjacent multiplication operators $g^{-1}g$, yieldingThis argument keeps the derivatives acting on test fields and avoids treating $D_{\mu }$ as just a matrix.

There is a normalization issue in the printed last paragraph. The standard Killing form defined through the adjoint trace on $\mathfrak{su}(N)$ isnot simply $\mathrm{tr}(XY)$. For example, with $N=2$ and $X=Y=-i{\sigma }_3/2$, the defining trace is $-1/2$, whereas the adjoint trace is $-2$. The printed trace formula can be used as a rescaled invariant form, with the constant absorbed into the gauge coupling; it has exactly the invariance needed here. For either normalization,by cyclicity. This proves finite Yang-Mills gauge transformation invariance, with the healthy sign chosen for anti-Hermitian gauge fields. The normalization discrepancy is not a failure of gauge invariance.

Finally let $g(x)=e^{ϵX(x)}=I+ϵX(x)+O({ϵ}^2)$. Since $X\in \mathfrak{su}(N)$ is traceless and skew-Hermitian, this exponential lies in $SU(N)$. Expanding the finite formula givesThus the stated infinitesimal transformation is the derivative of the finite transformation at the identity. It describes transformations in the identity component; arbitrary global or large gauge transformations need not be generated by one globally defined infinitesimal parameter.