{"id":"66e495fa-0c45-462d-a6c8-5efc0952d4ff","arxiv_id":"2411.15733","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Using standard Rayleigh-Rice and ABC spectral fitting, the authors produce a model scattering profile for KAGRA's core mirrors: dP/dΩ ≈ 2.5×10^2 [1+(θ/7.1×10^-5)^2]^{-1.6} per steradian.","lead":"This paper fits a standard scattering model to topographic surface maps of KAGRA's four core mirrors and derives an angular profile for scattered laser light. The result suggests the previous design estimate for scattered-light noise at KAGRA was conservative, so the existing baffle design does not need immediate changes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central profile Eq. (17) is not validated for multilayer-coated mirrors; the Rayleigh-Rice single-surface relation Eq. (2) is assumed, not established, and the cited multilayer treatment [40] is not applied.","rationale":"The paper is a useful compilation of standard scattering theory with an engineering case study. The Abel-transform derivations in Section 2 are standard and the cross-check between the 1D and 2D PSD routes in Figs. 4 and 5 supports the circular-symmetry assumption. The quantitative outcome, however, is only as good as the bridge from measured topography to scattering, and that bridge, Eq. (2), is a single-surface Rayleigh-Rice result whose applicability to coated mirrors is explicitly uncertain in Section 4.2. The paper cites Ref. [40] as a relevant multilayer treatment but does not use it, so the central claim is conditional rather than established. The reader's weakest_assumption identifies the same concern, so I agree with the reader's assessment. The PSD extrapolation beyond 1 mm^-1 is a real but secondary issue because the steep theta^-3.2 falloff means the integrated scattered power is dominated by the fitted band below about 1 mrad; still, a high-frequency PSD measurement would tighten the model. The normalization inconsistency between Eq. (1) and Eq. (17) is a genuine flaw in presentation but does not change the numerical profile for ray tracing, since the profile is used as a differential scattering coefficient. Overall, the verdict CONDITIONAL remains appropriate: the estimate is useful and honestly qualified, but it should not be treated as the validated scattering profile of the coated mirrors until a multilayer calculation or a BRDF measurement confirms Eq. (17).","tokens_in":12047,"tokens_out":6081,"duration_ms":58268,"concrete_test":"Implement the multilayer scattering calculation of Ref. [40] for the actual KAGRA coating stack, or equivalently measure the BRDF at 1064 nm of a coated witness sample made with the same IBS coating over theta_s from about 30 micro-rad to 10 mrad. Feed the same measured 2D PSD used in the paper into that multilayer calculation and compare the predicted dP/dOmega_s with Eq. (17) over the fitted band 0.03-1 mm^-1, corresponding to theta_s from about 32 micro-rad to 1.06 mrad. If the multilayer-predicted or directly measured profile deviates from Eq. (17) by more than the fit uncertainties, roughly 10-15% in amplitude, then Eq. (17) is not representative for the coated KAGRA core mirrors and the quoted noise implications require revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, Eq. (17), is a scattering profile for the KAGRA core mirrors, but it is obtained by inserting a measured top-surface PSD into Eq. (2), a Rayleigh-Rice relation quoted from Ref. [13] for a single uncoated surface. The KAGRA mirrors are multilayer-coated, and multilayer coatings can scatter light through buried interfaces, coating microstructure, and bulk inhomogeneities that are not represented by the top-surface height PSD alone. Section 4.2 explicitly concedes: 'Due to the multilayer, it is not obvious whether the actual scattering angular profile would follow the relation in Eq. (2).' The paper cites Zeidler et al. [40], which provides a calculation method for multilayer-coated mirrors, but does not apply it. Thus the fitted parameters in Eq. (13) and the resulting profile Eq. (17) are an unvalidated single-surface surrogate for the true coated-mirror scattering. The conclusion that KAGRA's scattered-light noise is smaller than the previous 2e-6/theta^2 assumption depends directly on this surrogate. A secondary internal inconsistency is that Eq. (1) defines dP/dOmega_s as a PDF normalized to unit integral, while Eq. (17) integrates to roughly 7e-6; the quantity is actually a differential scattering coefficient normalized to incident power. That does not invalidate the numerical profile for ray tracing, but it makes the 'PDF' label and any downstream use that assumes unit normalization incorrect.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits the traditional Rayleigh-Rice method for converting a measured 2D topographic height map of a mirror surface into an angular scattering profile. It provides self-contained derivations of the Abel transform between 2D and 1D power spectral densities and of the ABC (K-correlation) spectral model, then describes a practical data-processing pipeline (crop, detrend, window, FFT, circular mean, curve fit). As a case study, the authors apply the pipeline to measured surface maps of the four KAGRA core mirrors and obtain a parametric scattering profile, Eqs. (17) and (18), which they compare to the previously assumed 2×10^-6/θ_s^2 model. They conclude that the updated profile predicts less scattered-light noise at the angles relevant to KAGRA's arm cavities, so no immediate changes to the baffle design are needed.","tokens_in":125,"tokens_out":6250,"duration_ms":118403,"significance":"If the derived profile is valid, it provides the first topography-based scattering model for the KAGRA core mirrors and a simple parametric input for ray-tracing simulations of stray light. The paper also collects and derives standard formulas (Abel transforms, ABC model conversion) in a clear way that may be useful to practitioners. A strength is that the data-processing flow is explicitly described, and the cross-check between the 1D and 2D fits supports the assumed circular symmetry of the surface roughness. However, the central result rests on two explicitly acknowledged assumptions—the applicability of the Rayleigh-Rice single-surface relation to multilayer-coated mirrors and the extrapolation of the PSD fit beyond the measured spatial-frequency band—and on an internal inconsistency in the normalization of the quantity called 'scattering PDF'.","major_comments":[{"comment":"Section 4.2 concedes that 'Due to the multilayer, it is not obvious whether the actual scattering angular profile would follow the relation in Eq. (2)' and cites Ref. [40] without applying it. Since Eq. (17) is obtained by inserting the top-surface PSD into Eq. (2), the central claim is an unvalidated single-surface surrogate for the actual coated-mirror scattering. The authors should either implement the multilayer scattering calculation of Ref. [40] to quantify the deviation, or explicitly state in the abstract and conclusion that the profile is conditional on the single-surface approximation; as written, the phrase 'first model-based estimation' overstates the evidence.","section":"Section 4.2"},{"comment":"The measured PSD data in Figs. 4 and 5 cover spatial frequencies up to only ~1 mm^-1 (θ_s ≲ 1 mrad), while the text notes that full arm-cavity coverage requires angles up to ~0.1 rad (~100 mm^-1). The profile Eq. (17) and the noise comparison in Section 4.1 therefore depend on extrapolating the fitted ABC/power-law model by more than two orders of magnitude in frequency. The paper should provide a sensitivity analysis or a conservative upper bound for the extrapolated region (e.g., the Lambertian diffusive-scattering floor mentioned in Section 4.2) and discuss how the conclusions in Section 4.1 would change if the high-frequency PSD deviates from the fitted form.","section":"Section 4.2 / Eq. (17)"},{"comment":"Eq. (1) defines dP/dΩ_s as a probability density with ∫ dP/dΩ_s dΩ_s = 1, but the fitted profile Eq. (17) integrates to roughly 7×10^-6 over the sphere (using the small-angle approximation appropriate for the KAGRA geometry), not unity. The quantity in Eq. (17) is therefore a differential scattering coefficient normalized to incident power, not a normalized PDF. The label 'scattering PDF' in Fig. 4 and the text is misleading, and any downstream ray-tracing use that assumes unit normalization will be incorrect. The authors should either renormalize the profile or explicitly present it as (1/P_i)dP_s/dΩ_s, state that its integral equals the scattered-power fraction represented by the fitted angular range, and provide guidance on how it should be used in ray-tracing simulations.","section":"Eq. (1) vs. Eq. (17)"}],"minor_comments":[{"comment":"The Introduction contains a typo: 'ray-tracging' should be 'ray-tracing'.","section":"Abstract / Introduction"},{"comment":"The sentence 'the read data are the measured 2D maps of the KAGRA mirrors' is awkward; consider rephrasing to 'the input data are the measured 2D maps'.","section":"Section 3.2"},{"comment":"The right vertical axis label 'Scattering PDF dP/dΩ_s (1/sr)' is inconsistent with the normalization issue raised in the major comments; consider using a neutral term such as 'differential scattering coefficient'.","section":"Fig. 4"},{"comment":"The fitted parameters A', B, C are quoted with uncertainties in Section 3.2, but Eq. (17) gives only central values; propagation of the parameter uncertainties into the profile would be helpful for users.","section":"Eq. (17)"},{"comment":"Reference [10] is an internal LIGO document that is not publicly accessible; consider citing a published equivalent or indicating how interested readers can obtain it.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for the journal and the derivations of Eqs. (11) and (13) appear correct. The central concern is whether the case-study result, Eq. (17), is sufficiently supported given the acknowledged coatings and extrapolation issues; these are fixable with additional analysis or clearer caveats. I also note the normalization inconsistency between Eq. (1) and Eq. (17), which should be corrected for the paper to be self-consistent. No concerns about citation or novelty were identified beyond the points already raised."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a clean write-up of standard Rayleigh-Rice/ABC scattering theory, plus a genuinely useful engineering result for KAGRA. The new content is the fitted scattering profile for the four KAGRA core mirrors, Eqs. (17)-(18), derived from measured topography. The derivations in Section 2 are correct, the data-processing flow is sensible, and the cross-check between the 1D and 2D fit routes is a good idea. The paper also openly lists its biggest weaknesses, which is more than most do.\n\nWhere I part company with the reader's more favorable take: the central profile is a fit to a worst-case envelope over frequencies up to ~1 mm^-1, then extrapolated to the ~0.1 rad angles the arm baffles care about. That is a big extrapolation with no data. More importantly, Eq. (2) is Rayleigh-Rice for a single uncoated surface, and the KAGRA mirrors are multilayer-coated. The paper concedes it is not obvious that Eq. (2) applies, but then still uses it as the basis of Eq. (17) and the conclusion that the old 2e-6/theta^2 assumption was conservative. The paper cites Zeidler's multilayer-scattering calculation but does not apply it. I don't think the result is wrong; I do think calling it representative overstates what can be defended.\n\nThere is also a labeling problem. Eq. (1) defines dP/dOmega_s as a probability density with unit integral, but the fitted profile integrates to roughly 1e-5, not 1. That inconsistency won't matter for ray tracing, where only the relative angular shape is used, but it will confuse readers and downstream users who rely on the PDF normalization.\n\nThe limitations section is honest about coating and bandwidth, but the abstract and conclusion still call the result representative. Data and code are not public, so the reproducibility burden falls on the fit parameters alone. The stray textual artifact and missing error bars on the final profile are minor but should be cleaned up.\n\nWho benefits: anyone at KAGRA, LIGO, or Virgo doing stray-light ray tracing, and readers who want a compact derivation of the Abel/inverse Abel relations. I'd bring it to a reading group if the topic is on the table, but it is not a method breakthrough.\n\nVerdict: I'd send this to peer review rather than desk reject. The paper needs a major revision: either an explicit statement that Eq. (17) is a provisional upper-envelope model with unquantified systematic error from coating and bandwidth, or an actual multilayer calculation—even approximate—following Zeidler. Fix the PDF normalization, and it becomes a solid methods note with a useful case study.","headline":"Useful engineering estimate of KAGRA mirror scattering, but the headline profile is an extrapolated worst-case envelope, not a validated prediction for coated mirrors.","tokens_in":12916,"tokens_out":2439,"would_cite":true,"duration_ms":22840,"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":"From measured surface maps of KAGRA's core mirrors, the paper derives a scattering profile that puts the arm's scattered-light noise below the design-phase estimate.","keywords":["optical scattering","Rayleigh-Rice theory","power spectral density","ABC model","surface roughness","gravitational-wave detectors","KAGRA","ray tracing"],"falsifier":"Measure the actual scattering distribution of one KAGRA core mirror with a goniometer over $\\theta_s$ from roughly $4\\times10^{-5}$ to $10^{-3}$ rad; if the measured profile falls off more slowly than $\\theta_s^{-3.2}$ or lies above Eq. (17) beyond the fitted uncertainties, the central claim fails. A second, less direct check is measuring the surface PSD up to $\\sim100\\ \\mathrm{mm}^{-1}$ with a microscope; if the power-law slope changes beyond the measured $\\sim1\\ \\mathrm{mm}^{-1}$ band, the extrapolation underlying the full-angle profile is not supported.","tokens_in":11742,"feed_emoji":"🪞","tokens_out":13365,"duration_ms":100985,"temperature":0.7,"pith_summary":"Using topographic height maps measured on the KAGRA core mirrors, this paper derives a ready-to-use angular profile for the light those mirrors scatter. The central quantitative claim is that the scattering probability per unit solid angle is $dP/d\\Omega_s \\simeq 2.5\\times 10^{2}\\,[1+(\\theta_s/(7.1\\times10^{-5}\\,\\mathrm{rad}))^2]^{-1.6}\\ \\mathrm{sr}^{-1}$, with an equivalent power-law form $dP/d\\Omega_s \\simeq 1.3\\times 10^{-11}/\\theta_s^{3.2}\\ \\mathrm{sr}^{-1}$. That profile is the first model-based scattering estimate for the KAGRA arm mirrors, and it places the scattered-light noise below the earlier $2\\times10^{-6}/\\theta_s^2$ assumption used for baffle design. The authors also organize derivations of the Rayleigh-Rice, Abel-transform, and spectral-model equations so the whole flow can be reused for other mirrors.","feed_headline":"KAGRA mirror maps yield scattering profile below old noise estimate","feed_subtitle":"The profile falls off as angle to the -3.2 power, below the earlier stray-light curve.","key_machinery":"The engine is the Rayleigh-Rice relation, Eq. (2), which links the two-sided 2D power spectral density $S_2(f_x,f_y)$ of surface height errors to the scattering probability density per unit solid angle. With normal incidence, small scattering angles, and $P_s\\simeq P_i$, it reduces to $dP/d\\Omega_s \\simeq 16\\pi^2 S_2(f)/\\lambda^4$, where $f\\simeq\\theta_s/\\lambda$. To turn measured maps into a model, the paper circularly averages the raw 2D PSD, fits an ABC/K-correlation form $S_2(f)=A'/[1+(Bf)^2]^{(C+1)/2}$, and converts the fitted PSD into the angular profile analytically. An independent branch compresses the 2D PSD to a one-sided 1D PSD through the Abel transform, fits the corresponding spectral model, and cross-checks the 2D result. Together these pieces carry the argument from raw interferometric maps to Eq. (17).","core_discovery":"The paper's central result is Eq. (17): after cropping, detrending, windowing, and circularly averaging the two-sided 2D power spectral density (PSD) of the four KAGRA core mirror maps, a worst-case ABC/K-correlation fit gives $dP/d\\Omega_s \\simeq 2.5\\times 10^{2}\\,[1+(\\theta_s/(7.1\\times10^{-5}\\,\\mathrm{rad}))^2]^{-1.6}\\ \\mathrm{sr}^{-1}$. An equivalent inverse-power-law limit, Eq. (18), is $dP/d\\Omega_s \\simeq 1.3\\times10^{-11}/\\theta_s^{3.2}\\ \\mathrm{sr}^{-1}$. These formulas describe an isotropic scattering envelope that bounds the actual mirror behavior for ray-tracing studies. Compared with the design-phase profile $2\\times10^{-6}/\\theta_s^2$, the new profile scatters less light at $\\theta_s > 48\\,\\mu\\mathrm{rad}$, and also less than the design curve between 37 and 48 $\\mu$rad, so the scattered-light noise in the KAGRA arm is expected to be smaller than previously estimated. The consistency between the independent 1D-PSD Abel-transform fit and the 2D-PSD fit supports treating the roughness as circularly symmetric.","pith_inferences":["If the Rayleigh-Rice relation extends to multilayer-coated mirrors, the same flow would give comparable scattering envelopes for other interferometric detectors from their metrology data, enabling direct cross-detector stray-light comparisons.","The unresolved high-spatial-frequency behavior is the most promising place to test the model; a microscopic measurement of the coated surface roughness would sharpen or refute the large-angle tail of the profile.","A consistency check the paper does not perform is integrating Eq. (17) over the full sphere and comparing the total scattered fraction with the independently measured 50-100 ppm optical loss of the KAGRA core mirrors.","The circular-symmetry assumption could be tested by examining the azimuthal dependence of the raw 2D PSD; coating-induced anisotropy would break the azimuthal invariance of the scattering profile."],"forward_implications":["KAGRA's scattered-light noise is lower than the design-phase estimate, so the existing baffle system needs no immediate change.","Eq. (17) can be fed directly into ray-tracing simulations to set upper limits on stray-light contamination inside the arm cavity.","The same preprocessing-and-fit flow can be applied to any smooth mirror with a topographic height map, including optics for future gravitational-wave detectors.","The power-law form of Eq. (18) provides a simple conservative envelope for worst-case noise studies."],"supporting_citations":[{"why":"Supplies the Rayleigh-Rice scattering relation and the spectral-model formulas that turn a surface PSD into an angular scattering profile.","marker":"[13]"},{"why":"Provides the measured topographic maps of the four KAGRA core mirrors used as the case-study input.","marker":"[8]"},{"why":"Establishes the scattered-light noise formalism and reciprocity relation used to connect the scattering profile to noise estimates.","marker":"[10]"},{"why":"Defines the 37 microradian angular cut derived from KAGRA's arm geometry, which sets the lower spatial-frequency fitting limit.","marker":"[12]"},{"why":"Source of the earlier inverse-square scattering profile based on metrology of previous interferometer optics, the baseline the new model is compared against.","marker":"[38]"},{"why":"Gives the Abel-transform relations and ABC/K-correlation model used to convert between 1D and 2D PSDs and fit the measured spectra.","marker":"[15]"},{"why":"Introduces the inverse power-law spectral model used for the conservative upper-limit scattering profile.","marker":"[20]"},{"why":"Discusses light scattering from multilayer-coated mirrors, the caveat for whether the uncoated-surface Rayleigh-Rice relation applies.","marker":"[40]"}],"fun_headline_variants":["KAGRA mirror scatter drops below old noise curve","Mirror roughness yields lower scattered light for KAGRA","New scattering model cuts KAGRA stray-light estimate","KAGRA mirrors scatter less than design curve predicts","Scattering profile falls as angle^-3.2, below old limit"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result stands on the assumption that the Rayleigh-Rice formula, developed for bare surfaces, still describes scattering from the multilayer-coated KAGRA mirrors; the paper itself flags that this is not obvious.","fun_headline_variants_meta":{"raw":{"variants":["KAGRA mirror scatter drops below old noise curve","Mirror roughness yields lower scattered light for KAGRA","New scattering model cuts KAGRA stray-light estimate","KAGRA mirrors scatter less than design curve predicts","Scattering profile falls as angle^-3.2, below old limit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000299,"raw_usage":{"total_tokens":1701,"prompt_tokens":889,"completion_tokens":812,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":505,"completion_tokens_details":{"reasoning_tokens":731}},"tokens_in":505,"tokens_out":812,"duration_ms":7222,"temperature":1.0,"reasoning_tokens":731,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:58:51.230641+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the actual scattering distribution of one KAGRA core mirror with a goniometer over $\\theta_s$ from roughly $4\\times10^{-5}$ to $10^{-3}$ rad; if the measured profile falls off more slowly than $\\theta_s^{-3.2}$ or lies above Eq. (17) beyond the fitted uncertainties, the central claim fails. A second, less direct check is measuring the surface PSD up to $\\sim100\\ \\mathrm{mm}^{-1}$ with a microscope; if the power-law slope changes beyond the measured $\\sim1\\ \\mathrm{mm}^{-1}$ band, the extrapolation underlying the full-angle profile is not supported.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Rayleigh-Rice scattering relation and the spectral-model formulas that turn a surface PSD into an angular scattering profile."},{"cited_title":"Characterization of core optics in gravitational-wave detectors: Case study of KAGRA sapphire mirrors,","cited_arxiv_id":null,"evidence_quote":"Provides the measured topographic maps of the four KAGRA core mirrors used as the case-study input."},{"cited_title":"Noise due to light scattering in interferometric gravitational wave detectors. I: Handbook of Formulae, and their Derivations,","cited_arxiv_id":null,"evidence_quote":"Establishes the scattered-light noise formalism and reciprocity relation used to connect the scattering profile to noise estimates."},{"cited_title":"Vacuum and cryogenic compatible black surface for large optical baffles in advanced gravitational-wave telescopes,","cited_arxiv_id":null,"evidence_quote":"Defines the 37 microradian angular cut derived from KAGRA's arm geometry, which sets the lower spatial-frequency fitting limit."},{"cited_title":"Compilation of Metrology Data for the LIGO Large Optics,","cited_arxiv_id":null,"evidence_quote":"Source of the earlier inverse-square scattering profile based on metrology of previous interferometer optics, the baseline the new model is compared against."},{"cited_title":"The prediction of BRDFs from surface profile measurements,","cited_arxiv_id":null,"evidence_quote":"Gives the Abel-transform relations and ABC/K-correlation model used to convert between 1D and 2D PSDs and fit the measured spectra."},{"cited_title":"Fractal surface finish,","cited_arxiv_id":null,"evidence_quote":"Introduces the inverse power-law spectral model used for the conservative upper-limit scattering profile."},{"cited_title":"Calculation method for light scattering caused by multilayer coated mirrors in gravitational wave detectors,","cited_arxiv_id":null,"evidence_quote":"Discusses light scattering from multilayer-coated mirrors, the caveat for whether the uncoated-surface Rayleigh-Rice relation applies."}],"review_version":1}