Each state has exactly neighbours: choose one occupied coordinate and one unoccupied coordinate to swap. The normalized proposal distribution is therefore
The swap relation is symmetric and the number of neighbours is constant, so . The discrete uniform distribution used as the prior distribution is also constant on , with mass . Both factors cancel in the Metropolis–Hastings acceptance probability.
Let and . The displayed multivariate normal distribution presupposes that the covariance matrix is a positive-definite matrix, so is a symmetric matrix. Expanding the difference of quadratic forms gives
Consequently
The truncation at zero is essential: a log acceptance ratio can be positive, but the logarithm of an acceptance probability cannot.
For the particular proposed swap, let be the coordinate changing from zero to one and the coordinate changing from one to zero. Then , and the untruncated scalar simplifies further to
Only two rows and four entries of are needed for this evaluation. This is Metropolis sampling on fixed-cardinality subsets.
Solvent quality. The Flory interaction parameter measures the free-energy cost of polymer–solvent contacts relative to like contacts. In the simplest lattice model,
where is the coordination number and the values are contact energies. For lattice-site volume , the effective two-segment excluded volume is . Its sign determines the leading tendency of a long dilute flexible chain.
For , : a good solvent favors a swollen coil because segment encounters are effectively repulsive. For , the leading attraction and repulsion cancel, giving approximately ideal chain statistics. For , effective attraction favors contraction in a poor solvent, eventually producing a compact polymer globule under sufficiently poor conditions. This does not mean that every finite-chain solution with is automatically phase-separated; that also depends on composition and chain length.
The condition is the theta temperature condition for that polymer–solvent pair. For an energetic model , it occurs at . More generally, if , it occurs when . It is not the absence of all intermolecular interactions and does not require the polymer and solvent to be chemically identical.
Chain dimensions. Interpret as the root-mean-square end-to-end distance of a polymer chain. Let denote the supplied statistical step length, calibrated to the same count as the degree of polymerization. At the theta condition, random-walk statistics give
For a good solvent, the Flory approximation balances entropic stretching and two-segment interaction:
Minimization gives , hence
Thus the familiar good-solvent scaling is when its nonuniversal solvent-quality amplitude is suppressed. With the additional lattice estimate , the more explicit form is . This swollen-coil estimate is used away from the narrow theta crossover; it must not be extrapolated to predict as . The crossover occurs when is of order one. Nearer theta, the ideal-coil size supplies the appropriate limit.
Numerical sizes. The theta value is
For , the plain good-solvent scaling gives
If the solvent-quality amplitude is retained with , it instead gives
The distinction is an explicitly stated modeling convention: the question does not give a microscopic site volume or an absolute swollen-coil prefactor. Both estimates describe a coil larger than the theta coil. Here , so a swollen-coil estimate is reasonable. Neither result is an exact microscopic prediction from the supplied data alone.
Overlap concentration. The polymer overlap concentration corresponds to one chain per coil volume. Using the scaling convention gives
The first is chain number density; the second counts repeat units. Thus the repeat-unit molar concentration is approximately . If the lattice-site volume is taken as , this is equivalently
For a mass-concentration estimate, a polystyrene repeat unit has molar mass about , so this scaling convention gives
There is no sharp coil surface fixing an exact numerical prefactor. For example, using the often adopted geometric definition , with the ideal chain radius of gyration , gives . The invariant content is one chain per coil volume and the theta scaling ; the radius and coil-volume convention must accompany a numerical value.
Osmotic pressure. In the given lattice expression, is a dimensionless polymer volume fraction, not a dimensional number density, and is the mixing free energy per lattice site. Introduce to obtain a physical osmotic pressure. If is the number of sites and , differentiation of with respect to the solvent count at fixed chain count gives
Balancing this change in solvent chemical potential against the pressure contribution gives . For composition-independent ,
and therefore the osmotic pressure in Flory-Huggins solution theory is
Expanding the logarithm as gives
The first term is with now denoting chain number density: a whole chain supplies one translational contribution. The sign of the quadratic correction reproduces the solvent-quality classification. It is positive for a good solvent, zero at the theta temperature, and negative for a poor solvent. This is the same two-segment interaction that swells or contracts a single coil. At theta, a positive third-order contribution remains; the chain is approximately ideal in its large-scale conformation, not a gas of segments with no interactions at all.
First, planar Brownian motion started away from zero almost surely never hits zero. To see this without confusing it with neighborhood recurrence, stop in an annulus with inner radius and fixed outer radius . The preceding logarithmic calculation gives probability of reaching the inner circle first. A hit of zero before the outer exit would require all these inner hits, so its probability is zero as . On every finite path interval the continuous path is bounded; a countable union over integer outer radii and integer times proves the all-time assertion.
Consequently the integrand is locally square integrable: on every compact time interval the nonvanishing continuous path has a positive minimum modulus. Write with independent standard real coordinates. Complex quadratic covariation satisfies , so analytic Itô's lemma has no second-derivative drift. For one likewise has , and . The product rule therefore gives
Its initial value is , proving, simultaneously for all times by continuity,
Thus and is the continuous accumulated change of complex argument of the path about zero, measured in radians. It is a winding angle, without reduction modulo . A single-valued logarithm on the whole punctured plane is unnecessary; the path has a continuous stochastic logarithm of planar Brownian motion.
At the transverse coordinates are independent Brownian motions with diffusion coefficient . The cylinder is infinitely long, so the axial coordinate cannot affect its first contact time. By rotational invariance of Brownian motion, the mean first-passage time depends only on the transverse radius .
The Kolmogorov backward equation for a mean first-passage time is . Here , so the radial boundary value problem is
The absorbing boundary condition for a diffusion is imposed on the circular cross-section. At , the solution must extend smoothly over the axis. Multiplying by and integrating gives
Smoothness at the axis excludes the logarithm term, and fixes . Consequently
This also verifies the finiteness assumption behind the Kolmogorov backward equation. Let be the first contact time and apply the stopped Itô formula to the bounded nonnegative function just obtained. Its bounded derivatives on the disk make the stopped stochastic integral a martingale, giving
By monotone convergence, has finite expected value; continuity and dominated convergence then make the final term vanish as . Hence . In particular, the mean first-passage time is on the axis and tends to zero at the adsorbing wall.
For a radial function in three dimensions, the Laplacian is , while the imposed radial drift contributes . Thus the radial generator of an SDE diffusion is
The hitting probability for the inner sphere is harmonic for this diffusion generator in the shell, with value one on the inner sphere and zero on the outer sphere. At , the boundary value problem is therefore
Integrating gives . The two boundary conditions determine
Both arguments of the logarithms are dimensionless. The result decreases from one to zero across the shell and lies strictly between zero and one for .
For a direct probability interpretation, stop at the first contact with either sphere. The stopped process is a bounded martingale by Itô's formula. The shell has positive inner radius, so its drift is bounded and its diffusion is nondegenerate; exit occurs almost surely. Optional stopping and dominated convergence therefore give
The logarithm appears because the inward drift cancels one of the two geometric terms in the three-dimensional radial Laplacian, leaving the same radial diffusion generator as a planar Brownian motion with diffusion coefficient one.
Ensemble and applications. A grand canonical ensemble describes a system exchanging both energy and particles with a large reservoir. Its fixed control parameters are , while microstates with different have probabilities
Weak coupling allows the system's energy levels to remain meaningful; equilibrium fixes the reservoir temperature and chemical potential. The grand potential is , with now the mean population. It generates entropy, pressure and mean particle number by differentiation. This ensemble is convenient for gases in contact with particle reservoirs, independently occupied quantum levels and surface adsorption; its particle-number fluctuations are part of the description, not an error in counting.
Classical single-level limit. For the single-level grand partition function in the dilute limit, put . Empty and singly occupied levels contribute and . In the dilute classical limit , higher occupancies contribute only at order or beyond:
This truncation does not impose the Pauli principle on a classical gas. A fermion state has exactly , a boson state has , and Maxwell-Boltzmann occupancy weights give ; they all agree to this order. In particular, , with the classical Maxwell-Boltzmann result exactly .
grand potential and chemical potential. With , summing the single-level logarithms and using the stated density of states gives
Substitute , for which
Combining this with the density-of-states prefactor gives . Thus
Differentiating at fixed ,
Similarly , so . This is the grand-canonical chemical potential of a classical ideal gas. The condition states the classical regime. The supplied density of states counts one internal state; extra internal degeneracies would have to be included consistently rather than inserted at only one stage.
One adsorption site. Take the free-vapour energy reference as zero and the binding energy as . A site is either empty or singly occupied; the latter has vibrational partition function . Its adsorption-site grand partition function and mean occupation are
The chemical potential of an adsorbed species on equivalent sites follows by inverting the occupation ratio
The factor is the vacant-site probability. The binding reference and the energies included in must be chosen consistently, so that vibrational zero-point energy is not counted twice.
Isotherm. The vibrational factor in a Langmuir equilibrium constant is fixed by adsorption equilibrium, . Neglecting vapour internal degrees of freedom and writing its number density as gives . Therefore
This is the Langmuir isotherm: initially linear, increasing monotonically and concave down, approaching unit coverage. Half coverage occurs at . The coefficient has inverse-pressure units. It assumes equivalent independent sites, at most one molecule per site and no lateral interaction or multilayer adsorption. The molecular mass appropriate to the vapour determines .
Figure 1.
Original dimensionless Langmuir coverage curve with half-coverage and limiting behavior
. Details.
A real uptake measurement must separate adsorbed molecules from free gas in the apparatus. Langmuir's original mica experiment compared a mica-containing bulb with a similar blank and inferred uptake from calibrated gas volumes and pressures. His stated blank contribution and bath-temperature drift limited the inferred adsorption; a saturating curve alone does not establish ideal equivalent sites.[ref]
Classical translational statistics. Write the number density as and the thermal de Broglie wavelength as . Quantum exchange is negligible when wave packets rarely overlap:
The equivalent state-counting statement is that the mean occupation of each thermally relevant one-particle state is much smaller than one. The equations in this question count one internal species; with equally populated spin states the relevant density per species is .
Let be the dimensionless quantum gas fugacity, with the lowest one-particle energy chosen as zero. Both quantum distributions can be written
If , then for all , recovering Maxwell–Boltzmann statistics. Conversely, low occupation even at the spectrum's bottom requires for bosons or for fermions; either condition implies
A single high-energy state can be sparsely occupied even in a degenerate gas, so the condition must hold across the relevant spectrum, not just for one arbitrarily chosen energy. Also, classical translational statistics does not imply that all internal excitation energies or other quantum effects can be ignored.
Grand partition function, density and pressure. Noninteracting particles have independent occupations of each one-particle state. Put . Summing the allowed occupations gives
Therefore
The grand potential is for this homogeneous gas. Replacing the sum by and using gives
The mean number is , giving
For bosons, the grand sum requires when the ground energy is zero. A macroscopic condensate at the limiting must be separated from the continuum integral in the condensed phase. The following small- calculation is in the normal, nondegenerate regime, so that complication does not enter.
Density expansion and inversion. Expand the occupation as a geometric series:
The Gaussian wavevector integral for its th term is
It follows that
For the inverse series, put and try . Substitution gives , hence
The leading equality directly connects the density and occupancy criteria for classical behavior.
Pressure and its first quantum correction. Expanding the logarithm in the pressure expression gives one extra factor compared with the density series:
Insert the inverse series. To second order in ,
Thus
The first quantum correction at fixed density and temperature is
It is negative for bosons and positive for fermions, and vanishes in the classical density limit.
Virial coefficient. Define the number-density virial expansion by
The coefficient has units in this convention. Comparing terms gives
This is the exchange contribution to the second virial coefficient. It arises without a pair-interaction potential; the underlying gas is noninteracting. A molar-density convention would multiply this by the Avogadro constant. For equally populated internal species, exchange occurs within each species and the coefficient in terms of total density is smaller by a factor .
Why the two signs differ. Bosonic exchange symmetry enhances the probability of nearby identical bosons, whereas fermionic exchange symmetry suppresses nearby identical fermions in the same internal state. In the dilute regime, the normalized equal-time pair detection correlation is
The exchange part follows by squaring the Gaussian one-particle coherence factor . At zero separation this gives for bosons and for one-component fermions. Integrating the exchange pair weight provides a direct coefficient check:
Bunching therefore acts like an effective statistical attraction in the dilute equation of state, lowering the pressure needed at fixed density. Exclusion acts like an effective statistical repulsion, raising it. Neither description asserts an actual interparticle force. Different internal fermion species do not exclude each other from a common spatial location in the same way, explaining the importance of the species convention.
The logarithm requires , and is given. Set , so and the chain rule yields . Dividing the ordinary differential equation by therefore gives
Since , we obtain , including the case . Thus the general positive solution is . The initial condition gives :
For the separate decay requirement, precisely when : if , the result is , and if , it grows. Consequently the full family with the stated limit is
This latter family does not also satisfy ; the solution selected by that initial condition grows. The identically zero function is excluded by the original logarithm.
For , the chain rule gives and . Applying it twice to the radial function gives
Using , Laplace's equation becomes
Two integrations now give
These are the general solutions of the specified radial function form on an annulus or a punctured plane. If the solution must extend continuously to the origin, the logarithm term forces , leaving only a constant function.
The logarithm laws give
For , both and . The expression is therefore at most . Consequently
The Big O notation concerns eventual behaviour; the zero value of at is irrelevant.
The logarithm identity rewrites the expression as . Since , we have , and for ,
Thus the Big O notation bound holds with :
Repeated squaring computes by successively squaring. The binary representation of the exponent specifies which of these values to multiply. For modular exponentiation, reduce after every multiplication; only a number of steps proportional to the logarithm of the exponent is needed.
Here and . Successive squares give
Therefore
For sufficiently large positive , the logarithm identities give
After division by , the right-hand side tends to one. Thus gives the asymptotic order in Big theta notation and Big O notation.
The logarithm is strictly increasing, so . Since almost surely,
Thus the logarithm of a Pareto random variable has an exponential distribution with rate . Its probability density function value at the single endpoint can be chosen arbitrarily without changing the law.
For natural logarithms and sufficiently large , take the logarithm of the positive expression:
Set . Then . The exponential function grows faster than : for example, its power series gives for , so . In particular, eventually . Exponentiating yields
For logarithms to a fixed base , the natural logarithm of the left side is , still smaller than eventually. The big O notation conclusion is therefore independent of that base choice.
Use . Then , , and the secant method gives
Writing , the next step is
Thus the requested approximation is to three significant figures. All intermediate logarithms should be retained at higher precision before this final rounding.
True for every fixed . Extract from inside the logarithm:
Choose an integer . The nonnegative terms of the exponential series give for , whence
Since for , the displayed expression tends to zero and is therefore . The constant and cutoff may depend on the fixed ; no uniform assertion over unbounded is being made.
Eliminate the constraint by writing , with . The logarithms cancel, so the objective becomes . Its derivatives are and . Thus the unique stationary point is the unique global minimum, at . Restoring gives
As a direct global check, the arithmetic-geometric mean inequality gives , with equality exactly at the same point.
Integrate the solids’ molar volume contribution from bar and use ideal-gas water chemical fugacity . For pure solids,
The logarithm is dimensionless, and can be multiplied directly by pressure in bars. Therefore the dehydration equilibrium with ideal-gas fluid gives
Hence K, or C. Do not also add a constant gas-volume pressure term: the water contribution has already been integrated through . Ideal-gas behavior at kbar is a prescribed model approximation, not a realistic equation of state. The solid volume correction is distinct from it.