X-ray emission by sparks implies bremsstrahlung from a population of energetic electrons, but the details of this process remain a mystery. We present detailed statistical analysis of X-ray spectra detected by multiple detectors during sparks produced by 1 MV negative high-voltage pulses with 1 $\mu$s risetime. With over 900 shots, we statistically analyze the signals, assuming that the distribution of spark X-ray fluence behaves as a power law and that the energy spectrum of X-rays detectable after traversing $\sim$2 m of air and a thin aluminum shield is exponential. We then determine the parameters of those distributions by fitting cumulative distribution functions to the observations. The fit results match the observations very well if the mean of the exponential X-ray energy distribution is 86 $\pm$ 7 keV and the spark X-ray fluence power law distribution has index -1.29 $\pm$ 0.04 and spans at least 3 orders of magnitude in fluence.

We obtain the Maxwell-J\"uttner distribution function at first order in the post-Newtonian approximation within the framework of general relativity. Taking into account the aforesaid distribution function, we compute the particle four-flow and energy-momentum tensor. We focus on the search of static solutions for the gravitational potentials with spherical symmetry. In doing so, we obtain the density, pressure and gravitational potential energy profiles in terms of dimensionless radial coordinate by solving the aforesaid equations numerically. In particular, we find the parametric profile for the equation of state $p/\rho$ in terms of the dimensionless radial coordinate. Due to its physical relevance, we also find the galaxy rotation curves using the post-Newtonian approximation. We join two different kinds of static solutions in order to account for the linear regime near the center and the typical flatten behavior at large radii as well.