{"id":"2849003a-db59-4871-ae76-345fbf2de27e","arxiv_id":"2412.04299","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Low-mass companions and higher viscosity suppress Kozai-Lidov oscillations and disc tearing in Be star disc models, while radiative transfer predicts a first-ever triple-peaked H-alpha line from the oscillating disc.","lead":"This paper simulates Be star discs in misaligned binaries with new mass ratios and viscosities, finding that light companions or high viscosity suppress disc tearing and Kozai-Lidov oscillations, and it computes predicted spectral, polarimetric, and interferometric signatures.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq (2)'s alpha_ss-to-alpha_sph conversion is unquantified in warped/tearing regions; if effective alpha_ss is inflated, the claimed q/alpha_ss boundaries (including the q=0.1 null) shift.","rationale":"I read the paper in good faith. The central claim is a mapping between binary mass ratio/disc viscosity and the occurrence of KL oscillations and disc tearing in Be star discs, plus novel radiative-transfer predictions. The strongest load-bearing assumption is the fidelity of Eq (2) in converting Shakura-Sunyaev alpha_ss to SPH artificial viscosity. The conversion is designed to cancel h, but it uses the theoretical scale height H = c_s/Omega; in warped, tilted, tearing discs this prescription is not guaranteed. The manuscript itself flags this in Section 4.2, noting 'areas in the simulation where the disc does not follow this scale height prescription may have inflated artificial viscosity values.' No quantitative check is provided for the primary runs. Given that the paper's own alpha_ss=0.5 and 1.0 runs, and its constant-alpha_sph comparison, show that modest viscosity increases suppress tearing and KL, an unquantified inflation of the effective alpha_ss in the nominal 0.1 runs could shift the claimed boundaries. The q=0.1 null result, based on only two runs, is particularly sensitive because the same inflation could suppress the phenomena for viscosity reasons rather than mass-ratio reasons. The proposed resolution/convergence test would directly measure the realized alpha_ss. Other concerns (sparse parameter sampling, the 'first triple peak' novelty claim, circular-orbit assumption) are either explicitly qualified in the paper ('for the parameters examined') or do not affect the core dynamical mapping. Therefore the load-bearing concern aligns with the reader's weakest assumption, and the appropriate verdict remains CONDITIONAL.","tokens_in":19366,"tokens_out":10237,"duration_ms":91917,"concrete_test":"Re-run the q=0.5, beta=40, alpha_ss=0.1 tearing model (Model 2) with 4x the particle injection rate (2e6 particles per orbital period), and compute the radial profile of the effective alpha_ss = (1/10) alpha_sph h/H from the actual SPH particle data in the tearing region. If the disc still tears and the measured alpha_ss is within ~10% of 0.1 throughout the inner disc, the mapping is validated; if the effective alpha_ss is inflated or tearing disappears, the claimed q/alpha_ss boundaries are not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central mapping—where misaligned Be discs tear or undergo KL as a function of q and alpha_ss—rests on Eq (2), which converts a chosen Shakura-Sunyaev alpha_ss into the SPH artificial viscosity parameter alpha_sph using the theoretical isothermal scale height H = c_s/Omega. Equation (2) cancels the smoothing length h only if the disc's vertical structure actually follows H = c_s/Omega. Section 4.2 states the caveat explicitly: 'areas in the simulation where the disc does not follow this scale height prescription may have inflated artificial viscosity values (for example, when the disc becomes warped and tilted...).' The paper's main results occur precisely in warped, tilted, tearing discs, yet no measurement of the realized effective alpha_ss is presented for the primary (constant-alpha_ss) runs. Figure 6 only diagnoses the constant-alpha_sph comparison runs, for which H/h < 1 inflates alpha_ss. If the effective viscosity in the nominal alpha_ss=0.1 runs is inflated above 0.1 (e.g., in the outer disc, where H_theory grows as r^(3/2)), then the boundary for tearing/KL shifts to lower nominal alpha_ss. The q=0.1 null result is especially vulnerable: only two q=0.1 runs exist, and an inflated viscosity in those runs could suppress tearing/KL for reasons unrelated to mass ratio. Since the paper's own high-alpha_ss runs (0.5, 1.0) show that modest viscosity increases damp these phenomena, the quantitative boundaries are only as secure as the fidelity of Eq (2).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the authors' earlier SPH study of misaligned Be star binary discs (Paper I) by running eight new simulations that vary the binary mass ratio (q = 0.1, 0.5, 1.0) and the Shakura–Sunyaev viscosity parameter (alpha_ss = 0.1, 0.5, 1.0), and by post-processing selected runs with the radiative transfer code HDUST. The central dynamical claims are that q = 0.5 reproduces the previously seen disc-tearing (for 40 deg misalignment) and Kozai–Lidov (KL) oscillations (for 60 deg) with somewhat longer timescales, while q = 0.1 shows neither phenomenon in the parameters examined; that increasing alpha_ss to 0.5 or 1.0 damps these oscillations; and that a constant-SPH-viscosity prescription yields effectively higher alpha_ss values, explaining why tearing is suppressed in those comparison runs. On the observational side, the paper predicts that KL oscillations produce periodic changes in H-alpha equivalent width, V/R ratio, peak separation, polarization degree and position angle, including a triple-peaked H-alpha profile that the authors state is the first produced by radiative transfer modelling of a Be star disc. It also presents interferometric visibility and differential-phase signatures of gaps caused by KL oscillations and disc tearing, and suggests Pleione as a candidate for detecting such signatures.","tokens_in":19798,"tokens_out":4823,"duration_ms":51362,"significance":"If the dynamical claims hold, the paper offers a useful map of how misaligned Be discs transition between tearing, KL oscillation, and quiescent precession as functions of mass ratio and viscosity, which is directly relevant to Be/X-ray binaries and the growing population of Be stars with faint companions. The strengths of the work are that it uses a well-established SPH code, cross-checks its results against independent analytic scalings (Doğan et al. tearing radius, Eq. 6; Larwood et al. precession, Eq. 7; Lidov conservation, Eq. 8; Kiseleva et al. KL timescale, Eqs 9–10), and forward-models observables with HDUST rather than tuning the models to match data. The radiative-transfer and interferometric predictions, especially the triple-peaked H-alpha profile and the visibility-hump signature of a torn disc, are novel, falsifiable predictions that should guide future observations of systems like Pleione. The main reservations are that the viscosity conversion underlying the claimed q/alpha_ss boundaries is not quantitatively validated in the warped and tearing regions, and that the q = 0.1 null result rests on only one simulation per misalignment angle.","major_comments":[{"comment":"The mapping between the chosen alpha_ss and the SPH artificial viscosity parameter alpha_sph in Eq. (2) is load-bearing for the paper's central conclusion that q = 0.1 does not produce KL oscillations or disc tearing, and that alpha_ss = 0.5 or 1.0 damps them. The manuscript itself notes in §4.2 that 'areas in the simulation where the disc does not follow this scale height prescription may have inflated artificial viscosity values,' and Fig. 6 demonstrates H/h < 1 for a constant-alpha_sph run, inflating the realized alpha_ss. However, no analogous diagnostic is provided for the constant-alpha_ss runs on which the main scientific claims rest, even though those runs are precisely the ones that become warped, tilted, and torn. As a result, the effective viscosity in the main runs may be substantially larger than the nominal alpha_ss = 0.1, which could suppress tearing/KL in the q = 0.1 runs for reasons unrelated to mass ratio. I request that the authors report the radial and temporal distribution of the realized effective alpha_ss (e.g., from alpha_sph * h/H) for the constant-alpha_ss models, and, ideally, test whether the q = 0.1 null persists when either the resolution is increased (so that H/h approaches unity) or the nominal alpha_ss is lowered. Without such a test, the q and alpha_ss boundaries are not yet quantitatively secure.","section":"§4.2, Eq. (2), Fig. 6"},{"comment":"The null result for q = 0.1 is based on exactly one simulation for each of the two misalignment angles (40 deg and 60 deg), with no variation of random seed or independent realisation. Given the stochastic nature of SPH particle injection and the sensitivity of tearing and KL oscillations to the detailed disc structure at the end of the growth phase, a single realisation per parameter set is insufficient to establish a null result. The authors should either present multiple realisations with different noise seeds, or demonstrate robustness by varying some other parameter (e.g., lower alpha_ss, different mass-injection profile) while keeping q = 0.1. This is important because the paper's abstract and conclusions state that q = 0.1 does not cause KL oscillations or disc tearing, which is a strong claim for the current evidence.","section":"§3, Table 2, Figs 2–3"},{"comment":"The statement that applying the factor 2.7 from Martin & Franchini (2019) to the predictions of Eq. (10) 'results in the exact initial KL periods found in our simulations' is an overstatement. That factor is derived for a test particle with an initial eccentricity of 0.2, whereas the disc here is extended, has a changing surface-density profile, and is dissipating; moreover, the measured periods are read from the first eccentricity peak, which involves some judgement given the oscillations. The agreement is encouraging and worth reporting, but it should be presented as an approximate consistency check (within the factor-of-2.7 approximation), not as an 'exact' reproduction. This is a minor point for the overall conclusions, but it appears in the main body as a quantitative validation of the KL interpretation.","section":"§3, Eq. (10), Figs 2–3"}],"minor_comments":[{"comment":"The sentence 'we find we can recover the same qualitative disc evolution for a range of constant artificial viscosity as when alpha_ss is held constant' is immediately followed by the statement that disc tearing does not occur in the constant-alpha_sph = 1 models where it previously did. These two sentences are in tension; please rephrase to clarify that the qualitative agreement holds for some diagnostics (disc mass, inclination, precession) but not for tearing/KL, which depend sensitively on viscosity.","section":"§4.2, Fig. 6"},{"comment":"The interferometric predictions assume a distance of 100 pc and a specific orientation, but the text does not state whether the uncertainties in the HDUST grid resolution (50x50x50) or in the disc density have been checked for their effect on the visibility curves. A brief comment on numerical convergence of the image/visibility computation would strengthen the quantitative claims about the size of the visibility humps.","section":"§6, Figs 12–13"},{"comment":"The statement 'No new data was generated' is confusing given that eight new SPH simulations and their HDUST observables were produced. I suggest rewording it to 'No observational data were generated; the simulation outputs are available from the authors upon request,' which would be consistent with the models actually being new.","section":"§8, Data Availability"},{"comment":"The sentence 'Recently, simulations of Be stars in misaligned binary systems have revealed that misalignment between the disc and binary orbit can cause the disc to undergo Kozai-Lidov oscillations or disc-tearing' could benefit from a citation to the early papers that established these phenomena in Be-disc context (e.g., Martin et al. 2014; Suffak et al. 2022), rather than only to Paper I and the later sections.","section":"§1, Introduction"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid extension of the authors' prior work and will likely be of interest to the MNRAS readership, especially the radiative-transfer and interferometric predictions. My main editorial concern is the novelty claim of the 'first triple-peaked H-alpha profile' from a Be star disc model; it appears to be true within the cited literature, but I would ask the editor to ensure a broad enough search was carried out, since triple-peaked profiles have been discussed in several contexts. The more substantive issue is the viscosity-conversion problem, which is openly acknowledged in §4.2 but not resolved for the main runs; the authors should be pushed to quantify the effective alpha_ss in those runs rather than relying on the nominal values. I also note that the data-availability statement should be corrected, as it currently says no new data were generated despite the paper presenting new simulations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper does what a good follow-up should: it fills in the q and alpha_ss grid around the authors' own Paper I simulations, and it produces genuinely new observable predictions. The q=0.5 runs behaving like the equal-mass case but with delayed timescales, and the q=0.1 runs showing neither tearing nor KL, are clean and consistent with the Doğan tearing-radius estimate and the KL timescale scalings. The high-viscosity runs damping the oscillations are also expected and well-explained. The radiative-transfer work is the real payoff: the triple-peaked H-alpha profile, the polarization position-angle swings, and the interferometric visibility/phase signatures for gaps are new and will be useful for interpreting real Be binaries. I believe the claim of being the first RT model to produce a triple-peaked Be line, but I would want a referee to check that literature point carefully.\n\nThe soft spots are real but not fatal. The stress-test worry about Eq. (2) is legitimate: the alpha_ss-to-alpha_sph conversion assumes H = c_s/Omega, which is exactly what breaks down in warped, torn discs. The paper itself flags this in Section 4.2, and the constant-alpha_sph comparison runs show the effect—those runs damp tearing precisely because H/h < 1 inflates the effective alpha_ss. So the quantitative q/alpha_ss boundaries are less secure than the qualitative story. The q=0.1 null result rests on two simulations, one per inclination, which is thin for a load-bearing conclusion. Also, the circular-orbit assumption is never stated explicitly; this limits applicability to eccentric Be/X-ray binaries, though it may be acceptable for a first mapping.\n\nI'd send this to peer review. The core hydrodynamics is checked against independent analytic scalings, the observable predictions are forward-modeled, and the limitations are acknowledged. A good referee should ask for a direct measurement or estimate of the realized effective alpha_ss in the constant-alpha_ss runs, and ideally a few more q=0.1 realizations, but these are addressable without invalidating the main results. Worth citing for the triple-peaked line and the interferometric gap predictions.","headline":"A useful parameter-space extension with the first radiative-transfer observables for KL oscillations, but the viscosity mapping uncertainty and thin q=0.1 sampling mean the claimed boundaries should be treated as provisional.","tokens_in":20301,"tokens_out":1562,"would_cite":true,"duration_ms":17622,"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":"The paper identifies the binary mass ratio and disc viscosity as the controls for whether misaligned Be star discs undergo Kozai-Lidov oscillations and disc tearing, and shows these states produce distinct observable signatures including…","keywords":["Be stars","circumstellar discs","Kozai-Lidov oscillations","disc tearing","smoothed particle hydrodynamics","radiative transfer","triple-peaked H-alpha","stellar interferometry"],"falsifier":"Run the constant-$\\alpha_{\\rm sph}$ equal-mass, 40-degree misalignment model with particle resolution high enough that $H/h \\approx 1$ throughout the disc; if disc tearing then appears, the suppression of tearing in constant-$\\alpha_{\\rm sph}$ models is an artifact of the viscosity conversion rather than a physical effect of the higher viscosity.","tokens_in":19155,"feed_emoji":"🌀","tokens_out":11629,"duration_ms":98350,"temperature":0.7,"pith_summary":"This paper extends three-dimensional smoothed-particle-hydrodynamics simulations of Be star discs in misaligned binary systems to ask what controls two dramatic behaviours: Kozai-Lidov (KL) oscillations, where the disc's eccentricity and inclination repeatedly trade off while it dissipates, and disc tearing, where the disc breaks into separate precessing rings. The paper finds that lowering the binary mass ratio from 1 to 0.5 preserves both phenomena with a slightly longer timescale, while a mass ratio of 0.1 suppresses them for the parameters examined, and that raising the Shakura-Sunyaev viscosity parameter from 0.1 to 0.5 or 1.0 damps them as well. Using the HDUST radiative-transfer code, the paper produces the first triple-peaked H$\\alpha$ line profile from a Be disc model, tied to the asymmetric disc during a KL oscillation, and shows that polarization position angle, line peaks, and photometry oscillate with the disc's changing inclination. If these predictions hold, they give observers a concrete route to recognize KL oscillations and torn discs in real Be binaries, and they identify which binaries should be dynamically quiet.","feed_headline":"Binary mass and viscosity decide when Be discs tear","feed_subtitle":"Simulations map where Kozai-Lidov wobbles and torn discs appear, and how to spot them in lines and interferometry.","key_machinery":"The machinery is the three-dimensional smoothed-particle-hydrodynamics (SPH) code used for the disc models, with the Shakura-Sunyaev viscosity parameter $\\alpha_{\\rm ss}$ converted into the SPH artificial viscosity through $\\alpha_{\\rm sph} = 10\\,\\alpha_{\\rm ss} H/h$ (with $\\beta_{\\rm sph}=0$), paired with the HDUST non-local-thermodynamic-equilibrium Monte Carlo radiative-transfer code that converts the particle distributions into synthetic spectra, polarization, images, and interferometric visibilities. A second piece of machinery is the comparison between two viscosity prescriptions: constant $\\alpha_{\\rm ss}$ and constant $\\alpha_{\\rm sph}$. Because the simulations have $H/h < 1$, the constant-$\\alpha_{\\rm sph}$ runs correspond to an effective $\\alpha_{\\rm ss}$ above the nominal 0.1, which inflates the viscous torque and is why those runs fail to reproduce disc tearing; the constant-$\\alpha_{\\rm ss}$ runs are therefore the ones used to map the parameter boundaries.","core_discovery":"The central claim is that in a Be star with a misaligned binary companion, the binary mass ratio $q$ and the disc viscosity parameter $\\alpha_{\\rm ss}$ decide whether the disc tears apart or undergoes Kozai-Lidov oscillations. For $q=0.5$ the paper finds the same evolution as the equal-mass case, only delayed, whereas for $q=0.1$ neither phenomenon occurs for the examined parameters, and the disc simply grows larger and precesses; this is consistent with the analytical tearing radius moving to the disc's outer edge. Raising $\\alpha_{\\rm ss}$ from 0.1 to 0.5 or 1.0 damps the oscillations and suppresses tearing because the faster viscous communication keeps the disc intact. The paper then reports the first radiative-transfer-generated triple-peaked H$\\alpha$ line profile in a Be disc model, produced when the KL oscillation makes the disc strongly asymmetric, and shows that the observables—equivalent width, V/R ratio, peak separation, $V$-band magnitude, and polarization degree and position angle—oscillate in step with the changing disc inclination. Finally, it predicts that a gap in the disc, from either tearing or KL oscillations, creates hump-like bumps in the squared-visibility versus baseline curve and can increase the differential phase across H$\\alpha$ and Br$\\gamma$ lines by up to an order of magnitude, making the dynamical state recognizable with long-baseline interferometers.","pith_inferences":["The boundary between $q=0.1$ and $q=0.5$ suggests that many Be/X-ray binaries with low-mass neutron-star companions may lie below the tearing threshold, so superorbital variability from tearing should be rarer than the equal-mass case would suggest (an inference, not a paper claim).","Because the triple-peaked profile is attributed to disc asymmetry rather than to KL oscillations specifically, the same line-formation mapping may help interpret triple peaks seen in other Be stars with spiral waves or outbursts, if their discs have similar geometry (editorial extension).","The $H/h < 1$ issue implies that past constant-$\\alpha_{\\rm sph}$ simulations of Be discs may have effectively used a higher viscosity than intended; re-scaling their results with the effective $\\alpha_{\\rm ss}$ could shift which systems are expected to tear (editorial inference).","The predicted relationship between visibility hump location and gap radius could be inverted: a single observation with two position angles might estimate the gap's size and orientation, turning interferometry into a diagnostic of the tearing state (editorial suggestion)."],"forward_implications":["Be stars in misaligned binaries with a low-mass companion ($q \\sim 0.1$) are predicted not to show KL oscillations or disc tearing, so any observed such variability in those systems would point to a different mechanism or different parameters.","Discs with $\\alpha_{\\rm ss} \\gtrsim 0.5$ should dissipate smoothly without tearing or KL oscillations, which can be tested by comparing stars with known high-viscosity discs.","Triple-peaked H$\\alpha$ profiles can be produced by an asymmetric disc during a KL oscillation, and the line-formation mapping identifies the slow-moving, radially extended side and fast-moving, small side of the disc as the sources of the separate peaks.","A long-baseline interferometer whose projected baseline crosses a torn or eccentric gap should see a hump in the squared visibility versus baseline curve and an enhanced differential phase across H$\\alpha$ and Br$\\gamma$, a signature that could be searched for in Pleione.","The quantitative criteria (tearing radius and KL timescale) used to interpret the simulations can be applied to observed binaries to estimate whether their discs are in the tearing/KL regime or the quiet regime."],"supporting_citations":[{"why":"Paper I; the equal-mass simulations this work extends and compares against.","marker":"Suffak et al. 2022"},{"why":"Established the SPH-to-HDUST interface and observable predictions for a tearing disc, whose methods are reused here.","marker":"Suffak et al. 2024"},{"why":"Provides the tearing-radius criterion used to explain the absence of tearing at q=0.1.","marker":"Doğan et al. 2015"},{"why":"Supplies the precession rate formula used to check measured precession periods.","marker":"Larwood et al. 1996"},{"why":"Gives the global KL timescale estimate for discs compared to measured periods.","marker":"Martin et al. 2014b"},{"why":"Provides the factor of 2.7 that reconciles predicted and measured initial KL periods.","marker":"Martin & Franchini 2019"},{"why":"The HDUST radiative-transfer code used to compute all synthetic observables.","marker":"Carciofi & Bjorkman 2006"},{"why":"Adapted the SPH code to Be decretion discs and introduced the alpha_sph-alpha_ss conversion.","marker":"Okazaki et al. 2002"},{"why":"Defines the alpha_ss viscosity parameter that is varied across models.","marker":"Shakura & Sunyaev 1973"},{"why":"The artificial viscosity prescription that Eq. (2) tunes to mimic alpha_ss.","marker":"Monaghan & Gingold 1983"}],"fun_headline_variants":["Binary mass ratio and viscosity shape Be disc tearing","Be disc tearing depends on binary mass and viscosity","Mass ratio and viscosity set Be disc Kozai-Lidov fate","For Be stars, mass ratio and viscosity decide disc tears"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The main load-bearing assumption is that the formula converting the Shakura-Sunyaev viscosity into the simulation's artificial viscosity ($\\alpha_{\\rm sph} = 10\\,\\alpha_{\\rm ss} H/h$ with $\\beta_{\\rm sph}=0$) is faithful, since in these simulations $H/h < 1$ makes the effective viscosity higher than the nominal $\\alpha_{\\rm ss}$, and an error in that mapping would shift the mass-ratio and viscosity boundaries at which tearing and KL oscillations occur.","fun_headline_variants_meta":{"raw":{"variants":["Binary mass ratio and viscosity shape Be disc tearing","Be disc tearing depends on binary mass and viscosity","Mass ratio and viscosity set Be disc Kozai-Lidov fate","For Be stars, mass ratio and viscosity decide disc tears"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000741,"raw_usage":{"total_tokens":3402,"prompt_tokens":1137,"completion_tokens":2265,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":753,"completion_tokens_details":{"reasoning_tokens":2200}},"tokens_in":753,"tokens_out":2265,"duration_ms":16767,"temperature":1.0,"reasoning_tokens":2200,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T21:33:24.636543+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the constant-$\\alpha_{\\rm sph}$ equal-mass, 40-degree misalignment model with particle resolution high enough that $H/h \\approx 1$ throughout the disc; if disc tearing then appears, the suppression of tearing in constant-$\\alpha_{\\rm sph}$ models is an artifact of the viscosity conversion rather than a physical effect of the higher viscosity.","supporting_citations":[{"cited_title":"W., Jones C","cited_arxiv_id":null,"evidence_quote":"Established the SPH-to-HDUST interface and observable predictions for a tearing disc, whose methods are reused here."},{"cited_title":"C., Bjorkman J","cited_arxiv_id":null,"evidence_quote":"The HDUST radiative-transfer code used to compute all synthetic observables."}],"review_version":1}