{"id":"81a74225-c170-4361-9e7e-cde0f60a8427","arxiv_id":"2506.23441","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"LISA could detect the post-Newtonian (3,3) harmonic from the nearest close white dwarf binaries, with detection probability up to 98% if LISA, Taiji, and TianQin combine data over 10 years.","lead":"This paper calculates whether LISA and other space-based gravitational-wave detectors can see faint extra signals from nearby white dwarf pairs, beyond the main hum. These tiny signals encode how uneven the two stars' masses are, which affects how the pairs evolve.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 98%/75% detection probabilities are conditional on the nearest CWDB having the high mass asymmetries of two EM-selected systems; without marginalizing over Delta and inclination, the 'reasonably likely' claim is not yet established.","rationale":"The paper is a clean analytical estimate: standard 0.5PN waveform, LISA noise curve, disk distribution. I checked the arithmetic: Table II follows from Eq. (14) and the stated rho2; the distances for rho3=5 (0.338 and 0.228 kpc) and the combined-detector boost sqrt(1^2+1^2+2^2)~2.5 reproduce the quoted 98% and 75% via Eq. (18). So the central mathematical chain is sound. The load-bearing weakness is that the 'detection probability' is computed for a hypothetical nearest source with fixed Delta and I, not for the actual distribution of nearest sources. Since the HH amplitude is proportional to Delta, and the nearest-distance distribution is steep, this prior dominates the answer. The reader identified this same assumption, and my independent check confirms it is the correct weak point. A population-synthesis marginalization is the natural, feasible test. If it passes, the conditional verdict can be upgraded; if it fails, the abstract's 'reasonably likely' would need qualification. I do not see a basis for rejection, because the model dependence is stated in the text and the numbers are internally consistent.","tokens_in":9382,"tokens_out":8912,"duration_ms":97861,"concrete_test":"Run a Monte Carlo over the population: draw ~2000 f2>=4mHz Galactic CWDBs from a published population-synthesis catalog (e.g., Lamberts et al. 2019 or Toubiana et al. 2024), assign distances from the disk model of Eq. (16), inclinations from p(I) proportional to sin(I), and masses from the catalog; identify the nearest source and compute rho3 from Eq. (14) with 10-yr LISA/Taiji/TianQin sensitivities. Count the fraction of realizations with rho3>=5. If that fraction is substantially below 98%/75% (e.g., <50%), the headline claim is not robust to the Delta-prior assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equations (17)-(18) give the distribution of the minimum distance dmin among N=2000 f2>=4mHz CWDBs, but the detection probabilities quoted for LISA-Taiji-TianQin are then evaluated by setting the nearest source's parameters equal to those of HM Cancri (Delta=0.34, I=38deg) or ZTF J1539 (Delta=0.49, I=84deg) from Table I. The HH SNR is rho3 ~ (beta Delta E3(I)/E2(I)) rho2, so the distance at which rho3=5 scales as Delta E3(I)/E2(I); with F(<d)~d^2.5, the detection probability scales roughly as [Delta E3(I)/E2(I)]^2.5 in the small-probability regime. Halving Delta to 0.17 lowers the HM-Cancri-like combined 10-yr probability from ~98% to ~48%, and halving the ZTF-like Delta lowers ~75% to ~20%. The paper does not average over the unknown masses (or the published HM Cancri mass ambiguity) or over isotropic inclinations, so the headline numbers are an upper-tail scenario rather than a marginalized detection probability. The rest of the derivation is internally consistent; the weakness is the unstated conditioning on Table I's asymmetry values.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper estimates the detectability of 0.5PN higher-harmonic signals, mainly the (3,3) mode at f3 = 3 f_o and the f1 = f_o mode, from circular close white dwarf binaries in the mHz band. It uses published PN pattern functions and the LISA noise curve to express the harmonic SNRs rho3 and rho1 in terms of the quadrupole SNR rho2, the mass-asymmetry parameter Delta, and the inclination-dependent intensity functions E3 and E1. A Galactic disk model is then used to estimate the minimum distance d_min of f2 >= 4 mHz CWDBs, giving d_c ~ 0.5 kpc for N = 2000. For HM Cancri-like and ZTF J1539-like parameters placed at d = 0.5 kpc, the paper finds LISA-only 10-year detection probabilities (rho3 >= 5) of 31% and 13%, rising to 98% and 75% when LISA, Taiji, and TianQin data are coherently combined. The paper also discusses separation from eccentricity-induced harmonics, anomalous polarization tests, EM follow-up, and measurement of the two masses.","tokens_in":9637,"tokens_out":7438,"duration_ms":86282,"significance":"If the detection-probability claim holds, the paper identifies a genuinely new and observable channel: the 0.5PN (3,3) harmonic of nearby CWDBs, whose amplitude directly probes mass asymmetry and can break the chirp-mass degeneracy when combined with the quadrupole amplitude. The derivation is transparent and internally consistent, and the analytic scaling relations, especially rho3 = rho2 beta Delta E3/E2 and the F(<d) ~ d^2.5 population estimate, are useful tools. The paper does not introduce new waveform machinery or fit to a target result; the stress-test concern about circularity does not land. However, the headline probabilities are conditional on fixed representative values of Delta and inclination, and the abstract overstates the role of LISA alone, so the significance is real but more fragile than the current presentation suggests.","major_comments":[{"comment":"The 98% and 75% detection probabilities are conditional on the specific mass-asymmetry and inclination values of HM Cancri (Delta = 0.34, I = 38 deg) and ZTF J1539 (Delta = 0.49, I = 84 deg) from Table I. Because rho3 is proportional to Delta E3(I) and F(<d) is approximately proportional to d^2.5, the detection probability scales roughly as [Delta E3(I)/E2(I)]^2.5 in the small-probability regime. Halving Delta from 0.34 to 0.17 reduces the combined HM-Cancri-like probability from about 98% to about 50%, and halving the ZTF-like value reduces its probability from about 75% to about 20%. The paper does not marginalize over the mass-ratio distribution or over an isotropic inclination distribution, nor does it propagate the published mass ambiguity for HM Cancri. The claim that LISA is 'reasonably likely' to detect these modes therefore needs either a population-averaged estimate or a prominently stated conditional interpretation; as written, the headline numbers represent an upper-tail scenario rather than an expected detection probability.","section":"Nearby binaries; Eq. (18) and Table II"},{"comment":"The abstract and introduction say that 'LISA is reasonably likely' to detect higher-harmonic signals, but the quantitative support for the high probabilities comes from coherently combining LISA, Taiji, and TianQin, which gives an SNR gain factor of about sqrt(1^2+1^2+2^2) ~ 2.5. LISA alone yields only 31% and 13% for the two representative systems. This distinction should be made explicit in the abstract and conclusions; otherwise the reader is left with the impression that the 98%/75% numbers are LISA-only results.","section":"Abstract and Summary vs. Nearby binaries"},{"comment":"The estimate d_c ~ 0.5 kpc uses N = 2000 for the number of f2 >= 4 mHz CWDBs, but the text states that LISA is expected to detect about 2000 detached CWDBs and a similar number of AMCVn binaries in this band. If the relevant population is taken to be N ~ 4000, d_c decreases by a factor 2^0.4 ~ 1.32 and the quoted detection probabilities increase. The choice N = 2000 is conservative, but the paper should state this explicitly and ideally show how the probabilities depend on N, since the total population of f2 >= 4 mHz CWDBs is not known to a factor of two.","section":"Nearby binaries, Eq. (17)"}],"minor_comments":[{"comment":"The notation P(< dmin) is ambiguous; it should be written as the cumulative distribution function P(dmin < x) = 1 - exp[-N F(<x)] to make clear that the right-hand side is a function of the argument x, not of the random variable dmin.","section":"Eq. (18)"},{"comment":"The figure shows four combinations of (R_D, z_D) but has no legend; the caption lists the combinations without indicating which curve corresponds to which combination. Adding a legend or labeling the curves directly would improve readability.","section":"Fig. 1"},{"comment":"The SNR estimate uses the effective noise spectrum Sn(f) and the 0.89 averaging factor for linearly polarized radiation, but this is stated only in the text around Eq. (6). Since Table II is the main quantitative result, the paper should state explicitly that the tabulated rho_i values include that factor and use the extended 10-year mission.","section":"Eq. (6) and Table II"},{"comment":"The approximation Sn(f3)/Sn(f2) = 1 is reasonable in the flat part of the LISA noise curve, but the f1 mode falls at 3.1 mHz for HM Cancri and 2.4 mHz for ZTF J1539, where the noise curve is no longer flat. The paper acknowledges this, but it would be helpful to quote the actual rho1 values including the real noise curve rather than only the ratio.","section":"Section on PN corrections, text near Eq. (14)"}],"recommendation":"major_revision","confidential_remarks":"The core derivation is sound and the paper is clear and well-scoped. The main revision should center on properly conditioning the detection probabilities on Delta and inclination, and on reconciling the abstract's LISA-only language with the three-detector combination. If the authors add even a simple lognormal or flat mass-ratio marginalization and state the conditional nature of the quoted numbers, I would be supportive. The stress-test concern about conditional probabilities is the key unresolved point; the circularity concern is not."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the paper does something genuinely useful: it points out that the (3,3) post-Newtonian harmonic, normally ignored for LISA-band white dwarf binaries because beta ~ 1e-3, can actually be detectable for the nearest sources, and it carries a direct handle on mass asymmetry. The algebra in Eqs. (13)-(15) is a clean reduction of published 0.5PN pattern functions to a simple SNR scaling, and the numbers check out. The discussion of separating eccentricity-induced f3 contamination via apsidal precession is a nice, careful touch.\n\nSecond, the headline probabilities (98% and 75% for LISA+Taiji+TianQin) are conditional on an assumption that is load-bearing but not fully flagged: the nearest f2 >= 4mHz binary is assigned the mass asymmetry and inclination of either HM Cancri (Delta=0.34, I=38deg) or ZTF J1539 (Delta=0.49, I=84deg). The paper does not marginalize over these unknown parameters. Because the detection distance scales as Delta*E3(I)/E2(I) and the nearest-source distribution scales as d^2.5, the probability scales roughly as Delta^2.5. Halving Delta cuts the HM-Cancri-like probability from 98% to about 50%, and the ZTF-like from 75% to about 20%. The paper does note the HM Cancri mass ambiguity (0.55+0.27 vs 1.0+0.17 Msun) but does not propagate that spread into the detection estimates. So the 'reasonably likely' claim is really 'reasonably likely if the nearest source happens to look like these two known systems.' That is an upper-tail scenario, not a marginalized expectation.\n\nThe rest of the paper is solid. The spatial distribution estimate for d_c is standard, the disk-scale-length sensitivity is shown in Fig. 1, and the SNR values in Table II follow from the stated formulas. The combined-detector factor of ~2.5 assumes coherent stacking of LISA+Taiji+TianQin with design sensitivities; that is optimistic but not absurd, and it is clearly stated. The abstract slightly overstates LISA-only prospects, since LISA alone gives 31% and 13% for the two models.\n\nWho should read this: anyone planning LISA source catalogs or thinking about what can be learned from individual bright CWDBs. It is a good target-study paper, not a population synthesis. A serious referee should ask for a sensitivity plot over Delta and I, or at least a toy marginalization, before accepting the strong probability claims; the core qualitative result - that HH modes from nearby CWDBs are worth searching for - will survive that exercise.\n\nI'd send it to review. It is a legitimate, clearly written contribution with a fixable gap between the headline number and the evidence backing it.","headline":"Clean target study showing the (3,3) harmonic is within reach for the nearest CWDBs, but the headline detection probabilities are conditional on the nearest source having the large mass asymmetry of two EM-selected systems.","tokens_in":10203,"tokens_out":4349,"would_cite":true,"duration_ms":39001,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"LISA likely sees the (3,3) wave mode in nearby white dwarfs","keywords":["gravitational waves","white dwarf binaries","higher harmonics","post-Newtonian","LISA","AM CVn binaries","mass asymmetry"],"falsifier":"After ten years of combined LISA, Taiji, and TianQin data, look for the $(3,3)$ line in the nearest $f_2 \\ge 4$ mHz binary identified by electromagnetic surveys; if the binary lies within $0.3$ kpc with $\\Delta \\sim 0.4$ and no line appears at the predicted amplitude, the amplitude scaling or noise model is wrong. Alternatively, if an electromagnetic census shows the nearest such binary has $\\Delta < 0.2$, the 98% and 75% detection probabilities should be recomputed downward.","tokens_in":9129,"feed_emoji":"🔭","tokens_out":7537,"duration_ms":70627,"temperature":0.7,"pith_summary":"The paper argues that LISA has a good chance of detecting higher-harmonic gravitational-wave signals, particularly the $(l,|m|)=(3,3)$ mode, from a few nearby close white dwarf binaries even though their orbital velocities are only about $10^{-3}$ times the speed of light. The amplitude of this post-Newtonian mode is proportional to the binary's mass asymmetry, so a detection would give a direct, robust probe of how unequal the two white dwarf masses are, which matters for understanding mass-transferring AM CVn systems. The known template systems, HM Cancri and ZTF J1539, are too faint at their reference distances, but the paper estimates the nearest undiscovered binary sits near $0.5$ kpc, where the (3,3) mode becomes detectable. Combining LISA with Taiji and TianQin for ten years raises the detection probabilities to 98% and 75% for HM Cancri-like and ZTF J1539-like binaries, respectively. The claim makes the search a concrete observational program rather than a theoretical curiosity.","feed_headline":"LISA likely sees the (3,3) wave mode in nearby white dwarfs","feed_subtitle":"Close binaries should reveal mass asymmetry through post-Newtonian harmonics.","key_machinery":"The central object is the $(l,|m|)=(3,3)$ mode of the spin-weighted spherical harmonic decomposition, whose waveform pattern functions enter at 0.5PN order as $H_+^{(1/2)} \\propto \\beta\\Delta$. The argument's engine is the signal-to-noise formula $\\rho_3 \\simeq \\rho_2 \\beta \\Delta E_3(I)/E_2(I)$, where $\\rho_2$ is the loud quadrupole-mode signal-to-noise ratio and $E_3/E_2$ are inclination-dependent intensity functions; because $\\rho_2 \\sim 1000$ for the nearest binaries, the small factors $\\beta\\Delta \\sim 0.001$ are overcome. The probability estimate then uses the cumulative distance distribution $F(<d) \\propto d^{2.5}$ for the Galactic disk and the Poisson formula $P(<d_{\\min}) = 1 - \\exp[-N F(<d_{\\min})]$ to argue that the nearest $f_2 \\geq 4$ mHz binary is likely within $\\sim 0.5$ kpc.","core_discovery":"At its core the paper claims that the $(3,3)$ post-Newtonian harmonic, usually neglected for white-dwarf binaries because the PN parameter $\\beta$ is only $\\sim 0.004$, will be measurable in a handful of nearby systems. The mode amplitude is governed by the scaling factor $s = \\beta \\Delta$, where $\\Delta = (m_a - m_b)/m$ is the mass-asymmetry ratio, so measuring it yields mass information that the dominant quadrupole mode cannot provide. With a single LISA detector operating for ten years, an HM Cancri-like binary at $d = 0.5$ kpc gives $\\rho_3 \\simeq 3.4$, just below the detection threshold of $\\rho_3 \\geq 5$; combining LISA, Taiji, and TianQin boosts the signal-to-noise ratio by a factor of $\\sim 2.5$ and makes the detection probabilities 98% and 75% for the two representative templates. This is the paper's central result: the nearest binary, not the known ones, is the key target, and it should be close enough for the harmonic to show up in coordinated long-term observations.","pith_inferences":["If the detection succeeds, the measured $\\Delta$ values for nearby AM CVn binaries will bear directly on the classical stability criterion $\\Delta \\geq 1/5$, offering an empirical test of mass-transfer stability.","The same signal-to-noise framework could be applied to the $f_1$ mode to search for dipole gravitational radiation, since the paper shows the PN contribution can be predicted and subtracted using the strong quadrupole mode.","The probability estimate depends on the assumed disk scale lengths; future electromagnetic surveys that census the nearest $f_2 \\geq 4$ mHz binaries will directly test whether the $d^{2.5}$ scaling near 0.5 kpc is correct.","A null detection after ten years of combined data would not only challenge the noise and distance assumptions but would also tighten the upper bound on typical mass asymmetry of nearby white-dwarf binaries."],"forward_implications":["For several nearby close white dwarf binaries, the $(3,3)$ harmonic will be detectable, providing a direct measurement of the mass-asymmetry ratio $\\Delta$.","Coordinated 10-year operation of LISA, Taiji, and TianQin improves the harmonic signal-to-noise ratio by about a factor of 2.5, raising detection probabilities to 98% and 75% for the two representative systems.","The $f_1$ orbital-frequency mode can be separated from eccentricity-induced and anomalous-polarization signals, giving tests of gravitational-wave polarization and alternative gravity theories.","A nearby binary detected in harmonics becomes a prime electromagnetic follow-up target, with the gravitational-wave data predicting periodic variations to search in light curves and spectral lines.","Combining the harmonic amplitude with an electromagnetic parallax distance and the quadrupole amplitude can solve for the two component masses separately."],"supporting_citations":[{"why":"Supplies the expected number of Galactic CWDB detections above 4 mHz and the frequency-drift observational target.","marker":"[7]"},{"why":"Identifies the currently known EM-selected CWDBs above 4 mHz and notes the prospect of pre-launch identification of nearby binaries.","marker":"[14]"},{"why":"Provides the LISA noise curve used to set $S_n(f_3)/S_n(f_2)=1$ and to compute signal-to-noise ratios.","marker":"[21]"},{"why":"Supplies HM Cancri's basic parameters and quadrupole signal-to-noise ratio used as a template system.","marker":"[26]"},{"why":"Supplies ZTF J1539's parameters, distance, and quadrupole signal-to-noise ratio used as the second template.","marker":"[24]"},{"why":"Provides the distance lower bound and mass-model uncertainty for HM Cancri, affecting distance and mass-asymmetry assumptions.","marker":"[23]"},{"why":"Derives the post-Newtonian waveform pattern functions and quadrupole formula for circular binaries.","marker":"[2,3]"},{"why":"Gives Taiji's planned sensitivity used in the combined-detection signal-to-noise estimate.","marker":"[8]"},{"why":"Gives TianQin's design sensitivity used in the combined-detection signal-to-noise estimate.","marker":"[9]"}],"fun_headline_variants":["LISA may catch (3,3) harmonic from nearby white dwarfs","Mass asymmetry via (3,3) mode in nearby binaries","Close white dwarfs could reveal (3,3) wave with LISA","Post-Newtonian (3,3) signal possible for nearest systems"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The high detection probabilities rest on the premise that the nearest undiscovered $f_2 \\geq 4$ mHz binary has a mass-asymmetry ratio $\\Delta$ close to the 0.34 and 0.49 values of the two known template systems and a favorable inclination; if the typical nearby binary had $\\Delta \\sim 0.17$, the probabilities would drop by roughly a factor of six.","fun_headline_variants_meta":{"raw":{"variants":["LISA may catch (3,3) harmonic from nearby white dwarfs","Mass asymmetry via (3,3) mode in nearby binaries","Close white dwarfs could reveal (3,3) wave with LISA","Post-Newtonian (3,3) signal possible for nearest systems"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000207,"raw_usage":{"total_tokens":1378,"prompt_tokens":904,"completion_tokens":474,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":520,"completion_tokens_details":{"reasoning_tokens":395}},"tokens_in":520,"tokens_out":474,"duration_ms":5120,"temperature":1.0,"reasoning_tokens":395,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:42:28.977322+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"After ten years of combined LISA, Taiji, and TianQin data, look for the $(3,3)$ line in the nearest $f_2 \\ge 4$ mHz binary identified by electromagnetic surveys; if the binary lies within $0.3$ kpc with $\\Delta \\sim 0.4$ and no line appears at the predicted amplitude, the amplitude scaling or noise model is wrong. Alternatively, if an electromagnetic census shows the nearest such binary has $\\Delta < 0.2$, the 98% and 75% detection probabilities should be recomputed downward.","supporting_citations":[{"cited_title":"Expanding the ultracompacts: gravitational wave-driven mass transfer in the shortest-period binaries with accretion disks","cited_arxiv_id":"2411.12796","evidence_quote":"Identifies the currently known EM-selected CWDBs above 4 mHz and notes the prospect of pre-launch identification of nearby binaries."},{"cited_title":"Paczy´ nski, Acta Astronomica17, 287 (1967)","cited_arxiv_id":null,"evidence_quote":"Supplies HM Cancri's basic parameters and quadrupole signal-to-noise ratio used as a template system."},{"cited_title":"In Table I, we present the signal-to-noise ratios ρ2 for the quadrupole mode estimated in [26] and [24] for a 4yr operation of LISA","cited_arxiv_id":null,"evidence_quote":"Supplies ZTF J1539's parameters, distance, and quadrupole signal-to-noise ratio used as the second template."},{"cited_title":"Munday, T","cited_arxiv_id":null,"evidence_quote":"Provides the distance lower bound and mass-model uncertainty for HM Cancri, affecting distance and mass-asymmetry assumptions."},{"cited_title":"Hu and Y.-L","cited_arxiv_id":null,"evidence_quote":"Gives Taiji's planned sensitivity used in the combined-detection signal-to-noise estimate."}],"review_version":1}