Finite-frequency anomaly-induced electromechanical response of Dirac fermions in deformed graphene
Pith reviewed 2026-05-19 22:18 UTC · model grok-4.3
The pith
Mechanical deformations in graphene generate transverse electric currents via the parity-odd response of massive Dirac cones.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The coefficient of this mixed electromechanical response is the parity-odd current-current correlator of a massive Dirac cone. For an insulating cone the coefficient is the one-cone Chern-Simons value, while for a doped cone in the local regime it is reduced by the Berry curvature factor m/|μ|.
What carries the argument
Emergent phonon gauge field that couples to the Dirac current identically to the electromagnetic vector potential, whose strength is fixed by the parity-odd current-current correlator.
If this is right
- A traveling flexural wave generates a transverse second-harmonic current.
- A static ripple mixed with a dynamic phonon generates a transverse current at the drive frequency.
- Two non-collinear modes generate charge modulation through the emergent phonon flux.
- For sublattice-gapped graphene with valley-odd deformation gauge coupling the two valleys add rather than cancel.
- The response exhibits selection rules in direction, phase, frequency and gate-voltage dependence.
Where Pith is reading between the lines
- The same anomaly-driven channel could be realized in other two-dimensional Dirac systems whose deformations induce gauge-like couplings.
- Devices could convert mechanical waves directly into electrical signals at harmonic frequencies without external magnetic fields.
- Measurements at higher drive frequencies or with controlled intervalley scattering would test the range of validity of the low-energy approximation.
- Varying deformation amplitude independently of frequency could isolate the nonlinear contributions to the mixed response.
Load-bearing premise
The low-energy Dirac theory remains valid and the deformation-induced geometric vertices can be represented as an emergent phonon gauge field that couples to the Dirac current identically to the electromagnetic vector potential.
What would settle it
Detection of a transverse electrical signal at twice the frequency of an applied flexural wave whose amplitude is constant at small gate voltage then decays as one over chemical potential magnitude, with phase fixed by the sign of any sublattice gap.
Figures
read the original abstract
A deformation of a graphene sheet changes more than the positions of the atoms. In the low-energy Dirac theory it also produces geometric electron-phonon vertices. One of these vertices acts as an emergent phonon gauge field, $\calA_\mu$, which couples to the same Dirac current as the electromagnetic vector potential. This shared current vertex gives a direct route from mechanics to electronics: a moving deformation can generate a transverse electric current, and a deformation pattern with emergent phonon flux can bind electric charge. We show that the coefficient of this mixed electromechanical response is the parity-odd current-current correlator of a massive Dirac cone. For an insulating cone the coefficient is the one-cone Chern-Simons value, while for a doped cone in the local regime it is reduced by the Berry curvature factor $m/|\mu|$. We apply the response to explicit deformations. A traveling flexural wave generates a transverse second-harmonic current; a static ripple mixed with a dynamic phonon generates a transverse current at the drive frequency; and two non-collinear modes can generate charge modulation through the emergent phonon flux. We keep the spin and valley sum explicit, so the paper shows when the one-cone anomaly becomes a charge current in graphene and when it instead appears in a valley, spin, or spin-valley channel. For sublattice-gapped graphene with a valley-odd deformation gauge coupling, the two valleys add rather than cancel. The experimentally sharp signature is a transverse electrical signal at twice the flexural-wave frequency, with a phase fixed by the sign of the sublattice gap and a gate dependence that crosses over from a gap plateau to a $1/|\mu|$ decay. These direction, phase, frequency, and gate-voltage selection rules give clean tests of the anomaly-induced electromechanical channel in deformed graphene.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a low-energy Dirac theory of electromechanical response in deformed graphene. It identifies a geometric electron-phonon vertex that functions as an emergent phonon gauge field calA_μ coupling to the Dirac current identically to the electromagnetic vector potential. The coefficient of the resulting mixed response is equated to the parity-odd current-current correlator of a massive Dirac cone; this reduces to the one-cone Chern-Simons value for an insulating cone and to a Berry-curvature-suppressed factor m/|μ| for a doped cone in the local regime. Explicit applications are given to traveling flexural waves (producing transverse second-harmonic current), static ripples mixed with dynamic phonons, and non-collinear modes generating charge modulation via emergent phonon flux. Spin and valley sums are kept explicit, and experimental signatures (frequency doubling, phase fixed by sublattice-gap sign, gate-voltage crossover from plateau to 1/|μ| decay) are highlighted.
Significance. If the central identification is correct, the work supplies a direct, anomaly-based route from mechanical deformation to transverse electronic current and charge binding in graphene, extending known Chern-Simons and Berry-curvature responses to finite-frequency electromechanical settings. The explicit valley and spin accounting clarifies when the one-cone anomaly appears as a charge current versus a valley or spin channel, and the proposed selection rules (direction, phase, frequency, gate dependence) constitute clean, falsifiable tests. The approach is parameter-free once the Dirac mass and chemical potential are fixed, and the applications to concrete phonon modes are concrete.
major comments (2)
- [§2.2] §2.2 and the paragraph following Eq. (3): the claim that the deformation-induced geometric vertex can be rewritten exactly as minimal coupling to an emergent gauge field calA_μ with no additional scalar potentials or non-minimal terms is load-bearing for the parity-odd response. The continuum limit from the honeycomb lattice necessarily truncates higher-gradient and intervalley contributions; the manuscript must show explicitly that these corrections do not alter the transverse current at the order kept in the finite-frequency calculation.
- [§3] §3, Eq. (8): the reduction of the response coefficient by the factor m/|μ| for a doped cone in the local regime is obtained from the parity-odd correlator. The derivation assumes the local (q→0) limit remains valid at finite frequency and finite doping; an explicit check that finite-q or finite-ω corrections do not reintroduce additional terms that would change the quoted factor is required to support the doped-case claim.
minor comments (2)
- The symbol calA_μ is introduced without an immediate comparison table to the electromagnetic A_μ; a short side-by-side definition of the two couplings would remove potential confusion for readers.
- Figure 2 caption states 'transverse current at 2ω' but does not label the phase reference; adding the explicit phase relation to the sign of the sublattice gap would make the experimental signature easier to read.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. The two major comments raise valid points about the robustness of the continuum approximation and the local-limit derivation. We address each below and will incorporate the requested clarifications and checks into the revised version.
read point-by-point responses
-
Referee: [§2.2] §2.2 and the paragraph following Eq. (3): the claim that the deformation-induced geometric vertex can be rewritten exactly as minimal coupling to an emergent gauge field calA_μ with no additional scalar potentials or non-minimal terms is load-bearing for the parity-odd response. The continuum limit from the honeycomb lattice necessarily truncates higher-gradient and intervalley contributions; the manuscript must show explicitly that these corrections do not alter the transverse current at the order kept in the finite-frequency calculation.
Authors: We agree that an explicit demonstration is needed. The geometric vertex is obtained from the leading long-wavelength expansion of the honeycomb tight-binding model, and the parity-odd response is tied to the chiral anomaly, which is infrared-protected. In the revised manuscript we add a short paragraph at the end of §2.2 together with a new appendix. There we retain the next-order gradient terms and the leading intervalley matrix elements, compute their contribution to the current response, and show that they enter only the parity-even sector or appear multiplied by extra powers of (a q), where a is the lattice constant. Consequently they do not modify the transverse current at the order kept in the finite-frequency calculation. This justifies the minimal-coupling identification used for the anomaly-induced electromechanical response. revision: yes
-
Referee: [§3] §3, Eq. (8): the reduction of the response coefficient by the factor m/|μ| for a doped cone in the local regime is obtained from the parity-odd correlator. The derivation assumes the local (q→0) limit remains valid at finite frequency and finite doping; an explicit check that finite-q or finite-ω corrections do not reintroduce additional terms that would change the quoted factor is required to support the doped-case claim.
Authors: The factor m/|μ| is obtained from the exact analytic expression for the parity-odd current-current correlator of a massive Dirac cone at finite chemical potential in the q → 0 limit. To address the concern about finite-q and finite-ω corrections, the revised manuscript includes a supplementary calculation of the full momentum- and frequency-dependent bubble diagram. The leading term remains m/|μ|; relative corrections are O((q v_F/μ)^2, (ω/μ)^2). For the phonon wave-vectors and frequencies relevant to the applications discussed in the paper these corrections remain small. We will add this explicit check to §3 together with a brief discussion of the regime of validity. revision: yes
Circularity Check
No significant circularity; central identification follows from explicit low-energy mapping and independent correlator evaluation
full rationale
The derivation begins from the standard continuum limit of deformed graphene, where geometric vertices are expanded to yield an emergent phonon gauge field that couples minimally to the Dirac current (identical to the EM potential by construction of the low-energy theory). The mixed electromechanical coefficient is then identified with the parity-odd current-current correlator, whose value is obtained by direct computation for the massive Dirac cone rather than by fitting or redefinition. The insulating-case result recovers the standard one-cone Chern-Simons term, while the doped local-regime reduction by m/|μ| follows from the Berry curvature factor in the same correlator. No self-definitional loops, fitted inputs renamed as predictions, or load-bearing self-citations appear; the mapping and correlator evaluation remain independent of the final response formula. The paper keeps valley and spin sums explicit, allowing falsifiable selection rules without tautological reduction.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Low-energy effective Dirac theory remains valid for deformed graphene
- domain assumption Deformation produces geometric electron-phonon vertices that can be written as an emergent gauge field A_μ
invented entities (1)
-
emergent phonon gauge field A_μ
no independent evidence
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the coefficient of this mixed electromechanical response is the parity-odd current-current correlator of a massive Dirac cone... one-cone Chern-Simons value... reduced by the Berry curvature factor m/|μ|
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
A_μ which couples to the same Dirac current as the electromagnetic vector potential
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
The factor q3h2 0 is the product of a slopeqh 0 and a curvatureq 2h0. Thus the effect is strongest for fast, short-wavelength rip- ples while remaining within the small-slope regime. For the realistic sublattice-gapped, valley-odd-coupling case derived in Sec. VIII, this local signal is a true charge cur- rent rather than only a valley current. This makes...
-
[2]
Dy- namics of Dirac and Weyl fermions on a two-dimensional surface,
A. R. Kavalov, I. K. Kostov, and A. G. Sedrakyan, “Dy- namics of Dirac and Weyl fermions on a two-dimensional surface,” Phys. Lett. B175, 331 (1986)
work page 1986
-
[3]
Dirac and Weyl fermions coupled to two-dimensional surfaces: Determinants,
A. G. Sedrakyan and R. Stora, “Dirac and Weyl fermions coupled to two-dimensional surfaces: Determinants,” Phys. Lett. B188, 442 (1987)
work page 1987
-
[4]
Optical con- ductivity of graphene in the presence of random lattice deformations,
A. Sinner, A. Sedrakyan, and K. Ziegler, “Optical con- ductivity of graphene in the presence of random lattice deformations,” Phys. Rev. B83, 155115 (2011)
work page 2011
-
[5]
Deformation of a graphene sheet: Interaction of fermions with phonons,
A. Sedrakyan, A. Sinner, and K. Ziegler, “Deformation of a graphene sheet: Interaction of fermions with phonons,” Phys. Rev. B103, L201104 (2021)
work page 2021
-
[6]
M. A. H. Vozmediano, M. I. Katsnelson, and F. Guinea, “Gauge fields in graphene,” Phys. Rep.496, 109 (2010)
work page 2010
-
[7]
Novel effects of strains in graphene and other two dimensional materials,
B. Amorim, A. Cortijo, F. de Juan, A. G. Grushin, F. Guinea, A. Gutierrez-Rubio, H. Ochoa, V. Parente, R. Roldan, P. San-Jose, J. Schiefele, M. Sturla, and M. A. H. Vozmediano, “Novel effects of strains in graphene and other two dimensional materials,” Phys. Rep.617, 1 (2016)
work page 2016
-
[8]
Space dependent Fermi velocity in strained graphene,
F. de Juan, M. Sturla, and M. A. H. Vozmediano, “Space dependent Fermi velocity in strained graphene,” Phys. Rev. Lett.108, 227205 (2012)
work page 2012
-
[9]
Current-current correlation function in 3D massive Dirac theory with chemical potential,
E. Apresyan, Sh. Khachatryan, and A. Sedrakyan, “Current-current correlation function in 3D massive Dirac theory with chemical potential,” Mod. Phys. Lett. A30, 1550035 (2015)
work page 2015
-
[10]
E. Apresyan, and A. Sedrakyan, ”Transport properties of fermions with moat spectra”, Modern Physics Letters A34, 1950041 (2019)
work page 2019
-
[11]
Parity violation and gauge noninvariance of the effective gauge field action in three dimensions,
A. N. Redlich, “Parity violation and gauge noninvariance of the effective gauge field action in three dimensions,” Phys. Rev. Lett.52, 18 (1984)
work page 1984
-
[12]
Gauge noninvariance and parity noncon- servation of three-dimensional fermions,
A. N. Redlich, “Gauge noninvariance and parity noncon- servation of three-dimensional fermions,” Phys. Rev. D 29, 2366 (1984)
work page 1984
-
[13]
Condensed-matter simulation of a three-dimensional anomaly,
G. W. Semenoff, “Condensed-matter simulation of a three-dimensional anomaly,” Phys. Rev. Lett.53, 2449 (1984)
work page 1984
-
[14]
Fractional charge and zero modes for pla- nar systems in a magnetic field,
R. Jackiw, “Fractional charge and zero modes for pla- nar systems in a magnetic field,” Phys. Rev. D29, 2375 (1984)
work page 1984
- [15]
-
[16]
F. D. M. Haldane, “Model for a quantum Hall effect with- out Landau levels: Condensed-matter realization of the parity anomaly,” Phys. Rev. Lett.61, 2015 (1988)
work page 2015
-
[17]
Emergent Chern-Simons exci- tations due to electron-phonon interaction,
A. Sinner and K. Ziegler, “Emergent Chern-Simons exci- tations due to electron-phonon interaction,” Phys. Rev. B93, 125112 (2016)
work page 2016
-
[18]
Spontaneous mass generation due to phonons in a two-dimensional Dirac fermion sys- tem,
A. Sinner and K. Ziegler, “Spontaneous mass generation due to phonons in a two-dimensional Dirac fermion sys- tem,” Ann. Phys.400, 262 (2019)
work page 2019
-
[19]
Interplay of Coulomb and electron-phonon interactions in graphene,
D. M. Basko and I. L. Aleiner, “Interplay of Coulomb and electron-phonon interactions in graphene,” Phys. Rev. B 77, 041409(R) (2008)
work page 2008
-
[20]
The electronic properties 14 of graphene,
A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, “The electronic properties 14 of graphene,” Rev. Mod. Phys.81, 109 (2009)
work page 2009
-
[21]
Phonons and electron-phonon scattering in carbon nanotubes,
H. Suzuura and T. Ando, “Phonons and electron-phonon scattering in carbon nanotubes,” Phys. Rev. B65, 235412 (2002)
work page 2002
-
[22]
Symmetry-based approach to electron- phonon interactions in graphene,
J. L. Ma˜ nes, “Symmetry-based approach to electron- phonon interactions in graphene,” Phys. Rev. B76, 045430 (2007)
work page 2007
-
[23]
Energy gaps and a zero-field quantum Hall effect in graphene by strain engineering,
F. Guinea, M. I. Katsnelson, and A. K. Geim, “Energy gaps and a zero-field quantum Hall effect in graphene by strain engineering,” Nat. Phys.6, 30 (2010)
work page 2010
-
[24]
Generating quantizing pseudomagnetic fields by bending graphene ribbons,
F. Guinea, A. K. Geim, M. I. Katsnelson, and K. S. Novoselov, “Generating quantizing pseudomagnetic fields by bending graphene ribbons,” Phys. Rev. B81, 035408 (2010)
work page 2010
-
[25]
Strain-induced pseudomagnetic field for novel graphene electronics,
T. Low and F. Guinea, “Strain-induced pseudomagnetic field for novel graphene electronics,” Nano Lett.10, 3551 (2010)
work page 2010
-
[26]
Gauge fields from strain in graphene,
F. de Juan, J. L. Ma˜ nes, and M. A. H. Vozmediano, “Gauge fields from strain in graphene,” Phys. Rev. B 87, 165131 (2013)
work page 2013
-
[27]
Pseudo- magnetic field distribution and pseudo-Landau levels in suspended graphene flakes,
M. Mucha-Kruczy´ nski and V. I. Fal’ko, “Pseudo- magnetic field distribution and pseudo-Landau levels in suspended graphene flakes,” Solid State Commun.152, 1442 (2012)
work page 2012
-
[28]
Transport signatures of pseudomag- netic Landau levels in strained graphene ribbons,
D. A. Gradinar, M. Mucha-Kruczy´ nski, H. Schomerus, and V. I. Fal’ko, “Transport signatures of pseudomag- netic Landau levels in strained graphene ribbons,” Phys. Rev. Lett.110, 266801 (2013)
work page 2013
-
[29]
Strain-induced pseudo-magnetic fields greater than 300 tesla in graphene nanobubbles,
N. Levy, S. A. Burke, K. L. Meaker, M. Panlasigui, A. Zettl, F. Guinea, A. H. Castro Neto, and M. F. Crommie, “Strain-induced pseudo-magnetic fields greater than 300 tesla in graphene nanobubbles,” Science329, 544 (2010)
work page 2010
-
[30]
Visualizing strain-induced pseudomagnetic fields in graphene through an hBN mag- nifying glass,
Y. Jiang, J. Mao, J. Duan, X. Lai, K. Watanabe, T. Taniguchi, and E. Y. Andrei, “Visualizing strain-induced pseudomagnetic fields in graphene through an hBN mag- nifying glass,” Nano Lett.17, 2839 (2017)
work page 2017
-
[31]
Pseudo-magnetic field-induced slow carrier dynamics in periodically strained graphene,
D.-H. Kang, J. W. Park, S. J. Ha, G. Y. Choi, S. Lee, H. H. Park, C. J. Lee, K. Watanabe, T. Taniguchi, J. H. Kim, and H. Cheong, “Pseudo-magnetic field-induced slow carrier dynamics in periodically strained graphene,” Nat. Commun.12, 5087 (2021)
work page 2021
-
[32]
Protected fermionic zero modes in periodic gauge fields,
V. T. Phong and E. J. Mele, “Protected fermionic zero modes in periodic gauge fields,” Phys. Rev. B111, 125129 (2025)
work page 2025
-
[33]
Elastic screening of pseudogauge fields in graphene,
C. De Beule, R. Smeyers, W. N. Luna, E. J. Mele, and L. Covaci, “Elastic screening of pseudogauge fields in graphene,” Phys. Rev. Lett.134, 046404 (2025)
work page 2025
-
[34]
S. Heidari, S. Parsi, and P. Ghaemi, “Effects of strain- induced pseudogauge fields on exciton dispersion, trans- port, and interactions in transition metal dichalcogenide nanoribbons,” Phys. Rev. B112, 155149 (2025)
work page 2025
-
[35]
We useℏ= 1. Greek indicesµ, ν, ρ= 0,1,2 denote time and the two in-plane spatial directions; bold symbols such asqdenote two-dimensional spatial vectors. The tensor convention isϵ 012 = +1. The symbole >0 denotes the magnitude of the electron charge. The scalarµdenotes the chemical potential and should not be confused with the Greek spacetime index. These...
-
[36]
T. A. Sedrakyan, E. G. Mishchenko, and M. E. Raikh, “Smearing of the two-dimensional Kohn anomaly in a nonquantizing magnetic field: Implications for interac- tion effects,” Phys. Rev. Lett.99, 036401 (2007)
work page 2007
-
[37]
T. A. Sedrakyan and M. E. Raikh, ”Crossover from Weak Localization to Shubnikov–de Haas Oscillations in a High-Mobility 2D Electron Gas,” Phys. Rev. Lett.100, 106806 (2008)
work page 2008
-
[38]
T. A. Sedrakyan and M. E. Raikh, ”Magneto-Oscillations due to Electron-Electron Interactions in the ac Conduc- tivity of a Two-Dimensional Electron Gas,” Phys. Rev. Lett.100, 086808 (2008)
work page 2008
-
[39]
Persistent Friedel oscillations in graphene due to a weak magnetic field,
K. Wang, M. E. Raikh, and T. A. Sedrakyan, “Persistent Friedel oscillations in graphene due to a weak magnetic field,” Phys. Rev. B103, 085418 (2021)
work page 2021
-
[40]
Interac- tion effects in graphene in a weak magnetic field,
K. Wang, M. E. Raikh, and T. A. Sedrakyan, “Interac- tion effects in graphene in a weak magnetic field,” Phys. Rev. B104, L161102 (2021)
work page 2021
-
[41]
Ballistic magnetotrans- port in graphene,
K. Wang and T. A. Sedrakyan, “Ballistic magnetotrans- port in graphene,” Phys. Rev. B105, L121114 (2022)
work page 2022
-
[42]
A. A. Chaika, Y. Kulynych, D. O. Oriekhov, and S. G. Sharapov, “Density of states and differential entropy in Dirac materials in crossed magnetic and in-plane electric fields,” Phys. Rev. B111, 085426 (2025)
work page 2025
-
[43]
Quantization and quantum oscillations of the sublattice charge order in Dirac insulators,
A. Tarafdar and T. A. Sedrakyan, “Quantization and quantum oscillations of the sublattice charge order in Dirac insulators,” Phys. Rev. B112, 245106 (2025)
work page 2025
-
[44]
Tunable topological valley plasmons in graphene/hexagonal boron nitride moire superlattices,
Y. Du, L. Sun, C. Ding, Z. Zhang, H. Sun, and M. Zhao, “Tunable topological valley plasmons in graphene/hexagonal boron nitride moire superlattices,” Phys. Rev. B112, 155433 (2025)
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.