Pith. sign in

REVIEW 5 major objections 5 minor 15 references

Torsional Behavior of Carbon-Doped Ferrous Nanowires: Atomic-Scale Insights from MD Simulations

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

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

Pith's one-line read Carbon doping weakens iron nanowires under torsion, simulations show.

desk verdict First systematic MD map of Fe-C nanowire torsion, with a credible setup and expected trends, but the shear-stress metric is underdefined and the quantitative claims are not reproducible as written. read the letter →

arxiv 2507.04901 v2 pith:BCAULVVL submitted 2025-07-07 cond-mat.mtrl-sci

Charith L. Hirimuthugodage , Laalitha S.I. Liyanage , K. G. S. H. Gunawardana This is my paper · ORCID
classification cond-mat.mtrl-sci
keywords torsionnanowireiron-carbonmoleculardynamicsMEAMshearstressdislocationLAMMPS
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 uses molecular dynamics simulations to ask how carbon doping, temperature, and cross-sectional size change the way iron nanowires twist before they fail. It shows that adding carbon to a [001] ferrous nanowire lowers the maximum shear stress and the critical torsional angle at which elastic deformation gives way to plastic flow; at 10% carbon the elastic-plastic boundary becomes hard to identify. Raising temperature from 1 K to 900 K has the same weakening effect, while thicker nanowires need more shear stress to reach their smaller critical angle. The paper frames these as the main trends that define the torsional performance of Fe-C nanowires under the MEAM model.

What carries the argument

The central object is a set of circular BCC Fe and Fe-C nanowires with [001] orientation, twisted in LAMMPS with the Liyanage MEAM potential. The loading procedure rotates the two end blocks in opposite directions at 2e11 degrees per second in 1-degree steps with a 20 ps relaxation between steps; the response is tracked as shear stress versus torsional angle, with the critical angle defined as the point of maximum shear stress before plastic flow. Dislocation evolution is identified with Common Neighbor Analysis in OVITO. The MEAM potential is the load-bearing element: the claim that carbon weakens torsion rests entirely on how this potential describes Fe-C interactions under large strain.

What would settle it

Repeat the 1 K, 10%-carbon torsion runs with a different interatomic potential validated against high-strain Fe-C data (for example, one fitted to DFT shear and fracture energies) and with several independent random carbon distributions; if the maximum shear stress and critical angle do not decrease monotonically with carbon content, or if the 10% instability disappears, the claimed carbon-weakening trend is an artifact of the chosen potential rather than a property of Fe-C nanowires.

Watch

Extended reading notes

Core claim

The central claim is that carbon concentration, temperature, and cross-sectional diameter each move the shear stress versus rotation-angle curve of [001] Fe nanowires in a specific direction. Increasing carbon from 0% to 10% reduces the maximum shear stress and the critical torsional angle, and makes the elastic-plastic transition indistinct above 5%. Increasing temperature from 1 K to 900 K also lowers strength and critical angle. Increasing diameter from 10a to 15a raises the shear stress needed to reach the critical angle while reducing that angle. The paper attributes these trends to carbon-induced lattice distortion acting as stress concentrators, thermal vibration enhancing atomic mobility and dislocation nucleation, and the larger strained outer surface of thicker wires.

Load-bearing premise

The entire trend depends on the Liyanage MEAM potential being accurate for Fe-C interactions under torsional loading up to 10% carbon and through plastic failure, even though it was fitted to structural, elastic, and thermal properties, not to torsion or fracture.

Editorial extensions

If this is right

  • Fe-C nanowires intended for torsional components in NEMS or actuators should avoid carbon concentrations above 5% if a clean elastic-plastic transition is needed.
  • At high temperature and high carbon content, the critical torsional angle ceases to be a well-defined design parameter because the stress-angle curves become highly unstable.
  • Larger cross-sectional nanowires tolerate higher shear stress but reach their critical angle earlier, so size alone does not extend the usable elastic rotation range.
  • The same MEAM potential can be used to predict torsional strength of other Fe-C nanostructures, extending the earlier tensile-loading application of this potential.

Reading between the lines

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

  • Because the study uses a single random carbon distribution per concentration, the reported instability at 10% carbon could be sensitive to the specific placement of interstitial atoms; repeating with multiple random seeds or ordered carbon arrangements would test whether the trend is robust.
  • The weakening with carbon opposes the classical solid-solution strengthening seen in bulk steels, suggesting a possible size-dependent crossover: at the nanoscale, carbon-induced dislocation nucleation may dominate over dislocation pinning, which would have implications for miniaturized steel components.
  • The 20 ps relaxation after each 1-degree rotation acts as an effectively slower loading rate; varying this relaxation time could reveal whether the critical angles and stress values are rate-dependent, a natural next step given that the paper notes torsion rate effects in prior copper nanowire studies.
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.

Referee Report

5 major / 5 minor

Summary. The paper reports molecular dynamics simulations of torsional loading of [001]-oriented body-centered cubic iron and carbon-doped iron nanowires using the MEAM potential of Liyanage et al. The parameter sweep is over carbon concentration (0, 1, 5, 10%), temperature (nominally 1 K or 0 K to 900 K), and cross-sectional size (10a, 13a, 15a). The central claims are that increasing carbon concentration lowers the maximum shear stress and the critical torsional angle, that increasing temperature reduces torsional strength, and that larger cross-sections increase the shear stress required for deformation. The paper presents shear stress versus torsional angle curves, critical-angle values, and dislocation/strain visualizations from OVITO.

Significance. If the reported trends are quantitatively reproducible, the paper would provide a useful first mapping of carbon-doping, temperature, and size effects on the torsional response of Fe-C nanowires, with potential relevance to NEMS design. The work has strengths: it uses a previously published MEAM potential fitted to independent structural, elastic, and thermal data; it covers a systematic parameter matrix; and it includes atomic-scale visualization of dislocation and strain patterns. However, the quantitative content is currently not reproducible because the shear-stress observable is not defined, the curves come from single unseeded realizations without error bars, and the temperature labels are inconsistent. The qualitative directions of the trends may be correct, but the evidence as presented is insufficient to distinguish them from sampling noise or from artifacts of an ill-defined stress metric.

major comments (5)
  1. [Methodology, shear-stress calculation] The shear-stress metric used in Figures 3-17 is not defined. The text states only that it is "the sum of the per-atom stress components calculated by the Virial expression implemented in LAMMPS"; it does not state which tensor component is used (e.g., xy, the longitudinal shear components, von Mises, or the trace), how the per-atom values are aggregated, or whether the result is volume-normalized. LAMMPS per-atom virial stress has units of energy (or pressure-volume depending on output) and contains a kinetic term, so different choices change both the absolute magnitudes and the relative ordering of the curves. Since maximum shear stress is the basis for the critical-angle definition and for the carbon-weakening conclusion, the central quantitative claims cannot be tested until this quantity is specified and, ideally, validated against a torque-based shear-stress definition.
  2. [Methodology, simulation protocol and Figures 3-17] Each reported curve appears to come from a single unseeded MD run. Carbon atoms are said to be randomly distributed, but no random seed or ensemble of independent configurations is described, and no error bars or convergence tests are reported. Given the large fluctuations in the 10% carbon and high-temperature curves, the observed differences between concentrations and sizes could reflect sampling noise rather than systematic trends. Please provide multiple independent realizations per condition and report the mean and spread (or at least a convergence check for the key 0% versus 10% comparison at the low-temperature baseline).
  3. [Effect of Temperature on Torsion] The temperature axis is internally inconsistent. The abstract and Methodology state 1 K; the Effect of Temperature section states "Simulations were conducted at 0 K, 300 K, 600 K, and 900 K"; the figure captions and text report values at 1 K; and Figure 13 is described at 0 K. This must be reconciled because the low-temperature baseline anchors the temperature trend. If the simulations are at 1 K, relabel all text and figures; if some are at 0 K, the data and thermostat settings need to be clarified.
  4. [Methodology, critical torsional angle definition] The critical angle is defined as "the point at which maximum shear stress is reached before plastic deformation begins," but no operational rule is given for curves that have no clear maximum, as admitted for 10% carbon and high temperatures. Without an automated or at least reproducible peak/no-peak criterion, the reported critical-angle values (e.g., 137°, 140°, 100°, 114°) cannot be independently extracted from the raw curves. Please define the algorithm used to assign critical angles or state explicitly which curves were excluded from this determination.
  5. [Methodology, MEAM potential choice] The central results rely on extrapolating the Liyanage et al. MEAM potential to high carbon concentrations (10%) and to plastic failure under torsion, neither of which is part of the potential's fitted property set (structural, elastic, and thermal properties of BCC Fe and Fe-C phases). This is a correctness-risk concern, not a circularity claim. A concrete check would be to compare one representative stress-angle curve (e.g., 0% carbon, 10a, 1 K) against an alternative Fe potential or, if available, against experimental or DFT data on Fe shear behavior; at minimum, state the expected validity range of the potential for torsion and failure.
minor comments (5)
  1. [Abstract] The abstract states that "increasing carbon content weakens grain boundaries," but the simulated nanowires are single-crystal BCC without grain boundaries; the observed weakening is attributed to interstitial lattice distortion and should be reworded to avoid implying a grain-boundary mechanism.
  2. [Header] The preprint header says "Preprint submitted to Computational Material Scinece"; this should read "Computational Materials Science."
  3. [Methodology and Figure 6 caption] The cross-sectional size is described as "diameter" in the Methodology but as "radius" in the text and in the Figure 6 caption; the terminology should be made consistent throughout.
  4. [Effect of Temperature on Torsion] The section discusses potential-energy fluctuations and Figures 7 and 8 appear in that section, but the text never explicitly refers to those figures; please add the citations or remove the figures.
  5. [General] There are numerous grammar and spacing errors (e.g., "Figure 6, further illustrates," "di fferent," "Y , Z"); a careful proofreading pass is needed.

Circularity Check

0 steps flagged · score 2.0 of 10

No circular derivation: the torsion trends are MD outputs, not inputs; only a minor non-load-bearing self-citation lowers the score.

full rationale

The paper reports molecular dynamics simulation outcomes for Fe and Fe–C nanowires under torsion, and none of the claimed trends (carbon weakens torsional capacity, temperature lowers strength, larger cross-sections increase shear stress) is imposed by construction. The paper introduces no fitted parameter in this work that is later renamed as a prediction; the interatomic potential from Liyanage et al. [12] is a self-citation by the corresponding author, but that potential was fitted to structural, elastic, and thermal properties of BCC Fe and Fe–C phases, not to torsion, so the torsional response is not an input to the fit. The second same-group citation [3] is used only to state prior usage of the potential and does not carry any load in the derivation. The definitions of shear stress and critical torsional angle are operational and somewhat under-specified, but they do not encode carbon concentration, temperature, or cross-sectional size, so the reported trends are not true by definition. No equation in the manuscript reduces a claimed result to an input, and no uniqueness theorem or ansatz is imported through citation to forbid alternative outcomes. The genuine concerns—MEAM potential transferability to plastic torsion and the incomplete definition of the virial-based shear stress metric—are correctness and reproducibility risks, not circularity.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

No new particles, forces, phases, or conserved quantities are postulated; the simulation relies on an existing potential and standard MD protocol. The free parameters capture the fitted interatomic potential and the unreported random carbon seed, both of which the central trends depend on.

free parameters (2)
  • MEAM Fe-C potential parameters = not tabulated here; from Ref. [12]
    The potential parameters were fitted to independent Fe and Fe-C data, but every simulated stress value depends on them and the paper never re-validates them under torsion.
  • Random carbon distribution seed = not reported
    Carbon atoms are placed randomly, yet only one realization per concentration appears to be simulated; trends might change with a different seed.
assumptions (4)
  • domain assumption The initial nanowire is a perfect BCC Fe lattice with C atoms in interstitial sites.
    Real nanowires contain surfaces, possible oxides, pre-existing defects, and different carbon distributions; the idealized structure is assumed in the model construction (Methodology).
  • domain assumption The MEAM potential transfers to high-strain torsion and 10 percent carbon.
    The potential was fitted to structural, elastic, and thermal properties in Ref. [12]; accuracy under torsion is not tested, yet the entire study depends on it.
  • domain assumption NVT thermostat with stepwise 1-degree rotations and 20 ps relaxation represents quasi-static torsion.
    No convergence test is shown; the choice of loading protocol can shift critical angles and stress values.
  • domain assumption The sum of per-atom virial stresses equals the torsional shear stress.
    The paper does not specify which stress component is summed, so the reported shear stress may not correspond to the resolved shear stress on the cross-section.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Torsional Behavior of Carbon-Doped Ferrous Nanowires: Atomic-Scale Insights from MD Simulations." pith.science (2026). https://pith.science/paper/BCAULVVL

@misc{pith2026250704901,
  author       = {Pith},
  title        = {Pith review of: Torsional Behavior of Carbon-Doped Ferrous Nanowires: Atomic-Scale Insights from MD Simulations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/BCAULVVL}},
  note         = {Machine review of arXiv:2507.04901}
}
read the original abstract

This study investigates the torsional mechanical properties of pristine iron (Fe) and carbon-doped iron (FeC) nanowires with [001] orientation through molecular dynamics simulations utilizing the Modified Embedded Atom Method (MEAM) potential developed by Liyanage et al. for accurately modeling Fe-C interactions in body-centered cubic structures. Systematic analysis across carbon concentrations (0 - 10%), temperatures (1 - 900 K), and cross-sectional dimensions 10a, 13a, 15a, ( where a = 2.81 Angstrom represents the lattice constant ) within the LAMMPS environment reveals that increasing carbon content weakens grain boundaries, reducing the maximum shear stress required to reach the critical torsional angle, while higher temperatures promote phase transitions from elastic to plastic deformation due to enhanced atomic vibrations, and larger cross-sections exhibit higher shear stress resistance attributed to strengthening effects from the outer atomic layers. By elucidating these interrelationships between carbon content, temperature, and dimensional factors, this work provides fundamental insights into the mechanical behavior of FeC nanowires under torsional loading conditions, offering valuable guidance for their potential applications in nanoelectromechanical systems, nanorobotic actuators, and advanced structural materials where precise control of mechanical properties is essential.

Figures

Figures reproduced from arXiv: 2507.04901 by the authors.

Figure 1
Figure 1. BCC Fe unit cell with carbon impurity in interstitial position. [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Schematic of the Fe nanowire model. Yellow regions indicate torsional motion, and red regions deform naturally. [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Shear stress vs. torsional angle for 10a nanowire at 1 K. Each curve corresponds to a specific carbon concentration (0%, 1%, 5%, 10%) [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (15 more)
Figure 4
Figure 4. Figure 4: Shear stress vs. torsional angle for 13a nanowire at 1 K. Each curve corresponds to a specific carbon concentration (0%, 1%, 5%, 10%) [PITH_FULL_IMAGE:figures/full_fig_p006_4.png]
Figure 5
Figure 5. Figure 5: Shear stress vs. torsional angle for 15a nanowire at 1 K. Each curve corresponds to a specific carbon concentration (0%, 1%, 5%, 10%) [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]
Figure 6
Figure 6. Figure 6: Dislocation structures of the nanowire with a radius of 10 at di [PITH_FULL_IMAGE:figures/full_fig_p008_6.png]
Figure 7
Figure 7. Figure 7: Potential Energy vs. Rotation Angle for 0 % Carbon Content in a 10a Diameter Nanowire at 1 K [PITH_FULL_IMAGE:figures/full_fig_p008_7.png]
Figure 8
Figure 8. Figure 8: Potential energy vs. rotation angle for 0% carbon with a 10 [PITH_FULL_IMAGE:figures/full_fig_p009_8.png]
Figure 9
Figure 9. Figure 9: Shear stress vs. torsional angle for a 10 [PITH_FULL_IMAGE:figures/full_fig_p009_9.png]
Figure 10
Figure 10. Figure 10: Shear stress vs. torsional angle for a 10 [PITH_FULL_IMAGE:figures/full_fig_p010_10.png]
Figure 11
Figure 11. Figure 11: Shear stress vs. torsional angle for a 10 [PITH_FULL_IMAGE:figures/full_fig_p010_11.png]
Figure 12
Figure 12. Figure 12: Shear stress vs. torsional angle for a 10 [PITH_FULL_IMAGE:figures/full_fig_p011_12.png]
Figure 13
Figure 13. Figure 13: Dislocation structures of the model at a 137° torsional angle under varying temperatures. [PITH_FULL_IMAGE:figures/full_fig_p011_13.png]
Figure 14
Figure 14. Figure 14: Shear stress vs. torsional angle at 1 K for nanowires with 0% carbon concentration. Each curve corresponds to a di [PITH_FULL_IMAGE:figures/full_fig_p012_14.png]
Figure 15
Figure 15. Figure 15: Shear stress vs. torsional angle at 1 K for nanowires with 1% carbon concentration. Each curve corresponds to a di [PITH_FULL_IMAGE:figures/full_fig_p013_15.png]
Figure 16
Figure 16. Figure 16: Shear stress vs. torsional angle at 1 K for nanowires with 5% carbon concentration. Each curve corresponds to a di [PITH_FULL_IMAGE:figures/full_fig_p013_16.png]
Figure 17
Figure 17. Figure 17: Shear stress vs. torsional angle at 1 K for nanowires with 10% carbon concentration. Each curve corresponds to a di [PITH_FULL_IMAGE:figures/full_fig_p014_17.png]
Figure 18
Figure 18. Figure 18: Strain distribution across the cross-section of a 10 a nanowire. [PITH_FULL_IMAGE:figures/full_fig_p014_18.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

15 extracted references · 15 canonical work pages

  1. [1]

    Khalid, X

    K. Khalid, X. Tan, H. F. Mohd Zaid, Y . Tao, C. Lye Chew, D.-T. Chu, M. K. Lam, Y .-C. Ho, J. W. Lim, L. Chin Wei, Advanced in developmental organic and inorganic nanomaterial: a review, Bioengineered 11 (1) (2020) 328–355

    work page 2020

  2. [2]

    Kehoe, P

    M. Kehoe, P. Kelly, The role of carbon in the strength of ferrous martensite, Scripta Metallurgica 4 (6) (1970) 473–476

    work page 1970

  3. [3]

    Nadeesha, J

    T. Nadeesha, J. Liyanage, L. S. Liyanage, Prediction of mechanical properties of steel nanowires using molecular dynamics (2021)

    work page 2021

  4. [4]

    Sainath, B

    G. Sainath, B. Choudhary, Atomistic simulations on ductile-brittle transition in¡ 111¿ bcc fe nanowires, Journal of Applied Physics 122 (9) (2017)

    work page 2017

  5. [5]

    Z. Yang, M. Li, Y . Li, Y . Yang, J. Zhao, Molecular dynamics simulation on torsion deformation of copper aluminum core–shell nanowires, Journal of Nanoparticle Research 23 (2021) 1–10

    work page 2021

  6. [6]

    Z. Yang, M. Li, Y . Li, Y . Yang, X. Ding, Comparisons in mechanical properties of various copper nanowires under torsion: molecular dynamics simulation, Journal of Nanoparticle Research 24 (2) (2022) 42

    work page 2022

  7. [7]

    Y . Yang, Y . Li, Z. Yang, G. Zhang, X. Wang, J. Liu, Molecular dynamics simulation on the elastoplastic properties of copper nanowire under torsion, Journal of Nanoparticle Research 20 (2018) 1–10

    work page 2018

  8. [8]

    Y . Gao, F. Wang, T. Zhu, J. Zhao, Investigation on the mechanical behaviors of copper nanowires under torsion, Computational Materials Science 49 (4) (2010) 826–830

    work page 2010

Show all 15 references
  1. [9]

    Sung, C.-D

    P.-H. Sung, C.-D. Wu, T.-H. Fang, E ffects of temperature, loading rate and nanowire length on torsional deformation and mechanical prop- erties of aluminium nanowires investigated using molecular dynamics simulation, Journal of Physics D: Applied Physics 45 (21) (2012) 215303

    2012

  2. [10]

    S. Li, E. K. Salje, S. Jun, X. Ding, Large recovery of six-fold twinned nanowires of α-fe, Acta Materialia 125 (2017) 296–302

    2017

  3. [11]

    S. J. Plimpton, A. P. Thompson, Computational aspects of many-body potentials, MRS bulletin 37 (5) (2012) 513–521

    2012

  4. [12]

    L. S. Liyanage, S.-G. Kim, J. Houze, S. Kim, M. A. Tschopp, M. I. Baskes, M. F. Horstemeyer, Structural, elastic, and thermal properties of cementite (fe 3 c) calculated using a modified embedded atom method, Physical Review B 89 (9) (2014) 094102

    2014

  5. [13]

    HARDEN, R

    C. HARDEN, R. LEHOUCQ, Convergence verification of static solvers via energy minimization in lammps, CSRI SUMMER PROCEED- INGS 2008 (2008) 199

    2008

  6. [14]

    Miyazaki, Y

    N. Miyazaki, Y . Shiozaki, Calculation of mechanical properties of solids using molecular dynamics method, JSME international journal. Ser. A, Mechanics and material engineering 39 (4) (1996) 606–612

    1996

  7. [15]

    A. Stukowski, Visualization and analysis of atomistic simulation data with ovito–the open visualization tool, Modelling and simulation in materials science and engineering 18 (1) (2009) 015012. 13

    2009

Pith tools

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