Uniform expectations require nonempty finite sets; assume this throughout the analytic argument. For real-valued functions on , define the box norm by
All four variables are sampled as independent random variables with replacement, so repeated coordinates are included. The defining fourth power equals
If it vanishes, every summand vanishes. In particular, the terms give for each , so . Absolute homogeneity of a norm follows directly: .
For the triangle inequality, prove the mixed box Cauchy-Schwarz inequality. Write
Separating the two averages and applying the Cauchy-Schwarz inequality in gives
Expanding the square and reversing the order of finite sums,
by a second Cauchy-Schwarz inequality. Thus . Expand into its sixteen multilinear terms and apply this bound to each term. Their upper bounds sum to . Taking fourth roots yields the triangle inequality. Together with definiteness of a norm and absolute homogeneity of a norm, this proves the box norm is a norm. Complex functions require conjugates and are outside the printed real-valued convention.
For the bilinear correlation bound for the box norm, use the normalized L2 norm: let and , . The Cauchy-Schwarz inequality in gives
The remaining factor is . By the Cauchy-Schwarz inequality in , its magnitude is at most
Taking square roots proves
For the tripartite graph, assume its parts are nonempty and write for its adjacency indicator functions. Let , so . The normalized triangle count is
The constant-degree hypothesis is exactly for each . Therefore the contribution of the constant is
For fixed , apply the bilinear correlation bound for the box norm with and . Since these are indicator functions, and , the fraction of adjacent to . Averaging and applying the Cauchy-Schwarz inequality in gives the sharper triangle counting with one box-uniform pair and constant opposite degree estimate
Multiplying by yields
The middle bound also handles or , when the triangle count is zero.
To connect explicitly with the suggested expansion, put , . Both have overall mean zero, and for every . Consequently ; all terms with constant reduce to . The terms containing combine into the single error just estimated. This is where the exact degree assumption is used. If any part is empty the triangle count itself is zero, but the printed uniform expectations and edge density of a bipartite graph values would be undefined; the nonempty-parts convention must therefore be stated rather than dividing by zero.
Write and . A Jacobi field is a vector field along the geodesic satisfying
The curvature convention is that used in Question 1. A geodesic variation is a smooth map with whose -curves are affinely parametrized geodesics. Torsion-freeness gives . Differentiating and commuting covariant derivatives gives the Jacobi equation for .
For the converse, let and . Choose a curve with and . parallel transport along it identifies its tangent spaces with . Set , so and . The geodesics with initial data give a variation . Smooth dependence on initial conditions and compactness of the original time interval ensure that this family is defined for all after shrinking , even if is incomplete. Its variation field has the initial values , so uniqueness for the linear Jacobi equation identifies it with . This proves the realization of Jacobi fields by geodesic variations.
The endpoint-vanishing pointwise normal fields form a vector space, since their conditions and equation are linear. If is nonconstant, the map is injective because and zero initial derivative force the zero solution. Differentiating gives . Hence the endpoint-vanishing normal Jacobi fields have dimension at most n-1 bound is
For a constant geodesic, , and both zero endpoint values force , so the bound still holds.
On the unit round , let . Its speed is . The constant ambient vectors orthogonal to give independent parallel normal fields . Since on normal fields, the fields satisfy the Jacobi equation and vanish at both antipodal endpoints. They attain dimension .
For the last clause, interpret a closed geodesic as a nonconstant smoothly periodic geodesic. A constant loop has length zero and cannot be shortened. Parametrize on at constant nonzero speed. parallel transport once around it fixes . Because the manifold is orientable, this transport preserves orientation; on the normal space of dimension , which is odd, it lies in . An odd-dimensional special orthogonal transformation has a fixed vector: nonreal eigenvalues pair with their conjugates, real eigenvalues are , and determinant one in odd dimension forces an eigenvalue .
Transport such a nonzero normal fixed vector around the loop. It gives a periodic parallel normal field , with and . Its jets also agree at the seam, so it is smooth as a field on the parametrizing circle. For small , is a smooth variation through closed curves, providing their homotopy to .
For energy , the first variation vanishes at the closed geodesic. The permitted second variation of geodesic energy has no endpoint term for this periodic variation, so
Thus for small nonzero . The Cauchy-Schwarz inequality gives , while constant speed gives . Consequently the instability of a closed geodesic in positive even-dimensional curvature yields
The deformation need not remain a geodesic or an embedded curve.
A generalized conformal metric is locally , where is nonnegative and measurable and transforms as a length density under a holomorphic change of coordinate. Its area is . Use the usual convention of using locally rectifiable paths for a path family , and put . Then
This is extremal length. Zeros and isolated singularities of an admissible density are allowed; requiring a smooth strictly positive Riemannian metric would unnecessarily restrict the definition. Line integrals have their extended nonnegative values; if an arbitrary family is supplied, use its members that are locally rectifiable paths. An empty path family has infinite infimal length, whereas a family containing a constant path has extremal length zero.
Both numerator and denominator scale quadratically when is multiplied by a positive constant. The coordinate transformation of the area element makes the quotient unchanged under conformal equivalence.
The normalization relevant later is worth deriving. On the conformal cylinder , let contain the loops going once around it. For the horizontal loop at height , Cauchy-Schwarz inequality gives
Integrating in shows . The constant density achieves equality, since every winding-one loop has Euclidean length at least one. Hence
Here is the height divided by circumference, the conformal modulus of an annulus; the reciprocal is used for the family joining its boundary components.
Work on as a probability space, with Lebesgue measure of total mass one. Let be the filtration generated by the level- half-open dyadic cells together with the separate null cell . Define the dyadic slope martingale by
The endpoint may be assigned any value, since it is a null set for Lebesgue measure. The two child slopes average to their parent slope, by telescoping the two increments. Therefore almost everywhere, so this is a martingale. The Lipschitz condition gives everywhere except possibly at the freely chosen endpoint, where we take zero.
Apply the Lp martingale convergence theorem with . Its limit has almost everywhere and in L1 norm as well, by Cauchy-Schwarz inequality. Choose a measurable representative of and set it to zero on any exceptional null set; it is then a bounded measurable function on the entire interval.
Set . Telescoping at the grid points shows that is the linear interpolation of on the dyadic grid. The Lipschitz condition gives . Also
The two uniform limits coincide, giving the absolutely continuous function representation
The chosen bound holds for every after the null-set modification; the integral identity holds for every simultaneously.
Put and . The candidate is continuous, adapted, nondecreasing, starts at zero and has integrable terminal value. We prove is a martingale, rather than assuming the desired quadratic variation of a stochastic integral.
The assumed boundedness of also makes a true martingale on every finite horizon. Indeed, apply its local martingale identity at a localizing sequence; bounded stopped squares and the Fatou lemma show that is integrable. The stopped expressions are then dominated by , so the dominated convergence theorem gives the unrestricted conditional identity on that horizon.
First let be a bounded simple predictable process. On each interval with coefficient , an increment of is . The defining quadratic variation identity for , together with the martingale property of , makes
a martingale for . Expanding gives a term and the square term above; their conditional expectations yield the desired identity on that interval. Integrability follows because the variables are square-integrable; concatenating the finitely many intervals gives a martingale.
For general bounded predictable , choose bounded simple predictable processes in , which is possible since they generate the predictable sigma-algebra. Set and . The Itô isometry and Doob L2 maximal inequality imply in ; the Cauchy-Schwarz inequality gives in and
Pass the conditional expectation identity for to the limit. Thus is a martingale. Uniqueness in the defining quadratic variation property now gives
Using the usual zero-boundary Sobolev space convention,
where weak derivatives define and is the space of test functions. In the Gaussian free field convention, the same notation often denotes the Dirichlet energy space, the completion in the gradient norm alone. On bounded domains the Poincare inequality makes the two definitions equivalent; on unbounded domains one must distinguish them. The following gradient-pairing argument applies in either setting whenever the energy completion is realized as weak functions.
Identify with a vector subspace of by zero extension of H01. Approximating by test functions in shows that this is an isometric embedding in the inhomogeneous Sobolev norm and also in the Dirichlet inner product norm. In particular, arbitrary irregularity of causes no additional boundary term.
Define
These are weakly harmonic Sobolev functions. For , integration by parts in the weak sense gives . If , choose converging to in , or in energy for the homogeneous convention. The Cauchy-Schwarz inequality gives
Hence the orthogonality of supported and harmonic Dirichlet functions is
Both are linear vector subspaces. In the inhomogeneous convention they are closed: the first is the isometric image of a complete space, and the second is the intersection of the kernels of the linear functionals . No spanning assertion is needed. The PDF contains this orthogonality statement; the TeX has badly corrupted it into an assertion about openness.
Use the normalization
for the Lasso. A centered sub-Gaussian random variable with parameter satisfies for all ; for a noncentered variable apply this definition to . With fixed , independence of the errors gives
The exponential Markov bound, optimized over separately for the two signs, gives . Take and apply the union bound over columns:
The specified positive tuning parameter requires ; at its formula is zero. For small the probability lower bound can be negative and hence uninformative. Although centering makes the centered errors dependent, means , so the preceding proof uses the original independent errors, not independence after centering.
The Karush-Kuhn-Tucker conditions, using the subdifferential of the L1 norm, are
Consequently, for nonempty , on ,
This establishes the required active-set inequality without any rank condition on .
For , a sufficient Compatibility condition for the Lasso is
Indeed, for , comparison of the two Lasso objective values and the score bound on give the Basic inequality for the Lasso
It follows that satisfies the Lasso cone condition and . Thus , including the trivial case. This proves the slightly stronger , and in particular the requested . If , the same basic inequality forces on ; no nonempty-support compatibility constant is needed.
For the support-size argument, use the stated prediction bound with constant . Set . By the Cauchy-Schwarz inequality and the definition of the sparse maximum eigenvalue,
Combining with the lower bound and squaring gives for every nonempty . If and , choose of size ; this contradicts the strict inequality defining . If , proves the assertion directly. Hence in both cases.
When , minimality and monotonicity of the sparse maximum eigenvalues give
When , the first result gives , so the same final bound holds. Thus the finite-index conclusion is .
There is a genuine domain omission in the printed last assertion: the definition of only makes sense for , so is undefined when . This case can occur: take centered orthogonal columns with , , , and . Then for every available and none exceeds , so . A universally defined replacement, from and monotonicity, is
Alternatively, explicitly extend the definition by for , including infinity. Under that added convention the printed final expression is meaningful and follows from this replacement bound. Neither convention nor finiteness should be silently assumed.
Let and be Radon-Nikodym derivatives. Use the following normalizations of the Kullback-Leibler divergence, total variation distance and unnormalized Hellinger distance:
We use and assign infinite Kullback-Leibler divergence if , where denotes absolute continuity of measures. Common domination alone also suffices for the inequalities with the extended-value convention. The negative part of the logarithmic integral is integrable: on , . Thus the extended-value integral defining the Kullback-Leibler divergence is well-defined. These expressions do not depend on which common dominating measure is chosen. The Kullback-Leibler divergence is generally asymmetric and therefore fails a defining property of a metric.
For completeness, the integral formula for total variation distance follows by taking : since , the positive and negative parts of have the same integral, and this set attains the supremum.
Put , by the Cauchy-Schwarz inequality. A second application gives
The logarithmic inequality in the hint implies for . Apply it to on :
Multiplying by and integrating proves
If on a set of positive -measure, the right side is finite but the Kullback-Leibler divergence is infinite, so the inequality is immediate. Thus
This total variation–Hellinger–relative entropy inequality requires the stated Hellinger distance normalization for its first constant. With the alternative convention , the correct comparison is and . In particular, and have .
Here is a testing form of the Le Cam two-point lemma. Suppose two parameter values in a metric space have separation , and let be the corresponding observation laws. Every estimator satisfies
To prove it, classify by the closer of the two parameters, breaking a tie in favour of one. Let . On , the triangle inequality implies , while on it implies . Hence the sum of the two probabilities in the display is at least
At least one is at least half the sum, proving the claim. It also gives a lower bound for the maximum expected metric loss.
For absolute-error loss, or any loss given by a metric, there is a stronger expectation form. Take densities relative to and use the triangle inequality under their overlap:
Therefore the Le Cam lower bound under absolute-error loss is
Both arguments are valid for randomized estimators as well, by adjoining the independent randomization to the observation; doing so does not change total variation distance.
For the normal distribution model, assume the usual nondegenerate convention and use the full product observation laws. Choose
For one observation, its Gaussian log likelihood ratio is
Under the first law, and . Thus the Kullback-Leibler divergence between normal distributions is
For the independent and identically distributed random variables, the joint log likelihood ratio is the sum of the individual ones, so the Kullback-Leibler divergence of the two full statistical samples is . The comparison just proved implies . Applying the expectation form of the Le Cam two-point lemma gives
This is a Gaussian location minimax lower bound under absolute-error loss. The constant may depend on the fixed known , but neither on the estimator nor on . If a signed nonzero scale were used, the same result holds with in place of . If zero noise were allowed, would estimate exactly, so a positive lower bound would be false; positive variance is necessary.
For a finite irreducible Markov chain satisfying detailed balance, define the directed-edge capacity . For every ordered pair choose a path from to using only positive-capacity edges; take the empty path for . Let count occurrences of the directed edge . The canonical paths Poincare bound says that
is a valid Poincare inequality for a reversible Markov chain constant. Simple paths give ; allowing repetitions requires the multiplicities shown.
To prove the theorem, expand the variance using two independent samples from :
Telescoping along and applying the Cauchy-Schwarz inequality gives
Interchange the finite sums, then use the definition of :
Consequently . All edges here are directed, so the in the Dirichlet form of a Markov chain is retained consistently.
The reference kernel printed in the hint has total mass , so it is not a Markov kernel. We explicitly replace it by the normalized random transposition shuffle: choose two positions independently and uniformly and swap them. Its probabilities are
These sum to one, and its ordinary spectral gap is , the value intended in the hint. This value can also be certified by the Aldous spectral gap theorem: the single-label generator has rate between distinct positions, so on centered functions it acts as multiplication by . Both shuffles have the same uniform distribution on a finite set as stationary distribution.
We use the Canonical paths comparison theorem in the following precise form. For two finite reversible Markov chains with the same positive stationary distribution , route each directed transition of along a positive-capacity path of . If
then and . Indeed, telescope each difference along its path, apply the Cauchy-Schwarz inequality, and sum with weights . The load definition gives the Dirichlet form of a Markov chain inequality; the common variance denominator then gives the spectral gap inequality. Identity transitions use empty paths.
Put . A transposition is routed using the word
of length . Each generator occurs at most twice; write its occurrence count as . For a fixed directed adjacent edge and each occurrence of in a word, there is exactly one starting permutation whose translated word crosses this edge at that occurrence. This follows because right multiplication by the prefix is a bijection of the symmetric group. Hence the uniform stationary factors cancel in the comparison load, leaving
The crude bounds , , and imply . Therefore
Only pairs with contribute, and there are of these. Keeping this restriction gives the stronger and . Neither comparison bound needs the exact adjacent-shuffle spectral gap; the normalization repair of the reference kernel is essential.
The Cameron-Martin space of a Gaussian random variable in a Banach space is most conveniently defined using the first Gaussian chaos. Work over the real separable Banach space , with continuous dual space . For a centered Gaussian random variable in a Banach space, define
The set inside the closure is already a linear subspace. Every is a centered Gaussian random variable. The Bochner integral defining exists: Fernique's theorem states that for some , , and then Cauchy-Schwarz inequality gives .
The map is injective. Indeed, implies for all . Density of these observations in then implies . Define the Reproducing-kernel Hilbert space by
The injectivity makes the inner product unambiguous, and the completeness of makes a Hilbert space. The embedding into is continuous since
For put . The reproducing property is
This defines the Banach-space Gaussian RKHS even for a degenerate Gaussian measure.
For , write . There is a measurable coordinate on , obtained as a limit in L2 space of continuous linear observations, and has normal distribution . The Cameron-Martin theorem for a Gaussian measure, in its real separable Banach space form, states that if is the probability law of , then the probability law of is an equivalent probability measure to precisely when . For such a shift,
For the two laws are mutually singular measures. The coordinate need not be a continuous functional on , and the sample need not belong to .
Apply this theorem to the shift . With , it gives
Choose the measurable coordinate so that on a symmetric set of full -measure. This is possible by taking an almost-surely convergent subsequence of the approximating linear observations and intersecting its convergence set with its negative. Alternatively the joint law of is invariant under simultaneous negation, directly from those observations and their mean-square convergence. Since and the centered Gaussian measure is symmetric, the two exponential integrals with signs and agree. Consequently
because . The symmetric Gaussian translation lower bound is therefore
Only symmetry and Borel measurability of are required; convexity is not needed.
There is one degenerate flaw in the last printed claim. The separable Banach space with has dense in , and all positive-radius balls have probability one. Nevertheless for every . Thus the statement needs the additional hypothesis . If the course convention excludes the zero space, this hypothesis is already implicit.
Under this necessary hypothesis, density of implies that contains a nonzero vector. Given , scale that vector to obtain with , and choose . The triangle inequality gives
The centered closed ball is a symmetric Borel set. Apply the symmetric Gaussian translation lower bound to get
Hence the Gaussian norm distribution is strictly increasing when is nonzero. This proves the strict increase of a Gaussian norm distribution without assuming that its distribution has a density or that spheres have zero probability.
Write , , and , with as required by coercivity and by the printed inverse power. The Tikhonov regularization with a coercive penalty operator normal equation gives
This also holds for inconsistent data in , because and therefore . Taking the inner product with gives
For any with , put and use Cauchy-Schwarz inequality to bound the right-hand side by . Since , completing the square in yields
The positive root of this quadratic inequality is at most . Infimizing over proves the approximate source bound for coercive Tikhonov regularization, in fact slightly more strongly:
Injectivity gives . Given any positive tolerance, approximate by for some finite-norm ; for all sufficiently large it is admissible in the infimum. Thus
No attainment of that infimum is needed.
Because , a finite sum implies . Thus the printed weighted one-norm on l2 equals the extended nonnegative sum on every . Unbounded weights may make that sum infinite for some ; membership in alone is not a finiteness assertion.
For each , is continuous in the l2 sequence space, and . If in norm, then for every ,
Taking the supremum over proves
It is also weakly lower semicontinuous by the same finite-coordinate argument, and convex and proper because zero and finite sequences have finite penalty.
At a finite-penalty base point , coordinate variations give the subdifferential of a weighted one-norm on l2:
These conditions are sufficient by summing the coordinate supporting inequalities. Since , an infinitely supported would require infinitely many , contradicting . Thus a finite-penalty point can have empty subdifferential; a Bregman base with an actual subgradient must be finitely supported.
For any such , and the Bregman distance simplifies to
Each summand is nonnegative; the sum may be infinite. The pairing converges by Cauchy-Schwarz inequality. This definition does not silently assume a subgradient exists at every .
A pointed cone has no nonzero line through the origin: equivalently . If , its quadratic form vanishes on the nonnegative orthant. Testing the coordinate vectors gives . Testing then gives
Since is a symmetric matrix, every entry is zero. Thus .
The identity matrix gives an interior point. For a symmetric perturbation with operator norm , the Cauchy-Schwarz inequality gives
This open norm ball around lies even in the positive semidefinite cone, hence in . Consequently and the interior is nonempty.
More generally, the interior of the copositive cone consists exactly of strictly copositive matrices. Positivity on the compact nonnegative unit sphere has a positive minimum and persists under small perturbations; a zero there is destroyed by an arbitrarily small negative multiple of .
The Dirichlet Laplacian eigenfunctions and positive Dirichlet Laplacian eigenvalues on the interval of length two are
A Sturm-Liouville eigenfunction expansion gives , with . Thus all normal modes oscillate precisely when .
To justify a spatially uniform bound for arbitrary appropriate data, assume the usual wave equation energy class and , or smoother compatible data. The energy method conserves
The sharp Poincare inequality is . For , this makes coercive with constant ; for , coercivity is immediate. The one-dimensional Sobolev representative satisfies by the fundamental theorem of calculus and the Cauchy-Schwarz inequality. Hence the conserved energy bounds uniformly in both variables.
At , the smooth solution is unbounded. At , is unbounded. Therefore
The bound depends on the initial data; it is not a single bound for all arbitrarily rescaled solutions. The energy-class hypothesis supplies pointwise meaning and excludes undefined rough-data interpretations.
Use the real Hilbert space , the zero-boundary Sobolev space, with norm . The Poincare inequality makes this equivalent to its usual Sobolev space norm. Define the symmetric bilinear form and continuous linear functional
The Cauchy-Schwarz inequality and the sharp interval Poincare inequality imply
Thus is a bounded bilinear form and a coercive bilinear form, and is bounded because is square-integrable. The Lax-Milgram theorem supplies a unique weak solution:
This is obtained from the differential equation by integration by parts; the test functions have zero boundary trace. Conversely, the weak solution satisfies in the sense of distributions. Since , it belongs to and solves the equation almost everywhere.
Equivalently, the Ritz method minimizes
over . Indeed , since the mixed term is . This proves existence and uniqueness of the minimizer directly from the weak solution, as well as its equivalence to the variational problem.
For a fully discrete linear partial differential equation evolution, write in a specified discrete norm. Finite-time stability of a numerical method means
with independent of the allowed meshes. For a multistep method, augment the state with its time-history values and include a stable starter. Bounds may grow as when the continuous problem has growth; demanding decay for all time would be a stronger assertion. The treatment of initial data, forcing, Dirichlet boundary conditions, and the norm is part of the hypothesis, not a detail supplied by an interior calculation.
For constant coefficients on an infinite or periodic uniform grid, Von Neumann stability analysis substitutes the Fourier mode . The discrete Fourier transform and the Parseval identity turn a uniform multiplier estimate into a discrete L2 norm estimate. In a one-step scalar scheme, proves contractivity, while gives finite-time stability. In a multilevel scheme, every root of a polynomial in its amplification polynomial of a multilevel finite difference scheme matters, as does uniform control of the associated companion matrix.
For the Forward Euler diffusion scheme with ,
Thus for all frequencies exactly when
Sufficiency follows since ; necessity follows from the highest-frequency mode on even periodic grids. This gives the usual parabolic mesh restriction . For the Backward Euler diffusion scheme,
so every is stable. The Crank--Nicolson method has
again of modulus at most one for all , but stiff modes approach rather than zero. These are PDE manifestations of the A-stability and L-stability distinctions for time integration.
For advection , assume and set . The upwind finite difference scheme is , with
It is contractive when . This also has a direct maximum-norm proof: each newest value is a convex combination of two old values. The resulting Courant–Friedrichs–Lewy condition expresses that the numerical domain of dependence covers the physical one. For negative , the upwind direction must be reversed. By contrast, Forward Euler method time stepping with the centered first finite difference has and for nonzero . Under fixed nonzero refinement it is unstable, since a fixed nontrivial mode grows geometrically over steps. This is not a claim of instability under every imaginable mesh coupling: instead bounds its spurious finite-time growth, since .
The leapfrog advection scheme illustrates a multilevel subtlety. Its amplification polynomial of a multilevel finite difference scheme is
For , both roots of a polynomial have unit modulus and separation at least . A uniformly conditioned eigenbasis of the companion matrix then proves power boundedness of a two-level Fourier scheme for arbitrary history data. At and a grid admitting , the roots coincide on the unit circle, and a Jordan block produces growth. Thus the endpoint fails the ordinary arbitrary-history root condition for a multistep method; checking only that both moduli equal one misses the instability. Restricting the starter or filtering the parasitic mode is an additional hypothesis.
The energy method handles variable coefficients and finite boundaries for which Fourier stability analysis may be unavailable. For the Backward Euler diffusion scheme with homogeneous Dirichlet boundary conditions, define . The discrete summation by parts identity is
Take the inner product of with . The elementary identity gives
This proves unconditional contractivity including the actual boundary treatment. More generally, a dissipative operator gives of operator norm at most one: with , dissipativity implies , and the Cauchy-Schwarz inequality gives . In finite dimensions this also proves invertibility. For the second-order backward differentiation formula, the BDF2 discrete energy identity used in question 4 controls both time levels and proves unconditional diffusion stability. A repeated root strictly inside the disk is harmless here; the energy argument supplies the uniform bound without a singular eigenvector formula.
A second technique is eigenvalue stability analysis of a finite difference method. Under the method of lines, a time integrator advances by . If is a normal matrix, the scalar linear stability domain criterion on all controls its powers exactly in the corresponding L2 norm. For the centered Dirichlet discrete Laplacian, its eigenvalues lie between and zero; intersecting this interval with a time integrator's stability interval determines its mesh restriction. A uniform bound on the diagonalization of a matrix is needed for nonnormal diagonalizable systems. Eigenvalues alone can be misleading. For example,
has both eigenvalues inside the disk for , but the upper-right entry of is . Taking , , this grows like . Hence stable scalar eigenvalues do not give mesh-uniform stability. Boundary closures can introduce precisely the extra growth that an interior Fourier symbol overlooks, so boundary stability of a finite-difference method must be checked separately.
For nonlinear spatial discretizations, a useful replacement for Fourier stability analysis is a convexity argument. If the Forward Euler method map is nonexpansive in a chosen norm for , the second-order strong stability preserving Runge-Kutta method
is also nonexpansive under the same restriction. Indeed, for two inputs, the triangle inequality gives . Its Taylor expansion is , proving order two. Such strong stability preserving Runge-Kutta methods transfer suitable forward-step bounds without relying on a linear spectral calculation.
Finally, consistency of a numerical method explains what stability buys. For a linear Hadamard well-posed problem, the Lax equivalence theorem equates convergence of a consistent discretization with its stability, under the stated approximation-space and norm hypotheses. Directly, if the error satisfies , iteration gives the discrete Duhamel principle
Thus vanishing normalized local truncation error and initial error yield convergence. With a nonlinear Lipschitz step estimate, the discrete Gronwall inequality plays the same role. A stable but inconsistent stencil, such as the literal defective formulas in questions 3 and 4 under their usual refinement, is not rescued by any amplification bound. The decisive checks are a mesh-uniform evolution bound, a valid boundary closure, and consistency with the actual PDE.
Differentiating the grand potential gives
A reversible adiabatic change in an isolated gas keeps and fixed. To infer constant , one must check that this ratio is one-to-one, rather than merely assuming it.
Put , , and for . Integration by parts gives and , and
The strict inequality is the Cauchy-Schwarz inequality for the positive weight , since is not constant on its support. This proves strict monotonicity of Fermi gas entropy per particle and hence is constant. Then gives constant, and gives
This adiabatic scaling of a Fermi gas with power-law density of states concerns positive temperature; the zero-temperature limit is a separate limiting case.