{"id":"57c55f74-4192-4ce2-aecc-1ab7b4be5ded","arxiv_id":"2505.21056","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A scaling transformation for strained graphene preserves pseudogauge fields and enables quantum transport simulations showing valley-polarized pseudomagnetic focusing and snake-state oscillations.","lead":"This paper develops a scaled-down version of the graphene lattice that keeps the physics of strain-induced pseudomagnetic fields intact, making simulations of micrometer-sized devices possible. It predicts that bent and S-shaped graphene ribbons can focus valley-polarized currents and produce snake-state conductance oscillations with clear experimental fingerprints.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Snake-state transport is simulated at s=2 with Bs up to ~6 T, where the scaled strain is comparable to the s=4, Bs≈3 T regime already shown to break the scaling law; no s=1 benchmark is provided for this geometry.","rationale":"The central method is largely sound: for small scaled strains the derivation in SM S2 is consistent, the s=1 versus s=2 LDoS benchmarks match the analytic Landau levels, and the two-terminal conductance plateaus in SM S5 support the in-plane scaling for moderate fields. The reader's conditional verdict is appropriate. My stress test sharpens the reader's weakest assumption by translating the known s=4 failure into the dimensionless error parameter s·(strain) and applying it to the snake-state device. The parameters of Fig. 4 (s=2, Bs up to ~6 T, ribbon length ~350 nm) push the scaled strain close to or beyond the point where the linear bond expansion breaks down, and no direct unscaled benchmark exists for that geometry. Because the snake-state conductance oscillations are a headline prediction, the absence of such a check is the most load-bearing unresolved issue. I therefore keep the verdict CONDITIONAL: the paper should either provide an s=1 (or reduced-size) benchmark for the S-shaped device at the relevant Bs scale, or restrict the quantitative claims to the regime where the scaling is validated.","tokens_in":20448,"tokens_out":39687,"duration_ms":428849,"concrete_test":"Recompute the two-terminal conductance of the S-shaped ribbon at s=1 for a reduced-size copy (e.g., half the linear dimensions, with u0 adjusted so the physical PMF profile and Bs range of Fig. 4 are preserved) and overlay the G(u0,E) oscillations on the s=2 result. If the conductance peak positions in Bs shift by more than ~10% anywhere in the plotted range, the Fig. 4 snake-state predictions are not quantitatively reliable and the paper should either restrict the claims or provide a corrected scaling factor for the high-strain regime.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The scaling transformation (4) preserves the pseudomagnetic field only at linear order in the displacement gradients across a scaled bond, Eq. (S21). The relevant error is controlled by the product s × (local strain), not by s or Bs alone. The paper's own SM S5 benchmark shows that for s=4 deviations from the unscaled conductance steps set in for Bs ≳ 3 T, i.e. at s·Bs ≈ 12 T. The snake-state device in Fig. 4 uses s=2 and Bs up to ~6 T, giving the same s·Bs ≈ 12, and the local strains in the bent-ribbon displacement field (S39) grow as u0·L, reaching the order of 10% at the largest Bs shown; after scaling by s=2 the bond-level strain is roughly twice that. In this regime the linearized expansion u(r+s d0)−u(r)≈s(d0·∇)u, on which Eq. (4) is based, is not controlled, so the PMF profile actually implemented for the S-shaped ribbon can deviate from the intended uniform/tanh profile. The conductance oscillations in Fig. 4 are central predictions of the paper, yet no s=1 vs s=2 comparison is reported for this geometry or Bs range; the s=1 transport benchmarks in SM S5 use a smaller bent ribbon at lower fields. The high-Bs snake-state peaks could therefore be shifted or distorted by the scaling error rather than representing the unscaled physical device.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the scalable tight-binding model for graphene to strained systems. The proposed scaling transformation (Eq. 4) changes the lattice spacing, hopping amplitude, in-plane displacement field, and out-of-plane displacement field such that the pseudogauge field is approximately unchanged. The authors benchmark the transformation by comparing pseudo Landau levels for scaling factors s = 1 to 4, including coexisting real and pseudomagnetic fields. They then apply the method to quantum transport in two mesoscopic devices: a bent zigzag ribbon with transverse pseudomagnetic focusing and valley-polarized current (Fig. 3), and an S-shaped ribbon with a sign-changing pseudomagnetic field giving rise to pseudomagnetic snake states and conductance oscillations (Fig. 4). The paper argues that the method enables quantum transport simulations of micrometer-scale strained graphene devices with realistic geometries.","tokens_in":20746,"tokens_out":38321,"duration_ms":399106,"significance":"If the scaling method is fully validated, it is a practically valuable extension of the earlier scalable tight-binding model: it makes quantum transport calculations of micron-scale strained graphene devices computationally feasible and produces concrete, falsifiable transport signatures such as focusing peaks and snake-state oscillations. The analytic derivation of the scaling law and the s = 1 to 4 Landau-level benchmarks are genuine strengths, as is the authors' explicit acknowledgment of deviations for large displacements (Fig. 1g) and at s = 4 for Bs above about 3 T (SM S5). The principal concern is that the quantitative transport predictions, especially the high-field snake-state results, are made in regimes where this scaling validity has not been demonstrated.","major_comments":[{"comment":"The snake-state simulations in Fig. 4 are performed at s = 2 with Bs up to approximately 6 T. The scaling benchmark in SM S5 shows that for s = 4 the two-terminal conductance plateaus begin to deviate from the unscaled positions for Bs greater than about 3 T. Because the error in the linearized bond-vector expansion of Eq. (S21) grows with the product of the scaling factor and the local strain, the (s = 2, Bs approximately 6 T) regime used in Fig. 4 is comparable to the (s = 4, Bs approximately 3 T) regime in which the benchmark already fails. No s = 1 comparison is reported for the S-shaped geometry or for this Bs range; the s = 1 transport results in SM S5 use a smaller, uniformly bent ribbon at lower fields. The conductance oscillations in Fig. 4 are a central prediction, so the high-Bs peaks could be shifted or distorted by the scaling approximation rather than representing the physical device. The authors should either provide a scaling-convergence check for this geometry (for example, an s = 1 calculation on a smaller equivalent S-shaped device, or a direct comparison of the implemented PMF profile against the target profile over the full Bs range) or restrict the quantitative claims to the regime where the scaling has been validated.","section":"§4, Fig. 4; SM S5"}],"minor_comments":[{"comment":"The 1.1 conversion factor used to convert the strain parameter u0 to the pseudomagnetic field Bs is mentioned only in the Supplemental Material; since it sets the horizontal axes of Figs. 3 and 4, it should be stated in the main text and its sensitivity should be discussed.","section":"SM S5"},{"comment":"The abstract states that the scaling is valid as long as the atomic displacements vary slowly with respect to the scaled lattice, but the benchmarks in Fig. 1(g) and SM S5 show that large local strain also breaks the scaling; the stated validity condition should include a small-strain requirement.","section":"Abstract"},{"comment":"In the inset of Fig. 2(b) the plotted quantity is D/s^2, while the text says that the peak height scales with the area as s^2; a sentence clarifying the normalization would prevent confusion.","section":"Fig. 2(b)"},{"comment":"The conclusion mentions realistic device dimensions up to one micron, but the two main transport simulations use ribbons of width 600 nm and 206 nm; the micron-scale statement is supported only by the triaxial flake in SM S8 (Dh = 1500 nm), so the wording should be adjusted or the SM result should be cited.","section":"Conclusions"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about the snake-state simulations is legitimate and should be addressed before publication. The core scaling method appears sound in the benchmarked regimes, and the missing s = 1 validation for the S-shaped ribbon is fixable with additional calculations or a more cautious presentation. I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the scaling transformation for strained graphene, Eq. (4): scale the lattice spacing by s, the hopping by 1/s, in-plane displacements by s, and out-of-plane displacements by sqrt(s). I checked the derivation in SM S2; it follows from the linear expansion of the hopping and correctly rescales the pseudogauge field. The pseudo-Landau-level benchmarks for s=1..4 are credible, and the transport maps for transverse pseudomagnetic focusing and snake states are the first mesoscopic simulations of these effects with realistic geometry. That is a real contribution.\n\nThe soft spot is exactly where the stress-test note lands. The scaling law is controlled only at linear order in the displacement gradient across a scaled bond, Eq. (S21). The paper itself flags this: SM S5 shows s=4 conductance plateaus deviate for Bs ≳ 3 T, and Fig. 1(g) shows PMF profile mismatch for large displacements. But Fig. 4 pushes the snake-state device to s=2 and Bs up to about 6 T, which is a comparable s*Bs product (about 12 T) to the s=4, Bs=3 T failure point, and the bent-ribbon displacement field (S39) produces local strains that grow with device width. The snake device is roughly four times wider than the SM S5 benchmark ribbon, so the scaled strain there is likely larger than in the benchmark. The authors provide no s=1 versus s=2 comparison for the snake geometry, so the positions and strengths of the high-Bs conductance peaks in Fig. 4 could be shifted by the scaling error rather than reflecting the unscaled device. That needs to be addressed, either by a benchmark at s=1 for this device or by restricting the snake-state claims to the regime where the scaling is verified.\n\nAlso minor: the 1.1 calibration factor fitted to the s=1 pseudo Landau levels (SM S5) is a single free parameter. It does not make the method circular, but the transport maps in Figs. 3 and 4 inherit that fitted conversion, so the absolute Bs axis carries about 10% uncertainty. The neglect of the deformation potential is stated and justified by screening arguments, which is fine, though the qualitative valley-polarization conclusions are gauge-field-only results.\n\nWho should read this: condensed-matter theorists working on straintronics, mesoscopic transport, or valleytronics. The paper deserves a serious referee; the method is useful and the benchmarks are honest. Send it to review, and ask the authors to benchmark the snake-state geometry at s=1.","headline":"A genuinely useful scaling method for strained-graphene transport, with the snake-state results pushed into a regime where the scaling error is not yet benchmarked.","tokens_in":21311,"tokens_out":2266,"would_cite":true,"duration_ms":23262,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["72.80.Vp","73.43.-f"],"model":"deepseek-v4-flash","headline":"A scaling transformation for strained graphene preserves the pseudomagnetic field and low-energy Dirac physics, making micrometer-scale device transport simulations feasible.","keywords":["strained graphene","pseudomagnetic field","scalable tight-binding model","transverse pseudomagnetic focusing","valley polarization","snake states","Landau levels","straintronics"],"falsifier":"A direct check is to compare the local density of states or two-terminal conductance of a scaled strained ribbon at $s=4$ against $s=1$ for a pseudomagnetic field above 3 T; the authors report in the Supplemental Material that the conductance plateaus deviate from the predicted positions in exactly this regime. Experimentally, measuring the transverse-focusing peak positions in a bent graphene ribbon under strong bending would test whether the scaled model's predictions hold at large strain.","tokens_in":20253,"feed_emoji":"🧲","tokens_out":9918,"duration_ms":86526,"temperature":0.7,"pith_summary":"This paper argues that a simple scaling transformation makes quantum transport simulations of strained graphene devices tractable at micrometer scale. The transformation enlarges the lattice spacing by a factor $s$, reduces the hopping amplitude to $t_0/s$, and scales the in-plane displacement field by $s$ and the out-of-plane one by $\\sqrt{s}$; the pseudomagnetic field $B_s = \\nabla \\times \\mathbf{A}_s$ is then unchanged as long as the scaled strain remains small. The authors verify the invariance by matching pseudo Landau levels computed at scaling factors $s=1,2,3,4$. Using the method, they predict transverse pseudomagnetic focusing in a bent ribbon, generating a valley-polarized current with characteristic conductance oscillations, and pseudomagnetic snake states in an S-shaped ribbon, also marked by conductance oscillations. If correct, the approach gives experimenters a computational tool for designing graphene straintronics devices.","feed_headline":"Scaling law leaves graphene's pseudomagnetic fields intact","feed_subtitle":"A lattice-coarsening trick preserves strain-induced magnetic fields, making realistic graphene devices computable.","key_machinery":"The load-bearing object is the scaling transformation of Eq. (4): $a_0 \\to s a_0$, $t_0 \\to t_0/s$, $\\mathbf{u} \\to s\\mathbf{u}$, $h \\to \\sqrt{s}h$. It works because the pseudogauge field in lowest order depends on the strain tensor $u_{ij} = [\\partial_i u_j + \\partial_j u_i + (\\partial_i h)(\\partial_j h)]/2$, which scales linearly; the displacement-scaling part compensates for the fact that bond-vector differences are sampled over the scaled lattice vectors, $\\mathbf{u}(r+s d^0_n)-\\mathbf{u}(r) \\approx s (d^0_n \\cdot \\nabla)\\mathbf{u}(r)$. The pseudogauge-field relation of Eq. (3) connects the hopping modulations to $\\mathbf{A}_s$ and is what the scaling must preserve.","core_discovery":"The central claim is that strain-induced pseudogauge fields survive a coarse-graining of the graphene lattice. Starting from the tight-binding Hamiltonian with hopping amplitude $t_{ij}=t_0\\exp[-\\beta(|r_i-r_j|/a_0-1)]$, the strain enters through a modulation $\\delta t_n = -t_0\\beta \\sum_{ij} u_{ij} d^0_{n,i}d^0_{n,j}/a_0^2$, which produces the pseudogauge field via $ev_F e^{i\\tau\\theta}[A_{s,x}-i\\tau A_{s,y}] = -\\sum_n \\delta t_n e^{i\\tau K \\cdot d^0_n}$. Under the scaling $a_0 \\to s a_0$, $t_0 \\to t_0/s$, $\\mathbf{u} \\to s\\mathbf{u}$, $h \\to \\sqrt{s} h$, the strain tensor $u_{ij} \\to s u_{ij}$, so $B_s = \\nabla \\times \\mathbf{A}_s$ is unchanged and the low-energy Dirac physics is preserved. This is confirmed numerically by comparing pseudo Landau levels of a triaxially strained flake at $s=1,2,3,4$. Applied to a bent zigzag ribbon, the model yields transverse pseudomagnetic focusing with valley-polarized current; applied to an S-shaped ribbon, it yields pseudomagnetic snake states localized at the sign change of $B_s$, both producing experimentally visible conductance oscillations.","pith_inferences":["The same scaling could be applied to other Dirac materials with strain-tunable gauge fields, or to moiré systems where lattice relaxation creates pseudogauge fields, potentially enabling transport simulations of twisted multilayers at realistic moiré periods.","Because the scaling fails for large scaled strain, the practical sweet spot is moderate pseudomagnetic fields (below roughly 3 T at $s=4$); pushing to the hundreds-of-tesla regime of nanobubbles would need higher-order bond-length corrections.","The predicted valley-polarized focusing could be tested as a strain-only valley filter: a bent ribbon with two closely spaced contacts should show a nonlocal conductance signal that reverses when the bending direction is flipped, without any external magnetic field.","A natural numerical follow-up is to compute the full shot-noise or Fano factor of the focusing peaks, which would indicate whether the valley-polarized current is also phase-coherent, an aspect the conductance oscillations alone do not reveal."],"forward_implications":["Quantum transport in micrometer-sized strained graphene devices with nearly uniform pseudomagnetic fields becomes computationally feasible without full atomistic resolution.","A bent graphene ribbon acts as a strain-only valley splitter: transverse pseudomagnetic focusing sends one valley toward the collector and the other away, with conductance peaks that are asymmetric between inward and outward bending.","In an S-shaped ribbon, pseudomagnetic snake states bound to the sign change of $B_s$ propagate oppositely for the two valleys and produce conductance oscillations as the strain is tuned.","When a real magnetic field and a pseudomagnetic field coexist, the Landau-level ladder is $E_m(B_z \\pm B_s)$, and the zeroth Landau level changes sublattice support at $B_z = B_s$, giving a handle to distinguish strain effects from magnetic ones."],"supporting_citations":[{"why":"Introduces the scalable tight-binding model for pristine graphene that this paper generalizes to strain.","marker":"[26]"},{"why":"Provides the tight-binding treatment of uniaxial strain, giving the hopping modulation and pseudogauge-field connection used here.","marker":"[4]"},{"why":"Establishes the gauge-field description of strain in graphene that underlies the pseudomagnetic field concept.","marker":"[2]"},{"why":"Shows that bending graphene ribbons generates quantizing pseudomagnetic fields, motivating the bent-ribbon focusing device.","marker":"[7]"},{"why":"Supplies the triaxial displacement field that yields a uniform pseudomagnetic field, used to benchmark the scaling via pseudo Landau levels.","marker":"[8]"},{"why":"Source of the reduction factor kappa that accounts for optical displacement contributions to the pseudogauge field.","marker":"[30]"},{"why":"Also contributes the phonon-based derivation of the kappa reduction and elastic screening of pseudogauge fields.","marker":"[31]"},{"why":"Experimental demonstration of transverse magnetic focusing in graphene, the real-field analog whose conductance peaks the strained device mimics.","marker":"[39]"},{"why":"Shows snake trajectories in graphene p-n junctions, the analog for the pseudomagnetic snake states predicted in the S-shaped ribbon.","marker":"[44]"}],"fun_headline_variants":["Scaling law keeps graphene's pseudomagnetic fields intact","Coarse-graining graphene preserves strain-induced magnetic fields","Strain-induced magnetic fields survive graphene lattice coarsening","New scaling method models realistic graphene straintronics","Graphene pseudomagnetic fields persist under scaling trick"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method assumes the displacement field varies slowly enough that the bond-vector change can be linearized as $(d^0_n \\cdot \\nabla)\\mathbf{u}(r)$; when the scaled strain becomes large this breaks down, and the pseudomagnetic field is no longer preserved, as the authors demonstrate for large displacements and at $s=4$ for fields above 3 T.","fun_headline_variants_meta":{"raw":{"variants":["Scaling law keeps graphene's pseudomagnetic fields intact","Coarse-graining graphene preserves strain-induced magnetic fields","Strain-induced magnetic fields survive graphene lattice coarsening","New scaling method models realistic graphene straintronics","Graphene pseudomagnetic fields persist under scaling trick"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000201,"raw_usage":{"total_tokens":1405,"prompt_tokens":995,"completion_tokens":410,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":611,"completion_tokens_details":{"reasoning_tokens":332}},"tokens_in":611,"tokens_out":410,"duration_ms":4124,"temperature":1.0,"reasoning_tokens":332,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:37:45.488037+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct check is to compare the local density of states or two-terminal conductance of a scaled strained ribbon at $s=4$ against $s=1$ for a pseudomagnetic field above 3 T; the authors report in the Supplemental Material that the conductance plateaus deviate from the predicted positions in exactly this regime. Experimentally, measuring the transverse-focusing peak positions in a bent graphene ribbon under strong bending would test whether the scaled model's predictions hold at large strain.","supporting_citations":[{"cited_title":"Zhang, W","cited_arxiv_id":null,"evidence_quote":"Introduces the scalable tight-binding model for pristine graphene that this paper generalizes to strain."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the tight-binding treatment of uniaxial strain, giving the hopping modulation and pseudogauge-field connection used here."},{"cited_title":"Nigge, A","cited_arxiv_id":null,"evidence_quote":"Shows that bending graphene ribbons generates quantizing pseudomagnetic fields, motivating the bent-ribbon focusing device."},{"cited_title":"Guinea, A","cited_arxiv_id":null,"evidence_quote":"Supplies the triaxial displacement field that yields a uniform pseudomagnetic field, used to benchmark the scaling via pseudo Landau levels."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the reduction factor kappa that accounts for optical displacement contributions to the pseudogauge field."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Also contributes the phonon-based derivation of the kappa reduction and elastic screening of pseudogauge fields."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experimental demonstration of transverse magnetic focusing in graphene, the real-field analog whose conductance peaks the strained device mimics."},{"cited_title":"Chaves, L","cited_arxiv_id":null,"evidence_quote":"Shows snake trajectories in graphene p-n junctions, the analog for the pseudomagnetic snake states predicted in the S-shaped ribbon."}],"review_version":1}