{"id":"3ca59bff-1385-45d0-bf96-27eea62b2acf","arxiv_id":"2502.04705","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"SPH models of an eccentric Be/X-ray binary show the disc never reaches steady state and neutron star accretion is strongest in near-coplanar prograde geometries, decreasing by orders of magnitude for misaligned retrograde orbits.","lead":"SPH simulations of a Be star disc in a highly eccentric, short-period binary (A0538-66-like) sweep the neutron star's orbital tilt from 0 to 180 degrees and track the disc's mass, shape, and spin over 80 orbits. The result: all disc quantities pulse with the orbital period, and neutron star accretion is strongest for near-coplanar prograde orbits and much weaker for retrograde ones.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Single-orbit mean accretion rates in Table 2 lack error bars; with accretion concentrated in periastron bursts, the factor-1.7 gap between the 75° and 180° models could be closed by sampling noise.","rationale":"The reader's named weakest assumption is the Shakura-Sunyaev scale-height relation near the secondary (Section 4.3). I agree that is a real, acknowledged limitation, but it is not the most load-bearing for the headline claim: the authors' alpha_SPH = 5 control changes accretion rates by <30%, while the prograde/retrograde separation in Table 2 spans factors of 2–100, so the hierarchy would likely survive a corrected viscosity. The more dangerous assumption is that a single-orbit mean (or an unstated number of orbits) is a reliable estimator of the steady-state accretion rate. The accretion time series is strongly phase-locked, with essentially all mass arriving in bursts near periastron (Figs 8 and 15). For low-rate models, each orbit contributes tens of particles (or one burst), so the mean over one orbit is dominated by sampling noise. The 75° vs 180° gap is only a factor of 1.7, and a 30% systematic shift of the kind found in Section 4.3 plus burst noise could plausibly erode it. Without error bars or a statement of how many orbits were averaged, the central 'largest for <90°, smaller for >90°' claim is not yet quantitatively secure. This matches the reader's Conditional verdict, so I keep the verdict unchanged; I only sharpen the specific concern to the statistical robustness of Table 2 rather than the viscosity prescription.","tokens_in":21923,"tokens_out":11272,"duration_ms":99900,"concrete_test":"For the 75°, 105°, 120°, 150°, and 180° models, continue each simulation for at least 20 orbits after the nominal quasi-steady state, record the accreted mass in every orbital period, and report the per-orbit mean and standard deviation. Assess whether the 75° (prograde) and 180° (retrograde) distributions are separated by more than the combined 1σ uncertainties, and compare the 20-orbit mean with the Table 2 entry to test sensitivity to the averaging window. This directly settles whether the prograde > retrograde grouping is a robust result or an artifact of single-orbit sampling.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline accretion hierarchy in Table 2 rests on mean rates with no quoted uncertainty, despite strong phase-locked bursting (Figs 7–8, 14–15). Because accretion is confined to brief periastron encounters, the effective sample size is the number of orbits averaged, not the number of particles. For the retrograde models, per-orbit particle counts are small: with particle mass 4.6e-15 M_sun and P = 16.6409 d (0.0456 yr), the 120° rate of 2.2e-12 M_sun/yr corresponds to only ~22 particles per orbit, and the mass may arrive in one dominant burst. If Table 2's 'averaged over an orbital period' means a single orbit, the quoted rates are effectively single-burst draws. The smallest prograde-to-retrograde gap in Table 2 is between 75° (2.7e-11 M_sun/yr) and 180° (1.6e-11 M_sun/yr), only a factor of 1.7. The authors' own viscosity sensitivity test (Section 4.3) changes rates by up to ~30%, so a systematic shift of that size combined with burst-to-burst variability could erode or reverse this gap. The paper does not state how many orbits were averaged, the standard deviation, or the averaging window, so the central 'largest for <90°, smaller for >90°' ordering is not yet statistically established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents 3D SPH simulations of a Be star decretion disc in a highly eccentric, short-period Be/X-ray binary with parameters based on A0538-66, varying the neutron star orbital misalignment angle from 0 to 180 degrees. The authors track disc mass, angular momentum, eccentricity, inclination, and neutron star accretion rate, and report that all disc quantities vary with orbital phase, with prograde misalignments (<90 degrees) generally yielding higher accretion rates than retrograde ones (>90 degrees). The paper also tests an alternative viscosity prescription and reports that disc-scale evolution is robust, while accretion-rate differences can reach roughly 30% in some misaligned cases.","tokens_in":22078,"tokens_out":3858,"duration_ms":41956,"significance":"If the central accretion-rate hierarchy is robust, this is the first systematic SPH survey of misalignment effects in eccentric Be/X-ray binaries, and it would provide useful input for interpreting periodic Type I X-ray outbursts and for future observational predictions. The paper has clear strengths: it uses a broad and well-motivated parameter grid, all parameters except misalignment are fixed from prior literature, the headline quantities are emergent simulation outputs rather than fitted, and the phase-locked variability is convincingly illustrated in the figures. The main risk is quantitative rather than conceptual: the Table 2 accretion-rate hierarchy lacks uncertainty estimates, and the paper's own viscosity test changes rates by amounts comparable to the smallest margins in that hierarchy.","major_comments":[{"comment":"The central claim that accretion rates are largest for misalignment angles less than 90 degrees and smaller for angles greater than 90 degrees rests on single-orbit mean values in Table 2 with no error bars, no number of averaged orbits, and no run-to-run variance. The smallest prograde-to-retrograde margin is only a factor of 1.7 (75 degrees at 2.7e-11 M_sun/yr versus 180 degrees at 1.6e-11 M_sun/yr). Accretion is strongly burst-like (Figs. 7, 8, 14, 15), and at the quoted particle mass of 4.6e-15 M_sun and period of 0.0456 yr, the 120-degree rate corresponds to roughly 20 particles per orbit, so Poisson and burst-to-burst fluctuations could plausibly be comparable to or larger than the smallest reported gaps. The authors should report the averaging window, the scatter across orbits, and ideally multiple realizations or a phase-binned standard error, before the prograde/retrograde ordering can be considered established.","section":"Table 2; Sections 3.1.2 and 3.2.2"},{"comment":"The authors themselves note that Eq. (6) overestimates the scale height for particles bound to the neutron star, forcing the SPH artificial viscosity to compensate and affecting accretion timescales near the secondary. The corrective test with alpha_SPH = 5 changes accretion rates by up to roughly 30%, with the largest effects in the 45- and 105-degree models. Since the Table 2 hierarchy contains neighboring-model margins of only about 1.7-2.4, and since the viscosity sensitivity is largest on both sides of the 90-degree divide, this systematic uncertainty could alter or even invert specific orderings in the headline result. I ask the authors to present the viscosity-test accretion rates quantitatively (not only as time-series examples) and to state explicitly whether the '<90 vs >90' hierarchy is preserved under the alpha_SPH prescription in all cases.","section":"Section 4.3, Eq. (6)"}],"minor_comments":[{"comment":"The statement 'No new data were generated or analysed in support of this research' is inconsistent with the simulation outputs underlying Figs. 2-19; please clarify whether simulation outputs are available on request or through a repository.","section":"Data Availability"},{"comment":"The panels labeled 'i w.r.t. Secondary' mix angles with respect to the binary orbital plane and the secondary star's spin; please standardize the terminology and axis labels so the reader can distinguish these two reference frames.","section":"Figures 2 and 9"},{"comment":"The phrase 'accretion rates are strongly correlated with orbital phase' is qualitative; consider reporting a quantitative metric, such as the fraction of accreted mass within a given phase window or a phase-binned mean and standard deviation.","section":"Sections 3.1.2 and 3.2.2"},{"comment":"For the accretion-rate rows, please state explicitly how many orbital periods are included in the average and whether the quasi-steady-state interval is the same for all models; this information is needed to interpret the single quoted values.","section":"Table 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and the modelling effort is valuable. My recommendation is driven by the need to make the headline accretion-rate claim statistically and systematically robust; the disc-evolution results themselves appear solid and would survive a minor revision, but the quantitative hierarchy in Table 2 needs uncertainty quantification and a more quantitative reporting of the viscosity test. I would be willing to consider a revised version positively if these points are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the full 0-to-180 degree misalignment sweep in a highly eccentric (e=0.72), short-period Be/X-ray binary. Prior work covered circular orbits with misalignment up to 60 degrees, or coplanar eccentric cases, or a single misaligned eccentric case. This parameter study fills that gap systematically, and the phase-resolved disc diagnostics are a real step forward. The qualitative results – all disc quantities vary with orbital phase, and accretion is systematically larger for misalignment under 90 degrees – are internally consistent and survive the authors' own viscosity prescription switch, which is a good sign. The discussion of the two competing effects (plane overlap vs. relative velocity) is clear and physically sensible. I also appreciate that the authors flag the scale-height issue near the secondary star in Section 4.3 rather than hiding it.\n\nThe soft spots are mostly about the quantitative claims. Table 2 gives single-orbit mean accretion rates with no scatter, run-to-run variance, or averaging window, despite the accretion being strongly bursty and phase-locked. The stress-test note is on target: for the retrograde models, the per-orbit particle counts are tiny (tens of particles), so the quoted rates are effectively single-burst draws. The smallest gap in the prograde/retrograde ordering, between 75° and 180°, is only a factor of 1.7, and the authors' own viscosity test changes rates by up to ~30%. That means the specific ordering at the boundary is not yet statistically established. The broad trend (prograde far outpacing retrograde) is robust, but the fine structure needs multi-orbit statistics or at least standard deviations.\n\nThe viscosity caveat is acknowledged but not resolved: the overestimated scale height near the secondary affects exactly the region where accretion is measured, with the largest effects in the 45° and 105° models, which sit at the prograde/retrograde boundary. That tempers any quantitative reading of the hierarchy, though the sign of the effect is not obviously fatal. Also, the Data Availability statement says no data were generated, which is a miss for a numerical study; releasing accretion-rate time series and state dumps would let others check the burst statistics. Finally, the generalization from one A0538-66-like fiducial to all short-period eccentric Be/X-ray binaries is a stretch; this is a model family, not a survey.\n\nBottom line: this deserves a serious referee. The novelty is real, the physical reasoning is sound, and the fixable issues – error bars, data release, and tempering the generalization – are exactly what an effective review should ask for. I'd want those addressed before citing the numbers, but I'd probably cite the geometry sweep regardless.","headline":"A genuinely new geometry sweep for eccentric Be/X-ray binaries, with a plausible qualitative trend, but the headline accretion numbers need error bars before the prograde/retrograde ordering is taken as quantitative.","tokens_in":746,"tokens_out":1446,"would_cite":true,"duration_ms":30465,"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":"In simulations of the eccentric Be/X-ray binary A0538-66, neutron star accretion is strongest for prograde misalignments below 90° and weakest for retrograde misalignments above 90°, because tilt sets both particle encounters and…","keywords":["Be/X-ray binaries","decretion discs","smoothed particle hydrodynamics","misalignment angle","neutron star accretion","eccentric binary","A0538-66","X-ray outbursts"],"falsifier":"Rerun the 45° and 105° models with a viscosity that computes the local scale height from the actual particle distribution rather than from the relation $H(r) = c_s/v_{\\mathrm{crit}} (r/R_\\star)^{1.5}$: the central hierarchy survives only if the prograde 45° model still accretes faster than the retrograde 105° model. An observational counter-check would be a sample of eccentric Be/X-ray binaries in which bright periastron outbursts show no systematic preference for prograde-aligned systems.","tokens_in":21568,"feed_emoji":"🌠","tokens_out":10698,"duration_ms":94709,"temperature":0.7,"pith_summary":"This paper sets out to show that in short-period, highly eccentric Be/X-ray binaries the orientation of the neutron star's orbit relative to the Be star's disc controls how much matter the neutron star accretes. Using ten 3D smoothed particle hydrodynamics simulations of the A0538-66 system that sweep the misalignment angle from coplanar prograde (0°) to coplanar retrograde (180°), the authors find a systematic hierarchy: accretion is highest for prograde misalignments below 90° and lowest for retrograde misalignments above 90°, from about $3\\times10^{-10}\\,M_\\odot\\,\\mathrm{yr}^{-1}$ in the coplanar prograde case down to a few $\\times10^{-12}\\,M_\\odot\\,\\mathrm{yr}^{-1}$ in the retrograde family. The driving mechanism is geometric: the tilt of the orbit sets both how many disc particles the neutron star meets and how long it interacts with each one before relative motion separates them. The simulations also establish that the high eccentricity forces every disc property—mass, angular momentum, eccentricity, inclination—to oscillate with orbital phase, so accretion is phase-dependent even in a quasi-steady disc. A sympathetic reader would care because this connects binary geometry directly to the timing and strength of X-ray outbursts and to observable disc variability.","feed_headline":"Prograde tilts feed neutron stars; retrograde tilts starve them","feed_subtitle":"Simulations of A0538-66 show geometry, not chance, sets which misalignments accrete fastest","key_machinery":"The load-bearing machinery is a suite of ten three-dimensional smoothed particle hydrodynamics (SPH) simulations of a Be-star decretion disc interacting with a neutron star, all initialized with the parameters of A0538-66 ($e=0.72$, $P\\approx16.6$ d) and differing only in the misalignment angle of the neutron star's orbital plane, stepped from 0° (coplanar prograde) to 180° (coplanar retrograde). The argument is carried by two competing geometric effects: how much of the disc plane the neutron star's orbit overlaps, which sets the number of disc particles it encounters, and the relative velocity between the neutron star and those particles, which sets how long each encounter lasts. The disc transport is set by the Shakura–Sunyaev viscosity prescription with $\\alpha_{\\mathrm{SS}} = 0.5$, and the mapping to the SPH artificial viscosity uses the isothermal scale-height relation $H(r) = c_s/v_{\\mathrm{crit}} (r/R_\\star)^{1.5}$; particle splitting increases resolution near the neutron star so that accretion rates into the Eggleton Roche-lobe sink radius can be measured.","core_discovery":"The paper's central claim is that in a highly eccentric Be/X-ray binary like A0538-66, the misalignment angle between the neutron star's orbital plane and the Be star's equator sets a systematic ordering of accretion efficiency and disc response. For prograde misalignments (0°–75°) the neutron star's velocity roughly aligns with the disc particles, so overlapping the disc for longer and with smaller relative velocities yields the highest accretion rates; rates fall steadily as the angle grows. For retrograde misalignments (105°–180°), the same two effects compete rather than cooperate: the coplanar retrograde case maximizes particle encounters but with strongly antiparallel velocities, and misaligned retrograde orbits reduce both encounter number and interaction time, so accretion rates are systematically lower than in the prograde family. The paper also claims that in every model the disc's mass, angular momentum, eccentricity, and inclinations oscillate with orbital phase, with the largest disruption at periastron and a recovery afterward, and that the orbital phase of peak accretion shifts from periastron for near-coplanar cases to multiple post-periastron peaks for highly misaligned prograde cases.","pith_inferences":["If the same geometric competition operates in other short-period eccentric Be/X-ray binaries, the observed spread of Type I outburst fluences could be inverted to infer typical misalignment angles: the brightest accretors would be the most prograde-aligned systems.","A natural extension is to vary eccentricity: the paper's interaction-time argument implies the prograde/retrograde accretion gap should narrow at lower eccentricity, where periastron relative velocities are less extreme.","The angular momentum of captured retrograde material is opposite to the neutron star's spin direction, so retrograde accretors might show systematically different pulse-period changes than prograde accretors; this is not addressed in the paper.","Multi-wavelength monitoring of a single system could test the phase-locking corollary: if the simulated disc tilts are real, polarization angle and Balmer-line equivalent width should oscillate on the orbital period with amplitude growing as the misalignment approaches 90°."],"forward_implications":["Time-averaged neutron star accretion in eccentric Be/X-ray binaries should follow a prograde-favored ordering: for the same stellar and orbital parameters, misalignments below 90° out-accrete misalignments above 90°.","Every disc observable in these systems should vary with orbital phase, with mass and angular momentum dipping at periastron and rebuilding during the rest of the orbit, so phase-resolved observations should see periodic dips and recoveries.","For near-coplanar orbits the peak accretion occurs at or just after periastron; for highly misaligned prograde orbits the peak splits into multiple events as the neutron star re-crosses the disc plane or runs into its own induced spiral arms.","Retrograde orbits should produce weaker spiral arms and less disc disruption than prograde orbits, with coplanar retrograde accretion more efficient than any misaligned retrograde case because of the larger number of particles encountered.","The 30°–45° misaligned models develop accretion streams from the inner disc to the neutron star just after periastron, providing a geometric channel for the periodic Type I X-ray outbursts seen in these systems."],"supporting_citations":[{"why":"Supplies the SPH code adaptation for Be-star decretion discs and establishes the two-armed spiral and disc-truncation behaviour this work extends to misaligned eccentric orbits.","marker":"Okazaki et al. (2002)"},{"why":"Defines the $\\alpha_{\\mathrm{SS}}$ viscosity prescription (Eq. 2) that sets disc shear viscosity in all models.","marker":"Shakura & Sunyaev (1973)"},{"why":"Gives the isothermal scale-height relation (Eq. 6) used to convert $\\alpha_{\\mathrm{SS}}$ to the SPH artificial viscosity $\\alpha_{\\mathrm{SPH}}$.","marker":"Carciofi & Bjorkman (2006)"},{"why":"Provides the A0538-66 system parameters ($e=0.72$, $P\\approx16.6$ d, primary mass/radius/effective temperature) used to set the models.","marker":"Rajoelimanana et al. (2017)"},{"why":"Approximates the Roche lobe radius (Eq. 1) from which the neutron star accretion radius is defined.","marker":"Eggleton (1983)"},{"why":"Foundational SPH code that the simulations are built on.","marker":"Benz et al. (1990)"},{"why":"Refined the SPH code used here, including the sink-particle treatment of the two stars.","marker":"Bate et al. (1995)"},{"why":"Earlier coplanar prograde/retrograde Be/X-ray binary simulations that this work extends to a full misalignment grid.","marker":"Panoglou et al. (2016)"},{"why":"Supplies the misaligned-disc methodology and the equations for specific angular momentum and inclination that this paper reuses.","marker":"Suffak et al. (2022)"},{"why":"Circular misaligned binary simulations that provide the injection-radius and disc-truncation baseline for tilted discs.","marker":"Cyr et al. (2017)"}],"fun_headline_variants":["Prograde tilts feed neutron stars more than retrograde","A0538-66: prograde misalignment maximizes neutron star accretion","Eccentric Be/X-ray: prograde orbits outfeed retrograde","Geometry dictates feeding: prograde beats retrograde in A0538-66","Neutron star accretion ordered by misalignment angle"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that the viscous transport law that describes the Be disc also applies to gas once it becomes bound to the neutron star; the paper itself notes that its scale-height formula overestimates the thickness there, and using a fixed artificial-viscosity alternative changes the measured accretion rates by up to roughly 30 percent, with the largest shifts at 45° and 105°.","fun_headline_variants_meta":{"raw":{"variants":["Prograde tilts feed neutron stars more than retrograde","A0538-66: prograde misalignment maximizes neutron star accretion","Eccentric Be/X-ray: prograde orbits outfeed retrograde","Geometry dictates feeding: prograde beats retrograde in A0538-66","Neutron star accretion ordered by misalignment angle"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000311,"raw_usage":{"total_tokens":1835,"prompt_tokens":1075,"completion_tokens":760,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":691,"completion_tokens_details":{"reasoning_tokens":670}},"tokens_in":691,"tokens_out":760,"duration_ms":7630,"temperature":1.0,"reasoning_tokens":670,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T21:49:26.571580+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Rerun the 45° and 105° models with a viscosity that computes the local scale height from the actual particle distribution rather than from the relation $H(r) = c_s/v_{\\mathrm{crit}} (r/R_\\star)^{1.5}$: the central hierarchy survives only if the prograde 45° model still accretes faster than the retrograde 105° model. An observational counter-check would be a sample of eccentric Be/X-ray binaries in which bright periastron outbursts show no systematic preference for prograde-aligned systems.","supporting_citations":[{"cited_title":"T., Bate M","cited_arxiv_id":null,"evidence_quote":"Supplies the SPH code adaptation for Be-star decretion discs and establishes the two-armed spiral and disc-truncation behaviour this work extends to misaligned eccentric orbits."},{"cited_title":"E., Carciofi A","cited_arxiv_id":null,"evidence_quote":"Supplies the misaligned-disc methodology and the equations for specific angular momentum and inclination that this paper reuses."},{"cited_title":"H., Jones C","cited_arxiv_id":null,"evidence_quote":"Circular misaligned binary simulations that provide the injection-radius and disc-truncation baseline for tilted discs."}],"review_version":1}