you don't need to make an atmospheric correction for the MOD13Q1 product because it was generated by the Maximum Value Composite method,this means that among 16 rasters (16 days) the reflectance value for a pixel is the least contaminated by clouds, but also the one with the maximum value.
They provide information on the quality of data VI (vegetations index):
the MOD13Q1 Pixel Reliability, and MOD13Q1 VI Quality rasters provide informations about the quality of data VI (vegetations index):
Refer to TABLE 1: MOD13Q1 Pixel Reliability and TABLE 2: MOD13Q1 VI Quality in the https://lpdaac.usgs.gov/products/modis_products_table/mod13q1
If you use Daily products like MOD09GQ (MODIS Terra Surface Reflectance Bands 1–2 /250m), the atmospheric correction can be logical.