{"id":"b6c31e70-156e-4a1e-9e36-d954fa9aeed8","arxiv_id":"1908.04666","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"The paper derives an entropy-based shear stress law for twisted bilayer graphene, tau = 3 k T cot(theta/2) / (2 pi R^3), and argues configurational entropy is the origin of structural superlubricity.","lead":"This paper proposes that the vanishing friction in twisted bilayer graphene, structural superlubricity, comes from configurational entropy rather than force cancellation. It derives a shear stress formula depending on twist angle, temperature, and flake size, and claims it matches several experiments.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (2) asserts without derivation that the microstate ratio equals the moiré-cell area ratio; it is the sole source of the entropy behind Eq. (6), and no statistical-mechanical counting supports it.","rationale":"The reader's weakest assumption identifies Eq. (2), and my reading agrees: it is the single most load-bearing step. All subsequent results, especially Eq. (6), depend on converting the moiré-cell area ratio into a Boltzmann entropy difference. The paper provides no derivation, no counting argument, and no independent check that the number of microstates scales as the cell area. The internal-energy MD calculations only show dU/dθ ≈ 0, which is consistent with the standard incommensurability argument but does not establish the entropy term. In fact, if one takes the entropy as extensive and multiplies by the number of moiré cells, the torque acquires a different size dependence, so the paper's own thermodynamic framework is not internally justified. There is also a secondary inconsistency: Eq. (1) as written gives am decreasing monotonically on (0°, 60°), contradicting the stated symmetry about 30° that is used to motivate the negative-friction branch. This reinforces, rather than replaces, the concern about Eq. (2). The proposed free-energy calculation is a direct test: it measures the actual entropic torque in the same model and would settle whether the microstate-area identification is physically correct. Since the central claim lacks this support, the reader's REJECT verdict should stand.","tokens_in":6652,"tokens_out":10390,"duration_ms":125255,"concrete_test":"Use classical MD with the same AIREBO potential to compute the Helmholtz free energy F(θ) of a rigid circular graphene flake on an infinite substrate as a function of twist angle, sampling the flake's lateral displacement over the moiré cell via umbrella sampling or free-energy perturbation. Then compute the torque from M(θ) = −∂F/∂θ and compare with Eq. (5), M = kT cot(θ/2), for R = 1, 2, 5, and 10 nm at T = 300 K. If M deviates from Eq. (5) or depends on R, the Ω ∝ area assumption fails. A complementary check is to count the distinct translational registries sampled and test directly whether Ωθ/Ω0 equals Sθ/S0.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (2) is the only quantitative input that converts a geometric moiré-area ratio into an entropy difference, and it is asserted rather than derived: Ωθ/Ω0 = Sθ/S0. No partition function or microstate count is given, and for rigid layers at fixed θ the atomic configuration is deterministic; the number of relative registries is not obviously proportional to moiré-cell area. This ΔS is then substituted into dQ = T dS for the whole flake, so the derivation also uses a per-cell entropy as if it were the total entropy, omitting the number of moiré cells. That omission is what makes Eq. (5) size-independent. If the entropy were extensive, e.g. N_cells k ln(Sθ/S0), the predicted torque would have a different R and θ dependence. The authors' own MD result that U is nearly independent of θ removes the energy term, leaving Eq. (2) as the sole support for the cotangent law; if the area-to-microstate identification fails, Eq. (6) and the negative-friction prediction have no basis.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a thermodynamic origin for structural superlubricity in twisted bilayer graphene. The authors define a configurational entropy change on the basis of Eq. (2), where the ratio of microstates is set equal to the ratio of moiré unit-cell areas, leading to ΔS = 2k ln(a_m/a). Combining this with the first law under the assumption dU = 0, they derive a torque M = kT cot(θ/2) and a shear stress τ = 3kT cot(θ/2)/(2πR^3) in Eqs. (5) and (6). The paper claims this formula explains experimental observations of the twist-angle dependence of friction, size-dependent shear stress, temperature-dependent sliding, spontaneous reorientation, and persistence under high normal loads, and it reports MD simulations indicating that the internal energy is nearly independent of twist angle.","tokens_in":6930,"tokens_out":4142,"duration_ms":45109,"significance":"If the central relation Eq. (2) were derivable from statistical mechanics, the model would offer a remarkably concise analytic formula for interlayer shear stress in twisted bilayer contacts, and the predicted τ ~ 1/R^3 scaling and temperature dependence would be useful design rules for superlubric systems. The paper is clearly written and identifies a potentially important, nonstandard contribution (configurational entropy) to the energetics of twisted interfaces. The MD data in Fig. 1(d) and Fig. 4, showing near-constancy of the potential energy with twist angle for rigid layers, are useful supporting evidence for the dU ≈ 0 assumption. However, the load-bearing step, Eq. (2), is asserted without a microscopic counting argument, and the treatment of entropy as a per-unit-cell quantity without accounting for the number of moiré cells undermines the physical meaning of the derived torque. As a result, the central quantitative prediction is not established, and the experimental comparisons, though suggestive, do not independently validate the entropy mechanism.","major_comments":[{"comment":"Equation (2) is the central assumption of the paper, yet it is asserted without derivation: the ratio of microstate numbers is set equal to the ratio of moiré unit-cell areas, Ωθ/Ω0 = Sθ/S0. For rigid crystalline layers at a fixed twist angle, the atomic registry is deterministic, and it is not obvious why the number of microstates should scale with the moiré cell area. No partition function, phase-space argument, or counting procedure is provided. Because Eq. (2) is the sole source of the cotangent law in Eq. (5) and of the entropy term in Eq. (6), the central quantitative claim is unsupported unless this identification can be justified from a statistical-mechanical model.","section":"Eq. (2)"},{"comment":"The entropy ΔS in Eq. (2) is defined for a single moiré unit-cell, but Eq. (5) applies this entropy to the entire graphene flake without multiplying by the number of moiré cells. For a flake of radius R, the number of cells scales as N_cells ~ (R/a_m)^2, and the total entropy should be roughly N_cells times the per-cell entropy. If this extensivity is restored, the torque M would acquire an R-dependence and a different θ-dependence through the θ-dependence of N_cells; the size-independent torque in Eq. (5) is an artifact of treating the per-cell entropy as the system entropy. The statement in the text that the entropy is determined by one unit-cell 'as long as R > am' is not justified and is a critical flaw.","section":"Eq. (2) to Eq. (5): extensivity"},{"comment":"The derivation assumes dU = 0 during quasi-static rotation, supported only by MD simulations with the AIREBO potential for rigid layers. The authors themselves acknowledge in the paragraph following Eq. (5) that the model breaks down when energy contributions become comparable to entropy contributions, which is precisely the regime of standard registry-dependent interlayer interactions. Since real interfaces have finite stiffness and may undergo relaxation, the rigidity assumption is restrictive, and the extension of the dU ≈ 0 conclusion to normal loads up to 10 GPa in Fig. 4 is based on the same rigid-layer approximation. The validity of Eq. (6) beyond the idealized rigid case is therefore not established.","section":"Eq. (4) and dU = 0"},{"comment":"The prediction of a sign reversal of τ for 30° < θ < 60°, i.e., spontaneous rotation toward commensurate registry, is striking and falsifiable, but it rests entirely on the derivative of Eq. (3). The cited experimental support (Refs. [9] and [21]) is qualitative; no quantitative comparison with the predicted magnitude or angle dependence is attempted. In Fig. 2, the agreement with the experimental data of Ref. [8] is only approximate, and the statement that the theoretical curve is 'a little larger' is not a quantitative validation, especially because the plotted range and the experimental conditions (static vs. sliding friction) are not matched in detail.","section":"Negative friction prediction and experimental comparison"},{"comment":"The paper claims that Eq. (6) 'universally explains' a wide range of experiments, but most comparisons are qualitative. For the DWCNT pull-out (Ref. [13]), the twist angle is unknown, so the apparent agreement at R = 3663 nm is a consistency check with an adjustable geometric parameter, not an independent test of the entropy model. Given that Eq. (6) has no free parameters, a more stringent test would require experiments with known twist angles, sizes, and temperatures; the present set of comparisons does not achieve this.","section":"Comparisons with experiments"}],"minor_comments":[{"comment":"The abstract says the configuration entropy is 'directly derived from the Helmholtz free energy,' but the full text does not derive Eq. (2) from a free energy; it postulates the microstate-area proportionality. Please adjust the wording to avoid overstatement.","section":"Abstract and introduction"},{"comment":"The y-axis of Fig. 1(c) is missing units, and Fig. 1(d) uses 'internal energy' where the text describes potential energy; please clarify the notation.","section":"Fig. 1(c) and 1(d)"},{"comment":"The phrase 'micro-status' is nonstandard; the paper should use 'microstates' consistently.","section":"Terminology"},{"comment":"The text mentions a '9-mm-long inner shell,' but the original reference is about a micron-scale tube; this appears to be a typo and should be corrected.","section":"Reference [13]"},{"comment":"The Methods section is absent; the MD simulations lack details on system size, boundary conditions, equilibration, and whether the bottom layer is truly infinite. These details are necessary to assess the dU ≈ 0 conclusion.","section":"MD simulation details"}],"recommendation":"reject","confidential_remarks":"The central issue is not merely a presentation problem but a fundamental gap in the derivation: Eq. (2) is an unjustified postulate, and the non-extensive use of the configurational entropy leads to a torque that is not physically consistent. The paper's main contribution, a simple closed-form friction law, therefore rests on a broken foundation. Even if the authors were to replace Eq. (2) with a proper statistical-mechanical counting, the resulting formula would likely change, so the required revision goes beyond the scope of a normal major revision. I recommend rejection, though I note that the authors have identified an interesting and underexplored physical question that could be pursued in a future, more rigorous study."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nYou should know two things about this paper. The first is that it contains a clean, testable closed-form prediction that is genuinely new: tau = 3 k T cot(theta/2) / (2 pi R^3). The second is that the load-bearing step behind that formula, Eq. (2), is a postulate, not a derivation. Everything after it is algebra.\n\nWhat the paper does well: it reframes structural superlubricity in thermodynamic language, which is a fresh angle against the standard force-cancellation picture. The authors check their main assumption that internal energy is independent of twist angle with MD, and the result (dU ≈ 0) is credible. They also honestly flag limitations: perfect rigid lattices, no defects, and breakdown when chemical bonds form. The formula's symmetry with 60° and the 1/R^3 decay are elegant, and the comparison to the DWCNT experiment is suggestive.\n\nThe soft spot is where it hurts. In Eq. (2), the ratio of microstates is set equal to the ratio of moiré-cell areas without derivation. It is not obvious why the number of atomic configurations scales with the area of the superlattice cell; for rigid layers at fixed angle, the relative registries are constrained, and a statistical-mechanical counting is absent. More seriously, the paper uses the entropy per moiré cell as if it were the total entropy. That is why Eq. (5) is size-independent. If the entropy were extensive, with a factor of the number of cells, the torque would pick up an extra factor of (R/a_m)^2 and the 1/R^3 behavior would become 1/R or something else. The negative friction for 30°–60° is an artifact of that choice; it only appears because the area ratio is inverted. The experimental fits are qualitative, with no error bars and a static-vs-sliding mismatch that could hide a lot. The normal-load MD is a side point and doesn't rescue the derivation.\n\nI read the paper as a plausible hypothesis dressed as a derivation. The formula is new and worth discussing, but the current manuscript does not establish the entropy claim. If someone can replace Eq. (2) with a real partition-function argument, the paper could be valuable. As is, I would not cite it for the result, but I would bring it to a reading group to debate the extensivity issue.\n\nRecommendation: this deserves a serious referee, not a desk reject. The question is important, and the flaw is subtle enough that a good referee might help the authors fix it.","headline":"The new friction law rests entirely on an unproven entropy-area postulate, so the formula is not yet established, but the question is worth debating.","tokens_in":7402,"tokens_out":2960,"would_cite":false,"duration_ms":27684,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82C05","80A10"],"pacs":["81.40.Pq"],"model":"deepseek-v4-flash","headline":"A thermodynamic derivation traces structural superlubricity to configurational entropy, giving interlayer shear stress $\\tau = 3kT\\,\\cot(\\theta/2)/(2\\pi R^3)$ and a self-rotation branch beyond 30 degrees.","keywords":["structural superlubricity","twisted bilayer graphene","configurational entropy","moiré unit cell","interlayer shear stress","thermodynamic model","negative friction","size scaling"],"falsifier":"Measure, for a rigid, defect-free circular graphene flake of radius R twisted against graphite at fixed temperature, the torque as a function of twist angle from 0 to 60 degrees. The paper's claim predicts $M = kT\\,\\cot(\\theta/2)$: divergence near zero angle, linear temperature dependence, no dependence on R beyond the entropy area factor, and a sign reversal at 30 degrees. Observing no sign reversal, or a torque that is not linear in temperature, would falsify Eq. (6); a direct enumeration of finite-flake microstates that fails to reproduce $k\\ln(S_\\theta/S_0)$ would falsify Eq. (2).","tokens_in":6488,"feed_emoji":"🌀","tokens_out":9645,"duration_ms":91390,"temperature":0.7,"pith_summary":"This paper tries to establish a thermodynamic origin for structural superlubricity, the near-zero friction between atomically flat layers whose lattices are twisted out of commensurate registry. Its claim is that the effect comes from configurational entropy: a twisted bilayer has more accessible atomic arrangements inside its moiré unit cell, and that entropy difference dominates the free energy during quasi-static rotation. Working from this premise with rigid layers and constant internal energy, the paper derives a closed-form torque $M = kT\\,\\cot(\\theta/2)$ and an interlayer shear stress $\\tau = 3kT\\,\\cot(\\theta/2)/(2\\pi R^3)$. It argues that these formulas reproduce measured friction-versus-angle and size-dependent data, and that they imply superlubricity persists to large scales and low temperatures.","feed_headline":"Twist entropy sets the law of superlubric friction","feed_subtitle":"A closed-form torque kT cot(θ/2) from moiré entropy matches experiments and predicts self-rotation.","key_machinery":"The load-bearing machinery is the configuration-entropy ansatz of Eq. (2): the number of atomic microstates in a moiré unit cell is taken proportional to the cell area, $\\Omega_\\theta/\\Omega_0 = S_\\theta/S_0$. Combined with the moiré side-length law $a_m = \\sqrt{3}\\,a/(2\\sin(\\theta/2))$ and the quasi-static thermodynamic relation $dQ = T\\,dS$, it yields the entropy difference $\\Delta S = 2k\\ln(a_m/a)$ and then the torque. The argument also relies on the rigid-layer approximation, justified by the large ratio of in-plane stiffness to interlayer shear modulus, and on the observation that the interlayer potential energy is essentially unchanged during rotation, which sets $dU=0$ and makes the torque purely entropic.","core_discovery":"The central claim is that structural superlubricity is an entropic phenomenon. Relative to AA-stacked bilayer graphene, a twisted bilayer is assigned a configuration entropy $\\Delta S = 2k\\ln(a_m/a)$, where $a_m = \\sqrt{3}\\,a/(2\\sin(\\theta/2))$ is the moiré side length and $a$ is the carbon bond length; the entropy is taken to dominate the Helmholtz free energy. The rotation torque follows as $M = kT\\,\\cot(\\theta/2)$, and dividing by the flake area gives $\\tau = 3kT\\,\\cot(\\theta/2)/(2\\pi R^3)$. Thus the shear stress is nonzero but vanishingly small for large rigid flakes, grows linearly with temperature, and becomes negative for twist angles between 30 and 60 degrees, meaning the flake should spontaneously rotate toward a commensurate registry. The paper backs the derivation with molecular-dynamics evidence that interlayer potential energy is nearly independent of twist angle even under normal loads up to 10 GPa, and with comparison to existing graphite-flake friction, nanotube pullout, and torque-induced reorientation experiments.","pith_inferences":["An extension of the model is to measure the full torque-angle curve for other layered materials such as MoS2 or hexagonal boron nitride; if the $kT\\cot(\\theta/2)$ law holds unchanged, the entropic origin is generic rather than specific to graphene.","The predicted sign change near 30 degrees implies a thermodynamic tendency for twisted contacts to anneal toward commensurate registry; this could mean that preferred twist angles observed in devices partly reflect entropy maximization rather than preparation history alone.","A computational test of the weakest premise is direct enumeration of distinct rigid registries for finite flakes: if the microstate count does not track moiré cell area, the law would acquire finite-size corrections beyond the simple $R^{-3}$ dependence.","At very low temperatures quantum zero-point motion should eventually replace the classical accessible-volume estimate of entropy, so the strictly linear $T$ dependence may saturate; measuring that saturation would delimit the classical regime."],"forward_implications":["Shear stress in a twisted bilayer contact should vanish linearly as temperature approaches zero and grow linearly with temperature, so superlubricity should persist at low temperatures.","For fixed angle and temperature the stress scales as $R^{-3}$; ideal large-area rigid contacts should have immeasurably small shear stress, consistent with measured nanotube pullout strengths.","Between 30 and 60 degrees the torque reverses sign: a free flake should spontaneously rotate toward AB stacking, an entropy-driven self-rotation the paper identifies with observations of torque-induced reorientation.","The formula should remain valid under uniformly distributed normal loads up to 10 GPa, as long as potential energy stays nearly independent of twist angle.","Because the stress expression contains no material-specific interaction strength, the same law should describe other rigid layered contacts with the same geometry."],"supporting_citations":[{"why":"supplies the measured friction-versus-twist-angle data for a graphite flake that the model's curve is compared with at R=1 nm and 300 K.","marker":"[8]"},{"why":"states the original prediction that superlubricity exists for rigid, atomically flat, incommensurate surfaces, the baseline the entropic explanation is meant to replace.","marker":"[5]"},{"why":"provides the measured shear strength for pulling a long inner shell from a double-walled carbon nanotube, used to test the R^{-3} size prediction.","marker":"[13]"},{"why":"reports that superlubricity is eliminated by torque-induced reorientation, cited as experimental evidence for the spontaneous-rotation branch.","marker":"[21]"},{"why":"gives observations of a graphene nanoflake sliding on graphene, including self-termination by rotation and longer sliding at 5 K than at 77 K, supporting the temperature dependence.","marker":"[9]"},{"why":"supplies the moiré period formula a_m = sqrt(3) a / (2 sin(theta/2)) used to convert cell area into twist-angle dependence.","marker":"[25]"},{"why":"independently provides the same moiré side-length formula used in constructing the entropy difference.","marker":"[26]"},{"why":"justifies treating the layers as rigid during sliding, a necessary condition for the purely entropic torque calculation.","marker":"[24]"}],"fun_headline_variants":["Twist entropy dictates torque, not energy","Moiré entropy sets superlubric torque law","Entropy torque: kT cot(theta/2) predicts rotation","Superlubricity is an entropic phenomenon"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Equation (2) assumes that the number of possible atomic configurations inside a moiré unit cell is proportional to the area of that cell; if the microstate count does not scale with moiré cell area, the entropy difference and the derived cotangent law do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Twist entropy dictates torque, not energy","Moiré entropy sets superlubric torque law","Entropy torque: kT cot(theta/2) predicts rotation","Superlubricity is an entropic phenomenon"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000853,"raw_usage":{"total_tokens":3695,"prompt_tokens":921,"completion_tokens":2774,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":537,"completion_tokens_details":{"reasoning_tokens":2710}},"tokens_in":537,"tokens_out":2774,"duration_ms":20314,"temperature":1.0,"reasoning_tokens":2710,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:35:31.009391+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure, for a rigid, defect-free circular graphene flake of radius R twisted against graphite at fixed temperature, the torque as a function of twist angle from 0 to 60 degrees. The paper's claim predicts $M = kT\\,\\cot(\\theta/2)$: divergence near zero angle, linear temperature dependence, no dependence on R beyond the entropy area factor, and a sign reversal at 30 degrees. Observing no sign reversal, or a torque that is not linear in temperature, would falsify Eq. (6); a direct enumeration of finite-flake microstates that fails to reproduce $k\\ln(S_\\theta/S_0)$ would falsify Eq. (2).","supporting_citations":[{"cited_title":"Dienwiebel, G","cited_arxiv_id":null,"evidence_quote":"supplies the measured friction-versus-twist-angle data for a graphite flake that the model's curve is compared with at R=1 nm and 300 K."},{"cited_title":"Hirano and K","cited_arxiv_id":null,"evidence_quote":"states the original prediction that superlubricity exists for rigid, atomically flat, incommensurate surfaces, the baseline the entropic explanation is meant to replace."},{"cited_title":"Zhang, Z","cited_arxiv_id":null,"evidence_quote":"provides the measured shear strength for pulling a long inner shell from a double-walled carbon nanotube, used to test the R^{-3} size prediction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"reports that superlubricity is eliminated by torque-induced reorientation, cited as experimental evidence for the spontaneous-rotation branch."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives observations of a graphene nanoflake sliding on graphene, including self-termination by rotation and longer sliding at 5 K than at 77 K, supporting the temperature dependence."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the moiré period formula a_m = sqrt(3) a / (2 sin(theta/2)) used to convert cell area into twist-angle dependence."},{"cited_title":"Hod, Physical Review B 86, 075444 (2012)","cited_arxiv_id":null,"evidence_quote":"justifies treating the layers as rigid during sliding, a necessary condition for the purely entropic torque calculation."}],"review_version":1}