We introduce a pipeline including multifractal detrended cross-correlation analysis (MF-DXA) modified by either singular value decomposition or the adaptive method to examine the statistical properties of the pulsar timing residual ($PTR$) induced by a gravitational wave (GW) signal. We propose a new algorithm, the so-called irregular-MF-DXA, to deal with irregular data sampling. Inspired by the quadrupolar nature of the spatial cross-correlation function of a gravitational wave background, a new cross-correlation function, $\bar{\sigma}_{\times}$, derived from irregular-MF-DXA will be introduced. We show that, this measure reveals the quadrupolar signature in the $PTRs$ induced by stochastic GWB. We propose four strategies based on the $y$-intercept of fluctuation functions, the generalized Hurst exponent, and the width of the singularity spectrum to determine the dimensionless amplitude and power-law exponent of the characteristic strain spectrum as $\mathcal{H}_c(f)\sim\mathcal{A}_{yr}(f/f_{yr})^{\zeta}$ for stochastic GWB. Using the value of Hurst exponent, one can clarify the type of GWs. We apply our pipeline to explore 20 millisecond pulsars observed by Parkes Pulsar Timing Array. The computed scaling exponents confirm that all data are classified into a nonstationary class implying the universality feature. The value of the Hurst exponent is in the range $H\in [0.56,0.87]$. The $q$-dependency of the generalized Hurst exponent demonstrates that the observed $PTRs$ have multifractal behavior, and the source of this multifractality is mainly attributed to the correlation of data which is another universality of the observed datasets. Multifractal analysis of available $PTRs$ datasets reveals an upper bound on the dimensionless amplitude of the GWB, $\mathcal{A}_{yr}< 2.0\times 10^{-15}$.

Using Detrended Fluctuation Analysis (DFA) and box counting method, we test spacial correlation and fractality of Polarization Pulse Profiles (PPPs) of 24 millisecond pulsars (MSPs) which were observed in Parkes Pulsar Timing Array (PPTA) project. DFA analysis indicates that MSPs' PPPs are persistent and the results of box counting method confirm the fractality in the majority of PPPs. A Kolmogorov-Smirnov test indicates that isolated MSPs have more complex PPPs than binary ones. Then we apply our analysis on a random sample of normal pulsars. Comparing the results of our analysis on MSPs and normal pulsars shows that MSPs have more complex PPPs which is resulted from smaller angular half-width of the emission cone and more peaks in MSPs PPPs. On the other hand, high values of Hurst exponent in MSPs confirm compact emission regions in these pulsars.

In this paper for studying the anisotropic strange quark stars, we assume that the radial pressure inside the anisotropic star is a superposition of pressure in an isotropic case plus a Gaussian perturbation term. Considering a proportionality between electric charge density and the density of matter, we solve the TOV equation for different cases numerically. Our results indicate that anisotropy increases the maximum mass $M_{max}$ and also its corresponding radius $R$ for a typical strange quark star. According to our calculations, an anisotropy amplitude of $A=3\times10^{33}Nm^{-2}$ with a standard deviation of $\sigma=3\times10^{3}m$ leads to a neutron star of 1.97$M_{\odot}$. Furthermore, electric charge not only increases the maximum mass and its corresponding radius, but also raises up the anisotropy factor. We can see that the tangential pressure $p_{t}$ and anisotropy factor $\Delta$ unlike the radial pressure $p_{r}$ have a maximum on the surface and this maximum increases by adding electric charge effect. However, we show that anisotropy can be more effective than electric charge in rasing maximum mass of strange quark stars.