The coil-globule transition is the contraction of an isolated flexible polymer as solvent quality worsens. A swollen chain in a good solvent has size in the Flory approximation. At the theta temperature the effective two-segment excluded volume vanishes and . A poor solvent favors a dense polymer globule, with at fixed internal density.
The competition is visible in the scaling free energy
For , elastic cost balances two-segment repulsion to give the swollen-coil law. At , three-segment repulsion stabilizes ideal-chain scaling; for , attraction balanced by the positive three-segment term gives . The dense polymer globule also has an interfacial free energy. The scaling expression explains regimes, not all fluctuation corrections or a detailed finite-chain transition curve. A finite chain crosses over smoothly; the sharp theta-point limit refers to indefinitely long chains.
A dilute-solution light-scattering experiment can follow radius of gyration while changing temperature or solvent composition. Measure at several concentrations and extrapolate toward isolated chains; check molecular-weight distribution, reversibility and the absence of growing multi-chain aggregates. dynamic light scattering measures hydrodynamic size, which is related to but distinct from radius of gyration. A falling scattering-derived size accompanied by a rising aggregate population cannot by itself establish single-chain collapse. solvent quality controls conformation through the balance of segment interactions and conformational entropy.
A Couette rheometer measures the torque required to rotate one coaxial cylinder relative to another. With inner radius , outer radius , wetted length , inner angular velocity and stationary outer wall, laminar Newtonian flow has
The shear stress at the inner wall is . In a narrow gap , , so .
Measuring torque over a range of speeds gives a flow curve: constant viscosity identifies a Newtonian regime; decreasing or increasing apparent viscosity indicates shear thinning or thickening. A yield-stress fit needs a genuine low-rate regime and must not be confused with bearing drag or wall slip. Small oscillatory strain also probes viscoelasticity, with for imposed .
Control temperature, calibrate instrument torque and end effects, and check gap dependence or roughened surfaces for slip. Keep inertia small enough to avoid Taylor vortices. The exact Newtonian profile is a calibration model; a non-Newtonian material needs its own stress–rate relation.
Figure 1.
Original Couette geometry and dimensionless Newtonian velocity and shear-rate profiles
. Details.
An amphiphile contains a solvent-compatible head and a solvent-avoiding tail. In water, aggregation hides hydrocarbon tails while leaving heads exposed. This reduces the free-energy cost of tail–water contact, competing with loss of translational entropy, headgroup repulsion and chain packing. A micelle is a finite equilibrium aggregate, not necessarily a permanently bound particle.
Equality of molecular chemical potentials controls aggregation. For a large preferred aggregation number, approaches an aggregate value , giving
Above the critical micelle concentration, added surfactant mainly increases aggregate population while monomer activity changes relatively little. Finite aggregation numbers give a smooth crossover rather than an exact bulk phase transition. The amphiphile packing parameter compares tail volume, preferred head area and chain length; small favours spherical micelles, while larger values favour less curved aggregates. These geometric criteria are approximate, not universal shape phase boundaries.[ref]
Surface-tension or conductivity changes with concentration locate an operational CMC; scattering probes size and aggregation number. temperature, impurities and salt must be controlled, and independent methods help distinguish aggregation from precipitation. In an ionic surfactant, salt screens headgroup repulsion and can change both the CMC and aggregate shape. This discussion concerns generic synthetic surfactants only.
Figure 1.
Original generic micelle schematic and an explicitly illustrative aggregation crossover
. Details.
Interpret the stated free energy as per lattice site of volume . Flory–Huggins solution theory counts polymer chains, giving the entropy term , and solvent molecules, giving . The interaction term penalizes unlike contacts when ; a simple contact model has . Connectivity reduces polymer translational entropy relative to independent monomers. This is an incompressible mean-field model with composition-independent for the calculations below.[ref]
Expanding the dilute osmotic pressure identifies the effective second virial coefficient. A good solvent has and positive excluded volume, promoting swollen chains. A poor solvent has , favouring contraction and potentially phase separation. The theta-solvent condition is in this long-chain mean-field description: the second virial term vanishes, but higher terms remain. For finite , this condition is not the exact critical point, which instead has .
Let . The osmotic pressure in Flory-Huggins solution theory is
Thus . For mass concentration and chain molar mass ,
osmometry uses a membrane permeable to solvent but retaining every polymer chain, at known uniform temperature. Extrapolate to zero concentration to remove interactions; the intercept gives the number-average molar mass. Aggregation, leakage of small chains, ions with a Donnan contribution, or a poor-solvent two-phase sample invalidate the simple molecular-weight interpretation.
The measured ratios give an exactly linear reduced-pressure fit:
The reduced osmotic-pressure intercept is . The data determine ; the measurement temperature is needed to convert this into an absolute molar mass. For an explicitly assumed room temperature , . The boxed expression is the general answer for an arbitrary experimental temperature.
Figure 1.
Original reduced-osmotic-pressure regression and its zero-concentration intercept
. Details.
At fixed temperature and solvent permittivity, the Debye–Hückel screening length satisfies , where is ionic strength. Fully dissociated sodium chloride has two monovalent ion species, giving M. Fully dissociated magnesium sulfate has two divalent species, giving
The multivalent-electrolyte screening-length ratio is therefore
The squared valences are essential: counting only the total ion concentration misses a factor of four in the screening strength per formula unit. This is the comparable linear-screening estimate, assuming unchanged dielectric response and full dissociation. ion pairs and concentrated-electrolyte nonideality can make a real divalent-salt solution depart from it; no such correction is specified.
For a long ideal chain of statistical segments of length , . At fixed repeat chemistry and statistical segment length, molar mass is proportional to . The molecular-weight scaling of an ideal-chain gyration radius thus gives
The supplied size and ideal statistics are used as stipulated; an independently tabulated coil size would not replace the given baseline.
solvent quality changes the balance of segment–segment and segment–solvent contacts. In a good solvent, excluded-volume repulsion swells a sufficiently long isolated chain, with , in three dimensions; the Flory approximation gives . At theta conditions the effective two-segment excluded volume vanishes and the ideal exponent is recovered for long chains. A poor solvent favors a compact polymer globule, whose large-chain size scales approximately as at fixed internal density. solvent and temperature can also change local persistence length and hence the prefactor. These different conformations need not have the same radius at the same molar mass.
The Stokes–Einstein relation connects thermal diffusion with dissipative mobility. For a dilute suspension of isolated rigid spheres of hydrodynamic radius in an unbounded Newtonian fluid of viscosity ,
where is translational drag coefficient, diffusion coefficient and absolute temperature. Stokes drag supplies under no-slip, continuum, low-Reynolds-number conditions. The Einstein diffusion relation is more general than that particular spherical drag law.
A short equilibrium argument explains the connection. Let a weak external potential exert force , so drift velocity is . The number flux is . A Boltzmann distribution has ; requiring zero equilibrium flux fixes . thermal fluctuations and viscous dissipation are therefore linked, rather than independent adjustable effects. The original theory also relates diffusion to measured random displacements.[ref]
On times long compared with momentum relaxation, an isotropic diffusing sphere has
single-particle tracking can measure the slope. Alternatively, dilute single-scattering dynamic light scattering measures a field correlation ; the corresponding intensity-correlation excess is proportional to . The scattering wavenumber and the distinction between field and intensity decay must be retained.
A Stokes-Einstein experimental consistency test independently measures temperature and solvent viscosity, uses a narrow distribution of inert spheres, checks decay rate proportional to or displacement variance proportional to time, and extrapolates toward infinite dilution. Subtract drift, account for wall-enhanced drag, and verify that aggregation is absent. Tracking additionally needs spatial calibration and correction for finite exposure and localization precision. At finite concentration, collective diffusion is not automatically the single-particle self-diffusion in the equation. A hydrodynamic radius is not generally identical to a polymer radius of gyration; nonspherical particles, slip or viscoelastic solvents require a different mobility model.
A viscoelastic material stores part of the work of deformation and dissipates part. In the linear regime its stress depends on the strain history, not just the instantaneous strain or shear strain rate. For simple shear, a causal relaxation modulus gives . A sudden strain then probes stress relaxation, while a constant imposed stress probes material creep; these are different experiments.
For steady sinusoidal strain use the convention and complex amplitudes. The complex viscoelastic modulus is
so . The storage modulus is the in-phase elastic component. The loss modulus multiplies the component in phase with shear strain rate; stress leads strain by . The loss tangent is when .
Multiplying stress by shear strain rate and averaging gives . Hence the dissipated energy per oscillatory shear cycle is
The elastic contribution returns its work over the cycle. A passive material has nonnegative dissipation, which fixes the sign of the loss term in the chosen Fourier convention. This complex viscoelastic modulus is a response function, not the mathematical absolute value of a complex number.
A Maxwell model places a spring of shear modulus in series with a dashpot of viscosity , so the strains add at common stress:
The Maxwell fluid dynamic modulus is
At the response is predominantly viscous, and . At it approaches the elastic plateau . The two moduli cross, and reaches its maximum, at . A Kelvin-Voigt model instead has the spring and dashpot in parallel, giving and no fluid-like zero-frequency terminal flow. Real materials often need a spectrum of relaxation times rather than one Maxwell time.
Figure 1.
Original Maxwell storage and loss moduli, with a stress–strain cycle whose enclosed area measures dissipated work
. Details.
Oscillatory rheometry measures stress amplitude and phase while imposing a small strain. An oscillatory-shear amplitude sweep first identifies a linear regime; an oscillatory-shear frequency sweep then measures and . Repeat at controlled temperature, suppress evaporation, check wall slip and sample aging, and correct instrument inertia at high frequency. Frequency-dependent structure and a nonlinear response cannot be inferred reliably by fitting one spring and dashpot to one point.
An amphiphile has a solvent-affine part and a solvent-averse part within one molecule. In water these are commonly called hydrophilic head and hydrophobic tail. Two nonbiological examples are sodium dodecyl sulfate, with a twelve-carbon tail and an anionic sulfate head, and cetyltrimethylammonium bromide, with a sixteen-carbon tail and a cationic quaternary-ammonium head. Their counterions balance the ionic head-group charge; the head and tail do not have interchangeable interactions with water.
Aggregation reduces unfavorable exposure of hydrocarbon tails to water but loses molecular translational entropy. Interfacial free energy, head-group solvation and electrostatic repulsion, chain packing and conformational entropy determine the preferred aggregate. The hydrophobic effect is a solvent-mediated free-energy balance, rather than a literal attractive bond between tails.
The amphiphile packing parameter
compares tail volume , effective head-group area and maximum tail length . A spherical core of radius obeys , and a long cylindrical core obeys . The condition gives the conventional spherical window and cylindrical window . Less strongly curved bilayers become compatible for , with planar packing near one; values above one can favor inverse curvature in an appropriate nonpolar surrounding medium. These windows are geometric guides rather than exact phase boundaries.
A larger hydrated or electrostatically repelling head increases and favors positive curvature. Screening of charged heads or a smaller head can lower and favor rods or sheets. Multiple tails increase the tail volume relative to head area and can favor bilayers. Tail length, branching, temperature and solvent selectivity alter both geometry and interaction free energy, so one cannot predict a unique aggregate from hydrophobicity alone.[ref]
In the supplied chemical-potential expression, is the equilibrium chemical potential per amphiphile, the standard free energy per amphiphile within an aggregate of size , and the dimensionless molecular fraction belonging to that size class. The aggregate molecular fraction versus aggregate number fraction distinction explains the logarithm: the aggregate count fraction is , and dividing its translational chemical potential by gives a per-molecule value. Thus
Here is the Boltzmann constant, absolute temperature, and a positive integer. The logarithm acts on a dimensionless ideal activity; a concentrated nonideal solution would require activity coefficients.
equilibrium exchange of amphiphiles gives the same per-molecule for every size. For monomers, . Substitution yields the finite-micelle aggregation-number distribution
This formula by itself does not prescribe a preferred micelle size: the size dependence of must include the competing packing and interaction terms. If their balance favors a finite size , a monomer population and a peak around that size can coexist. Dividing by would instead plot the distribution of aggregate counts.
At very low , the factors suppress large aggregates and most molecules remain monomers. As rises, a favorable aggregate class grows extremely rapidly because of this high power. The critical micelle concentration from standard chemical potentials has the large-size estimate
More precisely, the illustrative crossover condition gives
Added material above the crossover enters aggregates predominantly. Monomer concentration then increases slowly rather than being exactly pinned: its logarithm enters , while aggregate populations respond with a large factor . For finite micelles this is a sharp association crossover, not an absolute prohibition on aggregates below one exact concentration.
Figure 1.
Original illustrative molecular-fraction distributions over aggregation number below and above a micellization crossover
. Details.
The drawing uses a declared illustrative standard-free-energy model favoring aggregates near . It is not an inferred distribution or critical concentration for either named surfactant. The separate monomer peak is at ; a distribution of molecular fractions weights each aggregate by its number of molecules.
For the general physics of a nearly planar fluid sheet, let be a small-slope displacement, membrane bending rigidity and mechanical tension. At fixed topology and symmetric reference curvature, the quadratic Helfrich energy is
The bending term is the small-slope curvature-elastic energy; additional spontaneous-curvature or boundary terms require extra information.[ref]
For Fourier amplitudes , equipartition gives the membrane fluctuation spectrum under tension,
Higher temperature increases fluctuations, while greater membrane bending rigidity or tension suppresses them. Short-wavelength modes are predominantly bending-controlled; longer wavelengths can be tension-controlled. Finite sheet size, supports, confinement and imposed area set boundary conditions and infrared cutoffs; continuum elasticity also needs a short-distance cutoff. solvent viscosity affects relaxation rates, although it does not enter this equilibrium amplitude at fixed elastic parameters.
TODO: user-deferred biology. The requested importance of these fluctuations for bodily function is deferred. The preceding result covers only the independent continuum and thermal physics.
In a regular solution, ideal combinatorial entropy competes with a nonzero unlike-contact energy. For equal molecular volumes and mole fraction of species A, write mixing free energy per site as
where
Here is the coordination number and a pair-contact energy counted consistently in the lattice model. A positive penalizes unlike contacts; a negative value favors them. The logarithmic terms lower free energy on mixing. The sign of the enthalpy of mixing alone therefore does not determine miscibility.
At fixed temperature and composition, a homogeneous phase is globally stable only if no division into two compositions lowers its free energy. The common-tangent construction selects coexisting compositions , equivalent to equality of both species' chemical potentials. Local stability to infinitesimal fluctuations instead requires
Zero curvature gives the spinodal, while the coexistence boundary is the binodal. Between them a homogeneous state can be metastable and require nucleation; inside the spinodal its small composition fluctuations grow. The symmetric regular solution has
Below this contact penalty, the free energy is convex at every composition. Above it a miscibility gap is possible around the central composition. If with , heating promotes mixing, but specific interactions or temperature-dependent noncombinatorial terms can produce other trends.
For two polymers, let count equal-volume segments per chain, and now let be the segment volume fraction of A. Connecting segments into chains reduces their independent translational entropy. The binary Flory-Huggins polymer blend has
The energetic contact term remains extensive in the number of segments, while the favorable entropy per segment is reduced by and . This is the key long-chain limit of polymer-blend miscibility: a small positive incompatibility can overcome the mixing entropy. The lattice treatment is a controlled model for this argument, rather than a claim that every real blend has composition-independent .[ref]
The local stability condition and critical interaction parameter become
For equal chain lengths, . Similar segment chemistry and matched cohesive interactions can reduce the unlike-contact penalty; favorable specific interactions, such as complementary hydrogen bonding, can make mixing favorable even for long chains. Lower molar mass increases mixing entropy. temperature can help when it lowers the effective , but heating is not universally sufficient. A suitable common solvent can facilitate processing or ternary mixing, and a block-copolymer compatibilizer can stabilize fine domains; neither guarantees that the final solvent-free binary material is one equilibrium molecular phase.
An oil–water emulsion illustrates emulsion stability versus molecular miscibility. A surfactant adsorbs with its water-affine head toward water and its oil-affine tail toward oil, lowering interfacial tension and the work needed to make additional interface. Breaking a fixed oil volume into spherical droplets of radius gives total area , so creating small droplets still requires substantial area and usually mechanical work.
Adsorbed surfactants can oppose droplet coalescence through charged double-layer repulsion, hydrated-head interactions or steric repulsion from sufficiently extended groups. Adequate coverage, adsorption time, ionic conditions and film drainage affect the lifetime. A lower interfacial tension does not itself remove the bulk miscibility gap, and a protected emulsion is commonly kinetically stabilized. A thermodynamically stable microemulsion is a distinct situation requiring an appropriate composition and free-energy balance. The two liquids must not be called molecularly miscible simply because droplets disperse visibly.
For the capillary, balance upward contact-line force against hydrostatic column weight, neglecting meniscus volume and gas density. For inner radius ,
The contact angle is measured through the liquid. complete wetting means , so the contact-angle change of capillary rise gives
This follows from the stipulated wetting change, holding radius, bulk liquid surface tension and density fixed. One should not replace that premise with a general claim that grease always improves glass wetting; if contamination also changes the liquid surface tension, the original height and new angle alone would not determine the new height.
For two wide vertical plates separated by gap , use unit width. There are two upward contact-line contributions and liquid cross-sectional area , giving the capillary rise between parallel plates:
Only one principal meniscus curvature contributes between wide plates, whereas a circular tube has two. For the same fluid and angle,
If the plate gap equals the tube diameter, , the rise is half as large: the initial and completely wetting heights would be cm and cm. If instead , the heights are equal. In either geometry, changing only the angle from to zero multiplies the height by . Since the question specifies no plate spacing, the general gap-dependent relation is the unambiguous result; the half-height answer requires the equal-gap-to-diameter convention.
The Flory interaction parameter measures the excess cost of unlike segment contacts. In a simple lattice with coordination number , contact energies and temperature ,
A positive Flory interaction parameter penalizes unlike contacts; a negative value favors them. Mixing still depends on the balance with entropy of mixing, rather than the sign of alone.
For two polymers of degrees of polymerization , let be the volume fraction of species . The binary Flory-Huggins polymer blend has mixing free energy per lattice site
The logarithmic terms favor mixing. Their factors and arise because an entire chain translates as one object; connecting many segments greatly reduces the mixing entropy per segment. Hence long polymers can demix for a much smaller positive Flory interaction parameter than small molecules.
For composition-independent , the local stability criterion is
Positive curvature means stability to small composition fluctuations, zero curvature locates the spinodal, and negative curvature means spontaneous growth of fluctuations. Actual two-phase coexistence is determined by the common-tangent construction, or equality of both species' chemical potentials; its boundary is the binodal. A locally stable but globally unfavourable homogeneous mixture can lie between the binodal and spinodal.
Setting the second and third derivatives of to zero gives
For the symmetric Flory-Huggins polymer blend, and . For a polymer solution with and long chains of species , , consistent with the theta solvent limit. An effective Flory interaction parameter can depend on temperature and composition; in the simple energetic case , heating encourages mixing, but that temperature trend is not universal.
Figure 1.
Increasing the Flory interaction parameter can destroy stability of a symmetric polymer blend
. Details.
The original illustration uses stipulated . At the mixing free energy is convex, whereas at a central unstable region appears. The horizontal common-tangent construction in the latter case joins the symmetric coexisting compositions; its contacts differ from the spinodal boundaries. These are model curves, not supplied experimental data.
The critical micelle concentration is the approximate surfactant concentration above which added amphiphile predominantly joins micelles. Isolated amphiphiles expose their tails to the surrounding solvent; aggregation can reduce the unfavorable exposure, while losing translational entropy and incurring head-group and chain-packing costs. The favored aggregate minimizes the total free energy.
A mass-action model of micellization makes the concentration dependence explicit. For aggregation of monomers into one micelle, equilibrium requires
The dimensionless thermodynamic activities therefore obey
For a large aggregation number , a modest increase in monomer activity produces a steep increase in micelle abundance. In an ideal dilute approximation with the same reference concentration for both species, the amphiphile balance is
Well below the crossover almost every amphiphile is isolated. Above it, most additional amphiphile enters micelles, so the free-monomer concentration rises only slowly. It is not exactly constant for a finite , and a small aggregate population can exist below the nominal critical micelle concentration. The experimentally sharp change is an association crossover, not an assertion that aggregation is absolutely absent below one concentration.
In the original illustrative model, and after choosing the arbitrary reference concentration. The normalization is a teaching choice, not a measured critical micelle concentration. The monomer curve flattens relative to total concentration, while the aggregated fraction rises rapidly.
An experimental critical micelle concentration can be estimated from a change of slope in surface tension against the logarithm of total surfactant concentration, as interfacial coverage approaches saturation and bulk aggregation becomes important. Changes in electrical conductivity for ionic surfactants, or in light scattering, supply other probes. Different observables can give slightly different operational values because the crossover has finite width. Molecular structure, temperature and ionic conditions affect the balance of chemical potentials; screening charged head groups can reduce repulsion and favor aggregation, but no universal sign of the temperature dependence should be assumed.
Entropic elasticity. A rubber contains flexible polymer strands tied into a network by cross-links. Stretching biases the strands toward relatively extended configurations, reducing their configurational entropy. The restoring entropic force follows from the increase in Helmholtz free energy . For a Gaussian strand with statistical steps of length , its deformation-dependent free energy is approximately . Its stiffness therefore increases with absolute temperature, even if the underlying bond stiffness changes little.
In the ideal affine-network model of rubber elasticity, the shear modulus is
where is the number density of elastically active network strands. Thus heating increases the modulus approximately in proportion to , provided the network density and rubbery state remain approximately unchanged. Thermal expansion, changing relaxation times or damage can modify that simple trend.
Examples are lightly crosslinked natural rubber, based on cis-rich polyisoprene, and crosslinked polydimethylsiloxane. The crosslinks prevent the long-time macroscopic flow of an unconnected chain solution.
Density estimate and the unspecified modulus convention. With room temperature , the Boltzmann constant gives . If the stated modulus is , then
The molar value counts moles of active strands. If “modulus” instead means Young's modulus , a nearly incompressible isotropic rubber has , and the inferred value is
Thus the robust order-of-magnitude estimate is . The cross-link density inferred from rubber elasticity is usually an effective strand density; the actual junction density is for junction functionality in an ideal network without inactive loops or dangling ends. For four-functional junctions this would be half the strand density. The paper supplies neither modulus type nor junction functionality, so a unique chemical junction count cannot be extracted.
For an ideally dissociated electrolyte, the Debye–Hückel screening length is
Here is ionic number density, the charge in units of the elementary charge, and is molar concentration. Conversion between and number density cancels in a ratio. At fixed relative permittivity, .
In sodium chloride, , giving . Complete dissociation of aluminium chloride gives and , so
The factor is sixty, not ten: the squared charge and the three counterions per formula unit both matter. Using the given reference length directly, and , ,
The modest temperature increase raises the screening length by only about at fixed concentration; the large increase in ionic strength shortens it by almost an order of magnitude overall.
Under charge stabilization of colloids, the electrostatic double-layer repulsion now has a much shorter range. Attractive Van der Waals forces can then dominate at closer separations, reducing the barrier to particle aggregation in DLVO theory. Thus an electrostatically stabilized colloid is likely to be destabilized, with aggregation more likely than in the original dispersion.
This is a conditional physical prediction, not a proof of instability from a length alone. No particle radius, surface potential, attractive interaction strength or steric coating is specified. Steric stabilization could preserve stability, and trivalent-ion adsorption can alter the surface charge. Moreover, ideal screening theory is only an estimate at , where ion association, hydrolysis and activities can matter. The numerical result uses the supplied reference length and the stipulated unchanged permittivity, rather than replacing either with an unprovided calibration.
Write the membrane height above its mean plane as a sum of Fourier modes, . The supplied average is a thermal mean-squared modal height; for complex Fourier coefficients it is understood as , with for a real height field. The terms are:
In the small-slope approximation, Helfrich membrane energy contains
The equipartition theorem, with the independent real modes counted consistently, gives the stated membrane fluctuation spectrum under tension. Its units are length squared: and both have units , so has units .
At small length scales, is large and the bending term dominates:
The fluctuation amplitude is controlled mainly by the membrane bending modulus and thermal energy, rather than tension. The crossover occurs when , hence
The bending-dominated regime extends up to the membrane elastic length when scale means , or up to when it means wavelength. At still larger scales tension gives a spectrum. At the opposite end, the continuum formula requires wavelengths larger than the microscopic membrane thickness or molecular cutoff. Without numerical and , no numerical crossover scale is determined.
Molecular area from the adsorption slope. The supplied Gibbs adsorption equation identifies the surface excess number density through
The logarithm is understood relative to a fixed concentration unit or reference concentration; only concentration ratios enter its differences. The pure-water point supplies a baseline surface tension, but cannot be included in .
An estimate near the later working concentration uses its two bracketing points:
Taking room temperature as ,
This corresponds to about adsorbed molecules per . The data are not perfectly linear in , so the inferred molecular area depends somewhat on the concentration range: a fit to all five positive-concentration points gives about , while the higher-concentration region gives roughly . An effective area of order – is justified; a concentration-independent exact value is not.
Figure 1.
Surface-tension interpolation and extrapolated wetting threshold, alongside the local Gibbs adsorption slope
. Details.
Contact angle. Linear interpolation between and gives . For an ideal equilibrium contact angle measured through the liquid, Young's equation is
The supplied solid-vapour and solid-liquid values have difference . Hence
This assumes clean homogeneous equilibrium interfaces, negligible contact-angle hysteresis, and the stated solid interfacial energies applicable to the solution.
Complete wetting. The spreading coefficient is . The onset of complete wetting occurs at , or , if the solid interfacial energies remain unchanged. The last measured value is still , so this requires extrapolation beyond the table.
Continuing the last local slope in concentration, per , gives
Continuing a constant local Gibbs slope in instead gives approximately . Thus the data suggest a threshold near under continued-trend assumptions, rather than fixing one exactly. A surface-tension plateau caused by micellization, or a change in , could prevent or shift this threshold; no such extra data are given.
Capillary rise. For a circular capillary of radius , vertical surface-tension force balances the weight . Thus capillary action gives
Using , and ,
This is the height above the reservoir level in the ideal static capillary approximation; meniscus-volume corrections and contact-angle hysteresis are neglected. The positive gives a rise, rather than a depression.
Meaning of the mixing free energy. The supplied expression is the Flory–Huggins solution theory mixing free energy per lattice site, or per solvent-sized reference cell. Although the question calls it “per particle,” it cannot mean per polymer molecule: one chain occupies sites. Let the site volume be , the total number of sites be , and the chain and solvent counts be , so
where
The first term is the polymer's translational entropy of mixing. Its factor expresses that connected segments translate as one molecule. The second is the solvent's translational mixing entropy. Both are negative in the mixing free energy for , favoring mixing.
The final term is a mean-field excess contact contribution. The Flory interaction parameter measures the relative preference for unlike polymer-solvent contacts compared with like contacts. In the simplest lattice model with coordination number and pair energies ,
Positive penalizes unlike contacts and promotes separation; negative favors them. Effective real-solution values can also include non-enthalpic effects, so their temperature dependence need not be exactly a constant divided by . The subsequent derivatives assume composition-independent as in the supplied model.
Osmotic pressure and the special value one half. At fixed chain count, . Differentiating with respect to volume, relative to the pure-solvent reference, gives
The resulting osmotic pressure in Flory-Huggins solution theory is
The leading term is the ideal osmotic pressure of chains, since their number density is . At the second-virial term vanishes: effective two-segment attraction and excluded-volume repulsion balance. This is the theta solvent condition for a sufficiently long dilute chain; the positive third-order contribution remains, so a finite-concentration solution is not a completely noninteracting gas. The site-volume factor is necessary to give physical pressure units.
For , the effective pair excluded volume is positive. In a good solvent, the isolated chain expands relative to a Gaussian coil, with , or in the Flory approximation. For , attraction wins in a poor solvent and a sufficiently long chain collapses toward a dense polymer globule, with asymptotic size . Close to the theta condition, finite-chain crossover effects matter. At , to leading scaling order. This theta criterion is not the exact finite- phase-separation critical point; the lattice model has and .
Overlap concentration. The polymer overlap concentration is reached when there is roughly one isolated coil in its own pervaded volume. If the coil size is and a repeat unit occupies volume ,
With , the general scaling law is
The usual single-segment-volume estimate gives : in a theta solvent, approximately in a good solvent, and order one for compact globules. A coil has no sharp surface, so choosing a radius of gyration or a geometric coil volume changes numerical prefactors.
For the stated polystyrene at and , the conventional scaling estimate is
This is a volume-fraction estimate, not a molar concentration. It uses the implicit site-volume approximation; the question does not provide an independent repeat-unit volume to determine a more precise prefactor.
Semidilute length scale. Beyond overlap, a semidilute polymer solution is characterized locally by a concentration blob size . Suppose one blob contains statistical units. At the theta condition the subchain is Gaussian, so , while matching its internal concentration to the bulk gives . Eliminating gives the concentration-blob size in a theta solvent:
Using the supplied statistical step length and ,
The same convention gives , and the whole ideal chain has . Thus the blob is smaller than the chain. Because is only modestly above overlap, the value is a scaling estimate in a crossover region, not a precision asymptotic measurement. A statistical step length is not automatically a physical monomer diameter: if differs from , the more general estimates are and , with the missing prefactor left explicit.
Two stabilization mechanisms. Steric stabilization coats each colloid with adsorbed or grafted polymer layers. Bringing two particles close overlaps or compresses those layers, increasing their free energy through loss of chain conformations and, in a suitable solvent, increased local osmotic pressure. The resulting repulsion can keep the cores apart. A well-solvated nonionic coating can remain effective when added salt would screen electrostatic forces. Coatings can fail through desorption, poor solvent conditions or bridging, so simply adding arbitrary free polymer does not guarantee stabilization.
Charge stabilization of colloids instead uses similarly charged surfaces and their diffuse electric double layers. Their approach costs electrostatic and ionic free energy. In DLVO theory, that repulsion competes with attractive Van der Waals forces. It is particularly sensitive to salt concentration, counterion valence and surface ionization.
Bjerrum length. The standard Bjerrum length is the separation at which the Coulomb interaction energy of two elementary charges equals thermal energy:
At and , this is
It compares ionic electrostatic coupling with thermal agitation. For ions of signed valences , . The standard reference length uses , not the whole colloid charge: assigning a colloid charge does not change the solvent's tabulated .
Screening and the printed normalization. In a symmetric electrolyte, is the ionic valence magnitude, the elementary electric charge, the bulk number density of each ionic species, the vacuum permittivity, the solvent's relative permittivity, the Boltzmann constant and absolute temperature. Mobile ions rearrange around a charged surface; the Debye–Hückel screening length sets the range of the resulting weak potential.
There is an SI normalization error in the original PDF. Combining the Boltzmann distribution with the SI Poisson equation and linearizing at gives
The extra in the PDF denominator is inconsistent with those SI charge, permittivity and density conventions. Keeping the printed algebra would make its screening length times the SI value at the same ion density. No numerical screening length is required here, and the qualitative dependence is the same in either normalization. Added salt shortens the range of repulsion, allowing the attractive potential to capture particles more readily.
Interaction-potential sketches. Let be the surface separation between two spheres of radius . A useful schematic DLVO theory potential has the form
away from the microscopic contact region. Here is a Hamaker constant; excluded-core repulsion and molecular details regularize contact. The exponential is the double-layer repulsion and the inverse-gap term the short-gap Van der Waals force contribution. A fixed illustrates the changing screening range; actual charge regulation can change the repulsive amplitude as well.
At low salt there is a high positive barrier between a deep primary minimum near contact and the dispersed state at . A shallow secondary minimum can occur outside the barrier and support weak reversible association. Increasing salt narrows the repulsive range and lowers the barrier. Thermally activated primary aggregation becomes appreciable when the barrier is only a few ; the idealized barrier-free condition defines the fast-aggregation threshold associated with the critical coagulation concentration. Above that threshold, collisions can result in rapid capture in the primary minimum.
Figure 1.
Original schematic DLVO potentials as screening increases
. Illustrative dimensionless amplitudes show the barrier and outer attractive minimum; the near-contact primary well extends below the plotted range. These curves do not determine a material-specific coagulation concentration. Details.
Aggregate structures. Well above the critical coagulation concentration, the sticking probability is close to one. In diffusion-limited cluster aggregation, the first encounter attaches particles or clusters at accessible outer branches before they can explore many alternative configurations. The result is an open, tenuous structure. Typical ideal three-dimensional cluster models have mass fractal dimension .
Near the threshold on its lower-salt, slow-aggregation side, a remaining barrier makes most encounters unsuccessful. Reaction-limited cluster aggregation then permits repeated approaches and orientations before attachment, commonly giving denser clusters with . Since particle count scales as , the larger exponent means more particles at a given large cluster size. This explanation does not require already formed bonds to rearrange: unsuccessful encounters themselves alter which configurations eventually bind.
The side of the threshold matters. At the ideal barrier-free threshold itself, or just above it, aggregation can already be diffusion limited; proximity to the critical coagulation concentration alone does not universally imply a denser reaction-limited structure. The quoted dimensions are typical of irreversible approximately isotropic cluster models, not exact values for every colloid. Strong bond rearrangement or anisotropic attractions can change this comparison.
There is a source error to resolve first. For a drag coefficient defined by , the velocity relaxation time is . The mass and drag coefficient in the PDF's exponential are interchanged: is dimensionless, whereas is not. The consistent mean relaxation law is therefore .
A mean relaxation law alone does not determine a variance. Use the usual thermal Langevin equation, with a bath at temperature and the fluctuation-dissipation relation for a Langevin particle:
Here the calculation is for one component in Cartesian coordinates. With fixed initial position and velocity, integration gives the mean displacement
Integrating the stochastic part of the velocity once more gives the centered displacement
The delta correlation reduces its variance to a single integral:
This is the integrated Ornstein-Uhlenbeck displacement for a prescribed initial velocity. Its short-time variance is ; at long times it is . Thus the diffusion coefficient is
If the requested mean-squared fluctuation instead refers to the conventional equilibrium mean squared displacement, the initial velocity is also averaged over its thermal distribution, with and . Its contribution adds , giving
This has the ballistic limit and the same diffusive limit . For fixed nonzero , the uncentered mean square instead equals the first boxed variance plus . In three independent spatial dimensions, the centered variances add, and the long-time result is . Stating the initial-velocity convention is necessary to distinguish these finite-time formulas.
Diffusion and mobility. For a dilute isolated sphere in a Newtonian fluid, Stokes's law gives the drag magnitude , where is the dynamic viscosity and is the hydrodynamic radius. The mechanical mobility is therefore ; the translational hydrodynamic mobility matrix is .
To relate mobility to diffusion, put the particles in a weak external potential energy . Their drift speed is . For number density , the sum of drift and Fick's first law fluxes is
At thermal equilibrium, the Boltzmann distribution gives , and the net flux is zero. Substituting into yields . Thus the Stokes–Einstein relation is
Here is the Boltzmann constant and is absolute temperature. The assumptions include a no-slip approximately spherical particle, a continuum fluid, a small particle Reynolds number, and negligible hydrodynamic interactions with other particles or boundaries.
Gravitational sedimentation. The buoyant force must be included, since the mass of displaced water is comparable to the protein mass. Using gives
so the buoyant mass in the stated spherical model is
The drag coefficient is . At terminal speed, , so
Ignoring buoyancy would give and would not use the supplied water density. This calculation treats the stated radius as both the hydrodynamic radius and the radius determining displaced volume; a real hydrated protein can require a separately measured partial specific volume.
Electrical drift. In the specified simple electrophoresis model, terminal force balance is . Therefore
where is the electrophoretic mobility and the elementary electric charge. For the given positive beta-lactoglobulin charge,
in the direction of the electric field. This exceeds the gravitational sedimentation speed by about .
Separation time. The two populations have equal drag coefficients because their radii are equal. Their difference in drift speed is therefore set by one elementary charge:
Starting with coincident band centres, their separation is . Hence
At that time the original band has travelled and the mutant . Their masses do not enter this ideal electrical force balance: equal-radius particles with different masses but the same charge would have equal electrophoretic mobility in this model.
Practical discrimination. A one-charge difference is enough for the stipulated separation, so the method is sensitive to mutations that change effective charge or friction. It cannot distinguish every mutation: substitutions that preserve charge and hydrodynamic size need not separate. Buffer pH changes protein ionization, and realistic electrophoresis also involves the ionic atmosphere, electro-osmosis, walls or gel interactions, and finite band widths.
For example, the Stokes–Einstein relation at gives . For an initially narrow band, diffusive broadening of electrophoretic bands gives
at the calculated separation time. Thus the centres are about diffusion standard deviations apart, rather than forming infinitely sharp bands. The stated centre-separation criterion is achievable in the model, but actual resolution also depends on the loading width and other dispersion mechanisms.
Figure 1.
Predicted band centres and diffusion widths after 508 seconds
. Original Gaussian drift-diffusion calculation for initially narrow bands, using the stipulated equal-radius model. Details.
Equivalent centrifugal acceleration. For the original protein to attain its electrical drift speed by sedimentation, its buoyant centrifugal force must equal :
Consequently
This is the acceleration at the protein's position in the rotating solution; no rotation rate can be assigned without a centrifuge radius. If the intended comparison is instead the relative separation speed of the two bands, use , since their displaced volumes are equal. Matching then gives . In that separation experiment the heavier original protein moves faster centrifugally, whereas the more highly charged mutant moves faster electrically.