{"id":"5bff757a-aeda-482c-8306-6360c47141d8","arxiv_id":"2501.11642","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"New analytic and numerical initial data are constructed for Kaluza-Klein spacetimes with a space-dependent compactification radius, covering black strings, naked or hidden KK bubbles, and multiple black strings.","lead":"This paper constructs initial data, valid snapshots for Einstein's equations, of five-dimensional Kaluza-Klein spacetimes in which the extra-dimension radius varies in space, including black strings and bubbles of nothing. The configurations are built to seed future numerical-relativity simulations of an expanding extra-dimensional bubble colliding with black objects.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The numerical data of Sec. V rest on a hand-chosen transition radius r1=3R∞ and outer boundary xout=10 with no sensitivity study or Hamiltonian-constraint residual check, so the reported physical properties of the black-string/KK-bubble slice are not yet established as robust.","rationale":"The analytic sections hold up under scrutiny: substituting the SO(3) forms (14)-(15) into Eq. (7) indeed yields only the relation (16), so the claimed one-parameter family is consistent with the Hamiltonian constraint; the ansatz (36) in Sec. IV correctly reduces Eq. (9) to Laplace's equation for harmonic Psi; and the ADM mass formulas (26), (34), and (41) follow from the stated asymptotics. The weakest point is the numerical construction of Sec. V, exactly as the reader identified. The hand-picked transition radius r1 and outer boundary xout, combined with the absence of any r1- or xout-dependence study and the absence of a constraint-residual diagnostic, mean that the reported physical properties of the black-string/KK-bubble data are not yet robustly established. The author's explicit Sec. VI caveat that evolving a non-SO(3) KK bubble surface is an unsolved challenge further limits the immediate usefulness of these data for the proposed collision simulations, but that limitation is acknowledged and does not invalidate the initial-data construction itself. Therefore the appropriate verdict remains CONDITIONAL: the numerical claim requires added sensitivity and residual checks before the data can be accepted as physically reliable starting slices, but there is no demonstrated internal inconsistency in the construction.","tokens_in":15756,"tokens_out":42287,"duration_ms":444700,"concrete_test":"Recompute the Sec. V initial data for the same physical parameters as Fig. 13 (r0=R∞, M/R∞=1,2,3,4,5, z0=5R∞) with r1/R∞ in {2,3,4} and with xout=5,10,20, holding all other choices fixed. For each run, compute the ADM mass via Eq. (63), the horizon shape from Eq. (45), and the maximum Hamiltonian-constraint residual over the domain. If MADM or the horizon location shifts by more than the grid truncation error, the reported quantities are sensitive to the arbitrary shell; if the residual does not converge to zero, the numerical solution is not a valid solution of Eq. (9).","verdict_should_be":"UNCHANGED","load_bearing_attack":"In the numerical construction (Sec. V), the compactification profile is fixed by Eq. (51) with Phi=1 for r>=r1, and the transition radius r1/R∞=3 and outer boundary xout=10 are chosen by hand (Sec. V B). The function S(theta) is then determined by the no-conical-singularity condition Eq. (60), and Delta-Psi is solved from Eqs. (55)-(56) with the bubble condition (58) and the Robin condition (62). The paper never tests whether the reported physical outputs—Delta-Psi (Fig. 11), S(theta) (Fig. 12), the total ADM mass MADM (Fig. 13), and the apparent horizon surfaces (Fig. 14)—change when r1 or xout is varied. Because r1 is not a gauge degree of freedom but part of the chosen initial compactification profile, different r1 generically select different physical configurations; without a sensitivity study, the specific data cannot be claimed to represent the intended collision scenario rather than an artifact of the shell. Moreover, the only numerical validation is the grid-convergence plot of Delta-Psi (Fig. 10); no residual of the Hamiltonian constraint (9) is reported, so it is not demonstrated that the discrete solution satisfies the continuum constraint beyond the observed convergence of Delta-Psi on one parameter set.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs time-symmetric vacuum initial data for five-dimensional Kaluza-Klein spacetimes with a space-dependent compactification radius. Three classes are presented: analytic SO(3)-symmetric data (black string without a bubble, black string hiding a KK bubble, naked KK bubble), analytic multi-black-string data generalizing Brill-Lindquist, and numerical axisymmetric data for a black string and a KK bubble at separated locations. The constraint equation (4)R=0 is solved using conformally flat ansätze; the analytic solutions are explicit, and the numerical part solves an elliptic problem for ΔΨ and S(θ) with a prescribed transition profile for the compactification radius.","tokens_in":16105,"tokens_out":12538,"duration_ms":120396,"significance":"If the results hold, these are useful first-step initial data for numerical relativity studies of extra-dimensional dynamics, and the analytic families are interesting in their own right. Strengths include the absence of any fitting to target quantities, the explicit reduction to the Einstein-Rosen bridge and Brill-Lindquist data in the uniform-radius limits, and the closed-form ADM mass and horizon formulas. The numerical section demonstrates grid convergence for ΔΨ, but as detailed below, the physical robustness of the numerical data is not yet established. The analytic sections are largely checkable algebraically, though one formula in Sec. III B needs correction.","major_comments":[{"comment":"The numerical initial data are constructed with a hand-chosen transition radius r1/R∞=3 and outer boundary xout=10, and no sensitivity study is reported. Since the profile (51) sets Φ=1 for r≥r1 and varies in r0≤r<r1, the radius r1 is part of the physical compactification profile rather than a numerical gauge parameter; changing r1 generically changes the physical configuration. The reported ADM mass, S(θ), and apparent-horizon shapes (Figs. 11-14) could therefore depend strongly on this choice. The paper should either demonstrate that these quantities are insensitive to r1 and xout over a range, or clearly characterize the intended physical regime selected by the chosen values.","section":"V.B, Eq. (51), Figs. 11-14"},{"comment":"The only numerical accuracy check is the grid-convergence plot of ΔΨ in Fig. 10; no residual of the Hamiltonian constraint (9) is reported. Convergence of ΔΨ on one parameter set does not by itself establish that the discrete solution satisfies the continuum constraint, particularly near the bubble surface r=r0 where the equation is singular unless the boundary condition (58) is enforced. The authors should report the L2 or maximum residual of Eq. (9) on the numerical grid and check its convergence, including the behavior of S(θ).","section":"V.B, Eq. (9), Fig. 10"},{"comment":"The displayed expression for Δχ in Eq. (33) appears inverted. From the no-conical-singularity condition Δχ=2πΩ²(rB)/F'(rB) and Eqs. (29)-(30), one obtains Δχ=2π(rB²+M rB+rmin²)(rB+M/2)/rB², not the reciprocal printed in the paper. As written, Eq. (33) has dimensions of inverse length and contradicts the statement immediately below it that Δχ/2πM is always greater than 1/2. This affects Fig. 3 and the quantitative discussion of the bubble data; the formula and figure should be corrected.","section":"III.B, Eq. (33), Fig. 3"}],"minor_comments":[{"comment":"The text says there are five parameters, but then lists six: Δχ, r0, r1, rout, M, and z0; the count should be corrected.","section":"V.B, parameter list"},{"comment":"The convergence study in Fig. 10 reports errors only in ΔΨ; an analogous error measure for S(θ) would make the numerical convergence statement more complete.","section":"V.B, Fig. 10"},{"comment":"For the results in Figs. 11-14 the resolution is stated as (Imax,Jmax)=(200,45), but it would be helpful to also state how this resolution relates to the convergence study in Fig. 10 and whether the reported physical quantities are those of the converged solution.","section":"V.B, resolution"}],"recommendation":"major_revision","confidential_remarks":"The main concern is concentrated in Sec. V, where the numerical data are not yet shown to be robust with respect to the chosen transition radius and outer boundary, and no Hamiltonian-constraint residual is reported. The analytic sections are publishable after correcting Eq. (33). I do not see grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a solid, workmanlike construction of time-symmetric initial data for 5D Kaluza-Klein spacetimes with a position-dependent compactification radius. The genuinely new things are the analytic families: the SO(3) data with R0/R_infinity not equal to 1, the multi-black-string data built from the Phi = C + D/Psi ansatz, and the numerical setup that puts a black string and a KK bubble at different positions. The Phi = C + D/Psi trick that reduces the Hamiltonian constraint to a Laplace equation for Psi is clean and worth remembering, and the limits back to the Einstein-Rosen bridge and ordinary Brill-Lindquist data check out. I did not find a circular step; the free parameters are genuinely free inputs and the masses and areas are derived, not fitted.\n\nThe soft spot is exactly the one the stress-test flags: the numerical section of Sec. V depends on a hand-chosen transition radius r1 = 3 R_infinity and an outer boundary xout = 10, with no test of how the ADM mass, the bubble shape, or the apparent horizon respond to those choices. That is not a trivial gauge choice because different r1 correspond to different physical compactification profiles. The paper also never reports a residual of the Hamiltonian constraint, only grid convergence of Delta-Psi. So the two-object collision data are not yet established as robust. That said, the analytic sections are not affected, and the author explicitly acknowledges the bigger limitation that evolving a non-SO(3) bubble is still an open problem. The paper is honest about what it does and does not do.\n\nCitation pattern looks fine; the author cites the exact-solution literature and previous initial-data work, and the self-citations to [35,36] are only contextual, not load-bearing.\n\nWhat this paper is for: people who build and run numerical relativity simulations of higher-dimensional spacetimes, and anyone who wants explicit initial data with a varying extra dimension. It is a useful starting point, not a final answer. I would send it to peer review, but the referee should insist that the numerical section either add a sensitivity study for r1 and xout, report a constraint residual, or explicitly frame the data as illustrative rather than ready-to-evolve. No code is shipped, which is a minor annoyance but not a blocker for a methods paper in this subfield.","headline":"A genuinely new family of KK initial data with varying radius, clean analytic parts, but the numerical section needs a sensitivity study before the collision-simulation claim is sold.","tokens_in":16646,"tokens_out":1717,"would_cite":true,"duration_ms":20182,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C05","83E15"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper constructs time-symmetric initial data for five-dimensional Kaluza-Klein spacetimes in which the compactification radius varies in space, including a black string and a KK bubble at different positions, all satisfying the…","keywords":["Kaluza-Klein bubble","black string","initial data","Hamiltonian constraint","Brill-Lindquist initial data","compactification radius","numerical relativity","extra dimensions"],"falsifier":"Repeat the Section V construction with $r_1/R_\\infty = 2, 3, 4, 5$ and outer boundaries $x_{\\rm out} = 5, 10, 20, 40$ for the same physical parameters, and check whether $\\Delta\\Psi$, $M_{\\rm ADM}$, the bubble shape where $\\Phi=0$, and the apparent-horizon area converge to common values; if they drift, the data are slicing artifacts rather than genuine two-object initial data.","tokens_in":15443,"feed_emoji":"🫧","tokens_out":7022,"duration_ms":69889,"temperature":0.7,"pith_summary":"The paper's aim is to build starting slices for numerical relativity studies of a dynamical extra dimension. In time-symmetric vacuum data the momentum constraint is automatic and the only requirement is that the four-dimensional Ricci scalar of the slice vanish, ${}^{(4)}R=0$; the paper constructs slices of this kind in five-dimensional Kaluza-Klein spacetimes where the radius of the compact circle depends on spatial position. It delivers three families: analytic spherically symmetric data (a black string alone, a black string with a bubble hidden inside its horizon, and a naked bubble), analytic multi-black-string data that generalize Brill-Lindquist initial data, and numerical data for a black string and a KK bubble placed at different positions. These are intended as the initial conditions for simulating an expanding bubble of nothing meeting black objects.","feed_headline":"Initial data built for black string and KK bubble collisions","feed_subtitle":"Three valid starting data sets satisfy Einstein's equations with a space-varying extra dimension.","key_machinery":"The central objects are the conformal metric ansatz $ds^2 = \\Psi^4(dx^2+dy^2+dz^2)+\\Phi^2 d\\chi^2$ and the reduced Hamiltonian constraint $\\Psi\\nabla^2\\Phi+2\\nabla\\Phi\\cdot\\nabla\\Psi+4\\Phi\\nabla^2\\Psi=0$ for time-symmetric data, with the equivalent $F$-$\\Omega$ form used in the SO(3) section. The argument is carried by exact ansätze that turn this elliptic equation into solvable form: the closed forms $F=1+a_1/r+a_2/r^2$ and $\\Omega=1+b_1/r+b_2/r^2$ obeying $2(a_2+b_2)=a_1b_1$, and the reduction $\\Phi=C+D/\\Psi$ that collapses the constraint to Laplace's equation. The numerical construction uses a split at radius $r_1$ with $\\Phi=1$ outside and a smooth transition profile $\\Phi=1-\\exp\\left(\\frac{r-r_0}{r-r_1}S(\\theta)\\right)$ inside, with $S(\\theta)$ solved so that the $\\chi$-period is constant and the bubble is free of a conical singularity.","core_discovery":"The central claim is that such slices exist and can be written down or computed explicitly. With the metric ansatz $ds^2 = \\Psi^4(dx^2+dy^2+dz^2)+\\Phi^2 d\\chi^2$, the Hamiltonian constraint ${}^{(4)}R=0$ reduces to the elliptic equation $\\Psi\\nabla^2\\Phi+2\\nabla\\Phi\\cdot\\nabla\\Psi+4\\Phi\\nabla^2\\Psi=0$. The SO(3)-symmetric ansatz $F=1+a_1/r+a_2/r^2$, $\\Omega=1+b_1/r+b_2/r^2$ satisfies the analogous constraint under the single algebraic relation $2(a_2+b_2)=a_1b_1$, producing black-string and bubble data parameterized by mass, minimal-sphere radius, and either the central or bubble compactification radius. The multi-black-string data follow from the substitution $\\Phi=C+D/\\Psi$, which reduces the constraint to Laplace's equation $\\nabla^2\\Psi=0$, exactly as in the Brill-Lindquist case. For a black string and a bubble at separate locations, the paper solves the constraint numerically by splitting the domain at a radius $r_1$, taking $\\Phi=1$ outside and a smooth transition profile inside with an auxiliary function $S(\\theta)$ that enforces a constant $\\chi$-period, and reports convergence of the solution with grid resolution.","pith_inferences":["The paper fixes the junction radius $r_1/R_\\infty=3$ and outer boundary $x_{\\rm out}=10$ without varying them; a natural next check is whether the ADM mass, bubble shape, and horizon geometry are stable under changes of these numerical boundaries.","The same split-domain construction should extend to non-axisymmetric or unequal-mass black-string/bubble systems, and to non-time-symmetric slices with nonzero extrinsic curvature, where the momentum constraint would also have to be solved.","If the naked-bubble slices with negative ADM mass are evolved, they offer a concrete arena to test cosmic censorship in five dimensions, since the paper's horizon-area comparison suggests a naked singularity may form.","A useful by-product of the analytic families is that they provide closed-form calibration cases for code tests of constraint solvers with position-dependent compactification radius."],"forward_implications":["The SO(3) data without a bubble give a two-parameter family of black-string slices whose ADM mass is always larger than $M/2$, and whose apparent-horizon area never exceeds that of the equal-mass Schwarzschild string.","The bubble data include naked bubbles with negative ADM mass; when the bubble is trapped, the horizon area can exceed the Schwarzschild-string value, a regime the paper connects to black-string instability and possible naked-singularity formation under evolution.","The Brill-Lindquist-type data show that the condition for a common horizon around two equal black strings depends strongly on the ratio $R_p/R_\\infty$ of the compactification radius at the puncture to that at infinity, with the largest critical separation (measured in ADM mass) at $R_p/R_\\infty=1$.","The numerical data for a separated black string and KK bubble satisfy the constraint to the achieved convergence, yield a common apparent horizon around the black string, and give total ADM mass $M_{\\rm ADM}=M+\\Delta M$ with $\\Delta M<0$.","Together, these slices are the input layer for evolving the nonlinear dynamics of a varying extra dimension, including bubble-black-string collisions."],"supporting_citations":[{"why":"Supplies the bubble-of-nothing solution that motivates the physical setup of these initial data.","marker":"[9]"},{"why":"Shows that KK bubbles can carry negative gravitational mass, a feature the paper's bubble data reproduce.","marker":"[10]"},{"why":"Provides the Brill-Lindquist multi-black-hole initial data that the multi-black-string construction generalizes.","marker":"[32]"},{"why":"Develops numerical treatments of KK bubble surfaces that the paper's bubble regularity and horizon analysis build on.","marker":"[13,14]"},{"why":"Documents black-string instability and naked-singularity formation, referenced when horizon area exceeds the Schwarzschild-string value.","marker":"[33,34]"},{"why":"Provides the higher-dimensional BSSN and Cartoon formalism that the new initial data are designed to feed into evolution codes.","marker":"[35,36]"}],"fun_headline_variants":["Space-varying extra dimension meets black strings in new data","Black strings and bubbles: initial data with varying compactification","First initial data for dynamical extra dimensions with black objects","Nonuniform extra dimension: analytic and numeric initial data","Ready for simulation: black string and KK bubble initial data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The numerical construction for a black string and a bubble at separate locations fixes the junction radius $r_1$ and outer boundary $r_{\\rm out}$ by hand, and the paper does not test whether physical quantities such as the ADM mass, bubble shape, or apparent horizon depend on these choices.","fun_headline_variants_meta":{"raw":{"variants":["Space-varying extra dimension meets black strings in new data","Black strings and bubbles: initial data with varying compactification","First initial data for dynamical extra dimensions with black objects","Nonuniform extra dimension: analytic and numeric initial data","Ready for simulation: black string and KK bubble initial data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000717,"raw_usage":{"total_tokens":3258,"prompt_tokens":1016,"completion_tokens":2242,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":632,"completion_tokens_details":{"reasoning_tokens":2164}},"tokens_in":632,"tokens_out":2242,"duration_ms":16838,"temperature":1.0,"reasoning_tokens":2164,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T18:03:07.711697+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the Section V construction with $r_1/R_\\infty = 2, 3, 4, 5$ and outer boundaries $x_{\\rm out} = 5, 10, 20, 40$ for the same physical parameters, and check whether $\\Delta\\Psi$, $M_{\\rm ADM}$, the bubble shape where $\\Phi=0$, and the apparent-horizon area converge to common values; if they drift, the data are slicing artifacts rather than genuine two-object initial data.","supporting_citations":[{"cited_title":"Witten, Nucl","cited_arxiv_id":null,"evidence_quote":"Supplies the bubble-of-nothing solution that motivates the physical setup of these initial data."},{"cited_title":"Brill and H","cited_arxiv_id":null,"evidence_quote":"Shows that KK bubbles can carry negative gravitational mass, a feature the paper's bubble data reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Brill-Lindquist multi-black-hole initial data that the multi-black-string construction generalizes."}],"review_version":1}