{"id":"a35eaf8e-10f0-4aab-b3f5-28940fb43bee","arxiv_id":"2411.14075","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Photon drag at metal-2D semiconductor contacts is predicted to be strong, with p-polarized fields giving a conductivity-independent momentum transfer and s-polarized fields scaling as η ln(1/η).","lead":"This paper calculates how light pressure on electrons, known as photon drag, behaves at the boundary between a metal contact and a two-dimensional semiconductor. It predicts the effect is much stronger than previously assumed and can dominate the usual thermoelectric response at terahertz frequencies.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central prediction rests on the local Ohm's law near a singular edge; for small η the drag is dominated by a near-edge region of width ~ηλ0 where nonlocal corrections become significant, so the numerically obtained α_x≈0.45 may not be universal.","rationale":"The most load-bearing step in the paper is the conversion of the exact (local-response) diffraction fields into a drag force using the same local conductivity. The reader identifies exactly this point, and I agree: the p-polarization drag is dominated by the near-edge region whose width l~ηλ0 shrinks as η→0. For the experimentally relevant graphene case η≈0.01 and λ0≈100 μm, l≈1 μm, comparable to the momentum relaxation mean free path at low temperature; the local approximation can fail by order-unity factors. The claimed constant α_x≈0.45 is obtained numerically from this local model, so it is not a protected quantity. The paper's own statement that the only assumption is locality (Sec. II.A) and the absence of an estimate of when it breaks make the central formula conditional. The s-polarization 'always dominates' claim depends on a local-heating approximation for the thermoelectric effect that likely underestimates the temperature at the contact when heat diffusion is included, but this is secondary compared with the validity of the drag force itself. A concrete nonlocal simulation would settle whether the constants and scaling survive.","tokens_in":13769,"tokens_out":16377,"duration_ms":163845,"concrete_test":"Perform a self-consistent numerical simulation of a metal-2DES junction solving the Boltzmann equation for the 2D electron gas (with a finite relaxation time τ) coupled to Maxwell's equations for the same half-plane geometry, for η=0.01 and η=0.001, with l_mfp = v_F τ such that l_mfp/λ0 = 0.1 and 1.0. Extract the drag photovoltage via Eq. (13) using the simulated fields and the nonlocal current, and compare the resulting α_x to the local value 0.45. If the deviation exceeds 20% in any of these cases, the local-Ohm's-law assumption is not justified and the reported α_x is model-dependent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The derivation of the drag photovoltage (Eqs. 4, 13, 18–19) assumes the local relation j_ω(x)=σ_ω(x)E_ω(x) everywhere, as stated in Sec. II.A. For p-polarization, α_x is dominated by the near-edge field, which grows as ~1/√η and varies over a length l~ηλ0 (Eq. 21 and the surrounding text). For small η (e.g., η≈0.01 in graphene at THz), l is in the range 0.1–1 μm, comparable to or shorter than the electron mean free path in high-mobility 2D systems at low temperature. In this regime the local conductivity is an idealization: the true current response is nonlocal, j(q)=σ(q,ω)E(q), and the electromagnetic field itself is modified at length scales below the mean free path. Since α_x is computed as an integral over the near-edge region (Eq. 19), a q-dependent σ(q,ω) can change the constant 0.45 and possibly the η→0 limiting value. The paper explicitly acknowledges the locality assumption but provides no estimate of the breakdown scale; this makes the central prediction of a universal, polarization-dependent α less secure. The s-polarization dominance claim is an additional concern, but the locality issue is more load-bearing because it affects the drag formula itself.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a theory of photon drag at the junction between a metal and a two-dimensional electron system (2DES). The authors combine exact Wiener-Hopf diffraction spectra for the local electromagnetic fields at a metal-2DES edge with a microscopic transport-theory expression for the ponderomotive force on the 2D carriers. Their central result, Eq. (18), is a photovoltage responsivity r_pd = -2/(omega e n_2d)(alpha_x + alpha_y), where the dimensionless momentum-transfer coefficients alpha_i depend only on the normalized 2D conductivity eta = 2*pi*sigma/c and on polarization. For p-polarized light, alpha_x approaches a finite constant (~0.45) as eta -> 0, which the authors attribute to the dynamic lightning-rod field singularity at the edge; for s-polarized light, alpha_y scales as eta ln(1/eta), arising from a long-range cylindrical wave emitted by the metal half-plane. The paper also compares the junction photon drag with the contact thermoelectric effect, concluding that for p-polarization the ratio is set by omega*tau_epsilon and for s-polarization photon drag dominates.","tokens_in":14048,"tokens_out":7844,"duration_ms":78517,"significance":"If the central prediction is correct, the paper identifies a qualitatively new route to terahertz photodetection: junction photon drag is not suppressed by the small free-photon momentum or by the small absorbance of typical 2D materials, and it can exceed the thermoelectric response in common contact geometries. The derivation from the ponderomotive force, Eqs. (3)-(13), is clean and internally consistent, and the use of parameter-free Wiener-Hopf diffraction spectra from Refs. 33 and 34, with no fitted parameters, is a genuine strength. The two scaling laws (alpha_x -> const and alpha_y ~ eta ln eta^{-1}) are falsifiable predictions that can be tested in photocurrent-mapping experiments. However, the universality of these predictions rests on assumptions whose breakdown scales are not quantified, in particular the locality of the Ohm's law in the singular near-edge region and the treatment of the logarithmic divergence in the s-polarization channel.","major_comments":[{"comment":"The locality assumption j_omega(x) = sigma_omega(x) E_omega(x) is load-bearing for the p-polarization result. The text states that for small eta the p-polarized field diverges at the edge as E_x(0) ~ E_0/sqrt(eta(1+eta)) and varies over a length l ~ eta*lambda_0. For eta ~ 10^-2 (graphene at THz frequencies) and lambda_0 = 30-100 micrometers, l is roughly 0.3-1 micrometer, which is comparable to or shorter than the electron mean free path in high-mobility 2D systems at low temperature. In this regime a q-dependent conductivity sigma(q,omega) modifies both the field profile and the integral in Eq. (19), so the numerical value alpha_x ~ 0.45 and its claimed eta -> 0 limit are not established. The manuscript acknowledges the locality assumption but provides no estimate of the breakdown scale. Please quantify the nonlocal correction, for example by comparing q ~ 1/l with the inverse mean free path, or clearly state the parameter regime in which the universal value of alpha_x is expected to hold.","section":"Section II.A and Eq. (21)"},{"comment":"The claim that alpha_y depends only on eta and polarization is complicated by the logarithmic divergence described for s-polarization. The text states that the integral (19) diverges logarithmically at large distances and is cut off by finite absorption in the 2DES at L ~ lambda_0/eta, but this cutoff is not derived. In a finite sample of length L_sd < lambda_0/eta, the coefficient will depend on ln(L_sd/lambda_0) or on the actual absorption length, so the abstract's statement that the momentum-transfer coefficient depends only on eta and polarization is not universal. A finite-length expression for alpha_y(L_sd), or an explicit statement that the scaling requires L_sd >> lambda_0/eta, is needed to make the s-polarization prediction usable for experiments.","section":"Section III, s-polarization scaling paragraph"}],"minor_comments":[{"comment":"In the sentence introducing Eqs. (12)-(13), \"Demoting the Seebeck coefficient\" should read \"Denoting the Seebeck coefficient\".","section":"Section II.B"},{"comment":"The word \"superliner\" should be \"superlinear\" in the discussion of the s-polarization scaling.","section":"Section III"},{"comment":"The second bullet in the caption describes the s-polarized case (E0 parallel to the junction), but the text labels it as \"For p-polarized incident wave\"; this should be corrected to s-polarized.","section":"Fig. 1 caption"},{"comment":"Equation (21) as typeset appears as E_x(0) = E_0 sqrt(eta(1+eta)), which contradicts the surrounding text stating that the field diverges as eta -> 0. The intended expression is presumably E_x(0) = E_0 / sqrt(eta(1+eta)); the typesetting should be corrected.","section":"Eq. (21)"}],"recommendation":"major_revision","confidential_remarks":"The numerical value alpha_x ~ 0.45 is inherited from the Wiener-Hopf solutions of Refs. 33 and 34, which are by the same group. This is not circular, since those solutions are parameter-free, but the referee should be aware that the present paper does not independently rederive that constant. The main technical risk is the locality issue in the near-edge region, which is the basis of my major_revision recommendation. The paper is otherwise well within the scope of Physical Review B."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something solid: it takes the exact Wiener–Hopf diffraction spectra at a metal–2DES junction and combines them with a transparent photon-drag transport derivation. The result is a clear, testable prediction that the junction drag photovoltage is not suppressed by the small photon momentum, with distinctive scaling laws: α_x stays finite (~0.45) as η→0 for p-polarization, and α_y ~ η ln(1/η) for s-polarization. The comparison with thermoelectric response, giving a ratio of order ωτ_ε for p-polarization, is a useful rule of thumb. The authors are also honest about the main assumptions, which is a plus.\n\nThe real soft spot is the locality assumption near the metal edge. The derivation uses j = σE everywhere, but the p-polarization drag is dominated by a near-edge region of width l ~ ηλ_0. For small η in realistic THz graphene, that width is comparable to or smaller than the electron mean free path. The paper acknowledges the assumption but does not estimate what nonlocal corrections do to α_x ≈ 0.45. This is not a fatal flaw, but it means the numerical value and possibly the η→0 limit are less secure than the text suggests. A second, milder issue is the claim that s-polarization photon drag always dominates thermoelectric effect. That rests on a local-heating approximation in which field at the junction sets the temperature rise; but heat generated over the larger strip where the s-polarized field is non-zero can diffuse to the junction, so the comparison deserves a more careful treatment.\n\nThe citation pattern is fine. The field spectra come from earlier work by the same group, but those are parameter-free exact solutions, not fitted hacks. The derivation is analytic and detailed enough for independent reimplementation. This is a good theory paper for people working on THz photodetection, photocurrent mapping, and photon drag in 2D materials. It deserves a serious referee; the concerns I raised are addressable in revision but do not invalidate the central idea.","headline":"Clean theory paper with a genuinely new prediction of strong junction photon drag in 2D materials; the main open question is whether nonlocal corrections near the metal edge change the headline numbers.","tokens_in":14592,"tokens_out":3289,"would_cite":true,"duration_ms":35593,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A metal contact to a 2D semiconductor converts incident light into a photon-drag photovoltage whose p-polarized response survives as the sheet conductivity tends to zero.","keywords":["photon drag","metal-2d semiconductor junction","terahertz photodetection","lightning rod effect","Wiener-Hopf diffraction","momentum transfer coefficient","photothermoelectric effect","two-dimensional electron gas"],"falsifier":"A direct check: measure the zero-bias photovoltage at a metal-graphene junction under focused p-polarized terahertz illumination while a back gate tunes the sheet density and hence $\\eta$ over at least a decade. The theory predicts a response that tracks $1/n_{2d}$ with $\\alpha_x\\simeq0.45$ essentially independent of $\\eta$; observation of a strong $\\eta$-dependent deviation, or of an $\\alpha_x$ that vanishes with $\\eta$, would rule out the local-field drag mechanism described here. For s-polarization, the signature is the superlinear $\\alpha_y\\propto \\eta\\ln(1/\\eta)$ growth, which should be absent if nonlocal corrections dominate.","tokens_in":13563,"feed_emoji":"⚡","tokens_out":8960,"duration_ms":79568,"temperature":0.7,"pith_summary":"The paper argues that photon drag—photocurrent produced when light transfers momentum to charge carriers—does not have to be weak at the contact between a metal and a two-dimensional semiconductor. Diffraction at the metal edge creates strongly nonuniform local fields, and the Lorentz force on induced charges in these fields does not average to zero. Combining an exact diffraction solution with a microscopic force-balance calculation, the authors derive a responsivity $r_{pd}=-2(\\alpha_x+\\alpha_y)/(\\omega e n_{2d})$, where the dimensionless momentum-transfer coefficients $\\alpha_x,\\alpha_y$ depend only on the normalized 2D conductivity $\\eta=2\\pi\\sigma/c$ and on polarization. For $p$-polarized light, $\\alpha_x$ remains finite ($\\simeq 0.45$) as $\\eta\\to0$, because the edge field diverges as $1/\\sqrt{\\eta(1+\\eta)}$; for $s$-polarized light, $\\alpha_y\\propto \\eta\\ln\\eta^{-1}$. If this is right, contact-induced photon drag can dominate the thermoelectric response at terahertz frequencies, with the ratio set by $\\omega\\tau_\\varepsilon$ for $p$-polarization and no such suppression for $s$-polarization.","feed_headline":"Photon drag at metal-2D junctions stays finite as conductivity vanishes","feed_subtitle":"P-polarized contact drag keeps a finite coefficient at vanishing 2D conductivity and dominates at terahertz.","key_machinery":"The calculation is carried by three ingredients: (i) the wave-mean-force expression $f_{pd}=\\langle\\rho E\\rangle+c^{-1}\\langle j\\times B\\rangle$, which after time-averaging and using the continuity equation and the local Ohm's law $j_\\omega(x)=\\sigma_\\omega(x)E_\\omega(x)$ becomes $\\frac{2}{\\omega}\\mathrm{Im}\\{\\sigma_\\omega[\\partial_x E_x E_x^*-\\partial_x E_y^*E_y]\\}$; (ii) the exact Wiener-Hopf diffraction spectra for the local fields, Eqs. (14) and (16), obtained from factorized 2D dielectric functions $\\varepsilon_p(q)=1+\\eta\\sqrt{k_0^2-q^2}/k_0$ and $\\varepsilon_s(q)=1+\\eta k_0/\\sqrt{k_0^2-q^2}$; and (iii) the momentum balance plus thermal balance equations that separate the photon-drag voltage from the thermoelectric voltage. The singular behavior enters through the edge field $E_x(0)=E_0/\\sqrt{\\eta(1+\\eta)}$, which keeps $\\eta|E_x|^2$ finite as $\\eta\\to0$, and through the long-range $e^{ik_0x}/\\sqrt{x}$ cylindrical wave whose logarithmic divergence in the drag integral is cut off at distance $\\sim\\lambda_0/\\eta$.","core_discovery":"At a straight interface between a perfect metal and a 2D electron system, the open-circuit photovoltage from photon drag is $V_{pd}=-(2/\\omega e n_{2d})\\int_0^\\infty \\mathrm{Im}\\{\\sigma_\\omega[\\partial_x E_x E_x^*-\\partial_x E_y^* E_y]\\}\\,dx$, and for normal incidence the responsivity factorizes into a universal $1/(\\omega e n_{2d})$ prefactor times the sum of momentum-transfer coefficients $\\alpha_x+\\alpha_y$. The $p$-polarized coefficient stays nonzero in the zero-conductivity limit because the singular edge field concentrates the drag force in a narrow region of width $\\sim\\eta\\lambda_0$; the $s$-polarized coefficient is dominated by the long cylindrical wave emitted by the half-plane and scales as $\\eta\\ln(1/\\eta)$. The same expression shows that the ordinary bulk photon drag is smaller by a factor $\\eta' k_x L$, so the junction effect removes both conventional smallness parameters: photon momentum and low absorbance.","pith_inferences":["Editorial inference: because the drag force is proportional to $\\mathrm{Im}\\,\\sigma_\\omega$ and the field profile is set by both $\\eta'$ and $\\eta''$, measuring the photovoltage with focused light at both polarizations gives separate access to the real and imaginary parts of the sheet conductivity without patterning a separate antenna.","Editorial inference: the same mechanism should contribute to the DC photoresponse of contacted 2D detectors usually fitted by plasma-wave or hot-carrier rectification, so polarization-resolved contact-illumination measurements could isolate this contact drag from those channels.","Editorial inference: the local-Ohm's-law assumption can be relaxed by replacing $\\sigma_\\omega(x)$ with a nonlocal conductivity kernel; since the $p$-polarized drag is generated within $\\sim\\eta\\lambda_0$ of the edge, the first nonlocal corrections are expected when this length becomes comparable to the electron mean free path or screening length."],"forward_implications":["At terahertz frequencies, junction photon drag can dominate the hot-carrier thermoelectric response for $p$-polarized light whenever $\\omega\\tau_\\varepsilon\\gtrsim1$, which is typical for wavelengths of about 30 $\\mu$m and shorter.","For $s$-polarized light the photovoltage is essentially pure photon drag, since the thermoelectric contribution at the edge is small and $\\alpha_y$ remains positive.","The responsivity contains no factor of photon momentum or small absorbance: the prefactor $2/(\\omega e n_{2d})$ alone gives tens of $\\mu$V cm$^2$/W at 1 THz and $n_{2d}=10^{11}$ cm$^{-2}$.","The theory applies to any 2D conductor describable by a local dynamic conductivity, including graphene, semiconductor quantum wells, and transition-metal dichalcogenides.","Plasmonic resonances do not resonantly enhance the drag: the plasmon quality factor multiplies the dissipative conductivity, so the momentum-transfer coefficients remain smooth across the complex-$\\eta$ plane."],"supporting_citations":[{"why":"supplies the exact Wiener-Hopf diffraction spectra for the local fields at metal-2D junctions used in Eqs. (14) and (16).","marker":"[33]"},{"why":"supplies the complementary exact scattering solution for s-polarized fields and transverse-electric plasmon launching at 2D junctions.","marker":"[34]"},{"why":"is the classic wedge-diffraction solution whose $x^{-1/2}$ edge singularity is the origin of the lightning-rod field enhancement.","marker":"[22]"},{"why":"gives the structured-light photon-drag formalism that the junction calculation extends to internally structured contact fields.","marker":"[11]"},{"why":"is the plasma-wave rectification theory whose local-field description the paper identifies with photon drag at contacts.","marker":"[26]"},{"why":"is the hot-carrier photocurrent theory at graphene-metal interfaces that the paper's photon-drag result is compared against.","marker":"[19]"},{"why":"provides experimental evidence of hot-carrier thermoelectric dominance that the new drag mechanism would compete with.","marker":"[37]"},{"why":"supplies the energy relaxation times used to estimate $\\omega\\tau_\\varepsilon$ and the drag-to-thermoelectric ratio.","marker":"[42]"}],"fun_headline_variants":["Lightning rod effect pins photon drag at metal-2D junction","P-polarized drag at metal-2D edge survives zero conductivity","Junction photon drag defies low photon momentum","Edge photon drag: p-mode finite as 2D conductivity vanishes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes that the current at each point responds locally to the electric field, $j_\\omega(x)=\\sigma_\\omega(x)E_\\omega(x)$; near the contact for small $\\eta$ the field changes over a length $\\sim\\eta\\lambda_0$, which can be shorter than the electron mean free path or screening length, so nonlocal spatial-dispersion corrections could shift the value of $\\alpha_x$ and the predicted scaling.","fun_headline_variants_meta":{"raw":{"variants":["Lightning rod effect pins photon drag at metal-2D junction","P-polarized drag at metal-2D edge survives zero conductivity","Junction photon drag defies low photon momentum","Edge photon drag: p-mode finite as 2D conductivity vanishes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001114,"raw_usage":{"total_tokens":4702,"prompt_tokens":1067,"completion_tokens":3635,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":683,"completion_tokens_details":{"reasoning_tokens":3564}},"tokens_in":683,"tokens_out":3635,"duration_ms":29076,"temperature":1.0,"reasoning_tokens":3564,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:34:17.710857+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct check: measure the zero-bias photovoltage at a metal-graphene junction under focused p-polarized terahertz illumination while a back gate tunes the sheet density and hence $\\eta$ over at least a decade. The theory predicts a response that tracks $1/n_{2d}$ with $\\alpha_x\\simeq0.45$ essentially independent of $\\eta$; observation of a strong $\\eta$-dependent deviation, or of an $\\alpha_x$ that vanishes with $\\eta$, would rule out the local-field drag mechanism described here. For s-polarization, the signature is the superlinear $\\alpha_y\\propto \\eta\\ln(1/\\eta)$ growth, which should be absent if nonlocal corrections dominate.","supporting_citations":[{"cited_title":"The difference in effective dielectric functions in s- and p- polarizations results in different plasmon dispersion laws","cited_arxiv_id":null,"evidence_quote":"supplies the exact Wiener-Hopf diffraction spectra for the local fields at metal-2D junctions used in Eqs. (14) and (16)."},{"cited_title":"(27) We observe that ’ordinary’ photon drag differs from the ’junction drag’ by factors of η′ (absorbance factor) and kxLsd (light momentum factor)","cited_arxiv_id":null,"evidence_quote":"gives the structured-light photon-drag formalism that the junction calculation extends to internally structured contact fields."},{"cited_title":"This results in αx = |Epl|2 |E0|2 { Reη × Qpl − Imη 2 } , (23) where Qpl = Reqp pl/Imqp pl is the quality factor of 2D plas- mons","cited_arxiv_id":null,"evidence_quote":"is the hot-carrier photocurrent theory at graphene-metal interfaces that the paper's photon-drag result is compared against."},{"cited_title":"In the far infrared range ( λ0 = 10 µm), we estimate ωτε ≫ 1 even for the shortest energy relaxation time τε = 0.1 ps","cited_arxiv_id":null,"evidence_quote":"supplies the energy relaxation times used to estimate $\\omega\\tau_\\varepsilon$ and the drag-to-thermoelectric ratio."}],"review_version":1}