In this paper, a detailed study of the multiplicities of charged particles produced in p–p collisions, at √s=13TeV, has been carried out. Analysis of UrQMD generated events along pseudo–rapidity (η), azimuthal angle (ϕ) and η-ϕ phase spaces is made with the help of scaled factorial moment method. The results have also been compared with the UrQMD events with finite impact parameter (b) and minimum bias (MB) conditions. From the intermittency exponent (αq), the anomalous fractal dimension dq is derived and the variations of dq with order q are investigated. The observed intermittent fluctuations are manifested by various parameters such as anomalous fractal dimension dq, degree of multifractality (r), critical exponent (ν), Levy index (μ) and multifractal specific heat. In the framework of Ginzburg–Landau theory, the second-order phase transition in the light of scaled factorial moment method has also been carried out to search for the quark-hadron phase transition. It is observed that the intermittent type of fluctuations are strong enough in η-ϕ space compared to the one-dimensional η space and ϕ space. The data reflects the signature of multifractality in ϕ space and η-ϕ space. Whereas in case of η space, it shows the evidence of monofractality. It is also interesting to note that the values of critical exponents for all impact parameters indicate the apparent existence of quark hadron phase transition except for MB events. However, no such behaviour is observed in ϕ space. From the knowledge of generalized fractal dimension Dq, multifractal specific heat is derived from scale factorial moment analysis in case of one-dimensional η space, ϕ space and two-dimensional η-ϕ space.