Recently, the covariant formulation of the Tolman-Oppenheimer-Volkoff (TOV) equations for studying the equilibrium structure of a spherically symmetric compact star in the presence of the pressure anisotropy in the interior of a star was presented in Phys. Rev. D \textbf{97} (2018) 124057. It was suggested there that the anisotropic solution of these equations can be obtained by finding, first, the solution of the common TOV equations for the isotropic pressure, and then by solving the differential equation for the anisotropic pressure whose particular form was established on the basis of the covariant TOV equations. It turns out that the anisotropic pressure determined according to this scheme has a nonremovable singularity $\Pi\sim\frac{1}{r^2}$ in the center of a star, and, hence, the corresponding anisotropic solution cannot represent a physically relevant model of an anisotropic compact star. A new scheme for constructing the anisotropic solution, based on the covariant TOV equations, is suggested, which leads to the regularly behaved physical quantities in the interior of a star. A new algorithm is applied to build model anisotropic strange quark stars with the MIT bag model equation of state.

Thermodynamic properties of strange quark matter (SQM) in a nonuniform magnetic field are studied within the phenomenological MIT bag model under the charge neutrality and beta equilibrium conditions, relevant to the interior of strange quark stars. The spatial dependence of the magnetic field strength is modeled by the dependence on the baryon chemical potential in the exponential and power forms. The total energy density, longitudinal and transverse pressures in magnetized SQM are found as functions of the baryon chemical potential. It is clarified that the central magnetic field strength in a strange quark star is bound from above by the critical value at which the derivative of the longitudinal pressure with respect to the baryon chemical potential vanishes first somewhere in the interior of a star under varying the central field. Above this upper bound, the instability along the magnetic field is developed in magnetized SQM. The change in the form of the dependence of the magnetic field strength on the baryon chemical potential between the exponential and power ones has a nonnegligible effect on the critical magnetic field strength while the variation of the bag pressure within the absolute stability window for magnetized SQM has a little effect on the critical field.

Spherically symmetric relativistic stars with the polytropic equation of state, which possess the local pressure anisotropy, are considered in the context of general relativity. The modified Lane-Emden equations are derived for the special ansatz for the anisotropy parameter $\Delta$ in the form of the differential relation between $\Delta$ and the metric function $\nu$. The analytical solutions of the obtained equations are found for incompressible fluid stars. The dynamical stability of incompressible anisotropic fluid stars against radial oscillations is studied.

Magnetized strange quark stars, composed of strange quark matter (SQM) and self-bound by strong interactions, can be formed if the energy per baryon of magnetized SQM is less than that of the most stable $^{56}$Fe nucleus under the zero external pressure and temperature. Utilizing the MIT bag model description of magnetized SQM under charge neutrality and beta equilibrium conditions, the corresponding absolute stability window in the parameter space of the theory is determined. It is shown that there exists the maximum magnetic field strength allowed by the condition of absolute stability of magnetized SQM. The value of this field, $H\sim3\cdot10^{18}$ G, represents the upper bound on the magnetic field strength which can be reached in a strongly magnetized strange quark star.

Spherically symmetric relativistic stars with the polytropic equation of state (EoS), which possess the local pressure anisotropy, are considered within the framework of general relativity. The generalized Lane-Emden equations are derived for the arbitrary anisotropy parameter $\Delta=p_t-p_r$ ($p_t$ and $p_r$ being the transverse and radial pressure, respectively). They are then applied to some special ansatz for the anisotropy parameter in the form of the differential relation between the anisotropy parameter $\Delta$ and the metric function $\nu$. The analytical solutions of the obtained equations are found for incompressible fluid stars and then used for getting their mass-radius relation, gravitational and binding energy. Also, following the Chandrasekhar variational approach, the dynamical stability of incompressible anisotropic fluid stars with the polytropic EoS against radial oscillations is studied. It is shown that the local pressure anisotropy with $p_t>p_r$ can make the incompressible fluid stars unstable with respect to radial oscillations, in contrast to incompressible isotropic fluid stars with the polytropic EoS which are dynamically stable.

The impact of a strong magnetic field, varying with the total baryon number density, on thermodynamic properties of strange quark matter (SQM) in the core of a magnetized hybrid star is considered at zero temperature within the framework of the Massachusetts Institute of Technology (MIT) bag model. It is clarified that the central magnetic field strength is bound from above by the value at which the derivative of the longitudinal pressure with respect to the baryon number density vanishes first somewhere in the quark core under varying the central field. Above this upper bound, the instability along the magnetic field is developed in magnetized SQM. The total energy density, longitudinal and transverse pressures are found as functions of the total baryon number density.

The stability of magnetized strange quark matter (MSQM) is studied in the MIT bag model with the density dependent bag pressure. In the consistent thermodynamic description of MSQM, the quark chemical potentials, the total thermodynamic potential and the anisotropic pressure acquire the corresponding additional term proportional to the density derivative of the bag pressure. The model parameter space is determined, for which MSQM is absolutely stable, i.e., its energy per baryon is less than that of the most stable $^{56}$Fe nucleus under the zero external pressure and vanishing temperature. It is shown that there exists the magnetic field strength $H_{u\,max}$ at which the upper bound $B_\infty^u$ on the asymptotic bag pressure $B_\infty\equiv B(\varrho_B\gg \varrho_0$) ($\varrho_0$ being the nuclear saturation density) from the absolute stability window vanishes. The value of this field, \hbox{$H_{u\,max}\sim$$(1$--$3)\cdot10^{18}$}~G, represents the upper bound on the magnetic field strength, which can be reached in a strongly magnetized strange quark star. It is clarified how the absolute stability window and upper bound on the magnetic field strength are affected by varying the parameters in the Gaussian parametrization for the density dependence of the bag pressure.

Spin polarized states in dense neutron matter with recently developed Skyrme effective interaction (BSk20 parametrization) are considered in the magnetic fields $H$ up to $10^{20}$ G at finite temperature. In a strong magnetic field, the total pressure in neutron matter is anisotropic, and the difference between the pressures parallel and perpendicular to the field direction becomes significant at $H>H_{th}\sim10^{18}$ G. The longitudinal pressure decreases with the magnetic field and vanishes in the critical field $10^{18}<H_c\lesssim10^{19}$ G, resulting in the longitudinal instability of neutron matter. With increasing the temperature, the threshold $H_{th}$ and critical $H_c$ magnetic fields also increase. The appearance of the longitudinal instability prevents the formation of a fully spin polarized state in neutron matter and only the states with moderate spin polarization are accessible. The anisotropic equation of state is determined at densities and temperatures relevant for the interiors of magnetars. The entropy of strongly magnetized neutron matter turns out to be larger than the entropy of the nonpolarized matter. This is caused by some specific details in the dependence of the entropy on the effective masses of neutrons with spin up and spin down in a polarized state.

Thermodynamic properties of strange quark matter in strong magnetic fields $H$ up to $10^{20}$ G are considered within the MIT bag model at zero temperature implying the constraints of total baryon number conservation, charge neutrality and chemical equilibrium. The pressure anisotropy, exhibiting in the difference between the pressures along and perpendicular to the field direction, becomes essential at $H>H_{th}$, with the estimate $10^{17}<H_{th}\lesssim10^{18}$ G. The longitudinal pressure vanishes in the critical field $H_c$, which can be somewhat less or larger than $10^{18}$ G, depending on the total baryon number density and bag pressure. As a result, the longitudinal instability occurs in strange quark matter, which precludes: (1) a significant drop in the content of $s$ quarks, which, otherwise, could happen at $H\sim10^{20}$ G; (2) the appearance of positrons in weak processes in a narrow interval near $H\sim2\cdot10^{19}$ G (replacing electrons). The occurrence of the longitudinal instability leaves the possibility only for electrons to reach a fully polarized state, while for all quark flavors the polarization remains mild even for the fields near $H_c$. The anisotropic equation of state is determined under the conditions relevant to the interiors of magnetars.

Spin polarized states in dense neutron matter with BSk20 Skyrme force are considered in magnetic fields up to $10^{20}$ G. It is shown that the appearance of the longitudinal instability in a strong magnetic field prevents the formation of a fully spin polarized state in neutron matter, and only the states with moderate spin polarization can be developed.

In dense neutron matter under the presence of a strong magnetic field, considered in the model with the Skyrme effective interaction, there are possible two types of spin ordered states. In one of them the majority of neutron spins are aligned opposite to magnetic field (thermodynamically preferable state), and in other one the majority of spins are aligned along the field (metastable state). The equation of state, incompressibility modulus and velocity of sound are determined in each case with the aim to find the peculiarities allowing to distinguish between two spin ordered phases.

The possibility of the appearance of the states with a nonzero average helicity in neutron matter is studied in the model with the Skyrme effective interaction. By providing the analysis of the self-consistent equations at zero temperature, it is shown that neutron matter with the Skyrme BSk18 effective force undergoes at high densities a phase transition to the state in which the degeneracy with respect to helicity of neutrons is spontaneously removed.

Spin polarization of neutron matter at finite temperatures and strong magnetic fields up to $10^{18}$ G is studied in the model with the Skyrme effective interaction. It is shown that, together with the thermodynamically stable branch of solutions for the spin polarization parameter corresponding to the case when the majority of neutron spins are oriented opposite to the direction of the magnetic field (negative spin polarization), the self-consistent equations, beginning from some threshold density, have also two other branches of solutions corresponding to positive spin polarization. The influence of finite temperatures on spin polarization remains moderate in the Skyrme model up to temperatures relevant for protoneutron stars, and, in particular, the scenario with the metastable state characterized by positive spin polarization, considered at zero temperature in Phys. Rev. C {\bf 80}, 065801 (2009), is preserved at finite temperatures as well. It is shown that above certain density the entropy for various branches of spin polarization in neutron matter with the Skyrme interaction in a strong magnetic field demonstrates the unusual behavior being larger than that of the nonpolarized state. By providing the corresponding low-temperature analysis, it is clarified that this unexpected behavior should be addressed to the dependence of the entropy of a spin polarized state on the effective masses of neutrons with spin up and spin down, and to a certain constraint on them which is violated in the respective density range.

Spin polarized states in neutron matter at a strong magnetic field are considered in the model with the Skyrme effective interaction (SLy4, SLy7 parametrizations). Analyzing the self-consistent equations at zero temperature, it is shown that a thermodynamically stable branch of solutions for the spin polarization parameter as a function of density corresponds to the negative spin polarization when the majority of neutron spins are oriented oppositely to the direction of the magnetic field. Besides, beginning from some threshold density being dependent on the magnetic field strength the self-consistent equations have also two other branches (upper and lower) of solutions for the spin polarization parameter with the positive spin polarization. The free energy corresponding to the upper branch turns out to be very close to the free energy corresponding to the thermodynamically preferable branch with the negative spin polarization. As a consequence, at a strong magnetic field, the state with the positive spin polarization can be realized as a metastable state at the high density region in neutron matter which under decreasing density at some threshold density changes into a thermodynamically stable state with the negative spin polarization. The calculations of the neutron spin polarization parameter and energy per neutron as functions of the magnetic field strength show that the influence of the magnetic field remains small at the field strengths up to $10^{17}$ G.

Spin polarized states in nuclear matter with an effective nucleon-nucleon interaction are studied for a wide range of isospin asymmetries and densities. Based on a Fermi liquid theory, it is shown that there are a few possible scenarios of spin ordered phase transitions: (a) nuclear matter undergoes at some critical density a phase transition to a spin polarized state with the oppositely directed spins of neutrons and protons (Skyrme SLy4 and Gogny D1S interactions); (b) at some critical density, a spin polarized state with the like-directed neutron and proton spins appears (Skyrme SkI5 interaction); (c) nuclear matter under increasing density, at first, undergoes a phase transition to the state with the opposite directions of neutron and proton spins, which goes over at larger density to the state with the same direction of nucleon spins (Skyrme SkI3 interaction).

The possibility of appearance of spin polarized states in nuclear matter is studied within the framework of a Fermi liquid theory with Skyrme effective forces in a wide range of isospin asymmetries and densities. There are a few possible scenarios of spin ordered phase transitions: (a) nuclear matter with SLy4 interaction undergoes at some critical density a phase transition to a spin polarized state with the oppositely directed spins of neutrons and protons; (b) in nuclear matter with SkI5 interaction a spin polarized state with the like-directed neutron and proton spins is formed; (c) nuclear matter with SkI3 interaction under increasing density, at first, undergoes a phase transition to the state with the opposite directions of neutron and proton spins, which goes over at larger density to the state with the same direction of nucleon spins. Spin polarized states at strong excess neutrons over protons are accompanied by the long tails in the density profiles of the neutron spin polarization parameter near the critical density, if the energy gain of the transition from the nonpolarized state to the polarized one is the decreasing function of isospin asymmetry (SLy4 force). In the opposite case, if the energy gain is increased with isospin asymmetry, there are no long tails in the density distribution of the neutron spin polarization parameter (SkI3, SkI5 forces).

The possibility of appearance of spin polarized states in symmetric and strongly asymmetric nuclear matter is analyzed within the framework of a Fermi liquid theory with the Skyrme effective interaction. The zero temperature dependence of the neutron and proton spin polarization parameters as functions of density is found for SkM$^*$, SGII (symmetric case) and SLy4, SLy5 (strongly asymmetric case) effective forces. By comparing free energy densities, it is shown that in symmetric nuclear matter ferromagnetic spin state (parallel orientation of neutron and proton spins) is more preferable than antiferromagnetic one (antiparallel orientation of spins). Strongly asymmetric nuclear matter undergoes at some critical density a phase transition to the state with the oppositely directed spins of neutrons and protons while the state with the same direction of spins does not appear. In comparison with neutron matter, even small admixture of protons strongly decreases the threshold density of spin instability. It is clarified that protons become totally polarized within a very narrow density domain while the density profile of the neutron spin polarization parameter is characterized by the appearance of long tails near the transition density.