We construct and analyze a novel family of exact higher-dimensional black hole solutions in Einstein-Power-Yang-Mills gravity minimally coupled to a scalar field multiplet supporting a global monopole and a surrounding anisotropic ''scalar gas'' whose stress-energy obeys a generalized Murnaghan equation of state. By adopting the Wu-Yang magnetic ansatz for the non-Abelian sector and enforcing the physically motivated radial condition pᵣ=-ρ, the field equations admit closed-form expressions for the matter density and a single-function lapse F (r) expressed in terms of elementary and Gauss hypergeometric functions. The Murnaghan fluid interpolates between a non-linear core and an effective vacuum at a large radius, producing a backreaction encoded by parameters that control the stiffness and scalar-backreaction scale. We perform a systematic survey of classical energy condition, identifying regions of parameter space where the null energy condition/dominant energy condition holds while the strong energy condition is generically violated (phantom-like behavior), and examine curvature invariants to demonstrate that the solutions possess a central curvature singularity for n>3. Thermodynamic properties are derived in the extended phase space: explicit formulas for the Hawking temperature, entropy, conjugate potentials, and a generalized Smarr relation are obtained; the heat capacity exhibits divergencies and sign changes that mark local stability boundaries and second-order phase transitions. Varying the Yang-Mills charge, the nonlinearity index, the monopole coupling, and the Murnaghan parameters generates a rich phase behavior structure and topology changes in the defect (φ-Duan) map. The results underscore the manner in which gauge nonlinearity, topological defects, and dual-polytropic matter collaboratively transform BH thermodynamics and the topology of phase spaces.
Sekhmani et al. (Sat,) studied this question.