Pith. sign in

REVIEW 3 major objections 4 minor 30 references

Angular Dependence of Four Mechanical Properties of Graphynes

T0 review · 3 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash

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

Pith's one-line read The paper shows that in asymmetric graphynes, the Young's modulus can change by a factor of about 10 with loading direction, and some structures exhibit negative and null linear compressibility at specific angles.

desk verdict A systematic but unvalidated extension of prior MD work; the angular maps are interesting, but the rotation equation is under-specified and all claims rest on AIREBO. read the letter →

arxiv 2412.05705 v1 pith:JFUSDRP5 submitted 2024-12-07 cond-mat.mtrl-sci

Guilherme B. Kanegae , Alexandre F. Fonseca This is my paper · ORCID
classification cond-mat.mtrl-sci
keywords graphyneelasticanisotropyYoung'smoduluslinearcompressibilityPoisson'sratioshearAIREBOpotentialtwo-dimensionalmaterials
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

The paper works out how the four main elastic properties of 70 graphyne structures change with the direction of applied stress. It finds that in asymmetric graphynes, especially those of family GnY5, the Young's modulus can change by a factor of about 10 depending on loading direction: it is near 4 N/m along the armchair direction and rises to roughly 40 N/m at 60 degrees, where the load falls directly on the acetylenic chains. GnY5 and GnY6 also show negative linear compressibility in some directions, with exactly zero linear compressibility at 60°, 120°, 240°, and 300°. If these results are correct, graphyne sheets are highly anisotropic two-dimensional mechanical metamaterials whose directional stiffness and compressibility can be tuned by chain length and connection pattern.

What carries the argument

The load-bearing object is the 2D elastic constant matrix C with components C11, C12, C22, C66 taken from the authors' earlier study [22]. The paper rotates it according to C′ = R(θ)·C·R(θ), where R(θ) is the plane rotation matrix, and then computes E, G, ν, and β from the rotated coefficients using standard orthotropic formulas. The rotation makes the angular dependence explicit; the independence of the shear modulus follows from a known relation among the Cij from Ref. [28]. The stretching-versus-hinging interpretation, following Ref. [30], explains why stiffness peaks along the acetylenic chain direction and why the polar plots become more asymmetric as n increases.

What would settle it

A density-functional-theory calculation of the elastic constants for GnY5 and GnY6 with n = 1 and n = 2, followed by the same rotation analysis, would settle the claim: if the computed maximum/minimum Young's modulus ratio is far below 10, or if the linear compressibility never changes sign, then the factor-10 anisotropy and the negative linear compressibility are artifacts of the force field.

Watch

Extended reading notes

Core claim

The central claim is that asymmetry in the graphyne lattice produces a strong angular dependence in the in-plane elastic response. Taking the elastic constants from a previous simulation study, the authors rotate the two-dimensional stiffness matrix and recompute Young's modulus, shear modulus, Poisson's ratio, and linear compressibility as functions of angle. They find the shear modulus is independent of orientation for all structures, but Young's modulus, Poisson's ratio, and linear compressibility vary strongly for the asymmetric families GnY2, GnY3, GnY5, and GnY6. In GnY5 the maximum-to-minimum Young's modulus ratio reaches about 10, and in GnY5 and GnY6 the linear compressibility reverses sign with direction, passing through zero at exactly 60°, 120°, 240°, and 300°. The paper interprets the angular behavior as a competition between stretching and hinging of the acetylenic chains: stretching dominates along the chain direction (60°) and hinging dominates along the armchair and zigzag directions, which is also why the asymmetry grows as the chain number n increases.

Load-bearing premise

All angular predictions are inherited from the elastic constants computed in the previous study using the AIREBO classical potential, and the paper provides no quantum-mechanical or experimental check of those constants.

Editorial extensions

If this is right

  • Graphyne membranes could serve as direction-tunable mechanical elements: stiff along the acetylenic chain direction and soft along armchair/zigzag directions, with the contrast set by the choice of n.
  • The directions of null linear compressibility mean that under hydrostatic-like in-plane pressure, GnY5 and GnY6 sheets would not contract along those directions, which could be exploited for dimensionally stable nanoscale components.
  • Because the shear modulus is independent of angle, graphyne-based composites would not need orientation control to have isotropic shear response, simplifying design.
  • Poisson's ratio values as high as about 2.5, combined with the large angular variation, imply that loading along one direction can produce very large perpendicular strain, a feature relevant for auxetic-type mechanical metamaterials.

Reading between the lines

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

  • If the null-compressibility angles are a geometric consequence of the triangular arrangement of acetylenic chains, other chain-based 2D carbon allotropes with different chain lengths might exhibit zero linear compressibility at the same 60-degree multiples.
  • A finite-temperature molecular dynamics study stretching these sheets along several directions could test whether the factor-of-10 anisotropy and the sign change of linear compressibility survive beyond the 0 K stiffness matrix.
  • The same rotation-and-recompute procedure could be applied to other stress-strain metrics, such as second-order elastic constants under uniaxial strain, to see whether the angular zeros shift with applied strain.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper reports the angular dependence of four in-plane elastic properties (Young's modulus, shear modulus, Poisson's ratio, and linear compressibility) for 70 graphyne structures, with emphasis on asymmetric families GnY2, GnY3, GnY5, and GnY6. Elastic constants Cij are taken from the authors' prior AIREBO-based molecular dynamics study (Ref. [22]), and the present work rotates the stiffness matrix using Eq. (2) to obtain polar plots. The headline results are a maximum-to-minimum Young's modulus ratio near 10 for GnY5, shear modulus independent of angle for all structures, and negative or null linear compressibility along some directions in GnY5 and GnY6.

Significance. If the angular predictions are correct, the paper provides a useful and simple extension of prior elastic-constant data, identifying graphynes as highly anisotropic 2D mechanical metamaterials with direction-dependent auxetic or negative-linear-compressibility behavior. The stretching-and-hinging interpretation is plausible and qualitatively connects the results to known mechanisms. However, the significance is currently limited because the new results are deterministic transformations of previously published Cij, with no independent validation, no error bars, and a questionable rotation formula. The strongest potential contribution is the identification of specific directions with null linear compressibility, but this claim depends critically on the correctness of the tensor rotation.

major comments (3)
  1. [Theory and Simulation Details, Eq. (2)] Equation (2), C' = R(θ)·C·R(θ), is not the standard transformation law for a fourth-rank stiffness tensor in Voigt notation. In 2D, the correct transformation requires a 3×3 Bond-type matrix M(θ), giving C' = M·C·M^T (or an equivalent form), not a similarity transformation with an ordinary rotation matrix. The authors do not define R(θ) or provide the Bond matrix. If Eq. (2) was actually implemented in the calculations, then all angular results in Figures 2–4—the factor-10 Young's modulus anisotropy, the negative/null linear compressibility, and the claimed shear-modulus independence—are unsupported artifacts. The authors must either correct Eq. (2) and show that their numerical implementation uses the proper tensor transformation, or supply the explicit Bond matrix and an internal consistency check such as equality of E(θ) computed from stiffness and from compliance. This is the central load-bearing point of the paper.
  2. [Results and Discussion, Figure 2 and shear-modulus independence] The claim that the shear modulus is independent of θ for all structures, symmetric or not, is asserted via an unspecified relation and Ref. [28], with no derivation or numerical check. For a general 2D orthotropic material, C66'(θ) is constant only under a specific relation among C11, C22, C12, and C66; the expression for G in Eq. (1) is not the condition that guarantees isotropy of the shear modulus under rotation. Since the shear-modulus plots in Figure 2 are presented as a central result, the authors should provide the explicit derivation (or reference containing it), state the required relation, and verify it against the Cij values from Ref. [22] for all 70 structures.
  3. [Entire manuscript, angular results inherited from Ref. [22]] All angular predictions are obtained by rotating the Cij values computed in Ref. [22] using the AIREBO classical potential. The paper provides no error bars on these constants, no comparison with DFT or experiment, and no sensitivity analysis. Given that the headline findings (factor-10 anisotropy, negative linear compressibility) are deterministic functions of these Cij, their physical relevance depends entirely on the accuracy of AIREBO for graphynes. The authors should include at least one benchmark (e.g., DFT calculation of Cij for a representative GnY5 or GnY6 structure) or a clear discussion of the known limitations of AIREBO for acetylenic carbon chains, and they should report uncertainties in the plotted quantities.
minor comments (4)
  1. [Theory and Simulation Details, Eq. (1)] Equation (1) is presented without citation or derivation; the definition of linear compressibility β_x should be stated explicitly (e.g., β = S11 + S12 in 2D) and consistent with the sign convention used in the plots.
  2. [Results and Discussion, Figure 4] The claim of null linear compressibility at exactly 60°, 120°, 240°, and 300° is based on visual inspection of polar and contour plots. The authors should provide a numerical tolerance (e.g., |β| < 10^-3 N/m per unit stress) or a table of β values at these angles to support the claim of exact zeros.
  3. [References] Reference [4] lists the journal as 'J. Chem. Phys. C'; the correct abbreviation is 'J. Phys. Chem. C'. Several reference titles contain typographical spacing issues (e.g., Ref. [10] 'Chem. Commun. 46, 3256' is fine, but Ref. [13] has a misplaced space in 'Nano Energy 43, 192').
  4. [Data Availability] The statement 'Data available on reasonable request from the authors' is weak for a computational paper. The authors should provide the Cij values for all 70 structures and, if possible, the code/scripts used to perform the rotation, to allow readers to reproduce Figures 2–4.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity identified; the angular curves are direct tensor rotations of the authors' previously published elastic constants, and no fitted quantity is relabeled as a prediction.

full rationale

The derivation chain is transparent: Cij are taken from the authors' prior study (Ref. [22]) through AIREBO/LAMMPS simulations, and Eq. (2) rotates the stiffness matrix to obtain angular-dependent values of E, G, beta, and nu using Eq. (1). The angular quantities are mathematical functions of Cij, not quantities that were used to fit or define Cij, so the present results do not reduce to their inputs by construction. The self-citation to Ref. [22] is load-bearing in the sense that the new results inherit the earlier elastic constants, but it is a normal transfer of previously computed numerical data with stated assumptions (the AIREBO potential), not an unverified uniqueness theorem or an ansatz smuggled in to force the present conclusions. No step in this paper defines Cij in terms of E(theta), G(theta), beta(theta), or nu(theta), and no parameter is fitted to the angular quantities. A substantive correctness concern lies outside circularity: Eq. (2) as written, C' = R*C*R, is not the standard fourth-rank stiffness rotation in Voigt notation, which requires a Bond-type matrix M with C' = M*C*M^T; if implemented literally, the plotted angular dependence would be erroneous. That issue, and the absence of independent validation of the AIREBO-based Cij, are correctness and robustness risks, not circularity. The paper does not invoke an external benchmark, but the absence of a benchmark does not by itself make the derivation circular.

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

The paper introduces no new free parameters itself. It inherits the AIREBO potential parameters and the Cij values from Ref 22, and it relies on standard but under-specified tensor-rotation formulas. The central assumptions are that the force field and the earlier MD results are accurate, and that the shear-modulus invariance relation holds. No new entities are postulated.

free parameters (1)
  • AIREBO C-C force-field parameters (inherited from Ref 25) = Not enumerated in this paper
    All elastic constants come from LAMMPS/AIREBO simulations in Ref 22; these empirical parameters were fitted to carbon and hydrocarbon properties and are not re-calibrated or validated for graphynes here.
assumptions (4)
  • standard math Rotation of the elastic stiffness matrix via Eq. (2) yields the correct fourth-rank tensor transformation
    The paper writes C' = R(theta)*C*R(theta) but does not define the explicit Bond rotation matrix; if the Voigt transformation is different, the angular maps change.
  • domain assumption The AIREBO potential accurately predicts elastic properties of graphynes
    No DFT or experimental benchmarks are provided; all angular predictions inherit this assumption from Ref 22.
  • domain assumption Cij values from Ref 22 are reliable and correctly represent equilibrium graphyne structures
    The current paper does not repeat simulations; it reuses published Cij without sensitivity analysis.
  • domain assumption The shear modulus invariance relation C66 = (C11 - 2C12 + C22)/4 holds for the computed structures
    This relation is cited from Ref 28 but not demonstrated for the present Cij data; it underlies the claim that the shear modulus is angle-independent.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Angular Dependence of Four Mechanical Properties of Graphynes." pith.science (2026). https://pith.science/paper/JFUSDRP5

@misc{pith2026241205705,
  author       = {Pith},
  title        = {Pith review of: Angular Dependence of Four Mechanical Properties of Graphynes},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JFUSDRP5}},
  note         = {Machine review of arXiv:2412.05705}
}
read the original abstract

Graphyne is a porous two-dimensional carbon allotrope of graphene that possesses interesting physical properties, including non-null bandgap. It is composed of carbon hexagonal rings or carbon-carbon bonds connected by acetylenic chains. The diverse forms of these connections yield a variety of graphyne structures. In a previous study, we have obtained the elastic properties of seven distinct families of graphyne structures as a function of the number of acetylene chains, from 1 to 10. The Young's modulus, shear modulus, Poisson's ratio and linear compressibility were predicted for the zigzag and armchair directions of all 70 graphyne structures. Here, we present some noteworthy findings regarding the angular dependence of these four elastic properties of asymmetric graphynes. Our results demonstrate that in a single structure, the minimum and maximum Young's modulus can vary by a factor of 10. Additionally, the directions of null linear compressibility in some asymmetric structures were determined.

Figures

Figures reproduced from arXiv: 2412.05705 by the authors.

Figure 1
Figure 1. Left: superposition of the graphyne structures belonging to the first family, GnY1, with different numbers of acetylenic chains, n. The cyan, orange, green and red colors correspond to n = 1, 2, 3 and 10, respectively. Right: the n = 1 members of the other six families of GYs as originally proposed by Baughman, Eckhardt and Kertesz [1]. For all structures, the GY’s armchair and zigzag directions are drawn along the … view at source ↗
Figure 2
Figure 2. Polar plots of the Young's modulus (left panel) of the symmetric GnY1 family and the shear modulus (right panel) of the asymmetric GnY5 family. The line colors correspond to structures with varying numbers of acetylene chains. RESULTS AND DISCUSSION As anticipated, the most intriguing findings emerge from the asymmetric GYs, namely families GnY2, GnY3, GnY5 and GnY6 [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Young’s modulus (top), Poisson’s ratio (middle) and linear compressibility (bottom) of GY structures from GnY3, GnY5 and GnY6 families. Positive and negative values of the linear compressibility of GYs from GnY5 and GnY6 families are drawn in full and dashed lines, respectively. The color code is the same as in [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

30 extracted references · 19 canonical work pages

  1. [22]

    G. B. Kanegae and A. F. Fonseca, Carbon Trends 7,100152 (2022) https://doi.org/10.1016/j.cartre.2022.100152

    arXiv 2022

  2. [28]

    H. B. Huntington, in SOLID STATE PHYSICS, Advances in Research and Applications, Volume 7, ed. by F. Seitz and D. Turnbull (Academic Press Inc. 1958), p. 221

    work page 1958

  3. [1]

    R. H. Baughman, H. Eckhardt and M. Kertesz, J. Chem. Phys. 87, 6687 (1987) https://doi.org/10.1063/1.453405

    doi:10.1063/1.453405 1987

  4. [2]

    A. L. Ivanovskii, Progress in Solid State Chemistry 41, 1 (2013) https://doi.org/10.1016/j.progsolidstchem.2012.12.001

    doi:10.1016/j.progsolidstchem.2012.12.001 2013

  5. [3]

    S. W. Cranford and M. J. Buehler, Carbon 49, 4111 (2011) https://doi.org/10.1016/j.carbon.2011.05.024

    work page doi:10.1016/j.carbon.2011.05.024 2011

  6. [4]

    J. E. Padilha, A. Fazzio and A. J. R. da Silva, J. Chem. Phys. C 118, 18793 (2014) https://doi.org/10.1021/jp5062804

    work page doi:10.1021/jp5062804 2014

  7. [5]

    Galhofo and N

    D. Galhofo and N. Silvestre, Mechanics of Advanced Materials and Structures 28, 495 (2021) https://doi.org/10.1080/15376494.2019.1578007

    arXiv 2021

  8. [6]

    X. –M. Wang, D. –C. Mo and S. –S. Lu, J. Chem. Phys. 138, 204704 (2013) https://doi.org/10.1063/1.4806069

    work page doi:10.1063/1.4806069 2013

Show all 30 references
  1. [7]

    S. A. Hernandez and A. F. Fonseca, Diamond and Related Materials 77, 57 (2017) https://doi.org/10.1016/j.diamond.2017.06.002

    2017 doi

  2. [8]

    G. B. Kanegae, M. L. Pereira Junior, D. S. Galvão, L. A. Ribeiro Junior and A. F. Fonseca, ACS Appl. Mater. Interfaces (2024) https://doi.org/10.1021/acsami.4c03302

    2024 doi

  3. [9]

    Narita, S

    N. Narita, S. Nagai, S. Suzuki and K. Nakao, Phys. Rev. B 58, 11009 (1998) https://doi.org/10.1103/PhysRevB.58.11009

    1998 doi

  4. [10]

    G. Li, Y. Li, H. Liu, Y. Guo, Y. Li and D. Zhu, Chem. Commun. 46, 3256 (2010) https://doi.org/10.1039/b922733d

    2010 doi

  5. [11]

    Matsuoka, R

    R. Matsuoka, R. Sakamoto, K. Hoshiko, S. Sasaki, H. Masunaga, K. Nagashio and H. Nishihara, J. Am. Chem. Soc. 139,3145 (2017) https://doi.org/10.1021/jacs.6b12776

    2017 doi

  6. [12]

    K. Khan, A. K. Tareen, M. Iqbal, Z. Shi, H. Zhang and Z. Guo, Nano Today 39, 101207 (2021) https://doi.org/10.1016/j.nantod.2021.101207

    2021

  7. [13]

    J. Gao, J. Li, Y. Chen, Z. Zuo, Y. Li, H. Liu and Y. Lia, Nano Energy 43, 192 (2018) https://doi.org/10.1016/j.nanoen.2017.11.005

    2018 doi

  8. [14]

    Q. Pan, S. Chen, C. Wu, F. Shao, J. Sun, L. Sun, Z. Zhang, Y. Man, Z. Li, L. He and Y. Zhao, CCS Chem 2, 1368 (2020) https://doi.org/10.31635/ccschem.020.202000377

    2020 doi

  9. [15]

    V. G. Desyatkin, W. B. Martin, A. E. Aliev, N. E. Chapman, A. F. Fonseca, D. S. Galvão, E. R. Miller, K. H. Stone, Z. Wang, D. Zakhidov, F. T. Limpoco, S. R. Almahdali, S. M. Parker, R. H. Baughman and V. O. Rodionov, J. Am. Chem. Soc. 144, 17999 (2022) https://doi.org/10.102...

    2022 doi

  10. [16]

    X. Liu, S. M. Cho, S. Lin, Z. Chen, W. Choi, Y. -M. Kim, E. Yun, E. H. Baek, D. H. Ryu and H. Lee, Matter 5, 2306 (2022) https://doi.org/10.1016/j.matt.2022.04.033

    2022 doi

  11. [17]

    Barua, A

    M. Barua, A. Saraswat, C. N. R. Rao, Carbon 200, 247 (2022) https://doi.org/10.1016/j.carbon.2022.08.061

    2022 doi

  12. [18]

    Z. Yang, Y. Zhang, M. Guo and J. Yun, Comput. Mater. Sci. 160, 197 (2019) https://doi.org/10.1016/j.commatsci.2018.12.033

    2019 doi

  13. [19]

    Lin and M

    S. Lin and M. J. Buehler, Nano 5, 11801 (2013) https://doi.org/10.1039/C3NR03241H

    2013 doi

  14. [20]

    Esmaeili, J

    C. Esmaeili, J. Song, Y. Li, L. Mao, H. Liu, F. Wu, Y. Chen, C. Liu and X. -E. Zhang, J. Electrochem. Soc. 168, 077520 (2021). https://doi.org/10.1149/1945-7111/ac139c

    2021 doi

  15. [21]

    Z. Zhu, Q. Bai, S. Li, S. Li, M. Liu, F. Du, N. Sui, W. W. Yu, Small 16, 2001440 (2020) https://doi.org/10.1002/smll.202001440

    2020 doi

  16. [23]

    G. B. Kanegae and A. F. Fonseca, MRS Advances 8, 355 (2023) https://doi.org/10.1557/s43580-023-00529-x

    2023 doi

  17. [24]

    A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, Comput. Phys. Commun. 271, 108171 (2022) https://doi.org/10.1...

    2022

  18. [25]

    D. W. Brenner, O. A. Shenderova, J. A. Harrison, S. J. Stuart, B. Ni and S. B. Sinnott, J. Phys.: Condens. Matter 14, 783 (2002) https://doi.org/10.1088/0953-8984/14/4/312

    2002 doi

  19. [26]

    P. V. Polyakova, R. T. Murzaev, D. S. Lisovenko and J. A. Baimova, Comput. Mater. Sci. 244, 113171 (2024) https://doi.org/10.1016/j.commatsci.2024.113171

    2024

  20. [27]

    K. J. Kotoko, K. Sodoga, Y. Shaidu, N. Seriani, S. Borah and K. Beltako, J. Phys. Chem. C 128, 17058 (2024) https://doi.org/10.1021/acs.jpcc.4c01233

    2024 doi

  21. [29]

    M. Liu, V. I. Artyukhov, H. Lee, F. Xu and B. I. Yakobson, ACS Nano 7, 10075 (2013) https://doi.org/10.1021/nn404177r

    2013 doi

  22. [30]

    J. N. Grima, D. Attard, R. Caruana-Gauci, R. Gatt, Scripta Materialia 65, 565 (2011) https://doi.org/10.1016/j.scriptamat.2011.06.011

    2011 doi

Pith tools

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