{"id":"b13dca76-73c9-4fa1-aa6d-449174bd9788","arxiv_id":"2411.13429","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"The authors derive exact transition amplitudes for a graphene Coulomb impurity driven by an ultrarelativistic ion and predict 2D ionization cross sections roughly 100 times larger than the 3D counterpart.","lead":"This paper calculates the probability that a fast, nearly light-speed ion passing over a graphene sheet excites or ionizes an artificial atom formed by a charged impurity. It reports ionization cross sections about two orders of magnitude larger than for ordinary three-dimensional atoms, which could make these table-top 2D 'relativistic atoms' observable.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Most load-bearing concern: Eq. (4) imports the 3D ultrarelativistic coupling (1-σz), but graphene quasiparticles have vF≈c/300, so minimal coupling gives an almost scalar potential; the exact amplitude and Table I cross sections are therefore not yet predictions for graphene.","rationale":"The reader's weakest-assumption analysis identifies exactly the load-bearing point: the form of the projectile potential in Eq. (4) is imported from 3D ultrarelativistic QED and is not justified for graphene quasiparticles with vF << c. My independent check of the derivation confirms that the same assumption underlies the light-cone integration: the projector (1-σz) makes the δ(z-ct) potential exactly integrable only when the quasiparticle velocity equals c, and the final amplitude Eq. (13) inherits this projector in the matrix element. For graphene, the correct minimal-coupling vertex is e(1 - (vF/c)σz), which is nearly scalar, so the transition amplitudes, probabilities, and Table I cross sections would change substantially. This is not merely a disagreement with an external convention; it is an internal mismatch between the model Hamiltonian H0, which explicitly contains vF≈c/300, and the perturbation V, which implicitly assumes vF=c. Therefore the exact-solution claim and the central quantitative prediction are not established for graphene. The reader's REJECT verdict is reasonable; my stress-test does not move the verdict. I also note the additional secondary issue that the comparison with the 3D cross section of Ref. [16] is not fully parameter-matched, but that is secondary to the incorrect effective potential. A concrete numerical recomputation with the scalar potential would settle whether the two-orders-of-magnitude claim survives or is an artifact of the (1-σz) projector.","tokens_in":8296,"tokens_out":9885,"duration_ms":117434,"concrete_test":"Re-derive the ion potential from minimal coupling to the boosted Coulomb field: V_eff = -δ(z-ct) α ZP [1 - (vF/c) σz] ln(1+x^2/b^2), using the vF in Eq. (1). Then repeat the derivation of Eq. (13) without the projector step, i.e., solve the jump condition across z=ct with the full matrix c - vF σz, and recompute the ionization cross section (17) for the Table I parameters (ZP=10, ZT αgr=0.4, M≈0.1 eV, vF=c/300). If the resulting cross section differs from 2.74×10^9 barn by a significant factor, or if the two-orders-of-magnitude enhancement over the 3D result disappears, then Eq. (13) is not a graphene prediction and the central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result, Eq. (13), is derived from the potential in Eq. (4), V(x,z,t) = -δ(z-ct) α ZP (1-σz) ln(1+x^2/b^2). This form is taken over from 3D ultrarelativistic ion-atom collisions, where the Dirac velocity operator along the beam is c αz and the combination e(A0 - c αz Az) becomes e(1-αz)A0 for A0≈Az. In graphene, the Hamiltonian (1) has velocity operator vF σz with vF≈c/300. Minimal coupling to the projectile's Liénard-Wiechert potentials gives V = eφ - e vF σ·A ≈ -δ(z-ct) α ZP [1 - (vF/c) σz] ln(1+x^2/b^2). Since vF/c ≈ 1/300, this is essentially the identity in sublattice space, not (1-σz). The (1-σz) factor in Eq. (13) is thus a model assumption, not a consequence of the graphene Dirac equation. Moreover, the light-cone integration leading to Eq. (11) uses the projector (1-σz) to exponentiate the δ(z-ct) potential; this step is exact only when the quasiparticle velocity along z equals c. For vF << c the jump condition across the light-cone delta contains c - vF σz rather than c(1-σz), so the claimed exact solution is exact for a different, artificially ultrarelativistic 2D Dirac fermion. Since the headline 'two orders of magnitude larger' ionization cross sections in Table I are computed from this amplitude, the central quantitative claim is not supported for graphene unless Eq. (4) is independently justified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a theoretical model for electronic transitions in a Coulomb impurity in graphene (a '2D artificial atom') driven by a fast ultrarelativistic ion. The authors write down a time-dependent 2D Dirac equation with a projectile potential taken as a light-cone delta function with a specific sublattice operator (1-σ_z). They solve the equation exactly using light-cone coordinates and derive transition amplitudes, survival/excitation/ionization probabilities, and ionization cross sections. The central quantitative claim is that the ionization cross section is two orders of magnitude larger than the corresponding 3D result.","tokens_in":8662,"tokens_out":10021,"duration_ms":106412,"significance":"If the model were correct, this would provide a table-top probe of 2D relativistic atomic physics and strong-field collision phenomena in graphene, with a clean analytic exact solution. The paper is internally consistent mathematically, and the use of unitarity to obtain ionization probabilities without continuum wavefunctions is elegant. However, the physical input potential is not derived from the graphene Dirac Hamiltonian, and the predictive power of the results for graphene is therefore not established.","major_comments":[{"comment":"The projectile potential in Eq. (4) is taken from the 3D ultrarelativistic calculation of Ref. [23], where the operator (1-α_z) appears because the electron's Dirac velocity along the beam is c α_z. For the 2D Dirac-Weyl Hamiltonian in Eq. (1), the velocity operator along the beam direction is v_F σ_z with v_F ≈ c/300. Minimal coupling to the Liénard-Wiechert potentials of the ion, with A_z ≈ A_0 in the ultrarelativistic limit, gives V = e φ σ_0 - e v_F σ_z A_z ≈ e φ (σ_0 - (v_F/c) σ_z), which is essentially proportional to σ_0, not to (1-σ_z). The paper does not derive Eq. (4) from the graphene Hamiltonian; it merely states 'Following the arguments in Ref. [23]'. Since the exact solution in Eq. (11) and the transition amplitude in Eq. (13) rely on the projector property of (1-σ_z), the replacement of (1-σ_z) by the correct near-scalar coupling would invalidate the exact solution. Consequently, the ionization cross sections in Table I are not predictions for graphene unless Eq. (4) is independently justified.","section":"Section 2, Eq. (4)"},{"comment":"The light-cone integration that produces Eq. (11) implicitly assumes that the quasiparticle propagates with velocity c in the z direction. For the graphene Hamiltonian (1), the characteristic velocity along z is v_F, so the jump condition across the moving delta function δ(z-ct) contains c - v_F σ_z rather than c(1-σ_z). Thus the derivation is internally consistent only for an artificially ultrarelativistic 2D Dirac fermion with v_F = c, not for graphene. This is a separate reason that the claimed exact solution does not apply to the stated physical system.","section":"Section 2, Eqs. (10)-(11)"}],"minor_comments":[{"comment":"The delta-function potential in Eq. (4) has implicit dimensions; please specify the natural units used and the origin of the dimensional factor in the coefficient so that the expression is dimensionally transparent.","section":"Section 2, Eq. (4)"},{"comment":"In Fig. 4 and the sum in Eq. (16), specify how many bound states are included in the numerical sum and demonstrate convergence of the result with respect to the number of states.","section":"Section 3, Eqs. (14)-(16)"},{"comment":"Provide the 3D cross-section values from Ref. [16] alongside the 2D values so that the 'two orders of magnitude larger' claim can be verified quantitatively.","section":"Table I"},{"comment":"There are several typos, for example 'formed offers' and 'extenstion' in the Conclusions; please proofread the manuscript carefully.","section":"Conclusions"},{"comment":"Reference [28] is missing the article title; please complete the bibliographic entry.","section":"References"},{"comment":"The caption contains the typo 'artifical' instead of 'artificial'; this should be corrected.","section":"Fig. 1 caption"}],"recommendation":"reject","confidential_remarks":"The paper is clearly written and the algebra is internally consistent, but the central physical input (the form of the projectile potential) is not appropriate for graphene. In my view, this is not a fixable local issue: the exact solution methodology relies on the (1-σ_z) projector, which is absent in the correct minimal-coupling potential for graphene. The authors would need to redo the calculation with a substantially different approach. I recommend rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper takes Baltz's exact ultrarelativistic ion–atom collision solution and transplants it to a gapped graphene Coulomb impurity. The algebra from Eq. (4) onward is clean, the use of Novikov's bound states is appropriate, and the combination is genuinely new. There is no fitting to data; the transition amplitudes are derived analytically from the stated model. Credit where earned: this is a serious, competently executed adaptation of a known method, and the survival, excitation, and ionization probabilities in Figs. 2–4 follow coherently from the model.\n\nThe soft spot is load-bearing. The potential in Eq. (4), V = −δ(z−ct) α ZP (1−σz) ln(1+x²/b²), is taken directly from the 3D case, where the (1−αz) factor arises because the electron's velocity along the beam is c. In graphene, the quasiparticle velocity is vF σz with vF ≈ c/300. Minimal coupling to the ion's Liénard–Wiechert fields gives essentially a scalar potential in sublattice space, not (1−σz). The exact solution in Eq. (13) is therefore exact for a different, artificially ultrarelativistic 2D Dirac fermion. The jump condition across the delta potential would contain c − vF σz, not c(1−σz), so the derivation does not carry over. Consequently, the Table I cross sections, including the 'two orders of magnitude larger than 3D' claim, are not predictions for graphene as it stands. The paper does not derive or justify the (1−σz) form for the 2D case, and that is not a minor omission.\n\nA secondary, lesser issue: the comparison with the 3D cross sections is not parameter-matched (different gap, different masses, different couplings), so the headline enhancement is not a controlled comparison. That is minor relative to the potential problem.\n\nWho is this for? A reader interested in exact methods for time-dependent Dirac problems in 2D might find the technique useful, but the physical conclusions should not be taken at face value. The paper is not incoherent—the flaw is a modeling assumption, not an internal contradiction—but the central quantitative claim is unsupported until the potential is rederived from minimal coupling. I would send it to a referee, because the method and the question are worth engaging and the issue is fixable. My own recommendation as a referee would be major revision: derive the correct potential, redo the calculation, and either show the (1−σz) form is justified at some level of approximation or retract the graphene-specific predictions.","headline":"A clean, internally consistent application of Baltz's exact collision solution to a gapped graphene Coulomb impurity, but the load-bearing (1−σz) potential is imported from 3D ultrarelativistic physics and is not justified for graphene quasiparticles with vF ≈ c/300.","tokens_in":9217,"tokens_out":1648,"would_cite":false,"duration_ms":19566,"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":"This paper claims that a fast ultrarelativistic ion sweeping past a graphene layer ionizes a Coulomb-impurity artificial atom with a cross section two orders of magnitude larger than the equivalent three-dimensional process, and it…","keywords":["graphene","Coulomb impurity","artificial atom","Dirac equation","ionization cross section","ultrarelativistic ion","exact solution","2D Dirac materials"],"falsifier":"Compare the predicted impact-parameter-resolved ionization probability $P(b)$ and the integrated cross sections (for example, about $2.7\\times 10^9$ barn for $Z_P=10$) against a controlled experiment in which a fast ion beam of known charge passes over a graphene sample with isolated Coulomb impurities; a measured cross section differing by orders of magnitude, or a survival probability that does not saturate to unity at large $b$, would rule out the assumed potential.","tokens_in":8070,"feed_emoji":"⚛️","tokens_out":5066,"duration_ms":50871,"temperature":0.7,"pith_summary":"The paper treats a charged impurity (adatom) in a graphene monolayer as a two-dimensional artificial atom and asks what happens when an ultrarelativistic ion passes parallel to the layer at impact parameter b. The authors solve the time-dependent 2D Dirac equation exactly for this setup, exploiting a light-cone coordinate trick that turns the ion's brief pulse into a simple phase factor. From that solution they obtain the survival probability of the ground state, excitation probabilities to bound states, and the total ionization probability via unitarity. Their central quantitative claim is that the ionization cross section is two orders of magnitude larger than the corresponding three-dimensional relativistic ion-atom collision. If correct, this makes collision-induced transitions in graphene Coulomb impurities observable with a table-top ion beam, offering a laboratory window into 2D relativistic atomic physics.","feed_headline":"Graphene impurity atoms ionize 100x more readily","feed_subtitle":"An exact 2D Dirac solution predicts huge ionization cross sections from a fast ion beam passing over graphene.","key_machinery":"The machinery is the reduction of the scattering problem to a light-cone instantaneous kick. The projectile potential is idealized as $V(x,z,t) = -\\delta(z-ct)\\,\\alpha Z_P (1-\\sigma_z)\\ln(1+x^2/b^2)$, a delta-function pulse along the ion's worldline with a sublattice projection. Integrating the Dirac equation across the light-cone coordinate $z_-$ yields the boundary condition $(1-\\sigma_z)\\Psi = (1-\\sigma_z) e^{-i\\theta(ct-z)\\alpha Z_P \\ln(1+x^2/b^2)} \\psi_i$, and substituting this into the interaction-picture equation gives the exact amplitude as a static matrix element. All probabilities follow from that one expression, and ionization is obtained without computing continuum wavefunctions by subtracting survival and excitation from unity.","core_discovery":"The central discovery is an exact closed-form transition amplitude for a Coulomb-impurity artificial atom in graphene struck by a fast ion. Working in light-cone coordinates, the authors reduce the time-dependent Dirac equation to a constraint that the spinor just before the pulse equals the initial spinor multiplied by a phase $e^{-i\\alpha Z_P \\ln(1+x^2/b^2)}$ on the sublattice selected by $(1-\\sigma_z)$. The amplitude from initial state $\\psi_i$ to final state $\\psi_f$ is then a single matrix element, Eq. (13), whose modulus squared gives survival and excitation probabilities; ionization follows from unitarity as $1 - P_{\\text{survival}} - P_{\\text{excitation}}$. Numerically, for parameters $\\alpha_{\\rm gr} Z_T = 0.4$, $M\\approx 0.1$ eV, the ionization cross sections range from $2.76\\times 10^7$ barn for $Z_P=1$ to $2.74\\times 10^9$ barn for $Z_P=10$, about two orders of magnitude above the 3D results from the corresponding exact Dirac calculation.","pith_inferences":["If the $(1-\\sigma_z)$ factor is replaced by a sublattice-diagonal coupling, the exact solution is lost and the cross sections may change by orders of magnitude; a first-principles derivation of the graphene projectile potential would settle whether the 'two orders of magnitude' advantage survives.","The same light-cone technique should apply to other two-dimensional Dirac materials with different Fermi velocities, with cross sections expected to scale with the effective fine-structure constant and with $v_F/c$, offering a tunable family of table-top experiments.","The predicted dominance of ionization over bound excitation suggests that charge-sensing or current-noise measurements in a graphene flake could detect individual ionization events, providing a direct experimental test."],"forward_implications":["A fast ion beam passing over a graphene layer should ionize Coulomb-impurity artificial atoms with high efficiency: cross sections of order $10^7$ to $10^9$ barn for projectile charge $Z_P$ between 1 and 10.","Ionization dominates over excitation at large impact parameters, so the dominant observable signature is a free electron rather than a bound excited state.","The unitarity relation in Eq. (16) lets one compute total ionization without explicit continuum wavefunctions, a shortcut that extends to any transition sum.","The exact amplitude in Eq. (13) can be reused to compute particle-hole pair creation and angular distributions of differential cross sections in the same setup."],"supporting_citations":[{"why":"Provides the 3D exact Dirac-equation ionization calculation whose cross sections the paper compares against, anchoring the two-orders-of-magnitude claim.","marker":"[16]"},{"why":"Supplies the exact Dirac-equation treatment of ultrarelativistic heavy-ion ionization and pair production that the paper adapts to 2D; the light-cone method and the form of the potential come from this line of work.","marker":"[23]"},{"why":"Gives the closed-form bound-state and scattering wave functions of the Coulomb impurity in graphene that the paper uses as basis states.","marker":"[27]"},{"why":"Provides the vacuum (supercritical) Dirac spectrum for graphene that underlies the target model and the discussion of pair creation.","marker":"[28]"},{"why":"Establishes the low-energy Dirac description of graphene and its effective fine-structure constant, the starting Hamiltonian of the paper.","marker":"[1]"},{"why":"Derives the Coulomb potential of a particle in uniform ultrarelativistic motion, the source of the projectile potential form.","marker":"[22]"}],"fun_headline_variants":["Exact Dirac solution: graphene impurity ionization 100x higher","Fast ions increase graphene impurity ionization by 100x","Exact 2D Dirac solution gives 100x graphene ionization","Graphene impurity atoms ionize 100x via fast ion beam","Coulomb impurity in graphene: 100x ionization cross-section"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire calculation rests on the assumed shape of the ion's potential, a delta-pulse with a logarithmic profile that acts only on one sublattice; if the real graphene coupling differs, the exact solution no longer applies to graphene.","fun_headline_variants_meta":{"raw":{"variants":["Exact Dirac solution: graphene impurity ionization 100x higher","Fast ions increase graphene impurity ionization by 100x","Exact 2D Dirac solution gives 100x graphene ionization","Graphene impurity atoms ionize 100x via fast ion beam","Coulomb impurity in graphene: 100x ionization cross-section"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000688,"raw_usage":{"total_tokens":3073,"prompt_tokens":857,"completion_tokens":2216,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":473,"completion_tokens_details":{"reasoning_tokens":2129}},"tokens_in":473,"tokens_out":2216,"duration_ms":19775,"temperature":1.0,"reasoning_tokens":2129,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:26:59.873616+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the predicted impact-parameter-resolved ionization probability $P(b)$ and the integrated cross sections (for example, about $2.7\\times 10^9$ barn for $Z_P=10$) against a controlled experiment in which a fast ion beam of known charge passes over a graphene sample with isolated Coulomb impurities; a measured cross section differing by orders of magnitude, or a survival probability that does not saturate to unity at large $b$, would rule out the assumed potential.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the 3D exact Dirac-equation ionization calculation whose cross sections the paper compares against, anchoring the two-orders-of-magnitude claim."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the exact Dirac-equation treatment of ultrarelativistic heavy-ion ionization and pair production that the paper adapts to 2D; the light-cone method and the form of the potential come from this line of work."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the closed-form bound-state and scattering wave functions of the Coulomb impurity in graphene that the paper uses as basis states."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives the Coulomb potential of a particle in uniform ultrarelativistic motion, the source of the projectile potential form."}],"review_version":1}