{"id":"715ea853-9158-481f-b89d-4d8bdcfd7f3d","arxiv_id":"2506.07745","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Germanium doping in graphene opens a bandgap and induces valley-contrasting Berry curvature and second harmonic generation, according to density functional theory calculations.","lead":"This computational study uses density functional theory to show that replacing some carbon atoms in graphene with germanium opens a bandgap, breaks inversion symmetry, and creates opposite Berry curvature at the K and K' valleys. The authors also report second harmonic generation and argue the material could be useful for valleytronics and nonlinear optics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"SHG pillar is unsupported: EPW (an electron–phonon code) is listed as extracting chi^(2), with no formula, projectors, or convergence tests; an independent check is required before the nonlinear-optical claim stands.","rationale":"The reader's weakest-assumption analysis correctly identified Wannier-based post-processing as a risk. I go further and flag a specific technical mismatch: EPW is an electron-phonon Wannier code, not an SHG code, and the text provides no formula or numerical details for chi^(2). This makes the nonlinear-optical pillar of the central claim the least secure. The bandgap and Berry-curvature parts are more credible: breaking inversion symmetry by substitutional Ge doping naturally opens a gap and produces opposite Berry curvature at K/K' under time-reversal symmetry, and the small Chern numbers are consistent with a trivial topological phase. However, the SHG results in Figure 9, which are used to confirm inversion-symmetry breaking and to support optoelectronic applications, cannot be assessed from the manuscript as written. Because this concern is addressable by the authors providing details and an independent reproduction, the appropriate verdict remains conditional rather than a rejection: the paper should not be accepted as-is, but no conclusion about the entire study should be drawn until the SHG methodology is verified.","tokens_in":10134,"tokens_out":5733,"duration_ms":73080,"concrete_test":"Ask the authors to provide the exact SHG workflow: the chi^(2) formula, Wannier projectors, disentanglement windows, k-grid sizes, and the precise role of EPW. Then independently recompute chi^(2)(omega) for the 12.5% and 2% Ge supercells using a dedicated nonlinear-optics code (e.g., ABINIT's DFPT nonlinear module or an independent Wannier-based shift-current routine) with the same relaxed geometries. If the main peaks in Figure 9 do not appear within roughly 0.1-0.2 eV, or if their signs or magnitudes change substantially, the SHG claim is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim has two quantitative pillars: valley-contrasting Berry curvature and tunable SHG. The nonlinear-optical pillar is the least secure. Section 2 ('Materials and computational methods') states that 'SHG coefficients were extracted using EPW and custom scripts.' EPW is a Wannier-based electron-phonon coupling package and does not implement second-order susceptibilities by itself; the accompanying mention of ph.x also points to phonons, not chi^(2). No expression for chi^(2) (velocity-gauge or length-gauge sum-over-states), no Wannier projectors or disentanglement windows, no interpolation k-mesh, no scissor correction, and no benchmark comparison are reported. Because Figure 9 is presented as evidence of broken inversion symmetry and application potential, the SHG results are currently unverifiable. This is a concrete methodological red flag, not merely a missing convergence test: if the authors cannot clarify and reproduce the workflow, the nonlinear-optical claim should not be accepted.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports first-principles density functional theory (DFT) calculations, using PBE and HSE06 functionals with and without spin-orbit coupling, on Ge-doped graphene monolayers at nominal concentrations of 12.5%, 5.5%, 3.125%, and 2%. The authors compute band structures, density of states, Berry curvature, Chern numbers, optical absorption, and second-harmonic generation, and claim that Ge doping opens a bandgap, breaks inversion symmetry while preserving time-reversal symmetry, produces valley-contrasting Berry curvature with a potential valley Hall effect, and yields finite second-order susceptibilities. The central claim is that the bandgap, valley polarization, and nonlinear optical response can be tuned by varying the Ge doping concentration, making Ge-doped graphene a candidate for valleytronic and optoelectronic applications.","tokens_in":10327,"tokens_out":7448,"duration_ms":79344,"significance":"If the reported results are correct and reproducible, this paper would provide a systematic study of a chemically accessible graphene derivative for valleytronics and nonlinear optics. The use of two exchange-correlation functionals, several doping concentrations, and spin-orbit coupling is a constructive approach, and the reported opposite-sign Berry curvature at K and K' is consistent with the expected broken-inversion-symmetry physics. However, the significance is heavily conditional on verification of the SHG workflow and the Wannier-interpolation parameters, neither of which is provided. As it stands, the nonlinear-optical and valley Hall effect pillars are not yet established to the standard required for the claimed application-level conclusions.","major_comments":[{"comment":"The text states that 'the width of bandgap increases with the increasing concentration of Ge doping' (paragraph near Figure 3), but Table 1 contradicts this monotonic claim: for the hybrid polarized calculations, the 3.125% doped supercell has a gap of 0.26 eV while the 5.5% doped supercell has a smaller gap of 0.168 eV; the hybrid unpolarized 4x4 supercell gap (0.273 eV) also exceeds the 3x3 gap (0.00 eV). Since the abstract presents doping-concentration tuning of the bandgap as a central result, the authors must either report the actual non-monotonic concentration dependence or explain the discrepancy and revise the associated discussion in Figures 3 and 4.","section":"Section 3, Table 1"},{"comment":"The SHG workflow is not reproducible as written. The text states that 'SHG coefficients were extracted using EPW and custom scripts,' but EPW is a Wannier-based electron-phonon package and does not implement second-order optical susceptibilities by itself. No expression for chi^(2) (velocity-gauge or length-gauge sum-over-states), no Wannier projectors, disentanglement windows, interpolation k-meshes, broadening parameters, or scissor corrections are reported, and no benchmark comparison is provided. Figure 9 therefore cannot be independently verified, and the central nonlinear-optical claim is unsupported. The authors should provide the exact chi^(2) formula, a step-by-step description of the post-processing workflow, and convergence tests, or the SHG pillar should be removed.","section":"Section 2, Materials and computational methods"},{"comment":"The Berry curvature results rely on Wannier interpolation via postw90, but the manuscript reports none of the required Wannier90 input parameters (projector functions, disentanglement windows, frozen windows, interpolation grids) and no convergence tests for the Berry curvature or the Chern number. Without these details, the magnitudes of Omega_z at K and K' and the small Chern numbers in Table 2 (on the order of 10^-4 to 10^-3) cannot be distinguished from interpolation artifacts. Please document the Wannierization parameters and add a convergence study with respect to the interpolation k-mesh and the Wannier subspace size.","section":"Section 2, Figure 7, Table 2"},{"comment":"The symmetry relations stated for the Berry curvature are mutually inconsistent. The text writes 'Omega_n(k) = Omega_n(-k). This clearly indicates that the inversion symmetry in the system is broken while the time reversal symmetry i.e. Omega_n(k) = -Omega_n(-k) is preserved.' The first relation is the signature of preserved inversion symmetry, not broken inversion; with time-reversal symmetry the correct relation for a system with broken inversion is Omega_n(-k) = -Omega_n(k). In addition, the conclusion that time-reversal symmetry is preserved 'as indicated by zero Chern number' is not logically valid: a zero Chern number can also occur in time-reversal-broken systems and does not by itself prove TRS. These statements should be corrected.","section":"Section 3, Berry curvature discussion and Conclusions"}],"minor_comments":[{"comment":"The manuscript contains numerous typographical errors and inconsistent notation, including 'valletronics', 'dopped', 'prinstine', 'SGH' instead of SHG, 'the energy band H amiltonian', and 'has is the corresponding geometrical phase'. A careful proofreading pass is needed.","section":"Throughout"},{"comment":"Equation (2) omits the reduced Planck constant and the vector character of the Pauli matrices; Equation (5) uses an undefined term 'ℏ l k_x sigma_0'. Since these equations are not used in the calculations, they can be streamlined or corrected to avoid confusion.","section":"Equations (1)-(5)"},{"comment":"The column headings 'Polarized' and 'Unpolarized' are confusing because Section 3 reports that the spin-up and spin-down DOS are identical and the magnetic moment is zero; please clarify what 'polarized' means in each case and why the hybrid polarized results differ from the unpolarized ones if there is no net spin polarization.","section":"Table 1"},{"comment":"The band structure plots do not include energy axis labels, Fermi-level markings, or consistent k-path tick labels in all panels, which makes the reported gaps difficult to read from the figures.","section":"Figures 2-4"},{"comment":"The text refers to 'a lower unspecified concentration' in the list of dopant concentrations; this should be identified as the 2% concentration shown in Figure 1(d).","section":"Section 2 and Figure 1"},{"comment":"The SHG and absorption plots do not specify the polarization geometry or which Cartesian tensor component of chi^(2) is shown; please state the component and the light polarization used in the calculation.","section":"Figures 9 and 10"}],"recommendation":"major_revision","confidential_remarks":"The core issue for me is the SHG methodology: EPW is not a code that computes chi^(2) on its own, and no formula or input details are given. This may reflect a misunderstanding of the software stack. I would ask the editor to require the authors to specify precisely how chi^(2) was obtained, including the working formula and the actual input files, and to add benchmark tests. If that cannot be provided, the SHG claims should be withdrawn. I also recommend requesting a correction of the monotonic bandgap claim and the Berry curvature symmetry relations before the paper is reconsidered."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Main thing you should know: the bandgap-opening part of this paper is probably fine, and the Berry curvature contrast is plausible, but the second-harmonic generation pillar is unsupported—the methods text attributes chi^(2) extraction to EPW, which is an electron-phonon package, and gives no formula, no Wannier parameters, no convergence tests. As it stands, the SHG results are unverifiable. That is not a nitpick; it is a load-bearing part of the abstract and conclusions.\n\nWhat is new: previous work (Denis, ref 24) reported bandgap opening for Ge-doped graphene but didn't compute Berry curvature, Chern numbers, or SHG. This paper adds a systematic scan of four concentrations with those properties, which is a legitimate extension. The bandgap opening itself is consistent with earlier work, and the Berry curvature maps show opposite signs at K and K' which is what you'd expect when inversion symmetry is broken. Table 1 does give a rough sense of gap sizes, though the text's claim that the gap increases with concentration is contradicted by the table: 5.5% gives 0.168 eV (HSE polarized) while 3.125% gives 0.26 eV and 2% gives 0.168 eV. Non-monotonic is not automatically wrong—different supercells, strain, and band-folding effects—but the authors don't acknowledge it.\n\nSoft spots, in order of importance. (1) SHG as described is a black box. EPW doesn't compute second-order susceptibilities; 'custom scripts' covering the whole derivation is not acceptable without at least citing the formula used and reporting Wannier/convergence details. An independent check is needed before that claim stands. (2) The title says 'Valley Quantum Hall Effect' but the paper computes neither a Hall conductivity nor quantized Chern numbers; the values in Table 2 are essentially zero, so the system is trivial. The phrase 'potential valley Hall effect' in the text is fine; the title overclaims. (3) Table 1 includes references [4] and [5] in the last column, which are about valleytronics in general, not about these bandgap values; those look like misplaced citation artifacts. (4) The paper never reports Wannier projectors or interpolation meshes for the Berry curvature, so the figures can't be reproduced as is.\n\nWho it's for: people working in 2D valleytronics and doping strategies. If the authors redo or at least fully document the SHG calculation, retitle the paper, and correct the bandgap trend discussion, this could be a modest but useful contribution. In current form, I would not cite it, and I'd send it back for major revision rather than reject out of hand—there's enough real DFT data here that a serious referee could help the authors fix it.\n\nMy recommendation: send to peer review with a clear request for the SHG details and a corrected title. The core band structure results look defensible.","headline":"A plausible bandgap-engineering study that overreaches with an unsubstantiated SHG claim and a misleading title.","tokens_in":10877,"tokens_out":3260,"would_cite":false,"duration_ms":37083,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Substitutional germanium in graphene breaks inversion symmetry and produces valley-antisymmetric Berry curvature and second harmonic generation, making the doped monolayer a candidate for valleytronic and nonlinear optoelectronic devices.","keywords":["Valleytronics","Valley Quantum Hall Effect","Berry Curvature","Second Harmonic Generation","Inversion Symmetry Breaking","germanium-doped graphene","bandgap engineering","density functional theory"],"falsifier":"Recompute the Berry curvature and SHG coefficients with explicitly converged Wannier interpolation on a fine k-grid (for example 24×24×1 or denser) with defined projectors and reported disentanglement windows; if the antisymmetric $\\Omega_z$ peaks at $K$ and $K'$, or the main $\\chi^{(2)}$ peak near 1.1 eV for the 2% Ge cell, shifts by more than a few percent when the grid is doubled, the reported valley Hall and nonlinear signatures are numerical artifacts rather than intrinsic properties of Ge-doped graphene.","tokens_in":9951,"feed_emoji":"⚛️","tokens_out":9420,"duration_ms":103777,"temperature":0.7,"pith_summary":"The paper claims that substitutional germanium doping converts gapless graphene into a bandgap semiconductor and activates the valley degree of freedom: at the $K$ and $K'$ valleys the Berry curvature takes opposite signs, which would drive carriers from the two valleys to opposite edges of a sample (a valley Hall effect) while time-reversal symmetry remains intact. It further claims that the bandgap, the valley-contrasting curvature, and second harmonic generation can all be tuned by changing the Ge concentration (12.5%, 5.5%, 3.125%, 2%), with gaps reaching about 1.06 eV and the sharpest nonlinear response near 1.1 eV at 2% doping. If these results hold, a simple one-atom-per-supercell doping recipe would give the most-studied 2D material both a valleytronic handle and a nonlinear optical response.","feed_headline":"Ge doping opens graphene's gap and activates valley transport","feed_subtitle":"A single substitution per supercell breaks inversion symmetry and tunes gap, valley Hall, and nonlinear response.","key_machinery":"The load-bearing object is the Berry curvature $\\Omega_n(\\mathbf{k})=\\nabla_\\mathbf{k}\\times A_n(\\mathbf{k})$ computed from the Bloch states, understood as an effective magnetic field in momentum space. Doping with Ge realizes an inversion-symmetry-breaking mass term $\\hbar m\\sigma_z$ in the Dirac Hamiltonian $H(k)=\\pm\\hbar v_F(\\sigma_x k_x+\\sigma_y k_y)+\\hbar m\\sigma_z$, which is exactly the term that allows $\\Omega_n$ to be nonzero and opposite at $K$ and $K'$; its sign pattern yields the valley Hall velocity $v_h=(e/\\hbar)E\\times\\Omega_n$, while the same broken inversion symmetry permits a second-order nonlinear susceptibility $\\chi^{(2)}$ and hence second harmonic generation. The near-zero Chern numbers then follow from time-reversal symmetry, leaving a trivial topological phase despite nontrivial valley physics.","core_discovery":"On the paper's own terms, the central discovery is that replacing one carbon atom in a graphene supercell with germanium breaks spatial inversion symmetry while preserving time-reversal symmetry, and that this symmetry splitting is enough to produce a bandgap, a nonzero and valley-antisymmetric Berry curvature $\\Omega_z(K)=-\\Omega_z(K')$, and a nonzero second-order susceptibility $\\chi^{(2)}$ that tracks optical absorption. The computed Chern numbers are near zero, so the system remains topologically trivial even though it shows valley-contrasting transport, and the authors read this combination as evidence of a potential (not quantized) valley Hall effect and of efficient second harmonic generation. The 12.5% and 3.125% doped supercells retain clear valley contrasts, while the 5.5% cell is reported as anomalous, with zero unpolarized gap and no valley polarization; the authors present this as part of the concentration dependence.","pith_inferences":["A direct convergence test with denser k-meshes and explicit Wannier projectors would settle whether the antisymmetric Berry-curvature peaks and the $\\chi^{(2)}$ peaks are intrinsic properties or interpolation artifacts; the near-zero Chern numbers make this check especially informative because small numerical noise could mimic or erase them.","The reported concentration dependence is not monotonic: 5.5% Ge gives zero unpolarized gap while 3.125% gives 0.26 eV in the polarized hybrid calculation, so if confirmed, local bonding geometry around the dopant, not just doping fraction, controls the valley physics; formation-energy or strain-field calculations would identify the controlling factor.","The same inversion-breaking mechanism should operate for other group-IV substituents such as Si and Sn, with the spin-orbit contribution to valley splitting expected to grow with atomic number, so a Ge-versus-Sn comparison would isolate the role of SOC in the valley response.","Because the paper does not compute valley-resolved optical selection rules, its suggestion that circularly polarized light can read out the valley polarization remains untested; a direct calculation of left- versus right-circular absorption at the $K$ and $K'$ transitions would close that gap."],"forward_implications":["Ge doping converts graphene into a semiconductor whose gap depends on dopant concentration, with 0.17 eV at 2% Ge and 1.06 eV at 12.5% Ge (HSE06), while the 5.5% cell behaves anomalously with zero unpolarized gap.","Carriers from opposite valleys acquire opposite Hall velocities, so a lateral electric field should accumulate $K$-valley carriers on one edge and $K'$-valley carriers on the other, giving a device-relevant valley Hall response.","The material should emit at $2\\omega$ when pumped at $\\omega$ near the gap, with the 2% Ge sample showing a clear second-harmonic feature around 1.1 eV that matches an absorption peak.","Because the Chern number is zero for every concentration, the predicted valley Hall effect is not a quantized anomalous Hall effect; any edge accumulation would be ordinary diffusive valley transport rather than topologically protected current.","Spin-orbit coupling from Ge lifts the orbital degeneracy of the bands, but the computed spin splitting is too small to select a single valley with circularly polarized light in the present doping range."],"supporting_citations":[{"why":"Prior finding that Ge, Ga, As, and Se substitution opens a 0.3-1.3 eV bandgap in graphene, the starting point the paper extends to valley and nonlinear response.","marker":"[24]"},{"why":"Introduces the valley filter and valley valve concept in graphene, which motivates the valley Hall interpretation used here.","marker":"[4]"},{"why":"Provides the valley-contrasting Berry curvature and topological transport formalism used to interpret the $K$/$K'$ asymmetry.","marker":"[5]"},{"why":"Underlies all structural, band-structure, and density-of-states calculations through the plane-wave DFT package used in the computational workflow.","marker":"[25]"},{"why":"First-principles analysis of spin-orbit coupling band-structure topologies in graphene, used to interpret the Ge-induced orbital degeneracy lifting.","marker":"[26]"},{"why":"Proposes second-harmonic spectroscopy to optically detect valley polarization, justifying the SHG calculations as a valley probe.","marker":"[27]"}],"fun_headline_variants":["Ge doping opens graphene bandgap and activates valley Hall effect","Broken symmetry from Ge dopant yields graphene valley Hall and SHG","Ge doping concentration tunes graphene's bandgap and valley polarization","One Ge atom per cell flips graphene into valley Hall material","Germanium doping breaks symmetry in graphene for valleytronics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire valley Hall and SHG picture depends on the accuracy of the numerical interpolation step that connects the plane-wave electronic states to the fine grids used for Berry curvature and nonlinear optical spectra, but the paper reports no projectors, disentanglement windows, interpolation grids, or convergence tests for that step; if it is unconverged, the reported $K$/$K'$ contrast and $\\chi^{(2)}$ peaks could be numerical artifacts rather than intrinsic material properties.","fun_headline_variants_meta":{"raw":{"variants":["Ge doping opens graphene bandgap and activates valley Hall effect","Broken symmetry from Ge dopant yields graphene valley Hall and SHG","Ge doping concentration tunes graphene's bandgap and valley polarization","One Ge atom per cell flips graphene into valley Hall material","Germanium doping breaks symmetry in graphene for valleytronics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000788,"raw_usage":{"total_tokens":3495,"prompt_tokens":982,"completion_tokens":2513,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":598,"completion_tokens_details":{"reasoning_tokens":2428}},"tokens_in":598,"tokens_out":2513,"duration_ms":20705,"temperature":1.0,"reasoning_tokens":2428,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:26:55.271292+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the Berry curvature and SHG coefficients with explicitly converged Wannier interpolation on a fine k-grid (for example 24×24×1 or denser) with defined projectors and reported disentanglement windows; if the antisymmetric $\\Omega_z$ peaks at $K$ and $K'$, or the main $\\chi^{(2)}$ peak near 1.1 eV for the 2% Ge cell, shifts by more than a few percent when the grid is doubled, the reported valley Hall and nonlinear signatures are numerical artifacts rather than intrinsic properties of Ge-doped graphene.","supporting_citations":[{"cited_title":"Giant and Controllable Valley Currents in Graphene by Double Pumped THz Light","cited_arxiv_id":null,"evidence_quote":"Prior finding that Ge, Ga, As, and Se substitution opens a 0.3-1.3 eV bandgap in graphene, the starting point the paper extends to valley and nonlinear response."},{"cited_title":"The reported results were calculated by employing the first -principle DFT methodology based on GGA-PBE potential and hybrid HSE06 functional with and without spin -orbit coupling","cited_arxiv_id":null,"evidence_quote":"Introduces the valley filter and valley valve concept in graphene, which motivates the valley Hall interpretation used here."},{"cited_title":"However, these results confirm that the doping of Ge in graphene facilitates to open a considerable bandgap","cited_arxiv_id":null,"evidence_quote":"Provides the valley-contrasting Berry curvature and topological transport formalism used to interpret the $K$/$K'$ asymmetry."},{"cited_title":"Advances and Trends in Chemically Doped Graphene","cited_arxiv_id":null,"evidence_quote":"Underlies all structural, band-structure, and density-of-states calculations through the plane-wave DFT package used in the computational workflow."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"First-principles analysis of spin-orbit coupling band-structure topologies in graphene, used to interpret the Ge-induced orbital degeneracy lifting."},{"cited_title":"QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials,","cited_arxiv_id":null,"evidence_quote":"Proposes second-harmonic spectroscopy to optically detect valley polarization, justifying the SHG calculations as a valley probe."}],"review_version":1}