We analyze different stages of magnetized quark star evolution incorporating baryon number conservation and using an anisotropic energy momentum tensor. The first stages of the evolution are simulated through the inclusion of trapped neutrinos and fixed entropy per particle, while in the last stage the star is taken to be deleptonized and cold. We find that, although strong magnetic fields modify quark star masses, the evolution of isolated stars needs to be constrained by fixed baryon number, which necessarily lowers the possible star masses. Moreover, magnetic field effects, measured by the difference between the parallel and perpendicular pressures, are more pronounced in the beginning of the star evolution, when there is a larger number of charged leptons and up quarks. We also show that having a spatially varying magnetic field allows for larger magnetic fields to be supported.

We calculate the number density, energy density, transverse pressure, longitudinal pressure, and magnetization of an ensemble of spin one-half particles in the presence of a homogenous background magnetic field. The magnetic field direction breaks spherical symmetry causing the pressure transverse to the magnetic field direction to be different than the pressure parallel to it. We present explicit formulae appropriate at zero and finite temperature for both charged and uncharged particles including the effect of the anomalous magnetic moment. We demonstrate that the resulting expressions satisfy the canonical relations, Omega = - P_parallel and P_perp = P_parallel - M B, with M = - d Omega/d B being the magnetization of the system. We numerically calculate the resulting pressure anisotropy for a gas of protons and a gas of neutrons and demonstrate that the inclusion of the anomalous magnetic increases the level of pressure anisotropy in both cases.

We calculate the equation of state for degenerate quark matter to leading order in hard-dense-loop (HDL) perturbation theory. We solve the Tolman-Oppenheimer-Volkov equations to obtain the mass-radius relation for dense quark stars. Both the perturbative QCD and the HDL equations of state have a large variation with respect to the renormalization scale for quark chemical potential below 1 GeV which leads to large theoretical uncertainties in the quark star mass-radius relation.