The bending energy of a semi-flexible polymer chain and the polygons of the polymer chain

Pramod Kumar Mishra

REVIEW 4 major objections 5 minor 12 references

The bending energy of a semi-flexible polymer chain and the polygons of the polymer chain

T0 review · 4 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The minimum energy to bend a lattice polymer chain into a closed loop is $k_B T \ln[2(d-1)g l_p]$, set by dimension, fugacity, and persistence length.

desk verdict The central formula is asserted and self-confirming; the only quantitative support contradicts it. read the letter →

arxiv 1908.05586 v1 pith:75DMPN62 submitted 2019-08-15 cond-mat.soft cond-mat.stat-mech

Pramod Kumar Mishra This is my paper · ORCID

classification cond-mat.softcond-mat.stat-mech PACS 36.20.Ey05.40.Fb

keywords semi-flexiblepolymerpolygonbendingenergypersistencelengthrandomwalklatticemodelpartitionfunctionloopclosure

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that for an ideal semi-flexible polymer chain on a square or cubic lattice, the minimum energy needed to bend the chain into a closed loop is a simple logarithmic function of space dimensionality $d$, step fugacity $g$, and persistence length $l_p$: $\epsilon_b = k_B T \ln[2(d-1) g l_p]$. The argument works by singling out 'polymer polygons' --- closed lattice walks whose four arms each contain $l_p$ monomers, so the perimeter is $4l_p$ --- and counting how many of them return to the grafted site: 8 in two dimensions and 12 in three. The paper states the probability of forming such a polygon as $8/(4l_p)^4$ in two dimensions and $12/(4l_p)^6$ in three. A sympathetic reader would care because if the relation is right, the bending cost of loop formation in DNA or protein folding is fixed by a few macroscopic parameters rather than by detailed enumeration of all chain shapes.

What carries the argument

The central object is the polymer polygon: a closed lattice walk on a square or cubic lattice whose four arms each have $l_p$ monomers, giving perimeter $4l_p$ and area $l_p^2$. The persistence length $l_p$ is defined as the average chain length per bend and is extracted from the logarithmic derivatives of the partition function with respect to step fugacity and bending weight. The formula $\epsilon_b = k_B T \ln[2(d-1)g l_p]$ is the result of combining the polygon counts (8 in 2D, 12 in 3D) with this persistence-length definition.

What would settle it

Enumerate all closed walks of length $4l_p$ that start and end at the grafted site on a square lattice, without imposing the four-straight-arm restriction; for $l_p=1$ the count should come out to 8 in two dimensions and 12 on the cubic lattice, and the directly computed bending energy should follow $k_B T \ln[2(d-1)g l_p]$ — if either check fails, the central claim is wrong.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that closed-loop ('polygon') conformations of an ideal lattice polymer determine the persistence length, and the bending energy per closed polygon is $$\epsilon_b = k_B T \ln[2(d-1)\, g\, l_p],$$ with $d=2$ on the square lattice and $d=3$ on the cubic lattice, $g$ the step fugacity, and $l_p$ the persistence length. The paper treats the chain as grafted at a point, keeps only walks whose first and last monomers lie on that site, and counts the closed polygons of perimeter $4l_p$ and area $l_p^2$: 8 in two dimensions and 12 in three dimensions. The probability of closing such a polygon is stated as $8/(4l_p)^4$ in 2D and $12/(4l_p)^6$ in 3D, and the logarithmic bending-energy formula is obtained by combining that counting with the persistence length extracted from the random-walk partition function.

Load-bearing premise

The load-bearing assumption is that the persistence length is controlled only by closed polygons of perimeter exactly $4l_p$, with exactly 8 such polygons in two dimensions and 12 in three; if closed walks can form other shapes, or if the side length of the polygon is not $l_p$, the enumeration and the logarithmic formula collapse.

Editorial extensions

If this is right

In two dimensions the bending energy is $k_B T \ln(2 g l_p)$; in three dimensions it is $k_B T \ln(4 g l_p)$; the lattice dimension changes only the constant $2(d-1)$.
The stated closure probabilities, $8/(4l_p)^4$ in 2D and $12/(4l_p)^6$ in 3D, make closed-loop conformations exponentially rarer as $l_p$ grows.
Given $d$ and $g$, a measured persistence length $l_p$ fixes the minimum bending energy of an ideal loop, and conversely the formula can be inverted to extract $l_p$ from a loop-closure energy.
Only 8 (2D) or 12 (3D) four-arm polygons contribute, so the bending energy does not depend on the detailed internal shape of the loop beyond its perimeter and area.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

On other lattice geometries (for example hexagonal or diamond lattices), the same four-arm counting would likely replace $2(d-1)$ by a coordination-number-dependent prefactor; the author does not state this extension.
Because the model treats the chain as ideal, excluded-volume interactions could change the polygon count; a direct test would be to repeat the enumeration with self-avoiding walks and see whether a logarithmic formula survives with a modified constant.
The formula could be used backwards as an experimental estimator: measure loop-closure probability or bending energy in single-molecule experiments on DNA or protein loops and infer $l_p$; this application is not in the paper.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Bottom line: the central formula, ε_b = k_B T ln[2(d-1) g l_p], is not derived. Section III says 'It has been found from these calculations' but shows no calculation. Section V gives counts of 8 and 12 polygons with probabilities 8/(4l_p)^4 and 12/(4l_p)^6; converting those probabilities to energies through a Boltzmann factor gives expressions with different l_p scaling (slope 4 in 2D, slope 6 in 3D), not the claimed logarithms. Figure 2 plots ε_b/k_B T against Log[2g l_p] or Log[4g l_p], which is just the formula plotted against its own argument, so it confirms nothing.

The paper does have a reasonable motivation: define persistence length as average chain length per bend, then ask what bending energy would be required to close the chain into a square polygon. That is a legitimate heuristic. And the target relation, if properly derived, would be a compact connection between dimensionality, step fugacity, and persistence length, which someone might want to cite. But the execution doesn't get there.

The soft spots are structural, not cosmetic. Eq. (2) defines l_p as an ensemble average, but that average is never evaluated; instead the polygon side length is set equal to l_p by fiat, and the perimeter 4l_p is imposed. The counts 8 and 12 are asserted without enumeration. The paper never shows how the probability of a polygon maps to a bending energy. The statements in Sections III and V are mutually inconsistent. There is also a missing Section IV, and the references are heavily self-cited, but that is a side issue; the missing derivation is the problem.

Who gets value from this? Someone looking for a heuristic formula might get a starting point, but it is not a trustworthy result. It should not be sent to a referee as it stands. If the author can actually derive Eq. (3)-(4) from the polygon enumeration, the result would be worth a short note; until then, a desk reject is appropriate.

Referee Report

4 major / 5 minor

Summary. The manuscript considers a random-walk (ideal chain) model of a semi-flexible polymer on square and cubic lattices and defines the persistence length as the average chain length per unit bend. The central claim, stated in the abstract and Section III, is that the minimum bending energy of closed-loop (polygon) conformations is epsilon_b = k_B T ln[2(d-1) g l_p], where d, g, and l_p denote the dimensionality, step fugacity, and persistence length. The paper introduces a partition function and a ratio definition for l_p, asserts that calculations yield Eqs. (3)-(4), and presents Figure 2 plotting epsilon_b/(k_B T) versus Log[2 g l_p] or Log[4 g l_p]. Section V states that there are 8 polygons in two dimensions and 12 in three dimensions, with occurrence probabilities 8/(4l_p)^4 and 12/(4l_p)^6. No explicit statistical-mechanical derivation connects these ingredients to the claimed logarithmic energy formula.

Significance. If the formula epsilon_b = k_B T ln[2(d-1) g l_p] were rigorously derived, it would provide a compact relation among bending energy, persistence length, step fugacity, and dimensionality for ideal lattice polymer loops, and the idea of relating bending energy to the logarithm of a polygon count is interesting. However, the present manuscript does not supply that derivation or a falsifiable numerical test: the only quantitative validation is a plot of the formula against its own argument. Because the central claim is unsupported, the paper does not currently make a significant contribution.

major comments (4)

[Section III, Eqs. (3)-(4)] The two bending-energy expressions are introduced with the phrase 'It has been found from these calculations,' but no calculation is shown in the preceding text. The reader is not told how the partition function in Eq. (1) or the ratio in Eq. (2) is evaluated, nor how the resulting averages lead to a logarithmic dependence on g and l_p. This is the central derivation of the paper, and its absence leaves the main claim unsupported.
[Section II, Eq. (2)] Equation (2) defines the persistence length as the ratio of two statistical averages, <L>/<N_B>. The manuscript never evaluates these averages; instead, the text immediately identifies l_p with the integer side length n of a square polygon and declares the perimeter to be 4n. This replaces the ensemble definition with an ad hoc geometric identification, and the final formula inherits this identification, so the claimed result is an input to the model rather than an output of a calculation.
[Section V] The polygon counts 8 (2D) and 12 (3D) and the probabilities 8/(4l_p)^4 and 12/(4l_p)^6 are asserted without enumeration. More importantly, converting these probabilities into an energy via a Boltzmann factor, -ln P, gives 4 ln(4l_p) - ln 8 in two dimensions and 6 ln(4l_p) - ln 12 in three dimensions. Neither expression equals ln[2(d-1) g l_p]; the three-dimensional expression has a slope of 6 with respect to ln l_p, not 1. Thus the stated counting, even if correct, does not yield Eqs. (3)-(4), and the missing intermediate step is the core of the claimed result.
[Figure 2] Figure 2 plots epsilon_b/(k_B T) against Log[2 g l_p] (2D) and Log[4 g l_p] (3D). With these axes, any relation of the form y = Log[a l_p] is a straight line with unit slope by construction. The figure therefore cannot validate Eq. (3) or (4); it simply redraws the formula being tested.

minor comments (5)

[Section II, Eq. (1)] Equation (1) is not legible as typeset; the summation limits and the integrand appear garbled. Please provide a clear definition of the partition function, including the fugacity weights for steps and bends.
[Section II] The statement 'The perimeter of the polymer polygon is 4n monomers (where n=l_p)' uses l_p both as a real-valued ensemble average from Eq. (2) and as an integer side length; the relation between these two uses should be clarified.
[Section V] The manuscript jumps from Section III to Section V; Section IV is absent, which breaks the organization described in the introduction.
[Figure 2] The figure caption states g=1, but Eqs. (3)-(4) include g as a free parameter; please state whether the figure is meant to show the dependence on l_p only and how g was set in the plot.
[References] The reference list is sparse and mostly cites the author's own prior work; the methods (e.g., recursion relations) are not attributed to a specific external source, making it difficult to place the contribution in context.

Circularity Check

2 steps flagged · score 8.0 of 10

The central bending-energy formula is not derived: lp is imposed as the polygon side length, the polygon counts do not produce Eqs. (3)-(4), and Fig. 2 plots the formula against itself.

self definitional [Sec. II, text immediately after Eq. (2), and Eqs. (3)-(4)]
"The persistence length (l p) of the polym er chain may be defined as the average length of the polymer chain per unit bend and it may be calculated from the following relations: Lp= < L>/< NB> ... The perimeter of the polymer polygon is 4n monomers (where n=l p) i. e. the n mon omers are along each of the direction so that the polymer polygon may be closed one, where n=1, 2, 3…. , lp monomers."

Eq. (2) defines l_p as an ensemble average <L>/<N_B>, but that average is never evaluated. Instead, the paper declares the perimeter of the polygon to be 4n with n = l_p, so the persistence length is simply renamed as the integer side length of the polygon. Eq. (3) then states that the bending energy is k_B T Log[2 g l_p]. The only dependence of this 'derived' energy on l_p is the dependence put in by the identification n = l_p; no calculation from the partition function or from bend statistics supplies the logarithmic form. The central result is therefore constructed from its own input variable rather than derived from the model.
other [Fig. 2 caption and Eqs. (3)-(4)]
"The variation of (Ԑb/kB*T) versus Log [2*g*lp] is shown for two dimensional case, while the variation of (Ԑb/kB*T) versus Log[4*g*lp] is also shown for the case of three dimensions in this figure (here g=1 unit)."

Equations (3) and (4) assert exactly that epsilon_b/k_B T = Log[2 g l_p] and epsilon_b/k_B T = Log[4 g l_p]. Figure 2 therefore plots the quantity epsilon_b/k_B T against the same expression that Eqs. (3)-(4) declare it to equal. The graph is a straight line y = x by construction, so it cannot test or validate the formula. The figure is presented as the numerical evidence for the central claim, but it merely restates the claim; it provides no independent calculation or simulation data.

full rationale

Eq. (2) defines l_p as the ensemble ratio <L>/<N_B>, but that average is never computed; instead the paper states 'The perimeter of the polymer polygon is 4n monomers (where n = l_p)', so the microscopic polygon size is declared equal to the persistence length. Eq. (3) then 'finds' epsilon_b = k_B T Log[2 g l_p]; no calculation connects the two. The only quantitative ingredients in Sec. V are the counts 8 and 12, which enter as probabilities 8/(4 l_p)^4 and 12/(4 l_p)^6; the paper never shows how these probabilities produce Eqs. (3)-(4). If one interprets the energy as -k_B T ln P, the result would scale as 4 ln l_p or 6 ln l_p, not ln l_p, so the enumeration cannot be the derivation. Finally, Fig. 2 plots epsilon_b/k_B T against Log[2 g l_p] and Log[4 g l_p], which is exactly Eqs. (3)-(4); the plotted line is y = x by construction. Thus the central claim is not an independent prediction: its l_p-dependence is inserted through the n = l_p identification, and its only graphical support is a restatement of the formula. The self-citations in Refs. [6-11] are not load-bearing here, so the circularity arises from the self-referential construction of the result rather than from citation practice.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The central claim rests entirely on a small set of modeling choices: ideal lattice chains, closed loops, a definition of l_p as average length per bend, and hand-counted polygon geometries. No fitted parameters appear, but the result is not independently grounded; it is determined by these assumptions.

assumptions (5)

domain assumption A random walk on a square or cubic lattice with step fugacity g and bend fugacity k represents a semi-flexible polymer chain.
Section II states 'random walk is used on a square and a cubic lattice to mimic the polymer chain conformations'. This neglects excluded volume, monomer interactions, and continuum elasticity.
domain assumption Only closed-loop conformations that begin and end at the graft site are considered, and the chain is treated as ideal.
Section II: 'only those conformations were taken into considerations which have end monomer on the site where the polymer chain is grafted' and 'the chain is modeled as an ideal polymer'.
ad hoc to paper The persistence length l_p is defined as average chain length per unit bend, and the polygon side length is set equal to l_p.
Eq. (2) defines l_p = <L>/<N_B>, and Section II states 'The perimeter of the polymer polygon is 4n monomers (where n=lp)'. The equality of the side length with the persistence length is imposed, not derived.
ad hoc to paper The relevant closed polygons are exactly 8 in two dimensions and 12 in three dimensions.
Section V states these counts with no exhaustive enumeration or derivation.
ad hoc to paper The bending energy is obtained from the logarithm of the fugacity-weighted polygon count or from the ratio defining l_p, without an explicit statistical-mechanical derivation.
The paper says 'It has been found from these calculations' but does not show how Eq. (3) or the final formula follows from the partition function.

how reviews work

0 comments

Cite this review

Pith. "Pith review of The bending energy of a semi-flexible polymer chain and the polygons of the polymer chain." pith.science (2026). https://pith.science/paper/75DMPN62

@misc{pith2026190805586,
  author       = {Pith},
  title        = {Pith review of: The bending energy of a semi-flexible polymer chain and the polygons of the polymer chain},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/75DMPN62}},
  note         = {Machine review of arXiv:1908.05586}
}

read the original abstract

We consider random walk model of a semi-flexible polymer chain on a square and a cubic lattice to enumerate conformations of the polymer chain in two and three dimensions, respectively. The bending energy of the chain is assumed as the key factor which controls the minimum average length of the chain in between two successive bends in the chain; and the average length of the chains per unit bend is defined as the persistence length of the polymer chain. It has been found that the minimum energy required to bend the chain is expressed in the form of simple relation which includes space dimensionality, step fugacity and persistence length.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

12 extracted references · 12 canonical work pages

[2]

in the case of two d imensions and while there are 12 polygons (each polygon has the perimeter 4*l p and the area is l p

work page
[3]

The probability of finding such polymer polygon in two dimensions is 8/(4*lp)4 and in three dimensions the probability is 12/(4*lp)6

which are found to have first as well as the last monomer on the site where the polymer chain is grafted. The probability of finding such polymer polygon in two dimensions is 8/(4*lp)4 and in three dimensions the probability is 12/(4*lp)6 . The value of lp is 1, 2, 3, …… , n monomers. The standard methods of the statistical Physics were used to calculate ...

work page 2019
[4]

N. M. Toan, G. Morrison, C. Hyeon and D. Thirumalai, J. Chem. Phys. B (2008) 112 6094

work page 2008
[5]

Jun and J

S. Jun and J. Bechhoefer, Physics in Canada (2003) 59 85-92

work page 2003
[6]

Tyagi, S

S. Tyagi, S. A. E. Marris and F. R. Kramer, Nature Biotech. (2000) 18 1191; R. Schleif, Annu. Rev. Biochem. Sci. (1992) 26 199; K. Rippe, Trends Biochem. Sci. (2001) 26 733

work page 2000
[7]

Kratky, G

O. Kratky, G. Porod, Recl. Trav. Chim. Pays-Bas 68 (1949) 1160; V. Privman, N.M. Svrakic, Directed Models of Polymers, Interfaces and Clusters: Scaling and Finite-Size Properties, Springer, Berlin, 1989; M. Doi, S.F. Edwards, Theory of Polymer Dynamics, Clarendon, Oxford, 1992

work page 1949
[8]

P. K. Mishra, S. Kumar and Y. Singh, Physica A (2003) 323 453-465; J. Zi Zhang, X. Y. Peng, S. Liu, B. P. Jiang, S. C. Ji and X. C. Shen, Polymers (2019) 11 295; P. K. Mishra, J Phys: Cond Matt. (2010) 22 155103

work page 2003
[9]

P. K. Mishra, Fizika A (2010) 19 145-152

work page 2010

Show all 12 references

[10]

P. K. Mishra, Phase Transitisions (2011) 84 291-298

2011
[11]

P. K. Mishra, Condens Matter Phys. (2014) 17 23001

2014
[12]

P. K. Mishra, Phase Transitisions (2015) 88 593-604

2015
[13]

P. K. Mishra, Ind J Phys (2017) 91 1297-1304

2017

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [2]

in the case of two d imensions and while there are 12 polygons (each polygon has the perimeter 4*l p and the area is l p

work page

[2] [3]

The probability of finding such polymer polygon in two dimensions is 8/(4*lp)4 and in three dimensions the probability is 12/(4*lp)6

which are found to have first as well as the last monomer on the site where the polymer chain is grafted. The probability of finding such polymer polygon in two dimensions is 8/(4*lp)4 and in three dimensions the probability is 12/(4*lp)6 . The value of lp is 1, 2, 3, …… , n monomers. The standard methods of the statistical Physics were used to calculate ...

work page 2019

[3] [4]

N. M. Toan, G. Morrison, C. Hyeon and D. Thirumalai, J. Chem. Phys. B (2008) 112 6094

work page 2008

[4] [5]

Jun and J

S. Jun and J. Bechhoefer, Physics in Canada (2003) 59 85-92

work page 2003

[5] [6]

Tyagi, S

S. Tyagi, S. A. E. Marris and F. R. Kramer, Nature Biotech. (2000) 18 1191; R. Schleif, Annu. Rev. Biochem. Sci. (1992) 26 199; K. Rippe, Trends Biochem. Sci. (2001) 26 733

work page 2000

[6] [7]

Kratky, G

O. Kratky, G. Porod, Recl. Trav. Chim. Pays-Bas 68 (1949) 1160; V. Privman, N.M. Svrakic, Directed Models of Polymers, Interfaces and Clusters: Scaling and Finite-Size Properties, Springer, Berlin, 1989; M. Doi, S.F. Edwards, Theory of Polymer Dynamics, Clarendon, Oxford, 1992

work page 1949

[7] [8]

P. K. Mishra, S. Kumar and Y. Singh, Physica A (2003) 323 453-465; J. Zi Zhang, X. Y. Peng, S. Liu, B. P. Jiang, S. C. Ji and X. C. Shen, Polymers (2019) 11 295; P. K. Mishra, J Phys: Cond Matt. (2010) 22 155103

work page 2003

[8] [9]

P. K. Mishra, Fizika A (2010) 19 145-152

work page 2010