{"id":"9d1f5d3a-09c1-4891-9c22-228fa927f368","arxiv_id":"2608.13367","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"An autocorrelation-tail metric makes time-domain topology optimization suppress in-band spectral ripple in broadband dielectric mirrors.","lead":"Optical designers who optimize broadband devices from short-pulse simulations can overlook narrow color-specific defects because they barely affect total reflected energy. This paper adds an autocorrelation-based penalty that suppresses such spectral ripple, yielding inverse-designed mirrors with flatter broadband response at no measurable loss in average reflectance.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Finite FDTD window may truncate the long-lag autocorrelation tail that T_AC is designed to measure, so very narrow spectral dips can evade the constraint by placing their tail beyond T_sim.","rationale":"The reader's conditional verdict is appropriate. The finite-window issue is an operational failure mode of the very premise the reader identified: that the long-lag autocorrelation tail is a faithful and sufficient proxy for in-band spectral defects. The paper has independent support in the form of open code, a 100-run ensemble for the constrained formulation, and binarization checks, so I do not recommend rejection. However, the central claim's generality for arbitrarily narrow in-band features is not fully established until the autocorrelation tail is shown to be completely captured, or until T_AC is shown to be insensitive to further increases in T_sim. A single rerun at longer simulation time would settle this concern without changing the current conditional verdict.","tokens_in":9991,"tokens_out":16125,"duration_ms":168409,"concrete_test":"Take the unconstrained optimized design from Fig. 4(a), which has a narrow spectral feature and T_AC = 8.628. Rerun the FDTD simulation with T_sim equal to 25,000, 50,000, 100,000, and 200,000 time steps, holding all other parameters fixed, and extend each run until the residual energy in the domain is below 1e-6 of the injected energy. Recompute eta_r, T_AC, and the reflectance spectrum near the dip for each T_sim. If T_AC and the dip depth or shape stabilize for longer windows and T_AC remains above 1, the reported detection is physical. If T_AC decreases substantially, for example toward or below 1, or if the dip shrinks, the constraint's apparent smoothing is an artifact of truncating the autocorrelation tail.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires Eq. (3) to integrate the true long-lag autocorrelation tail of Er(t). In the FDTD runs of Figs. 4-5 (Table II: 25,000 steps at 63.38 as, i.e., T_sim ~ 1.58 ps; frequency resolution ~ 0.63 THz), the autocorrelation is computed from a signal truncated at T_sim, so all lags tau > T_sim contribute zero. A spectral dip of width delta ~ 0.1 THz has an autocorrelation tail that decays on a timescale ~ 1/(pi*delta) ~ 3.2 ps, longer than the recording window. For such a defect, the computed E_AC underestimates the true value, and the constraint T_AC <= 1 can be satisfied not because the spectrum is smooth but because the offending tail lies beyond the recorded interval. The same truncation affects the objective eta_r in Eq. (4), so a high-Q resonance can be underrepresented in the objective and later appear as a narrow spectral feature. The paper asserts T_sim is long enough for the fields to decay but does not report a convergence check, such as residual energy or T_AC versus T_sim. Therefore, the demonstrated benefit of the constraint is not yet separated from the finite-window artifact for the very narrow features the metric is intended to catch.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript introduces a time-domain metric T_AC, defined as the weighted long-lag autocorrelation energy of the reflected field normalized by that of the incident broadband pulse (Eq. (3)), and argues that this quantity is sensitive to narrow in-band spectral features that conventional Parseval-type integrated objectives miss. The metric is incorporated into FDTD-based density topology optimization of 1D dielectric Bragg mirrors, either as a hard constraint (P_c, Eq. (6b)) or as a penalty term (Eq. (7)). Numerical results over 100 random initial designs show that the constrained ensemble exhibits visibly reduced sharp spectral features with essentially unchanged median reflectance (0.8096 vs. 0.8093), while the penalty formulation produces a Pareto trade-off and recovers apodized and chirped grating profiles (Sec. III). The paper includes FDTD and adjoint-method details in Appendix A and states that the code and dataset are publicly available.","tokens_in":10317,"tokens_out":9161,"duration_ms":98376,"significance":"If the central claim holds, the paper makes a useful contribution to time-domain inverse design, where quadratic integrated objectives are known to be blind to spectral ripple and localized defects. The Fourier derivation in Sec. II is clear and essentially correct, and the 100-run ensemble comparison is a genuine strength because it addresses the non-convexity of the design problem. The open code and dataset are also valuable for reproducibility. However, the finite FDTD recording window may compromise the very long-lag tails that T_AC is designed to measure, and the paper does not yet provide the convergence check needed to separate the effect of the constraint from a finite-window artifact. For this reason the demonstrated benefit is not yet fully established.","major_comments":[{"comment":"The finite FDTD recording window is not separated from the claimed effect of the constraint. For the runs in Figs. 4 and 5, T_Sim = 25,000 × 63.38 as ≈ 1.58 ps, so the FFT-based autocorrelation in Eq. (3) is evaluated only for lags |τ| < T_Sim and all longer lags contribute zero. A narrow spectral dip of width Δ ≈ 0.1 THz produces an autocorrelation tail that decays on a timescale of order 1/Δ ≈ 3.2 ps, longer than the window; such a defect would be essentially invisible to T_AC even though it is exactly the class of feature the metric is meant to detect. The same truncation affects η_r in Eq. (4), because late-time reflected energy is omitted. The statement in Appendix A that T_Sim is \"sufficiently long for the fields to decay\" is not supported by any reported convergence check. Please add a T_Sim-convergence study, for example by doubling T_Sim for representative constrained and unconstrained designs and reporting η_r, T_AC, and residual field energy, and verify that the sharp-feature suppression in Fig. 4(b.1) persists when the spectra are computed from longer runs.","section":"Appendix A, Table II and Eq. (3)"},{"comment":"The metric's parameter choices, specifically the 1% Hilbert-envelope threshold defining τ_min and the quadratic lag weight in E_AC, are introduced without any sensitivity or robustness analysis. Since the central claim is that T_AC is a faithful proxy for narrow in-band spectral defects, the paper should demonstrate that the reported smoothing and the T_AC = 1 bound are not artifacts of these particular values. A minimal check would be to repeat one constrained optimization with τ_min varied by a factor of two and with a linear instead of quadratic lag weight, and to report how the optimized spectra and T_AC values change.","section":"Section II and Appendix A"}],"minor_comments":[{"comment":"The text reports the unconstrained median T_AC = 1.854 but does not quote the median or spread of T_AC for the constrained ensemble; please report these values so the reader can confirm that the constraint is actually satisfied across the ensemble.","section":"Sec. III-A, Fig. 4(b.2)"},{"comment":"The expression \"Mnd − → 0\" appears to be a typographical error; it should read \"M_nd → 0\".","section":"Sec. III-B"},{"comment":"The notation is inconsistent between T_AC and TAC; please use a single subscripted form consistently throughout the manuscript.","section":"Eq. (3) and surrounding text"},{"comment":"Please state explicitly that E_AC is a heuristic estimator rather than an exact spectral-derivative norm, because the integral is restricted to |τ| > τ_min and uses a normalized autocorrelation; this would clarify the relationship between Eq. (3) and the Parseval-type identity that motivates it.","section":"Sec. II"},{"comment":"The fitted chirp parameters a1, a2, and D_0 are reported without uncertainties or residuals; adding these would make the comparison between the fitted values and the source-reference Bragg vectors G_hi and G_lo quantitative.","section":"Sec. III-C and Table I"},{"comment":"The Pareto plot would be easier to read if the axis labels and the meaning of the color or marker coding for different α values were explained in the caption.","section":"Fig. 5(a)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of IEEE TAP and the open-code policy is a plus. The revision should focus on the finite-window convergence check; without it, the central comparison between constrained and unconstrained ensembles remains ambiguous. The sensitivity of τ_min and the lag-weight exponent is a secondary concern that can be addressed with a small additional experiment."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this one before deciding what to do with it: it is a solid, incremental methods paper, not a field-reshaping result. The new thing is T_AC, a normalized weighted long-lag autocorrelation energy of the reflected field, introduced as a constraint or penalty in FDTD-based topology optimization. The underlying Wiener–Khinchin idea is textbook, and the authors say so. What is theirs is the specific definition, the normalization to the incident pulse, and the demonstration that it works on 1D Bragg mirror design.\n\nCredit where it is due. The math in Sec. II is clean, the 100-run ensemble comparison is a real step beyond the usual single-start demo, and the median reflectance is essentially unchanged while T_AC drops. They also ship code and data, which is what makes the paper reviewable. The chirp reconstruction in Sec. III-C is a nice bonus, though not load-bearing.\n\nNow the soft spots. The finite-window worry is real. With T_sim around 1.58 ps (25,000 steps at 63.38 as), any spectral dip narrower than roughly 0.5–1 THz has an autocorrelation tail that extends beyond the recorded window. The computed E_AC then undercounts that tail, so the constraint can be satisfied not because the spectrum is smooth but because the offending feature lives outside the measurement. The paper says T_sim is long enough for fields to decay but gives no convergence check, like E_AC or T_AC versus T_sim. In the demonstrated mirror case the visible feature is wide enough that the benefit still shows clearly, so this is a limitation of the method's claimed range, not a fatal flaw. I would ask for a supplementary T_sim sweep.\n\nTwo smaller things. The adjoint derivative of T_AC is described in prose, not written out; the code fills the gap but an explicit expression would help. And the weak-penalty conclusion in Sec. III-B rests on single-start runs, unlike the constraint result, which has the ensemble. That makes the Pareto story anecdotal. Both are minor.\n\nWho is this for: people doing time-domain inverse design in photonics, particularly topology optimization. It merits a serious referee, and in fact would benefit from one who knows signal processing and FDTD—to push on the truncation issue and on how generally the metric maps to spectral quality. Yes, send it out.","headline":"A useful, incremental tool paper: a time-domain autocorrelation metric that gives gradient-based broadband photonic design sensitivity to in-band ripple that energy-only objectives miss, deserving a normal peer-review round.","tokens_in":10812,"tokens_out":1977,"would_cite":true,"duration_ms":23066,"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":"This paper claims that a weighted long-lag autocorrelation metric turns narrow in-band spectral defects, which integrated energy objectives miss, into a visible time-domain penalty, and that constraining it yields broadband dielectric…","keywords":["autocorrelation metric","broadband reflectance","inverse design","topology optimization","FDTD","time-domain adjoint method","Bragg grating","spectral smoothness"],"falsifier":"Run the unconstrained optimization of Sec. III-A and re-evaluate $T_{AC}$ with progressively longer simulation times so a narrow phase-defect dip is resolved at finer frequency resolution; if the dip deepens while $T_{AC}$ stays near the reference value of 1, then the long-lag autocorrelation metric is insensitive to that class of spectral features.","tokens_in":9819,"feed_emoji":"🪞","tokens_out":8796,"duration_ms":80559,"temperature":0.7,"pith_summary":"Broadband inverse design in optics usually optimizes time-integrated energies, which by Parseval's theorem measure only the total response over a band and cannot see how that response is distributed in frequency. The paper proposes a time-domain metric $T_{AC}$, the weighted long-lag autocorrelation energy of the reflected field relative to the incident pulse, that is Fourier-dual to spectral sharpness and therefore exposes narrow dips, ripple, and sidelobes. Adding $T_{AC}\\le 1$ as a constraint, or a weak penalty, removes these defects in one-dimensional dielectric mirrors while leaving median reflectance essentially unchanged: $0.8096$ unconstrained versus $0.8093$ constrained over 100 runs. If this works as claimed, designers get a cheap way to enforce spectral flatness in time-domain inverse design without prescribing a target waveform or adapting the excitation pulse.","feed_headline":"Autocorrelation tails expose hidden spectral dips","feed_subtitle":"A time-domain metric turns invisible in-band ripple into a constraint, keeping broadband mirrors flat without losing reflectance.","key_machinery":"The central object is $T_{AC}[s;s_0] = E_{AC}[s]/E_{AC}[s_0]$, with $E_{AC}[s] = 2\\int_{\\tau_{\\min}}^{\\infty}(\\tau/\\tau_{\\min})^2 |c_s(\\tau)|^2\\,d\\tau$. It measures weighted long-lag autocorrelation energy of the reflected signal relative to the incident pulse; the reference delay $\\tau_{\\min}$ is the first lag where the Hilbert envelope of the reference autocorrelation has decayed to 1% of its peak. By the Fourier-dual relation between $\\partial_f S$ and $\\tau C(\\tau)$, this machinery converts sharp spectral features into long autocorrelation tails that the optimizer can see, quantify, and penalize either as a constraint or a weighted objective.","core_discovery":"The central discovery is that the long-lag tail of the autocorrelation of the reflected field, weighted by lag squared and compared with a broadband reference pulse, is a computable time-domain surrogate for narrow-band spectral nonuniformity. Because the energy spectral density $S(f)=|\\hat{s}(f)|^2$ and the autocorrelation $C(\\tau)$ form a Fourier pair, a localized spectral feature produces a slowly decaying oscillatory tail; the quadratic weight makes the metric sensitive to $\\partial_f S$, hence to sharp features and ripple. Adding $T_{AC}\\leq 1$ to the topology-optimization problem removes a phase-defect-like resonance in the reflectance spectrum, and across 100 random initializations the constrained ensemble keeps median reflectance at $0.8093$ versus $0.8096$ unconstrained while reducing $T_{AC}$ outliers. A weak penalty formulation with weighting $\\alpha=10^{-2}$ even reaches higher reflectance than the unconstrained problem while driving $T_{AC}$ near 1, and strong penalization recovers apodized and chirped grating morphologies whose near-unity reflectance matches a parametric linear-chirp reconstruction.","pith_inferences":["Editorial inference: because $T_{AC}$ is defined on the reflected signal alone, the same construction should apply to transmission or absorption spectra, making it a general-purpose spectral-flatness regularizer for any transfer function once a reference pulse is chosen.","Editorial inference: the reference delay $\\tau_{\\min}$ is derived from the incident pulse, so the metric's sensitivity can be tuned per spectral region; a natural extension is a multi-scale version with several threshold lags that targets ripple at different bandwidths simultaneously.","Editorial inference: the automatic emergence of a chirped grating suggests the metric could serve as an unsupervised probe for discovering structured dispersive designs, and comparing its Pareto front with hand-built apodized and chirped templates would isolate what the autocorrelation constraint adds beyond known grating recipes."],"forward_implications":["Maximizing the Parseval-type energy fraction alone is insufficient: one unconstrained optimization produced a reflectance spectrum with a narrow phase-defect-like resonance and $T_{AC}=8.628$ even at a local optimum.","Enforcing $T_{AC}\\le 1$ removes such defects without statistically reducing reflectance, with median $\\eta_r$ equal to 0.8093 against 0.8096 for the unconstrained ensemble.","A weak penalty with $\\alpha=10^{-2}$ in $J=(1-\\alpha)\\eta_r-\\alpha T_{AC}$ finds designs with higher reflectance than no penalty while keeping $T_{AC}$ close to 1.","Strong penalization trades reflectance for spectral flatness and yields apodized and chirped grating morphologies with near-unity flat reflectance, with $\\eta_r=0.978$ for the chirped design.","The same metric can be used as a post-evaluation diagnostic and can be carried into other broadband inverse-design problems, including near-field coupling, achromatic focusing, and broadband absorption."],"supporting_citations":[{"why":"Supplies the Parseval–Plancherel and Wiener–Khinchin identities that convert long-lag autocorrelation energy into sensitivity to spectral sharpness.","marker":"[10]"},{"why":"Provides the weighted integrated sidelobe level concept from signal processing that the metric adapts.","marker":"[11]"},{"why":"Identifies phase-shifted-grating phase-defect resonances as the class of narrow spectral features the metric targets.","marker":"[17]"},{"why":"Describes the adaptive spectral-weighting alternative whose need to modify the excitation pulse motivates the new approach.","marker":"[5]"},{"why":"Supplies the Bragg-grating design background, including apodized and chirped structures used as comparison.","marker":"[12]"},{"why":"Provides the time-domain adjoint method used to compute gradients of the objective and constraint with respect to the density.","marker":"[9]"},{"why":"Defines the FDTD discretization that simulates the Maxwell equations and evaluates the reflected field.","marker":"[14]"},{"why":"Shows the stronger pulse-shaping alternative that prescribes a target Hilbert envelope, against which the weaker autocorrelation constraint is contrasted.","marker":"[22]"}],"fun_headline_variants":["Autocorrelation constraint flattens broadband spectra","Time-domain metric kills spectral ripple in inverse design","Weighted autocorrelation smooths optical spectra","Inverse design via autocorrelation penalty removes ripple","Autocorrelation-aware inverse design polishes mirrors"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim stands on the assumption that the weighted long-lag autocorrelation tail of the reflected signal, measured against the incident pulse's coherence time, faithfully captures every narrow in-band spectral defect that matters, so that suppressing that tail is enough to guarantee a spectrally smooth response.","fun_headline_variants_meta":{"raw":{"variants":["Autocorrelation constraint flattens broadband spectra","Time-domain metric kills spectral ripple in inverse design","Weighted autocorrelation smooths optical spectra","Inverse design via autocorrelation penalty removes ripple","Autocorrelation-aware inverse design polishes mirrors"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000455,"raw_usage":{"total_tokens":2276,"prompt_tokens":923,"completion_tokens":1353,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":539,"completion_tokens_details":{"reasoning_tokens":1283}},"tokens_in":539,"tokens_out":1353,"duration_ms":10931,"temperature":1.0,"reasoning_tokens":1283,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:45:46.458169+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the unconstrained optimization of Sec. III-A and re-evaluate $T_{AC}$ with progressively longer simulation times so a narrow phase-defect dip is resolved at finer frequency resolution; if the dip deepens while $T_{AC}$ stays near the reference value of 1, then the long-lag autocorrelation metric is insensitive to that class of spectral features.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Parseval–Plancherel and Wiener–Khinchin identities that convert long-lag autocorrelation energy into sensitivity to spectral sharpness."},{"cited_title":"Sequence design to minimize the weighted integrated and peak sidelobe levels,","cited_arxiv_id":null,"evidence_quote":"Provides the weighted integrated sidelobe level concept from signal processing that the metric adapts."},{"cited_title":"Phase-shifted fiber bragg gratings and their application for wavelength demultiplexing,","cited_arxiv_id":null,"evidence_quote":"Identifies phase-shifted-grating phase-defect resonances as the class of narrow spectral features the metric targets."},{"cited_title":"Time-domain adjoint optimization for metalens design toward enhanced broadband efficiency and unifor- mity,","cited_arxiv_id":null,"evidence_quote":"Describes the adaptive spectral-weighting alternative whose need to modify the excitation pulse motivates the new approach."},{"cited_title":"Fiber bragg gratings,","cited_arxiv_id":null,"evidence_quote":"Supplies the Bragg-grating design background, including apodized and chirped structures used as comparison."},{"cited_title":"Time-domain topology optimization of power dissipation in dispersive dielectric and plasmonic nanostructures,","cited_arxiv_id":null,"evidence_quote":"Provides the time-domain adjoint method used to compute gradients of the objective and constraint with respect to the density."},{"cited_title":"Taflove and S","cited_arxiv_id":null,"evidence_quote":"Defines the FDTD discretization that simulates the Maxwell equations and evaluates the reflected field."},{"cited_title":"Topology optimization of pulse shaping filters using the hilbert transform envelope extraction,","cited_arxiv_id":null,"evidence_quote":"Shows the stronger pulse-shaping alternative that prescribes a target Hilbert envelope, against which the weaker autocorrelation constraint is contrasted."}],"review_version":1}