{"id":"cde85124-fec1-4e43-b575-8aa95e1d5e9f","arxiv_id":"2607.04646","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"A non-Hermitian Dirac-vortex model with complex-mass winding and infinite imaginary potential yields closed-form TCSEL frequencies, thresholds, and tunable vector-beam polarizations, confirmed experimentally.","lead":"The authors solve a non-Hermitian Dirac vortex with complex mass and an absorbing boundary, yielding closed-form laser frequencies, thresholds, and vector-beam polarizations for topological-cavity surface-emitting lasers. The formulas match their experiments and give a design rule for single-mode TCSELs.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"Continuum Dirac + infinite-imaginary barrier + μ-perturbation hierarchy is the softest point for the claimed closed-form TCSEL design formulas.","rationale":"The Reader correctly isolates the continuum/infinite-barrier/μ-perturbation hierarchy as the weakest load-bearing premise; the supplied manuscript supplies only scale-separation arguments and visual far-field agreement (Figs. 6–7), with no quantitative error bars or direct comparison of predicted versus measured losses. No stronger internal inconsistency appears in the analytic structure (Jackiw–Rossi extension, Whittaker radial solutions, PT pairing, asymptotic scalings). Because that single concern already justifies CONDITIONAL rather than unconditional ACCEPT, and because the concrete numerical check above would settle it, the Reader’s verdict and confidence need no revision.","tokens_in":13525,"tokens_out":604,"duration_ms":18036,"concrete_test":"Run a 3-D CWT or full-wave simulation of a realistic triangular-lattice TCSEL (finite absorption 600 cm^{-1} outside a circular pumped region) at mR=1.5 and 2.0; extract the complex eigenvalues of the zero mode, unbound singlet and bound doublet. Check whether (i) zero-mode α∥R matches the asymptotic 16(mR)^{2}e^{-4mR} within 20 %, (ii) α⊥ remains mode-independent to better than the analytic ΔαR, and (iii) the ordering and ΔαR agree with Fig. 3. Failure of any check by more than the design margin falsifies quantitative reliability of the closed forms.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (closed-form ω, α∥, α⊥ and θ0-tunable vector beams from the non-Hermitian Dirac vortex capture essential TCSEL single-mode physics) rests on three linked approximations that are only scale-justified, not error-bounded: (1) continuum Dirac description of a discrete C3 photonic-crystal slab, (2) replacement of the finite-absorption exterior (~600 cm^{-1}) by an infinite imaginary potential at r=R (justified by 600 ≫ 1/R ~20 cm^{-1}), and (3) solving the Hermitian vortex first then treating the imaginary-mass winding μ≪m purely perturbatively so that α⊥=2μ is identical for every mode (Eq. 5) and the radiation operator (Eq. 8) is accurate. If lattice-scale corrections, soft-boundary leakage, or residual μ-induced mode mixing alter the relative boundary losses by an amount comparable to the predicted ΔαR peak near mR≈1.5 (Fig. 3), the analytic mode ordering, threshold margins and single-mode stability window cease to be quantitatively predictive even while remaining qualitatively useful.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript constructs a non-Hermitian Dirac-vortex model that combines a complex-mass winding (real mass m with winding w and imaginary mass μ with winding −2w) with an infinite-imaginary-potential boundary of radius R. This extends the Jackiw–Rossi vortex and neutrino-billiard models into the dissipative regime and is presented as a minimal continuum theory for topological-cavity surface-emitting lasers (TCSELs). From the 4\times4 Hamiltonian (Eq. 1) and the absorbing boundary condition (Eq. 2), the authors obtain closed-form modal frequencies and boundary losses (via Whittaker radial functions after neglecting μ for the eigenvalue problem), show that vertical radiation loss α⊥ = 2μ is mode-independent to first order (Eq. 5), and derive the radiation operator that yields θ0-tunable cylindrical vector beams (Eq. 8). Asymptotic scalings of α∥R (zero-mode ∼16(mR)^{2}e^{-4mR}, bound doublet exponentially suppressed, unbound singlet power-law) and the single-mode stability window near mR ≈ 1.5 are given in Fig. 3. End Matter solves the massless non-Hermitian billiard and reports experimental far-field polarization control and the zero-mode/unbound-singlet crossover.","tokens_in":13799,"tokens_out":1295,"duration_ms":9628,"significance":"If the continuum-plus-perturbation hierarchy holds at the claimed level of accuracy, the work supplies a rare, analytically closed non-Hermitian topological design theory for a practical large-area laser. Closed-form loss scalings, the identification of an optimal mR window, and the explicit θ0 control of vector-beam polarization are directly usable for device engineering and go beyond purely numerical coupled-wave theory. The experimental matches for polarization patterns (Fig. 6) and the lasing-mode crossover (Fig. 7) provide independent, falsifiable checks. The construction also cleanly unifies Jackiw–Rossi, neutrino-billiard and non-Hermitian optics, which is of broader interest to topological photonics.","major_comments":[{"comment":"The central claim that the closed-form thresholds and mode ordering are quantitatively predictive for real TCSELs rests on three linked approximations that are only scale-justified: continuum Dirac description of a discrete C3 photonic-crystal slab, replacement of finite exterior absorption (~600 cm^{-1}) by an infinite imaginary potential (justified by 600 ≫ 1/R ~ 20 cm^{-1}), and first-order treatment of μ ≪ m so that α⊥ is identical for every mode (Eq. 5) and the radiation operator (Eq. 8) is accurate. No error bound or direct comparison of the analytic ΔαR peak near mR ≈ 1.5 (Fig. 3) against full CWT or finite-absorption numerics is provided. If lattice-scale or soft-boundary corrections reorder the lowest-loss modes by an amount comparable to that peak, the claimed single-mode stability window ceases to be quantitative. A short numerical validation (or an explicit statement of the e","section":null},{"comment":"The radiation operator (Eq. 8) and the assertion that α⊥ = 2μ is strictly mode-independent both rely on the same first-order perturbation in μ after the Hermitian vortex has been solved. Residual μ-induced mixing between the zero mode and the nearby unbound singlet/bound doublet is not estimated. Because the threshold margin ΔαR is itself of order the boundary-loss differences plotted in Fig. 3, even a modest second-order correction could shift the optimal operating point. A brief estimate of the size of these corrections (or a statement that they remain negligible throughout the recommended mR window) would strengthen the load-bearing claim.","section":null}],"minor_comments":[{"comment":"The factor of 1/2 that converts the complex eigenvalue into the intensity decay rate α is introduced without derivation in Eq. (1); a one-sentence reminder of the e^{iωt} convention would help non-specialist readers.","section":null},{"comment":"Fig. 3(a) black dashed curve (minor-to-major component ratio) is useful but its definition appears only in the caption; a short inline definition would improve readability.","section":null},{"comment":"End Matter Table I compares Hermitian and non-Hermitian billiards; the boundary-condition operators are written with slightly different conventions from the main-text Eq. (2). Aligning the notation would avoid confusion.","section":null},{"comment":"Several key experimental references (e.g., the original TCSEL papers) are self-citations; a brief pointer to independent experimental realizations of related Dirac-vortex cavities would help place the work in the broader literature.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is a natural and high-quality follow-up to the authors’ own experimental TCSEL series. The analytic advance is genuine and the experimental checks are clean. The only substantive risk is over-claiming quantitative design accuracy without an error bar on the continuum/μ-perturbation hierarchy; once that is addressed (even with a short numerical appendix), the paper is suitable for a high-impact optics/photonics venue. No novelty or citation-pattern concerns beyond the expected self-citation of prior device work."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This paper is the analytical design theory the TCSEL community has been missing. They take the Jackiw–Rossi vortex, put complex mass windings on it, and close it with an infinite imaginary potential; the resulting Whittaker/Bessel spectrum, asymptotic losses (zero-mode α∥R ∼ 16(mR)2 e−4mR, unbound-singlet power-law, bound-doublet exponential), and the radiation operator that produces θ0-tunable +1/−2 vector beams are all new and closed-form. That is real progress over pure CWT numerics for devices thousands of lattice periods across.\n\nWhat they do well is keep the hierarchy clean: solve the Hermitian vortex first (μ ≪ m), get identical vertical loss α⊥ = 2μ by first-order perturbation, then read boundary loss and far-field from the minor spinor components induced by the absorbing wall. Anti-PT pairing, the mR > 1/4 zero-mode threshold, and the ΔαR peak near mR ≈ 1.5 all fall out cleanly. The End Matter massless “neutrino-billiard” limit and the experimental polarization rotation plus unbound-singlet crossover (Figs. 6–7) are honest checks, not fits.\n\nThe soft spot is exactly the continuum + infinite-absorption + μ-perturbation stack. They justify it by scale separation (600 cm−1 ≫ 1/R ∼ 20 cm−1) and μ ≪ m, but they never bound the error on relative losses near the design window mR ≈ 1.5. If lattice-scale or soft-boundary corrections scramble the ordering by an amount comparable to ΔαR, the analytic thresholds become qualitative only. That is a real limitation for quantitative engineering, not a fatal hole in the math. Self-citation of their own TCSEL papers is heavy but expected; the new analytic content stands on its own.\n\nThis is for people who design large-area single-mode surface emitters or who care about solvable non-Hermitian Dirac models. It deserves a serious referee. I would send it out.","headline":"Closed-form non-Hermitian Dirac-vortex formulas that actually give usable TCSEL design rules and match the polarization and mode-crossover experiments.","tokens_in":14459,"tokens_out":534,"would_cite":true,"duration_ms":5008,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"A solvable non-Hermitian Dirac vortex gives closed-form frequencies, thresholds, and polarizations for topological-cavity surface-emitting lasers.","keywords":["non-Hermitian Dirac vortex","topological-cavity surface-emitting laser","TCSEL","Jackiw-Rossi model","vector-beam polarization","complex-mass winding","absorbing boundary","photonic crystal laser"],"falsifier":"Fabricate a series of devices with fixed radius while continuously varying the mass amplitude so that mR sweeps through 0.85; if the lasing far-field does not switch from the three-lobe zero-mode pattern to the unbound-singlet pattern exactly as predicted, or if measured threshold margins deviate strongly from the analytic ΔαR curves, the continuum model fails.","tokens_in":14399,"feed_emoji":"🔆","tokens_out":1026,"duration_ms":7867,"temperature":0.7,"pith_summary":"Topological-cavity surface-emitting lasers (TCSELs) are large photonic-crystal devices that lase from a Dirac-vortex zero mode. Full-wave simulation of cavities thousands of periods across is impractical, so designers have relied on numerical coupled-wave theory. This paper constructs a minimal continuum model that extends the Jackiw–Rossi vortex and the neutrino-billiard boundary into the non-Hermitian regime: a complex-mass winding encodes vertical radiation loss while an infinite imaginary potential defines the absorbing edge of the active region. The model admits closed-form modal frequencies, boundary losses that scale differently for the zero mode and competing modes, and far-field vector-beam polarizations that can be rotated by the initial mass phase. Experiments on optically pumped devices confirm both the predicted polarization tuning and the crossover from zero-mode to unbound-singlet lasing when the normalized mass is reduced. The result supplies an analytical design tool that explains single-mode stability margins without heavy computation.","feed_headline":"Closed-form Dirac vortex solves TCSEL modes and polarizations","feed_subtitle":"Analytic losses and tunable vector beams match experiment, giving a design tool for large topological lasers","key_machinery":"Non-Hermitian Dirac vortex: the 4×4 continuum Hamiltonian with real mass m e^{i(wθ+θ0)}, imaginary mass −μ e^{−2i(wθ+θ0)}, and the absorbing boundary condition −τ_z(σ·n̂)|ψ⟩ = |ψ⟩ at r = R. Separation of variables plus first-order perturbation in μ produces the closed-form spectrum and radiation patterns.","core_discovery":"A non-Hermitian Dirac Hamiltonian with complex-mass winding number w and an infinite-imaginary-potential disk of radius R yields analytic eigenvalues and spinors. After nondimensionalization by R, the zero-mode boundary loss scales as 16(mR)^2 exp(−4mR), the next modes scale more slowly, and the radiation operator maps the spinor to two cylindrical vector beams of charges +1 and −2 whose relative weight and common polarization angle are fixed by mR and the initial mass phase θ0. These expressions reproduce the measured three-lobe far fields and the observed mode crossover.","pith_inferences":["Because the analytic thresholds depend only on the dimensionless product mR, the same formulas can be reused for electrically pumped or mid-infrared TCSELs once the effective mass is extracted from band structure.","The anti-PT pairing that isolates the zero mode may be portable to other non-Hermitian topological lasers that combine a mass defect with an absorbing boundary.","If lattice-scale corrections mix the valleys enough to spoil the shared vertical loss α⊥, the single-mode advantage would shrink even while the continuum topology remains intact—an effect testable by comparing devices with different supercell sizes."],"forward_implications":["Designers can choose mR ≈ 1.5 to maximize the analytic threshold margin while keeping free spectral range and boundary loss favorable.","Output polarization (radial, azimuthal, or spiral) is set by the single geometric parameter θ0 without redesigning the cavity shape.","The same closed-form loss hierarchy applies to any C3-symmetric lattice that realizes a Dirac vortex, independent of microscopic lattice details.","The massless non-Hermitian billiard spectrum supplies a universal lower bound αR ≥ 0.5 that any competing whispering-gallery mode must respect."],"fun_headline_variants":["Non-Hermitian Dirac vortex yields closed-form TCSEL modes","Analytic Dirac-vortex model solves TCSEL thresholds and beams","Minimal non-Hermitian Dirac theory captures TCSEL polarizations","Complex-mass winding gives exact TCSEL frequencies and vector beams","Closed-form Dirac eigenvalues match TCSEL losses and far fields"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The continuum Dirac description plus the infinite-absorption boundary remain accurate for a real photonic-crystal slab whose absorption is large but finite and whose lattice is discrete.","fun_headline_variants_meta":{"raw":{"variants":["Non-Hermitian Dirac vortex yields closed-form TCSEL modes","Analytic Dirac-vortex model solves TCSEL thresholds and beams","Minimal non-Hermitian Dirac theory captures TCSEL polarizations","Complex-mass winding gives exact TCSEL frequencies and vector beams","Closed-form Dirac eigenvalues match TCSEL losses and far fields"]},"model":"grok-4.5","effort":"low","cost_usd":0.00347,"raw_usage":{"total_tokens":1091,"prompt_tokens":721,"num_sources_used":0,"completion_tokens":90,"cost_in_usd_ticks":34700000,"prompt_tokens_details":{"text_tokens":721,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":280,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":721,"tokens_out":90,"duration_ms":2811,"temperature":1.0,"reasoning_tokens":280,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-11T15:50:35.805948+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Fabricate a series of devices with fixed radius while continuously varying the mass amplitude so that mR sweeps through 0.85; if the lasing far-field does not switch from the three-lobe zero-mode pattern to the unbound-singlet pattern exactly as predicted, or if measured threshold margins deviate strongly from the analytic ΔαR curves, the continuum model fails.","supporting_citations":[],"review_version":1}