In order to understand how the doping with self-assembled nanorods of different sizes and concentrations as well as applied magnetic fields affect the critical current anisotropy in YBa 2 cu 3 o 7−x (YBCO) thin films close to YBCO c-axis, we present an extensive and systematic computational study done by molecular dynamics simulation. the simulations are also used to understand experimentally measured J c (θ) curves for BaHfO 3 , BaZrO 3 and BaSno 3 doped YBCO thin films with the help of nanorod parameters obtained from transmission electron microscopy measurements. our simulations reveal that the relation between applied and matching field plays a crucial role in the formation of J c (θ)-peak around YBCO c-axis (c-peak) due to vortex-vortex interactions. We also find how different concentrations of different size nanorods effect the shape of the c-peak and explain how different features, such as double c-peak structures, arise. In addition to this, we have quantitatively explained that, even in an ideal superconductor, the overdoping of nanorods results in decrease of the critical current. our results can be widely used to understand and predict the critical current anisotropy of YBco thin films to improve and develop new pinscapes for various transport applications.High temperature superconductors (HTS) are expected to have large number of applications in different fields of technology and power industry in the future 1-3 . Since all known HTS are of type II, the critical current passed through them is highly dependent on the surrounding magnetic field due to the movement of vortices. Thus, to enhance and widen the usability of HTS, the dynamics of vortices need to be well understood.Among the high temperature superconductors, YBa 2 Cu 3 O 7−x (YBCO) seems the most practical choice when thinking for the applications 1 . The intrinsic anisotropy of the critical current, in thin films and coated conductors, can be modified by adding impurities within the lattice of YBCO which pin the vortices restricting their movement. Based on growth conditions and lattice mismatch between the YBCO and the dopant as well as their elastic properties 4,5 , impurities such as Y 2 O 3 6 , BaCeO 3 7-9 and BaZrO 3 (BZO) 10,11 can form uncorrelated randomly distributed nanoparticles within the YBCO lattice. Under optimized deposition conditions, via a spontaneous phase-separation and strain-driven self-assembly process during film deposition 12 , self-assembly of nanorods of BaHfO 3 (BHO) 1 , BaZrO 3 (BZO) 4,13,14 , BaSnO 3 (BSO) 15,16 , Ba 2 YTaO 3 (BYTO) 17 or Ba 2 YNbO 6 (BYNO) 18 within the YBCO lattice can be realized.Recently, a topic of interest has been to add both point-like nanodots and nanorods within the YBCO lattice simultaneously. This has been achieved by doping YBCO simultaneously with both BYTO and BYNO (referred as BYNTO) with an additional rare earth oxide, leading to continuous niobiate/tantalate nanorods and rare-earth oxide nanoparticles 19 . A lot of experimental research has been done in order to understand the...