pith. machine review for the scientific record. sign in

arxiv: 2604.20584 · v1 · submitted 2026-04-22 · ✦ hep-th · nlin.PS

Recognition: unknown

Generalized BPS magnetic monopoles in inhomogeneous Yang-Mills-Higgs models

Filipe Rodrigues da Silva , Azadeh Mohammadi
Authors on Pith no claims yet

Pith reviewed 2026-05-09 23:59 UTC · model grok-4.3

classification ✦ hep-th nlin.PS
keywords BPS monopolesinhomogeneous mediaYang-Mills-Higgspower-law profilesshell monopolesintegrable models
0
0 comments X

The pith

Generalized Yang-Mills-Higgs models admit regular BPS monopoles for power-law inhomogeneities in a determined region of parameter space.

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

The paper constructs a non-Abelian model for magnetic monopoles where the medium has spatially varying couplings in the gauge and scalar sectors. These couplings are constrained to preserve the BPS bound, allowing first-order equations to be solved for static spherically symmetric fields. For a specific class of permeabilities depending on a power-law radial profile, the authors map out where smooth solutions exist and find an integrable line where exact solutions can be written down. This produces a variety of monopole profiles including compact cores, hollow regions, and multiple concentric shells in the energy density.

Core claim

By imposing the product of the generalized permittivity and permeability to equal unity, the model retains a BPS bound. For permeabilities of the form f(r)/H^α with f(r) = r^β, regular BPS monopoles exist in a specific domain of the (α, β) plane, with exact closed-form solutions available when α equals 1.

What carries the argument

The constraint that the product of the generalized permittivity P and permeability M equals one, which keeps the energy bounded below by the topological charge and reduces the equations to first order.

If this is right

  • Regular monopole solutions exist for a range of exponents α and β in the power-law profile.
  • Along the line α=1 the equations integrate exactly, giving closed-form expressions for the monopole fields in the inhomogeneous background.
  • The energy density can exhibit single peaks, multiple concentric peaks, or be concentrated near the origin or in shells.
  • Numerical solutions fill in the parameter space away from the integrable line.

Where Pith is reading between the lines

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

  • This approach could be used to model monopoles in materials with radially varying properties, such as graded dielectrics or superconductors.
  • Multi-shell configurations might correspond to bound states or resonances in the inhomogeneous medium.
  • The exact solutions for α=1 provide a benchmark for testing numerical methods in more general inhomogeneous settings.

Load-bearing premise

The model requires the product of the generalized permittivity and permeability to be tuned to exactly one so that the BPS bound remains saturated.

What would settle it

Failure to find regular solutions numerically in the claimed domain of the (α, β) plane, or explicit substitution showing that the closed-form expressions for α=1 do not satisfy the first-order equations.

Figures

Figures reproduced from arXiv: 2604.20584 by Azadeh Mohammadi, Filipe Rodrigues da Silva.

Figure 1
Figure 1. Figure 1: Regularity map of the BPS equations (20) in the (α, β) plane, indicating the singular behavior of the system at the origin across different regions of the plane. In the interval 0 ≤ α ≤ 1, the system contains no negative powers of H, and the singularity structure is controlled entirely by the powers of r. In this region, the singularities are therefore of pole type, and their order is determined by the val… view at source ↗
Figure 2
Figure 2. Figure 2: Admissible region in the (α, β) plane. The blue curve corresponds to the parabola α(β) = 1 + (β + 1)2/8 and the dashed line to α = 1. A. The line α = 1 For α = 1, the general solution derived in Sec. II can be evaluated analytically. Using Eqs. (17) and (18), the integrals reduce to I1(r) = Z r 0 t β dt = r β+1 β + 1 , (22) which is finite provided β + 1 > 0, precisely the condition already identified in t… view at source ↗
Figure 3
Figure 3. Figure 3: BPS monopole solutions along the analytical line [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Field profiles and energy density for α = 1 and β ∈ [2, 6]. From left to right: scalar profile H(r), gauge profile K(r), and energy density E(r) [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: A shell-like monopole obtained along the analytical line [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: The cavity radius r⋆ as a function of β. Thus, even within the analytically solvable line α = 1, the model displays a wide variety of monopole configurations. Depending on the value of β, one finds compact core-like monopoles, broadened monopoles, and hollow shell-like monopoles. The parameter β therefore controls how the energy is redistributed in space, interpolating between configurations concentrated n… view at source ↗
Figure 7
Figure 7. Figure 7: BPS monopole solutions along the parabola [PITH_FULL_IMAGE:figures/full_fig_p014_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Monopole profiles and energy densities along the line [PITH_FULL_IMAGE:figures/full_fig_p015_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Comparison of energy densities along constant- [PITH_FULL_IMAGE:figures/full_fig_p016_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Energy densities for two examples of multi-shell monopoles, one with [PITH_FULL_IMAGE:figures/full_fig_p016_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: The standard Prasad-Sommerfield monopole at ( [PITH_FULL_IMAGE:figures/full_fig_p017_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Monopole profiles and energy densities along the line [PITH_FULL_IMAGE:figures/full_fig_p018_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Emergence and growth of hollow-monopole configurations along the line [PITH_FULL_IMAGE:figures/full_fig_p018_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Monopole solutions along the line α = −2 for negative values of β. The observed results show that as α decreases further below unity while β remains above the limiting value −1, the monopole configurations become increasingly spatially extended and less localized. Consequently, progressively larger values of β are required for shell-like monopoles to emerge. This behavior reinforces the global picture tha… view at source ↗
Figure 15
Figure 15. Figure 15: Left: surface m+(α, β), which grows towards α = 1 and negative values of β, but due to the boundary β = −1 acquires the form of a cusp accumulating near the vertex (1, −1). Right: surface m−(α, β), which grows towards α = 1 and positive values of β. For the numerical implementation, we regularize the problem near ξ = 0 by introducing a small cutoff ξa ≪ 1 and replacing the exact boundary conditions at the… view at source ↗
Figure 16
Figure 16. Figure 16: Analytical behavior of the parameter a0(β) for α = 1, showing its divergence as β → −1 + and its asymptotic approach to unity for large β [PITH_FULL_IMAGE:figures/full_fig_p028_16.png] view at source ↗
read the original abstract

We present a non-Abelian model for magnetic monopoles in inhomogeneous media, based on a generalization of the standard 't~Hooft-Polyakov model. The medium is described by spatially dependent couplings in the gauge and scalar sectors, constrained by $P(|\Phi|,r)M(|\Phi|,r)=1$ so that the Bogomol'nyi-Prasad-Sommerfield (BPS) bound is preserved. For static spherically symmetric configurations, we study the first-order monopole equations for the class of generalized permeabilities $M(H,r)=f(r)/H^\alpha$. For the power-law profile $f(r)=r^\beta$, we determine the domain in the $(\alpha,\beta)$ plane where regular BPS solutions exist. On the line $\alpha=1$, the system becomes exactly integrable, with closed-form monopole solutions in an inhomogeneous background. Away from this analytical sector, the solutions are constructed numerically. The model supports a rich spectrum of configurations, including effectively point-like monopoles, compact-core monopoles, hollow monopoles, shell-like structures, and multi-shell monopoles characterized by multiple concentric peaks in the energy density.

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

0 major / 2 minor

Summary. The manuscript presents a generalization of the 't Hooft-Polyakov monopole model to inhomogeneous media by introducing spatially dependent couplings in the Yang-Mills and Higgs sectors, subject to the constraint P(|Φ|, r) M(|Φ|, r) = 1 that preserves the BPS bound. For static spherically symmetric configurations, the authors derive the first-order BPS equations for the class of permeabilities M(H, r) = f(r)/H^α with power-law f(r) = r^β. They identify the region in the (α, β) parameter plane where regular solutions exist, provide exact analytic solutions when α = 1, and construct numerical solutions otherwise. The resulting configurations include point-like, compact, hollow, shell-like, and multi-shell monopoles.

Significance. If the results hold, the work provides a controlled extension of BPS monopoles to inhomogeneous backgrounds while retaining first-order equations and the energy bound. The exact integrability along α=1 yields closed-form solutions that serve as benchmarks, and the numerical survey maps a concrete domain in parameter space supporting a spectrum of distinct profiles. This construction is valuable for exploring topological defects in media with spatial variation and supplies explicit examples that can be used for further analytic or phenomenological study.

minor comments (2)
  1. The abstract refers to 'effectively point-like monopoles' and 'hollow monopoles'; a short parenthetical definition in terms of the radial energy-density profile (e.g., location of the peak relative to the core radius) would improve immediate readability.
  2. [§4] In the numerical section, the boundary conditions imposed at r=0 and r→∞ to enforce regularity and finite energy should be stated explicitly, together with the integrator and tolerance used to delineate the (α,β) domain.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive and accurate summary of our manuscript, as well as for recommending acceptance. Their assessment correctly identifies the key features of the generalized BPS monopoles in inhomogeneous Yang-Mills-Higgs models.

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained

full rationale

The paper defines its model by explicitly imposing the algebraic constraint P(|Φ|,r)M(|Φ|,r)=1 as part of the construction to enforce preservation of the BPS bound and reduce to first-order equations. It then substitutes the specific ansatz M(H,r)=f(r)/H^α with power-law f(r)=r^β, integrates the resulting ODEs (analytically on α=1, numerically elsewhere), and reports the domain of regular solutions together with the observed profile types. No output quantity is used to define or fit an input parameter, no self-citation supplies a load-bearing uniqueness theorem, and the claimed spectrum follows directly from solving the stated equations rather than from any renaming or re-derivation of the inputs themselves.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The model starts from the standard Yang-Mills-Higgs Lagrangian and adds only position-dependent factors whose product is fixed to one. The only free choices are the exponents α and β that label the family being studied; no new particles or forces are postulated.

free parameters (2)
  • α
    Exponent controlling the Higgs-field dependence in the permeability M(H,r)=f(r)/H^α; defines the class of models examined.
  • β
    Exponent in the radial power-law f(r)=r^β; scanned to locate the region of regular solutions.
axioms (2)
  • domain assumption P(|Φ|,r) M(|Φ|,r)=1 preserves the BPS bound
    Imposed so that the energy is bounded by the topological charge and first-order equations exist.
  • domain assumption Static spherically symmetric field configurations
    Reduces the PDE system to ODEs in the radial coordinate.

pith-pipeline@v0.9.0 · 5505 in / 1824 out tokens · 49088 ms · 2026-05-09T23:59:17.376703+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

34 extracted references · 24 canonical work pages

  1. [1]

    Since eH(0) = 0 and eK(0) = 1, we may write eH(ξ) =δ eH(ξ), eK(ξ) = 1−δ eK(ξ)

    Nearξ= 0 Near the origin, it is convenient to introduce the deviations δ eH(ξ)≡ eH(0)− eH(ξ) , δ eK(ξ)≡ eK(0)− eK(ξ) ,(A5) which vanish atξ= 0 and grow toward unity asξ→1, consistent with (A4). Since eH(0) = 0 and eK(0) = 1, we may write eH(ξ) =δ eH(ξ), eK(ξ) = 1−δ eK(ξ). For regular monopole solutions, both eH(ξ) and eK(ξ) must approach their boundary va...

  2. [2]

    ∞X k=0 β k (−t)k # ∞X n=0 ( ¯H1−α)(n)(0) n! tn ! dt    = ∞Y k,n=0 exp ( β k (−1)k ( ¯H1−α)(n)(0) n!

    Nearξ= 1 To analyze the behavior of the solutions near spatial infinity, we introduce the variable τ≡1−ξ,(A12) so thatτ→0 corresponds tor→ ∞. We further define ¯H(τ)≡ eH(ξ(τ)), ¯K(τ)≡ eK(ξ(τ)).(A13) 25 In terms of these variables, the BPS equations can be written as ¯K(τ) = exp Z τ 1 t−(β+2)(1−t) β ¯H 1−α(t) dt ,(A14a) ¯K2 = 1 +τ −β(1−τ) β+2 ¯H −α ¯H ′,(A...

  3. [3]

    (A11), up tor=r b

    One way to proceed is defining this parameter as the root of a function shooting(a0), constructed as follows: for a given trial value ofa 0, the BPS system is integrated from the initial conditions (H a, Ka) atr=r a, specified by Eq. (A11), up tor=r b. The function shooting(a0) then returns the mismatch shooting(a0) =H tar −H b,(A24) whereH b denotes the ...

  4. [4]

    Armend´ ariz-Pic´ on, T

    C. Armend´ ariz-Pic´ on, T. Damour, and V. Mukhanov, Physics Letters B458, 209 (1999), ISSN 0370- 2693, URLhttps://www.sciencedirect.com/science/article/pii/S0370269399006036

    1999

  5. [5]

    Tinto and J

    E. Babichev, Phys. Rev. D74, 085004 (2006), URLhttps://link.aps.org/doi/10.1103/PhysRevD. 74.085004

    work page doi:10.1103/physrevd 2006

  6. [6]

    C. Adam, J. S´ anchez-Guill´ en, and A. Wereszczy´ nski, Journal of Physics A: Mathematical and Theo- retical40, 13625 (2007), URLhttps://doi.org/10.1088/1751-8113/40/45/009

    work page doi:10.1088/1751-8113/40/45/009 2007

  7. [7]

    C. Adam, J. Sanchez-Guillen, and A. Wereszczynski, Journal of Physics A: Mathematical and Theo- retical42, 089801 (2009), URLhttps://doi.org/10.1088/1751-8121/42/8/089801

    work page doi:10.1088/1751-8121/42/8/089801 2009

  8. [8]

    Manifestly crossing-invariant parametrization of n-meson amplitude,

    G. ’t Hooft, Nuclear Physics B79, 276 (1974), ISSN 0550-3213, URLhttps://www.sciencedirect. com/science/article/pii/0550321374904866

    work page arXiv 1974

  9. [9]

    A. M. Polyakov, JETP Lett.20, 194 (1974), [Pisma Zh. Eksp. Teor. Fiz. 20, 430 (1974)]

    1974

  10. [10]

    E. B. Bogomolny, Sov. J. Nucl. Phys.24, 449 (1976), [Yad. Fiz. 24, 861 (1976)]

    1976

  11. [11]

    M. K. Prasad and C. M. Sommerfield, Phys. Rev. Lett.35, 760 (1975), URLhttps://link.aps.org/ doi/10.1103/PhysRevLett.35.760

    work page doi:10.1103/physrevlett.35.760 1975

  12. [12]

    Casana, M

    R. Casana, M. M. Ferreira, and E. da Hora, Phys. Rev. D86, 085034 (2012), URLhttps://link. aps.org/doi/10.1103/PhysRevD.86.085034

    work page doi:10.1103/physrevd.86.085034 2012

  13. [13]

    A. N. Atmaja and I. Prasetyo, Advances in High Energy Physics2018, 7376534 (2018), URLhttps: //doi.org/10.1155/2018/7376534

    work page doi:10.1155/2018/7376534 2018

  14. [14]

    A. N. Atmaja, Eur. Phys. J. C82, 602 (2022), URLhttps://doi.org/10.1140/epjc/ s10052-022-10569-6

    work page doi:10.1140/epjc/ 2022

  15. [15]

    Mulyanto, E. S. Fadhilla, and A. N. Atmaja, Physics Letters B874, 140267 (2026), ISSN 0370-2693, URLhttps://doi.org/10.1016/j.physletb.2026.140267

    work page doi:10.1016/j.physletb.2026.140267 2026

  16. [16]

    Bazeia, R

    D. Bazeia, R. Casana, M. Ferreira, E. da Hora, and L. Losano, Physics Letters B727, 548 (2013), ISSN 0370-2693, URLhttps://www.sciencedirect.com/science/article/pii/S0370269313008939

    2013

  17. [17]

    Bazeia, M

    D. Bazeia, M. A. Marques, and R. Menezes, Phys. Rev. D97, 105024 (2018), URLhttps://link. aps.org/doi/10.1103/PhysRevD.97.105024

    work page doi:10.1103/physrevd.97.105024 2018

  18. [18]

    Bazeia, M

    D. Bazeia, M. A. Marques, and G. J. Olmo, Phys. Rev. D98, 025017 (2018), URLhttps://link. aps.org/doi/10.1103/PhysRevD.98.025017

    work page doi:10.1103/physrevd.98.025017 2018

  19. [19]

    Bazeia, M

    D. Bazeia, M. A. Marques, and R. Menezes, Phys. Rev. D98, 065003 (2018), URLhttps://link. aps.org/doi/10.1103/PhysRevD.98.065003

    work page doi:10.1103/physrevd.98.065003 2018

  20. [20]

    Bazeia, M

    D. Bazeia, M. A. Liao, and M. A. Marques, The European Physical Journal C81, 552 (2021), URL https://doi.org/10.1140/epjc/s10052-021-09352-w

    work page doi:10.1140/epjc/s10052-021-09352-w 2021

  21. [21]

    Beneˇ s and F

    P. Beneˇ s and F. Blaschke, Phys. Rev. D107, 125002 (2023), URLhttps://link.aps.org/doi/10. 1103/PhysRevD.107.125002

    2023

  22. [22]

    Beneˇ s and F

    P. Beneˇ s and F. Blaschke, Progress of Theoretical and Experimental Physics2026, 023B04 (2026), 30 ISSN 2050-3911, URLhttps://doi.org/10.1093/ptep/ptag004

    work page doi:10.1093/ptep/ptag004 2026

  23. [23]

    C. Adam, J. M. Queiruga, and A. Wereszczynski, J. High Energ. Phys.07, 164 (2019), URLhttps: //doi.org/10.1007/JHEP07(2019)164

    work page doi:10.1007/jhep07(2019)164 2019

  24. [24]

    C. Adam, T. Romanczukiewicz, and A. Wereszczynski, J. High Energ. Phys.03, 131 (2019), URL https://doi.org/10.1007/JHEP03(2019)131

    work page doi:10.1007/jhep03(2019)131 2019

  25. [25]

    C. Adam, K. Oles, J. M. Queiruga, T. Romanczukiewicz, and A. Wereszczynski, J. High Energ. Phys. 07, 150 (2019), URLhttps://doi.org/10.1007/JHEP07(2019)150

    work page doi:10.1007/jhep07(2019)150 2019

  26. [26]

    J. a. G. F. Campos and A. Mohammadi, Phys. Rev. D102, 045003 (2020), URLhttps://link.aps. org/doi/10.1103/PhysRevD.102.045003

    work page doi:10.1103/physrevd.102.045003 2020

  27. [27]

    S lawi´ nska, Phys

    K. S lawi´ nska, Phys. Rev. E111, 014228 (2025), URLhttps://link.aps.org/doi/10.1103/ PhysRevE.111.014228

    2025

  28. [28]

    Tong and K

    D. Tong and K. Wong, J. High Energ. Phys.01, 090 (2014), URLhttps://doi.org/10.1007/ JHEP01(2014)090

    2014

  29. [29]

    Cockburn, S

    A. Cockburn, S. Krusch, and A. A. Muhamed, J. Math. Phys.58, 063509 (2017), URLhttps://doi. org/10.1063/1.4984980

    work page doi:10.1063/1.4984980 2017

  30. [30]

    Almeida, R

    V. Almeida, R. Casana, E. da Hora, and S. Krusch, Phys. Rev. D106, 016010 (2022), URLhttps: //link.aps.org/doi/10.1103/PhysRevD.106.016010

    work page doi:10.1103/physrevd.106.016010 2022

  31. [31]

    Bazeia, J

    D. Bazeia, J. G. F. Campos, and A. Mohammadi, J. High Energ. Phys.12, 108 (2024), URLhttps: //doi.org/10.1007/JHEP12(2024)108

    work page doi:10.1007/jhep12(2024)108 2024

  32. [32]

    N. H. Gonzalez-Gutierrez, R. Casana, and A. C. Santos, J. High Energ. Phys.01, 005 (2026), URL https://doi.org/10.1007/JHEP01(2026)005

    work page doi:10.1007/jhep01(2026)005 2026

  33. [33]

    Bazeia, M

    D. Bazeia, M. A. Marques, and R. Menezes, Eur. Phys. J. C85, 836 (2025), URLhttps://doi.org/ 10.1140/epjc/s10052-025-14582-3

    work page doi:10.1140/epjc/s10052-025-14582-3 2025

  34. [34]

    Bazeia, M

    D. Bazeia, M. A. Liao, and M. A. Marques, arXiv preprint (2025), 2505.11215v2, URLhttps://arxiv. org/abs/2505.11215v2

    work page arXiv 2025