{"id":"4e81d81d-117b-4792-a4ec-5f15f5e7a419","arxiv_id":"1908.11258","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A potential-based method constructs stationary, differentially rotating, non-barotropic neutron stars in general relativity, with dynamical evolutions supporting that they are equilibrium configurations.","lead":"This paper builds the first general-relativistic models of neutron stars that rotate differentially and have a non-barotropic (temperature and entropy dependent) equation of state. It wraps the Euler equation into a potential and checks the resulting stars with two independent numerical codes, which matters because real newborn neutron stars are hot and not barotropic.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Method constructs only the entropy profile implied by a chosen Q; it cannot match a prescribed thermal profile, so the claimed potential-formulation extension overstates its scope.","rationale":"The paper is careful and technically sound within its stated construction: Q is chosen first, h and F are obtained as partial derivatives, and the entropy is derived from the EOS inversion. The Euler residuals and the BAM stability tests are appropriate and give real support for the narrow claim that XNS can produce stationary, differentially rotating, non-barotropic configurations in the inverse sense. The reader's conditional verdict is therefore appropriate. The most load-bearing concern is scope, not internal inconsistency: the title and abstract can be read as promising a method for models with prescribed non-barotropic thermal profiles, whereas the construction only yields the thermal profile implied by a chosen Q. This is acknowledged in Sec. III B and Sec. VI B, but the abstract does not carry the qualification. I do not see a fatal mathematical error in the derivation: the Maxwell relation is automatically satisfied along any actual solution because the Euler equation is exactly dQ = dp/h + F dΩ in the two-dimensional stellar domain; the restriction appears only when one tries to prescribe s instead of Q. The concrete test above would settle whether the method can be pushed toward realistic prescribed profiles; absent that, the central claim should be read narrowly as an inverse construction for a restricted family.","tokens_in":21486,"tokens_out":15582,"duration_ms":158152,"concrete_test":"Take a target entropy profile s_target(r,θ) from a BNS merger simulation, discretize it on the XNS grid, and search over Q of the form Eq. (48) with an increasing number of modes for the equilibrium that best reproduces s_target, recording both the L2 error in s and the Euler residual ⟨log|δ|⟩ of Eq. (46). If the achievable L2 error cannot be made small while keeping the Euler residual near the paper's 10^-7 level, the method cannot realize prescribed non-barotropic thermal profiles, and the abstract should be qualified to the inverse family actually constructed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central limitation is stated by the authors themselves in Sec. III B: for fixed pointwise entropy, Eqs. (35)-(37) are three equations in two unknowns and 'in general has no solution'; therefore s(r,theta) is left free and determined after the fact from Eq. (37). The entropy profile is an output of the chosen Q, not an input. Consequently, the abstract's claim that the potential formulation 'can be extended to the non-barotropic case' is not a general extension: for a generic non-barotropic EOS and a desired thermal profile, the integrability condition Eq. (34), ∂(1/h)/∂Ω|_p = ∂F/∂p|_Ω, will fail, so no such Q exists. For any actual equilibrium solution the relation is automatic because dQ = dp/h + F dΩ with Q = -ln(α/γ); the real restriction is prescribability of the physical profile, not existence of a potential. Sec. VI B offers only a proof-of-principle fitting procedure, explicitly limited to planar configurations and stated to be cumbersome. The numerical evidence — Euler residuals and BAM evolutions — supports the narrower inverse construction, not a general method for prescribed non-barotropic models.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a new method for constructing stationary, axisymmetric, differentially rotating relativistic neutron stars with non-barotropic equations of state. The key idea is to write the relativistic Euler equation as the differential of a potential Q(p,Omega), with partial derivatives equal to 1/h and F, generalizing the barotropic first-integral method. After defining the potential Q in Sec. III A and a simple non-separable model Q = Q0 + H(p) + F(Omega) + b H(p) F(Omega) in Sec. III B, the authors implement the scheme in the XNS code and construct seven models with an analytic two-parameter EOS. They validate the models in two ways: the Euler-equation residuals of the consistent models are about 1e-7, two to three orders of magnitude smaller than in the inconsistent 'control' models, and BAM dynamical evolutions show that non-convective consistent models oscillate at a level comparable to a cold rigidly rotating star. The paper also derives consequences such as the relativistic von Zeipel theorem (non-barotropic stars must be differentially rotating) and discusses extensions to more general equations of state, Legendre-transformed potentials, and Newtonian gravity.","tokens_in":21746,"tokens_out":14302,"duration_ms":140214,"significance":"This is a genuinely useful contribution to relativistic stellar structure. The potential formulation Q(p,Omega) is elegant and, to my knowledge, new in the GR context, and the numerical implementation with careful controls provides strong evidence that the constructed configurations are stationary. The paper is unusually honest about limitations: Sec. III B states that the system (35)-(37) is overdetermined for a prescribed entropy profile and that s(r,theta) is derived after the fact, and Sec. VI B presents entropy fitting only as a proof of principle. The Euler-residual test and the BAM stationarity test are well designed: the consistent models give residuals around 1e-7 while the inconsistent C_Omega and C_p controls are two to three orders worse, and the non-convective NN model behaves like the cold rigid rotator in evolution. The paper is a solid proof of concept rather than a tool for arbitrary prescribed thermal profiles.","major_comments":[{"comment":"The abstract claims that the potential formulation 'can be extended to the non-barotropic case' without stating the inverse nature of the construction. As the authors themselves note in Sec. III B, for a prescribed pointwise entropy the system (35)-(37) is three equations in two unknowns and in general has no solution; the entropy profile is an output of the chosen Q, not an input. For a generic tabulated EOS and a desired thermal profile, the Maxwell relation (34) will generally fail, so no Q exists. I recommend rewording the abstract and introduction to say that the method constructs the unique thermal profile consistent with a chosen potential Q, and to state explicitly that matching a prescribed s(r,theta) is not possible in general. The current wording overstates the scope of the extension.","section":"Abstract and Sec. III B"}],"minor_comments":[{"comment":"The symbol h is used both for the specific enthalpy (e.g., Eq. (10)) and for the enthalpy density (e.g., Eq. (4) and Eq. (42)); this is confusing, and the footnote in Sec. II does not cover this distinction. Please use separate symbols or state clearly which quantity is meant in each equation.","section":"Sec. II and Sec. IV A"},{"comment":"Below Eq. (32), the statement that 'given the pair p and Omega, we must be able to determine the pair r and theta' should be qualified as a local condition. The paper mentions the two-hemisphere degeneracy but does not discuss points where grad p and grad Omega are parallel or where Omega is constant on the symmetry axis; please state the assumption that the Jacobian of the map (r,theta) to (p,Omega) is nonzero except on a set of measure zero and explain how the construction behaves near such points.","section":"Sec. III A"},{"comment":"The caption contains what appears to be a typo: 'p, r,'. Please fix the caption.","section":"Fig. 3"},{"comment":"The asterisk footnote says that b was included in a non-consistent way in C_Omega and C_p. This is explained in Sec. IV C, but a parenthetical reference to Eq. (40) versus Eq. (42) would make the construction of the control models clearer.","section":"Table II"},{"comment":"In Eq. (47), the index i in delta_i denotes the direction of differentiation while the average is taken over points j. This is fine, but writing the two explicit averages <log|delta_r|> and <log|delta_theta|> in the text would avoid possible confusion.","section":"Sec. V B"},{"comment":"There are several minor grammatical slips, for example 'the remaining density oscillations is likely' in Sec. IV B, and the phrase 'We have first defined' in Sec. III B. A careful proofread is recommended.","section":"General"}],"recommendation":"minor_revision","confidential_remarks":"The novelty claim for the Newtonian case should be scrutinized by the editor against the baroclinic star literature cited in Refs. [21-26], since some of those methods may already use a potential of the form Q(p,Omega). The main result is nonetheless a useful relativistic extension with solid numerical evidence, and the fit to the journal's scope is good."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the headline: this is the first tool for stationary, differentially rotating, non-barotropic neutron-star models in GR, and the numerical validation is genuinely good. The paper deserves a serious referee. But the abstract oversells the potential-formulation novelty; the construction is an inverse method that chooses Q and then derives the entropy profile, not a method for prescribing a thermal profile.\n\nWhat is actually new: the paper produces the first GR equilibrium models of this class, implements them in XNS, and checks them two independent ways. The Euler-equation residuals for consistent models are ~1e-7, two to three orders better than the deliberately inconsistent control models. The BAM evolutions show that the convectively stable non-barotropic model (NN) oscillates at the same level as the cold rigid model (CR), while controls drift. That is solid evidence the configurations are stationary. The authors also derive the Maxwell relation correctly and are transparent that for fixed pointwise entropy the system is overdetermined (Sec III B), and that fitting a prescribed profile is only a proof of principle, planar, and cumbersome (Sec VI B).\n\nThe soft spots are in proportion. The abstract and conclusion claim the potential formulation 'can be extended to the non-barotropic case'; that is too broad. The extension only works for the entropy profiles that are integrable with the chosen Q. If you want a specific thermal profile from a finite-temperature EOS, the method does not give you it, except through the pointwise fitting procedure the authors themselves say is cumbersome. The same is true for the claim of novelty in the Newtonian case: non-barotropic (baroclinic) stellar models are well established in Newtonian gravity, and the potential Q is a local rewriting of the Euler equation, so the novelty should be restricted to the GR equilibrium construction and its numerical realization. Two minor issues: no code or data are released, and the paper does not discuss degeneracies where grad p and grad Omega are parallel (the map to (p,Omega) is assumed nondegenerate).\n\nThe math, data, and citation pattern look honest. The limitations are stated, not hidden. My recommendation: send it to peer review, and in the revision ask the authors to recalibrate the novelty claims in the abstract and conclusions.","headline":"First GR non-barotropic rotating star models with solid numerical validation, but the advertised potential-formulation extension is narrower than the abstract suggests.","tokens_in":22267,"tokens_out":2862,"would_cite":true,"duration_ms":27800,"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":"The relativistic Euler equation can be cast in potential form for non-barotropic fluids, giving stationary, differentially rotating neutron star models.","keywords":["non-barotropic neutron stars","differentially rotating stars","potential formulation of the Euler equation","relativistic stellar structure","thermal entropy profiles","baroclinic stars","numerical relativity hydrodynamics","von Zeipel theorem"],"falsifier":"Take one of the constructed non-barotropic equilibrium models and evolve it with a fully general-relativistic hydrodynamics code that does not rely on the conformal-flatness approximation for several tens of dynamical timescales; if the central density, entropy map, and angular-velocity profile drift systematically away from their initial values beyond the small oscillations seen in the paper, the configurations are not true stationary solutions. A cheaper check is to evaluate the Euler-equation residual $\\delta_i=\\partial_i Q-\\partial_i p/h-F\\,\\partial_i\\Omega$ at increasing resolution and require it to converge to zero.","tokens_in":21241,"feed_emoji":"🔄","tokens_out":11052,"duration_ms":96594,"temperature":0.7,"pith_summary":"The paper shows that the potential formulation of the relativistic Euler equation, long used to build rotating stellar models under the barotropic assumption, survives when the fluid is non-barotropic: one can introduce a potential $Q(p,\\Omega)$ whose partial derivatives reproduce the enthalpy and angular-momentum terms. With this potential the paper constructs stationary, axisymmetric, differentially rotating neutron stars whose entropy and specific angular momentum depend on both pressure and angular velocity, not on a single variable. The authors validate the construction by evolving the configurations with a dynamical general-relativistic hydrodynamics code: consistently built models stay close to the initial equilibrium, while deliberately inconsistent control models drift away. The result matters because supernova and merger remnants are born hot and non-barotropic, and such equilibrium models give a starting point for studying their structure and evolution.","feed_headline":"Euler equation gains a potential form for non-barotropic neutron stars","feed_subtitle":"A single potential makes the Euler equation integrable for hot stars; dynamical runs confirm the equilibria.","key_machinery":"The central object is the potential $Q(p,\\Omega)$ defined as $-\\ln(\\alpha/\\gamma)$, with pressure $p$ and angular velocity $\\Omega$ as independent variables, in exact analogy to a thermodynamic potential. Its partial derivatives give the specific volume $1/h$ and the angular-momentum function $F$, and at each grid point the construction enforces the three equations $Q(p,\\Omega)=-\\ln(\\alpha/\\gamma)$, $\\partial_\\Omega Q=F$, and $\\partial_p Q=1/h$. The cross term $bH(p)\\mathcal{F}(\\Omega)$ in the model potential is what breaks barotropicity, because it makes $\\partial_p\\partial_\\Omega Q$ nonzero; the Maxwell-like relation then forces the entropy to depend on $\\Omega$ and the specific angular momentum to depend on $p$. This machinery converts the Euler equation from a differential relation into an algebraic system for $(p,\\Omega)$, with the entropy obtained afterward from the equation of state.","core_discovery":"The central claim is that the Euler equation $\\partial_i p/h+\\partial_i\\ln(\\alpha/\\gamma)+F\\partial_i\\Omega=0$ can be integrated through a potential $Q(p,\\Omega)=-\\ln(\\alpha/\\gamma)$ provided $1/h=\\partial Q/\\partial p|_\\Omega$ and $F=\\partial Q/\\partial\\Omega|_p$, which requires the Maxwell-like relation $\\partial_\\Omega(1/h)|_p=\\partial_p F|_\\Omega$. Choosing a potential such as $Q=Q_0+H(p)+\\mathcal{F}(\\Omega)+bH(p)\\mathcal{F}(\\Omega)$, the paper solves for pressure and angular velocity at every point and then derives the entropy profile that makes the one-form $\\mathrm{d}p/h+F\\,\\mathrm{d}\\Omega$ integrable. The resulting stars have $s=s(p,\\Omega)$ and $l=l(p,\\Omega)$, so they are genuinely non-barotropic, and the same potential construction is stated to be new even in the Newtonian limit.","pith_inferences":["A practical fitting algorithm is the natural next step: expand $Q(p,\\Omega)$ in a series of separable terms and adjust coefficients to match a target entropy at selected grid points; the paper only sketches this as a proof of principle.","Because the Maxwell-like relation is a solvability condition, it could be used as a diagnostic on dynamical data: compute $\\partial_\\Omega(1/h)|_p-\\partial_p F|_\\Omega$ on a merger remnant to test whether a stationary non-barotropic equilibrium is a good local approximation, a check the paper does not perform.","The non-barotropicity parameter $b$ need not be constant; promoting it to a function of $p$ and $\\Omega$ would give extra freedom to fit realistic entropy gradients, an extension left implicit.","Since only $(p,h)$ are needed from the Euler equation, the method should combine with tabulated finite-temperature equations of state by inverting $(p,h)$ to temperature and composition, provided the inversion is unique; the paper demonstrates the inversion only for a simple analytic equation of state."],"forward_implications":["A stationary non-barotropic star must be differentially rotating; uniform rotation forces the star to be barotropic, a relativistic version of the von Zeipel theorem.","If the angular-momentum function depends only on $\\Omega$, then the entropy is a function of pressure alone and the star is an effective barotrope; the converse also holds, so non-barotropic thermal structure and a genuinely two-dimensional rotation law are inseparable.","The Euler equation fixes only pressure and enthalpy density directly; every other thermodynamic quantity must come from inverting the equation of state, which the paper demonstrates for a two-dimensional equation of state and outlines for one with an additional variable such as electron fraction.","The potential construction is independent of the conformal-flatness approximation and can be adapted to the full stationary metric or to Newtonian gravity, and the second independent variable could in principle be something other than $\\Omega$, such as a magnetic-field-related coordinate.","The method can produce equilibrium configurations that are convectively unstable, so stability is a separate question to be checked for each model."],"supporting_citations":[{"why":"Supplies the standard equations for stationary axisymmetric rotating relativistic stars, including the Euler equation form the paper starts from.","marker":"[30]"},{"why":"The stationary stellar-structure code that is extended to solve the new non-barotropic matter equations.","marker":"[14, 15]"},{"why":"The modified and validated version of the stationary code on which the implementation is based.","marker":"[16]"},{"why":"The dynamical general-relativistic hydrodynamics code used to evolve the configurations and test their stationarity.","marker":"[39–44]"},{"why":"Supernova and merger simulations showing that realistic remnants are non-barotropic, providing the physical motivation and comparison data.","marker":"[4, 5]"},{"why":"The von Zeipel theorem whose relativistic version underpins the result that non-barotropic stationary stars must be differentially rotating.","marker":"[56, 57]"},{"why":"Provides the differential-rotation law adopted to build the model configurations.","marker":"[31]"},{"why":"Newtonian non-perturbative baroclinic star construction, the baseline showing the potential formulation is new even in the Newtonian case.","marker":"[21]"}],"fun_headline_variants":["Potential form tames Euler equation for non-barotropic stars","Non-barotropic neutron stars: Euler equation now integrable","Barotropic assumption bypassed: new Euler potential","Euler equation integrates for rotating hot neutron stars","First potential method for non-barotropic relativistic stars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a potential $Q(p,\\Omega)$ can be found whose mixed partial derivatives satisfy the Maxwell-like relation $\\partial_\\Omega(1/h)|_p=\\partial_p F|_\\Omega$; for a generic non-barotropic equation of state and a prescribed thermal profile this integrability condition will fail, and then the method cannot reproduce that profile—it can only produce the entropy profile that makes the potential integrable.","fun_headline_variants_meta":{"raw":{"variants":["Potential form tames Euler equation for non-barotropic stars","Non-barotropic neutron stars: Euler equation now integrable","Barotropic assumption bypassed: new Euler potential","Euler equation integrates for rotating hot neutron stars","First potential method for non-barotropic relativistic stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000799,"raw_usage":{"total_tokens":3532,"prompt_tokens":977,"completion_tokens":2555,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":593,"completion_tokens_details":{"reasoning_tokens":2480}},"tokens_in":593,"tokens_out":2555,"duration_ms":21155,"temperature":1.0,"reasoning_tokens":2480,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:21:08.532134+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take one of the constructed non-barotropic equilibrium models and evolve it with a fully general-relativistic hydrodynamics code that does not rely on the conformal-flatness approximation for several tens of dynamical timescales; if the central density, entropy map, and angular-velocity profile drift systematically away from their initial values beyond the small oscillations seen in the paper, the configurations are not true stationary solutions. A cheaper check is to evaluate the Euler-equation residual $\\delta_i=\\partial_i Q-\\partial_i p/h-F\\,\\partial_i\\Omega$ at increasing resolution and require it to converge to zero.","supporting_citations":[{"cited_title":"Disk formation in the collapse of supramassive neutron stars","cited_arxiv_id":"1806.07775","evidence_quote":"Provides the differential-rotation law adopted to build the model configurations."}],"review_version":1}