After the load is removed, a reversibly stretched DNA molecule relaxes towards a thermally fluctuating coil, releasing stored free energy as heat through the surrounding fluid. There are many more disordered chain conformations than strongly aligned ones, so entropy favours retraction; viscous drag sets how rapidly it occurs. At long times an unconstrained ideal chain has zero mean end-to-end vector but a nonzero typical size, ; it does not contract to a point.
At finite temperature the chain continues fluctuating, and an anchored or otherwise constrained molecule need not regain an isotropic ensemble. If stretching broke attachments, denatured the DNA or caused irreversible structural changes, recovery can be incomplete or hysteretic. Those outcomes cannot be decided from the ideal FJC alone.
Dense DNA packaging opposes three main physical effects. The negatively charged DNA backbones repel through electrostatic interactions, moderated by electrostatic screening and counterion condensation. At close approach, hydration forces and excluded volume become important. Bending a semiflexible DNA polymer on a radius comparable with its persistence length stores elastic energy. Confinement also removes configurational entropy, and concentrating mobile counterions can impose an additional osmotic pressure and free energy cost. Their relative importance depends on geometry and ionic conditions.
A schematic lytic bacteriophage life cycle is illustrated below.
Figure 1.
Conceptual lytic bacteriophage cycle and active DNA packaging
. Details.
Attachment recognizes a suitable host; delivery brings the genetic information into the host; host resources support production of components; assembly creates progeny that are released. Some bacteriophages instead have a latent integrated phase, but that is not a necessary stage of this lytic schematic.
In many double-stranded DNA bacteriophages, a portal-associated molecular motor uses ATP hydrolysis to move DNA into a preformed capsid. Chemical free energy pays for bending, repulsion and confinement, so the opposing force rises as filling proceeds. This is active packaging rather than spontaneous equilibrium compaction. Optical tweezers can infer that opposing force from displacement of a trapped bead and monitor packaging through the changing external DNA extension. Such measurements establish the mechanical principle without requiring any virus-handling procedure. The inverse spool model used below is an idealized packing geometry, not an exact structural claim for every bacteriophage.
The bare packed-sphere estimate exceeds by approximately ; subtracting the estimated extended-rod self-energy still leaves a factor of about . The dominant missing effect is electrostatic screening: mobile positive counterions accompany the negatively charged DNA, greatly reducing long-range repulsion. Counterion condensation can further reduce its effective electric charge. Treating the genome's bare charge as uncompensated throughout and outside a sphere therefore grossly overestimates its long-range electrostatic energy.
An optical tweezers measurement integrates an imposed force over a displacement. A reversible extraction or packaging comparison concerns a free energy difference relative to the unconfined DNA and ionic reservoirs. It can include entropy, bending and hydration contributions; dissipative pulling can add further work. Extraction can also be assisted by stored internal free energy, so the sign of externally supplied work requires the experimental convention. Here the supplied positive value is compared as an energy magnitude.
Adding monovalent counterions makes the combined genome-plus-ion electric charge nominally zero. It does not make the local electric field or free energy vanish: opposite charges have spatial distributions, and confining mobile ions costs entropy.
A definite screened continuation of b replaces the bare Coulomb potential by , with and . In the Debye–Hückel approximation, a uniformly distributed fixed charge density satisfies inside, with zero fixed charge outside and matching potential and radial field at . Its interior solution is
The screened charging free energy is therefore
For , the bracket approaches one: only charges within roughly a screening length contribute appreciably, yielding the local estimate . With the same as before, gives bracket and
about times smaller than the unscreened packed-sphere self-energy. This is a specified continuum estimate of packed-state charging free energy; a full packaging free energy also subtracts the screened unpacked reference and includes ion redistribution. It should not be presented as a unique exact packaging work fixed by the problem's few inputs.
The supplied counterion count also permits an independent estimate of the entropic contribution. For ideal confined ions, , so compression from an accessible reference volume to costs
For example, approximate the original screening atmosphere as an annulus of thickness around a rod. With , . Taking gives . This useful order-of-magnitude ion-entropy model is distinct from the screened charging model; the two estimates must not simply be added as independent terms without a consistent ionic free energy functional. The rod cutoff, fixed DNA volume and screening-atmosphere geometry are rough geometric assumptions, not a microscopic reconstruction of the source sketch.
There is a further consistency limit. The nominal ion number density is . For a uniform ideal monovalent screening gas, would give at , rather than ; added salt would shorten it further. Thus the stated ion inventory and screening length cannot both be interpreted as a self-consistent dilute uniform screening gas with this dielectric. Bound counterions, nonuniform distributions and an effective phenomenological screening length can change that interpretation. Linear screening is itself only an estimate near highly charged DNA.
The specified screened-sphere charging model gives , about times the supplied experimental magnitude. It is much closer than the unscreened result. The separate ideal counterion-compression estimate gives , within about of that magnitude, but this numerical proximity does not validate its rough geometry.
Remaining differences can arise from the unpacked reference self-energy, counterion condensation, ion exchange with the surrounding reservoir, nonuniform charge density, nonlinear electrostatic screening, finite DNA thickness and the actual capsid shape. DNA bending energy, hydration forces and configurational entropy also contribute. The experiment measures a process-dependent work or inferred free energy, whereas the simple screened sphere computes just one model contribution. Neither a zero-net-charge argument nor agreement of one approximate number establishes an exact total energy.
Each Cartesian component of one link has zero mean and variance , with zero cross-component covariance. By the multivariate central limit theorem, for large on the coil scale the vector probability density is
It is normalized with respect to ; the radial probability density is . The Gaussian limit of the end-to-end distribution of a freely jointed chain is an approximation, not the exact finite- law near .
At prescribed vector separation the configurational entropy is . The ideal chain has no extension-dependent internal energy, so its Helmholtz free energy is
The externally required force is , while the chain's restoring force has the opposite sign. Thus
The last equation also follows by multiplying by the Boltzmann factor and completing the square. This entropic spring resists extension because fewer link configurations realize a stretched chain. At the same small applied force, increasing temperature reduces extension; at the same extension it increases the required force. Do not differentiate to obtain the fixed-vector Hooke law, since that radial measure includes an additional geometric multiplicity. The low-force condition is , equivalently .
In atomic force microscopy, one end of the polyacrylic acid chain is attached to the substrate and the other to a compliant cantilever tip. Controlled separation displaces the cantilever; calibrated cantilever stiffness converts its deflection into force. The molecular extension is the tip–surface separation after correcting for cantilever deflection and instrumental offsets. A quasi-static interpretation requires relaxation to be fast compared with the pulling motion.
At low force a coiled polymer has many available configurations and extends mainly by changing segment orientations. This is the entropic spring regime. As the available contour is aligned, the loss of configurational entropy becomes steep. For an inextensible three-dimensional freely jointed chain,
where is the Langevin function. The limiting extension is the polymer contour length ; the ideal model approaches it only as . Real polymer backbone stretching permits extension beyond the fitted rigid-link contour, explaining why the experimental curve lies on the more extensible side at high force.
The original trace has a sharp force drop after its peak. Finite contour length explains the steep preceding rise, but does not itself cause that drop: release of an attachment, desorption or mechanical failure ends tensile loading. The graph alone does not distinguish these possibilities. If the intact chain is released from tension, thermal motion restores a coil with characteristic end-to-end distance of a polymer chain of order , because this macrostate has many more configurations. Recoil is reversible for an intact equilibrating chain; irreversible chain scission or attachment loss can prevent repetition of the same pulling curve.
For the adsorbed lattice gas, the Helmholtz free energy is
Use Stirling's approximation with large, then differentiate at fixed . Writing surface coverage as gives the adsorption chemical potential
The occupation ratio comes from site-mixing entropy; it is essential near saturation.
Let be the number of free atoms and . In the field-free gas, both spin states have zero energy, so the one-atom partition function is , where is the thermal de Broglie wavelength. For a dilute ideal gas,
Therefore
The last equality uses . The bulk spin degeneracy is two even though the bulk magnetic field vanishes: degeneracy is not the absence of spin states.
For fixed degeneracies and ionization energy,
Therefore increases with temperature in this model. Ionization costs positive energy, so it is suppressed at low temperature by its Boltzmann factor. At higher temperature that penalty becomes less important.
Ionization also creates an additional independently translating particle and more ways to distribute the particles and their momenta. This favours the ionized mixture through its translational and mixing entropy, consistent with the prefactor and increasing ionization at high temperature.
Assume a dilute classical ideal gas of indistinguishable, noninteracting particles, with translational states dense enough to replace their sum by a continuum integral. Separate translation from all internal degrees of freedom. Then
Using Stirling's approximation in the Helmholtz free energy,
Differentiate at fixed . Since ,
Only after this derivative set and for one mole at standard pressure:
This is the translational Sackur-Tetrode equation. The factor is essential for the correct extensive entropy; replacing by a temperature-dependent standard molar volume before taking the fixed-volume derivative would give the wrong constant term.
For a linear rigid rotor in the high-temperature regime, . Its molar rotational entropy is
The rotational symmetry number must already be included in .
For identical, non-interacting molecules, the translational canonical partition function is . Stirling's approximation gives the Helmholtz free energy
At fixed , , so differentiation gives . For one mole at the standard-state molar volume,
Here the quoted molar translational quantity is per mole, and its ratio to is dimensionless. The derivative defining entropy is taken at fixed volume before inserting the standard-state volume; differentiating along a constant-pressure path would give the wrong additive constant.
The substituted molecule has the larger entropy. Increased mass and moment of inertia reduce translational and rotational spacings and make more states thermally accessible. The largest contribution here is the loss of the twofold rotational indistinguishability of the oxygen ends, contributing . It is a rotational symmetry number correction for one isotope-defined species, not an entropy of mixing between different molecules.
Use the same standard pressure and temperature, the same bond length, and the same electronic and vibrational-ground-state approximations. The translational partition function scales as mass to the power , giving the factor . For the original isotopologue, and the rotational symmetry number is two. For the isotope-substituted molecule, the two oxygen ends are distinguishable, so the symmetry number is one. Hence
The change in standard molar entropy is
so
This uses a convention omitting a separate nuclear-spin entropy. If nuclear-spin degeneracies are instead included explicitly, the spin- carbon-13 nucleus adds a further relative to the spin-zero original isotopologue; that convention must be stated rather than silently mixed with the quoted partition functions.
Under the stated vibrational-ground-state approximation and a nondegenerate electronic ground state, only translation and rotation contribute. Therefore
A temperature-independent ground-state energy, including zero-point energy, contributes no entropy. This estimate deliberately omits thermal vibrational excitation; it need not equal a tabulated value that includes all modes.
As a first estimate, neglect the small mass change and changes in geometry and vibrational populations. Removing an electron creates an unpaired spin, giving at least the spin-only electronic entropy increase
For a more specific estimate, the actual cation ground term is , rather than an orbitally nondegenerate doublet. Ignoring spin-orbit coupling, its two orbital components and two spin components give electronic degeneracy four and
The ground-state assignment is confirmed by state-resolved photoionization measurements.[ref] Spin-orbit splitting and changed vibrational/rotational constants refine this estimate; the question's supplied partition functions alone do not determine those corrections. The answer is a spin-only estimate, while additionally retains the orbital degeneracy.
The strongest donor belongs to the acid, , and the strongest acceptor to the aminopyridine, . The cocrystal pairs the strongest donor with the strongest acceptor, while the reciprocal amino-NH···carbonyl contact closes a complementary two-contact motif.
Using the simple favorable-contact score , one heterodimer gives
The self-associated comparison, normalized to the same two molecules and two hydrogen bonds, is half an acid dimer plus half an aminopyridine dimer:
Thus
In a model where bond free energy contains , this is a favorable energy change. No bonds have been artificially added to one side of the comparison: both alternatives contain two contacts per acid/aminopyridine pair.
The calculation rationalizes the hydrogen-bond preference. It is not a complete crystal free energy: packing, additional contacts, entropy and possible solvent effects also matter. Without a specified calibration, the numerical score should not be promoted to a measured cocrystallization free energy.
In water, separated benzene molecules require separate hydration environments around nonpolar surface. Bringing them together buries part of that surface and reorganizes the surrounding water. The resulting solvent contribution favours association relative to cyclohexane, where benzene was already well solvated by dispersion. This solvent-mediated contribution is the hydrophobic effect, not formation of a new covalent bond or a claim that benzene has become a conventional hydrogen-bond acceptor like oxygen.
The supplied constants imply a thirtyfold enhancement. With matched standard states,
at an illustrative 298 K. The water value is still slightly below , so enhancement relative to cyclohexane does not mean a negative standard association free energy for a 1 M standard state, nor does it imply that dilute benzene is all dimer. No universal “entropy only” explanation is needed: solvent enthalpy and entropy both contribute and depend on temperature.
Association in solution includes replacement of solute–solvent contacts, not just formation of the isolated dimer. cyclohexane solvates the benzene surface well through dispersion. Much of the London dispersion gained on bringing two benzene molecules together replaces comparable benzene–cyclohexane contacts, leaving only a modest net gain before the translational and orientational entropy penalty. cyclohexane cannot provide the complementary donor/acceptor interactions available when two waters hydrogen-bond, so water association gains a more distinctive interaction relative to its separated solvated monomers.
Using the same concentration standard state, and . For illustration at 298 K and , the ratio of the two supplied association constants corresponds to
Thus the water-molecule pair can associate more favourably in this solvent even though the stipulated isolated-pair binding energies are similar. Gas-phase potential energy alone omits both differential solvation and the entropy of association.
The clearest comparison is between the two amines, which have the same molecular formula and molar mass. isopropylamine can donate N–H hydrogen bonds to neighbouring nitrogen lone pairs; trimethylamine has no N–H donor. The resulting intermolecular association raises the primary amine's boiling point, allowing it to remain liquid at room temperature, whereas trimethylamine is gaseous. Independent phase-change data place their normal boiling points on opposite sides of room temperature.[ref]
acetone has a substantial carbonyl dipole and London dispersion forces sufficient for a room-temperature liquid, even though it cannot build a strong donor–acceptor hydrogen bond network by itself. isobutane has mainly London dispersion forces and a compact shape and remains gaseous. acetamide combines carbonyl acceptors with two N–H donors, allowing an extended, directionally bonded crystal lattice; favourable packing as well as intermolecular association stabilizes its solid state.
These compounds have similar molar masses, so mass alone cannot explain the differences. The hydrogen-bond acceptor parameter also cannot predict the phase on its own: a pure substance needs complementary donors to exploit strong acceptor sites, and freezing additionally depends on crystal packing and entropy. A stronger acceptor does not automatically imply a higher melting point for arbitrary molecules.