By using a first-principles pseudopotential method based on the density functional theory and Vienna ab initio Simulation Package (VASP), we investigate the multiple trapping of C by Ni vacancy (VNi) and its temperature effects in NiAl intermetallics. A single C atom is energetically and favorably situated at the Ni-rich octahedron interstitial site that surrounds Ni vacancy, which is shown via calculating the formation energy of C atom in NiAl with Ni vacancy system. Single C atom prefers to interact with neighboring Ni atom and Al atom to form a covalent bond. In NiAl intermetallics, C atoms prefer to be trapped in the Ni vacancy in the sequential way, thus easily forming the CnVNi (n=1, 2, 3, 4) clusters, in which the C4VNi clusters are most energetically favorable. It is interesting to find that when C atoms are trapped by Ni vacancy, all the C atoms themselves prefer to be combined with each other to form a bond, surrounding Ni vacancy. With the C atoms further added, both the charge density and the deformation charge prefer to bind with each other despite the Ni or Al environment and the intrinsic bonding properties of CC bond contain obvious covalent characteristics. Furthermore, using first-principles calculations combined with statistical model, we quantitatively predict point defect concentration as a function of temperature in NiAl intermetallics. It is concluded that the concentration of intrinsic Ni vacancies (VNi) will obviously increase as temperature increases. With the increase of temperature, the concentration of C atoms in the CnVNi cluster is higher than that at the intrinsic position. Besides, it indicates that most of C atoms in NiAl intermetallics are trapped by Ni vacancy, which is due to the larger binding energy of the CnVNi clusters and most of the C atoms are trapped directly by vacancies at room temperature or high temperature to form CnVNi clusters. Since the formation of CnVNi clusters is a process of heat releasing which will further increase the temperature of the NiAl system and produce more and more Ni vacancies, we can conclude that much more vacancies are created as a result of the presence of C impurity in NiAl intermetallics. However, the Ni vacancies exist in the form of CnVNi clusters from our calculation in a certain temperature range (less than 700 K). The existence of this kind of CnVNi cluster can effectively restrain the generations of cracks in the vacancies, which will produce some influences on the mechanical properties of NiAl intermetallic compound. Consequently, our results will provide a valuable reference for understanding the effects of C and vacancy on the mechanical properties of the NiAl intermetallics.