For the real coefficients of the stated self-adjoint Sturm-Liouville problem, define the weighted inner product
Since , a nonzero eigenfunction has positive squared norm. Let . Integration by parts and the Dirichlet boundary conditions give
Taking to be an eigenfunction shows that its eigenvalue is real, since is real. For two eigenfunctions with eigenvalues , the same identity gives
Hence distinct eigenvalues imply orthogonality in the weighted inner product. This is the orthogonality of Sturm-Liouville eigenfunctions. The calculation assumes the regularity needed for the stated self-adjoint boundary problem, so all boundary terms and integrals are defined.
Set and , so . By the chain rule, the equation becomes . The substitution reduces it further to
Multiplication by and integration by parts gives for a nonzero solution, ruling out zero and negative eigenvalues. For , the first boundary condition gives , and the second gives , . Thus all eigenpairs are
with arbitrary nonzero scalar multiples of each eigenfunction. These are the complete list since all solutions of the constant-coefficient differential equation have been considered; the transformed sine system is also complete in the corresponding weighted Hilbert space.
Use the Fourier transform . For a rapidly decaying smooth function, integration by parts gives , and differentiation under the integral gives . Both are justified here by Gaussian decay. Since ,
The Gaussian integral supplied in the hint gives . Solving this first-order equation yields
The calculation states the transform normalization and the two differentiation properties explicitly.
Write . Assume and have right limits at zero and sufficient regularity and exponential-order bounds for the following integration by parts; take large enough to kill the boundary terms at infinity. First,
Apply the same calculation to :
This is the Laplace transform of a derivative formula with both initial-data terms retained.
Take and a positive prescribed distance ; if , the zero-speed journey has zero cost. Minimize subject to and the isoperimetric constraint . Introduce a constant Lagrange multiplier and the integrand . For a variation with , integration by parts gives
The interior Euler-Lagrange equation is
When terminal speed is free, the same calculation also imposes the natural boundary conditions for a free endpoint:
One cannot drop this boundary term when interpreting an extremal as a genuine minimum. If , the interior condition makes constant; with a continuous zero initial speed and it cannot be attained. In that case Cauchy-Schwarz inequality gives the infimum , approached by speeds nearly constant except in a short initial layer. For , the derivative cost controls such layers and the natural terminal condition matters in parts (b) and (c).
Here because and the two cost coefficients are not both zero. The interior Euler-Lagrange equation makes quadratic. Imposing each displayed pair of endpoint data and the distance constraint gives the following formal profiles and their costs:
Thus choice (3) is the best strategy, and its true minimum cost is
The qualification from part (b) applies to row (1): its formal quadratic is not a minimizer with free terminal speed, since . Row (2) is a genuine minimizer for its fixed zero terminal speed, and row (3) is the genuine free-terminal-speed minimizer, satisfying its natural condition .
A direct lower bound makes the comparison rigorous. For every admissible speed with ,
by integration by parts. The Cauchy-Schwarz inequality gives . Equality requires to be proportional to , giving exactly row (3). Row (1)'s zero initial acceleration excludes equality; however short initial layers, as in part (b), approach that same infimum. Hence choice (1) has infimum but does not attain it, while choice (3) does. This distinction resolves the boundary-condition issue in the printed minimization request.
The function is even, so its Fourier series has . In the convention ,
Twice applying integration by parts gives, for ,
Therefore the periodic Fourier representation is
with periodic continuation. The periodic extension is continuous, including the identified endpoints, and piecewise smooth, so the Fourier series converges to it. The coefficients also give absolute convergence and uniform convergence by comparison with .
At , the cosine cancels the alternating coefficient sign; at , it is identically one. These evaluations yield
The first identity solves the Basel problem.
Write . Smooth compactly supported functions lie in and are dense in L2 space, so is densely defined.
For , use the form
on the one-dimensional harmonic oscillator form domain . It is bounded and coercive in the norm. The Lax-Milgram theorem gives a unique satisfying for every . Thus in distributions.
We must verify the full specified operator domain, not just the form domain. The equation first gives . For compactly supported functions, integration by parts gives the graph estimate for the harmonic oscillator, obtained from
Apply this to , with on , supported in , and derivatives bounded by and . Its right-hand source is , uniformly bounded in L2 space. The estimate and the Fatou lemma give ; the original equation then gives , and hence . Thus , proving that is onto. The energy identity proves injectivity and the resolvent of the shifted harmonic oscillator estimate
To prove closedness, suppose and in L2 space. For one fixed ,
Therefore belongs to the range of the resolvent, namely , and .
All hypotheses of the Hille-Yosida theorem now hold; the resolvent bound is stronger than . Let be its contraction semigroup. The exponentially shifted semigroup
has generator . It is continuous in L2 space, takes the initial value , and satisfies . For , its time derivative exists in L2 space and equals , so the equation holds as a strong L2 space solution. For general initial data it is the semigroup solution, without an unproved pointwise differentiability claim.
A bounded-domain hypothesis is needed for the requested global construction and uniqueness. The printed question does not include it. For example, the upper half-plane satisfies the stated exterior cone condition, but both and are harmonic functions with zero Dirichlet boundary condition. Thus part (iv), in its printed unrestricted solution class, is false. In this question's solutions we add that is bounded; the exterior cone condition and all other data remain as printed.
The half-plane also rules out the printed global separated barrier. In its angular interval , varying on a fixed ray forces if for every . We would then have and . Put and . Twice integration by parts gives
contradicting the positive integrand. Thus a source qualification is necessary for part (i) as well as part (iv).
Fix and its exterior cone of half-angle . Choose and set
On the complement of the cone, unwrap the angle in polar coordinates as , measured from the cone axis. Then is the angle from the opposite axis and on . Define the power barrier for an exterior cone by
This has the requested separated form . In the printed signed-angle convention, the same function is outside the cone. It is smooth across the negative axis: near that axis the angular expression is the even function .
Since , put . The cosine is at least , so away from , and is continuous at . The Laplacian in polar coordinates gives
Let . Since , the barrier for the Dirichlet problem satisfies
Both the uniform lower bound and the later global comparison use boundedness. In particular, is positive at every other boundary point, and is bounded away from zero on boundary sets staying a positive distance from .
Choose a smooth cutoff function supported in , equal to one on , with and . For , put . The zero Dirichlet boundary condition implies that the tangential boundary derivative vanishes on the flat face. Differentiating the Poisson equation gives . Test this equation with . The boundary terms vanish because on the flat face and near the curved face. Integration by parts in the tangential direction moves off , yielding
This is the essential Caccioppoli inequality step: it involves , rather than . Expanding and using the Cauchy-Schwarz inequality and the Young inequality gives
The term is bounded by . Absorbing the first term and summing over proves the hinted estimate:
All mixed second partial derivatives are now controlled, since . The remaining normal derivative is recovered from the equation:
The Cauchy-Schwarz inequality bounds its squared L2 norm by times the sum of the squared L2 norms on the right. Including the already controlled zeroth and first partial derivatives, and taking square roots, gives the second-derivative estimate at a flat Dirichlet boundary:
The full second-order Sobolev norm is controlled without any norm of .
Use a signed power test for a Laplacian eigenfunction. For , test the Dirichlet Laplacian eigenfunction equation with . Its derivative is ; at the test function is just . Integration by parts gives the exact energy estimate
Set . It belongs to the zero-boundary Sobolev space and obeys the Sobolev chain rule:
For , the power map is with bounded derivative on the bounded range of . For , the absolute value map is Lipschitz, and almost everywhere; on the zero set, use the fact that a gradient of a Sobolev function vanishes on a level set. These facts justify the identity even when changes sign.
Apply the given Sobolev inequality with the Sobolev conjugate exponent :
Testing with a signed power raises the integrability exponent from to . Absorbing the factor into gives precisely the requested inequality.
Multiply the Dirichlet Laplacian eigenfunction equation by and use integration by parts. The zero Dirichlet boundary condition removes the boundary term and gives
The denominator is positive because the eigenfunction is nonzero. The numerator is also positive: if it were zero, the given Sobolev inequality for would force . Therefore the Dirichlet Laplacian eigenvalue has the positive Rayleigh quotient
Every nonzero Dirichlet eigenfunction in this problem has a positive eigenvalue.
Choose a constant positive Hermitian form on . Its underlying flat Kähler metric is translation-invariant and descends through the Euclidean lattice to the compact complex torus . The Dolbeault Hodge decomposition gives a unique Dolbeault Laplacian-harmonic representative of every Dolbeault cohomology class. By the Kähler Laplacian identity, these are also Hodge Laplacian-harmonic.
Choose a constant coframe on . A lifted harmonic form is a sum of constant wedge basis forms with smooth periodic coefficients. For the flat metric, the Hodge Laplacian acts coefficientwise by the scalar constant-coefficient Laplacian. Each coefficient is a harmonic function on a compact flat real torus, hence is constant: integration by parts gives . Conversely every constant-coefficient form is harmonic. Thus the harmonic representatives are exactly the translation-invariant differential forms of type (p, q).
The constant covectors form ; the constant covectors form the space of antilinear maps . Sending a constant wedge form to its cohomology class gives
The map is intrinsic and does not depend on the auxiliary flat metric: that metric proves every class has exactly one constant representative, while the inclusion of constant forms is defined directly by the quotient. The Dolbeault cohomology of a complex torus consists of constant forms, and in particular
For indices outside , both the exterior-power expression and the cohomology group vanish.
For a continuous vector semimartingale and , the Itô formula is
The Itô product rule, or stochastic integration by parts, is
The additional second-order term is the quadratic covariation.
For general semimartingales with càdlàg paths, use left limits in the integrands. The Itô formula for semimartingales with jumps, for time-independent , is
Here denotes the continuous local martingale part; the continuous bracket is that of these parts. For time-dependent add and insert in the other evaluations of . The integration by parts formula becomes
For continuous paths the jump sum vanishes, recovering the first two formulas.
Return to the mostly-minus Minkowski metric of Question 3, and write . The field is assumed real and nonzero wherever its inverse occurs. Under the usual Abelian gauge transformation , with inert, the electromagnetic field tensor is invariant but . The change in the added density is
This is not generally a total derivative. The action is not invariant under arbitrary gauge transformations. It retains residual Lorenz gauge symmetry for transformations satisfying , with the same boundary conditions imposed before and after the transformation. The action itself is undefined at ; that value can only be considered as a limiting gauge.
To vary the electromagnetic four-potential, use and . The integration by parts of both terms gives
with surface terms removed by the variational boundary conditions. Hence
For variable , the derivative must act on as well as on :
Since enters algebraically, its Euler-Lagrange equation is
For a real field this gives . Thus a dynamical gauge-fixing parameter imposes a constraint equation in field theory; it is not an ordinary propagating scalar. On configurations satisfying both equations, , , and therefore . There is no independent kinetic equation that determines .
For the momentum-space equation at a prescribed constant , use . Then , and the linear equation is
Equivalently, . If the separately varied equation is also imposed on a real classical solution, then and its nonzero modes lie on the massless mass shell. For constructing a full-field Green function, however, invert the fixed-background quadratic operator before imposing that on-shell constraint. Holding fixed and integrating over is a Gaussian functional integral; integrating over the original variable as well is a different constrained problem.
For , the Lorentzian versions of the transverse projector of a vector field and longitudinal projector of a vector field are
They satisfy , , and . The mixed-index operator is
Its inverse follows by inverting these two scalar eigenvalues. Define by . The positive-sign covariant gauge-fixing inverse is
This is the algebraic inverse of the kinetic operator. A vacuum photon propagator instead has the factor and a Feynman i-epsilon prescription for its simple and double poles. If the name is used for the vacuum two-point function, include that factor consistently in its defining source equation. A retarded Green function would use different boundary conditions at the same poles. The original plus sign corresponds to the usual covariant gauge parameter ; in particular, gives Feynman gauge and .
Contracting gives the longitudinal gauge propagator contraction
Thus the contraction is not identically zero as a function of four-momentum for any , and it is a nonzero vector at every nonzero non-null four-momentum. A particular component may vanish when . At , the displayed inverse has poles and must be interpreted using the chosen Green function prescription, rather than as a pointwise finite matrix. In the limiting Landau gauge, , the longitudinal part vanishes. A covariant photon Green function may have a longitudinal component even though physical photon polarizations are transverse.
Use Lie algebra generators with , and define the adjoint gauge covariant derivative by . Varying the gauge field strength gives . Antisymmetry followed by integration by parts therefore gives
For variations vanishing at the boundary, the Yang-Mills equations are
The Jacobi identity for the commutators of gauge covariant derivatives, using , gives the gauge-theory Bianchi identity:
In the second form of the gauge-theory Bianchi identity . The identity is a geometric consequence of the definition of curvature, not a second dynamical field equation.
For the usual non-Abelian Yang-Mills theory, the quantum running coupling becomes strong at low energies. Dimensional transmutation generates a scale absent from the classically scale-invariant equations. The expected confining dynamics and massive colour-singlet spectrum involve confinement and a mass gap, rather than freely propagating weakly coupled coloured waves. These nonperturbative effects are not captured by simply solving the classical equations; this is not a claim of a mathematical proof of the Yang-Mills mass gap.
Write the gauge condition as to distinguish it from the structure constants. Under an infinitesimal gauge transformation, . Define the Faddeev-Popov operator by . In a Euclidean convention, the gauge-fixing and ghost additions can be chosen as
A nonlocal kernel would require a double integral; a local functional instead makes a local differential operator. Choosing a local gauge functional thus retains a local gauge-fixed action and the ordinary local interaction structure of perturbative quantum field theory.
The fields and are independent Grassmann-valued fields in the Adjoint representation, Lorentz scalars with ghost numbers and . They are the Faddeev-Popov ghost fields: their functional integral produces , the Faddeev-Popov determinant correcting the volume element along a gauge orbit. They are internal fields, not physical asymptotic particles. Their statistics supply a minus sign for every closed ghost loop. No auxiliary field is needed in the displayed gauge-fixing representation.
In axial gauge, , so
With the Fourier transform convention , the ghost Feynman rules in the original normalization of are
Momentum conservation accompanies the vertex; each ghost loop has the extra minus sign. With a canonically normalized field , the vertex is instead . The gauge fixing term is quadratic and introduces no further interaction vertex.
The Euclidean quadratic kernel of is . For , its inverse is
This follows by multiplication with the kernel: the piece gives , and the piece supplies the missing . The strict axial-gauge propagator is therefore
It obeys . Every ghost attachment to an internal gauge propagator therefore vanishes in the strict limit; equivalently, on the gauge slice. Thus the interacting Faddeev-Popov ghosts decouple. Residual gauge transformations and the poles at require compatible boundary conditions and an axial-gauge pole prescription.
After Wick rotation, the corresponding Minkowski propagator has replaced by and the usual overall factor , with a compatible causal prescription. The fixed makes the gauge-fixed propagator nonmanifestly Lorentz covariant. Nevertheless, gauge invariance makes physical observables independent of this gauge-choice vector, so their perturbative predictions are Lorentz invariant. Already at tree level, contraction with conserved external currents removes all terms containing or , leaving . The quantum Ward identities give the analogous cancellation in complete physical amplitudes; arbitrary gauge-dependent Green functions need not be independent of .
The odd symmetries of a zero-dimensional polynomial model preserve both the action and the flat integration measure: their coefficients have zero superdivergence. The confining polynomial weight removes boundary terms at infinity. Thus the integral of an odd-symmetry derivative is zero, a supersymmetric Ward identity for this finite-dimensional integral in this zero-dimensional supersymmetric field theory.
Since is a holomorphic function, the barred symmetry obeys . Applying the supersymmetric Ward identity gives the requested unnormalized insertion:
One can also verify the result after Berezin integration, without invoking the symmetry terminology:
The final integration by parts uses and rapid decay. The derivative on in the insertion is essential.
A spherically symmetric stellar model has no preferred direction in the tangent plane at a fixed radius. Rotations about the radial axis exchange its two transverse velocity directions, so . This uses spherical symmetry of the velocity distribution, not just a spherical mass density.
The augmented density treats radius and relative potential as independent variables, agreeing with the physical density only along . For any regular bound spherical galactic distribution function , write and . At fixed , and . Integration by parts therefore gives
For , the velocity measure is , and . Integrating the transverse derivative gives , while integrating the radial derivative in gives . Consequently
where is fixed. These identities show that the augmented representation holds for any such spherical distribution, not only a monomial ansatz. With radial pressure vanishing at the zero-binding boundary, set
Along the physical curve, the chain rule gives . The Spherical Jeans equation then requires
Thus the Jeans moments from an augmented density are
The partial derivative holds fixed. Substitution cancels the integral of and leaves precisely in the Spherical Jeans equation. This is a representation by a chosen augmented density, not a unique solution determined by the one-variable density: different extensions off the physical curve encode different anisotropies, and not every formal extension necessarily gives a nonnegative galactic distribution function.
For the hypervirial density-potential family, differentiate the potential and use the spherical shell theorem:
Differentiating and dividing by gives
The enclosed mass tends to at infinity and to zero at the centre for every . The same result follows from , with the stated sign convention.
From now on set . The separable augmented density is . Integrating it gives
and hence
With these are the requested radial expressions. The velocity-anisotropy parameter is , and the total mean square speed is .
The local kinetic-energy density is . With the potential zero at infinity, the gravitational energy density is , where . Therefore
This is the local virial relation of the hypervirial model, a special property stronger than the global virial theorem. To compute the global energies as well, set . The substitution yields the global binding integral of the hypervirial model
The beta function integral converges for all . Thus
Restoring dimensions multiplies by .
To derive the hypervirial distribution function, let be the positive relative energy and the magnitude of the specific angular momentum. Seek for and zero otherwise. Its density is
The angular integral is . Substitution makes the radial integral
Matching the augmented-density power gives . Matching its coefficient gives the normalization of the hypervirial distribution function
All integrals are finite in velocity for , and . Because and are orbital integrals, this is a steady solution by the Jeans theorem. As checks, gives for the Hernquist model, and gives for the isotropic Plummer model.
The density and potential do not uniquely fix the stellar distribution. The preceding coefficient is fixed within the chosen separable power-law model, but nonuniqueness of spherical galactic distribution functions remains without that extra restriction. An explicit positive counterexample uses the same density. Define the Osipkov-Merritt distribution function variable and take
Rescale the transverse velocity by . Then and the velocity Jacobian contributes . Direct integration gives
exactly the same Plummer model density. This anisotropic Plummer distribution with unit anisotropy radius has , so it is distinct from the isotropic distribution while remaining positive and spherical. Its second moments differ, so it does not share the extra augmented-density and local-virial restrictions of the chosen hypervirial model.
A small-angle gravitational encounter between a test star and a field star of mass gives a transverse kick , where is the impact parameter and the relative speed. For number density , encounters in occur at rate . Independent kicks add in mean square, giving the random-walk derivation of stellar relaxation
The Coulomb logarithm in stellar dynamics has of order system size , and of order the strong-deflection scale . For a virialized -star system, , so is of order . The stellar relaxation time is the time for accumulated velocity variance to become comparable with :
Using gives . With a diameter-crossing convention , this is , conventionally rounded to
The coefficient is approximate: spatial profile, velocity averages and crossing-time convention change it by factors of order unity. The robust result is crossing times.
For and years, the estimate is about years, enormously longer than a Hubble time. Most galaxies therefore behave as collisionless stellar systems. Dense nuclei and star clusters can relax more rapidly. Negligible stellar encounters do not mean negligible collective gravitational instabilities.
Define the mass-weighted galactic distribution function by , so . Hamiltonian gravitational motion preserves phase space volume by the Liouville theorem in Hamiltonian mechanics. In the collisionless regime no encounter term redistributes stars between neighbouring phase space trajectories, so . With acceleration , the Collisionless Boltzmann equation is
For a steadily rotating galactic bar, write the inertial polar angle as . The inertial radial and angular equations are and . Substitution gives the equations of motion in a rotating frame:
Here the velocities are measured in the bar frame. Set , and define the positive-force rotating-frame relative effective potential
The centrifugal sign is positive in this relative-potential convention. The characteristics are
Consequently the rotating-frame collisionless Boltzmann equation is
For these planar equations, integrate over and interpret as the planar or vertically integrated mass density. Define . Boundary terms in velocity vanish for a sufficiently decaying distribution.
Let the velocity accelerations above be . Their velocity divergence is . Thus integration by parts in the zeroth moment contributes , producing the cylindrical continuity equation for a rotating stellar system
This is conservation of mass: change of density balances flux through the radial and azimuthal sides of a small cylindrical element. The factor is geometric; the physical phase space measure is .
For the radial first moment, multiply the rotating-frame collisionless Boltzmann equation by . Its velocity integrations are and . Hence the radial cylindrical Jeans equations in a rotating frame are
For the azimuthal first moment, the velocity integral vanishes, while . This gives
The Jeans equations are local momentum-balance equations. The second moments include streaming momentum flux and random-velocity stress, while gravity, centrifugal acceleration and Coriolis acceleration supply the frame-dependent forces. They do not by themselves close the full distribution: a stress prescription or a galactic distribution function is also needed.