The clustering of local extrema including peaks and troughs will be exploited to assess Gaussianity, asymmetry and the footprint of cosmic strings network on the CMB random field observed by {\it Planck} satellite. The number density of local extrema reveals some non-resolved shot noise in {\it Planck} maps. The \texttt{SEVEM} data has maximum number density of peaks, $n_{pk}$, and troughs, $n_{tr}$, compared to other observed maps. The cumulative of $n_{pk}$ and $n_{tr}$ above and below a threshold, $\vartheta$, for all {\it Planck} maps except for the 100GHz band are compatible with the Gaussian random field. The unweighted Two-Point Correlation Function (TPCF), $\Psi(\theta;\vartheta)$, of the local extrema illustrates significant non-Gaussianity for angular separation $\theta\le 15'$ for all available thresholds. Our results show that to put the feasible constraint on the amplitude of the mass function based on the value of $\Psi$ around the {\it Doppler peak} ($\theta\approx 70'-75'$), we should consider $\vartheta\gtrsim+1.0$. The scale independent bias factors for peak above a threshold for large separation angle and high threshold level are in agreement with that of expected for a pure Gaussian CMB. Unweighted TPCF of local extrema demonstrates a level of rejecting Gaussian hypothesis in \texttt{SMICA}. Genus topology also confirms the Gaussian hypothesis for different component separation maps. Tessellating CMB map with disk of size $6^{\circ}$ based on $n_{pk}$ and $\Psi_{pk-pk}$ demonstrate statistical symmetry in {\it Planck} maps. Combining all maps and applying the $\Psi_{pk-pk}$ puts the upper bound on the cosmic string's tension: $G\mu^{(up)} \lesssim 5.00\times 10^{-7}$.

The Cold Spot is a statistically significant anomaly in the Cosmic Microwave Background (CMB) sky. In this work we assess whether a huge void or a cosmic texture could have produced such an anomaly through searching for their gravitational redshift (Rees-Sciama effect) and lensing signatures on the {\it Planck} CMB sky. In the flat sky approximation, we find the amplitudes for the corresponding templates for both candidates consistent with zero, leaving little room, according to {\it Planck}, for those structures as the main source of the Cold Spot anomaly.

Accurate and precise measurements of masses of galaxy clusters are key to derive robust constraints on cosmological parameters. Rising evidence from observations, however, confirms that X-ray masses, obtained under the assumption of hydrostatic equilibrium, might be underestimated, as previously predicted by cosmological simulations. We analyse more than 300 simulated massive clusters, from `The Three Hundred Project', and investigate the connection between mass bias and several diagnostics extracted from synthetic X-ray images of these simulated clusters. We find that the azimuthal scatter measured in 12 sectors of the X-ray flux maps is a statistically significant indication of the presence of an intrinsic (i.e. 3D) clumpy gas distribution. We verify that a robust correction to the hydrostatic mass bias can be inferred when estimates of the gas inhomogeneity from X-ray maps (such as the azimuthal scatter or the gas ellipticity) are combined with the asymptotic external slope of the gas density or pressure profiles, which can be respectively derived from X-ray and millimetric (Sunyaev-Zeldovich effect) observations. We also obtain that mass measurements based on either gas density and temperature or gas density and pressure result in similar distributions of the mass bias. In both cases, we provide corrections that help reduce both the dispersion and skewness of the mass bias distribution. These are effective even when irregular clusters are included leading to interesting implications for the modelling and correction of hydrostatic mass bias in cosmological analyses of current and future X-ray and SZ cluster surveys.

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}$.

We examine linear perturbation theory to evaluate the contribution of viscosity coefficient in the growing of dark matter perturbations in the context of the bulk viscous dark energy model inspired by thermodynamical dissipative phenomena proposed by [B. Mostaghel et al., 2017, The European Physical Journal C, 77, 541]. As the cosmological implementations, we investigate the {\it Integrated Sachs-Wolfe} (ISW) auto-power spectrum, the ISW-galaxy cross-power spectrum and derive limits on $f\sigma_8$. The dimensionless bulk viscosity coefficient ($\gamma$) in the absence of interaction between dark sectors, modifies the Hubble parameter and the growth function, while the Poisson equation remains unchanged. Increasing $\gamma$ reduces the dark matter growing mode at the early epoch while a considerable enhancement will be achieved at the late time. This behavior imposes non-monotonic variation in the time evolution of gravitational potential generating a unique signature on the CMB photons. The bulk viscous dark energy model leads to almost a decreasing in ISW source function at the late time. Implementation of the Redshift Space Distortion (RSD) observations based on "Gold-2017" catalogue, shows $\Omega^0_{\rm m} =0.303^{+0.044}_{-0.038}$, $\gamma=0.033^{+0.098}_{-0.033}$ and $\sigma_8= 0.769^{+0.080}_{-0.089}$ at $1\sigma$ level of confidence. Finally, tension in the $\sigma_8$ is alleviated in our viscous dark energy model.

We explore the use of random forest and gradient boosting, two powerful tree-based machine learning algorithms, for the detection of cosmic strings in maps of the cosmic microwave background (CMB), through their unique Gott-Kaiser-Stebbins effect on the temperature anisotropies.The information in the maps is compressed into feature vectors before being passed to the learning units. The feature vectors contain various statistical measures of processed CMB maps that boost the cosmic string detectability. Our proposed classifiers, after training, give results improved over or similar to the claimed detectability levels of the existing methods for string tension, $G\mu$. They can make $3\sigma$ detection of strings with $G\mu \gtrsim 2.1\times 10^{-10}$ for noise-free, $0.9'$-resolution CMB observations. The minimum detectable tension increases to $G\mu \gtrsim 3.0\times 10^{-8}$ for a more realistic, CMB S4-like (II) strategy, still a significant improvement over the previous results.

We propose a multi-scale edge-detection algorithm to search for the Gott-Kaiser-Stebbins imprints of a cosmic string (CS) network on the Cosmic Microwave Background (CMB) anisotropies. Curvelet decomposition and extended Canny algorithm are used to enhance the string detectability. Various statistical tools are then applied to quantify the deviation of CMB maps having a cosmic string contribution with respect to pure Gaussian anisotropies of inflationary origin. These statistical measures include the one-point probability density function, the weighted two-point correlation function (TPCF) of the anisotropies, the unweighted TPCF of the peaks and of the up-crossing map, as well as their cross-correlation. We use this algorithm on a hundred of simulated Nambu-Goto CMB flat sky maps, covering approximately $10\%$ of the sky, and for different string tensions $G\mu$. On noiseless sky maps with an angular resolution of $0.9'$, we show that our pipeline detects CSs with $G\mu$ as low as $G\mu\gtrsim 4.3\times 10^{-10}$. At the same resolution, but with a noise level typical to a CMB-S4 phase II experiment, the detection threshold would be to $G\mu\gtrsim 1.2 \times 10^{-7}$.

In this paper we investigate ghost dark energy model in the presence of non-linear interaction between dark energy and dark matter. We also extend the analysis to the so called generalized ghost dark energy (GGDE) which $\rho_D=\alpha H+\beta H^2$. The model contains three free parameters as $\Omega_D, \zeta(=\frac{8\pi G \beta}{3})$ and $b^2$ (the coupling coefficient of interactions). We propose three kinds of non-linear interaction terms and discuss the behavior of equation of state, deceleration and dark energy density parameters of the model. We also find the squared sound speed and search for signs of stability of the model. To compare the interacting GGDE model with observational data sets, we use more recent observational outcomes, namely SNIa from JLA catalog, Hubble parameter, baryonic acoustic oscillation and the most relevant CMB parameters including, the position of acoustic peaks, shift parameters and redshift to recombination. For GGDE with the first non-linear interaction, the joint analysis indicates that $\Omega_D=0.7192\pm0.0062$, $b^2=0.146^{+0.030}_{-0.026}$ and $\zeta=0.104\pm0.047$ at 1 optimal variance error. For the second interaction, the best fit values at $1\sigma$ confidence are $\Omega_D=0.72091\pm0.0065$, $b^2=0.0395\pm0.0080$ and $\zeta\le0.0173$. According to combination of all observational data sets considered in this paper the best fit values for third non-linearly interacting model are $\Omega_D=0.7287\pm0.0062$, $b^2=0.0109\pm0.0023$ and $\zeta\le0.00764$ at $1\sigma$ confidence interval. Finally we found that the presence of interaction is compatible in mentioned models via current observational data sets.

Inspired by thermodynamical dissipative phenomena, we consider bulk viscosity for dark fluid in a spatially flat two-component Universe. Our viscous dark energy model represents Phantom crossing avoiding Big-Rip singularity. We propose a non-minimal derivative coupling scalar field with zero potential leading to accelerated expansion of Universe in the framework of bulk viscous dark energy model. In this approach, coupling constant ($\kappa$) is related to viscosity coefficient ($\gamma$) and energy density of dark energy at the present time ($\Omega_{\rm DE}^0$). This coupling is bounded as $\kappa\in [-1/9H_0^2(1-\Omega_{\rm DE}^0), 0]$ and for $\gamma=0$ leads to $\kappa=0$. To perform robust analysis, we implement recent observational data sets including Joint Light-curve Analysis (JLA) for SNIa, Gamma Ray Bursts (GRBs) for most luminous astrophysical objects at high redshifts, Baryon Acoustic Oscillations (BAO) from different surveys, Hubble parameter from HST project, {\it Planck} data for CMB power spectrum and CMB Lensing. Joint analysis of JLA$+$GRBs$+$BAO$+$HST shows that $\Omega_{\rm DE}^0=0.696\pm 0.010$, $\gamma=0.1404\pm0.0014$ and $H_0=68.1\pm1.3$ at $1\sigma$ confidence interval. {\it Planck} TT observation provides $\gamma=0.32^{+0.31}_{-0.26}$ at $68\%$ confidence limit for viscosity coefficient. Tension in Hubble parameter is alleviated in this model. Cosmographic distance ratio indicates that current observed data prefer to increase bulk viscosity. Finally, the competition between Phantom and Quintessence behavior of viscous dark energy model can accommodate cosmological old objects reported as a sign of age crisis in $\Lambda$CDM model.

In this paper, we propose crossing statistics and its generalization, as a new framework to characterize the anisotropy in a 2D field, e.g. height on a surface, extendable to higher dimensions. By measuring $\nu^+$, the number of up-crossing (crossing points with positive slope at a given threshold of height ($\alpha$)), and $N_{tot}$ (the generalized roughness function), it is possible to distinguish the nature of anisotropy, rotational invariance and Gaussianity of any given surface. For the case of anisotropic correlated self- or multi-affine surfaces (even with different correlation lengths in various directions and/or directional scaling exponents), we analytically derive some relations between $\nu^+$ and $N_{tot}$ with corresponding scaling parameters. The method systematically distinguishes the directions of anisotropy, at $3\sigma$ confidence interval using P-value statistics. After applying a typical method in determining the corresponding scaling exponents in identified anisotropic directions, we are able to determine the kind and ratio of correlation length anisotropy. To demonstrate capability and accuracy of the method, as well validity of analytical relations, our proposed measures are calculated on synthetic stochastic rough interfaces and rough interfaces generated from simulation of ion etching. There are good consistencies between analytical and numerical computations. The proposed algorithm can be mounted with a simple software on various instruments for surface analysis and characterization, such as AFM, STM and etc.

We consider models of inflation with U(1) gauge fields and charged scalar fields including symmetry breaking potential, chaotic inflation and hybrid inflation. We show that there exist attractor solutions where the anisotropies produced during inflation becomes comparable to the slow-roll parameters. In the models where the inflaton field is a charged scalar field the gauge field becomes highly oscillatory at the end of inflation ending inflation quickly. Furthermore, in charged hybrid inflation the onset of waterfall phase transition at the end of inflation is affected significantly by the evolution of the background gauge field. Rapid oscillations of the gauge field and its coupling to inflaton can have interesting effects on preheating and non-Gaussianities.

We consider the Weyl-Dirac theory within the framework of the weak field approximation and show that the resulting gravitational potential differs from that of Newtonian by a repulsive correction term increasing with distance. The scale of the correction term appears to be determined by the time variation rate of the gravitational coupling. It is shown that if the time variation rate of gravitational coupling is adopted from observational bounds, the theory can explain the rotation curves of typical spiral galaxies without resorting to dark matter. To check the consistency of our theoretical model with observation we use Likelihood analysis to find the best-fit values for the free parameters. The mean value for the most important free parameter, $\beta\times 10^{14} (1/yr)$, using the Top-Hat and Gaussian priors are $6.38^{+2.44}_{-3.46}_{-6.71}^{+6.18}$ and $5.72_{-1.18}^{+1.22}_{-2.69}^{+2.90}$, respectively. Although the interval for which $\beta$ is defined is wide, our results show that the goodness of the fit is, by and large, not sensitive to this quantity. The intergalactic effects and gravitational lensing of clusters of galaxies are estimated and seem to be consistent with observational data.