Subbarao–Warren Problem: Multidisciplinary Views

• The Subbarao–Warren Problem is a collective term for distinct research challenges in thermostatistics, phase-field modeling, arithmetic congruences, partition theory, and unitary perfect numbers.
• It emphasizes ensemble-specific criteria and structured analytical reductions, such as verifying Gibbs temperature consistency and deriving nonlinear jump cost models.
• Each branch leverages rigorous methods—from microcanonical equilibrium proofs to prime classification and partition bijections—to address domain-specific puzzles.

Searching arXiv for the cited Subbarao/Warren-related papers to ground the article. The Subbarao–Warren Problem does not denote a single universally standardized problem across mathematics and physics; rather, the label is used for several distinct research threads associated with Subbarao-, Warren-, or Subbarao–Warren-type questions. In the supplied literature, it refers to at least four technically separate domains: a thermostatistical dispute over whether the Gibbs temperature $T_G$ correctly characterizes microcanonical equilibrium (Dunkel et al., 2014); sharp-interface limits of Kobayashi–Warren–Carter phase-field energies leading to a jump-sensitive total-variation-type functional (Giga et al., 2024); arithmetic congruence problems involving $\varphi(n)$ and $\sigma(n)$, including a variation of a congruence of Subbarao (Bujačić, 2016); and partition-theoretic generalizations of Subbarao’s finitization and Bressoud–Subbarao weighted identities (Nyirenda et al., 2022, Agarwal et al., 2022). A more recent use concerns the Subbarao–Warren problem for unitary perfect numbers, formulated through a structured multiplicative balance and a dependency-graph analysis (Maciejewski, 19 May 2026). The term therefore functions as an umbrella for a family of problems rather than a single theorem statement.

1. Thermostatistical formulation: Gibbs temperature and thermal equilibrium

In statistical mechanics, the central question is explicitly stated as follows: What temperature correctly describes thermal equilibrium for an isolated system? For a confined classical system with Hamiltonian

$H(\zeta,A)=E,$

where $\zeta=(\zeta_1,\ldots,\zeta_N)$ are canonical coordinates and $A=(A_\mu)$ are external control parameters, the answer advanced by Dunkel and Hilbert is that the Gibbs temperature does (Dunkel et al., 2014).

The paper fixes the thermodynamic objects by defining the microcanonical density operator $\rho_M$, supported on the shell $H=E$, together with the density of states and integrated density of states,

$$\omega(E,A)=\mathrm{Tr}[\delta(E-H)], \qquad \Omega(E,A)=\mathrm{Tr}[\Theta(E-H)].$$

With $k_B=1$, the entropies and temperatures are

$\varphi(n)$0

$\varphi(n)$1

The canonical density operator and Shannon entropy are also defined: $\varphi(n)$2 These definitions organize the dispute between Boltzmann, Gibbs, and Shannon entropy assignments in different ensembles (Dunkel et al., 2014).

The key benchmark is the thermodynamic consistency condition

$\varphi(n)$3

Three exact facts are then highlighted. First,

$\varphi(n)$4

so the Gibbs temperature satisfies microcanonical equipartition for all finite $\varphi(n)$5, whereas $\varphi(n)$6 does not (Dunkel et al., 2014). Second, the pair $\varphi(n)$7 satisfies the consistency relation for all finite $\varphi(n)$8, whereas $\varphi(n)$9 does not. Third, $\sigma(n)$0 satisfies the same consistency relation for all finite $\sigma(n)$1.

The rebuttal to Frenkel and Warren is that the consistency of $\sigma(n)$2 does not undermine the consistency of $\sigma(n)$3; these are ensemble-specific statements, not mutually exclusive alternatives. This suggests that the “problem” in this formulation is one of ensemble-appropriate thermodynamic description, rather than a competition to identify a single entropy valid in all settings (Dunkel et al., 2014).

2. Coupled isolated systems and the equilibrium criterion

The claim that $\sigma(n)$4 does not characterize thermal equilibrium is answered through a weak-coupling setup involving two isolated systems with Hamiltonian

$\sigma(n)$5

where $\sigma(n)$6. If the systems initially have energies $\sigma(n)$7 and $\sigma(n)$8, then after coupling the total energy is

$\sigma(n)$9

and the joint state is microcanonical: $H(\zeta,A)=E,$0 In the weak-coupling limit, microcanonical equipartition implies

$H(\zeta,A)=E,$1

The equilibrium statement is therefore that when two isolated systems are brought into thermal contact and described microcanonically, they equilibrate to a common Gibbs temperature (Dunkel et al., 2014).

Within this framework, the authors conclude that $H(\zeta,A)=E,$2 does correctly characterize thermal equilibrium and that $H(\zeta,A)=E,$3 does not correctly characterize equilibrium for finite systems. They further note that the common claim that “energy per particle becomes equally distributed” can be wrong for some finite systems; the quantity that equalizes is the virial/equipartition observable $H(\zeta,A)=E,$4, not necessarily a simple energy-per-particle share (Dunkel et al., 2014).

The discussion extends to the thermodynamic limit and to bounded-spectrum systems. The exact facts are said to hold for every finite $H(\zeta,A)=E,$5, so if different limits yield different answers, that indicates ensemble inequivalence, not failure of Gibbs temperature. The paper emphasizes that microcanonical and canonical ensembles are not equivalent in systems with bounded spectra, especially for population-inversion systems, spins, and ultracold gases. The correct ensemble must be matched to the experimental situation: canonical for a system coupled to an effectively infinite heat bath, microcanonical for an isolated system with fixed energy (Dunkel et al., 2014).

3. Phase-field and free-discontinuity formulation: the Kobayashi–Warren–Carter limit

A separate use of the label concerns a variational model tied to the Kobayashi–Warren–Carter energy. The functional studied is of Rudin–Osher–Fatemi type, but with a regularizer that is not classical total variation. For $H(\zeta,A)=E,$6,

$H(\zeta,A)=E,$7

and the fidelity-augmented functional is

$H(\zeta,A)=E,$8

Here the regularizer measures jump discontinuities through a nonlinear cost $H(\zeta,A)=E,$9, rather than penalizing total variation linearly (Giga et al., 2024).

The paper explains that $\zeta=(\zeta_1,\ldots,\zeta_N)$0 arises by minimizing the diffuse KWC energy over the order parameter $\zeta=(\zeta_1,\ldots,\zeta_N)$1 in the sharp-interface limit $\zeta=(\zeta_1,\ldots,\zeta_N)$2: $\zeta=(\zeta_1,\ldots,\zeta_N)$3 with

$\zeta=(\zeta_1,\ldots,\zeta_N)$4

In the limit, minimizing over the phase-field object $\zeta=(\zeta_1,\ldots,\zeta_N)$5 yields the effective jump cost

$\zeta=(\zeta_1,\ldots,\zeta_N)$6

where

$\zeta=(\zeta_1,\ldots,\zeta_N)$7

For the model choice $\zeta=(\zeta_1,\ldots,\zeta_N)$8, $\zeta=(\zeta_1,\ldots,\zeta_N)$9, the paper gives explicitly

$A=(A_\mu)$0

This establishes a rigorous passage from a diffuse-interface KWC model to a free-discontinuity functional with a nonlinear jump term (Giga et al., 2024).

The main one-dimensional theorem states that if $A=(A_\mu)$1 and $A=(A_\mu)$2 satisfies the structural assumptions, then any minimizer $A=(A_\mu)$3 of $A=(A_\mu)$4 is piecewise constant with finitely many jumps. More precisely, $A=(A_\mu)$5 is piecewise constant on $A=(A_\mu)$6, satisfies $A=(A_\mu)$7, and its number of jumps $A=(A_\mu)$8 obeys

$A=(A_\mu)$9

with a sharper estimate for monotone $\rho_M$0,

$\rho_M$1

The proof uses the coincidence set

$\rho_M$2

together with comparison arguments showing that outside $\rho_M$3 the minimizer cannot vary smoothly and is forced to flatten into constant pieces (Giga et al., 2024).

The paper emphasizes the contrast with classical ROF: for

$\rho_M$4

continuous data produce a minimizer with no jumps, whereas in the KWC-derived model the minimizer may develop jumps even when $\rho_M$5 is continuous. This suggests that the variational branch of the Subbarao–Warren label concerns phase separation and interface generation, not ordinary smoothing (Giga et al., 2024).

4. Arithmetic congruence problems involving $\rho_M$6 and $\rho_M$7

A number-theoretic branch concerns congruences involving Euler’s totient function and the divisor-sum function. One paper studies the congruence

$\rho_M$8

described as a variation of a congruence of Subbarao, because Subbarao studied the analogous congruence

$\rho_M$9

The roles of $H=E$0 and $H=E$1 are interchanged, so the work addresses a dual variation rather than the classical congruence itself (Bujačić, 2016).

The main theorem gives a complete classification in the family

$H=E$2

The only such integers satisfying

$H=E$3

are

$H=E$4

This is proved by using the multiplicative formulas

$H=E$5

when both exponents are positive, and

$H=E$6

The proof treats prime cases, pure powers of $H=E$7, pure powers of $H=E$8, and the mixed case separately (Bujačić, 2016).

In the mixed case, setting

$H=E$9

the congruence leads to

$$\omega(E,A)=\mathrm{Tr}[\delta(E-H)], \qquad \Omega(E,A)=\mathrm{Tr}[\Theta(E-H)].$$0

Further modular arguments show that $$\omega(E,A)=\mathrm{Tr}[\delta(E-H)], \qquad \Omega(E,A)=\mathrm{Tr}[\Theta(E-H)].$$1 and $$\omega(E,A)=\mathrm{Tr}[\delta(E-H)], \qquad \Omega(E,A)=\mathrm{Tr}[\Theta(E-H)].$$2 must both be even and that $$\omega(E,A)=\mathrm{Tr}[\delta(E-H)], \qquad \Omega(E,A)=\mathrm{Tr}[\Theta(E-H)].$$3. The problem is then rewritten as

$$\omega(E,A)=\mathrm{Tr}[\delta(E-H)], \qquad \Omega(E,A)=\mathrm{Tr}[\Theta(E-H)].$$4

with

$$\omega(E,A)=\mathrm{Tr}[\delta(E-H)], \qquad \Omega(E,A)=\mathrm{Tr}[\Theta(E-H)].$$5

and diagonalized into the Pell-type equation

$$\omega(E,A)=\mathrm{Tr}[\delta(E-H)], \qquad \Omega(E,A)=\mathrm{Tr}[\Theta(E-H)].$$6

Congruence restrictions yield

$$\omega(E,A)=\mathrm{Tr}[\delta(E-H)], \qquad \Omega(E,A)=\mathrm{Tr}[\Theta(E-H)].$$7

and continued-fraction methods of Worley–Dujella and Dujella–Jadrijević reduce the possibilities to a finite set of values of $$\omega(E,A)=\mathrm{Tr}[\delta(E-H)], \qquad \Omega(E,A)=\mathrm{Tr}[\Theta(E-H)].$$8, all of which are eliminated. Hence there are no solutions with $$\omega(E,A)=\mathrm{Tr}[\delta(E-H)], \qquad \Omega(E,A)=\mathrm{Tr}[\Theta(E-H)].$$9 (Bujačić, 2016).

This branch of the topic belongs to a broader family of arithmetic characterization problems similar in spirit to Lehmer’s totient problem and Subbarao’s congruence questions. A plausible implication is that the expression “Subbarao–Warren Problem” is being used here to denote a family of totient/divisor-sum congruence classification problems, rather than one fixed conjecture.

5. Partition-theoretic generalizations: finitization and weighted identities

Another major branch concerns partition theory. One paper generalizes Subbarao’s finitization of Andrews’ theorem. It recalls the chain of results from MacMahon’s theorem and Andrews’ extension to Subbarao’s finitized form

$k_B=1$0

where $k_B=1$1 and $k_B=1$2 encode multiplicity and residue restrictions on partitions (Nyirenda et al., 2022).

The generalized finitization is stated as Theorem 3.1. For positive integers $k_B=1$3 and $k_B=1$4 with $k_B=1$5, the classes $k_B=1$6 and $k_B=1$7 are defined so that

$k_B=1$8

The multiplicity conditions in $k_B=1$9 are

$\varphi(n)$00

for multiplicities congruent to $\varphi(n)$01. The residue conditions in $\varphi(n)$02 require that parts divisible by $\varphi(n)$03 are not divisible by $\varphi(n)$04, while parts not divisible by $\varphi(n)$05 are congruent to

$\varphi(n)$06

The paper explicitly notes that setting $\varphi(n)$07 and $\varphi(n)$08 reduces Theorem 3.1 to Subbarao’s finitization (Nyirenda et al., 2022).

A central feature is the bijection

$\varphi(n)$09

whose construction uses base-$\varphi(n)$10 expansion in one case and base-$\varphi(n)$11 expansion in the other. Passing to the limit $\varphi(n)$12, the paper defines

$\varphi(n)$13

and proves

$\varphi(n)$14

The authors state that this extends Sellers’ bijection and the Sellers–Fu bijection by treating all possible residue classes modulo $\varphi(n)$15 rather than fixing two classes (Nyirenda et al., 2022).

Arithmetic consequences are also derived. Among them are parity congruences such as

$\varphi(n)$16

whenever $\varphi(n)$17 is a quadratic nonresidue modulo $\varphi(n)$18, and

$\varphi(n)$19

whenever $\varphi(n)$20 is a quadratic nonresidue modulo $\varphi(n)$21, for the prime ranges stated in the paper (Nyirenda et al., 2022).

A related direction concerns Bressoud–Subbarao type weighted partition identities for the generalized divisor function

$\varphi(n)$22

The central theorem states that for any $\varphi(n)$23 and complex $\varphi(n)$24,

$\varphi(n)$25

Setting $\varphi(n)$26 yields the generalized divisor-function identity originally due to Bressoud and Subbarao,

$\varphi(n)$27

and the specialization $\varphi(n)$28, $\varphi(n)$29 gives

$\varphi(n)$30

The proof uses a sign-reversing pairing on $\varphi(n)$31 and extends the original argument to complex $\varphi(n)$32 and the extra parameter $\varphi(n)$33 (Agarwal et al., 2022).

The same paper also derives identities from Ramanujan’s $\varphi(n)$34-series, Uchimura’s identity, Dilcher-type formulas, and an Andrews–Garvan–Liang identity. For example,

$\varphi(n)$35

and

$\varphi(n)$36

These results place the partition-theoretic Subbarao–Warren usage in a broader framework of weighted partition identities for divisor functions (Agarwal et al., 2022).

6. Unitary perfect numbers and the modern bounded-box reduction

A recent formulation explicitly speaks of the Subbarao–Warren problem for unitary perfect numbers. A unitary perfect number is a positive integer $\varphi(n)$37 satisfying

$\varphi(n)$38

where $\varphi(n)$39 sums unitary divisors. Writing

$\varphi(n)$40

with odd primes $\varphi(n)$41, unitary multiplicativity gives

$\varphi(n)$42

so the defining equation becomes

$\varphi(n)$43

This exact multiplicative balance is the organizing principle of the paper (Maciejewski, 19 May 2026).

The analysis uses the recursive notion of a 3-Higgs prime. The structural fact stated is that every prime divisor of a unitary perfect number must be 3-Higgs. This motivates the set

$\varphi(n)$44

and its even part

$\varphi(n)$45

An odd dependency graph is then introduced, with vertices equal to odd 3-Higgs primes and edges $\varphi(n)$46 whenever $\varphi(n)$47 for some admissible exponent $\varphi(n)$48. Strongly connected components of this graph are interpreted as feedback kernels in the prime-forcing cascade (Maciejewski, 19 May 2026).

Within the bounded box

$\varphi(n)$49

every admissible source kernel is either one of the two kernels occurring in the known nonsquarefree examples,

$\varphi(n)$50

or one of five impostor kernels: $\varphi(n)$51 The paper gives a reproducible three-filter certificate eliminating those impostor kernels for all relevant seed classes with $\varphi(n)$52 (Maciejewski, 19 May 2026).

The three filters are:

Filter	Description
Z	Zsigmondy-type exponent obstructions
N	inherited non-3-Higgs witnesses
O	deterministic 2-adic budget overshoot

The exact split among 2119 impostor candidates up to $\varphi(n)$53 is reported as 495 killed by Z, 1614 killed by N, 10 killed by O, and 0 unresolved (Maciejewski, 19 May 2026).

The remaining obstruction is the auxiliary set $\varphi(n)$54. A structural lemma shows that if

$\varphi(n)$55

then every prime divisor of $\varphi(n)$56 is 3-Higgs, each such prime appears with exponent at most $\varphi(n)$57, and every odd divisor $\varphi(n)$58 gives $\varphi(n)$59. The paper proves that $\varphi(n)$60 is finite if and only if the subset

$\varphi(n)$61

is finite, and if there are $\varphi(n)$62 such prime-branch elements then

$\varphi(n)$63

This reduces the open problem to the prime branch $\varphi(n)$64 (Maciejewski, 19 May 2026).

Certified frontier bounds are then established: $\varphi(n)$65 The paper also uses Ford’s theorem for downward-closed prime sets to show that the set of 3-Higgs primes is thin,

$\varphi(n)$66

for some $\varphi(n)$67, and in particular

$\varphi(n)$68

However, this gives thinness rather than finiteness. The unresolved branch is reformulated as a divisor-level problem for the cyclotomic value $\varphi(n)$69, equivalently for the factors in

$\varphi(n)$70

The remaining target is to show that for all sufficiently large odd primes $\varphi(n)$71, $\varphi(n)$72 has at least one prime divisor that is not 3-Higgs (Maciejewski, 19 May 2026).

7. Scope, ambiguity, and recurring structural themes

Across these sources, the phrase Subbarao–Warren Problem is not attached to a single invariant mathematical object. Instead, it denotes a family of problems unified only loosely by naming conventions and by the appearance of Subbarao- or Warren-linked structures. The supplied literature supports at least the following usages:

Domain	Representative problem	arXiv id
Statistical mechanics	Whether $\varphi(n)$73 correctly characterizes microcanonical equilibrium	(Dunkel et al., 2014)
Phase-field / free discontinuity	Sharp-interface limit of KWC energy and piecewise-constant minimizers	(Giga et al., 2024)
Arithmetic congruences	Classification for $\varphi(n)$74 in $\varphi(n)$75	(Bujačić, 2016)
Partition theory	Generalized finitization and weighted divisor-partition identities	(Nyirenda et al., 2022, Agarwal et al., 2022)
Unitary perfect numbers	Bounded-box reduction via source kernels and 3-Higgs primes	(Maciejewski, 19 May 2026)

Several recurring themes nonetheless appear. One is the replacement of a naive global principle by an ensemble-specific, model-specific, or class-specific criterion: Gibbs vs. Shannon entropy in thermostatistics (Dunkel et al., 2014), nonlinear jump cost $\varphi(n)$76 vs. classical total variation in KWC limits (Giga et al., 2024), restricted prime-support families in arithmetic classification (Bujačić, 2016), residue-class/multiplicity frameworks in partition theory (Nyirenda et al., 2022), and source-kernel plus dependency-graph reductions in the unitary perfect-number setting (Maciejewski, 19 May 2026).

A second recurring feature is reduction to structured exact statements. Examples include the three exact facts E1–E3 in the Gibbs-temperature dispute (Dunkel et al., 2014), the gamma-limit jump-energy formula in the KWC setting (Giga et al., 2024), the Pell-type reduction in the congruence problem (Bujačić, 2016), explicit bijections and generating-function identities in partition theory (Nyirenda et al., 2022, Agarwal et al., 2022), and the three-filter elimination plus cyclotomic reduction in the unitary perfect-number problem (Maciejewski, 19 May 2026).

A plausible implication is that the expression “Subbarao–Warren Problem” is best treated as a context-dependent research label. In one context it concerns microcanonical thermodynamics; in another it names a phase-field singular-limit problem; in others it refers to congruence classification, partition identities, or the structure of unitary perfect numbers. Any precise use therefore requires the surrounding domain and the defining equations to be stated explicitly.