pith. sign in

arxiv: 2604.21325 · v1 · submitted 2026-04-23 · 🌀 gr-qc · math-ph· math.MP

Gravitational Collapse of a Chiellini Integrable Scalar Field

Mohamed Aarif A , Soumya Chakrabarti This is my paper

Pith reviewed 2026-05-09 21:30 UTC · model grok-4.3

classification 🌀 gr-qc math-phmath.MP
keywords gravitational collapsescalar fieldChiellini integrableasymptotic collapseMilne-Pinney equationenergy conditionsapparent horizonVaidya exterior
0
0 comments X

The pith

A spatially homogeneous scalar field with an extended Higgs potential undergoes asymptotic gravitational collapse without reaching zero volume at finite time.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper studies the gravitational collapse of a non-interacting mixture of perfect fluid and a homogeneous scalar field. By adopting an extended Higgs-type potential that permits reduction of the Klein-Gordon equation to a generalized damped Milne-Pinney equation, the authors obtain closed-form solutions for the scale factor and scalar field. These solutions show that the proper volume decreases monotonically but approaches zero only as time tends to infinity. The analysis also covers energy conditions, possible multiple apparent horizons, and smooth matching of the interior to a generalized Vaidya exterior.

Core claim

We find that it exhibits an asymptotic collapse in which the proper volume decreases monotonically but never reaches zero at finite time. Closed-form analytical solutions are derived for the scalar field and scale factor in the collapsing branch, with the scalar field remaining canonical while the perfect fluid can violate the null energy condition, and the apparent horizon condition yielding no trapped surface or multiple horizons depending on parameters.

What carries the argument

Reduction of the Klein-Gordon equation to a generalized damped Milne-Pinney class of differential equation via the Chiellini-integrable framework with an extended Higgs-type self-interaction potential.

If this is right

  • The proper volume decreases monotonically toward zero only as time tends to infinity.
  • The perfect fluid component can violate the null energy condition while the scalar field remains canonical.
  • The apparent horizon condition admits either no trapped surface or multiple apparent horizons depending on the parameter values.
  • The interior homogeneous solution matches smoothly to a generalized Vaidya exterior metric through the Israel-Darmois junction conditions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This class of integrable models may provide explicit examples of collapse that avoids finite-time singularities.
  • The parameter-dependent horizon behavior could be checked against numerical relativity simulations of similar setups.
  • The smooth exterior matching offers a concrete starting point for studying the emitted gravitational radiation or the final state as time goes to infinity.

Load-bearing premise

The scalar field is spatially homogeneous and the self-interaction potential is chosen as an extended Higgs-type form that permits reduction to a Chiellini-integrable damped Milne-Pinney equation.

What would settle it

A direct numerical solution of the Einstein-scalar field equations for the same potential without imposing the Chiellini integrability condition, showing whether the volume reaches zero in finite time.

Figures

Figures reproduced from arXiv: 2604.21325 by Mohamed Aarif A, Soumya Chakrabarti.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: The Apparent Horizon Condition as a function of [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: The Apparent Horizon Condition as a function of [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
read the original abstract

We study the gravitational collapse of a non-interacting mix of perfect fluid and a spatially homogeneous scalar field within a Chiellini-integrable framework. We choose an extended Higgs-type self-interaction potential and reduce the Klein-Gordon equation into a generalized damped Milne-Pinney class of differential equation. We derive a closed-form analytical solution for the scalar field, the scale factor and explore the collapsing branch of the same. We find that it exhibits an asymptotic collapse in which the proper volume decreases monotonically but never reaches zero at finite time. We analyze the energy conditions for the constituent elements of the collapsing sphere. While the scalar field remains canonical in nature, we find that the perfect fluid can violated the Null Energy Condition. We also study the formation of apparent horizon condition and find multiple possibilities depending on the parameter space : either no trapped surface or the formation of multiple apparent horizons. We match the interior homogeneous solution to a generalized Vaidya exterior via the Israel-Darmois junction conditions, yielding the corresponding boundary mass function, ensuring a smooth collapse scenario.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The paper claims to study the gravitational collapse of a non-interacting mixture of perfect fluid and a spatially homogeneous scalar field in a Chiellini-integrable framework. By choosing an extended Higgs-type self-interaction potential, the Klein-Gordon equation is reduced to a generalized damped Milne-Pinney differential equation, for which a closed-form analytical solution is derived for the scalar field and the scale factor. The collapsing branch exhibits asymptotic collapse where the proper volume decreases monotonically but never reaches zero at finite time. Energy conditions are analyzed, with the scalar field remaining canonical and the perfect fluid capable of violating the null energy condition. The formation of apparent horizons is studied, showing possibilities of no trapped surface or multiple apparent horizons depending on the parameter space. The interior solution is matched to a generalized Vaidya exterior using Israel-Darmois junction conditions, providing the boundary mass function for a smooth collapse.

Significance. If the closed-form solution and its properties hold, this work offers a valuable exact model for asymptotic gravitational collapse in the presence of scalar fields, contributing to the understanding of singularity avoidance and horizon formation in general relativity. The strengths include the derivation of an exact integrable solution and the explicit matching to an exterior spacetime, which allows for concrete predictions about energy conditions and trapped surfaces. This could be useful for testing numerical simulations or exploring cosmic censorship in scalar field models.

major comments (2)
  1. [Apparent horizon analysis] In the apparent horizon analysis, the statement that there are 'multiple possibilities' for apparent horizon formation depending on the parameter space is made without providing the specific conditions or ranges of parameters that lead to no trapped surface versus multiple horizons; this makes the claim difficult to verify and is central to the horizon matching discussion.
  2. [Energy conditions section] While it is stated that the perfect fluid can violate the Null Energy Condition, the explicit expressions for the energy density and pressure of the fluid (or their ratio) are not given in terms of the derived scale factor, preventing assessment of the extent and consistency of the violation during collapse.
minor comments (3)
  1. The abstract and introduction should include a brief definition or reference for 'Chiellini integrability' to make the paper more accessible.
  2. Some figures (if present) showing the scale factor evolution would benefit from labels indicating the specific parameter values used for the plotted curves.
  3. [Matching section] The Israel-Darmois conditions are invoked, but a short recap of the continuity requirements for the metric and extrinsic curvature would improve clarity for readers.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We agree that the two major points raised require clarification and will revise the manuscript accordingly to improve verifiability of the results.

read point-by-point responses

  1. Referee: In the apparent horizon analysis, the statement that there are 'multiple possibilities' for apparent horizon formation depending on the parameter space is made without providing the specific conditions or ranges of parameters that lead to no trapped surface versus multiple horizons; this makes the claim difficult to verify and is central to the horizon matching discussion.

    Authors: We agree that the specific conditions on the parameter space were not stated explicitly. The apparent horizon condition is given by the vanishing of the expansion of outgoing null geodesics, which reduces to an algebraic equation involving the scale factor a(t) and its derivative at the boundary. Since we have a closed-form expression for a(t) from the integrable Milne-Pinney equation, the number of real positive roots (corresponding to horizons) can be analyzed via the discriminant or by examining the monotonicity of the relevant function. In the revised manuscript we will derive and tabulate the explicit inequalities on the integration constants and potential parameters that separate the regimes of no trapped surface from those permitting one or multiple apparent horizons. This will be added as a new subsection or appendix to make the claims directly verifiable. revision: yes

  2. Referee: While it is stated that the perfect fluid can violate the Null Energy Condition, the explicit expressions for the energy density and pressure of the fluid (or their ratio) are not given in terms of the derived scale factor, preventing assessment of the extent and consistency of the violation during collapse.

    Authors: We acknowledge the omission. The energy density and pressure of the perfect fluid follow directly from the Friedmann and acceleration equations once the total energy density and pressure (scalar field plus fluid) are known. With the closed-form scale factor a(t) and scalar field φ(t) in hand, we can isolate ρ_f(t) = 3H² - ρ_φ and p_f(t) = -2Ḣ - 3H² - p_φ. In the revised version we will write these explicit expressions in terms of a(t) and its derivatives, then evaluate the combination ρ_f + p_f throughout the collapse to identify the intervals where the null energy condition is violated and to quantify the duration and strength of the violation for representative parameter choices. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained

full rationale

The paper imposes spatial homogeneity and selects an extended Higgs-type potential specifically to reduce the Klein-Gordon equation to a Chiellini-integrable damped Milne-Pinney equation, then obtains an exact closed-form solution for the scalar field and scale factor. The central result—an asymptotic collapse in which proper volume decreases monotonically but never reaches zero at finite time—follows from direct analysis of this exact solution and is not equivalent to the input assumptions by construction. No parameters are fitted to data and then relabeled as predictions; no self-citation chain is invoked to justify uniqueness or load-bear the collapse behavior; energy conditions and apparent-horizon formation are likewise computed from the derived solution. The setup satisfies the Einstein equations within the chosen ansatz, and the integrability condition is an external mathematical property of the chosen potential rather than a tautological redefinition of the output. The derivation therefore remains independent of its inputs.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The model rests on the assumption of spatial homogeneity for the scalar field and the selection of an extended Higgs potential that forces the equation into an integrable form; no new particles or forces are introduced.

free parameters (1)
  • potential parameters
    Parameters in the extended Higgs-type potential are chosen by hand to satisfy the Chiellini integrability condition.
axioms (2)
  • domain assumption Spatially homogeneous scalar field
    Invoked to simplify the Klein-Gordon equation to an ordinary differential equation.
  • domain assumption Non-interacting perfect fluid and scalar field
    Stated in the abstract as the physical setup.

pith-pipeline@v0.9.0 · 5481 in / 1376 out tokens · 40285 ms · 2026-05-09T21:30:15.231112+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

53 extracted references · 53 canonical work pages

  1. [1]

    A similar reduction can be found in the theory of solitons and non-linear wave equation

    showed that a hybrid Ermakov-Painleve system may arise from a reduction of a coupled nonlinear Schrodinger equation system. A similar reduction can be found in the theory of solitons and non-linear wave equation. In extension to these works, there is a re- newed interest in Chiellini-type integrability conditions [17]. These conditions establish a specifi...

    work page 2023

  2. [2]

    B. Datt, Z. Phys.108, 314 (1938)

    work page 1938

  3. [3]

    J. R. Oppenheimer and H. S. Snyder, Phys. Rev.56, 455 (1939)

    work page 1939

  4. [4]

    Penrose, Phys

    R. Penrose, Phys. Rev. Lett.14, 57 (1965)

    work page 1965

  5. [5]

    Penrose, Nuovo Cimento Rivista Serie1(1969)

    R. Penrose, Nuovo Cimento Rivista Serie1(1969)

    work page 1969

  6. [6]

    S. M. C. V. Goncalves and I. G. Moss, Class. Quant. Grav.14, 2607 (1997)

    work page 1997

  7. [7]

    M. W. Choptuik, Phys. Rev. Lett.70, 9 (1993)

    work page 1993

  8. [8]

    P. R. Brady, C. M. Chambers, and S. M. C. V. Goncalves, Phys. Rev. D56, R6057 (1997)

    work page 1997

  9. [9]

    Gundlach, Phys

    C. Gundlach, Phys. Rev. Lett.75, 3214 (1995) ; Liv. Rev. Rel.2, 1 (1999)

    work page 1995

  10. [10]

    Christodoulou, Comm

    D. Christodoulou, Comm. Math. Phys.109, 591 (1987) ; Comm. Math. Phys.109, 613 (1987) ; Ann. Math.140, 607 (1994)

    work page 1987

  11. [11]

    D. S. Goldwirth and T. Piran, Phys. Rev. D36, 3575 (1987)

    work page 1987

  12. [12]

    W. E. Milne, The numerical determination of character- istic numbers, Phys. Rev. 35 (1930) 863-867

    work page 1930

  13. [13]

    Pinney, Proc

    E. Pinney, Proc. A. M. S. 1 (1950) 681

    work page 1950

  14. [14]

    V. P. Ermakov, Second-order differential equations. Con- ditions of complete integra- bility, Univ. Isz. Kiev Series III 9, 1-25 (1880) ; Appl. Anal. Discrete Math. 2 (2008) 123-148

    work page 2008

  15. [15]

    Haas, Approximate solution for a damped Pinney equation, math-ph/0712.4083

    F. Haas, Approximate solution for a damped Pinney equation, math-ph/0712.4083

    work page arXiv

  16. [16]

    Cari˜ nena and J

    J.F. Cari˜ nena and J. de Lucas J, Int. J. Geom. Methods Mod. Phys. 6 (2009) 683-99

    work page 2009

  17. [17]

    Rogers, A novel Ermakov-Pinalev´ eII system, N+ 1-dimensional coupled NLS and elastodynamic reductions, Stud

    C. Rogers, A novel Ermakov-Pinalev´ eII system, N+ 1-dimensional coupled NLS and elastodynamic reductions, Stud. Appl. Math.133(2014) 214-231

    work page 2014

  18. [18]

    A.Chiellini, Sull’integrazione dell’equazione differeziale ˙y+P y 2 +Qy 3 = 0, Bollettino dell’Unione Matematica Italiana, 10, 301-307 (1931)

    work page 1931

  19. [19]

    Lienard,Revue generale de l’electricite, 23, 901-912 and 946-954(1928)

    A. Lienard,Revue generale de l’electricite, 23, 901-912 and 946-954(1928)

    work page 1928

  20. [20]

    Ivan Bandic, Bollettino dell’Unione Matematica Italiana, 16, 59-67 (1961)

    work page 1961

  21. [21]

    Ivan Bandic, Bollettino dell’Unione Matematica Italiana, 17, 81-91 (1962)

    work page 1962

  22. [22]

    M.K.Mak, T.Harko, Comput.Math.Appl.43,91-94(2002)

    work page 2002

  23. [23]

    M.K.Mak, H.W.Chan, T.Harko, Comput.Math.Appl.41,1395-1401(2001)

    work page 2001

  24. [24]

    T.Harko,M.K.Mak, Comput.Math.Appl.46,849- 853(2003)

    work page 2003

  25. [25]

    A.V.Yurov, V.A.Yurov, J.Math.Phys.51,082503(2010)

    work page 2010

  26. [26]

    T.Harko,F.S.N.Lobo, M.K.Mak,arXiv:1302.0836v2

    work page arXiv

  27. [27]

    Mancas, and P

    H.C.Rosu, S. Mancas, and P. Chen, Phys Lett. A 379 (2015) 882-887

    work page 2015

  28. [28]

    Mancas, H.C

    S.C. Mancas, H.C. Rosu, Phys. Lett. A, 377 (2013), 1434

    work page 2013

  29. [29]

    Mancas, H.C

    S.C. Mancas, H.C. Rosu, Integrable equations with Ermakov–Pinney nonlinearities and Chiellini damping, Appl. Math. Comp., 259 (2015) 1-11

    work page 2015

  30. [30]

    Banerjee and S

    N. Banerjee and S. Chakrabarti, Phys. Rev. D.95, 024015 (2017)

    work page 2017

  31. [31]

    Polyanin and V.F.Zaitsev,Handbook of Exact So- lutions for Ordinary Differential Equations, Chapman & Hall/CRC, 2003

    A.D. Polyanin and V.F.Zaitsev,Handbook of Exact So- lutions for Ordinary Differential Equations, Chapman & Hall/CRC, 2003

    work page 2003

  32. [32]

    R. A. Jimenez and J. Sola, Phys. Lett. B,389(1): 53 (1996)

    work page 1996

  33. [33]

    H. E. Camblong, L. N. Epele, H. Fanchiotti and C. A. Garcia Canal, Phys. Rev. Lett.,85: 1590 (2000)

    work page 2000

  34. [34]

    Lopez-Val and J

    D. Lopez-Val and J. Sola, Phys. Rev. D,81: 033003, (2010)

    work page 2010

  35. [35]

    Hinterbichler, J

    K. Hinterbichler, J. Khoury, A. Levy and A. Matas, Phys. Rev. D,84(10)(2011)

    work page 2011

  36. [36]

    J. Sola, E. Karimkhani and A Khodam-Mohammadi, Class. Quant. Gravit.,34(2): 025006 (2016)

    work page 2016

  37. [37]

    R-G Cai and S-J Wang, Phys. Rev. D,103: 023502 (2021)

    work page 2021

  38. [38]

    Chakrabarti, K

    S. Chakrabarti, K. Dutta and J. L. Said, Mon. Not. Roy. Astron. Soc.,514(1): 427 (2022)

    work page 2022

  39. [39]

    P. G. C. Rufo, A. Kumar, C. Sabin and A. Mazumdar, Phys. Rev. D111: 115005 (2025)

    work page 2025

  40. [40]

    G. D. Moore, Phys. Rev. D,93: 065043, (2016)

    work page 2016

  41. [41]

    Barrabes and G

    C. Barrabes and G. F. Bressange, Class. Quant. Grav. 14: 805 (1997)

    work page 1997

  42. [42]

    Viqar Husain, Phys. Rev. D53, 4 (1996)

    work page 1996

  43. [43]

    Anzhong Wang, Wu Yumei, Gen. Relativ. Grav.31, 1 (1999)

    work page 1999

  44. [44]

    P. C. Vaidya, Proc. Indian Acad. Sci. A33, 264 (1951)

    work page 1951

  45. [45]

    P. C. Vaidya, Phys. Rev.83, 1 (1951)

    work page 1951

  46. [46]

    P. C. Vaidya and I.M. Pandya, Prog. Theor. Phys.35, 1 (1966)

    work page 1966

  47. [47]

    P. C. Vaidya, Phys. Rev.174, 5 (1968)

    work page 1968

  48. [48]

    Vaidya, Gen

    P.C. Vaidya, Gen. Relativ. Grav.31, 1 (1999) 8

    work page 1999

  49. [49]

    Deruelle, M

    N. Deruelle, M. Sasaki, and Y. Sendouda, Prog. Theor. Phys.119, 237 (2008)

    work page 2008

  50. [50]

    J. M. M. Senovilla, Phys. Rev. D.88, 064015 (2013)

    work page 2013

  51. [51]

    N. O. Santos, Mon. Not. R. Astron. Soc.216, 403 (1985)

    work page 1985

  52. [52]

    Chan, Mon

    R. Chan, Mon. Not. R. Astron. Soc.316, 588 (2000)

    work page 2000

  53. [53]

    C. W. Misner and D. Sharp, Phys. Rev.136, B571 (1964)

    work page 1964