[Paper Review] Statistical Mechanics of Semiflexible Chains: A Meanfield Variational Approach

D. Thirumalai, Bae‐Yeun Ha|ArXiv.org|May 20, 1997

Protein Structure and Dynamics7 references21 citations

TL;DR

This paper introduces a meanfield variational approach to model semiflexible polymer chains, replacing the local tangent vector constraint with a global average constraint to enable tractable statistical mechanics calculations. It accurately predicts end-to-end distance distributions and elastic response under tension, showing quantitative agreement with simulations and experiments on DNA, though it underestimates persistence length in strong nematic fields due to failure to capture symmetry breaking.

ABSTRACT

We describe a simple meanfield variational approach to study a number of properties of intrinsically stiff chains which are appropriate models for a large class of biopolymers. We present the calculation of the distribution of end-to-end distance and the elastic response of stiff chains under tension using this approach. In the former example we find that the simple expression almost quantitatively fits the results of computer simulation. For the case of the stiff chain under tension we recover analytically all the known limits. We obtain quantitative agreement with recent experiments on the stretching of DNA. The limitations of our approach are also discussed.

Motivation & Objective

• To develop a tractable statistical mechanics framework for intrinsically stiff polymers like DNA and biopolymers.
• To address the intractability of enforcing the local tangent vector constraint in the wormlike chain model.
• To calculate conformational properties such as end-to-end distance distribution and elastic response under tension.
• To test the method’s accuracy against simulations and experimental data on DNA stretching.
• To identify limitations in cases involving broken symmetry, such as strong nematic fields.

Proposed method

• Replace the local constraint ||u(s)|| = 1 with a global average constraint on the tangent vector magnitude.
• Use a variational free energy functional to approximate the true statistical mechanics of the system.
• Apply the stationary phase approximation to derive effective persistence lengths and correlation functions.
• Incorporate external fields (e.g., nematic or tensile) via quadratic interaction terms in the action.
• Solve the resulting equations analytically in various limits (weak and strong field) to extract physical observables.
• Validate results against known analytical limits and simulation data for end-to-end distance and elasticity.

Experimental results

Research questions

• RQ1 How can the statistical mechanics of semiflexible chains be simplified while preserving key physical properties like persistence length and end-to-end distance distribution?
• RQ2 What is the elastic response of a stiff chain under tensile force, and how well does the meanfield approach reproduce known analytical limits?
• RQ3 How does a nematic field affect the conformational properties of a semiflexible chain, particularly the persistence length?
• RQ4 Why does the meanfield variational approach fail to capture the exponential increase in persistence length under strong nematic fields?
• RQ5 Can the method quantitatively describe experimental DNA stretching data using only contour length and persistence length as parameters?

Key findings

• The meanfield variational approach yields a simple expression for the end-to-end distance distribution that agrees almost quantitatively with simulation results.
• The method analytically recovers all known limits of the elastic response of stiff chains under tension.
• For DNA stretching, the model achieves quantitative agreement with experiments using only contour length and persistence length as adjustable parameters.
• In the weak nematic limit, the method predicts a linear increase in persistence length along the field direction by a factor of (1 + 1/3 gl₀), consistent with previous work.
• In the strong nematic limit, the approach incorrectly predicts a finite increase in persistence length (factor of 2), while exact results show an exponential growth as exp(√(glₚ)), due to failure to capture instanton effects.
• The method breaks down in broken-symmetry regimes because it cannot describe tunneling between degenerate minima in the effective potential, highlighting its limitations in such cases.