In this work, we have analyzed the $\gamma$-ray data in the energy range 100 MeV to 500 GeV from the \emph{Fermi}-Large Area Telescope (LAT) observations for the period August 4, 2008 to March 31, 2019. The $\gamma$-ray emission from AR Scorpii over the last decade is not statistically significant and therefore 2$\sigma$ upper limit on the integral flux above 100 MeV has been estimated. We reproduce the non-thermal broadband spectral energy distribution of AR Scorpii using an emission model having two synchrotron components due to the relativistic electrons in very high magnetic fields. The first component (Synchrotron-1) broadly describes the emissions at radio to high energy X-rays through the synchrotron radiation originating from a spherical region of radius $\sim$ 1.8$\times$10$^{10}$ cm and a magnetic field strength of $\sim$ 10$^3$ Gauss. The second component (Synchrotron-2) which reproduces the X-ray emission at lower energies and predicts the $\gamma$-ray emission, originates from another spherical region with radius $\sim$ 1.4$\times$10$^{10}$ cm and a magnetic field strength of $\sim$ 10$^6$ Gauss. The relativistic electron populations in both the emission regions are described by a smooth broken power law energy distribution. The $\gamma$-ray emission predicted by the Synchrotron-2 model is below the broadband sensitivity of the \emph{Fermi}-LAT and is also consistent with the 95$\%$ confidence level upper limit on the integral flux above 100 MeV derived from more than 10 years of observations. According to our model, the binary system AR Scorpii could be a $\gamma$-ray source, although its emission level must be below the current detection limit of the \emph{Fermi}-LAT.

The discovery of the white dwarf binary system AR Scorpii (AR Sco) with its fascinating non-thermal dominated multi-frequency emission has sparked renewed interest in potential high energy gamma-ray emission from white dwarf pulsars. The Spectral Energy Distribution (SED) below and above optical shows evidence of non-thermal synchrotron emission, with pulsed emission in optical and X-ray bands at the white dwarf spin period ($P_{*} = 117 \, \mbox{s}$) as well as a beat period ($P_{\rm b} = 118.1 \, \mbox{s}$) with the binary period. From an energy perspective, the highly magnetic rotating white dwarf can accelerate particles to TeV energies. In this study, a search for high energy gamma-ray emission was conducted between 100 MeV - 500 GeV by analysing the newly available Fermi-LAT Pass 8 data with the new Fermi 1.0.1 Science Tools. Binned likelihood analysis was done using power law, broken power law and log parabola models. From the selected Region of Interest (ROI) centred on AR Sco's position, we calculated a significance of $\sqrt{TS}\leq$ 3.87$\sigma$ for the integrated gamma-ray activity between 100 MeV - 500 GeV at a photon flux level of 0.486 $\pm$ 0.261 x 10$^{-8}$ photons cm$^{-2}$ s$^{-1}$ using the broken power law model. This resulted in a 3$\sigma$ upper-limit detection from the position of AR Sco. The location of AR Sco inside the Rho Ophiuchi (Rho Oph) molecular cloud complex combined with the poor spatial resolution of Fermi-LAT, complicates any positive identification of low-level gamma-ray activity at the location that coincides with the position of AR Sco.