{"id":"cb8323b1-a969-49bc-bc42-5955e83cef52","arxiv_id":"2411.17346","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Topological insulator surface states in a magnetic field produce strong, tunable second-harmonic generation governed by hexagonal-warping-modified Landau level selection rules.","lead":"This paper calculates how topological insulator surface states convert light to double its frequency when a perpendicular magnetic field is applied, finding sharp resonances from Landau level transitions. The predicted second-harmonic response is extremely strong and tunable by magnetic field and doping, which could make these surfaces useful for terahertz and infrared nonlinear optics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The '10^7 pm/V' susceptibility claim is not an invariant result: it follows from dividing a 2D sheet conductivity by an assumed surface-state thickness d≈9.4 Å, with no validation against the actual field profile or radiation boundary conditions, so the bulk-material comparison is not secured.","rationale":"I read the paper as a model calculation of SHG from TI surface Landau levels with hexagonal warping. The strong points are the microscopic conductivity formula, the perturbative selection rules, and the exact diagonalization spectra. The selection rules for one- and two-photon resonances are internally consistent: they follow from the Δn=3 coupling of the Fu Hamiltonian and the Berry connection selection rule |s1|-|s2|=3l+1. The numerical spectra show peaks whose positions match the Landau-level transition energies. The reader's verdict CONDITIONAL is appropriate. The weakest link is the conversion from 2D sheet conductivity to an effective 3D susceptibility. A 2D sheet response is characterized by the sheet conductivity; converting it to a bulk χ by dividing by a heuristic penetration depth is a convention that is not validated in the paper. The quoted 10^7 pm/V therefore is not a settled physical result and should not be compared to bulk crystals without a boundary-condition calculation. My proposed check would settle whether the d-division is equivalent to the actual SHG radiation. I agree with the reader's assessment; no change in verdict is needed (remain CONDITIONAL).","tokens_in":13691,"tokens_out":24255,"duration_ms":218509,"concrete_test":"For the peak transition at B=5 T, compute the reflected SHG amplitude from the 2D sheet conductivity σ_2D using the exact boundary conditions for a nonlinear surface current at the TI interface (e.g., Sipe's surface-current method). Then compute the reflected SHG from a homogeneous slab of thickness d=ℏv_F/Δ_gap with χ_eff=iσ_2D/(2ωε0 d) and identical linear optical constants. If the two reflected amplitudes differ by more than ~30%, or if recomputing χ_eff with d taken as the actual wavefunction decay length (e.g., 2-3 nm) changes the quoted 10^7 pm/V by more than an order of magnitude, the bulk-comparison claim should be removed or explicitly redefined as a 2D sheet response.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—effective surface SHG susceptibility up to ~10^7 pm/V—is obtained by converting the 2D sheet conductivity to a bulk-like χ using d=ℏv_F/Δ_gap≈9.4 Å. This conversion presumes the surface response can be represented as a homogeneous 3D slab of thickness d, ignoring (i) the true z-profile of the surface-state wavefunction (the current is a sheet, not a uniform slab), (ii) local-field and dielectric-screening effects at the vacuum/TI interface, and (iii) the fact that the measurable quantity is the radiated SHG from a surface current, which is not directly a bulk χ. If the appropriate nonlinear thickness is a few times larger, or the field distribution is nonuniform, χ_eff drops by orders of magnitude, and the comparison to conventional bulk nonlinear crystals (and to twisted bilayer graphene) loses meaning. The spectral/selection-rule part of the paper is robust and independently supported by the perturbative analysis; the magnitude headline is the load-bearing weak point.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies second harmonic generation (SHG) from the surface states of three-dimensional topological insulators in a perpendicular quantizing magnetic field. Using a Fu-type surface Hamiltonian with hexagonal warping and Zeeman coupling, the authors construct Landau levels numerically and evaluate the second-order magneto-optical conductivity given in Eq. (4). A first-order perturbative treatment of the warping yields approximate selection rules: one-photon inter-Landau-level transitions require |s1|-|s2| = τ or -2τ, while two-photon transitions require |s1|-|s2| = -τ or 2τ, where τ labels circular polarization. The numerical spectra are then interpreted by assigning resonant peaks to these one- and two-photon transitions, with additional peaks appearing for different chemical potentials via Pauli blocking and intraband transitions. The paper's main quantitative claim is that the effective surface SHG susceptibility can reach about 10^7 pm/V in the terahertz/infrared range, tunable by magnetic field and doping.","tokens_in":13922,"tokens_out":5857,"duration_ms":54686,"significance":"If the spectral and selection-rule results are correct, the paper provides a concrete route to tunable THz/IR second-harmonic generation from topological surface Dirac electrons and establishes hexagonal warping as the essential symmetry-breaking ingredient. The strengths of the manuscript are its analytical selection rules, the internal consistency between the perturbative transition amplitudes and the numerically computed peak positions, the finite-temperature Fermi-Dirac treatment, and the independent check against the gapped-graphene limit at small magnetic field. The paper also makes clear which response tensor components are allowed by the reduced C3 symmetry. The main quantitative claim, however, rests on an unvalidated conversion from a 2D sheet conductivity to a bulk-like χ, and the numerical diagonalization lacks convergence documentation; both points need addressing before the headline result can be accepted.","major_comments":[{"comment":"The quantitative headline that the effective SHG susceptibility reaches about 10^7 pm/V is obtained by dividing the 2D sheet conductivity by a thickness d = ℏv_F/Δ_gap ≈ 9.4 Å. This conversion is not derived in the paper: the surface-state current is a 2D sheet with a material-specific wavefunction decay profile, not a homogeneous slab, and the vacuum/TI interface is subject to local-field, dielectric-screening, and radiation-boundary effects that are absent from Eq. (4). The measurable quantity is the SHG radiated by a surface current, which is not automatically expressible as a bulk χ. Because the comparison to conventional bulk nonlinear materials and to twisted bilayer graphene in Refs. [46]-[52] is a central selling point, the authors should either carry out a proper electromagnetic calculation of the reflected SHG from the 2D nonlinear sheet or explicitly reframe the result as a 2D nonlinear conductivity with an honest discussion of the assumed effective thickness. As written, the claim is inversely proportional to an arbitrary d; an order-of-magnitude larger effective thickness would erase the advertised enhancement.","section":"Section 3.2"},{"comment":"","section":"Section 3.2"},{"comment":"","section":"Section 3.1"}],"minor_comments":[{"comment":"","section":"Equation (4) and following text"},{"comment":"","section":"References [19] and [20]"},{"comment":"","section":"Figure 4 caption"},{"comment":"","section":"Figure 2 caption"},{"comment":"","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":""},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core of this paper is a careful calculation of second harmonic generation from Landau-quantized topological insulator surface states with hexagonal warping. I find the perturbative selection rules convincing: for one-photon transitions |s1|-|s2| = τ or -2τ, for two-photon -τ or 2τ, where τ is the circular polarization. The exact numerical spectra show peaks at exactly the transition energies identified from the S and T amplitudes, so the internal consistency is strong. The weak-field limit also returns the expected gapped-graphene-like response. That is a genuinely new and useful result: prior nonlinear TI work used magnetic impurities or circular polarization, while an external field gives a continuously tuneable response.\n\nWhat the paper does not secure is the headline magnitude. The claim of effective susceptibility up to 10^7 pm/V comes from dividing the 2D sheet conductivity by d = ℏv_F/Δ_gap ≈ 9.4 Å. That is a unit conversion, not a model of the surface nonlinearity. The actual SHG radiation is produced by a surface current with a finite penetration depth, not a homogeneous slab, and local-field and dielectric-confinement effects at the vacuum/TI interface are ignored. A factor of a few in the effective thickness changes the conclusion by orders of magnitude. So the comparison to bulk crystals and twisted bilayer graphene is not established. This is the load-bearing weak point for the practical-device narrative.\n\nTwo smaller issues: the \"numerically exact\" Landau level calculation is not documented with truncation or convergence data, so the exactness claim can't be checked; and the linewidth Γ = 1.3 meV is introduced phenomenologically without justification. Neither affects the spectral assignments, but they should be addressed in revision.\n\nOverall I would send this to a knowledgeable referee. The core calculation is novel, reproducible in principle, and likely to matter for THz/IR nonlinear studies of TIs. The magnitude claim needs to be either softened or given a proper treatment of surface radiation. A referee should ask for that, but this is the kind of paper that deserves referee time.\n\nWho's it for: people working on nonlinear optics of Dirac materials, magneto-optical responses of TIs, and anyone planning THz SHG experiments on Bi2Se3/Bi2Te3 under a magnetic field.\n\nRecommendation: accept with major revision focus on the thickness conversion and numerical convergence.","headline":"Solid and internally consistent theory for SHG selection rules in TI surface Landau levels, but the 10^7 pm/V susceptibility claim rests on an unvalidated unit conversion.","tokens_in":14457,"tokens_out":4715,"would_cite":true,"duration_ms":45153,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["78.20.Ls","42.65.Ky"],"model":"deepseek-v4-flash","headline":"This paper argues that a perpendicular magnetic field gives the Dirac surface states of a three-dimensional topological insulator a strong, tunable second-harmonic response through the combined action of hexagonal warping and Landau…","keywords":["second harmonic generation","topological insulator surface states","Landau levels","hexagonal warping","magneto-optical conductivity","selection rules","terahertz nonlinear optics","circular polarization"],"falsifier":"At $B=5$ T with the chemical potential set so the highest occupied level is $s=0$, measure the SHG spectrum for circularly polarized normal incidence: the selection-rule picture predicts a first two-photon peak at $\\hbar\\omega=(\\varepsilon_2-\\varepsilon_0)/2$ (the $T^+_{2,0}$ transition) and a first one-photon peak at $\\hbar\\omega=\\varepsilon_1-\\varepsilon_0$, with the output in the opposite circular polarization. Observing no such peaks at these energies, or seeing the same polarization as the input, would refute the claim; a cheaper check is that all peak positions should scale as $\\sqrt{B}$ at fixed level index.","tokens_in":2306,"feed_emoji":"🧲","tokens_out":9979,"duration_ms":146625,"temperature":0.7,"pith_summary":"This paper argues that the Dirac surface states of a three-dimensional topological insulator, placed in a perpendicular quantizing magnetic field, develop a strong and tunable second harmonic generation (SHG) response. The mechanism is that the hexagonal warping of the surface band structure becomes decisive once the continuum of Dirac states collapses into discrete Landau levels, providing the dipole matrix elements and selection rules that allow a circularly polarized photon to be up-converted. The paper derives these selection rules perturbatively and verifies them against numerically exact Landau-level spectra. The payoff is a predicted effective surface susceptibility up to about $10^7$ pm/V in the far-infrared and terahertz range, controlled by magnetic field and chemical potential. If correct, this makes magnetic-field-tuned topological insulator surfaces a practical nonlinear medium rather than a theoretical curiosity.","feed_headline":"Magnetic field unlocks 10^7 pm/V harmonics on topological surfaces","feed_subtitle":"Quantized Landau levels plus hexagonal warping make the surface SHG tunable by field and doping in the THz/IR.","key_machinery":"The central object is the second-order magneto-optical conductivity $\\sigma^{\\tau\\alpha\\beta}(\\omega)$ written in the circular-polarization basis, following from the density-matrix expression used in the paper: $$\\$sigma^{{\\tau\\alpha\\beta}}$(\\omega) = -\\frac{i $e^{4}$ B}{2\\pi\\$hbar^{2}$}\\sum_{s_1s_2s} \\frac{\\hbar\\omega_{s_2s_1}\\$xi^{{\\bar\\tau}}$_{s_2s_1}(\\xi^\\alpha_{s_1s}\\xi^\\beta_{ss_2}+\\xi^\\beta_{s_1s}\\xi^\\alpha_{ss_2})}{2\\hbar\\omega-\\hbar\\omega_{s_1s_2}+i\\Gamma} \\left(\\frac{f_{s_2s}}{\\hbar\\omega-\\hbar\\omega_{ss_2}+i\\Gamma}-\\frac{f_{ss_1}}{\\hbar\\omega-\\hbar\\omega_{s_1s}+i\\Gamma}\\right).$$ The machine that carries the argument is the Berry connection $\\xi^+_{s_1s_2}=-il_c\\int dx\\,\\Phi^\\dagger_{s_1}\\hat a^\\dagger\\Phi_{s_2}$ between Landau levels. Without warping it is nonzero only for $|s_1|-|s_2|=1$; the warping term $\\sqrt{2}\\lambda(\\hbar/l_c)^3[(\\hat a^\\dagger)^3+\\hat a^3]\\sigma_z$ mixes levels differing by 3, so the connection becomes nonzero for $|s_1|-|s_2|=3l+1$. Inserting these into the $\\Gamma\\to0$ resonant decomposition splits the conductivity into one-photon ($S^\\tau_{s_1s_2}$) and two-photon ($T^\\tau_{s_1s_2}$) amplitudes, from which the selection rules and the assignment of every spectral peak follow.","core_discovery":"The paper's central claim is that placing the surface states of a three-dimensional topological insulator in a perpendicular quantizing magnetic field produces a finite, highly tunable second harmonic response, even for normally incident light, and that the effect is controlled by the interplay of hexagonal warping and Landau quantization. With $C_{3v}$ symmetry reduced to $C_3$ by the field, the allowed conductivity components satisfy $\\sigma_{xyy}=\\sigma_{yyx}=\\sigma_{yxy}=-\\sigma_{xxx}$ and $\\sigma_{yxx}=\\sigma_{xxy}=\\sigma_{xyx}=-\\sigma_{yyy}$, so circularly polarized input generates a second-harmonic current in the opposite circular polarization. Computing the second-order magneto-optical conductivity from the density-matrix expression and diagonalizing the warped Landau Hamiltonian numerically, the paper finds sharp imaginary-part peaks at one-photon ($\\hbar\\omega=\\hbar\\omega_{s_1s_2}$) and two-photon ($\\hbar\\omega=\\hbar\\omega_{s_1s_2}/2$) inter-Landau-level resonances. A perturbation treatment in the warping parameter yields the selection rules $|s_1|-|s_2|=\\tau$ or $-2\\tau$ for one-photon and $-\\tau$ or $2\\tau$ for two-photon transitions, and identifies every peak with a specific transition (e.g. $T^+_{2,0}$, $T^+_{1,-2}$, $S^+_{s+1,-s}$). In weak fields the response becomes continuous with peaks at $\\hbar\\omega\\approx|\\mu|$ and $2|\\mu|$, $\\sigma_{xxx}\\propto B$, and $\\sigma_{yyy}$ roughly field-independent. The paper concludes that the effective surface susceptibility reaches about $10^7$ pm/V in the THz/IR, controlled by magnetic field and chemical potential.","pith_inferences":["Beyond the paper, the $10^7$ pm/V number is a conversion, not a direct observable; the robust prediction is the sheet conductivity itself, and the effective bulk susceptibility will move by orders of magnitude if the nonlinear thickness is set by the surface-state wavefunction penetration rather than by $d=\\hbar v_F/\\Delta_{\\rm gap}$.","A clean experimental signature that avoids the thickness issue is the $\\sqrt{B}$ scaling of the first few peak positions at fixed level index, which follows from the Landau-level energies and is independent of the susceptibility normalization.","The same $\\tau$ and $-2\\tau$/$2\\tau$ selection-rule pattern should appear in any two-dimensional Dirac or Rashba system with cubic warping under a quantizing field; the specific transition labels will change, but the circular-polarization-dependent up-conversion mechanism is generic."],"forward_implications":["For a circularly polarized input at normal incidence, the SHG output has the opposite circular polarization; the allowed components reduce to $\\sigma_{+--}$ and $\\sigma_{-++}$, and their half-sum gives $\\sigma_{xxx}$.","The SHG spectrum is discrete and magnetically tunable: as $B$ increases, all resonant peaks shift to higher photon energy and grow in amplitude because the Landau-level degeneracy increases.","Tuning the chemical potential through the Landau levels controls the response: Pauli blocking removes interband peaks, while intraband transitions around the highest occupied level add new peaks with large amplitudes.","In the limit of small magnetic field, the Landau levels merge into a continuum, the spectrum peaks at $\\hbar\\omega\\approx|\\mu|$ and $2|\\mu|$, and $\\sigma_{xxx}(\\omega)$ is linear in $B$ while $\\sigma_{yyy}(\\omega)$ is essentially $B$-independent.","With the surface-state thickness estimated at $d\\approx9.4$ Å, the effective SHG susceptibility reaches about $10^7$ pm/V in the THz/IR, exceeding many bulk nonlinear materials and comparable to twisted bilayer graphene."],"supporting_citations":[{"why":"Supplies the hexagonal warping Hamiltonian $\\lambda/2\\hbar^3(p_+^3+p_-^3)\\sigma_z$ that breaks the continuous rotational symmetry and generates the SHG.","marker":"[25]"},{"why":"Gives the density-matrix expression for the second-order magneto-optical conductivity, Eq. (4), from which all spectra and selection rules are computed.","marker":"[42]"},{"why":"Provides the experimental and theoretical framework for separating surface-state SHG in Bi2Se3 and supplies the thickness estimate $d=\\hbar v_F/\\Delta_{\\rm gap}\\approx9.4$ Å used to convert sheet conductivity to susceptibility.","marker":"[13]"},{"why":"Establishes that the surface-state SHG response is tunable by gating and doping, the experimental handle the predicted field and doping tunability builds on.","marker":"[21]"},{"why":"Reports terahertz SHG from Bi2Se3 Dirac surface states, the frequency regime the predicted resonances target.","marker":"[22]"},{"why":"Computes optical nonlinearity in doped and gapped graphene and underlies the weak-field limit comparison with peaks at $\\hbar\\omega\\approx|\\mu|$ and $2|\\mu|$.","marker":"[45]"}],"fun_headline_variants":["Magnetic fields make topological surfaces emit 10^7 pm/V harmonics","Landau levels + warping: topological SHG hits 10^7 pm/V","Quantizing field turns topological surfaces into tunable THz emitters","Topological insulator surfaces: field-tunable second harmonic to 10^7 pm/V","Hexagonal warping unlocks giant SHG on topological surfaces in B-field"],"cache_read_input_tokens":16640,"weakest_assumption_plain":"The paper assumes the surface SHG response can be summarized as an effective bulk susceptibility by dividing the two-dimensional sheet conductivity by a fixed thickness $d=\\hbar v_F/\\Delta_{\\rm gap}\\approx9.4$ Å, with no local-field correction, surface profile, or bulk screening; if the actual nonlinear thickness is much larger or the field distribution is nonuniform, the headline $10^7$ pm/V number changes by orders of magnitude.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic fields make topological surfaces emit 10^7 pm/V harmonics","Landau levels + warping: topological SHG hits 10^7 pm/V","Quantizing field turns topological surfaces into tunable THz emitters","Topological insulator surfaces: field-tunable second harmonic to 10^7 pm/V","Hexagonal warping unlocks giant SHG on topological surfaces in B-field"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000206,"raw_usage":{"total_tokens":1496,"prompt_tokens":1147,"completion_tokens":349,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":763,"completion_tokens_details":{"reasoning_tokens":247}},"tokens_in":763,"tokens_out":349,"duration_ms":3826,"temperature":1.0,"reasoning_tokens":247,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:13:43.376436+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"At $B=5$ T with the chemical potential set so the highest occupied level is $s=0$, measure the SHG spectrum for circularly polarized normal incidence: the selection-rule picture predicts a first two-photon peak at $\\hbar\\omega=(\\varepsilon_2-\\varepsilon_0)/2$ (the $T^+_{2,0}$ transition) and a first one-photon peak at $\\hbar\\omega=\\varepsilon_1-\\varepsilon_0$, with the output in the opposite circular polarization. Observing no such peaks at these energies, or seeing the same polarization as the input, would refute the claim; a cheaper check is that all peak positions should scale as $\\sqrt{B}$ at fixed level index.","supporting_citations":[{"cited_title":"A.; Muryumin, E.; Gaiduk, E","cited_arxiv_id":null,"evidence_quote":"Supplies the hexagonal warping Hamiltonian $\\lambda/2\\hbar^3(p_+^3+p_-^3)\\sigma_z$ that breaks the continuous rotational symmetry and generates the SHG."},{"cited_title":"Nonlinear response in a noncentrosymmetric topological insulator","cited_arxiv_id":null,"evidence_quote":"Provides the experimental and theoretical framework for separating surface-state SHG in Bi2Se3 and supplies the thickness estimate $d=\\hbar v_F/\\Delta_{\\rm gap}\\approx9.4$ Å used to convert sheet conductivity to susceptibility."},{"cited_title":"Observation of terahertz second harmonic generation from Dirac surface states in the topological insulator B i _2 S e _3","cited_arxiv_id":null,"evidence_quote":"Establishes that the surface-state SHG response is tunable by gating and doping, the experimental handle the predicted field and doping tunability builds on."},{"cited_title":"Topological insulators in B i _2 S e _3 , B i _2 T e _3 and S b _2 T e _3 with a single D irac cone on the surface","cited_arxiv_id":null,"evidence_quote":"Reports terahertz SHG from Bi2Se3 Dirac surface states, the frequency regime the predicted resonances target."},{"cited_title":"M.; Alencar, T","cited_arxiv_id":null,"evidence_quote":"Computes optical nonlinearity in doped and gapped graphene and underlies the weak-field limit comparison with peaks at $\\hbar\\omega\\approx|\\mu|$ and $2|\\mu|$."}],"review_version":1}