PaperlayerView exact revision on arXiv

Source-linked AI summary

Flory theory for Polymers

Somendra M. Bhattacharjee, Achille Giacometti, Amos Maritan

arXiv:1308.2414v2cond-mat.stat-mech

TL;DR

This review revisits simple analytical theories for homopolymers, centering Flory theory to clarify its nuances and place it in current context. It surveys Flory predictions and extensions, including dimensional behavior, crossover, and force response.

  • Problem

    The review aims to clarify the nuances of Flory theory and place it in the current context of other mean-field theories.

    Source: §I · p. 3

  • Method

    The review synthesizes Flory theory and its variants, including scaling functions, blob-based modifications, temperature crossover, and force response.

    Source: §I · p. 3, B. Scaling functions and interpolation · p. 12, C. Modification of the Flory estimate · p. 19

  • Results

    Flory predictions distinguish swollen, theta, and attractive regimes, with critical dimensions dc = 4 for excluded volume and dc = 3 for the theta point.

    Source: A. Summary of Flory predictions · p. 17

  • Takeaways & Limitations

    Flory theory is presented as applicable to systems with long-range correlations or no relevant length scale beyond the size set by the system.

    Source: §IX · p. 37

    The review is selective rather than exhaustive, leaving many important aspects and contributions uncovered or uncited.

    Source: §IX · p. 37

Abstract

from arXiv · show

We review various simple analytical theories for homopolymers within a unified framework. The common guideline of our approach is the Flory theory, and its various avatars, with the attempt of being reasonably self-contained. We expect this review to be useful as an introduction to the topic at the graduate students level.

I. INTRODUCTION

This review centers Flory theory to clarify its nuances and place it in the current context of homopolymer theories. It introduces polymer representations and size scaling, then considers excluded-volume effects and collapse.

I. INTRODUCTION

  • Flory theory predicts critical-like power laws with dimensionality dependence without fine-tuning, and the review aims to clarify its nuances and current context.The review also covers modifications, crossover behavior, and force responses.Source: §I · p. 3

1. Freely Jointed chain

  • A freely jointed chain has independent bonds that rotate freely, and its beads and bonds may cross without penalty.The review also introduces bead-spring and continuum representations.Source: §1 · p. 4, Figure 1 · p. 4

2. Size of a polymer

  • For a freely jointed chain, the rms size exponent is ν = 1/2, corresponding to fractal dimension dF = 2.The monomers remain correlated even when the bonds are uncorrelated.Source: §2 · p. 5

B. Non-Gaussian Behaviour

  • The review asks whether repulsive excluded-volume interactions change the polymer’s size exponent, and therefore its fractal dimension, rather than merely its size amplitude..Source: B. Non-Gaussian Behaviour · p. 6

1. Simple Flory theory

  • The elementary Flory estimate balances excluded-volume repulsion against chain entropy by minimizing the free energy with respect to polymer size.The free energy is estimated for a chain constrained to have radius R.Source: §1 · p. 6
  • The Flory exponent is ν = 1 for d = 1, ν = 3/4 for d = 2, ν = 1/2 for d > 4, and ν ≈ 0.588 for d = 3.The d = 3 estimate is close to the cited best estimate.Source: §1 · p. 6

2. Collapse

  • At T > Tθ the polymer is swollen, at T < Tθ it is compact, and at T = Tθ it is Gaussian-like at least in d > 3.Attraction favors compact packing with monomer density O(1) in N.Source: Figure 2 · p. 7, §2 · p. 7

III. THE EDWARDS CONTINUUM MODEL

The review moves from a Gaussian-chain description to continuum paths and effective interactions, then uses scaling to connect polymer size exponents with Flory theory.

A. Discrete Gaussian model

  • ⟨R2⟩=db2N for the discrete Gaussian chain, whose harmonic-spring representation simplifies calculations.Its extensible springs allow rare configurations larger than Nb, so the model is most useful when highly stretched states are not central.Source: A. Discrete Gaussian model · pp. 7–8

B. Continuum model

  • The continuum model replaces discrete bonds with a curvilinear coordinate and a formal sum over polymer paths, interpreted as a limit of the discrete sum.The path integral is ill-defined if taken literally as the b→0 limit; the continuum string also need not be a unit-speed space curve, and a lower cutoff may be needed.Source: B. Continuum model · p. 8

C. Interactions: The Edwards model

  • Effective monomer interactions model solvent quality: good solvents use repulsive pair contacts, while collapse uses attractive pairs stabilized by repulsive three-body contacts.The Edwards Hamiltonian combines connectivity with excluded-volume penalties; the minimal models use L, u, and v.Source: C. Interactions: The Edwards model · p. 9

D. Green Functions

  • Green functions normalize the interacting polymer’s end-to-end distribution, distinguishing a free endpoint from configurations constrained to a fixed end-to-end vector.The fixed-R and fixed-force ensembles are conjugate, but polymer results need not be equivalent across them.Source: D. Green Functions · p. 10

IV. FLORY THEORY IN A MODERN PERSPECTIVE

  • The review frames Flory theory through scale invariance, universality, and crossover behavior, linking it to ideas developed in modern critical phenomena.It notes Flory theory’s role in revealing power laws and the role of dimensionality.Source: §IV · p. 10

A. Scaling analysis

  • Continuous scale invariance implies power laws, and polymer scaling transformations keep microscopic lengths fixed while rescaling large lengths.The scaling form suppresses small-scale dependence and describes the large-distance end-to-end distribution.Source: A. Scaling analysis · pp. 10–11
  • At the Gaussian value ν=1/2, the interaction grows under rescaling for d<4; balancing the Hamiltonian terms gives the Flory exponent νF=3/(d+2).Scale invariance of the free energy, rather than the Hamiltonian alone, is required; the vertex correction is small for polymers.Source: A. Scaling analysis · p. 11

B. Scaling functions and interpolation

The scaling function ψ describes interpolation from Gaussian to swollen-chain behavior, while the review frames the exponent as independent of interaction parameters for nonzero interactions. The mean-field free energy is obtained through a density-factorization approximation.

1. Why exponent?

B. Scaling functions and interpolation

A. Mean Field Approximation

  • The mean-field Helmholtz free energy takes the standard Flory form after density correlations are neglected through factorization.The density profile is assumed to depend on the polymer size R at large scales, not microscopic scales.Source: A. Mean Field Approximation · p. 14
  • The usually subdominant logarithmic term encodes loop-formation probability, the reunion exponent, and can control transitions such as DNA melting [37].Source: A. Mean Field Approximation · p. 15

B. Solution through steepest descent method

The paper uses steepest descent to estimate the partition integral around a free-energy minimum, then derives size behavior across good-solvent, theta, and compact regimes. The predicted exponents depend on the regime and spatial dimension.

B. Solution through steepest descent method

  • Steepest descent estimates the partition integral by expanding around the free-energy minimum R* and imposing the saddle-point condition.The method is described as systematically improvable.Source: B. Solution through steepest descent… · p. 15

1. Flory regime(α > 0, T > Tθ)

  • The exponent is ν=max(νF, 1/2); for d<4, νF=3/(d+2), and d=4 is special because u is dimensionless there.Source: a) Case 2 < d < 4 · p. 16

2. Theta regime (α = 0, T = Tθ)

  • At the theta point, ν=1/2 in d=3, while in d=2 Flory predicts ν=2/3 rather than the known exact ν=4/7.Source: §2 · p. 16
  • At the theta point, Cθ=(2γ)^(1/(d+1)) for d<3 and Cθ=1 for d≥3.Source: §2 · p. 17

3. Compact regime (α < 0, T < Tθ)

  • For α<0, matching interaction terms gives the compact exponent ν=1/d, consistent with O(1) density; v>0 helps prevent complete collapse.The Gaussian term remains necessary for convergence of the integral even though it is subdominant at large N.Source: §3 · p. 17

VI. ADDITIONAL REMARKS ON FLORY THEORY

The section reviews Flory predictions and refinements to the interpolation formula, including a fractal-counting argument and an explicit two-dimensional scaling-function calculation. It also discusses how the interpolation formula compares with numerical and experimental results.

A. Summary of Flory predictions

  • For swollen and theta regimes, interactions cease to change Gaussian behavior above their critical dimensions; below them, nonzero positive interaction strength does not affect the size exponents.The critical dimensions are dc = 4 for excluded-volume interactions and dc = 3 for the theta point; Table II summarizes predicted exponents across regimes and dimensions.Source: A. Summary of Flory predictions · p. 17, Table II · p. 17

B. Modification of the entropic term and Flory interpolation formula

  • The modified free energy adds an endpoint-placement entropy term −d ln R; Eq. (6.4) shows a loop at low 2γ, and experiments found its form for polystyrene but a constant-term variant for polymethyl methacrylate.The figure uses 2γ values 0.01, 0.038, 0.1, and 1.Source: B. Modification of the entropic term… · p. 18, Figure 4 · p. 19

C. Modification of the Flory estimate

  • Fractal counting suggests an energy term R^(2dF−d), but this conflicts with the N²/R^d repulsive-energy estimate; a refined blob-contact argument recovers the Flory estimate.The blob argument estimates contacts using the intersection dimension of two fractals.Source: C. Modification of the Flory estimate · p. 19

D. Explicit computation of the scaling function in d = 2

  • In d = 2, an exact integral calculation yields a scaling-function result consistent with the large-z behavior of Eq. (6.3).The authors note that comparison of this crossover function with experiments on polymers adsorbed at surfaces or planar interfaces would be important.Source: D. Explicit computation of the scaling… · p. 20

E. The Uniform expansion method

The section develops a uniform-expansion route to the Flory interpolation formula and sketches extensions to directed and branched polymers. It then frames temperature-driven behavior as a finite-chain crossover and treats the scaling function on either side of the theta temperature for 2 ≤ d ≤ 3.

E. The Uniform expansion method

  • Choosing an effective Gaussian elastic constant makes the first-order scheme nonperturbative; its uniform-expansion condition reproduces Eq. (6.4) when γ = 0.The condition sets the correction term in brackets to zero.Source: E. The Uniform expansion method · pp. 20–21

F. Extension of Flory theory to more complex systems

  • For directed polymers, the approach uses separate transverse and longitudinal radii and gives an upper critical dimension dc = 3; for dendrimers, minimization yields an expansion-factor equation analogous to Eq. (6.4).The dendrimer discussion specializes the repulsive term to a three-dimensional system and cites applications to branched polymers.Source: F. Extension of Flory theory to more… · pp. 21–22

VII. TEMPERATURE INDUCED TRANSITION AND EXTERNAL FORCE

  • At large but fixed N, changing temperature can describe crossover among the swollen, theta, and compact regimes; on the high-temperature side, the scaling function must recover the swollen-phase Flory exponent.The exponent relation stated is νθ + qφ = νF.Source: §VII · p. 22

A. Case 2 ≤d ≤3

  • For 2 ≤ d ≤ 3 near Tθ, the scaling function has distinct solutions Ψ+ above and Ψ− below the theta temperature.The stated scaling form generalizes Eq. (4.6).Source: A. Case 2 ≤d ≤3 · p. 23

1. Case T > Tθ.

Polymer-size scaling separates swollen, compact, and theta regimes, with crossover near Tθ. For d>3, the Flory correction matters at smaller chain lengths above Tθ, while below Tθ the γ term remains relevant.

3. Phase diagram

  • The phase diagram separates swollen (T > Tθ), compact (T < Tθ), and theta regimes according to N relative to |T − Tθ|^-1/φ, with a finite-chain crossover region around θ.The theta point is at T = Tθ and 1/N = 0.Source: §3 · p. 23, Figure 5 · p. 24

B. Case d > 3

  • For d > 3, Eq. (7.12) cannot be expressed using a single scaling variable because the terms multiplying α, γ, and higher terms also depend on N.The analysis instead assumes z = |α|/N^(d−4)/2.Source: B. Case d > 3 · p. 25

1. Case T > Tθ (α > 0)

  • For 3 < d < 4 above Tθ, N^1/2 scaling holds below the crossover length and N^νF scaling above it; for d > 4, the large-N scaling is N^1/2.The Flory correction matters at relatively small N, while self-avoidance dominates at large N; the γ term is irrelevant.Source: §1 · p. 25
  • Below Tθ, 3 < d < 4 has N^1/2 scaling at smaller N and N^νc at larger N, while d > 4 has N^1/2 scaling for N ≫ 1.In this phase, both the Flory correction and γ term are relevant depending on N; the γ term is needed for convergence and leaves Ψ N-dependent.Source: §2 · p. 25

C. Inclusion of an external force

The review extends Flory theory to external pulling forces, polymer solutions, and semiflexible chains under tension. It describes force-driven phase behavior, concentration crossovers, and the limits of the Gaussian-chain description.

C. Inclusion of an external force

  • A first-order transition separates swollen and compact phases at a critical force, which Flory theory predicts varies linearly with reduced temperature.The force–temperature phase diagram is schematic, and the approach has been used to infer DNA unzipping transitions.Source: C. Inclusion of an external force · p. 27, Figure 7 · p. 27

D. Polymer Solution

  • Interpenetrating chains form a semidilute solution through a smooth crossover, and interactions screen long-range self-repulsion so sufficiently long chains become Gaussian.The crossover concentration is where individual polymer coils begin to touch.Source: D. Polymer Solution · p. 28
  • For polymer solutions, finite chains have phase-separation critical points rather than a true theta point, while the infinite-chain theta point is tricritical.In the semidilute regime, Flory scaling gives osmotic pressure Π ∼ c^5/4 in three dimensions, consistent with reported experiments.Source: D. Polymer Solution · pp. 28–29

VIII. SEMIFLEXIBLE CHAIN UNDER TENSION

  • A stretched state competes with entropy and can arise from pulling force or bending rigidity; its configurations lie in the Gaussian distribution’s tail, where the continuum Gaussian model is inadequate.The review considers semiflexible polymers under tension and derives an interpolation between small- and large-force regimes.Source: §VIII · p. 29

A. Discrete approach

The discrete-chain treatment handles bending rigidity and force while preserving fixed bond lengths. The review then motivates a distinct continuum limit for semiflexible chains and checks the large-force approximation against an earlier asymptote.

A. Discrete approach

  • The discrete semiflexible chain enforces fixed bond length with a delta-function constraint and adds a bending-energy penalty between neighboring bond directions.A force on the endpoint is equivalent to an orientational force on each bond.Source: A. Discrete approach · pp. 29–30
  • At small force, the chain’s extension is treated by linear response using zero-force orientational correlations, which decay exponentially in the one-dimensional model.The susceptibility scales as N times the persistence length when endpoint effects are ignored.Source: A. Discrete approach · p. 30

B. Continuum limit: A detour

  • Unlike the freely-jointed-chain procedure, the continuum limit must retain the unit-tangent constraint along the space curve, replacing the Edwards Gaussian term with a bending-energy description.The discrete approach is presented as simpler while yielding the same results as the continuum formulation.Source: B. Continuum limit: A detour · pp. 30–31

C. Large f limit: detailed calculations

  • In the large-force calculation, assuming small transverse fluctuations and bf ≫ 1, the leading terms agree with the earlier large-force asymptote.The derivation reduces the transverse fluctuations to a Gaussian integral involving a tridiagonal matrix.Source: C. Large f limit: detailed calculations · pp. 31–32

D. Small force limit: detailed calculations

The zero-force chain’s mean-square size distinguishes rodlike, semiflexible, and Gaussian length regimes. The section also obtains the small-force response through an expansion about the zero-force partition function.

D. Small force limit: detailed calculations

  • The small-force partition function is expanded around zero force, where rotational invariance makes the first-order term vanish and supports the expected linear response.The expansion uses zero-force averages and neglects subdominant terms in the stated large-chain separation limit.Source: D. Small force limit: detailed… · pp. 32–33, §1 · p. 33
  • For N ≪ lp/b the chain is rodlike with ⟨R2⟩ ∼ N^2, whereas for N ≫ lp/b it is Gaussian with ⟨R2⟩ ∼ N.The persistence length remains finite for finite bending strength, and the semiflexible regime lies between these limits.Source: §1 · p. 34
  • The zero-force mean-square size grows linearly with chain length for N → ∞, consistent with Gaussian behavior at lengths large relative to the persistence length.The persistence length is the tangent-correlation length in the chain description.Source: §1 · pp. 33–34

E. An interpolation formula

The review gives an interpolation formula for semiflexible-chain extension under force and relates structure-factor calculations to end-to-end size. Its outlook presents Flory theory as a useful but selective graduate-level introduction.

E. An interpolation formula

  • An interpolation formula combines the small- and large-force asymptotes for semiflexible-chain extension, with the large-force derivation taken in the lp ≪ Lc regime.The continuum limit keeps persistence length and contour length fixed as bond length tends to zero.Source: E. An interpolation formula · pp. 34–35

F. Structure factor and end-to-end distance

  • The semiflexible-chain structure factor is expanded in powers of wavevector, with odd powers vanishing by symmetry.The second-order expansion yields an end-to-end distance consistent with the direct result and with Ref. after a mean-field translation 2lp/3 → lp.Source: F. Structure factor and end-to-end… · p. 36

IX. OUTLOOK

  • The review aims to introduce polymer-physics methods and results self-containedly for graduate students, using Flory theory as a unifying framework.It connects the Flory approach with the Edwards continuum model and field-theoretical methods.Source: §IX · p. 36
  • Although Flory theory ignores fluctuations and can give coarse predictions of phase transitions, it captures dimensionality-dependent signatures and often gives results close to correct ones.The authors note that few approaches outside rare exact solutions and full renormalization-group calculations provide dimension-dependent results.Source: §IX · p. 37
  • The review is selective rather than exhaustive, leaving out important contributions and covering polymer solutions only briefly.Its chosen emphasis includes steepest descent, interpolation formulas, finite-size crossover, and tricritical-point effects.Source: §IX · p. 37
  • The authors present Flory theory and its variants as continuing tools for studying polymer systems and other scale-invariant problems.They point to its use for systems with long-range correlations or no relevant length scale beyond system size.Source: §IX · p. 37

Appendix A: Gaussian integrals and the Hubbard-Stratonovich transformation

The appendices use Gaussian identities and the Hubbard–Stratonovich transformation, then derive the large-N end-to-end-distance distribution and compare bond-length models. The Gaussian limit applies when the bond distribution has finite variance.

Appendix A: Gaussian integrals and the Hubbard-Stratonovich transformation

  • For symmetric Gaussian interactions, correlations ⟨φ_iφ_j⟩ are related to the inverse interaction matrix, and the generalized Gaussian identity underlies the Hubbard–Stratonovich transformation.The transformation is stated for a symmetric matrix with positive eigenvalues.Source: Appendix A · p. 37

Appendix B: Distribution of the end-to-end distance in d dimensions

  • For N ≫ 1 and kσ ≪ 1, expanding the bond-distribution Fourier transform gives the large-scale end-to-end-distance distribution; in d = 3 it reduces to Eq. (2.5).The expansion uses the variance σ² of the bond distribution.Source: Appendix B · p. 38

1. Examples

  • The freely jointed chain has fixed bond length, whereas the Gaussian bond model has no fixed length; both yield the expansion with σ² = b²/d.For the Gaussian distribution, σ = b; the fixed-length example is also shown to match the expansion.Source: §1 · pp. 38–39

2. Non-gaussian case

The Gaussian approximation to polymer end-to-end statistics requires finite bond variance, while perturbation theory examines small interaction strength. The appendix presents direct and Laplace–Fourier procedures for first-order corrections, with divergences handled by analytic continuation or renormalization.

2. Non-gaussian case

  • The Gaussian large-N limit requires finite variance; a Cauchy bond distribution instead retains Cauchy end-to-end statistics, so the review excludes distributions dominated by very large sizes.Its large-|x| tail produces divergent mean and variance.Source: §2 · p. 39

Appendix C: Perturbation Theory

  • Perturbation theory expands in small z, but the resulting series is often asymptotic; its first-order two-body term represents a single monomer contact.The end-to-end-distance correction requires expanding both numerator and denominator in the interaction parameter.Source: Appendix C · pp. 39–40

1. Direct evaluation

  • The direct evaluation encounters divergent contour-length integrals for d ≤ 4 and uses analytic continuation; the perturbative end-to-end distance diverges as d → 4.Poles at d = 4 and d = 6 are identified as responsible for divergences.Source: §1 · pp. 40–41

2. Laplace-Fourier approach

  • The Laplace–Fourier approach is a grand-canonical alternative to the canonical calculation, decoupling the R and s convolutions for the two-body correction.The discussion limits itself to two-body interactions, though additional terms can also be considered.Source: §2 · p. 41

a. First order correction from perturbation theory

The first-order perturbative calculation regularizes a divergent integral and quotes its final result in Eq. (6.9). The appendix also relates polymer size scaling to spin-system susceptibility and identifies the Gaussian-chain structure factor’s scale.

a. First order correction from perturbation theory

  • The first-order perturbative calculation uses a regulator whose divergence is handled by renormalization, and its final result is quoted in Eq. (6.9).The text states that the divergence is irrelevant to computing ⟨R²⟩.Source: a. First order correction from… · pp. 42–43

Appendix D: Issue of thermodynamic limit

  • For an interacting discrete polymer, the corresponding susceptibility per spin has finite-size behavior linked to the polymer size exponent; requiring extensive susceptibility would enforce Gaussian polymer behavior.The mapping identifies polymer end-to-end distance with total spin and mean-square distance with susceptibility.Source: Appendix D · p. 43

Appendix E: The structure factor of a Gaussian chain

  • The Gaussian-chain structure factor depends on the dimensionless variable kR0, so its wavevector scale is set by the polymer’s overall size rather than microscopic scales.Source: Appendix E · p. 44

Appendix F: Exact solution of the freely jointed chain model with external force

The appendices give an exact force response for the freely jointed chain, derive a Green function for the semi-flexible polymer, and verify a mean-square displacement result using spherical harmonics.

Appendix F: Exact solution of the freely jointed chain model with external force

  • The freely jointed chain partition function can be solved exactly with an external force, yielding the Langevin-function response and a leading-order result in the small-force limit.The physical force is defined as f_phys = f/β.Source: Appendix F · p. 44

Appendix G: Derivation of the Green function for semi-flexible polymer

  • The semi-flexible-polymer Green-function derivation recovers Eq. (8.57) using Lc = Nb and the wormlike-chain relation b = 2lp.Source: Appendix G · p. 45

Appendix H: Calculation of ⟨[rz(s) −rz(s′)]2⟩

  • The spherical-harmonic calculation of the z-component displacement coincides with Eq. (8.53) when the x and y components are included.Source: Appendix H · p. 45
Loading 1308.2414v2…