{"id":"f043fd8d-bcbe-4a64-a9f1-1ab0b79fb45f","arxiv_id":"2608.08490","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"For stochastic gravitational wave backgrounds containing simultaneous tensor, vector, and scalar modes, the more inclined LISA-TAIJIm network separates polarization components substantially better than LISA-TAIJIp, although Fisher forecasts for power-law spectra are broadly similar.","lead":"This paper compares two proposed LISA-TAIJI detector networks for their ability to separate the tensor, vector, and scalar polarization components of a stochastic gravitational wave background. It finds that the configuration with a larger angle between the two detector triangles, LISA-TAIJIm, suffers less from polarization confusion and can measure the components more cleanly.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (43) inverts the per-frequency Fisher matrix before integrating, but the stated model has global power-law amplitudes; the TVS factor-2–3 SNR claim may not survive an integrated Fisher-matrix recalculation.","rationale":"The reader focused on the absence of released ORF code. Reproducibility is a real issue, but the more targeted concern is internal to the SNR formalism. The qualitative geometric argument is standard and plausible: a larger relative inclination decorrelates polarization responses. The quantitative factor 2–3, however, rests on Eq. (43). The per-frequency inverse versus integrated-inverse distinction is a well-known pitfall in multi-parameter SGWB analyses when amplitudes are global. The paper's own Fisher analysis in Sec. IV uses global parameters and Eq. (47), so the integrated construction is available; applying it to the profiled SNR would settle the issue directly. I therefore recommend keeping a conditional verdict, with the condition expanded to include this SNR recalculation in addition to code or data release for the ORF computations.","tokens_in":25153,"tokens_out":19041,"duration_ms":229294,"concrete_test":"Recompute the TVS SNR maps of Fig. 7 with the frequency-integrated Fisher matrix for the three global amplitude parameters: define F_ab = Σ_ij 2T_obs ∫_0^{fmax} df [ (∂S_h,tot/∂A_a)(∂S_h,tot/∂A_b) / M_ij(f) ], with α fixed at 2/3, and set the profiled SNR for component P equal to A_P / sqrt[(F^{-1})_{PP}]. Then regenerate the three SNR-ratio panels of Fig. 7. If the LISA-TAIJIp/TAIJIm ratios move outside the stated 0.3–0.6 band by more than about 20%, the headline factor-2–3 advantage is an artifact of the per-frequency profiling convention; if the ratios remain within the band, the per-frequency formula is an acceptable numerical proxy.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The quantitative headline is the SNR ratio in Fig. 7: LISA-TAIJIp/TAIJIm ≈ 0.3–0.6 in the TVS case, i.e. a factor 2–3 in favor of LISA-TAIJIm. That number comes from Eq. (43), ρ²(P) = Σ_ij 2T ∫ [S_h^P(f)/σ_P(M(f))]² df, where σ_P is the diagonal element of the inverse of the frequency-by-frequency Fisher matrix F(f) in Eq. (33). The text in Sec. III.C says the SNR calculation treats the amplitudes A_P of the power-law model as free parameters and marginalizes over the remaining polarization amplitudes. If those amplitudes are global parameters common to all frequencies, the correct profiled SNR for component P is A_P² / [(∫ F(f) df)^{-1}]_{PP}, not ∫ (A_P t(f)/σ_P(f))² df. Integration and inversion do not commute for a non-diagonal matrix; the per-frequency expression discards the information that the nuisance amplitudes are identical across frequency bins, so it is not the profile likelihood for the stated global model. The discrepancy is configuration-dependent because the frequency structure of the ORF degeneracies differs between LISA-TAIJIp and LISA-TAIJIm, so the reported factor of 2–3 is not robust until this is checked. If Eq. (43) was intended as two-step nonparametric component separation followed by matched filtering, that interpretation should be stated explicitly; as written, Sec. III.C describes a global profiled likelihood for power-law amplitudes.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compares two proposed LISA-TAIJI network geometries, LISA-TAIJIp and LISA-TAIJIm, for detecting and separating tensor, vector, and scalar polarization components of an isotropic stochastic gravitational-wave background. Using the PD4L time-delay interferometry response model, the authors compute overlap reduction functions, effective sensitivities, power-law integrated sensitivities, profiled signal-to-noise ratios, and Fisher-matrix parameter forecasts. The central claim is that the larger relative constellation inclination of LISA-TAIJIm substantially weakens polarization degeneracies and yields recovered SNRs that are factors of roughly two to three larger than LISA-TAIJIp in the tensor-vector-scalar case, while the Fisher-forecast parameter constraints for the two geometries are broadly comparable.","tokens_in":25513,"tokens_out":7775,"duration_ms":91026,"significance":"The question addressed is timely and relevant for design studies of future space-based gravitational-wave networks. The paper's strength is its systematic application of the standard ORF and Fisher-matrix formalism to a concrete pair of mission geometries, and its explicit separation of model-independent component separation from model-dependent parameter estimation. If the headline numerical comparison survives scrutiny, the work would provide a useful design argument in favor of large-inclination LISA-TAIJI configurations for SGWB polarization studies. However, the main quantitative claim rests on an SNR formula that is not the profile likelihood for the stated global power-law model, and the numerical ORF pipeline is not independently checkable from the manuscript. The qualitative distinction between component separation and parameter estimation is likely robust, but the factor-of-2-3 SNR advantage requires correction and recomputation before it can be accepted.","major_comments":[{"comment":"Equation (43) does not compute the profiled SNR for the global power-law model stated in Sec. III.C. In that model A_P is a single amplitude common to all frequencies and the nuisance amplitudes A_Q are also global parameters. The correct profile for A_P is A_P^2 / [(2T ∫ W(f) df)^{-1}]_{PP}, where W_ab(f) is the per-frequency Fisher matrix for the global amplitudes, related to F(f) in Eq. (33). Equation (43), by contrast, evaluates ∫ [S_P(f)/σ_P(f)]^2 df with σ_P(f) = (F(f)^{-1})_{PP}; this is the SNR obtained by estimating A_P independently in every frequency bin and then combining the bins incoherently. Since integration and matrix inversion do not commute for a non-diagonal F(f), the two expressions differ, and the discrepancy is configuration-dependent because the frequency structure of the ORF degeneracies differs between LISA-TAIJIp and LISA-TAIJIm. The factor 2-3 ratio in Fig. 7 is therefore not the profile-likelihood ratio for the stated global model. Please recompute Fig. 7 with the integrated-Fisher profile, or explicitly redefine Eq. (43) as a two-step per-frequency component-separation SNR and adjust the claims accordingly.","section":"Sec. III.C, Eq. (43)"},{"comment":"The central quantitative results depend entirely on the numerically computed overlap reduction functions from the PD4L TDI response model, yet the manuscript provides no convergence checks, no comparison with known low-frequency or geometric limits, and no code or data release. Because the factor-2-3 advantage is the headline claim, a reader cannot distinguish a true geometric effect from an integration or pipeline artifact. Please add numerical validation, such as convergence with sky resolution, comparison with analytic zero-frequency ORF values, or release of the ORF computation code, before publication.","section":"Sec. II.B"},{"comment":"When the SGWB self-noise terms are retained, Eqs. (29) and (43) treat the nine cross-correlation baselines as statistically independent, each with variance M(f). That is not correct: the same SGWB realization contributes to all baselines, so the cross-spectra C_ij(f) have non-vanishing cross-covariances once the signal term in Eq. (29) is non-negligible. The likelihood in Eq. (31) implicitly requires the full 9x9 covariance matrix, not a per-baseline scalar. Since some plotted SNRs in Figs. 7-10 are O(100) at the largest amplitudes, the weak-signal limit is not valid over the full parameter space, and the per-baseline sum may overestimate or misorder the two networks. Please either restrict the SNR calculation to the weak-signal regime or implement the full multi-baseline covariance.","section":"Sec. III.C"}],"minor_comments":[{"comment":"The notation in Eq. (27) is ambiguous: the numerator should be written as the squared expectation of the cross-correlation estimator, and the denominator as its variance, e.g., ρ² = ⟨Ĉ⟩² / Var(Ĉ).","section":"Eq. (27)"},{"comment":"The conclusion contains a typo: 'nontrivial spectral features structures' should read 'nontrivial spectral features' or 'spectral structures'.","section":"Sec. V"},{"comment":"The Fisher analysis excludes all auto-correlation information within each mission by design; this should be stated more prominently, since including auto-correlations could change the relative constraints of the two geometries and the interpretation of the 'broadly comparable' result.","section":"Sec. IV"},{"comment":"The identical acceleration-noise assumption for LISA and TAIJI in Eq. (22) is a strong modeling choice; a brief justification or a sensitivity check would help the reader assess how much of the reported differences could come from this assumption.","section":"Sec. II.B"}],"recommendation":"major_revision","confidential_remarks":"The skeptic's concern about Eq. (43) is mathematically correct and should be the centerpiece of the revision. The paper's reliance on the authors' own prior papers is not a circularity problem, but the lack of ORF code or numerical validation is a practical reproducibility issue for a numerically driven claim. The distinction between component separation and parameter estimation in Sec. IV is a useful contribution and should be retained."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read: the paper's central distinction—network geometry matters more for model-independent component separation than for parametric Fisher forecasts—is well argued and likely correct. The per-frequency sensitivity curves (Figs. 3–6) cleanly show LISA-TAIJIm retaining more pure-mode sensitivity in the TVS case, and the paper is honest about the VS case reversing. That part is solid and is a genuinely new comparison.\n\nThe soft spot is Eq. (43). The text says the SNR is computed after marginalizing over the remaining polarization amplitudes, with the power-law amplitudes global across frequencies. But Eq. (43) inverts the per-frequency Fisher matrix F(f), takes the diagonal element, and integrates. That is the SNR for a per-frequency free amplitude, not for a global amplitude. For a single global amplitude, the profiled SNR is A_P^2 / [(∫F(f)df)^{-1}]_{PP}. Inversion and integration do not commute, and the discrepancy is frequency-structure dependent, so the reported ratio contours in Fig. 7 could shift. This is not fatal for the qualitative conclusion, because the per-frequency sensitivity analysis in Sec. III.B is independent of Eq. (43), but the abstract's \"factors of two to three\" rests on the questionable calculation. Either the authors should redo the SNR with the integrated Fisher matrix, or explicitly frame Eq. (43) as a two-step nonparametric separation followed by matched filtering.\n\nOther soft spots are more minor. The ORF computation code isn't released, so the load-bearing numerical chain (orbits, PD4L response, noise PSDs) can't be independently checked. No Monte Carlo validation of the inverse covariance, and assuming identical LISA/TAIJI acceleration noise is optimistic but standard. The self-citations for the alternative network configurations and PD4L are to published results, so that is not a problem.\n\nWho this is for: anyone making mission-geometry arguments for LISA-TAIJI or doing SGWB polarization forecasts. It deserves a serious referee. The referee should ask for the integrated-Fisher recalculation and code release, but the core insight is worth publishing.","headline":"The qualitative claim about LISA-TAIJIm's weaker polarization degeneracies is supported by the per-frequency analysis, but the headline factor-2–3 SNR advantage rests on a per-frequency Fisher inversion that is not the profile likelihood for the stated global power-law model.","tokens_in":26000,"tokens_out":2161,"would_cite":true,"duration_ms":23942,"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":"A larger relative inclination between the LISA and TAIJI constellations, LISA-TAIJIm, separates tensor, vector, and scalar stochastic background components with recovered signal-to-noise ratios up to two to three times higher than the…","keywords":["stochastic gravitational wave background","gravitational wave polarizations","LISA-TAIJI network","overlap reduction function","time-delay interferometry","polarization separation","component separation","Fisher matrix"],"falsifier":"An independent computation of the tensor, vector, and scalar overlap reduction functions using a different TDI combination or published orbital ephemerides could settle the claim: if the determinant of the 3 x 3 component-separation matrix for LISA-TAIJIp is within a few percent of that for LISA-TAIJIm, then the claimed factor-of-two-to-three SNR improvement cannot be correct.","tokens_in":1687,"feed_emoji":"🌌","tokens_out":9345,"duration_ms":146137,"temperature":0.7,"pith_summary":"Future space-based gravitational-wave observatories LISA and TAIJI could cross-correlate their measurements to detect the stochastic background and, crucially, to separate its tensor, vector, and scalar polarization components. The paper compares two proposed network geometries that differ only in the tilt of the TAIJI constellation plane relative to LISA's, changing the angle between the two planes from roughly 34.5 to 71 degrees. It finds that the larger-tilt configuration, LISA-TAIJIm, produces substantially weaker correlations among the three polarization responses, so its sensitivity to a mixed tensor-vector-scalar background degrades far less than that of the smaller-angle LISA-TAIJIp. In a full three-component fit, the recovered signal-to-noise ratios for each component are higher by factors of about two to three for LISA-TAIJIm. This advantage persists for tensor-scalar and tensor-vector mixtures but reverses for a pure vector-scalar background, so the geometry's benefit is specific to tensor-containing polarization mixtures.","feed_headline":"Tilting LISA-TAIJI wider separates polarizations 2-3x better","feed_subtitle":"The larger the angle between the LISA and TAIJI planes, the cleaner the tensor-vector-scalar decomposition of the stochastic background.","key_machinery":"The load-bearing object is the set of overlap reduction functions Gamma_T, Gamma_V, Gamma_S defined as sky-averaged cross-correlations of the detector responses to each polarization pair, computed with the PD4L TDI observable, a second-generation combination of laser links with an effective duration of 4L. These overlap reduction functions enter the component-separation matrix F = M^dagger $N^{{-1}}$ M, whose inverse gives the covariance of the recovered polarization amplitudes; off-diagonal terms encode the polarization degeneracies. The larger relative constellation inclination in LISA-TAIJIm makes the three overlap-reduction-function vectors more nearly orthogonal, and that increased orthogonality is the mechanism that reduces the variance of each recovered component and raises the recovered signal-to-noise ratio.","core_discovery":"The central claim is that the linear independence of the detector network's responses to different gravitational-wave polarizations is controlled by the relative inclination of the two constellations. Using the PD4L time-delay-interferometry response model, the authors compute overlap reduction functions Gamma_T, Gamma_V, Gamma_S for each polarization sector and assemble them into the component-separation matrix F = M^dagger $N^{{-1}}$ M. The off-diagonal entries of F measure polarization leakage; when the TAIJI plane is flipped from +60 to -60 degrees, the angle between the LISA and TAIJI planes grows from about 34.5 to about 71 degrees, the overlap-reduction-function vectors become more orthogonal, and the off-diagonal covariances drop. As a result, the variances of the recovered tensor, vector, and scalar amplitudes shrink substantially, and the profiled signal-to-noise ratios in a joint tensor-vector-scalar fit rise by factors of roughly two to three relative to the smaller-angle configuration. The same matrix explains why the vector-scalar case reverses: at the frequencies where the fiducial power-law spectrum contributes most to the SNR, the smaller-angle configuration has better individual component sensitivities, so the larger tilt is not universally beneficial. The paper also shows that a Fisher forecast for power-law spectral parameters yields nearly identical uncertainties for both geometries, because the model already treats the polarization sectors as distinct, reducing the role of geometric degeneracies.","pith_inferences":["Beyond the paper: the same matrix-orthogonality argument suggests that an even larger relative inclination than 71 degrees, or non-equilateral constellation shapes, could push tensor-vector-scalar separation further, though possibly at some cost to overall sensitivity; the paper does not explore this trade-off.","Beyond the paper: if real background spectra deviate from power laws and contain spectral features, the Fisher-forecast comparability would probably shrink, and the model-independent advantage of LISA-TAIJIm would become the dominant factor; the paper hints at this but does not quantify it.","Beyond the paper: a natural testable extension is to apply the same component-separation calculation to anisotropic or parity-violating backgrounds, where the larger inclination is already known to help; whether the polarization-separation gain persists when sky direction is also reconstructed remains open."],"forward_implications":["For a model-independent search that lets tensor, vector, and scalar amplitudes float simultaneously, the LISA-TAIJIm geometry should recover each component with roughly two to three times the signal-to-noise ratio of LISA-TAIJIp.","In mixed tensor-scalar and tensor-vector backgrounds, the larger-inclination network remains better but by a smaller margin; in a pure vector-scalar background, the smaller-inclination network can win by about a factor of 1.3 to 1.4 in SNR.","The sensitivity loss caused by simultaneously fitting several polarization sectors is largest at low frequencies, where the responses are hardest to distinguish, and the large-inclination geometry mitigates that loss.","Under a power-law model with known component spectra, both geometries deliver nearly equal parameter uncertainties; the geometric advantage shows up mainly when the polarization sectors are not assumed a priori.","These results give a concrete design criterion for future space-based networks: for polarization separation, maximize the relative inclination of the detector planes rather than the raw overlap."],"supporting_citations":[{"why":"Defines the tensor, vector, and scalar overlap reduction functions and the component-separation formalism this paper adopts.","marker":"[16]"},{"why":"Introduces the alternative LISA-TAIJI configurations, including the two constellation geometries compared here.","marker":"[47]"},{"why":"Computes the SGWB overlap reduction functions for LISA-TAIJI networks, establishing the cross-correlation behavior the present work extends to polarization separation.","marker":"[48]"},{"why":"Provides the PD4L second-generation TDI observable whose link responses enter every overlap reduction function.","marker":"[54]"},{"why":"Supplies the unified cross-correlation and maximum-likelihood component-separation framework used for the sensitivity and SNR calculations.","marker":"[63]"},{"why":"Derives the optimal-filter signal-to-noise ratio for stochastic background cross-correlation, the basis of Eq. (28).","marker":"[68]"},{"why":"Defines the power-law integrated sensitivity curves used to summarize broadband detectability.","marker":"[72]"},{"why":"Supplies LISA mission parameters and the acceleration-noise spectrum assumed for both LISA and TAIJI.","marker":"[32]"}],"fun_headline_variants":["Wider LISA-TAIJI angle improves polarization separation by 2-3x","LISA-TAIJIm tilt reduces degeneracy, triples SNR for polarizations","Large relative inclination in LISA-TAIJI unlocks cleaner polarization maps","Tilting TAIJI relative to LISA sharpens tensor-vector-scalar separation","How a wider LISA-TAIJI angle boosts stochastic background polarimetry"],"cache_read_input_tokens":28160,"weakest_assumption_plain":"The entire comparison rests on the numerical evaluation of the overlap reduction functions from the assumed orbits, TDI combination, and noise power spectral densities; if any of these inputs are inaccurate, the claimed factor-of-two-to-three geometric advantage could be misestimated.","fun_headline_variants_meta":{"raw":{"variants":["Wider LISA-TAIJI angle improves polarization separation by 2-3x","LISA-TAIJIm tilt reduces degeneracy, triples SNR for polarizations","Large relative inclination in LISA-TAIJI unlocks cleaner polarization maps","Tilting TAIJI relative to LISA sharpens tensor-vector-scalar separation","How a wider LISA-TAIJI angle boosts stochastic background polarimetry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000555,"raw_usage":{"total_tokens":2709,"prompt_tokens":1076,"completion_tokens":1633,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":692,"completion_tokens_details":{"reasoning_tokens":1527}},"tokens_in":692,"tokens_out":1633,"duration_ms":12058,"temperature":1.0,"reasoning_tokens":1527,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:33:45.628470+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An independent computation of the tensor, vector, and scalar overlap reduction functions using a different TDI combination or published orbital ephemerides could settle the claim: if the determinant of the 3 x 3 component-separation matrix for LISA-TAIJIp is within a few percent of that for LISA-TAIJIm, then the claimed factor-of-two-to-three SNR improvement cannot be correct.","supporting_citations":[],"review_version":1}