Galaxy clusters are known to be reservoirs of Cosmic Rays (CRs), as inferred from theoretical calculations or detection of CR-derived observables. CR acceleration in clusters is mostly attributed to the dynamical activity that produces shocks. Shocks in clusters emerge out of merger or accretion, but which one is more effective in producing CRs? at which dynamical phase? and why? To this aim, we study the production or injection of CRs through shocks and its evolution in the galaxy clusters using cosmological simulations with the {\sc enzo} code. Particle acceleration model considered here is primarily the Diffusive Shock Acceleration (DSA) of thermal particles, but we also report a tentative study with pre-existing CRs. Defining appropriate dynamical states using the concept of virialization, we studied a sample of merging and non-merging clusters. We report that the merger shocks (with Mach number $\mathcal{M}\sim2-5$) are the most effective CR producers, while high-Mach peripheral shocks (i.e. $\mathcal{M}>5$) are mainly responsible for the brightest phase of CR injection in clusters. Clusters once merged, permanently deviate from CR and X-ray mass scaling of non-merging systems, enabling us to use it as a tool to determine the state of merger. Through a temporal and spatial evolution study, we found a strong correlation between cluster merger dynamics and CR injection. We observed that the brightest phase of X-ray and CR injection from clusters occur respectively at about 1.0 and 1.5 Gyr after every mergers, and CR injection peaks near to the cluster virial radius (i.e $r_{200}$). Delayed CR injection peaks found in this study deserve further investigation for possible impact on the evolution of CR-derived observables from galaxy clusters.

We present the discovery of ringlike diffuse radio emission structures in the peripheral regions of the Bullet cluster 1E 0657$-$55.8. Ring formations are spanning between 1--3 Mpc away from the center of the cluster, significantly further away from the two already reported relics. Integrated fluxes of four of the sub-regions in the inner `ring' from 4.5 to 10 GHz have also been reported. To understand the possible origin of these structures, here we present a maiden attempt of numerical modelling of a 3D and realistic `bullet' like event in a full cosmological ($\Lambda$CDM) environment with N-body plus hydrodynamics code. We report a simulated `bullet' found inside a (128 Mpc)$^3$ volume simulation with a speed of 2700 km s$^{-1}$, creating a high supersonic bow shock of Mach $M=3.5$ and a clear evidence of temporal separation of dark matter and baryons, assuring no challenge to $\Lambda$CDM cosmology from the bullet event as of now. We are also able to unveil the physics behind the formation of these observed multiple shock structures. Modelled radio emissions in our simulation support a complex combination of merger-associated processes that accelerates and re-accelerates fossil and cosmic-ray electrons. With a time evolution study and the computed radio emissions, we have shown that the ring like formation around the bullet is originated due to the interaction of the strong merger shocks with the accretion shocks at the periphery. The multiple shock structures observed are possibly originated from multiple mergers that have taken place at different times and much before the bullet event.

Observed steep mass scaling of radio power from the available high mass clusters has ruled out the prospect of detection of 'galaxy groups'. But, the available simulations and observations of thermal emissions show that the groups are merger prone, thus non-virialised, indicating better visibility in the radio waves. Detection of radio emissions from them would help us to understand the scale-dependent particle acceleration mechanisms also groups can be a unique laboratory to test the models of cosmic magnetism and canbe the potential source of WHIMs. So, we have modelled radio emissions from the simulated structures using {\sc{ENZO}}. We present a model for computing magnetic field and for the first time, used the electron energy spectrum from both the Fermi I (DSA) and Fermi II (TRA) mechanisms to compute radio emissions. Computed radio power from more than 200 simulated objects, mass ranging $\geq 10^{13}$ to $2\times 10^{15} M_{\odot}$ show a new mass scaling of $M_{500} \propto P_{1.4\;GHz}^{2.17\pm 0.08}$ and a strong correlation scale of $L_X \propto P_{1.4\;GHz}^{1.08\pm 0.05}$. Both magnetic field and radio power are shown to have adequately replicated the available observations at high mass, allowing us to extend the results to further smaller masses. We report that groups below $10^{14}\;\rm{M_{\odot}}$ show the existence of 10s of nano to a sub-$\mu$G magnetic field and about 10$^{19-23}$ W Hz$^{-1}$ of radio power, much higher than what existing mass scaling predicts. We found that the combined radio power from TRA and DSA electrons can only fit very well to all the observed `radio halos'. Finally, we have implemented this model on a real data set obtained from the Sloan Digital Sky Survey (SDSS). It predicts about 10s to 100s $\mu$Jy/(10$\arcsec$ beam) of radio flux in groups indicating their detectability with existing and aplenty with the future radio telescopes.

`Galaxy groups' have hardly been realised as a separate class of objects with specific characteristics in the structural hierarchy. The presumption that the self-similarity of dark matter structures is a valid prescription for the baryonic universe at all scales has rendered smaller structures undetectable by current observational facilities, leading to lesser dedicated studies on them. Some recent reports that indicate a deviation from $\rm{L_x}$-T scaling in groups compared to clusters have motivated us to study their physical properties in depth. In this article, we report the extensive study on physical properties of groups in comparison to the clusters through cosmological hydrodynamic plus N-body simulations using ENZO 2.2 code. As additional physics, radiative cooling, heating due to supernova and star motions, star formation and stellar feedback has been implemented. We have produced a mock sample of 362 objects with mass ranging from $5\times10^{12}\; \rm{M_{\odot}}$ to 2.5$\times 10^{15}\; \rm{M_{\odot}}$. Strikingly, we have found that objects with mass below $\sim$ $8\times 10^{13}\;\rm{M_{\odot}}$ do not follow any of the cluster self-similar laws in hydrostatics, not even in thermal and non-thermal energies. Two distinct scaling laws are observed to be followed with breaks at $\sim$ $8\times 10^{13}\;\rm{M_{\odot}}$ for mass, $\sim$1 keV for temperature and $\sim$1 Mpc for radius. This places groups as a distinct entity in the hierarchical structures, well demarcated from clusters. This study reveals that groups are mostly far away from virialization, suggesting the need for formulating new models for deciphering their physical parameters. They are also shown to have high turbulence and more non-thermal energy stored, indicating better visibility in the non-thermal regime.