Molecular simulation of melting in tetracosane (C24H50) monolayers and bilayers on graphite

Cary Pint

REVIEW

Molecular simulation of melting in tetracosane (C24H50) monolayers and bilayers on graphite

Cary Pint

Soft Condensed Matter · Submitted 21 Feb 2006 · 1 reference

Pith's verdict: Layer-by-layer melting signal is consistent and worth publishing, but the headline monolayer-vs-bilayer Tm claim is confounded by commensurability and needs a control simulation.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

T0 AI review An AI read the paper and answered a fixed list of questions.
T1 Audit Automatic checks that the references and files the paper relies on exist and match.
T2 Compute The paper's code or calculations re-run with the same results.
T3 Formal statement The claims written out in a formal math language a computer can read.
T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →

Challenge this review Comment (0) Read on arXiv ↗

Cite ▾

arXiv IDcond-mat/0602478 (version 1)

Pith IDGQS7HSZU, Pith's permanent record of this paper and its review

CitationCary Pint. "Molecular simulation of melting in tetracosane (C24H50) monolayers and bilayers on graphite." arXiv:cond-mat/0602478 (2006). https://pith.science/paper/GQS7HSZU

Cite this reviewBibTeX for Pith's review of this paper

Claim paper

Why Pith says this

Highlight any text to comment on it.

Editor's note

read the letter

Worth a look if you care about alkane/graphite melting. The genuinely new contribution is the bilayer simulation with separate top/bottom layer analysis, plus four monolayer variations that test mechanisms. The layer-by-layer ordering difference is a solid observation: order parameter, methyl-methyl distance, bond-roll parameter, and out-of-plane fluctuations all point in the same direction, with the top layer losing order ~20 K before the bottom layer. That consistency is the best evidence in the paper, and I buy it as a real effect within the model.\n\nThe problem is the abstract's claim that the bilayer melts significantly hotter than the monolayer. That comparison is not controlled. The monolayer is uniaxially incommensurate (a=4.53, b=4.66); the bilayer is fully commensurate (a=4.26, b=5.36). In Section 3.4 the author runs a fully commensurate monolayer using the bilayer's lattice parameters and finds it melts at ~380 K, about the same as the bilayer's bottom layer (~375 K). So the 'bilayer increase' could simply be commensurability, not coverage. The paper even says the phase diagram is delicate, but the headline doesn't carry that caveat. The bottom layer may still get extra stability from the top layer, but the data don't isolate it. A lattice-matched monolayer control would settle it.\n\nOther soft spots are the usual: 1.2-1.4 ns runs, no error bars, no finite-size scaling, small cells. The 20 K layer separation could in principle be a sampling artifact, though the multi-observable consistency makes that less likely.\n\nCredit where due: the variations are the most useful part. Stiffening the dihedrals moves melting to ~475 K, isolating molecules lets them kink easily, truncating interactions drops Tm to ~310 K. Those cleanly separate the roles of intramolecular stiffness, packing, and intermolecular interactions, and they make sense.\n\nBottom line: this is a credible simulation study with one confounded headline claim. It deserves a serious referee; I'd send it back for major revision, asking for the matched-lattice monolayer/bilayer comparison and at least some convergence or block-averaging analysis.

Referee report

Significance

The topic is relevant to surface freezing, quasi-2D melting, and the behavior of physisorbed alkane films, and the layer-resolved analysis is a useful approach. The paper ships a substantial set of order parameters (OP2, methyl-methyl distance, bond-roll parameter, out-of-plane fluctuations, density profiles) and multiple control variations, which are valuable diagnostics. If the central claims held, the work would provide a clear molecular picture of how an additional molecular layer alters melting relative to a monolayer, with implications for understanding the crossover from 2D to bulk behavior. However, the principal quantitative claim about the bilayer melting temperature is undermined by a confounded comparison, and the absence of statistical error bars and finite-size checks leaves the 20 K layer-separation estimate unsupported.

Major objections (3)

[§3.1 and Table I] The central claim that the bilayer melts significantly hotter than the monolayer is not a controlled comparison. Table I gives different lattice parameters for the monolayer (a=4.53 Å, b=4.66 Å, uniaxially incommensurate) and the bilayer (a=4.26 Å, b=5.36 Å, fully commensurate). Section 3.4 reports that a fully commensurate monolayer, using the same lattice parameters as the bilayer, melts at ca. 380 K, which is comparable to the bilayer bottom layer melting temperature of ca. 375 K reported in Section 3.1. Thus the melting-temperature increase attributed to the bilayer can be reproduced by a monolayer alone when the lattice is changed; the comparison conflates coverage with commensurability. A control simulation with matching lattice parameters (for example, a bilayer on a uniaxially incommensurate lattice, or a direct baseline comparison against the fully commensurate monolayer) is needed to isolate the coverage effect, and the claims in the abstract and Section 5 should be reframed accordingly.
[§2 and §3.1] No statistical error bars, time-convergence tests, or finite-size scaling are reported for any of the order parameters from which melting temperatures are read. The monolayer and bilayer cells contain only 46 and 96 molecules, respectively, and trajectories are 1.2-1.4 ns total including equilibration. The 20 K separation between the loss of strict rectangular-centered order in the top layer (above ca. 355 K) and the bottom layer (above ca. 375 K) could be within the sampling noise of such short, small runs, particularly since the transition in the monolayer is described as gradual. The authors should provide block averages or other error estimates for OP2, the methyl-methyl distance, and the bond-roll parameter near the transition, and should demonstrate that the reported order-parameter values have converged over the averaging window.
[§3.4 and §4.1] The fully commensurate monolayer variation in Section 3.4, which melts at ca. 380 K, is directly relevant to the paper's headline comparison, yet it is not integrated into the discussion in Section 4.1 of why the bilayer bottom layer melts at ca. 375 K. This variation shows that the lattice commensurability, rather than the presence of a second layer, can account for the elevated melting temperature. The discussion should explicitly acknowledge that the monolayer-versus-bilayer comparison is confounded and should either present a control that isolates the coverage dependence or substantially soften the coverage-related conclusion.

Minor objections (5)

[§3.1] The layer-separation criterion is stated as "the molecular centers of the molecules in the top layer be greater 1.56 Å less than the average height of all molecules in the top layer," which is internally inconsistent as written and appears to contain a typographical error; the intended criterion should be clarified.
[§3.2] The equation for the bond-roll order parameter βroll is garbled in the manuscript; the summation limits and the vector dot/cross product structure should be written out clearly so the definition is unambiguous.
[General] The manuscript contains numerous typographical and formatting errors (e.g., "seperately," "analagous," "perpindicular," "psuedoatoms," "equilbration," stray "6" characters throughout the text, and garbled figure labels). A thorough copyedit is needed before publication.
[§4.2] The comparison to experiment for the bilayer relies on an analogy to hexane on graphite, as the authors state that no direct experimental bilayer melting temperature is available; this limitation should be stated prominently in the abstract or conclusions rather than appearing only in the discussion.
[References] References [8] and [36] both cite a Ph.D. dissertation by B. Matthies but with different years (1999 and 1989); one of these is likely a typo and should be checked.

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

Bilayer tetracosane melts at a higher temperature than the monolayer, and the top layer disorders first.

The reading

This paper uses molecular dynamics to ask how melting changes when a physisorbed film of tetracosane (C24H50) on graphite grows from one layer to two. It reports that the bilayer melts at a significantly higher temperature than the monolayer, and that the monolayer itself stays ordered above the bulk melting point of tetracosane. The simulations also show layer-by-layer melting: the top layer forms chain-kinking gauche defects and rolls its chains at lower temperatures, while the bottom layer keeps its rectangular-centered order until roughly 20 K higher. If this picture is right, a second layer stabilizes the interfacial layer, and coverage, not just chain length, controls how alkane films melt.

What carries the argument

The argument is carried by layer-resolved observables computed from the same trajectories. The structural order parameter $OP_2 = \frac{1}{N_m}\sum_i \cos(2\phi_i)$ tracks the in-plane azimuthal alignment of molecular long axes; the average methyl–methyl distance reports trans versus gauche chain conformations; the bond-roll parameter $\beta_{\mathrm{roll}}$ measures whether chains lie flat on the substrate or have rolled perpendicular to it; and center-of-molecule density profiles along the long axis show lamellae broadening and merging. Averaging these quantities separately for bottom and top layers is what lets the paper attribute different melting behavior to each layer and identify the coupling between them.

Extended reading notes

Core claim

On its own terms, the paper claims that in constant-NVT molecular dynamics the C24H50 monolayer loses rectangular-centered order starting near 340 K and is essentially disordered by 360 K, while the bilayer's bottom layer keeps that order until about 375 K and the top layer begins losing it near 355 K. Gauche-defect formation and perpendicular chain rolling appear in the top layer at lower temperatures than in the bottom layer, yet the top layer retains significant order until the bottom layer gives way, so the two layers are coupled. The author interprets this as layer-by-layer melting in which the top layer wants to disorder first but cannot fully do so while the interfacial layer remains solid, and attributes the bottom layer's extra stability to confinement by the second layer. Supporting variations show that suppressing gauche defects pushes monolayer melting up to about 475 K, that a fully commensurate monolayer melts near 380 K (above the incommensurate monolayer), and that truncating intermolecular interactions to first neighbors lowers melting to about 310 K. The paper further claims these results are consistent with neutron diffraction and STM observations on similar alkane/graphite systems.

Load-bearing premise

The load-bearing premise is that 1.2–1.4 ns trajectories on cells of 46–96 molecules sample the slow collective motions that control melting; if those runs are too short or too small, the reported roughly 20 K difference between layers and the melting temperatures themselves could be artifacts.

Editorial extensions

If this is right

If the bilayer truly melts above the monolayer, then adding a second layer raises the temperature at which an alkane boundary film loses solid order, so film stability on graphite is not a monotonic function of reducing confinement.
The layer-by-layer picture predicts that a partially disordered top layer still cannot become fluid until the interfacial layer melts, so any experiment or application that pins the bottom layer should also suppress total film melting.
The variation with suppressed gauche defects, which delays melting to about 475 K, implies that chain kinking is the rate-limiting step in monolayer melting, while lamellar fluctuations alone are not enough.
The commensurability variation, where the fully commensurate monolayer melts near 380 K, implies that lattice matching with the substrate can raise the melting point by about 40 K relative to the uniaxially incommensurate layer.
The shorter-interaction variation, which lowers melting to about 310 K, implies that longer-range intermolecular packing is what sustains order above the bulk melting temperature.

Reading between the lines

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

A natural extension, not stated in the paper, is that with three or more layers the melting temperature should decrease layer by layer away from the graphite surface, eventually approaching bulk behavior; the bottom layer would remain the last to disorder.
A direct numerical extension would be to run the same bilayer at several independent initial configurations and larger cells to see whether the roughly 20 K separation between top- and bottom-layer order loss is a thermodynamic feature or a short-trajectory artifact.
The alternating parallel/perpendicular chain orientation hinted by the bond-roll parameter suggests a concrete experimental signature: polarized infrared or neutron reflectivity measurements on a well-characterized two-layer film should detect the top-layer chains standing up before the bottom layer disorders.
For applications, the confinement effect implies that a strongly adsorbing substrate plus a capping layer—not just a single coating—should be the design route for alkane boundary films that must keep their solid order at elevated temperatures.

Original abstract

This work reports an investigation into the solid phase behavior and melting behavior of tetracosane (C24H50) monolayers and bilayers physisorbed onto the graphite basal plane using molecular dynamics simulations performed in the constant (N,V,T) ensemble. This study is conducted in order to perform an in-depth analysis of the melting process in the quasi-2D monolayer, and compare the molecular behavior in the monolayer and the bilayer at temperatures near and at the melting transition. The bilayer is studied by looking at quantities for the bottom and top layer separately in order to break apart the evolution of differing structural and molecular properties that are normally associated with melting in the monolayer. Simulations suggest that the melting temperature in the bilayer is increased significantly from that of the monolayer- which already occurs at temperatures higher than the bulk tetracosane melting temperature. In addition, a layer-by-layer melting effect is observed where significant factors associated with the melting transition- the formation of molecular gauche defects and the perpendicular rolling of the molecule chains- occur in the top layer of the bilayer at significantly lower temperatures than they occur in the bottom layer. Variations are conducted to better understand the role of different interactions in melting of the crystalline monolayer, and the results of simulations are discussed in comparison to recent STM and diffraction experiment and found to be in good agreement.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 5,240 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.08 · human participants 0.01 · controlled comparison 0.52 · randomized assignment 0.01 · states quantitative result 0.08 · compares with prior results 0.95 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; show every role.

This work reports an investigation into the solid phase behavior and melting behavior of tetracosane (C24H50) monolayers and bilayers physisorbed onto the graphite basal plane using molecular dynamics simulations performed in the constant (N,V,T) ensemble. This study is conducted in order to perform an in-depth analysis of the melting process in the quasi-2D monolayer, and compare the molecular behavior in the monolayer and the bilayer at temperatures near and at the melting transition. The bilayer is studied by looking at quantities for the bottom and top layer separately in order to break apart the evolution of differing structural and molecular properties that are normally associated with melting in the monolayer. Simulations suggest that the melting temperature in the bilayer is increased significantly from that of the monolayer- which already occurs at temperatures higher than the bulk tetracosane melting temperature. In addition, a layer-by-layer melting effect is observed where significant factors associated with the melting transition- the formation of molecular gauche defects and the perpendicular rolling of the molecule chains- occur in the top layer of the bilayer at significantly lower temperatures than they occur in the bottom layer. Variations are conducted to better understand the role of different interactions in melting of the crystalline monolayer, and the results of simulations are discussed in comparison to recent STM and diffraction experiment and found to be in good agreement.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#	sentence	role	the authors' own work	states the main conclusion	strength
s0	This work reports an investigation into the solid phase behavior and melting behavior of tetracosane (C24H50) monolayers and bilayers physisorbed onto the graphite basal plane using molecular dynamics simulations performed in the constant (N,V,T) ensemble.	contribution 0.92	0.95	0.04	asserted 0.97
s1	This study is conducted in order to perform an in-depth analysis of the melting process in the quasi-2D monolayer, and compare the molecular behavior in the monolayer and the bilayer at temperatures near and at the melting transition.	contribution 0.96	0.68	0.04	asserted 0.71
s2	The bilayer is studied by looking at quantities for the bottom and top layer separately in order to break apart the evolution of differing structural and molecular properties that are normally associated with melting in the monolayer.	method 0.90	0.84	0.04	asserted 0.94
s3	Simulations suggest that the melting temperature in the bilayer is increased significantly from that of the monolayer- which already occurs at temperatures higher than the bulk tetracosane melting temperature.	result 1.00	0.92	0.57	hedged possibility 0.99
s4	In addition, a layer-by-layer melting effect is observed where significant factors associated with the melting transition- the formation of molecular gauche defects and the perpendicular rolling of the molecule chains- occur in the top layer of the bilayer at significantly lower temperatures than they occur in the bottom layer.	result 1.00	0.93	0.25	asserted 0.98
s5	Variations are conducted to better understand the role of different interactions in melting of the crystalline monolayer, and the results of simulations are discussed in comparison to recent STM and diffraction experiment and found to be in good agreement.	result 0.34	0.89	0.05	asserted 0.98

What kind of paper the classifier reports this is: simulation study 1.00 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read

arXiv listing abstract v1

not examined

full text
figures and tables
references and citing sentences
relations to other papers

levels run

P abstract

not run

S sections
T body sentences
C citations
R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures

**Figure 3.** Snapshots of (a) the low-temperature crystalline monolayer C24 solid phase in simulations (260 K), (b) the C24 monolayer at temperatures approaching melting (310 K), and (c) the top and bottom layers of the C24 bilayer at 320 K

**Figure 4.** Height profiles for the (a) C24 monolayer and (b) C24 bilayer at temperature points near and after the melting transition. 200 220 240 260 280 300 320 340 360 380 400 0 5 10 15 20 25 30 35 <rCH3-CH > /Å 3 Temperature/K (a) 200 220 240 260 280 300 320 340 360 380 400 420 0 5 10 15 20 25 30 35 (b) bottom layer top layer <rCH3-CH > /Å 3 Temperature/K

**Figure 6.** Bond-roll order parameter, βroll, for the C24 monolayer and bilayer at simulated temperature points. 200 240 280 320 360 400 440 0 1 2 3 4 C24 monolayer C24 bilayer <∆z>/Å 2 T/K

**Figure 7.** Total mean square out-of-plane fluctuations at temperatures near and after the melting transition for the C24 monolayer and bilayer. Fig. 8a. Center-of-molecule density distributions for the C24 monolayer in the direction of the molecule long axes,defined at labelled temperature points

**Figure 9.** Snapshots of the C24 monolayer at 300 K with (left panel) no gauche defects present in the layer, and (right panel) with only nearest-neighbor pair interactions allowed

Showing 5 of the paper's 9 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:cond-mat/0602478 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.

[1] [1]

1M.J.W.6Povey,6S.A.6Hindle,6A.6Aarflot,6and6H.6Hoiland,6 Cryst.JGrowthJDes.J6,62976(2006).6 2S.M.6Clarke,6A.6Inaba,6T.6Arnold,6and6R.K.6Thomas,6J.JTherm.JAnal.JCal. 657 ,66436(1999).6 3M.A.6 Castro,6 S.M.6 Clarke,6 A.6 Inaba,6 T.6 Arnold,6 and6R.K.6 Thomas,6 J.J Phys.J Chem.J SolidsJ60 ,6 14956 (1999).6 4M.A.6Castro,6S.M.6Clarke,6A.6Inaba,6T.6Arnold,6and6...

View on Pith 2006