Pith. sign in

REVIEW 3 major objections 6 minor 1 cited by

Casimir effect with an unstable mode

T0 review · 3 major / 6 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read In a scalar field with an unstable mode, the condensate itself contributes a repulsive Casimir force.

desk verdict A genuine first treatment of the unstable-mode Casimir effect with a condensate; the central stability claim is asserted rather than proved, so it deserves refereeing with a request for the missing proof. read the letter →

arxiv 2502.09407 v1 pith:ZBCYNHAU submitted 2025-02-13 quant-ph

M. Bordag , I.G. Pirozhenko This is my paper · ORCID
classification quant-ph MSC 81T5535Q55 PACS 03.70.+k03.75.Lm
keywords CasimireffectunstablemodecondensateGross-PitaevskiiequationellipticJacobifunctionsvacuumenergyRobinboundaryconditionzetaregularization
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper argues that the Casimir effect has a second source when the vacuum is unstable: a static condensate that forms through the nonlinear Gross-Pitaevskii equation. The standard vacuum-fluctuation force is joined by a force coming from the condensate energy, and in the models studied the condensate force is always repulsive. In a (1+1)-dimensional scalar field with self-interaction, the paper computes both contributions for a finite interval with Robin and Dirichlet boundaries and shows they compete. If the claim holds, Casimir-force calculations for systems with critical modes cannot stop at the fluctuation spectrum; the condensate contributes on the same footing.

What carries the argument

The load-bearing object is the condensate solution $\phi_0(x)$ of the Gross-Pitaevskii equation $(-\partial_x^2 + m^2 + V(x) + \lambda \phi_0^2)\, \phi_0 = 0$, obtained exactly in terms of elliptic Jacobi functions for a delta potential, a Robin-Dirichlet interval, and a potential hole. After the field shift $\phi = \phi_0 + \text{fluctuations}$, the fluctuation operator is $-\partial_x^2 + m^2 + V(x) + 3\lambda \phi_0^2$; the paper asserts it is stable by construction for an exact GP solution. The vacuum energy is computed from the mode-generating function, built as a Wronskian, through a contour representation regularized by a zeta function, with heat-kernel coefficients subtracted to renormalize. The condensate energy simplifies to $-\tfrac{\lambda}{4}\int \phi_0^4\, dx$, and its derivative with respect to the interval length gives the condensate force, which the paper finds to be always repulsive.

What would settle it

Compute the lowest eigenvalue of the fluctuation operator $-\partial_x^2 + m^2 + V(x) + 3\lambda \phi_0^2$ for one of the exact elliptic solutions (for example the Robin-Dirichlet interval with $m=1$, $\lambda=1$, $\kappa=2$, and $L$ in the region $L_0 < L < L_1$) by direct numerical solution of the boundary-value problem, and look for a bound state with energy below $-m^2$. If any such state exists, the asserted stability by construction fails and the vacuum-energy integral misses a contribution. Alternatively, check whether the mode-generating function $\Phi(i\xi)$ has a zero for real $\xi > m$ in the critical case.

Watch

Extended reading notes

Core claim

The paper's central claim is that in a (1+1)-dimensional real scalar field theory whose potential $V(x)$ is strong enough to create an imaginary-frequency bound state, the correct ground-state description requires splitting the field into a static condensate $\phi_0(x)$ and quantum fluctuations. The condensate solves the real Gross-Pitaevskii equation, and after the shift the fluctuation operator has only real eigenfrequencies, so the vacuum energy can be computed by standard zeta-function methods. For the finite-interval example with a Robin condition at one end and Dirichlet at the other, the total Casimir force splits into a fluctuation part and a condensate part. The fluctuation part is repulsive in this asymmetric setup, and the condensate part is repulsive in all cases considered, so the two can either reinforce or compete. The paper also shows that a mean-field approximation based on the linear bound state reproduces the exact condensate energy near the threshold of instability, with a small error, and identifies a parameter region where the approximate method gives a false positive.

Load-bearing premise

The entire calculation rests on the assertion that after shifting by a Gross-Pitaevskii solution, the fluctuation operator has no negative or imaginary eigenvalues; the paper guarantees this 'by construction' for exact solutions but does not prove it, and for approximate solutions it is explicitly not guaranteed.

Editorial extensions

If this is right

  • In any system with an unstable mode, the total Casimir force is the sum of a condensate force and a vacuum-fluctuation force, and both must be included when comparing with experiment.
  • The condensate force is always repulsive in the models treated, so a configuration whose vacuum force is attractive (for example equal Robin conditions on both sides) could have a total force that changes sign as the parameters vary.
  • Near the threshold of instability, the linear bound-state wave function is a quantitatively reliable approximation to the exact condensate for computing the condensate energy.
  • The k-gap shows that the existence of a condensate solution can depend sensitively on interval length and coupling, and a mean-field approximation can falsely predict a solution there.
  • The vacuum energy is much more sensitive to boundary conditions than to the shape of the condensate potential inside the interval.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the repulsive condensate force persists in higher-dimensional geometries, the usual expectation that Casimir forces between identical plates are attractive could be reversed in systems with critical modes, a testable prediction for condensed-matter or superconducting analogues.
  • The near-independence of the vacuum energy from the condensate profile suggests that measurements of the Casimir force in such systems would constrain boundary parameters much more strongly than the self-interaction strength $\lambda$.
  • The k-gap structure may have an analogue in the allowed parameter space of trapped Bose-Einstein condensates, where interval length plays the role of a tunable cavity size.
  • Generalizing the calculation to a complex field or to higher dimensions would clarify whether 'always repulsive' is a structural feature of the GP condensate contribution or an artifact of the one-dimensional elliptic solutions.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper studies the Casimir effect in a (1+1)-dimensional real scalar field theory with a quartic self-interaction, in the presence of background potentials that support an unstable (imaginary-frequency) mode. The authors shift the field by a static condensate satisfying the Gross-Pitaevskii (GP) equation and then identify two contributions to the Casimir energy: the classical energy of the condensate, Eq. (19), and the vacuum energy of the quantum fluctuations, computed via the zeta-function method, Eq. (96). Three solvable models are presented: a delta-function potential, a finite interval with Robin and Dirichlet boundary conditions, and a potential hole. Exact GP solutions are given in terms of elliptic Jacobi functions, and a bound-state approximation scheme is developed and tested near the critical threshold. For the Robin–Dirichlet cavity the authors compute both contributions to the force and find a repulsive condensate force competing with a repulsive vacuum-fluctuation force.

Significance. If the central result holds, the paper extends the Casimir effect to a new class of systems—those with an unstable mode—and demonstrates a mechanism by which a condensate contributes a repulsive force that can compete with or compensate the vacuum-fluctuation force. The work is largely self-contained: the exact GP solutions are checked against the nonlinear equations, the bound-state approximation is compared with exact results within the paper, and no external parameters are fitted. The numerical evaluation of the vacuum energy via a mode-generating function and the subtraction of heat-kernel asymptotics follows standard methods and is clearly presented. The main weakness is that the spectral stability of the fluctuation operator after the GP shift is asserted rather than proved; this is a load-bearing issue because the vacuum-energy calculation is valid only for a positive spectrum. The paper also states a global claim about the repulsiveness of the condensate force on the basis of only a few examples.

major comments (3)
  1. The claim that the fluctuation operator in Eq. (82) 'has no imaginary eigenvalues by construction' is not proved. The associated quadratic form is Q(φ) = ∫_0^L (φ')² dx + ∫_0^L V0(x) φ² dx - κ φ(0)², where V0 = 3λϕ0² and the boundary term is negative. Solving the GP equation (18) is a stationarity condition on the energy functional, not a sufficient condition for the second variation to be nonnegative. If a negative eigenvalue were present, the contour representation (68) and the renormalized vacuum energy (96) would require substantial revision, and the central sum-of-forces result would be unreliable. The authors should provide a proof of spectral positivity for the exact GP solution, or at least a systematic numerical search for negative eigenvalues over the parameter range used in Figs. 8 and 9.
  2. The statement that the condensate force 'is always repulsive' is a global claim, but it is supported only by the numerical examples in Figs. 6 and 8. No analytic argument or comprehensive parameter scan is given. Either prove this monotonicity property from the expression for Econd (e.g., by showing ∂Econd/∂L ≤ 0 for all admissible parameters) or rephrase the claim as 'in the examples considered'.
  3. In the approximate-solution case of Section 6.3, the potential (98) is used to compute the vacuum energy. The paper itself admits in Section 7 that the approximate solution 'has no guarantee' of stability and notes only that no instability was observed numerically. Since the vacuum-energy calculation depends on the positivity of the fluctuation spectrum, a numerical check of the lowest eigenvalue of Eq. (82) with the approximate potential should be reported, especially because the exact and approximate potentials differ notably (Fig. 9).
minor comments (6)
  1. The phrase 'the module |x|' should read 'the modulus |x|'.
  2. The text 'on the function Φ 0(x)' should read 'on the function ϕ0(x)'.
  3. The word 'accsin' is a typo and should read 'arcsin'.
  4. The derivation of the Wronskian equality W = ξ u(L) is compressed to 'from the x-independence of W and the boundary conditions'; a one-line explanation would improve readability.
  5. The legends in the upper panels overlap with the curves, and the different curves are not distinguished by line styles. Please use clear labels or distinct markers to improve readability.
  6. Ref. [15] is cited for the Casimir effect in λϕ⁴ theory, but the relation of that earlier work to the present condensate-plus-fluctuations decomposition could be stated more explicitly in the introduction.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the condensate stability assertion follows from the GP solution rather than being used as an input, and all computed energies/forces are model-derived.

full rationale

The derivation chain is self-contained. The only passage that could look like a definitional shortcut is Section 6.2: 'By construction of the Lagrangian 9, the equation 82 has no imaginary eigenvalues.' This is an assertion, not a circular reduction: the exact condensate φ0 is obtained by solving the Gross-Pitaevskii equation and the matching condition (Eq. 53), and the fluctuation operator in Eq. 82 is then a Schrödinger operator with potential 3λφ0². For a positive GP solution satisfying the same Robin/Dirichlet boundary conditions, the absence of negative eigenvalues is a theorem (by comparing with the positive solution via a Wronskian argument), not the definition of the prediction. The paper suppresses the proof, but that is a rigor/correctness issue, not circularity. The condensate energy Eq. 19 follows algebraically from the GP equation; the exact elliptic-function solutions are checked against the field equation and boundary conditions; the perturbative parameter µ is fixed by minimizing the bound-state energy (Eq. 29), not by the energies later compared; and the condensate force is computed as -∂_L E_cond. The approximate solution is explicitly benchmarked against exact solutions in the same paper (e.g., Eq. 50 vs. Eq. 43 and Figures 3, 5, 7), so no fitted parameter is renamed as a prediction. Self-citations (Refs. 1, 15, 16, 17) supply standard methods and background context; none is invoked as a uniqueness theorem or as the sole justification of the central claim. Therefore no step reduces by construction to its own input, and the paper has no significant circularity.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

No free parameters are fitted to data: the couplings m, λ, the boundary strength κ, and the potential parameters L, U0, R are inputs, while the elliptic module k and shift x1 are fixed by boundary conditions. The main unproved input is the stability of the fluctuation spectrum after the GP shift, and the use of the classical GP equation as the condensate description is standard but not derived from the quantum theory. No new particles or forces are introduced.

assumptions (5)
  • domain assumption The condensate is correctly described by the real Gross-Pitaevskii equation (12), derived from the shifted Lagrangian (9) with λ > 0.
    Standard mean-field description; higher-order terms in (9) are dropped without a quantitative estimate of their effect.
  • domain assumption After shifting to a GP solution, the fluctuation operator -∂² + V + 3λφ0² has no imaginary eigenvalues.
    Invoked in Sections 2 and 6.2 as 'stable by construction'; not proven, only numerically checked in the examples, and explicitly not guaranteed for approximate solutions.
  • standard math Zeta-function regularization with heat-kernel subtraction (72) is the correct renormalization prescription for the vacuum energy.
    Standard method from Refs. [1,19]; the paper applies it without re-deriving the subtraction.
  • standard math The mode generating function is represented by a Wronskian, and its large-ξ asymptotic expansion determines the divergent part of the vacuum energy.
    Standard contour-integral technique, Eqs. (66)-(71), based on Ref. [1].
  • standard math Jacobi elliptic function identities from the NIST handbook [18] are used to construct exact GP solutions.
    The delta-function solution (41) can be verified directly; the Robin and hole solutions rely on machine algebra and are stated to satisfy the GP equation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Casimir effect with an unstable mode." pith.science (2026). https://pith.science/paper/ZBCYNHAU

@misc{pith2026250209407,
  author       = {Pith},
  title        = {Pith review of: Casimir effect with an unstable mode},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZBCYNHAU}},
  note         = {Machine review of arXiv:2502.09407}
}
read the original abstract

We consider the Casimir effect in a (1+1)-dimensional model with a critical mode. Such a mode gives rise to a condensate described by the nonlinear Gross-Pitaevskii equation. In the condensate, there are two sources of the Casimir force; one is the conventional one resulting from the fluctuations, the other follows from the condensate. We consider three simple models that allow for condensate solutions in terms of elliptic Jacobi functions. We also investigate a method for obtaining approximate solutions and show its range of applicability. In all three examples we compute the condensate energy. In one example with a finite interval with Robin boundary conditions on one side and Dirichlet conditions on the other side, we calculate the vacuum energy and the Casimir force. There is a competition between the forces from the condensate and the fluctuations. We mention that the force from the condensate is always repulsive.

Figures

Figures reproduced from arXiv: 2502.09407 by the authors.

Figure 1
Figure 1. The energy 36 of the bound state solution and the energy 43 of the exact solution of [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. The left panel shows the left-hand side of equation 53 as a function of [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. The energy [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figures from the paper (6 more)
Figure 4
Figure 4. Figure 4: The left-hand side of condition 53 as a function of [PITH_FULL_IMAGE:figures/full_fig_p012_4.png]
Figure 5
Figure 5. Figure 5: The energy Ebs of the bound state (upper curve) and the energy Econd of the condensate (lower curve) are shown in the left panel. The latter is shown in the region L > L2, 57 with κ = 2, m = 1. For L0 < L < L1 see figure 3. The right panel shows the Li , 56 and 57, as …
Figure 6
Figure 6. Figure 6: In the left panel, for the Robin boundary condition, the force 54 resulting from the [PITH_FULL_IMAGE:figures/full_fig_p013_6.png]
Figure 7
Figure 7. Figure 7: The energy Ebs, 63, of the bound state solution (upper curve) and the energy Econd, 11, of the exact solution (lower curve). In both cases the parameters are m = 1, λ = 1 and R = 1. The energy of the approximate and the exact solution can be calculated explicitly. For …
Figure 8
Figure 8. Figure 8: figure 8. The vacuum energy can have either sign, the force is repulsive, as known from previous [PITH_FULL_IMAGE:figures/full_fig_p017_8.png]
Figure 9
Figure 9. Figure 9: In the left panel, the background potentials [PITH_FULL_IMAGE:figures/full_fig_p020_9.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Effects of Curvature-Scalar Coupling on Vacuum Energy in Flat (3+1)-Dimensional Space-Time

    hep-th 2026-07 conditional novelty 3.0 of 10

    In flat (3+1)-dimensional spacetime, the vacuum energy induced by an impenetrable magnetic flux tube with Robin boundary conditions depends on the curvature coupling ξ; only Dirichlet and Neumann boundary conditions m...

Reference graph

Works this paper leans on

19 extracted references · 18 canonical work pages · cited by 1 Pith paper

  1. [1]

    Bordag, G

    M. Bordag, G. L. Klimchitskaya, U. Mohideen, and V. M. Mostepanenko. Advances in the Casimir Effect. Oxford University Press, Oxford, 2009

    work page 2009

  2. [2]

    Radiative corrections as the origin of spontaneous symmetry breaking

    Sidney Coleman and Erick Weinberg. Radiative corrections as the origin of spontaneous symmetry breaking. Phys. Rev. D , 7:1888–1910, 1973

    work page 1910

  3. [3]

    A.B. Migdal. Stability of vacuum and limiting fields. Soviet Physics JETP-USSR, 34(6):1184– &, 1972. 21

    work page 1972

  4. [4]

    Lattimer

    James M. Lattimer. Neutron star equation of state. New Astron. Rev., 54(3-6):101–109, 2010

    work page 2010

  5. [5]

    Voskresenskii

    D.N. Voskresenskii. On the possibility of bose condensation of pions in ultrarelativistic colli- sions of nuclei. SovPhys-JETP, 105(6):793–805, 1994

    work page 1994

  6. [6]

    Variations of casimir energy from a super- conducting transition

    G Bimonte, E Calloni, G Esposito, and L Rosa. Variations of casimir energy from a super- conducting transition. Nucl. Phys. B , 726(3):441–463, 2005

    work page 2005

  7. [7]

    Towards measuring variations of casimir energy by a superconducting cavity

    Giuseppe Bimonte, Enrico Calloni, Giampiero Esposito, Leopoldo Milano, and Luigi Rosa. Towards measuring variations of casimir energy by a superconducting cavity. Phys. Rev. Lett., 94:180402, 2005

    work page 2005

  8. [8]

    Lally, Lawrence K

    Diego Perez-Morelo, Alexander Stange, Richard W. Lally, Lawrence K. Barrett, Matthias Imboden, Abhishek Som, David K. Campbell, Vladimir A. Aksyuk, and David J. Bishop. A system for probing casimir energy corrections to the condensation energy. Microsystems & Nanoengineering, 6(1):115, 2020

    work page 2020

Show all 19 references
  1. [9]

    Stationary solutions of the one-dimensional nonlinear schrodinger equation

    LD Carr, CW Clark, and WP Reinhardt. Stationary solutions of the one-dimensional nonlinear schrodinger equation. i. case of repulsive nonlinearity. Phys. Rev. A , 62(6):063610, 2000

    2000

  2. [10]

    Stationary solutions of the one-dimensional nonlinear schrodinger equation

    LD Carr, CW Clark, and WP Reinhardt. Stationary solutions of the one-dimensional nonlinear schrodinger equation. ii. case of attractive nonlinearity. Phys. Rev. A , 62(6):063611, 2000

    2000

  3. [11]

    Bose-einstein condensates in a one-dimensional double square well: Analytical solutions of the nonlinear schrodinger equation

    KW Mahmud, JN Kutz, and WP Reinhardt. Bose-einstein condensates in a one-dimensional double square well: Analytical solutions of the nonlinear schrodinger equation. Phys. Rev. A , 66(6):063607, 2002

    2002

  4. [12]

    Ragan, Asaad R

    Robert J. Ragan, Asaad R. Sakhel, and William J. Mullin. The gross-pitaevskii equation for an infinite square well with a delta-function barrier. J. Low Temp. Phys. , 216(5-6):814–838, 2024

    2024

  5. [13]

    Nonlinear Schr¨ odinger flow past an obstacle in one dimension

    Vincent Hakim. Nonlinear Schr¨ odinger flow past an obstacle in one dimension. Phys. Rev. E , 55:2835–2845, 1997

    1997

  6. [14]

    Espich´ an Carrillo, J. A. and Maia Jr., A. and Mostepanenko, V.M. Jacobi elliptic solutions of λϕ4 theory in a finite domain. Int. J. Mod. Phys. A , 15(17):2645–2659, 2000

    2000

  7. [15]

    Vacuum Energy for a Scalar Field With Self-Interaction in (1 + 1) Dimen- sions

    Michael Bordag. Vacuum Energy for a Scalar Field With Self-Interaction in (1 + 1) Dimen- sions. Universe, 7:55, 2021

    2021

  8. [16]

    Casimir effect for scalar fields under robin boundary conditions on plates

    A Romeo and AA Saharian. Casimir effect for scalar fields under robin boundary conditions on plates. J. Phys. A: Math. Gen. , 35(5):1297–1320, 2002

    2002

  9. [17]

    Bordag, H

    M. Bordag, H. Falomir, E. M. Santangelo, and D. V. Vassilevich. Boundary dynamics and multiple reflection expansion for Robin boundary conditions. Phys. Rev.D, 65:064032, 2002

    2002

  10. [18]

    Olver, D.W

    F.W.J. Olver, D.W. Lozier, R.F. Boisvert, and Ch.W. Clark. NIST Handbook of mathematical functions: with formulas, graphs, and mathematical tables . Cambridge university press, 2010

    2010

  11. [19]

    Vassilevich

    D.V. Vassilevich. Heat kernel expansion: User’s manual. Phys. Rept., 388:279–360, 2003. 22

    2003

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.