{"id":"76999462-448f-43d8-9493-f5b58629a066","arxiv_id":"2411.09112","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A two-star extension of the earlier compact-component tidal theory shows that near critical curves with zero time-averaged apsidal precession, the apsidal angle librates and spin-orbit inclinations can undergo significant oscillations.","lead":"This paper derives equations for how the orbits and spin axes of two non-compact stars in a close binary change over time when the spins are tilted relative to the orbit. It identifies special critical curves where the apsidal precession rate vanishes, near which the apsidal angle can librate and the spin-orbit inclination can oscillate.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The two-component torque model assumes, without derivation, that IP's one-sided torque expressions and the perpendicularity dS_k/dt⊥S_k carry over independently to each non-compact component.","rationale":"I read this as a theory paper whose central claim is conditional: near critical curves the parallel tidal torque drives libration of the apsidal angle and oscillations of the spin-orbit inclinations, with the pendulum equation (82) describing the idealized regime satisfying Eq. (84). The internal algebra from Eq. (25)-(30) to Eq. (82) is coherent, and Case A provides a numerical check of the pendulum picture where its assumptions hold. The weakest point is not the algebra but the physical input: the torque model. Section 2.4 states Eq. (17) as an assumption, and Appendix B imports IP's torque expressions without re-deriving them for two deformable components. In IP the companion is a point mass, so there is only one non-spherical body and no question of superposing two tidal responses or of cross-torques between the two stars' deformation fields. Here both components are non-compact, and the claim that the same expressions apply independently is an extra assumption. If it is wrong, the entire equation system on which the critical curves and libration results are built is wrong. I agree with the reader's identification of this as the weakest assumption. I do not think the torque model is necessarily incorrect; leading-order additivity is plausible in the weak-tide limit, and the mutual quadrupole interaction may indeed be second order. But the paper does not show that, and the central claim depends on it. A direct linearized calculation of the two-body torque would settle the issue. The other concern raised by the reader, the violation of Eq. (84) in Case B, is explicitly acknowledged in Section 4.4.3 and limits the analytic description rather than the full numerical system, so I treat it as secondary. My verdict remains CONDITIONAL as the reader had it; I do not move it because the concern reinforces an already-conditional recommendation rather than overturning the paper's internal consistency.","tokens_in":31965,"tokens_out":24010,"duration_ms":349633,"concrete_test":"Perform an independent linearized calculation of the tidal torque on one component in a two-star system where both stars have finite tidal and rotational deformations, using the same normal-mode expansions as Appendix A for both stars, and compare the resulting T∥,k and T⊥,k with Eqs. (B1)-(B2) together with Eq. (17). Specifically check (i) whether terms proportional to Q_eq,1 Q_eq,2 appear at the same order as the retained terms, and (ii) whether S_k · dS_k/dt vanishes identically. If cross terms appear at leading order or if the along-spin component is nonzero, the central claim fails; if both checks pass, the inherited torque model is justified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.4 adopts Eq. (17), dS_k/dt = T_k with T_k perpendicular to S_k, so |S_k| is conserved, and the torque components T∥,k, T⊥,k are taken verbatim from IP (Appendix B, Eqs. (B1)-(B2)), where the companion was a compact point-like object. In the present two-non-compact system, the paper assumes that these one-sided IP torques apply independently to each star and that there are no additional spin-spin or cross-tidal contributions. This is not derived. If the tidal response of star 1 depends on the deformation of star 2 beyond the point-mass monopole, or if the torque on a star has a component along S_k so that |S_k| changes, then Eqs. (25)-(30), the critical-curve condition, and the pendulum Eq. (82) are not the correct dynamical system. The paper provides no calculation isolating the two-body linear response and no test of the perpendicularity assumption in the two-deformable-body context; it only cites the compact-companion derivation. Since every subsequent result, including the claimed libration near critical curves, is a consequence of this assumed torque structure, this is the most load-bearing unverified step. The timescale limitation in Eq. (84) is real but explicitly acknowledged and affects only the analytic pendulum description, not the full system; the torque foundation affects everything.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the non-dissipative tidal torque formalism of Ivanov & Papaloizou (IP, IP1, IP2) from binaries with one compact component to binaries in which both components are non-compact and carry spin angular momentum. The authors write down a system of ordinary differential equations for the orbital angular momentum, the spin vectors, and the apsidal angle, using torque components taken from IP for each component. In the limit of small spin-to-orbit angular momentum ratio, they reduce the system to a smaller set of equations, identify critical curves in parameter space where the time-averaged apsidal precession rate vanishes, and show that near these curves the apsidal angle can librate, described by a forced pendulum equation. They compute critical curves for parameters appropriate to DI Her using MESA/GYRE stellar models and present numerical integrations for two representative cases, finding agreement with the analytic theory in one case (A) and qualitative libration/circulation behavior in the other (B).","tokens_in":32190,"tokens_out":25910,"duration_ms":257428,"significance":"If the underlying torque model is accepted, this is a valuable extension because real close binaries such as DI Her have two non-compact components. The paper recovers the compact limit of IP1, provides transparent derivations of the angle evolution equations, and uses independent stellar-structure calculations for the parameters. The analytic pendulum description gives a falsifiable prediction: near critical curves, libration of the apsidal angle and oscillations of spin-orbit inclinations should occur, and in the small-spin limit such curves require a retrograde component. The numerical comparison in Case A is a good consistency check of the secular reduction. The main caveat is that the torque expressions are assumed to carry over independently from the compact-companion calculation, and the simulations do not test this assumption.","major_comments":[{"comment":"The paper assumes that the torque on each component has the form dS_k/dt = T_k with T_k perpendicular to S_k, and that the components T_parallel,k and T_perpendicular,k are exactly the expressions derived in IP for a binary with one compact companion. This is not derived for the two-non-compact case. If the tidal response of one star is affected by the deformation of the other beyond the point-mass monopole, or if the torque has a component along S_k so that |S_k| changes, then Eqs. (25)-(30), the critical-curve condition, and the pendulum Eq. (82) are not the correct dynamical system. I request an explicit ordering argument or an order-of-magnitude estimate showing that the neglected cross-tidal and along-spin torque contributions are higher order in (R/a)^5 or (Omega_r/Omega_*)^2, or a statement of the approximation regime in which the IP expressions are valid for each component.","section":"Section 2.4, Eq. (17), Appendix B"},{"comment":"The numerical simulations integrate the same assumed torque model, so they validate the reduction from the full ODE system to the secular and pendulum descriptions, but they cannot validate the torque model itself. The agreement in Case A is evidence that the secular theory correctly represents the assumed equations, not that the assumed equations are the true two-deformable-body dynamics. The abstract and conclusions should be phrased as conditional on the torque model, and the paper should explicitly state that the torque additivity and perpendicularity are inputs rather than outputs of the numerical comparison.","section":"Sections 3.3 and 4.4"}],"minor_comments":[{"comment":"The statement that the semi-major axis is 'conserved' is not consistent with Eq. (25), which allows dL/dt to be nonzero. In the small-spin limit the change is negligible, but the statement should be qualified as an approximation.","section":"Section 5.1"},{"comment":"The discussion of the possibility that only the secondary of DI Her is retrograde should be connected to the condition derived in Section 3.3, which requires the component with the larger perpendicular spin projection to be retrograde for a critical state to exist.","section":"Section 4.3.1"},{"comment":"The phrase 'at least one component with retrograde rotation' is imprecise; the derived condition is that the component that dominates the averaged non-inertial apsidal precession must be retrograde, not merely that some component is retrograde.","section":"Abstract"},{"comment":"In the typeset version, Eq. (26) and Eq. (27) appear to have unbalanced parentheses in the right-hand sides; please re-check the formulas to ensure they are typeset correctly.","section":"Equations (26) and (27)"}],"recommendation":"major_revision","confidential_remarks":"The main risk to the paper is the unproven additivity of the IP torque expressions. I do not think this is a fatal flaw, because a leading-order argument likely suffices, but the authors should provide it. The application to DI Her is speculative but clearly labeled. The paper is a good fit for MNRAS, and the numerical comparison in Case A is convincing within the assumed model."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper is IP plus a second spinning star, done with the same level of care as the earlier papers. The genuinely new content is the two-spin generalization: both spin-orbit inclinations evolve under parallel and perpendicular torques, the relative azimuthal angle nu_1 - nu_2 enters as a dynamical variable, and the critical-curve/pendulum analysis is carried out for two non-compact components. There is a real result here: near critical curves the tidal torque drives apsidal libration and inclination oscillations; in the small-spin limit the critical curves exist only when at least one component rotates retrograde.\n\nWhat is done well: the derivation of the governing equations is transparent. The compact limit reduces to IP1, checked explicitly. The analytic secular theory is compared with direct integration in Case A and the libration/circulation boundary tracks E/E_crit as predicted. The MESA/GYRE inputs are described in enough detail that the stellar parameters are reproducible, and the convergence of the mode sums with N_r is shown.\n\nThe soft spots are real but smaller than the stress-test note implies. The torque on each star is taken from IP, where the companion was a point mass. Applying those expressions independently to each component assumes that the two tidal responses add linearly and that each spin magnitude is conserved. The conserved spin magnitude is standard for conservative non-dissipative torques; the additivity is the natural leading-order weak-tide approximation, since quadrupole-quadrupole terms are higher order. But the paper never states this, and a referee should ask for a short justification or an estimate of the neglected cross-tidal terms. The more concrete limitation is the timescale separation |Delta_omega| >> Omega_parallel behind the pendulum description, eq. (84). Case B violates it, and the authors are upfront about that; the analytic model is tested in Case A, where it works.\n\nMinor: no code or data is shipped, only 'available on reasonable request.' That is acceptable for a theory paper but it would help a referee check the torque additivity assumption.\n\nBottom line: this deserves a serious referee and, after a modest revision, publication. If I were working on misaligned binaries, I would cite it.","headline":"A careful two-spin extension of the IP framework; the main unproven step is real but likely not fatal, and the paper deserves a serious referee.","tokens_in":32765,"tokens_out":3138,"would_cite":true,"duration_ms":39056,"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":"Tidal torque can make the apsidal angle librate instead of precess in a binary of two ordinary stars with misaligned spins.","keywords":["tidal interaction","apsidal precession","spin-orbit misalignment","non-dissipative torque","critical curves","apsidal libration","binary stars","DI Her"],"falsifier":"A direct numerical computation of the tidal response of two extended, deformable stars in a misaligned eccentric binary, without assuming each star feels only the point-mass field of the other, would settle whether the torque decomposition is valid; a mismatch in the magnitude or phase of the parallel torque would shift the critical curves and change the libration predictions.","tokens_in":31711,"feed_emoji":"💫","tokens_out":6026,"duration_ms":57497,"temperature":0.7,"pith_summary":"This paper extends the theory of non-dissipative tidal evolution in misaligned binaries from the case of one compact, spinless companion to two ordinary non-compact stars, each carrying its own spin angular momentum. It argues that the parallel component of the tidal torque, arising from inertial and Coriolis forces on the tidal response, comes to dominate the orbital dynamics near critical curves in parameter space where the time-averaged apsidal precession rate vanishes. Near these curves the apsidal angle librates rather than circulating, and the spin-orbit inclination angles undergo bounded oscillations whose amplitude is set by the distance from the critical curve. In the small-spin limit, critical curves exist only for binaries with at least one retrograde component. The paper provides the full governing equations, a forced-pendulum approximation for the near-critical motion, and numerical simulations for a DI-Her-like system that reproduce the analytic libration-circulation transition.","feed_headline":"Tidal torque makes apsidal angle librate in misaligned binaries","feed_subtitle":"Near critical curves, the orbit can swing between prograde and retrograde states while spin inclinations oscillate.","key_machinery":"The central object is the parallel tidal torque $T_{\\parallel,k}$ acting on each component, defined through eq. (17) with the explicit form in eq. (B1). This torque depends on the apsidal angle $\\hat\\varpi_k=\\varpi+\\gamma_k$ and arises from inertial and Coriolis forces acting on the tidal bulge. The critical curves are the loci in parameter space where the time-averaged apsidal precession rate $\\dot\\varpi_T+\\dot\\varpi_E+\\dot\\varpi_R+\\dot\\varpi_{NI}$ vanishes; near these curves the system reduces to the forced simple pendulum of eq. (82), whose pendulum frequency $\\Omega_\\parallel$ is given by eq. (83), predicting both librating and circulating solutions separated by a separatrix.","core_discovery":"The central claim is that for a binary with two non-compact components, the long-term evolution of the apsidal angle and of the spin-orbit inclinations is controlled by the parallel tidal torque, and that the most interesting behaviour occurs on critical curves where the time-averaged total apsidal precession rate is zero. On those curves the apsidal angle satisfies a forced-pendulum equation, so it can librate between prograde and retrograde orientations instead of steadily precessing, while the obliquities oscillate with amplitude bounded by the distance from the critical curve. In the small-spin limit such critical states require at least one retrograde component, because the non-inertial spin-precession contribution to the apsidal rate changes sign for retrograde spins. Numerical integrations for parameters similar to DI Her confirm the expected transition between librating and circulating solutions, and show that near-polar, rapidly rotating systems can still exhibit libration even when the timescale separation assumed by the pendulum model breaks down.","pith_inferences":["Beyond the paper: the same pendulum mechanism should apply to star-planet systems where the star is non-compact and the planet can be treated as a point mass, so the critical-curve condition may explain observed spin-orbit misalignments in close-in exoplanet systems that linger near apsidal resonances.","Beyond the paper: a direct hydrodynamical simulation of two extended, deformable stars with misaligned spins would test whether the torque decomposition into two independent point-mass-companion expressions remains valid when the two tidal responses can interact.","Beyond the paper: if apsidal libration near critical curves is common, the distribution of apsidal angles in a large sample of misaligned eclipsing binaries should cluster near the librating phase rather than being uniformly distributed.","Beyond the paper: the requirement of at least one retrograde component in the small-spin limit creates a selection effect: binaries with anomalously slow measured apsidal precession may be those currently undergoing near-critical libration, so their observed $\\dot\\varpi$ is systematically smaller than the steady-precession prediction."],"forward_implications":["Binary systems whose parameters lie near a critical curve will show libration of the apsidal line rather than steady precession, so the orbit can oscillate between prograde and retrograde apsidal states.","Spin-orbit inclination angles undergo bounded oscillations whose amplitude is controlled by the distance to the critical curve, and near the separatrix high-frequency forcing can generate chaotic motion.","In the small-spin regime, a critical state is possible only with at least one retrograde component, so the predicted apsidal libration is a signature of retrograde rotation in the binary.","For a close, eccentric binary like DI Her, the measured apsidal precession rate and spin inclinations can be used to assess proximity to a critical curve; if the secondary is slightly retrograde, libration about a critical curve is a realistic outcome.","The closed system of equations (25)-(30) provides a directly integrable non-dissipative model for arbitrary eccentricity and two spinning non-compact stars, replacing the earlier compact-companion restriction."],"supporting_citations":[{"why":"Derived the non-dissipative tidal torques, including the parallel torque from inertial and Coriolis forces, whose expressions the paper adopts for each component.","marker":"Ivanov & Papaloizou (2021) (IP)"},{"why":"Introduced the critical-curve concept and the pendulum description for the one-compact-component case that this paper generalises to two non-compact stars.","marker":"Ivanov & Papaloizou (2023a) (IP1)"},{"why":"Provided the preliminary numerical simulations of libration near critical curves that the present paper extends and compares with.","marker":"Ivanov & Papaloizou (2023b) (IP2)"},{"why":"Supplies the general precession formalism for the orbital angular momentum and eccentricity vector that underlies the apsidal precession contributions.","marker":"Barker & O'Connell (1975)"},{"why":"Gives the classical tidal-distortion apsidal precession contribution adopted as $\\dot\\varpi_T$.","marker":"Sterne (1939)"},{"why":"Provides the spin-orbit inclination estimates for DI Her used to set the system parameters for the numerical study.","marker":"Philippov & Rafikov (2013)"},{"why":"Supplies updated DI Her parameters, including masses, radii, rotation periods, inclinations, and the apsidal angle, used for the stellar models and case B.","marker":"Liang, Winn, & Albrecht (2022)"}],"fun_headline_variants":["Tidal torque swings apsidal angle in misaligned binaries","Critical curves flip orbit orientation in misaligned binaries","Apsidal libration emerges near critical curves in binaries","Orbit can flip prograde-retrograde under tidal torque","Misaligned binaries show orbit flips near critical curves"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the torque expressions derived for a binary with one compact, spinless companion apply independently to each non-compact star, so that the total non-dissipative torque is simply the sum of the two identical functional forms applied separately.","fun_headline_variants_meta":{"raw":{"variants":["Tidal torque swings apsidal angle in misaligned binaries","Critical curves flip orbit orientation in misaligned binaries","Apsidal libration emerges near critical curves in binaries","Orbit can flip prograde-retrograde under tidal torque","Misaligned binaries show orbit flips near critical curves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00114,"raw_usage":{"total_tokens":4770,"prompt_tokens":1020,"completion_tokens":3750,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":636,"completion_tokens_details":{"reasoning_tokens":3667}},"tokens_in":636,"tokens_out":3750,"duration_ms":21697,"temperature":1.0,"reasoning_tokens":3667,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T21:00:47.720162+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct numerical computation of the tidal response of two extended, deformable stars in a misaligned eccentric binary, without assuming each star feels only the point-mass field of the other, would settle whether the torque decomposition is valid; a mismatch in the magnitude or phase of the parallel torque would shift the critical curves and change the libration predictions.","supporting_citations":[{"cited_title":"B., Papaloizou J","cited_arxiv_id":null,"evidence_quote":"Derived the non-dissipative tidal torques, including the parallel torque from inertial and Coriolis forces, whose expressions the paper adopts for each component."},{"cited_title":"A., Rafikov R","cited_arxiv_id":null,"evidence_quote":"Provides the spin-orbit inclination estimates for DI Her used to set the system parameters for the numerical study."},{"cited_title":"N., Albrecht S","cited_arxiv_id":null,"evidence_quote":"Supplies updated DI Her parameters, including masses, radii, rotation periods, inclinations, and the apsidal angle, used for the stellar models and case B."}],"review_version":1}