{"id":"7e869823-a41d-4f63-bc5a-ecbfb8454c50","arxiv_id":"2507.00818","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Stable three-component vortex solitons, whose third-harmonic part is an anti-vortex with charge -1, are predicted in a dual-period quasi-phase-matched lithium-niobate photonic crystal.","lead":"This paper predicts a new class of stable composite vortex solitons: light beams carrying orbital angular momentum that survive long propagation in a lithium-niobate photonic crystal, thanks to a dual-period poling design that phase-matches both second- and third-harmonic generation. The third-harmonic part forms an anti-vortex, a phase pattern spinning opposite to the fundamental, which the authors propose as a new handle for vortex-beam control.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reduced model drops a residual carrier phase e^{-2iOmega_a z} in the u1^2 term of Eq. (15) without a validity condition, so the detuning stability maps in Figs. 3(b) and 5 may not correspond to the physical system (1)-(3).","rationale":"The reader's weakest_assumption correctly identifies the rotating-wave reduction as the soft spot, and I agree that dropping residual phase factors without a validity condition is dangerous. However, close re-derivation shows the reader's specific example of a phase factor in the u2*u3 drive of Eq. (14) is wrong: that term, like the u1*u3 and u1*u2 drives, is phase-free with the paper's gauge choice in Eq. (9). The only nontrivial residual phase is in the u1^2 drive of Eq. (15), e^{-2iOmega_a z}. This is sufficient to undermine the nonzero-detuning results because the u1^2 term is the primary SHG drive and its phase error affects the relation between propagation constants and detuning. The central claim of stable composite vortex solitons at exact phase matching (Omega_a=Omega_b=Omega_3=0) is not threatened by this concern, since all residual phases vanish there. The paper also has internal inconsistencies (existence vs stability statements and power ranges) and an overstatement about 1 m propagation in lossless LiNbO3, but these are secondary to the model-reduction issue for the detuning scans. The reader's CONDITIONAL verdict remains appropriate: the authors should either state a validity condition for the RWA (e.g., |Omega_a| << 1 over the formation length) or verify the nonzero-detuning stability predictions against the full physical model. My concrete test would settle the question directly.","tokens_in":1085,"tokens_out":2624,"duration_ms":285395,"concrete_test":"Directly simulate the physical system (1)-(3) with the full poling function d(Z) of Eq. (6), or with the residual phase term e^{-2iOmega_a z} added to the u1^2 term in Eq. (15), for parameters on the stability boundaries of Fig. 3(b), e.g., (P,D)=(50,4), Omega_a=0.5, Omega_b=0.5, Omega_3=0.5. Compare the soliton evolution and stability with the prediction of the reduced model (14)-(16). If the soliton becomes unstable, fails to form, or its propagation constant shifts as beta_2 = 2beta_1 - 2Omega_a, the reduced model is invalid at nonzero detuning and the detuning stability maps must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing concern is the validity of the rotating-wave reduction at nonzero detuning, which underlies the stability maps in Figs. 3(b) and 5. Re-deriving the reduction from the physical three-wave system (1)-(3) with the dual-period poling (6) shows that the gauge choices in Eq. (9) cancel the carrier phases in all nonlinear terms except the u1^2 drive in Eq. (15), which retains a factor e^{-2iOmega_a z} in scaled units. The paper's Eq. (15) omits this factor, justified only as 'neglect rapidly oscillating terms' with no validity condition. This is not rapidly oscillating when Omega_a is O(1) over the soliton-formation distance: the phase advances by ~2Omega_a per unit z, and the soliton forms over z ~ 6 (Figs. 6,7). Thus for the scanned ranges in Fig. 3(b), the dropped phase materially changes the nonlinear coupling. If included, the stationary condition for u2 would shift (beta_2 = 2beta_1 - 2Omega_a instead of beta_2 = 2beta_1), altering existence and stability boundaries. The exact-matching result (Omega_a=Omega_b=Omega_3=0) is unaffected, so the abstract's headline survives, but the detuning scans - important for experimental robustness - are not supported by the reduced model.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a quasi-phase-matched (QPM) photonic-crystal design, with a checkerboard transverse modulation and a dual-period modulation along the propagation direction, to support three-component (FF/SH/TH) composite vortex solitons in a purely quadratic medium. Starting from the physical three-wave system (1)-(3), the authors derive the reduced scaled model (14)-(16), solve it for stationary solutions, and report two families of four-peak solitons: rhombic and square-shaped. The FF component carries topological charge s=1, the SH component s=2, and the TH component s=-1, with the TH anti-vortex produced by a cascade through the SH quadrupole. Stability is checked by the Vakhitov-Kolokolov criterion and by direct perturbed propagation to z=1000-1300, and LiNbO3 parameters are used to estimate a physical propagation distance of about 1 m.","tokens_in":14200,"tokens_out":13320,"duration_ms":136940,"significance":"If the results hold, this is a substantial step beyond earlier QPM vortex-soliton work: it adds the third-harmonic component and shows that its anti-vortex structure arises naturally from the cascaded quadratic nonlinearity, while the checkerboard lattice can stabilize the composite object over many diffraction lengths. The paper gives direct numerical evidence for existence and stability at exact phase matching, and the experimental parameter estimates are useful and concrete. The main weakness is that the detuning scans, which are presented as control-parameter results, rely on a rotating-wave reduction whose validity at nonzero detuning is not established; this limits the support for the detuning-dependent stability maps but does not undermine the exact-phase-matching headline result.","major_comments":[{"comment":"The reduction from the physical three-wave system (1)-(3) to the scaled model (14)-(16) is load-bearing for all results at nonzero detuning, but it is not fully derived and the rotating-wave approximation is applied without a quantitative validity condition. With the gauge defined in Eq. (9), the nonlinear drives generated by the m=+1 and l=-1 Fourier harmonics of Eq. (7) do not all lose their carrier phases when the detunings are nonzero; in particular, the u1^2 drive in Eq. (15) retains a residual phase factor of the form e^{-2iΩ_a z} (up to sign, depending on convention). The statement in Section II that 'we neglect rapidly oscillating terms' is not sufficient: for Ω_a of order unity and soliton formation over z≈6, as seen in Figs. 6 and 7, this phase advances by ~2Ω_a per unit z and is not rapidly oscillating on the soliton-formation scale. Consequently, the stability diagrams in the (Ω_a, Ω_b) plane (Fig. 3(b)) and the H and β curves versus detuning (Fig. 5) may not correspond to the original physical system. The Ω_a=Ω_b=Ω_3=0 results are not affected, but the detuning branches of the paper need either an explicit smallness condition on Ω_a and Ω_b (showing that the dropped phases are negligible) or a re-simulation that retains the residual phases.","section":"Section II, Eqs. (9), (14)-(16)"},{"comment":"The detuning stability maps are the only support for the claim that the solitons can be controlled by the detuning parameters, and they inherit the uncertainty described above. As written, the manuscript asks the reader to accept the reduced model (14)-(16) as exact for all Ω_a, Ω_b in the scanned ranges, but this is not established. At minimum, the authors should display the full intermediate reduction, identify precisely which Fourier harmonics are kept, and state the parameter range in which the RWA is controlled. If the residual phases are genuinely negligible for the scanned parameters, a quantitative check (e.g., comparing a few full-model simulations of Eqs. (1)-(3) with the reduced model at representative Ω_a, Ω_b) would resolve the concern.","section":"Section III, Figs. 3(b) and 5"}],"minor_comments":[{"comment":"The transformation is written only as u_j = (ω_1/(n_1 I_0))^{-1/2} A_j e^{i(...)Z}, but the inverse transformation (A_j in terms of u_j) is not stated; this makes the sign convention ambiguous and complicates verification of the linear terms in Eqs. (14)-(16).","section":"Section II, Eq. (9)"},{"comment":"The perturbation protocol for the 'perturbed evolution' used to test stability is not specified; please state the perturbation type and amplitude, and the number of independent realizations, so that the stability region boundaries in Fig. 3 are reproducible.","section":"Section III, stability simulations"},{"comment":"The sentence 'Outside of the stability regions, unstable solutions have not been found' is potentially confusing in view of the later statement that unstable solitons exist outside the boundaries in Fig. 3(b); clarify that the former refers to the (P,D) plane only.","section":"Section III, text near Fig. 3(a)"},{"comment":"There are several typographical and grammatical errors, e.g., 'as a chains of rectangles' in the abstract, 'once again' in Section III, and 'shaoe'/'sahped' in the caption of Fig. 5.","section":"Abstract and Introduction"},{"comment":"The nonlinear coefficient d33 is cited to a Wikipedia article; please replace this with a primary or standard reference for lithium niobate's nonlinear tensor.","section":"Reference [66]"},{"comment":"The conversion from z=1300 to a physical distance of 1 m is not stated explicitly in Table I (which gives z=1000 as 76.5 cm); please add the z=1300 conversion to make the 'up to ~1 m' claim directly traceable.","section":"Section IV and Table I"}],"recommendation":"major_revision","confidential_remarks":"The exact-phase-matching results appear sound and are likely publishable, but the paper should not be accepted with the detuning scans in their current form unless the rotating-wave reduction is made explicit and its validity region is quantified. A representative comparison with the full model (1)-(3) at a few nonzero detunings would be a convincing fix; restricting the detuning claims to the range where the RWA is controlled would also be acceptable. The novelty relative to Ref. [33] is clear through the TH component and the anti-vortex cascade, so the manuscript is worth a revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a serious numerical extension of the group's PRL 2023 SHG vortex-soliton work to a full FF/SH/TH cascade. The genuinely new pieces are the dual-period poling that phase-matches both cascade steps and the topological cascade: FF vortex (s=1) to SH quadrupole (s=2) to TH anti-vortex (s=-1). At exact phase matching the reduced model (14)-(16) is standard, the VK criterion holds where the paper says, and direct evolution over z=1000 supports the stability claim. The LiNbO3 parameter estimates are plausible, though the '1 m propagation' statement ignores losses and photorefraction.\n\nThe soft spot is the nonzero-detuning analysis. The gauge in Eq. (9) does not cancel the residual phases in the nonlinear terms. I re-derived the u1^2 drive in Eq. (15) and find a residual e^{-3iΩ_a z} (from the m=+1 poling term), and the u1 u2 drive in Eq. (16) carries e^{-iΩ_3 z}. These are not rapidly oscillating when Ω is O(1). The paper's 'we neglect rapidly oscillating terms' comes with no validity condition, so the stability boundaries in Fig. 3(b) and Fig. 5 are not tied to the physical system (1)-(3). The exact-matching results are safe because those phases vanish. The stress-test's specific exponent may be off, but the concern is sound.\n\nThe text also contradicts itself: 'outside of the stability regions, unstable solutions have not been found' but Fig. 4 shows unstable rhombic branches; and Pmax≈350 in Fig. 3 vs P≤100 at D=4 in Fig. 4. The perturbation protocol for the stability runs is not specified. These are fixable but need attention.\n\nThis paper deserves a serious referee. The central existence claim at exact matching is defensible, and the TH anti-vortex is a nice addition. A referee should push for (i) either a validity condition for the RWA or a direct check of the detuning maps against the original equations, (ii) reconciliation of the existence/stability statements, and (iii) numerical details.","headline":"Solid exact-matching extension to THG, but the detuning stability maps rest on an unjustified rotating-wave reduction.","tokens_in":14995,"tokens_out":14160,"would_cite":true,"duration_ms":128109,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.65.Tg","42.65.Ky","42.70.Qs"],"model":"deepseek-v4-flash","headline":"The paper reports stable composite vortex solitons in a three-dimensional quasi-phase-matched photonic crystal, whose three color components carry topological charges $s=1$, $s=2$, and $s=-1$.","keywords":["vortex solitons","third harmonic generation","quasi-phase matching","photonic crystals","quadratic nonlinearity","cascaded nonlinear interactions","lithium niobate","anti-vortex"],"falsifier":"Integrate the original physical equations (1)-(3) with the complete poling function $d(Z)$ from Eq. (6) at a detuned point such as $(\\Omega_a,\\Omega_b)=(1,0)$ and compare the soliton profiles and stability boundaries with the reduced model (14)-(16); a visible mismatch would falsify the reduction and the claimed stability maps.","tokens_in":13674,"feed_emoji":"🌀","tokens_out":12760,"duration_ms":122376,"temperature":0.7,"pith_summary":"The paper argues that a purely quadratic ($\\chi^{(2)}$) medium can host stable composite vortex solitons once third-harmonic generation is included through a two-period quasi-phase-matched poling. The soliton is a three-color bound state: a fundamental-frequency vortex with topological charge (winding number) $s=1$, a second-harmonic quadrupole with $s=2$, and a third-harmonic anti-vortex with $s=-1$, created by cascaded quadratic interactions. In the reduced three-wave model, two families of four-peak solitons are found, rhombic and square, with the rhombic family stable over a broader range of total power and detuning. In lithium niobate parameters the same solitons are estimated to propagate stably for about one meter, at peak intensities below the threshold where cubic nonlinearity would matter. If the claim is right, it supplies a route to stable vortex solitons and higher-order optical states without invoking cubic nonlinearities.","feed_headline":"Three-color vortex solitons survive a full meter","feed_subtitle":"FF vortex, SH quadrupole, and TH anti-vortex form stable solitons in a checkerboard crystal.","key_machinery":"The load-bearing reduction is the near-resonance approximation that turns the physical three-wave system (1)-(3) into the scaled model (14)-(16): only the $m,l=\\pm1$ Fourier harmonics of the two-period poling are kept and the remaining oscillatory factors are dropped. In that model, the checkerboard sign $\\sigma(x,y)$ multiplies all nonlinear couplings, giving the transverse lattice structure that localizes the four peaks. The phase-matching identities (23)-(24) tie the component phases to the local sign of $\\sigma$, and the lattice's $D_4$ symmetry caps the phase circulation at charge magnitude one, which is what converts the cascade FF($s=1$) to SH($s=2$) to TH($s=-1$) into an anti-vortex rather than a higher-charge state. Stability is then asserted via the slope criterion $d\\beta/dP>0$ and by direct perturbed propagation over $z=1000$.","core_discovery":"On the paper's terms, the central discovery is that adding third-harmonic generation to a checkerboard quasi-phase-matched quadratic crystal produces a new class of stable composite vortex solitons with a nontrivial charge pattern. In the scaled system (14)-(16), stationary solutions $u_p=\\phi_p e^{i\\beta_p z}$ are built from four intensity peaks arranged as a rhombus or a square, selected by the relative phase $\\phi_d=0$ or $\\pi$ in the phase-matching conditions (23)-(24). The FF component is a vortex with $s=1$; the SH component has $s=2$, read as a quadrupole after phase folding; the TH component has $s=-1$, called an anti-vortex. Perturbed simulations keep both shapes stable up to $z=1000$, about 100 diffraction lengths, and to $z=1300$, about one meter in lithium niobate. The paper also shows that launching a vortex beam with unit winding number in the FF channel excites the solitons, with broader input beams favoring square shapes and narrower beams favoring rhombic shapes.","pith_inferences":["If the reduced model is accurate only near exact phase matching, then the stability maps computed at nonzero $\\Omega_a$ and $\\Omega_b$ may shift when the dropped oscillatory terms are restored; a direct comparison with the full three-wave equations would settle this.","The same two-period QPM cascade could be exported to higher harmonic orders: whenever the lattice symmetry caps vorticity, the charge sequence of doubling then subtracting one should generically produce anti-vortex components.","The demonstrated dependence of soliton shape on input beam width hints at an all-optical control channel, where the same crystal could switch between rhombic and square solitons by changing the launch profile."],"forward_implications":["Both rhombic and square composite vortex solitons remain stable in the model up to $z=1000$, about 76.5 cm in lithium niobate, and simulations still show stability at $z=1300$, about 1 m.","The rhombic family has a broader stability region, with total power stable up to $P\\approx 350$, versus $P\\approx 140$ for the square family, and lower Hamiltonian values for the same parameters.","A single vortex input with winding number 1 is enough to excite the three-color solitons; the output shape can be selected by the input beam width.","The third-harmonic anti-vortex is produced by the cascade FF($s=1$) $\\to$ SH($s=2$) $\\to$ TH($s=-1$), so no direct seeding of the TH component is required.","In lithium niobate, the peak intensities of all three components stay below about 1 GW/cm², so cubic nonlinearities can be neglected and the QPM structure is within current fabrication capabilities."],"supporting_citations":[{"why":"Gives the checkerboard QPM platform where stable FF/SH vortex solitons were previously predicted; this paper extends it to include THG.","marker":"[33]"},{"why":"Demonstrates quasi-phase-matched third-harmonic generation in a quasi-periodic superlattice, the mechanism this paper builds on for THG.","marker":"[34]"},{"why":"Shows a three-dimensional lithium niobate nonlinear photonic crystal can be fabricated, underpinning the experimental-feasibility claim.","marker":"[40]"},{"why":"Establishes the vorticity cutoff in nonlinear photonic crystals, used to explain why the TH component is an anti-vortex with charge magnitude one.","marker":"[62]"},{"why":"Supplies the slope criterion $d\\beta/dP>0$ that the paper applies to test stability of both soliton families.","marker":"[65]"},{"why":"Provides the $d_{33}=27$ pm/V value for lithium niobate used in the feasibility estimate and unit conversions.","marker":"[66]"},{"why":"Provides the Sellmeier refractive indices of lithium niobate used to convert scaled units to physical lengths and powers.","marker":"[67]"},{"why":"Provides the third-order susceptibility of lithium niobate, used to argue that cubic nonlinearity is negligible at the relevant intensities.","marker":"[68]"}],"fun_headline_variants":["Three-color vortex solitons stable for a full meter","Checkerboard crystal hosts stable composite vortices","Vortex solitons with charge trio survive 1 m","Stable three-wave vortex solitons in photonic crystals","Quasi-phase-matched crystal yields stable vortex solitons"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the reduced three-wave model, which drops rapidly oscillating terms, stays faithful to the physical crystal when the phase mismatches are not exactly zero; if those terms contribute, the predicted stability regions do not belong to the original model.","fun_headline_variants_meta":{"raw":{"variants":["Three-color vortex solitons stable for a full meter","Checkerboard crystal hosts stable composite vortices","Vortex solitons with charge trio survive 1 m","Stable three-wave vortex solitons in photonic crystals","Quasi-phase-matched crystal yields stable vortex solitons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000834,"raw_usage":{"total_tokens":3704,"prompt_tokens":1071,"completion_tokens":2633,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":687,"completion_tokens_details":{"reasoning_tokens":2553}},"tokens_in":687,"tokens_out":2633,"duration_ms":22503,"temperature":1.0,"reasoning_tokens":2553,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:11:38.909926+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Integrate the original physical equations (1)-(3) with the complete poling function $d(Z)$ from Eq. (6) at a detuned point such as $(\\Omega_a,\\Omega_b)=(1,0)$ and compare the soliton profiles and stability boundaries with the reduced model (14)-(16); a visible mismatch would falsify the reduction and the claimed stability maps.","supporting_citations":[{"cited_title":"Zhang, Z","cited_arxiv_id":null,"evidence_quote":"Gives the checkerboard QPM platform where stable FF/SH vortex solitons were previously predicted; this paper extends it to include THG."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows a three-dimensional lithium niobate nonlinear photonic crystal can be fabricated, underpinning the experimental-feasibility claim."},{"cited_title":"Ferrando, Vorticity Cutoff in Nonlinear Photonic Crystals, Phys","cited_arxiv_id":null,"evidence_quote":"Establishes the vorticity cutoff in nonlinear photonic crystals, used to explain why the TH component is an anti-vortex with charge magnitude one."},{"cited_title":"Generally, d0 = (2 /π )d33","cited_arxiv_id":null,"evidence_quote":"Provides the $d_{33}=27$ pm/V value for lithium niobate used in the feasibility estimate and unit conversions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Sellmeier refractive indices of lithium niobate used to convert scaled units to physical lengths and powers."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the third-order susceptibility of lithium niobate, used to argue that cubic nonlinearity is negligible at the relevant intensities."}],"review_version":1}