We simulated three different thunderstorm events over India using a regional model with two different multi-moments cloud microphysical schemes: (i) Morrison double moment (Morrison-2 M) and (ii) National Taiwan University (NTU) triple moment (NTU-3 M). The results show that the choice of microphysical scheme significantly impacts storm structures, cloud formations, lightning intensity, and rainfall patterns. The lightning flash counts from the in-situ lightning detection network (LDN) are also used to compare the simulated storms. The Lightning Potential Index (LPI) was computed for both Morrison-2 M and NTU-3 M microphysics schemes and compared it with the LDN observation. Differences in simulations are attributed to the simulation of ice crystals, snow, and graupel by the microphysical schemes. The effects on the size distributions of cloud hydrometeors between the two microphysical schemes are also found to be important. The inclusions of ice crystal shapes likely accounted for the key differences between the two microphysics schemes. The results demonstrate that the simulation of lightning event characteristics is sensitive to the choice of microphysical parameterization schemes in numerical weather prediction models. In comparison to both the microphysical schemes, some aspects are better in the 2 M microphysics scheme, and some are better in the 3 M microphysics scheme. Morrison-2 M simulates more mixing ratio for the majority of cloud hydrometeors across all cases than the NTU-3 M whereas the Morrison-2 M has less rainfall bias over the ocean for the majority of the cases. Overall, the present study suggests that the choice of microphysical scheme can greatly influence the accuracy of thunderstorm and lightning event predictions in operational weather forecasting models.