{"id":"e1ab67e7-d2bb-4383-bd03-94c9d5f00cf1","arxiv_id":"2607.28920","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A simulation-based line shape code that alternates between atomic and molecular bases yields broader quasi-H2+ satellites that better match a white dwarf's ultraviolet spectrum.","lead":"This paper introduces a new way to compute spectral line shapes in white dwarf star atmospheres by switching between atomic and molecular descriptions when atoms collide. The resulting hydrogen Lyman profiles with quasi-molecular satellites match an ultraviolet spectrum of a white dwarf better than the standard theory.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Infinite-separation overlap approximation in Eq. (33) is untested at small R and could distort the quasi-molecular satellites that drive the claimed spectral agreement; a recomputation with R-dependent overlaps is needed.","rationale":"The reader identified the infinite-separation overlap approximation in Sec. 4.3 as the weakest assumption, and I agree. The central claim is not just that the method works in principle, but that the resulting line shapes are broader than ULBT and agree with observed FUSE spectra. The quasi-molecular satellites are wing features originating from close collisions at R ~ 10 a0. Eq. (33) projects the molecular dipole into an atomic basis using infinite-separation overlaps at all R, although the paper explicitly says this is invalid at small R. Because the projection determines how much quasi-molecular opacity is assigned to each atomic channel, it directly shapes the satellite profiles used in the Wolf 1346 comparison. The paper tests several other approximations (screening, velocity distribution, simultaneous broadening, critical radius), but it does not test the R-dependence of the overlaps in Eq. (33). That is a gap between the method's physical claims and its numerical implementation. The proposed test—rerunning one representative profile with R-dependent overlaps—would settle whether this approximation is benign or whether it is partly responsible for the claimed broadening. The reader's CONDITIONAL verdict remains appropriate: the paper is promising and internally consistent, with a useful validation in Fig. 6 showing it can reproduce ULBT when the relevant approximations are imposed, but the untested dipole projection prevents full acceptance. I recommend no change to the reader's verdict because the identified concern is the same one already used to justify the conditional status, and the needed test is clearly specified.","tokens_in":23584,"tokens_out":5497,"duration_ms":62280,"concrete_test":"Recompute the T = 10,000 K, n_e = 10^17 cm^-3 ion-only Lyα line shape with Eq. (33) replaced by explicit R-dependent overlaps ⟨a|m(R)⟩, using the same finite n=1–3 atomic basis and the Zammit H2+ molecular wavefunctions, keeping r_crit = 45 a0 and all other Xenomorph settings unchanged. Compare the 1400 Å satellite peak strength, width, and red-wing slope against the infinite-overlap result. If the satellite changes by more than the factor-of-two agreement margin quoted in Sec. 5.1, the infinite-separation dipole projection is a load-bearing distortion; if it is unchanged, the approximation is benign and the central claim is substantially strengthened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that the new multi-basis method produces quasi-H2+ Lyman features that are broader than ULBT and agree with the Wolf 1346 FUSE spectrum. The load-bearing assumption is the infinite-separation overlap approximation used in Eq. (33) to transform the time-dependent molecular dipole D_m(t) into the atomic basis at every instant of a close collision. The paper concedes in Sec. 4.3 that this approximation 'is no longer valid as R can be arbitrarily small,' because a molecular state at small R overlaps many atomic states and exact R-dependent overlaps would 'split the quasi-molecular opacity across countless different atomic line profiles.' This matters because the 1060/1080 Å Lyβ and 1400 Å Lyα satellites are formed at R ~ 10 a0 (Fig. 1), not at infinite separation. The strength, width, and shape of these satellites are set by how the molecular dipole is projected onto atomic channels before the Fourier transform in Eq. (11). Using infinite-R overlaps at small R could systematically mis-assign oscillator strength among atomic final states, producing broader or misplaced features that are an artifact of the projection rather than of the physical multi-basis dynamics. The U(t) basis change at r_crit (Eq. 26) is justified because R is large there, but Eq. (33) applies the same infinite-R overlaps throughout the collision. No test in Sec. 5 varies this assumption; Fig. 10 only tests convergence with r_crit. The paper's other approximations are individually tested (Secs. 5.2–5.6), but this dipole decomposition—the one most directly tied to the satellite features—is not. Until it is checked, the claimed broader agreement with FUSE rides on an untested bookkeeping choice.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a multi-basis method for simulation-based Stark-broadening line shape codes, in which the radiating system is evolved in a single-center atomic basis when the nearest-neighbor perturber is far, and in a two-center molecular basis during close collisions (Sec. 4.1). The time-evolution operator is transformed between bases using infinite-separation overlap integrals (Eqs. 16, 21, 26), and the molecular dipole is similarly projected onto atomic states after time evolution (Eq. 33). The method is implemented in the Xenomorph code and applied to hydrogen Lyman-alpha and Lyman-beta lines with quasi-H2+ satellites at white-dwarf photosphere conditions. The resulting profiles are compared with ULBT profiles (Figs. 4–6), and a grid of profiles is inserted into Tlusty to generate model spectra, which are compared qualitatively with a FUSE spectrum of the DA white dwarf Wolf 1346 (Fig. 15). The paper reports broader quasi-molecular features than ULBT and good qualitative agreement with the observed spectrum.","tokens_in":23944,"tokens_out":3897,"duration_ms":42216,"significance":"If the method holds up, it is a substantial advance: it provides the first simulation-based line shape implementation of quasi-molecular structure, removes several ULBT approximations (single-velocity, no screening, separate ion/electron broadening, no directional correlations), and reproduces ULBT profiles in a simplified limit (Fig. 6). The model-spectrum comparison to Wolf 1346 is suggestive and could help resolve known UV/optical discrepancies in DA white dwarfs. However, the central validation rests on at least one approximation that the authors explicitly concede is not valid in the regime where the quasi-molecular satellites form, and at least one additional ad hoc symmetry assumption in the appendices. These need to be tested or justified quantitatively before the observational agreement can be taken as evidence for the method.","major_comments":[{"comment":"The transformation of the time-dependent molecular dipole D_m(t) into the atomic basis uses infinite-separation overlap integrals at every internuclear separation R. The paper concedes in Sec. 4.3 that this approximation is 'no longer valid as R can be arbitrarily small.' This is precisely the regime that produces the quasi-H2+ satellites: Sec. 5.3 notes the 1400 Å Ly-alpha feature forms at R ~ 10 a0, and the 1060/1080 Å Ly-beta features are governed by similar close collisions. Because the line shape is obtained by Fourier-transforming D_fi(t) (Eq. 11), any systematic misassignment of molecular oscillator strength among atomic final states at small R will directly modify the width, strength, and position of the satellites that drive the agreement in Fig. 15. The paper provides no test of this approximation; Fig. 10 varies only r_crit, not the overlap prescription. A recomputation with R","section":"Sec. 4.3, Eq. (33)"},{"comment":"The inter-atomic contribution to the line shape is retained through the symmetry approximations F{U^†_pp D_p U_pr} ≈ F{U^†_rr D_r U_rp} and F{U^†_pr D_p U_pr} ≈ F{U^†_rp D_r U_rp}. These are asserted on the basis of H2+ symmetry and a 'time-averaged power spectrum' equivalence, but no test is shown. They are load-bearing: Sec. 5.8 and Fig. 11 demonstrate that the inter-atomic term is comparable to the single-site term in the line wings, exactly where the quasi-molecular satellites appear. Replacing these two terms with a symmetry assumption, without quantitative verification against a two-center simulation that tracks both atoms, re-introduces an uncontrolled approximation into a method whose stated goal is to remove ad hoc ULBT-style truncations.","section":"App. B, Eqs. (B11)-(B12)"}],"minor_comments":[{"comment":"The text says 'observational data shown later in Ch. 6' but the comparison is in the same section; 'Ch. 6' should be 'Sec. 6'.","section":"Sec. 6"},{"comment":"The legend uses 'dot-dash' for the Pelisoli+ (2015) profile, but the line style in the figure is not obvious in the printed version. Consider making the line styles more distinct and matching the caption wording.","section":"Fig. 4"},{"comment":"The r_crit convergence test in Fig. 10 is shown only for the Ly-beta red-wing features. Since the 1400 Å Ly-alpha satellite is also used in the spectral comparison, it would be useful to show r_crit convergence for that feature as well.","section":"Sec. 5.7"},{"comment":"The new line shapes are said to be 'available from the corresponding author upon request.' For reproducibility, a permanent archive (e.g., Zenodo or a journal repository) would be preferable, especially because the model comparison in Fig. 15 depends on the exact numerical profiles.","section":"Data Availability"},{"comment":"Author affiliation 7 contains the typo 'Tuscon' and should read 'Tucson'.","section":"Affiliation"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a genuinely novel and useful approach, and the acknowledgment of the Eq. (33) limitation is honest. However, the untested infinite-R overlap approximation sits exactly at the point where the paper claims improved agreement with observations, so I cannot support acceptance without a dedicated sensitivity test. The symmetry assumptions in App. B are also weaker than the main derivation and should be either substantiated or removed in a revised version. The observational comparison is only qualitative, as the authors state; I would not treat Fig. 15 as a strong validation until the underlying line-shape approximations are independently checked."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a genuinely new method—switching between an atomic and a molecular basis inside a simulation-based line-shape code—and the paper is honest about its own approximations. The main unresolved question is whether the infinite-separation overlap approximation used to project the molecular dipole onto the atomic basis throughout a close collision actually holds at the small R where the satellites form. The authors admit it does not, but they use it anyway for tractability and do not test it. That is the load-bearing soft spot.\n\nWhat is actually new: the multi-basis procedure itself (Sec. 4), the implementation in Xenomorph, and the demonstration that the resulting Lyα/Lyβ profiles reproduce the ULBT when the ULBT's approximations are artificially reintroduced (Fig. 6), while also removing some of those approximations (thermal velocity distribution, directional correlations, simultaneous ion+electron broadening). The comparison to the FUSE spectrum of Wolf 1346 is qualitative—one star, per-model normalization, no error bars—but the improvement in the Lyβ red wing is visible, and the authors present it as initial evidence rather than a finished fit. That is the right framing.\n\nThe soft spots are real but proportionate. The infinite-separation overlap approximation is used exactly at the R values where the quasi-molecular features form (R ≈ 10 a0), and the only basis-related test (Fig. 10) varies r_crit, not the overlap treatment. The arbitrary critical radius, the µ-ion symmetry assumption in App. B, and the low-density rescaling are also acknowledged but not stress-tested; the µ-ion simplification in particular rests on a symmetry argument that is plausible but not demonstrated. The code and line-shape data are not public, which limits independent verification. None of these are fatal, but together they mean the paper's central claim—broader quasi-H2+ features that agree with observations—rides on an approximation that is plausible but unverified.\n\nWho this is for: anyone working on white dwarf atmosphere modeling, line broadening theory, or H2+ quasi-molecular physics. The paper deserves a serious referee; the method is important and the authors are transparent. My recommendation: send it to peer review, but require the authors to test—or at least bound—the infinite-separation overlap approximation before publication. An R-dependent overlap calculation for a subset of trajectories, or a simplified one-dimensional test, would go a long way. The observational comparison should also be expanded beyond one star, even if only to a small sample, before strong claims about resolving UV/optical discrepancies are made.","headline":"Genuinely new multi-basis method for quasi-molecular line shapes, honestly written, but the untested infinite-separation overlap projection is a load-bearing approximation that needs checking before the claimed spectral agreement can be trusted.","tokens_in":24513,"tokens_out":2954,"would_cite":true,"duration_ms":28494,"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":"By switching between atomic and molecular bases in a line-shape simulation, the new method yields Lyman-series quasi-H2+ satellites that are broader than standard profiles and closer to observed white-dwarf ultraviolet spectra.","keywords":["quasi-molecules","H2+ satellites","Lyman series","Stark broadening","white dwarf atmospheres","simulation line shapes","multi-basis method","ultraviolet spectra"],"falsifier":"Compute the multi-basis profiles using R-dependent overlap integrals in the dipole transformation and check whether the 1400 angstrom Ly-alpha and the 1060/1080 angstrom Ly-beta satellites move or change width. If they do, the qualitative match to the observed white-dwarf spectrum is not decisive; a controlled laboratory measurement of hydrogen plasma at roughly 10,000 K and 10^17 cm^-3 would settle which profile family is correct.","tokens_in":23455,"feed_emoji":"🔭","tokens_out":5377,"duration_ms":58064,"temperature":0.7,"pith_summary":"The paper tries to establish that quasi-molecular structure—the transient H2+ states formed when a hydrogen atom and a proton collide without binding—can be built into simulation-based Stark broadening by switching between an atomic basis and a molecular basis during the closest part of each encounter. The authors implement this multi-basis procedure, compute hydrogen Lyman-alpha and Lyman-beta profiles with quasi-H2+ satellites, and insert them into a stellar atmosphere code. They find that the resulting satellites are broader and smoother than those from the standard semi-analytic unified broadening theory, and that model spectra built from them agree well with an observed ultraviolet spectrum of a hydrogen-atmosphere white dwarf. If right, this supplies a missing ingredient in the line shapes used to fit white-dwarf spectra, potentially easing known discrepancies between ultraviolet and optical estimates of temperature and mass.","feed_headline":"Quasi-molecule line shapes match white-dwarf UV spectra","feed_subtitle":"A line-shape simulation that switches between atomic and molecular bases broadens the hydrogen satellites that standard models get wrong.","key_machinery":"The load-bearing mechanism is a critical-radius basis switch: when the separation R between the radiating atom and the nearest perturbing ion drops below a chosen radius, the time-evolution operator U(t) and the dipole operator are transformed into a precomputed H2+ molecular basis; when R rises again, they are transformed back. The molecular Hamiltonian uses Born-Oppenheimer potential energy curves and dipole moments as functions of R, while distant perturbers are included through the same multipole electric-field interaction used in the atomic phase. Unitarity is retained by enlarging the atomic basis to include perturber states and using infinite-separation overlap integrals, and the off-","core_discovery":"The paper's central claim is that a single simulation line shape code can account for close-collision quasi-molecules by time-evolving the system in a molecular basis whenever the nearest-neighbor ion falls inside a critical radius, and in the usual atomic basis otherwise. To keep the two representations equivalent, the atomic basis is extended to include both radiator and perturber states and the transformation uses overlap integrals evaluated at infinite separation; this preserves unitarity and, after the collision, produces nonzero 'inter-atomic' blocks in the time-evolved dipole that encode charge exchange. The resulting hydrogen Lyman profiles have quasi-H2+ resonances that are broader","pith_inferences":["A systematic fit to many ultraviolet white-dwarf spectra would test whether the broader satellites actually resolve the known ultraviolet-optical discrepancies in effective temperature and mass; the paper only demonstrates qualitative agreement with one star and explicitly defers such fits.","The approximation the paper itself flags in Section 4.3—infinite-separation overlap integrals used at arbitrarily small R—means the very features being compared to observations could shift if R-dependent overlaps were used; recomputing with R-dependent overlaps is the natural next test.","Because the same dipole basis transformation is central, the method should transfer to Balmer lines, heavier perturbers, and two-electron quasi-molecules; the paper lists these as future work, and the mechanism has no evident Lyman-specific barrier.","A laboratory spectrum of a hydrogen plasma at white-dwarf photosphere conditions could discriminate between the broader multi-basis satellites and the standard semi-analytic satellites, since the paper's comparison to one star is qualitative."],"forward_implications":["Quasi-molecular resonances no longer have to be bolted onto Stark-broadened profiles afterward; they emerge from the same time-dependent simulation that handles ordinary Stark broadening.","Simultaneous ion and electron broadening replaces the usual separate treatment, so line cores are no longer distorted by convolving or adding independent ion- and electron-only profiles.","Inclusion of inter-atomic transitions is required: without them, the far wing is too weak relative to the core.","The broader, smoother satellites change the model flux between Ly-alpha and Ly-beta and improve the slope of the Ly-beta red wing compared with an observed ultraviolet white-dwarf spectrum.","The approach is the first simulation line shape implementation to include quasi-molecular structure, bringing Lyman-series profiles into the same simulation machinery used for ordinary Stark broadening."],"fun_headline_variants":["Quasi-molecule line shapes match white dwarf UV spectra","New line shapes fix white-dwarf spectral mystery","Atomic-molecular switch broadens hydrogen satellites","Simulation reproduces white dwarf hydrogen UV lines","Multiple-basis method improves stellar line fitting"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The assumption that the dipole-moment transformation between atomic and molecular bases can be done with overlap integrals evaluated at infinite separation, even while the collision distance R is arbitrarily small, is the load-bearing premise; the paper itself states this is no longer valid as R becomes small.","fun_headline_variants_meta":{"raw":{"variants":["Quasi-molecule line shapes match white dwarf UV spectra","New line shapes fix white-dwarf spectral mystery","Atomic-molecular switch broadens hydrogen satellites","Simulation reproduces white dwarf hydrogen UV lines","Multiple-basis method improves stellar line fitting"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000104,"raw_usage":{"total_tokens":867,"prompt_tokens":742,"completion_tokens":125,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":486,"completion_tokens_details":{"reasoning_tokens":54}},"tokens_in":486,"tokens_out":125,"duration_ms":2419,"temperature":1.0,"reasoning_tokens":54,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T17:09:54.004766+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the multi-basis profiles using R-dependent overlap integrals in the dipole transformation and check whether the 1400 angstrom Ly-alpha and the 1060/1080 angstrom Ly-beta satellites move or change width. If they do, the qualitative match to the observed white-dwarf spectrum is not decisive; a controlled laboratory measurement of hydrogen plasma at roughly 10,000 K and 10^17 cm^-3 would settle which profile family is correct.","supporting_citations":[],"review_version":1}