{"id":"d8a6e588-249f-4a7f-b394-a3afdc8c767f","arxiv_id":"2512.19542","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Combining Buchdahl and Fonarev transformations yields the first exact non-vacuum, asymptotically FLRW, axisymmetric black/white-hole solution of the Einstein–scalar system, with horizons located by the mean curvature vector.","lead":"This paper constructs the first exact, axisymmetric, time-dependent black/white-hole solution of general relativity with a scalar field, embedded in an expanding or contracting cosmology. The authors combine two known solution-generating methods and use the mean curvature vector to locate the dynamical trapping horizons.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Extended Fonarev proof has an algebraic inconsistency at (3.45b)-(3.48): the Buchdahl profile gives ξ0ξ1P=ξ0, not 1, so the theorem as printed is not established; direct substitution is needed to confirm the new solution.","rationale":"The reader's weakest assumption correctly identifies the algebraic spot in the proof of the extended Fonarev theorem. My own reduction confirms that the printed condition ξ0ξ1P=1 is inconsistent with the Buchdahl profile: for φ=ξ0 ln(g^aa), the proportionality constant is P=ξ0/β, so ξ0ξ1P=ξ0 (using the stated ξ1=β/ξ0), not 1. The equation (3.45b) itself, when reduced carefully, gives ˙μ=κ(ξ0ξ1/β)˙Ψ=κ˙Ψ, so the intended result Ψ=μ/κ is recovered if the condition is corrected to ξ1P=1. This is therefore a typographical/algebraic presentation error at the load-bearing step, not a demonstrated invalidity of the final solution. The exactness of (4.4)-(4.6) should nevertheless be confirmed by direct substitution, since the paper does not provide such an independent check and the proof as printed is not logically complete. The horizon-locus analysis in §4.2.3 also relies on undocumented numerics, but that is secondary to the existence of the solution. Given the supporting limits and the likely fixability of the proof step, the appropriate verdict remains CONDITIONAL: the authors should supply the corrected derivation or an explicit computer-algebra verification.","tokens_in":30927,"tokens_out":29079,"duration_ms":246089,"concrete_test":"Use a computer algebra system (e.g., xAct or GRTensor) to substitute the metric (4.4)-(4.6) and scalar (4.5) into the full Einstein-scalar equations (3.12)-(3.13) with the parameter constraints (3.26), for generic M, δ, β, C, and check that all field equations and the Klein-Gordon equation vanish identically. If they do, the proof typo is harmless and the solution is exact; if not, the central construction fails. As a second check, re-derive (3.45b) with the corrected condition ξ1P=1 and confirm that (3.48)-(3.51) follow.","verdict_should_be":"UNCHANGED","load_bearing_attack":"At §3.2.2, after (3.45b), the proof defines ∂_kφ = P(g^aa)^{-β}∂_k(g^aa)^β and claims (3.45b) becomes ˙μ=κξ0ξ1P˙Ψ, then imposes ξ0ξ1P=1 to get Ψ=μ/κ. For the Buchdahl seed φ=ξ0 ln(g^aa), one has ∂_kφ=(ξ0/β)(g^aa)^{-β}∂_k(g^aa)^β, so P=ξ0/β. The equation (3.45b) actually reduces to ˙μ=κ(ξ0ξ1/β)˙Ψ, and with the stated ξ1=β/ξ0 this is ˙μ=κ˙Ψ, yielding Ψ=μ/κ. Thus the condition should be ξ1P=1 (or ∂_kφ=ξ0P(...) with P=1/β); the printed ξ0ξ1P=1 is inconsistent because for the Buchdahl profile ξ0ξ1P=ξ0≠1. As written, the proof of the central solution-generating theorem does not go through; since (4.4)-(4.6) is obtained by applying this theorem, the exactness of the new solution is not established by the text. The massless HMN limit and δ=1 Schwarzschild limit provide supporting evidence, but the central proof step needs correction or independent verification.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a solution-generating technique for the self-interacting Einstein-scalar system, combining Buchdahl-type transformations with a generalized Fonarev map to produce non-stationary, axisymmetric solutions. The central result is the extended Fonarev theorem of §3.2.1, whose proof is presented in §3.2.2. The method is then applied to the static Zipoy–Voorhees seed, yielding the time-dependent metric (4.4)–(4.6) with the scalar (4.5), claimed to be asymptotically FLRW and to describe a dynamical axisymmetric black or white hole. The paper also reviews Anco's mean-curvature vector as a generalization of the Kodama vector and uses its zero-locus to identify trapping and anti-trapping horizons. The novelty rests on the exactness of the new solution and on the reliability of the horizon classification.","tokens_in":31272,"tokens_out":13938,"duration_ms":121026,"significance":"If the construction is valid, this is a useful and nontrivial step: it provides an explicit, non-spherically-symmetric, time-dependent solution of the Einstein-scalar system motivated by primordial black hole physics, and it demonstrates a practical method for locating dynamical horizons beyond spherical symmetry. Strengths of the manuscript include the fully explicit metric and scalar field, parameter constraints derived from the field equations rather than fitted, and the fact that the spherical δ=1 limit and the massless HMN limit reduce to previously known solutions. Credit is also explicitly given to Anco for the MCV construction. However, the significance is contingent on repairing the proof of the generating theorem and on documenting the numerical horizon classification; as it stands, the central exactness claim is not fully established.","major_comments":[{"comment":"The proof of the extended Fonarev theorem contains a factor-ξ0 error. For the Buchdahl seed φ=ξ0 ln(g^aa), one has ∂_kφ=(ξ0/β)(g^aa)^{-β}∂_k(g^aa)^β, so the constant P defined by ∂_kφ=P(g^aa)^{-β}∂_k(g^aa)^β equals ξ0/β. Equation (3.45b) therefore reduces to ˙μ=κξ1P ˙Ψ=κ(ξ0ξ1/β)˙Ψ, not κξ0ξ1P ˙Ψ. With the stated ξ1=β/ξ0, the correct conclusion is ˙μ=κ˙Ψ and Ψ=μ/κ. The printed condition ξ0ξ1P=1 is inconsistent with the given P: for the Buchdahl profile it gives ξ0ξ1P=ξ0, not 1. This step is load-bearing because the new solution (4.4)–(4.6) is generated through this theorem. The theorem may be salvageable by removing the superfluous ξ0 or redefining P, but as written the proof does not establish the exactness of the solution. Please correct the proof or provide an independent direct substitution of (4.4)–(4.6).","section":"§3.2.2, Eqs. (3.45b)–(3.48)"},{"comment":"The stated condition R<0 for past anti-trapping horizons is inconsistent with Eq. (4.36). For the C<0 branch with t<0 and ξ2<0, Eq. (4.36) reads 2ξ2/t=+R. Since the left-hand side is positive, one must have R>0, not R<0. The text's statement that 'such a horizon is defined only for R(r_h,θ_h)<0' corresponds instead to the θ+=0 condition with ξ2<0. This sign error affects the white-hole branch analysis and should be corrected.","section":"§4.2.2, Eq. (4.36) and following paragraph"},{"comment":"The identification of black-hole versus cosmological horizon segments relies on the numerical sign of L_{l_-}θ_+ evaluated along the horizon-locus curves, but the numerical procedure is not described. No algorithm, grid parameters, precision, or representative values of the critical points (t_1,r_1) and (t_2,r_2) are provided, and the figure alone is not sufficient to reproduce the computation. Since the physical interpretation — black hole vs contracting cosmological horizon, and horizon production/annihilation — depends on these numerical signs, please provide reproducible numerical data or analytic criteria for the sign of (4.42).","section":"§4.2.3, Fig. 1 and Eq. (4.42)"}],"minor_comments":[{"comment":"The expression for ξ1 appears to be misprinted: β/ξ0 = sqrt(2κβ^2/(1-β^2)), not sqrt(2κβ/(1-β^2)).","section":"Eq. (3.52)"},{"comment":"The acronym for the Husain–Martinez–Nuñez solution is written both as HNM and HMN; please use one consistently.","section":"Throughout"},{"comment":"The Buchdahl transform is applied with a=t, where the seed component is ¯g_tt=-f^δ. Raising a negative metric component to a real power β is branch-dependent, and the scalar profile φ=ξ0 ln(¯g_tt) formally involves the logarithm of a negative quantity. The paper silently writes f^{δβ} and δξ0 ln f. This is a standard convention in FJNW-type solutions, but it should be stated explicitly, e.g., by working with |¯g_aa|, to make the theorem mathematically precise for non-integer β.","section":"§4.2, Eqs. (4.4)–(4.5)"}],"recommendation":"major_revision","confidential_remarks":"The central construction is promising and the manuscript is generally careful about attributing prior work, but the proof of the extended Fonarev theorem needs a concrete algebraic fix and the horizon classification needs better documentation. I do not see this as a fundamental obstruction, hence major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: this paper almost certainly delivers what it claims—a new exact, non-spherical, time-dependent scalar-FLRW compact object—but the proof of the generating theorem contains a fixable algebraic typo, and the horizon analysis relies on numerics the authors don't describe. Worth refereeing, not desk-rejecting.\n\nWhat's genuinely new: combining Buchdahl's static-axisymmetry map with Fonarev's time-dependence map to build an 'extended Fonarev theorem' that takes a static axisymmetric vacuum seed to a dynamical axisymmetric Einstein-scalar solution with Liouville potential. The explicit example, the time-dependent Zipoy-Voorhees-FJNW metric (4.4)-(4.6), is new and the limits check out: β=1 gives ZV vacuum, δ=1 gives Fonarev-FJNW, and the massless HMN limit is consistent. The MCV review in section 2 is a useful advertisement of Anco's 2007 construction, properly credited. The discussion of limitations—contracting vs expanding branch, horizon formation only after critical time—is honest and grounded.\n\nSoft spots, in proportion. The proof of the extended Fonarev theorem has a real typo. At (3.45b) they define ∂_k φ = P(g^aa)^{-β}∂_k(g^aa)^β. For the Buchdahl seed φ=ξ0 ln(g^aa), P=ξ0/β. They then impose ξ0ξ1P=1 to get Ψ=μ/κ. But ξ0ξ1P = ξ0² ξ1/β, and with ξ1=β/ξ0 that equals ξ0, not 1. The correct condition is ξ1P=1, which does give Ψ=μ/κ. So the theorem is not false, but the text as written does not prove it. The fix is small; a referee should ask for a corrected step or an independent substitution check of (4.4)-(4.6).\n\nSecond, the horizon-locus analysis in §4.2.3 relies on 'numerically computed' critical points in Fig. 1 with no procedure described. For a paper that advertises analytic control, that's a gap. The sign of L_{ℓ-}θ_+ flips between critical points—fine, but show the method or at least give the equations solved.\n\nThird, the 'foliation-independent' phrasing in the abstract is stronger than what is shown. For a fixed 2-surface S, the MCV norm is indeed independent of the normal-frame choice; but S itself is defined through a time slicing and a timelike boundary, so the locus of zero norm can depend on that choice. The paper does say 'embedding-dependent' in the abstract, so this is a soft spot, not a fatal flaw.\n\nBottom line: this is a solid paper for the exact-solutions and PBH-modeling communities, with a correctable proof issue and under-documented numerics. I'd send it to a serious referee and ask for the fix, not reject it.","headline":"New exact axisymmetric dynamical scalar-FLRW black hole solution, worth refereeing; the extended Fonarev proof has a fixable factor typo, and the horizon numerics need more detail.","tokens_in":31837,"tokens_out":4507,"would_cite":true,"duration_ms":40801,"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":"This paper constructs the first exact axisymmetric, non-vacuum, asymptotically-FLRW black/white hole solution of the Einstein–scalar system, and locates its dynamical horizons with the mean curvature vector.","keywords":["primordial black holes","exact solutions of Einstein equations","Einstein-scalar system","axisymmetric spacetimes","dynamical horizons","mean curvature vector","Kodama vector","Zipoy-Voorhees metric"],"falsifier":"Substitute the metric (4.4)–(4.6) and the scalar (4.5) directly into the Einstein and scalar field equations with the Liouville potential, verifying every component symbolically under the constraints (3.26); any non-vanishing component away from the known singularities would falsify the solution. Separately, an independent numerical computation of θ+θ− from the null expansions (4.22)–(4.23) should reproduce the reported S-shaped horizon curves and the critical times (t1, r1) and (t2, r2).","tokens_in":30763,"feed_emoji":"🕳️","tokens_out":9166,"duration_ms":86299,"temperature":0.7,"pith_summary":"This paper sets out to establish the first exact, non-vacuum, axisymmetric solution of the Einstein–scalar system that is asymptotically FLRW and contains a dynamical black (or white) hole. To achieve it, the authors extend the Fonarev solution-generating map so it can act on static axisymmetric seeds obtained through the Buchdahl transformation, producing a time-dependent conformal factor and a shifted scalar profile that solve the field equations with a Liouville potential. Their main example is a dynamical Zipoy–Voorhees geometry, and they show that its trapping and anti-trapping horizons are located by the zero-norm surfaces of the mean curvature vector, the natural generalization of the Kodama vector beyond spherical symmetry. A sympathetic reader would care because exact analytic benchmarks like this are rare and are needed to test numerical collapse simulations, model primordial black holes, and develop tools for non-spherical dynamical horizons.","feed_headline":"New exact solution embeds an axisymmetric black hole in a cosmos","feed_subtitle":"It also supplies a coordinate-free way to locate the hole's horizons as the universe expands or contracts.","key_machinery":"The load-bearing object is the extended Fonarev theorem, which combines two transformations: the Buchdahl map turns a static axisymmetric vacuum seed into an Einstein-scalar seed by raising the ignorable-coordinate metric component to powers and taking the scalar proportional to its logarithm; the Fonarev step then multiplies the whole metric by e^{2μ(a)}, with μ(a)=ξ2 ln(Ca+B), and shifts the scalar by (ξ1/κ) μ(a), so the pair solves the field equations with a Liouville potential V=V0 e^{ξ3 φ} subject to the constraints (3.26). The second tool is the mean curvature vector, built from the traces of the extrinsic curvatures of the two normals to a closed 2-surface; its norm is proportional to","core_discovery":"The central claim is that the time-dependent Zipoy–Voorhees metric (4.4)–(4.6) with scalar (4.5) is an exact solution of the Einstein–scalar system with a Liouville potential under the constraints (3.26), and that it is the first exact non-vacuum asymptotically FLRW axisymmetric black/white hole solution of that system. The solution reduces to FLRW at large r, has a time-like singularity at r=2M (ring-like for δ>1) and a space-like singularity at t=0, and its apparent horizons are given by the zero-locus of the mean-curvature-vector norm (4.20), equivalently θ+θ−=0. In the contracting branch (C<0, t<0, ξ2>0) the geometry contains a future trapping/black-hole horizon; in the expanding branch","pith_inferences":["If the theorem is correct, the same two-step construction should work for other static axisymmetric vacuum seeds in the Weyl family, offering a way to build a catalog of dynamical compact objects with different multipole structures; this is not demonstrated in the paper.","The MCV zero-norm criterion suggests a practical numerical recipe — compute the norm from the two-surface normals and track its zero-set — which may be more robust than conventional foliation-by-foliation apparent-horizon searches in non-spherical collapse.","An independent symbolic substitution of (4.4)–(4.6) and (4.5) into the field equations would settle the construction, since the proof's integration step fixes the scalar profile through a proportionality condition that the paper does not exhibit explicitly.","A natural testable extension is to use the MCV's zero-locus as the basis for a quasi-local compaction function and to compute primordial-black-hole formation thresholds beyond spherical symmetry; the paper announces this as a companion project."],"forward_implications":["The map generates an entire family of asymptotically FLRW axisymmetric scalar-sourced spacetimes: any static axisymmetric vacuum seed satisfying the Buchdahl conditions gives a dynamical solution, with the Zipoy–Voorhees geometry as one representative.","For the explicit solution, horizon trajectories are given analytically by 2ξ2/t = ±R(r,θ), and the distinction between black-hole and cosmological horizons is fixed by the sign of ξ2 together with the Lie derivative of the null expansions.","The black-hole branch lives in a contracting FLRW background; the time-reversed branch is a white hole in an expanding background, so the solution provides an exact analogue of black-hole formation in a contracting universe.","The mean curvature vector provides a foliation-independent diagnostic for (anti-)trapped regions that extends the Kodama construction and applies to non-spherical dynamical geometries generally.","The metric is an analytic benchmark that can be used to test numerical collapse codes and to study dynamical-horizon thermodynamics and semi-classical evaporation."],"fun_headline_variants":["First exact axisymmetric black hole in an expanding universe","Exact solution for a spinning black hole in a dynamic cosmos","New solution: black hole embedded in FLRW cosmology","Dynamical axisymmetric black hole solution found","Axisymmetric black hole exact solution in a evolving universe"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The whole construction depends on the proof of the extended Fonarev theorem being exactly right: the step that fixes the scalar profile in terms of the metric deformation must be a genuine solution of the field equations, because if the constraints (3.26) are not satisfied by an independent substitution, the metric (4.4)–(4.6) does not solve Einstein's equations and the black-hole interpretation fails.","fun_headline_variants_meta":{"raw":{"variants":["First exact axisymmetric black hole in an expanding universe","Exact solution for a spinning black hole in a dynamic cosmos","New solution: black hole embedded in FLRW cosmology","Dynamical axisymmetric black hole solution found","Axisymmetric black hole exact solution in a evolving universe"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000164,"raw_usage":{"total_tokens":1135,"prompt_tokens":849,"completion_tokens":286,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":593,"completion_tokens_details":{"reasoning_tokens":209}},"tokens_in":593,"tokens_out":286,"duration_ms":4335,"temperature":1.0,"reasoning_tokens":209,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T14:41:31.280322+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Substitute the metric (4.4)–(4.6) and the scalar (4.5) directly into the Einstein and scalar field equations with the Liouville potential, verifying every component symbolically under the constraints (3.26); any non-vanishing component away from the known singularities would falsify the solution. Separately, an independent numerical computation of θ+θ− from the null expansions (4.22)–(4.23) should reproduce the reported S-shaped horizon curves and the critical times (t1, r1) and (t2, r2).","supporting_citations":[],"review_version":1}