{"id":"5494a16a-4bad-4e0e-aecd-3de5a26b54e1","arxiv_id":"2507.09525","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"Tilted polarization stripes in a quadratic crystal create an effective gauge field, and numerical simulations report stable chiral solitons whose fundamental and second-harmonic parts carry opposite energy currents.","lead":"This paper proposes making an artificial gauge field in a nonlinear crystal by tilting the polarization stripes used for quasi-phase matching, and reports numerical chiral solitons that carry opposite energy flows in the fundamental and second-harmonic waves. A generalist might read it because it suggests a new, all-optical route to chiral transport without linear coupling.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The stationary phase-gradient ansatz (20)-(22) is inconsistent: for real envelopes it forces ηu=ηv=0 and hence α=0, so Eq. (36) for chiral flow is not derived for α≠0.","rationale":"The reader's REJECT verdict is based on exactly this inconsistency, and I agree. The ansatz is not merely an unproven approximation; it is internally contradictory for localized functions. Since Eq. (36) is built directly from ηu and ηv obtained under that ansatz, the central semi-analytical result is unsupported. The paper's numerical simulations may still indicate stable localized structures, and those could be real, but the text does not provide code or data, and the analytic explanation of the chiral current and optimal inclination is what the abstract and conclusion advertise. The additional factor-of-2 error in the current algebra reinforces the rejection but is secondary; even a corrected factor would not repair the imaginary-part inconsistency. If the authors could supply an explicit stationary solution with nonlinearly varying phase and re-derive Sc, the verdict could change, but as written REJECT is appropriate.","tokens_in":18797,"tokens_out":11671,"duration_ms":133762,"concrete_test":"Symbolically substitute ψ_u=φu e^{iηu x}, ψ_v=φv e^{iηv x} with real φ into Eqs. (14)-(15) and separate real and imaginary parts before using Eq. (22). If the imaginary equations reduce to ηu φu'=0 and ηv φv'=0, then any localized solution of this ansatz has ηu=ηv=0 and α=0, so the ansatz cannot describe the reported α≠0 solitons; this directly invalidates Eqs. (30)-(36). A complementary numerical check is to evaluate the L2 norm of Im[β1ψu + (1/2)ψu'' + Ωψu + 2e^{iαx}ψu*ψv] for the ITP-converged profile; for a true stationary solution of the assumed form it must vanish, and for the profiles shown it will not.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Substituting ψ_u=φu e^{iηu x} and ψ_v=φv e^{iηv x} with real φu,φv into the stationary equations (14)-(15), and using 2ηu-ηv=α, makes the nonlinear phases match, but the kinetic terms generate imaginary parts -ηu φu' and -(ηv/2)φv' with no counterpart on the left-hand side. For localized profiles φu' and φv' do not vanish identically, so ηu=ηv=0; Eq. (22) then implies α=0. Thus for tilted stripes (α≠0) no localized stationary soliton of the assumed form exists, and the derivation of ηu,ηv in Eqs. (30)-(31) and of Sc in Eq. (36) is not supported. There is also an independent algebraic inconsistency: from P=Pu+2Pv, current conservation ηu Pu+ηv Pv=0 and Eq. (22) give ηu=αPv/P, not αPv/(2P); recomputing Su and Sv changes Sc by a factor of 2. The numerical existence or stability of some chiral soliton is not disproved by this, but the claimed semi-analytic chiral-flow formula and the associated optimal-inclination argument rest on this invalid ansatz.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a quasi-phase-matched (QPM) quadratic nonlinear photonic crystal in which tilted ferroelectric stripes encode an x-dependent phase that acts as an artificial gauge field in the nonlinear coupling. The authors report stable chiral solitons whose fundamental-frequency (FF) and second-harmonic (SH) components carry opposite energy currents, and they derive a semi-analytical formula for the chiral energy flow, Sc = -αβ1PuPv/P (Eq. 36). They further study the dependence of Sc on power, detuning, and stripe angle, identify an optimal inclination that maximizes Sc, show that kicked solitons move and collide nearly elastically, and provide an experimental parameter estimate. The numerical exploration is extensive, but the analytic derivation underlying Eq. (36) contains several specific technical errors that are load-bearing for the paper's central claims.","tokens_in":19097,"tokens_out":12967,"duration_ms":133310,"significance":"If correct, the proposal would be an interesting way to engineer synthetic gauge fields through the QPM geometry rather than through linear couplings, and the reported stable chiral solitons with coupled FF-SH currents would be a useful addition to nonlinear optics. The numerical stability maps (Fig. 4), the mobility and collision simulations (Figs. 9-10), and the experimental parameter estimates (Sec. V) are potentially valuable. However, the paper's headline analytic result is not sound: the stationary phase-ansatz used to derive Eq. (36) is inconsistent for nonzero α, and the derivation of the current formula contains algebraic and conceptual errors. Since the optimal-inclination claim is presented as following from this formula, the main conceptual conclusions are not established by the present analysis. The numerical phenomenology might still be publishable, but not with the current analytic scaffolding.","major_comments":[{"comment":"The stationary ansatz with real envelopes and constant phase slopes is inconsistent for α≠0. Substituting ψu=φu e^{iηu x} and ψv=φv e^{iηv x} with real φu,φv into Eq. (14) produces an imaginary term -ηu φu′ (and similarly -ηv φv′/2 in Eq. (15)) with no balancing imaginary source once Eq. (22) enforces the nonlinear phase. For localized non-constant envelopes this forces ηu=ηv=0 and hence α=0. Thus the assumed form cannot describe the tilted-stripe solitons, and the derivation of Eqs. (30)-(31) and Eq. (36) is unsupported.","section":"§II, Eqs. (20)-(22) and (14)-(15)"},{"comment":"Eq. (28) is not a valid identity. Integrating Eq. (27) over x gives a boundary value of ju+2jv, not the local current density; the vanishing of the integral does not imply the local current vanishes. Consequently Eq. (29) does not follow as a local statement. In addition, direct calculation from Eqs. (7)-(8) shows that the nonlinear terms cancel in ∂z(|u|^2+2|v|^2), so the right-hand side of Eq. (27) should be zero, not -8Im(e^{iαx}ψu*^2ψv).","section":"§II, Eqs. (27)-(29)"},{"comment":"The quoted values ηu=αPv/(2P) and ηv=-αPu/P are algebraically inconsistent with Eq. (22): substituting them gives 2ηu-ηv=α(Pu+Pv)/P, not α. Solving Eq. (22) together with the integrated current conservation ηuPu+ηvPv=0 gives ηu=αPv/P and ηv=-αPu/P. With these corrected values, Su and Sv are still equal and opposite, but Sc=Su-Sv becomes -2β1αPuPv/P, a factor of 2 different from Eq. (36). Thus the chiral-flow formula is not derived as stated.","section":"§II, Eqs. (30)-(31) and (36)"},{"comment":"Even if Eq. (36) were algebraically correct, the claimed semi-analytical prediction is not independent: β1, Pu, and Pv are not input parameters but are quantities read off the numerically obtained solitons. The optimal inclination α(OI) is therefore a property of the numerically determined product αβ1PuPv/P, not a parameter-free prediction. The text's statement that Sc is \"a quantity determined by the input parameters\" is misleading in this context.","section":"§III, Eq. (36) and Fig. 7"}],"minor_comments":[{"comment":"The manuscript contains multiple encoding artifacts, such as \"dieresis.ts\" in the references and \"£º\" after Eq. (27); these should be cleaned before any resubmission.","section":"Throughout"},{"comment":"The text refers to panels (c1,c2) in Figs. 2 and 6, but the corresponding figure captions and panel layouts use (a3,a4,b3,b4) and (a,b); the panel labels and text references should be reconciled.","section":"§III, Fig. 2 and Fig. 6 captions"},{"comment":"The conclusion mentions verification with the \"RTP method,\" whereas the introduction and Sec. III describe the ITP method and direct simulations; the acronyms and methods should be used consistently.","section":"§VI, Conclusion"}],"recommendation":"reject","confidential_remarks":"The numerical study may contain useful phenomenology, but the analytic framework that gives the paper its central claim is invalid at several independent points: the stationary ansatz does not admit nonzero α, the current derivation conflates local and integrated quantities, and the algebra leading to Eq. (36) is inconsistent. These are not local presentation issues; they affect the main result and the optimal-inclination interpretation. I would not invite a standard major revision because repairing the derivation would require replacing the core ansatz and reworking the paper's central argument. The numerical results, if separated from the invalid analytic claims, might be suitable for a different presentation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read. The setup is genuinely neat: tilting the QPM stripes creates an x-dependent phase in the quadratic coupling, which you can view as an artificial gauge field acting on the nonlinear interaction. That specific mechanism I haven't seen before. And the numerics look like real work: ITP convergence, BdG stability, direct simulation, stability maps in the (Ω, α) plane, collisions, and a plausible experimental parameters section. If the existence of stable chiral solitons with counter-propagating FF/SH currents is correct, that is a valid subfield result.\n\nThe soft spots are, unfortunately, at the center. The semi-analytic derivation of the chiral flow doesn't hold. Substitute the ansatz (20)-(22) into the stationary equations (14)-(15). The kinetic terms generate imaginary parts -iηu φu' and -(i/2)ηv φv' that nothing balances; for localized real envelopes this forces ηu=ηv=0 and, from (22), α=0. So the derivation of ηu,ηv and formula (36) is not valid for tilted stripes. There is also a standalone algebraic error: integrating the current-conservation expression (29) gives ηu Pu + ηv Pv = 0, which with (22) yields ηu = α Pv/P, not α Pv/(2P). That changes Su, Sv and Sc by a factor of two. And the 'semi-analytical solution' is not a prediction—β1, Pu, Pv are read off the same numerically found solitons, so the optimal-inclination curve is a property of that numerical data, not an independent derivation.\n\nI want to be clear: none of this disproves the numerical chiral solitons themselves. The stability maps and the observed counter-propagating currents may well be correct. But the paper's central analytic claim—the chiral flow formula and the optimal inclination argument—rests on the invalid ansatz. As written, I would not accept it. The authors could fix this by deriving the phase-gradient properly (allowing x-dependent η) or by dropping the semi-analytic claim and reporting the numerics alone. The topic is good enough that I'd send it to a knowledgeable referee rather than desk reject; hopefully they can untangle whether the numerical states survive scrutiny once the derivation is corrected.","headline":"Neat idea and plausible numerics, but the semi-analytic derivation of the chiral flow is wrong as written; the optimal-inclination argument rests on it.","tokens_in":19603,"tokens_out":6493,"would_cite":false,"duration_ms":65609,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.65.-k","42.65.Ky"],"model":"deepseek-v4-flash","headline":"Tilted crystal stripes create chiral light solitons.","keywords":["chiral solitons","quasi-phase-matching","artificial gauge field","quadratic nonlinear optics","second-harmonic generation","energy flow","photonic crystals"],"falsifier":"Take the numerically converged chiral soliton from the paper's imaginary-time method, substitute it into the stationary equations (14)-(15), and compute the imaginary part of the FF equation; if the imaginary residue is not zero, or if the integrated energy current of the exact solution differs from Eq. (36), the central claim is contradicted. Alternatively, solve for stationary solutions without the constant-slope ansatz and check whether any localized solution with $\\eta_u \\neq 0$ exists.","tokens_in":18609,"feed_emoji":"🌀","tokens_out":5971,"duration_ms":67575,"temperature":0.7,"pith_summary":"The paper proposes that tilting the ferroelectric polarization stripes in a quasi-phase-matched quadratic crystal acts as an artificial gauge field, but one carried by the nonlinear coupling rather than by linear tunnelling or coupling. In this setting, the fundamental-frequency (FF) and second-harmonic (SH) components of a soliton lock into a chiral state: their energy currents are equal in magnitude and opposite in direction, so the two waves circulate around each other. The paper derives a semi-analytic expression for the chiral energy flow, $S_c = -\\alpha\\beta_1 P_u P_v/P$, and shows numerically that stable solitons exist, that an optimal stripe inclination maximizes the current, and that the solitons survive kicks and collisions. If correct, this gives a way to generate chiral optical currents using only standard poled crystals and a single pump beam.","feed_headline":"Tilted crystal stripes create chiral light solitons","feed_subtitle":"Fundamental and second-harmonic waves carry equal, opposite currents, offering chip-scale chiral transport.","key_machinery":"The load-bearing mechanism is the tilted-stripe quasi-phase-matching phase $\\phi_d = (2\\pi/\\Lambda)\\tan\\theta\\,X$, which converts the poling pattern into the phase factor $e^{i\\alpha x}$ in the dimensionless coupled equations. Applying the rotating-wave-like reduction to the $m=\\pm1$ Fourier components, this phase becomes a constant artificial vector potential $\\alpha$ acting through the nonlinear terms. The current formula then follows from the stationary ansatz $\\psi_u = \\varphi_u e^{i\\eta_u x}$, $\\psi_v = \\varphi_v e^{i\\eta_v x}$ with the resonance condition $2\\eta_u - \\eta_v = \\alpha$, which forces the two currents to balance.","core_discovery":"On the paper's own terms, the central discovery is that a spatially linear phase ramp imposed on the quasi-phase-matching modulation, achieved by tilting the domain stripes by an angle $\\theta$, behaves as an equivalent magnetic field in the coupled FF-SH equations. Under this field the stationary soliton develops opposite linear phase slopes for the two components, $\\eta_u = \\alpha P_v/(2P)$ and $\\eta_v = -\\alpha P_u/P$, so the FF and SH currents $S_u$ and $S_v$ are equal and opposite and the total chiral flow is $S_c = -\\alpha\\beta_1 P_u P_v/P$ (Eq. 36). The paper verifies by direct simulation and Bogoliubov-de Gennes analysis that these chiral solitons are stable over wide parameter ranges, finds an optimal inclination angle that maximizes $|S_c|$ for fixed power and detuning, and demonstrates kicked motion and near-elastic collisions.","pith_inferences":["A natural extension is to repeat the tilted-stripe construction in two-dimensional poling patterns, where the synthetic field could support chiral edge transport or vortex solitons rather than single-peaked currents.","Measuring $S_c$ versus tilt angle in a short periodically poled sample would test the optimal-inclination prediction directly; the peak location should track the paper's $\\alpha^{(\\mathrm{OI})}(P,\\Omega)$ curves.","An exact stationary-solution check, solving Eqs. (14)-(15) without the constant-slope ansatz, would show whether the analytic formula survives beyond the assumed phase-ramp form.","The same tilted-stripe gauge field may also couple to other nonlinear excitations in the crystal, such as vortices or dissipative solitons, giving a general tool for chiral nonlinear optics."],"forward_implications":["Stable chiral energy transport can be engineered into a bulk quadratic crystal by geometry alone, using only the poling tilt angle and a single pump beam.","The closed-form current $S_c = -\\alpha\\beta_1 P_uP_v/P$ provides a design rule for maximizing chiral flow by choosing the optimal inclination for given power and detuning.","Because the solitons respond to kicks and collide nearly elastically, they could serve as movable, reconfigurable carriers in nonlinear-photonic routing.","The stability regions in the $(\\Omega,\\alpha)$ plane map out where such chiral solitons can be observed in practice, guiding future experiments in lithium-niobate-style crystals."],"supporting_citations":[{"why":"Supplies the geometric-phase formalism that identifies the QPM phase as an equivalent vector potential.","marker":"[34]"},{"why":"Provides the rescaling and phase-matching reduction used to reach the two dimensionless coupled equations.","marker":"[37]"},{"why":"Establishes vortex solitons in quasi-phase-matched photonic crystals, the prior demonstration of soliton solutions and stability methods in this system.","marker":"[49]"},{"why":"Gives the earlier chiral-soliton-in-gauge-field result that this paper contrasts with its nonlinear-coupling gauge field.","marker":"[69]"},{"why":"Foundational quasi-phase-matching theory for second-harmonic generation that the tilt-phase mechanism extends.","marker":"[47]"},{"why":"Documents that the poling modulation leaves the linear refractive index unchanged, so the gauge field acts only through the nonlinear coupling.","marker":"[60]"},{"why":"Experimental observation of an all-optical Stern-Gerlach effect in QPM crystals, the closest prior demonstration of QPM-induced equivalent magnetic fields.","marker":"[30]"}],"fun_headline_variants":["Tilted stripes create chiral solitons with opposing currents","Nonlinear gauge field yields chiral solitons in photonic crystals","Optimal tilt boosts chiral flow in quadratic crystals","Chiral solitons move and collide elastically in tilted crystals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole analytic result depends on the assumption that a localized soliton can be written as a real profile times a plane-wave phase ramp with constant slope for each component; if the stationary equations forbid such a form, the current formula does not follow.","fun_headline_variants_meta":{"raw":{"variants":["Tilted stripes create chiral solitons with opposing currents","Nonlinear gauge field yields chiral solitons in photonic crystals","Optimal tilt boosts chiral flow in quadratic crystals","Chiral solitons move and collide elastically in tilted crystals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000317,"raw_usage":{"total_tokens":1788,"prompt_tokens":936,"completion_tokens":852,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":552,"completion_tokens_details":{"reasoning_tokens":782}},"tokens_in":552,"tokens_out":852,"duration_ms":9864,"temperature":1.0,"reasoning_tokens":782,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:55:10.478606+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the numerically converged chiral soliton from the paper's imaginary-time method, substitute it into the stationary equations (14)-(15), and compute the imaginary part of the FF equation; if the imaginary residue is not zero, or if the integrated energy current of the exact solution differs from Eq. (36), the central claim is contradicted. Alternatively, solve for stationary solutions without the constant-slope ansatz and check whether any localized solution with $\\eta_u \\neq 0$ exists.","supporting_citations":[{"cited_title":"Karnieli and A","cited_arxiv_id":null,"evidence_quote":"Supplies the geometric-phase formalism that identifies the QPM phase as an equivalent vector potential."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the rescaling and phase-matching reduction used to reach the two dimensionless coupled equations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes vortex solitons in quasi-phase-matched photonic crystals, the prior demonstration of soliton solutions and stability methods in this system."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the earlier chiral-soliton-in-gauge-field result that this paper contrasts with its nonlinear-coupling gauge field."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Foundational quasi-phase-matching theory for second-harmonic generation that the tilt-phase mechanism extends."},{"cited_title":"Zhang, Y","cited_arxiv_id":null,"evidence_quote":"Documents that the poling modulation leaves the linear refractive index unchanged, so the gauge field acts only through the nonlinear coupling."},{"cited_title":"Yesharim, A","cited_arxiv_id":null,"evidence_quote":"Experimental observation of an all-optical Stern-Gerlach effect in QPM crystals, the closest prior demonstration of QPM-induced equivalent magnetic fields."}],"review_version":1}