pith. sign in

arxiv: 2604.11367 · v1 · submitted 2026-04-13 · ❄️ cond-mat.mes-hall · cond-mat.mtrl-sci

Strain-Induced Curvature in Monolayer Graphene: Effects on Electronic Structure, Phonon Dynamics, and Lattice Thermal Conductivity

M. C. Santos , E. Lora da Silva , D. S. Baptista , T. Santos , M. Molinari , F. J. Manj\'on , Yin Cui , Xidong Lin
show 1 more author
Tao Yang
This is my paper

Pith reviewed 2026-05-10 15:59 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall cond-mat.mtrl-sci
keywords graphenestraincurvaturelattice thermal conductivityphonon dispersionelectronic structureflexural mode
0
0 comments X

The pith

Applying x-y strain to monolayer graphene induces curvature that lowers lattice thermal conductivity by changing phonon behavior.

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

The paper examines how controlled x-y strain on a single layer of graphene creates a curved shape that makes the structure more stable energetically. This curvature changes the electronic states so that peaks in the density of states move closer to the energy where conduction begins, and the bands show both flat and straight-line behaviors near that level. Calculations of atomic vibrations reveal that one key low-energy mode switches from a quadratic to a linear pattern, which increases scattering of heat-carrying vibrations and therefore reduces the material's ability to conduct heat. The authors find that the curved sheets stay stable across the range of strains they tested, leading to the conclusion that strain can be used to adjust heat flow by producing different amounts of curvature.

Core claim

Imposing an x-y strain constraint on monolayer graphene produces a curvature that energetically stabilizes the sheet and shifts the flexural acoustic phonon mode from quadratic to linear dispersion, which increases phonon scattering and thereby decreases the lattice thermal conductivity. At the same time the electronic density of states develops Van Hove singularities that approach the Fermi energy, while the highly curved sheets exhibit coexisting flat and linear dispersions near the Fermi level. Phonon dispersion curves confirm dynamic stability for the studied strain and curvature values.

What carries the argument

The x-y strain constraint that forces topological curvature in the graphene sheet, altering both phonon mode dispersions and the electronic density of states.

If this is right

  • Lattice thermal conductivity decreases as curvature increases because of enhanced phonon scattering from the linear flexural mode.
  • Electronic density of states peaks move closer to the Fermi energy with larger curvature.
  • Highly curved sheets develop mixed flat and linear band dispersions near the Fermi level.
  • All examined strain values produce dynamically stable phonon dispersions.
  • The curvature approach offers a route to adjust heat transport while preserving structural stability.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same strain-curvature route might be tested on other two-dimensional sheets such as boron nitride or transition-metal dichalcogenides to control their heat flow.
  • If the stabilization holds in experiment, it could allow free-standing curved graphene shapes without additional supports.
  • Device designs that combine local strain with heat management could exploit the lowered conductivity and shifted electronic states together.
  • Direct comparison of measured versus calculated thermal conductivity on suspended strained samples would test the phonon-scattering mechanism.

Load-bearing premise

The imposed strain constraints create stable curved geometries whose phonon and electronic properties are captured accurately by the computational models.

What would settle it

Experimental measurements of lattice thermal conductivity on graphene sheets subjected to controlled x-y strain that show no reduction, or phonon spectra that lack the predicted linear flexural mode, would falsify the tuning claim.

Figures

Figures reproduced from arXiv: 2604.11367 by D. S. Baptista, E. Lora da Silva, F. J. Manj\'on, M. C. Santos, M. Molinari, Tao Yang, T. Santos, Xidong Lin, Yin Cui.

Figure 1
Figure 1. Figure 1: FIG. 1. The 5 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Potential-energy surface of the S5 graphene system [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Density of States of the 5 [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Electronic band structure of the S1 (left) and S5 (right) graphene systems along the high symmetry [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Isosurfaces of the S1 (left) and S5 (right) enlarged supercells, for better visualization of the charge distribution [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Phonon band dispersion of the five 5 [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. In-plane lattice thermal conductivity of the 5 [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Calculated phonon lifetimes of the S1 (left) and S5 (right) systems at 300 K. The color shades represent the phonon [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Potential energy surfaces of the four graphene systems for different induced-strain. [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. C-C bond length differences of the S5 graphene system. [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Electronic band structure of the S5 graphene system along the high symmetry [PITH_FULL_IMAGE:figures/full_fig_p012_11.png] view at source ↗
read the original abstract

We present a comprehensive set of calculations to investigate the effect of strain-induced x-y topological perturbation in the monolayer graphene sheet. We show that the induced curvature with the defined strain constraint, energetically stabilizes the systems. The electronic properties are modified when the amplitude of the curvature of the sheet increases, which induces Van Hove singularities of the electronic Density of States to approach the Fermi energy. The highly curved system exhibits coexisting flat and linear dispersions close to the Fermi level, which is a promising feature for thermoelectric applications. We also demonstrate, through the phonon dispersion curves, that respective systems are dynamically stable within the studied range of strains/curvatures. Moreover, the flexural acoustic mode transitions from quadratic to linear dispersion under strain, mimicking the 3D behavior and enhancing phonon scattering. The increase of phonon scattering will therefore decrease the value of the lattice thermal conductivity, $\kappa_L$. Such results allows us to conclude that it is possible to tune $\kappa_L$ by applying x-y strain to the monolayer sheet, and inducing different topological curvatures.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper claims to show that strain-induced curvature in monolayer graphene stabilizes the system energetically, modifies the electronic structure by shifting Van Hove singularities in the DOS towards the Fermi energy and creating flat and linear bands, maintains dynamic stability as per phonon dispersions, and alters the flexural acoustic phonon branch from quadratic to linear dispersion. This change is argued to increase phonon scattering, thereby reducing the lattice thermal conductivity κ_L, allowing it to be tuned by x-y strain and different curvatures.

Significance. Should the central claims be substantiated with quantitative calculations, the work would offer a pathway to engineer thermal transport in graphene via curvature, with potential relevance for thermoelectric devices. The electronic structure modifications also suggest possible applications in electronics. However, the current presentation relies on qualitative inferences for the key thermal conductivity result.

major comments (2)
  1. [Abstract] Abstract: The claim that 'the increase of phonon scattering will therefore decrease the value of the lattice thermal conductivity, κ_L' and that 'it is possible to tune κ_L' is not backed by any direct computation of κ_L for the curved systems. Only phonon dispersion curves are referenced, without solving the phonon Boltzmann transport equation, calculating third-order anharmonic force constants, or applying the relaxation time approximation to obtain numerical values of κ_L under different strains.
  2. [Abstract] Abstract: Details on the computational methodology are absent from the abstract and appear insufficiently described. This includes the choice of exchange-correlation functional, plane-wave cutoff, k-point mesh, supercell size for phonon calculations, and convergence criteria, which are essential to assess the reliability of the reported stability, DOS, and dispersion results.
minor comments (2)
  1. [Abstract] The manuscript should clarify the precise definition of the 'x-y strain constraint' used to induce curvature and how the curvature amplitude is quantified.
  2. Figures showing phonon dispersions would benefit from explicit labeling of the flexural acoustic branch and comparison to unstrained graphene to highlight the quadratic-to-linear transition.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment point by point below, acknowledging where revisions are needed to strengthen the presentation.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The claim that 'the increase of phonon scattering will therefore decrease the value of the lattice thermal conductivity, κ_L' and that 'it is possible to tune κ_L' is not backed by any direct computation of κ_L for the curved systems. Only phonon dispersion curves are referenced, without solving the phonon Boltzmann transport equation, calculating third-order anharmonic force constants, or applying the relaxation time approximation to obtain numerical values of κ_L under different strains.

    Authors: We agree that the statement on the reduction and tunability of κ_L is an inference based on the transition of the flexural acoustic phonon branch from quadratic to linear dispersion, which is expected to enhance scattering according to established phonon transport theory in 2D materials. No explicit solution of the phonon Boltzmann transport equation or calculation of third-order anharmonic force constants was performed in this work. We will revise the abstract and discussion sections to present this as a qualitative expectation rather than a computed result, and we will add a statement noting that quantitative confirmation would require additional anharmonic calculations. This revision will be incorporated in the next version of the manuscript. revision: yes

  2. Referee: [Abstract] Abstract: Details on the computational methodology are absent from the abstract and appear insufficiently described. This includes the choice of exchange-correlation functional, plane-wave cutoff, k-point mesh, supercell size for phonon calculations, and convergence criteria, which are essential to assess the reliability of the reported stability, DOS, and dispersion results.

    Authors: Abstracts do not conventionally include detailed methodological parameters. However, we acknowledge that the main text description of the computational setup can be expanded for completeness. In the revised manuscript, we will add explicit details on the exchange-correlation functional, plane-wave energy cutoff, k-point sampling mesh, supercell sizes employed for the phonon calculations, and all convergence criteria for total energy and forces. These additions will be placed in an expanded Methods section to ensure full reproducibility. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation relies on independent computational results

full rationale

The paper reports DFT-based calculations of electronic DOS, band structures, and phonon dispersions for strain-induced curved graphene configurations. The central claim that κ_L can be tuned follows from an inference that the observed transition of the flexural acoustic branch to linear dispersion increases scattering rates. This inference draws on standard phonon transport physics rather than any self-definition, fitted parameter renamed as prediction, or load-bearing self-citation. No equation or step reduces to its own input by construction, and the computational outputs (stability, dispersions) are obtained independently of the final tunability statement.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The abstract does not specify any free parameters, axioms, or invented entities; all details are at a high level with no explicit model assumptions listed.

pith-pipeline@v0.9.0 · 5534 in / 1192 out tokens · 34019 ms · 2026-05-10T15:59:33.316795+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages

  1. [1]

    for the system with largest deformation (S5). By interpolating the curved surface onto a planar ge- ometry, and subsequently following the transformation path toward the system with inverse curvature topology, we observe that the potential well of the S5 system is considerably deep. The energy difference between the FIG. 2. Potential-energy surface of the...

    work page doi:10.54499/- 2022

  2. [2]

    FCT/Mobility/1303210959/2024-25 projects

    and UIDB/04968/2020 (doi.org/10.54499/UIDB/- 04968/2020); and the FCT - Mobility outgoing with Ref. FCT/Mobility/1303210959/2024-25 projects. ELdS further acknowledges the High Performance Computing Chair - a R&D infrastructure (based at the University of ´Evora; PI: M. Avillez), endorsed by Hewlett Packard Enterprise, and involving a consortium of higher...

    work page doi:10.54499/uidb/- 2020

  3. [3]

    A. J. Chaves, T. Frederico, O. Oliveira, W. de Paula, and M. C. Santos, Optical conductivity of curved graphene, J. Phys. Condens. Matter.26, 185301 (2014)

    work page 2014

  4. [4]

    A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, The electronic properties of graphene, Rev. Mod. Phys.81, 109 (2009)

    work page 2009

  5. [5]

    N. Levy, S. A. Burke, K. L. Meaker, M. Panlasigui, A. Zettl, F. Guinea, A. H. C. Neto, and M. F. Crommie, Strain-induced pseudo–magnetic fields greater than 300 tesla in graphene nanobubbles, Science329, 544 (2010)

    work page 2010

  6. [6]

    Lora da Silva, S

    E. Lora da Silva, S. Galmarini, L. Maurizi, M. J. C. San- tos, T. Yang, D. J. Cooke, and M. Molinari, Solid state chemistry: Computational chemical analysis for materi- als science, inComputational Techniques for Analytical Chemistry and Bioanalysis(The Royal Society of Chem- istry, 2020)

    work page 2020

  7. [7]

    Silva, M

    E. Silva, M. Santos, J. Skelton, T. Yang, T. Santos, S. Parker, and A.Walsh, Electronic and phonon insta- bilities in bilayer graphene under applied external bias, Materials Today: Proceedings20, 373 (2020), 11th In- ternational Conference on Advanced Nano Materials

    work page 2020

  8. [8]

    Santos, E

    M. Santos, E. L. da Silva, T. Yang, A. Lopes, and J. Ara´ ujo, Strain-induced effects of topological deformed graphene, J. Magn. Magn. Mater.540, 168429 (2021)

    work page 2021

  9. [9]

    E. L. da Silva, A. Leonardo, T. Yang, M. C. Santos, R. Vilaplana, S. Gallego-Parra, A. Bergara, and F. J. Manj´ on, As2Te3: Pressure-induced three-dimensional Dirac semimetal with ultralow room-pressure lattice thermal conductivity, Phys. Rev. B104, 024103 (2021)

    work page 2021

  10. [10]

    Hussain, H

    W. Hussain, H. Shahzad, P. A. Garc´ ıa-Salaberri, H. Ashraf, and N. A. Shah, Comprehensive review on thermoelectric materials principles, synthesis, nanostruc- turing and advance applications, Int. Commun. Heat Mass Transf.169, 109672 (2025)

    work page 2025

  11. [11]

    N. D. Wood, L. J. Gillie, D. J. Cooke, and M. Moli- nari, A review of key properties of thermoelectric com- posites of polymers and inorganic materials, Materials 15, 10.3390/ma15238672 (2022)

    work page doi:10.3390/ma15238672 2022

  12. [12]

    N. D. Wood, J. S. Tse, J. M. Skelton, D. J. Cooke, L. J. Gillie, S. C. Parker, and M. Molinari, First principles insights into oxide/polymer composites: Sr- tio3/polyaniline/graphene, J. Mater. Sci. Technol.166, 250 (2023)

    work page 2023

  13. [13]

    J. D. Baran, C. Eames, K. Takahashi, M. Molinari, M. S. Islam, and S. C. Parker, Structural, electronic, and trans- port properties of hybrid SrTiO 3-graphene and carbon nanoribbon interfaces, Chem. Mater.29, 7364 (2017)

    work page 2017

  14. [14]

    Y. Cao, V. Fatemi, K. Fang, S. Watanabe, T. Taniguchi, E. Kaxiras, and P. Jarillo-Herrero, Unconventional super- conductivity in magic-angle graphene superlattices, Na- ture556, 43 (2018)

    work page 2018

  15. [15]

    K. Bera, P. Chatterjee, P. Mohan, and A. Saha, Ther- moelectric properties of magic angle twisted bilayer graphene–superconductor heterojunction: Effect of val- ley polarization and trigonal warping, Phys. Rev. B112, 165416 (2025)

    work page 2025

  16. [16]

    pudding mold

    A. J. Thakur, M. Thees, F. Fortuna, E. Frantzeskakis, D. Shiga, H. Kuriyama, M. Nohara, H. Takagi, H. Kumi- gashira, and A. F. Santander-Syro, Large seebeck coeffi- cient driven by “pudding mold” flat band in hole-doped CuRhO2, Phys. Rev. Mater.9, L032401 (2025)

    work page 2025

  17. [17]

    A. A. A. Elameen, D. Dutta, S. Duman, M. Rosmus, G. D’Olimpio, B. G¨ urbulak, D. W. Boukhvalov, A. Agar- wal, and A. Politano, Electronic band structure of gal- lium sulfide (gas) with thickness reduction unveiling parabolic and pudding mold band dispersion, J. Phys. Chem. C.129, 3967 (2025)

    work page 2025

  18. [18]

    G. K. Madsen, J. Carrete, and M. J. Verstraete, Boltz- TraP2, a program for interpolating band structures and calculating semi-classical transport coefficients, Comput. Phys. Commun.231, 140 (2018)

    work page 2018

  19. [19]

    Cheng, Z

    Y. Cheng, Z. Fan, T. Zhang, M. Nomura, S. Volz, G. Zhu, B. Li, and S. Xiong, Magic angle in thermal conductiv- ity of twisted bilayer graphene, Mater. Today Phys.35, 101093 (2023)

    work page 2023

  20. [20]

    Hohenberg and W

    P. Hohenberg and W. Kohn, Inhomogeneous electron gas, Phys. Rev.136, B864 (1964)

    work page 1964

  21. [21]

    Kresse and J

    G. Kresse and J. Furthm¨ uller, Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set, Comput. Mat. Sci.6, 15 (1996)

    work page 1996

  22. [22]

    J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, and K. Burke, Restoring the density-gradient expansion for exchange in solids and surfaces, Phys. Rev. Lett.100, 136406 (2008)

    work page 2008

  23. [23]

    J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, and K. Burke, Erratum: Restoring the density-gradient expansion for exchange in solids and surfaces [Phys. Rev. Lett. 100, 136406 (2008)], Phys. Rev. Lett.102, 039902 (2009)

    work page 2008

  24. [24]

    H. J. Monkhorst and J. D. Pack, Special points for Brillouin-Zone integrations, Phys. Rev. B13, 5188 (1976)

    work page 1976

  25. [25]

    Oba, and I

    A., F. Oba, and I. Tanaka, First-principles calculations of the ferroelastic transition between rutile-type and CaCl2- type SiO 2 at high pressures, Phys. Rev. B78, 134106 (2008)

    work page 2008

  26. [26]

    Chaput, A

    L. Chaput, A. Togo, I. Tanaka, and G. Hug, Phonon- phonon interactions in transition metals, Phys. Rev. B 84, 094302 (2011)

    work page 2011

  27. [27]

    A. Togo, L. Chaput, and I. Tanaka, Distributions of phonon lifetimes in brillouin zones, Phys. Rev. B91, 094306 (2015)

    work page 2015

  28. [28]

    S. Y. Davydov, Effect of van hove singularities on ther- moelectric properties of graphene, Fizika Tverdogo Tela 65, 2024 (2023)

    work page 2024

  29. [29]

    Born, On the stability of crystal lattices

    M. Born, On the stability of crystal lattices. i, Proc. Cam- bridge Philos. Soc.36, 160 (1940)

    work page 1940

  30. [30]

    M. T. Dove,Introduction to Lattice Dynamics(Cam- bridge University Press, 1993)

    work page 1993

  31. [31]

    M. T. Dove,Structure and Dynamics: An Atomic View of Materials(Oxford Master Series in Physics, 2003)

    work page 2003

  32. [32]

    M. T. Dove, Review: Theory of displacive phase transi- tions in minerals, Am. Mineral.82, 213 (2015)

    work page 2015

  33. [33]

    Venkataraman, Soft modes and structural phase tran- sitions, Bull

    G. Venkataraman, Soft modes and structural phase tran- sitions, Bull. Mater. Sci.1, 129 (1979)

    work page 1979

  34. [34]

    Di Gennaro, S

    M. Di Gennaro, S. Saha, and M. Verstraete, Role of dy- namical instability in the ab initio phase diagram of cal- cium, Phys. Rev. Lett.111, 025503 (2013)

    work page 2013

  35. [35]

    Ochoa, E

    H. Ochoa, E. V. Castro, M. Katsnelson, and F. Guinea, Scattering by flexural phonons in suspended graphene under back gate induced strain, Physica E Low Dimens. Syst. Nanostruct.44, 963 (2012)

    work page 2012

  36. [36]

    Amorim and F

    B. Amorim and F. Guinea, Flexural mode of graphene on a substrate, Phys. Rev. B88, 115418 (2013)

    work page 2013

  37. [37]

    Y. W. Sun, W. Liu, I. Hernandez, J. Gonzalez, F. Ro- driguez, D. J. Dunstan, and C. J. Humphreys, 3D strain in 2D materials: To what extent is monolayer graphene graphite?, Phys. Rev. Lett.123, 135501 (2019)

    work page 2019

  38. [38]

    Manunza, R

    L. Manunza, R. Dettori, A. Cappai, and C. Melis, Thermal conductivity of graphene Moir´ e superlattices at small twist angles: An approach-to-equilibrium molec- ular dynamics and Boltzmann transport study, C11, 10.3390/c11030046 (2025)

    work page doi:10.3390/c11030046 2025

  39. [39]

    Han and X

    Z. Han and X. Ruan, Thermal conductivity of monolayer graphene: Convergent and lower than diamond, Phys. Rev. B108, L121412 (2023)

    work page 2023

  40. [40]

    B. D. Kong, S. Paul, M. B. Nardelli, and K. W. Kim, First-principles analysis of lattice thermal conductivity in monolayer and bilayer graphene, Phys. Rev. B80, 033406 (2009)

    work page 2009

  41. [41]

    A. S. Nissimagoudar and N. S. Sankeshwar, Significant reduction of lattice thermal conductivity due to phonon confinement in graphene nanoribbons, Phys. Rev. B89, 235422 (2014)

    work page 2014