{"id":"42f13e1e-cac2-4aef-bfe9-99a7ef3f958c","arxiv_id":"2411.17954","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A symmetry-reduced eigenfunction matching method yields frequency- and time-domain scattering solutions for a rectangular junction of four finite-width waveguides.","lead":"Scientists solve wave scattering at a four-way junction of rectangular waveguides by splitting the problem into four symmetric pieces and matching eigenfunctions at the seams. The result is a scattering matrix and time-domain movies that could help design acoustic or optical waveguide networks.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Appendix C's right-arm reconstruction for the symmetric case uses odd sine modes, contradicting Eq. (11) and the even symmetry of the recovered field.","rationale":"The reader's weakest assumption concerns the unstated zero-extension in the derivative matching at (14a), integrated via H_mn in (42). That is a legitimate presentation gap, but the final matrix system in §3 is consistent with the zero-extension if one reads (14a) together with (41)–(42), so an expert implementer can recover the intended method. The more decisive problem is in Appendix C, which is the paper's explicit construction of the recovered full solution. For a symmetric incident wave, the right-arm field must be even in y. Eq. (11) gives a cosine combination; Appendix C writes sine modes. This is not a matter of interpretation: the sine formula violates the symmetry imposed by the incident wave and gives ∂_yφ≠0 on the interior line y=0, which is impossible for the smooth Helmholtz solution. Because the abstract and Sections 3/4 rely on the full solution being recovered correctly, the manuscript as written contains a concrete mathematical inconsistency in its central claim. The fix is straightforward (replace sin γ_n and γ̄_n by cos β_n and β̄_n, with the same coefficient combination), and the numerical figures may well have been produced with the correct formula; nevertheless, the stated reconstruction must be corrected and verified before the central claim can be accepted. This does not change the overall CONDITIONAL verdict, but it sharpens the reason for it.","tokens_in":20988,"tokens_out":15003,"duration_ms":120896,"concrete_test":"Run the following check for a symmetric incidence case (e.g., Fig. 4a: a1=a2=3, b1=b2=5, k=5, p=2). Assemble the right-arm field φ_N(x,y) for x>b2 in two independent ways: (i) using the Appendix C formula with sin γ_n(y); (ii) using Eq. (11) directly, evaluating φ_NN(−x,y) and φ_DN(−x,y) from the computed quadrant expansions. Compare the two fields pointwise and test the symmetry condition φ_N(x,y)=φ_N(x,−y) and the condition ∂_y φ_N(x,0)=0. The Appendix C formula will fail these tests (its ∂_y at y=0 is generically nonzero), whereas the Eq. (11) field passes. This settles whether the stated reconstruction is internally consistent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Appendix C (second displayed equation, last line) gives the full field φ_N(x,y) for a symmetric incident wave in the right arm x≥b2 as ∑_n ((A_n^NN − A_n^DN)/2) e^{−iγ̄_n(−x+b2)} sin(γ_n y), with γ_n=(2n+1)π/(2a1). This is not the function defined by the paper's own reconstruction (11). For a wave incident from x=−∞ that is symmetric about y=0, (11) requires φ(x,y)=φ(x,−y) everywhere. In the right arm the domain is y∈(−a1,a1), so the field must be even in y. Directly from (11), φ(x,y)=½[φ_NN(−x,|y|) − φ_DN(−x,|y|)], and both quadrant solutions have y-dependence cos(nπy/a1) on Ω1. The Appendix C formula instead uses sin((2n+1)πy/(2a1)), which is odd, and evaluates ∂_y φ(x,0) to the generically nonzero value ∑ ((A_n^NN−A_n^DN)/2) γ_n e^{iγ_n(x−b2)}. Hence the stated reconstructed solution violates the required symmetry and cannot satisfy the Helmholtz problem. Since the abstract's claim 'the solution to the problem on the full region is recovered' rests directly on these reconstruction formulas, this is a load-bearing inconsistency in the central claim as written.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper considers two-dimensional linear wave scattering at a junction of four rectangular waveguides with Neumann walls. The authors split the problem into four quadrant subproblems using reflection symmetry, solve each subproblem by eigenfunction matching, and reconstruct the full field from the quadrant solutions. They also construct an 8×8 block scattering matrix for the junction and compute time-domain pulses as superpositions of frequency-domain solutions. The central claims are that the full-region solution is recovered from the subproblems and that a scattering matrix for the junction is presented.","tokens_in":51,"tokens_out":23510,"duration_ms":213477,"significance":"If correct, the paper would supply a useful semi-analytic benchmark for multi-mode scattering at a four-waveguide junction and a building block for lattice-type multiple-scattering studies. The strengths are that no physical parameters are fitted, the quadrant matching framework is standard, an energy-conservation identity is checked for the Neumann-Neumann case, and the time-domain evaluation is reduced to matrix multiplication. However, the current manuscript contains load-bearing internal inconsistencies: the reconstructed right-arm field in Appendix C does not follow from the paper's own reconstruction formula (11), and the derivative matching derivation in Section 3 hides an essential zero-extension assumption. These issues must be corrected before the central claim can be accepted.","major_comments":[{"comment":"The right-arm formula for the symmetric case is internally inconsistent. For x≥b2, Eq. (11) reconstructs the field as (φNN(−x,y)−φDN(−x,y))/2 for y>0 and (φNN(−x,−y)−φDN(−x,−y))/2 for y<0. Both φNN and φDN in Ω1 have y-dependence cos(β_n y) from Eq. (18), and mirroring gives e^{−iβ_n(−x+b2)} = e^{iβ_n(x−b2)}. Appendix C instead gives ∑ ((A_n^NN−A_n^DN)/2) e^{−iγ̄_n(−x+b2)} sin(γ_n y), whose transverse modes are odd in y and cannot represent an even field. In particular, ∂_yφ_N(x,0) is generically nonzero, violating the symmetry imposed by (11). This directly undermines the abstract's claim that the full-region solution is recovered, and the formula must be corrected to use the proper even modes and wavenumbers.","section":"Appendix C"},{"comment":"The derivative matching step is not stated correctly. Condition (14a) is piecewise: on (0,a1) the Ω1 derivative equals the Ω2 derivative, while on (a1,b1) the derivative is zero on the rigid wall. Equation (41) instead writes the Ω1 derivative as equal to the Ω2 derivative for all y∈(0,b1), which is not the physical condition. Equation (42) then integrates the left side only over (0,a1), through H_mn, while the right side is integrated over (0,b1); this is a valid weighted-residual formulation only if the Ω1 derivative is implicitly zero-extended over (a1,b1). That zero-extension is never stated. Since all computed coefficients come from this and the analogous systems in Appendix A, the derivation must state the extension explicitly or present the weak form directly.","section":"§3, Eqs. (41)–(42)"},{"comment":"The odd-mode expansions are inconsistently indexed. Eq. (26) writes the reflected field as a sum from n=1, while γ_n is defined for n=0,1,2,... and the matching equation (49) sums from n=0. These cannot both be correct. If the n=0 term is truly omitted, the fundamental odd mode sin(πy/(2a1)) is missing from the expansion, which would violate completeness and affect all antisymmetric results. The same inconsistency appears in the scattering-matrix count q̃ in §4, which appears to exclude the lowest odd mode. The indexing must be corrected and made uniform across the definitions, the matching equations, and the scattering matrix.","section":"§2.3 and Appendix A"},{"comment":"The naming of the four subproblems is inconsistent. The enumeration in §2 defines Neumann-Dirichlet as Dirichlet at x=0 and Neumann at y=0, and Dirichlet-Neumann as Neumann at x=0 and Dirichlet at y=0. In contrast, §2.4 and §2.5 define ND with Neumann at x=0 and Dirichlet at y=0, and DN with Dirichlet at x=0 and Neumann at y=0. Since Eqs. (11)–(12) and Appendix C rely on which condition holds on y=0, the labels must be harmonized before the reconstruction and the scattering-matrix notation can be unambiguously interpreted.","section":"§2 and §2.4/§2.5"}],"minor_comments":[{"comment":"The sentence introducing Eq. (12) says the incident wave is antisymmetric about the line x=0, but both Eq. (12) and Appendix C use antisymmetry about y=0; please correct this.","section":"§2, before Eq. (12)"},{"comment":"The quadrature parameters Nk and Δk used in the time-domain computations for Figures 5–8 are not reported, so the time-domain visualisations are not reproducible as written.","section":"§5 and figure captions"},{"comment":"The caption of Figure 8 is garbled: it reads \"b1 = 4, b2 = a2 = 3, a2\" with a duplicated a2; please correct the parameter list.","section":"Figure 8 caption"},{"comment":"There are numerous typographical errors, including \"donated\" for \"denoted\", \"team\" for \"term\", \"donate\" for \"denote\", and \"probl em\" in the abstract; a careful copyedit is needed.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The Appendix C inconsistency is severe but local: the quadrant eigenfunction-matching framework is standard, and the errors appear correctable by reworking the reconstruction formulas and the matching derivation. I therefore do not recommend rejection, but the authors must verify every reconstruction formula against the original PDE and boundary conditions before resubmission."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"New here is the geometry: scattering by a rectangular junction of four finite-width waveguides, solved by symmetry reduction to four quadrant problems and eigenfunction matching. That specific configuration is not in the cited literature, and the assembled method is standard but competently laid out. The energy identity (45) and the matching plots are reasonable internal consistency checks, and the time-domain visualizations are a useful extra.\n\nThe main problem is Appendix C. For the symmetric incident wave, the reconstructed field in the right arm is written as a series in sin(γ_n y) with γ_n=(2n+1)π/(2a1). The reconstruction formula (11) forces the full field to be even in y in the right arm, and both quadrant solutions (NN and DN) have cos(β_n y) y-dependence there. So the appendix formula contradicts the paper's own symmetry reconstruction and cannot solve the Helmholtz problem. This looks like a copy-paste from the antisymmetric case, but as written it invalidates the abstract's claim that the full-region solution is recovered.\n\nSecond, Eq (41) states the derivative matching as an equality on the whole interval (0,b1), but condition (14a) requires the left-hand side to be zero on the wall (a1,b1). The subsequent inner products integrate only over (0,a1), so an implicit zero-extension is doing work that is never stated. A reader implementing (41) literally would get the wrong system. This is fixable by rewriting the matching equation with the zero-extension made explicit.\n\nThird, there is no convergence study, no independent benchmark, and no code. The truncation N=100 is justified by the energy identity and linewidth matching, which are self-consistency checks, not validation. The scattering-matrix section is a sketch: only two sub-matrix columns are given, and the general entries are described rather than listed.\n\nNone of this is fatal to the method. The symmetry decomposition and matching are sound, and the Appendix C error is very likely a fixable slip. But the paper as written has a load-bearing inconsistency in its central reconstruction formula, and the lack of verification means the numerics are only as credible as the authors' self-checks.\n\nThe paper is for people working on mode-matching for multi-port waveguides or on lattice models with finite-width arms. I would send it to a serious referee, and expect the report to demand a corrected reconstruction formula, an explicit statement of the zero-extension, and at least a convergence check or independent comparison. Until then, I wouldn't cite the full-solution formula as it stands.","headline":"Useful new geometry for mode-matching, but the Appendix C reconstruction contradicts the paper's own symmetry and needs correction before the results can be trusted.","tokens_in":21788,"tokens_out":9392,"would_cite":false,"duration_ms":75794,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35J05","35P25","78A50"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper establishes that wave scattering at a rectangular four-waveguide junction is recovered exactly from four quadrant problems, produced by splitting the domain with Neumann and Dirichlet mirror conditions, each solved by…","keywords":["four-waveguide junction","eigenfunction matching","symmetry decomposition","scattering matrix","time-domain scattering","Helmholtz equation","Neumann-Dirichlet quadrant problems","wave propagation"],"falsifier":"Compute the normal derivative of the reconstructed Neumann-Neumann solution on the rigid-wall segment y in (a1,b1) at x=-b2 from the coefficients obtained by (39)-(44); if that derivative is not zero to within the truncation error, the implicit zero-extension of the Omega_1 field is not supplying the boundary condition stated in (14a).","tokens_in":20766,"feed_emoji":"🌊","tokens_out":5653,"duration_ms":50571,"temperature":0.7,"pith_summary":"The paper solves a two-dimensional linear-wave scattering problem at a rectangular junction where four waveguides meet, with sound-hard (zero normal derivative) boundaries and an incident plane wave from one arm. It shows that the full solution is recovered from four quadrant problems, each carrying a different combination of Neumann and Dirichlet artificial conditions along the two mirror lines. Each quadrant is solved by eigenfunction matching, expanding the field in channel eigenfunctions and matching value and normal derivative across the interfaces. From these quadrant solutions the authors assemble a scattering matrix for all four ports and construct time-domain pulses as quadrature superpositions of frequency-domain fields. The central claim, read sympathetically, is that this symmetry-plus-matching construction yields converged coefficients that satisfy energy conservation.","feed_headline":"Four-waveguide junction solved from one quadrant","feed_subtitle":"Neumann/Dirichlet mirror symmetry plus eigenfunction matching yields a full scattering matrix and time-domain pulses.","key_machinery":"The load-bearing mechanism is the reflection-symmetry decomposition of the rectangular junction into four quadrants, combined with the eigenfunction matching method. Each quadrant problem is solved on three rectangles by separation of variables, and matching the expansions across the two interior interfaces produces a linear system of algebraic equations for the mode amplitudes. A named piece is the scattering matrix S (equation 46) with its even and odd blocks; the reconstruction identities (11)-(12) are what convert quadrant solutions back to the physical field.","core_discovery":"The central claim is that the Helmholtz solution on the full four-waveguide domain need not be computed directly. The two mirror symmetries of the junction split the problem into four quadrant problems: Neumann-Neumann, Dirichlet-Dirichlet, Neumann-Dirichlet, and Dirichlet-Neumann, named for the artificial boundary conditions imposed on the coordinate axes. Each quadrant field is written piecewise in three rectangular subregions, expanded in cosine or sine eigenfunctions with complex propagation constants, and the unknown mode amplitudes are fixed by matching the field and its normal derivative at the two interior interfaces x=-b2 and y=b1. Formulas (11) and (12) recombine the quadrant solutions into the full-domain field for even and odd incident waves. The paper packages the outcome as an 8-by-8 block scattering matrix with even and odd modes separated and symmetry-enforced zero blocks, and illustrates time-domain pulses obtained by summing frequency-domain solutions.","pith_inferences":["For equal waveguide widths (a1=a2), the even and odd block sizes of S differ by at most one propagating mode, which may produce near-zero reflection zeros at specific frequencies; plotting reflection coefficients against k would test this directly.","The zero-extension subtlety in the derivative matching suggests that a direct finite-element comparison at one or two parameter sets would provide an independent numerical check beyond the paper's internal matching plots.","The wide-spacing approximation invoked for the scattering matrix ignores evanescent coupling between nearby junctions; at lattice spacings comparable to the junction size, a transfer-matrix or full-domain calculation would be needed instead of S alone.","The same quadrant decomposition could be applied to junctions with higher rotational symmetry, but the artificial boundaries would become wedges whose corner singularities require the matching basis to be enriched beyond simple cosine and sine modes."],"forward_implications":["The 8-by-8 scattering matrix, with even and odd modes in separate blocks, gives outgoing amplitudes for any incident mode combination and can drive multiple-scattering simulations of a lattice of junctions.","Once the frequency-domain field matrix is computed for a given geometry, time-domain responses for many different incident spectra cost only one matrix multiplication per time step.","Because the reconstruction formulas and the matching systems are written for general a1, a2, b1, b2, the method covers junctions with unequal channel and rectangle widths, not only the symmetric case used for the scattering matrix.","The numerical results at N=100 satisfy the conservation-of-energy identities (45), which is the paper's internal check that truncation has converged.","A junction of this type can be embedded in a square lattice whose band structure could be studied through the junction scattering matrix, connecting to quantum-graph-inspired metamaterial design."],"supporting_citations":[{"why":"Supplies the symmetry-and-scattering-matrix approach used here for time-domain waveguide scattering.","marker":"[5]"},{"why":"Supplies the mode-matching without root-finding technique that the quadrant solution relies on.","marker":"[8]"},{"why":"The multifurcated-duct mode-matching analysis that this paper extends to four waveguides.","marker":"[17]"},{"why":"Justifies the convergence of the discrete mode-matching truncation used in the numerical systems.","marker":"[23]"},{"why":"Precedent for using symmetry to solve wave scattering in a waveguide-resonator geometry.","marker":"[34]"},{"why":"Provides the time-domain construction as a superposition of frequency-domain solutions via quadrature.","marker":"[40]"},{"why":"Supplies the generalized eigenfunction expansion method for two-dimensional acoustic time-domain scattering used here.","marker":"[41]"}],"fun_headline_variants":["Mirror symmetry splits four-waveguide junction into one quadrant","Quadrant solutions reconstruct full scattering matrix for junction","Eigenfunction matching solves junction via Neumann-Dirichlet quadrants","Symmetry reduces four-waveguide scattering to quarter-domain problem","Time-domain pulses from quadrant eigenfunction solutions at junction"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that, at each interface, matching the normal derivative over the open part (with the field implicitly zero over the rigid wall) is enough to enforce the true boundary condition; this implicit zero-extension is never stated and every computed coefficient depends on it.","fun_headline_variants_meta":{"raw":{"variants":["Mirror symmetry splits four-waveguide junction into one quadrant","Quadrant solutions reconstruct full scattering matrix for junction","Eigenfunction matching solves junction via Neumann-Dirichlet quadrants","Symmetry reduces four-waveguide scattering to quarter-domain problem","Time-domain pulses from quadrant eigenfunction solutions at junction"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000148,"raw_usage":{"total_tokens":1130,"prompt_tokens":824,"completion_tokens":306,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":440,"completion_tokens_details":{"reasoning_tokens":226}},"tokens_in":440,"tokens_out":306,"duration_ms":3200,"temperature":1.0,"reasoning_tokens":226,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:42:07.505901+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the normal derivative of the reconstructed Neumann-Neumann solution on the rigid-wall segment y in (a1,b1) at x=-b2 from the coefficients obtained by (39)-(44); if that derivative is not zero to within the truncation error, the implicit zero-extension of the Omega_1 field is not supplying the boundary condition stated in (14a).","supporting_citations":[{"cited_title":"ul Hassan, M","cited_arxiv_id":null,"evidence_quote":"Supplies the symmetry-and-scattering-matrix approach used here for time-domain waveguide scattering."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the mode-matching without root-finding technique that the quadrant solution relies on."},{"cited_title":"ul Hassan, M","cited_arxiv_id":null,"evidence_quote":"The multifurcated-duct mode-matching analysis that this paper extends to four waveguides."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Justifies the convergence of the discrete mode-matching truncation used in the numerical systems."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Precedent for using symmetry to solve wave scattering in a waveguide-resonator geometry."},{"cited_title":"Meylan, Fabien Montiel, and Sarah Wakes","cited_arxiv_id":null,"evidence_quote":"Provides the time-domain construction as a superposition of frequency-domain solutions via quadrature."},{"cited_title":"Wilks, M","cited_arxiv_id":null,"evidence_quote":"Supplies the generalized eigenfunction expansion method for two-dimensional acoustic time-domain scattering used here."}],"review_version":1}