For discrete random variables, write for the joint probability mass function and for its marginal. The conditional entropy is the average of the entropies of the conditional distributions:
Conditional distributions at zero-probability values of can be chosen arbitrarily; they contribute nothing. Here denotes Shannon entropy, with one common logarithm base and .
For , use and sum with weight . Summing out in the first term gives
For countable alphabets, the entropy contributions are nonnegative, so the equality remains valid with extended values by summing nonnegative terms. Applying this identity successively to the prefix and proves the chain rule for information entropy:
For , the conditioning tuple is empty and the term is .
For the subset inequalities, first suppose the marginal entropies are finite, so all joint entropies below are finite by the chain rule and the assumed conditioning reduces entropy property. Put , and let be the joint entropy with coordinate omitted. Apply the two-variable chain rule in the order to get
The inequality holds because the prefix conditioning set is contained in the set of all other coordinates. Summing over and using the chain rule gives
This proves Han's entropy inequality. In terms of the normalized subset entropy, and , hence for .
Apply the same result to each -element coordinate set . It gives
When this is summed over all -subsets, each fixed -subset occurs once for each possible added coordinate, namely times. Therefore
Substitute and the analogous identity at . The binomial coefficient identity cancels the common positive factor and gives
This is the monotonicity of normalized subset entropies, obtained from the local -coordinate inequality by explicit counting.
If some marginal entropy is infinite, each subset average contains at least one subset including that coordinate. The entropy of such a subset is at least that marginal entropy, by the nonnegative conditional-entropy chain rule, so every is infinite. The inequalities then hold in the extended sense. Otherwise the finite-entropy proof applies. For independent random variables the averages are all ; repeated copies of one variable instead give , illustrating the entropy reduction caused by redundancy.
Let be the number of isolated vertices and set . For fixed , the falling factorial counts ordered -tuples of distinct isolated vertices. A specified tuple is isolated exactly when all edges incident to its vertices are absent. There are such edges, so
Since and , taking logarithms gives
Thus every fixed factorial moment tends to . We use the standard factorial-moment criterion for Poisson convergence: nonnegative integer-valued random variables whose factorial moments converge, for every fixed order, to converge in distribution to . Consequently has that limiting distribution, and its mass at each fixed integer converges. In particular,
This establishes the limiting probability of at least one isolated vertex, rather than merely a bound based on its expectation.
The Poisson distribution mass function is
It is an exponential dispersion family of order one with
Indeed, and . A generalized linear model specifies independent random variables from a response exponential family, together with a linear predictor and a link function relating it to the response mean. The Poisson canonical link is .
For the insurance counts, let class and merit be the reference levels in a regression factor. The fitted Poisson regression assumes
The insured exposure is known; its logarithm is an offset, not a coefficient to estimate. The model has eight free mean parameters and fixes the dispersion parameter at one. It assumes additive class and merit effects on the log-rate scale, hence multiplicative effects on the rate scale without a class-by-merit interaction.
Apart from constants independent of the parameters, the log-likelihood is
Differentiating gives the Poisson regression margin-matching score equations
Thus the fitted totals equal the observed totals in each class and merit category, including the reference categories by subtraction. In design matrix notation these equations are . Fisher scoring solves them iteratively; the reported three iterations describe the numerical fit, not an additional statistical result.
The regression intercept gives the baseline fitted claim rate per insured car year. The merit rate ratios, relative to merit at the same class, are
The class rate ratios, relative to class at the same merit, are
Each predicted cell rate is the baseline rate times its class and merit multipliers. The output's coefficient-to-standard error ratios are asymptotic normal Wald test statistics under the fitted Poisson regression, despite the software label “t value”; there is no estimated Gaussian residual scale in this model. All nonreference effects have very large absolute ratios.
The null Poisson deviance is on nineteen statistical degrees of freedom; adding the seven factor coefficients reduces it by . But the final Poisson deviance is still on twelve statistical degrees of freedom, and the deviance residuals include values near and . These are far beyond what a well-fitting unit-dispersion model would normally produce. The additive Poisson model fits substantially better than a common rate, but remains grossly inadequate. A class-by-merit interaction, overdispersion, or statistical dependence between claims could contribute; the summary alone does not distinguish these explanations. Its very small nominal standard errors and formal Wald tests rely on the rejected model, so improved mean structure and error assumptions are needed before treating them as reliable uncertainty assessments.
For , an exterior polynomial root of persists by continuity for sufficiently small negative real in the stability polynomial . Hence these parameters cannot be A-stable.
For , the principal polynomial root at is simple. Let denote its analytic logarithm near zero. Using the expansion in part (b), and , solve the polynomial root equation to obtain
Indeed the residual at is , and the correction in cancels it by multiplication with . Choose . This lies strictly in the left half-plane, but
for small positive . Thus , proving linear instability within the A-stability domain. This is a direct principal-root obstruction to A-stability, rather than an appeal to an order-barrier theorem alone.
Finally, at , the full stability polynomial is
It retains a double unit polynomial root for every negative , allowing growing parasitic solutions. Therefore
The canceled Backward Euler method recurrence at is A-stable, but its admissible solutions and starting relations differ from those of the original method.
For each fixed , continuity of the logarithm gives the pointwise convergence . But at , which belongs to for ,
Consequently the uniform error cannot tend to zero. In fact is unbounded as for every fixed . The pointwise limit is , and convergence is not uniform.
Let count primes for which is also prime. We prove using a quadratic Selberg sieve, then apply Abel summation.
For , each odd prime has exactly two roots modulo that prime. The Chinese remainder theorem gives the polynomial root density in a sieve
for odd squarefree integers . Put and choose the optimizing real Selberg sieve weights , supported on these , with . If are prime and , their weight is one. It is nonnegative for every other integer. Expanding the square therefore gives
The Selberg sieve diagonalization makes the main quadratic form , where
The optimal Selberg weights have modulus at most one. Since , the remainder is at most .
For completeness, the lower bound on is where the second logarithm enters. Restrict to primes and put . The full product is , since its logarithm is and the Mertens second theorem applies. Give the associated squarefree divisor probability . Prime inclusion probabilities are , so the Mertens first theorem gives
for large . The Markov inequality shows that at least half the weight has . This proves the truncated Euler-product lower bound . Consequently the twin-prime upper bound from a quadratic sieve is
Finally, Abel summation gives
The tail is dominated by . Thus reciprocal-sum convergence from a counting bound proves convergence of the first reciprocal sum. The second is smaller termwise. We have proved Brun's theorem:
Counting each twin prime just once changes at most a finite repetition and gives the same conclusion.
Use base-two logarithms, so Von Neumann entropy is measured in bits, and set . Write for the quantum relative entropy, the quantity denoted in the paper. It is unless the support of a positive operator lies within that of .
We first establish nonnegativity of quantum relative entropy, including its equality condition. Let and , and suppose the support condition holds. Put . Concavity of the scalar logarithm gives
Terms with have zero weight here. Working first with natural logarithms, the inequality yields
since the sum to at most . Equality forces for every positive-weight eigenvector, no remaining weight of outside their span, and equality in the logarithm's strict concavity. The latter makes each an eigenvector of with eigenvalue . Thus exactly when .
For the bipartite density operator,
Indeed, a projector onto the kernel of , tensored with , has zero expectation in . Positivity implies that annihilates its range: . The same argument applies to . Hence the quantum relative entropy below is finite. On this support,
Using the partial trace to evaluate the two local terms gives
This proves Subadditivity of Von Neumann entropy:
The proof also identifies equality exactly with the product state .
A nondegenerate distribution function is a max-stable distribution if, for every integer , there are and such that
Equivalently, the suitably centered and positively scaled sample maximum of independent and identically distributed random variables with law again has law . The extremal types theorem states that every nondegenerate limiting law of normalized sample maxima is a max-stable distribution. Thus this property identifies all possible limiting types, although a particular starting distribution function need not admit any such limiting normalization.
The three standard types are the Gumbel distribution
the Fréchet distribution, with ,
and the negative Weibull distribution, with ,
They are max-stable distributions, as follows respectively by taking , , and . The necessary and sufficient classification is: is max-stable exactly when for , , and one of these three types . A positive affine map sends a max-stable distribution to another max-stable distribution.
For the requested threshold equivalence, set . By independence,
If , then and the expansion gives , since . Exponentiation proves the claimed limit. Conversely, if , then : any subsequence with would have its powers tending to zero. Taking logarithms now gives . Since , including the continuous extension at , it follows that . Hence
The exceedance count has the binomial distribution with parameters . For and each fixed nonnegative integer , its probability mass function satisfies
where and . Indeed , , and . If , the bound gives the degenerate Poisson distribution at zero. Summing the finitely many terms proves, in both cases,
This derives the needed instance of the Poisson limit theorem directly.
Finally fix with and put . The assumed limit for the normalized sample maximum and the proved equivalence give . The second-largest order statistic is at most exactly when at most one observation is strictly greater than . This event identity holds even if observations have ties. Apply the preceding Poisson limit theorem with to obtain
The endpoint corresponds to and is already included. This is the case of the exceedance count limit for extreme order statistics.
For a pro-p group, the lower p-series is the descending sequence of closed characteristic subgroups
The bar means closure of the subgroup generated by the displayed powers and group commutators. Throughout this question is odd. For a finitely generated powerful pro-p group, group commutator collection gives , each is powerful, and . All these subgroups are open. A uniform pro-p group is such a group with all subgroup indices equal to . We use the standard powerful-group torsion criterion: a finitely generated powerful pro-p group is uniform exactly when it is a torsion-free group. For a uniform group of generator number , the indices are , and is a homeomorphism from onto .
We record explicitly the collection facts used for the operations below. With , one has
The second statement is the power-congruence property: the th-power map is an isomorphism on adjacent uniform layers, and collection gives for . Thus the first nonzero layer of a difference moves down exactly one step under a th power. Iterating gives the asserted equivalence and also justifies taking roots of congruences. These facts follow from powerful-group group commutator collection, using divisibility by of and the successive-layer power isomorphisms; they are the uniform-group structure results being used here.
The transported addition on a uniform pro-p group is
The product is in and has a unique th root in . For the comparison at define . The numbered proofs below establish associativity, approximate commutativity and compatibility of these operations.
The intrinsic Lie algebra of a uniform pro-p group has underlying set , zero equal to the group identity, additive inverse , and operations
Here the group commutator is . For , is the limit of for integers ; the pro- topology makes it well defined. The addition sequence is Cauchy by the compatibility proved below. The group commutator estimate puts in , so its displayed root exists and lies in . Collection of group commutators shows that successive bracket approximants differ in layers tending to the identity, giving the other limit. The group commutator identities give additivity of the limiting bracket, and the Hall-Witt identity gives its Jacobi identity. Approximate commutativity makes limiting addition commutative. The resulting lattice is a free module over of rank , and .
A powerful Lie algebra over p-adic integers, in the convention relevant here, is a Lie algebra whose underlying module is finite free over and whose bracket satisfies . For the condition is instead . For odd , define multiplication on its underlying set by the convergent Baker--Campbell--Hausdorff formula
An iterated bracket of degree lies in . The BCH coefficient bound loses at most powers of , so all higher-degree terms are integral, and their valuations tend to infinity. This explains both convergence and why multiplying a lattice with small bracket stays inside the lattice. Formal associativity of BCH passes to the convergent series, the identity is , the inverse is , and every group power satisfies . The quotients by are finite p-groups. Hence the resulting compact group is a finitely generated torsion-free group and a pro-p group, and its group commutator lies in . Thus is a uniform pro-p group, with and successive index .
The standard uniform-group/Lie-lattice correspondence identifies these constructions as mutually inverse: and . To see why the inverse limits recover the lattice operations, the addition approximant for a BCH group is : its linear term is and every higher term tends to zero. Similarly its group commutator approximant has limiting term . In the other direction, a local logarithm in the complete normed group algebra identifies the intrinsic operations with algebra addition and commutator and identifies group multiplication with BCH; the uniform power maps extend the identification from a sufficiently small power subgroup to all of . Continuous group homomorphisms respect powers, unique roots and limits, while Lie algebra homomorphisms respect every BCH term, so the equivalence also preserves morphisms.
The objective is a continuous faithful representation of a uniform pro-p group in some . The proof proceeds from the intrinsic Lie lattice to a small matrix group, then extends the resulting representation to all of . This last step matters: a representation defined only on an open subgroup does not yet establish linearity of the original group.
Let be the Lie algebra of a uniform pro-p group and put . The uniform-group/Lie-lattice correspondence identifies the group law on with BCH. The substantial finite-dimensional Lie-theoretic input is the Ado theorem: every finite-dimensional Lie algebra over a characteristic-zero field admits a faithful finite-dimensional representation. Apply it to obtain
An Adjoint representation would not suffice here, since it kills the centre; Ado supplies injectivity even in central directions.
Choose a p-adic lattice spanning . Since has finitely many lattice generators, their representing matrices have bounded denominators. There is therefore an integer such that
On this matrix lattice the series and converge and are inverse. For example , so the valuations of the exponential terms tend to infinity with this choice of ; the logarithm terms do as well. They take values in and in respectively. The formal identity holds for the convergent series.
Under the correspondence, has Lie lattice and is an open characteristic subgroup. Define
where denotes the identification of with . A Lie homomorphism preserves BCH, so the exponential identity proves . If , the matrix logarithm gives ; faithfulness of gives . Thus is a continuous faithful integral representation of .
To extend it, take the module
Values on a transversal of determine , so is a free module over of rank . These functions are continuous, since on each open coset they are given by the continuous representation . Define . The equivariance condition is preserved, and . In a transversal basis each is a permutation of blocks followed by matrices from , so it belongs to and depends continuously on .
Faithfulness can be checked directly. If , right multiplication moves the coset to a different coset, so moves a function supported on one block and cannot be the identity. If and , evaluation at gives for every , whereas equivariance gives . Since is arbitrary, , and hence . Therefore
Compactness of makes this continuous injection a homeomorphism onto its closed image. This proves the linearity of a uniform pro-p group for the full group, including its centre, and explains why taking a sufficiently small power lattice permits matrix integration before finite-index induction restores the whole group.
The Prime number theorem says as . Bertrand's postulate says that for every integer there is a prime number strictly between and .
To count S-smooth numbers, split each prime exponent into its parity and an even part. This gives the unique expression , where is a squarefree product of a subset of . There are possible , while whenever . The squarefree-part bound for smooth numbers therefore gives .
For an integer , take to be all prime numbers at most . Every positive integer at most is then S-smooth, so . Taking logarithms yields
Use base-two logarithms, so entanglement entropy is measured in bits. Write for the binary entropy, with . Taking a partial trace over any two parties removes the off-diagonal terms: the two traced strings are orthogonal. Consequently all three one-party reduced density matrices of the first pure state are
Each division into two groups is one party versus the other two. The Schmidt coefficients are and , so the entanglement entropy across each of , and is . For the second pure state, set : every one-party reduced density matrix is , and all three entanglement entropies are one bit.
For every pair , the first pure state instead gives the reduced density matrix
This is a separable quantum state, explicitly a convex combination of two product states. Its relative entropy of entanglement is therefore zero: in , choose to attain zero, while nonnegativity of quantum relative entropy excludes a negative value. The same argument with equal weights applies to the second pure state. Thus
The pairwise density matrices contain classical correlations, despite the entanglement of the global pure state; discarding a party loses the coherence between its two branches.
The entanglement entropy of independent source copies across any cut is , whereas copies of the target GHZ state have entanglement entropy . Average monotonicity of pure-state entanglement entropy under LOCC gives for exact successful forward branches. In a reversible asymptotic conversion the reverse direction similarly prevents a first-order deficit, since . For asymptotically faithful outputs, a continuity bound for quantum conditional entropy with trivial conditioning changes these inequalities only by . Hence the reversible leading-order rate is
The zero pairwise relative entropies of entanglement are consistent with this conversion: no independent pairwise resource has to be supplied or removed.
There is a qualification to the literal finite-block formulation in the source. For , the Schmidt rank of the exact state of source copies is , whereas target GHZ states have Schmidt rank . Every refined LOCC branch acts by a product of local linear operators and cannot increase this Schmidt rank. For , the displayed rate has for sufficiently large blocks when , so an exact successful reverse conversion is impossible, irrespective of its claimed probability. The exact dilution obstruction from Schmidt rank is avoided in the usual asymptotic formulation by allowing outputs approaching the desired pure state. The leading-order rate above and the exact forward protocol below remain valid; neither establishes the source's stronger exact round-trip claim. At the two states are identical. At or , the source is a known product state, the rate is zero, and local preparation supplies any desired number of source copies.
Divide metal enrichment by gas depletion in the closed-box model of galactic chemical evolution:
The logarithm appears on integrating from and :
With the gas fraction of a galaxy , this is
It gives at and increasing gas-phase metallicity as gas is consumed, independently of the time dependence of the star formation rate. For this closed initially star-free system, , so this initial-reservoir fraction equals the present gas-to-total baryonic fraction. That equality fails in an open system.
First, for , the derivative of is , while . Hence .
Choose large enough that the given ratio bound holds for every . Positivity permits taking logarithms, and the preceding inequality gives
Summing, and using the integral test for convergence from part (b), gives for
Thus . Since , the P-series converges, so the comparison test for series proves
The finitely many initial terms do not affect convergence. This power-law comparison from a refined ratio bound supplies information when the ordinary ratio test can be inconclusive.
Put , take without loss of generality, and write . The Newtonian gravitational potential is singular at the origin; all pointwise claims below concern . The zero of the Newtonian gravitational potential is the one in the paper, with an implicit fixed length unit inside the logarithm.
The Poisson equation for Newtonian gravity in cylindrical coordinates gives
The two contributions to the Laplacian are
Consequently the mass density of the axisymmetric logarithmic gravitational potential is
The coefficient of is positive. The coefficient of determines the density positivity for an axisymmetric logarithmic potential: strict positivity at every noncentral point requires and . Nonnegativity permits , where the mass density vanishes on the noncentral symmetry axis. For it is negative near that axis. Either sign of a nonzero real gives the same model, and is undefined. Setting yields vacuum away from the singular origin, not strictly positive mass density. The central cusp is locally integrable and contains no hidden point mass: the Newtonian gravitational field flux through a shrinking sphere is of order its radius.
An equatorial circular orbit requires inward gravitational acceleration equal to . Thus the circular speed is
This is a flat galaxy rotation curve. The singular center does not define an additional circular orbit at .
Figure 1.
Flat equatorial rotation curve of the singular logarithmic galaxy; the central point is excluded
. Details.
The specific orbital energy is conserved because the Newtonian gravitational potential is time independent. The axial specific angular momentum is conserved because the Newtonian gravitational potential is axisymmetric. Hence the Jeans theorem, or directly in the Collisionless Boltzmann equation, permits a stationary two-integral galactic distribution function. Using makes this galactic distribution function even under reversal of the azimuthal velocity, so it describes a nonstreaming population. A general axisymmetric population can have an odd part in , or dependence on a third integral of motion; neither is required for this construction.
At each position define . The two-integral galactic distribution function is invariant under swapping and , since only their sum of squares occurs in . It is also even in each velocity component. Swapping integration variables proves meridional velocity isotropy of a two-integral distribution, while reflecting one variable makes each mixed integrand odd. Therefore
All mean velocities vanish, so these are also the corresponding statements for the velocity ellipsoid and its centered covariance tensor. They do not require to equal either meridional second moment.
To normalize the proposed logarithmic-potential two-integral distribution, first separate the mass density into terms with the same spatial factors as the two exponentials:
Since and , the Gaussian integral and its derivative give
Thus velocity integration of the galactic distribution function yields . Matching the independent spatial factors determines
This verifies both the stationary Collisionless Boltzmann equation and the required Poisson equation for Newtonian gravity. A shift changes the normalizations to and ; the physical mass density is unchanged.
There is an important distinction between positive mass density and a nonnegative galactic distribution function. For , both coefficients above are nonnegative. For a prolate model , and one must test the sum rather than reject the model just because one term is negative. On accessible phase space,
with equality in the equatorial plane at , . The nonnegativity bound for a prolate logarithmic-potential distribution is consequently
For larger , the displayed algebraic expression still integrates to the positive mass density, but is negative at those equatorial circular orbit phase points and is not a physical galactic distribution function. This extra restriction is distinct from the earlier density-only answer.
The locally integrable logarithm has distributional derivative . By the Fourier transform of a derivative and the previous result,
For , this gives . At zero the expression must be extended as a Hadamard finite-part integral, not treated as a locally integrable function. For example, define the finite-part inverse-absolute-value distribution by
Then . Different reference cutoffs change only by a multiple of the Dirac delta distribution. The general solution of the transform equation is consequently
Indeed multiplication by annihilates only a delta term among distributions supported at zero; delta derivatives would generate additional nonzero terms. This supplies the precise interpretation of the paper's formula, and the convention-dependent constant need not be found.
For the binomial branching process, the offspring probability generating function is . The Galton-Watson extinction fixed point gives .
The printed finite-binomial upper bound is false without an additional asymptotic qualification. For and , solving the fixed-point equation exactly gives
For instance gives . More generally, at every fixed the Taylor expansion gives , again contradicting that upper bound for sufficiently small positive .
Here is the corrected binomial branching survival correction. The Poisson branching process of mean has branching survival probability satisfying . Expanding the logarithm gives
Thus , which implies the requested lower estimate . Since on , iteration of the two offspring probability generating functions gives .
For the actual binomial branching process, expansion of its exact fixed-point equation gives
If , all coefficients are nonnegative. Keeping the first term proves
In particular, for and this gives . The printed bounds also hold for the binomial branching process in an explicit large- regime: and . Indeed these conditions give and . In the preceding nonnegative series, evaluation at gives
The series is increasing, and its value at the actual branching survival probability is , so . Together with the lower bound already proved, this recovers in that regime. In particular it applies eventually to part (ii). The counterexample shows why a regime condition is needed for a literal finite- statement.
Figure 1.
Finite-binomial survival exceeds the printed upper bound, while Poisson survival lies between the corrected comparison bounds
. Details.
Set . The chain rule gives , so the transformed equation is the forced harmonic oscillator , with and . Therefore
This real solution is defined only where its logarithm has positive argument. Put and . Then , and the connected interval containing zero is .
Since , a stationary point requires . The positive branch cannot contain the oscillator minimum . It contains the maximum , attained first at . Hence
At the endpoints where , tends to minus infinity; the oscillator solution cannot be continued through those zeros as a finite real .
For the normalized second-order equation, an ordinary point is one where both coefficient functions are analytic. A singular point of a second-order linear ODE is a point that is not ordinary. It is a regular singular point if and extend analytically to ; a singular point failing that condition is irregular.
Here , . Thus zero is singular but regular, because and are analytic. Write the analytic solution as . Its constant equation gives , and for ,
Using gives
The factorial denominators show convergence for every finite .
For the other branch, integrating the first correction equation twice and imposing its two conditions gives and
Then , so and
Each correction tends to zero at the origin and its derivative vanishes at one. Although the second solution has a finite limiting value one, it has a logarithmically divergent endpoint derivative. Indeed the differential equation gives , and integration gives
A finite initial derivative cannot be prescribed for this second branch. The logarithm is consistent with the integer separation of the two exponents in the Frobenius method; this is not the repeated-exponent case.