The catalytic properties of zeolites, which are primarily determined by the framework topology and active centers (e.g. Ti 4+ , Sn 4+ or Al 3+ ), remain challenging to be controlled in zeolite synthesis. Here, we combined first-principal and classical molecular simulations to investigate how Al siting, and OSDA orientation impacts the energy of a zeolite supercell. 36 T-site CHA zeolite with TMAda + (OSDA) was chosen as the model system. By applying a Boltzmann factor to each Al configuration and as a function of TMAda + orientation, we came to the conclusion Al pairs prefer to locate in 8-MRs compared to 6-MR, 4-MR and D6R, which is consistent with our previous experimental finding. We also found that the potential energy was governed by the distance between the anionic AlO − 4 tetrahedra and the cationic quaternary ammonium