The condensation of silicic acid with aluminate in alkaline environment, the essential reaction of zeolite synthesis, is studied using the density functional theory, with the hybrid functional B3LYP in conjunction with the 6-311++G(d, p) basis set. The Si(OH) 4 monomer and Al(OH) 4 -anion are used as the reactant models to study the condensation pathway in basic solution. The solvent effect is included by the COSMO-RS model. The study includes the complete geometry optimization and frequency calculation of reactants, products, reaction intermediates, and transition states, as well as the calculation of the activation energy of the different pathways involved. The intrinsic reaction coordinate method is used to verify the reactant and product corresponding to the transition state. The calculation shows that the formation of Si-O-Al linkage can proceed via two possible reaction pathways. The first is a single-step process, in which the formation of SiOÁÁÁAl bond and removal of water are synchronous, with the activation energy of 83.7 kJ/mol. The second is a stepwise route, in which the AlOÁÁÁSi bond is first formed to give a 5-coordianted Si intermediate, and then water is removed to yield a dimer aluminosilicate, with the barriers of 62.7 and 69.3 kJ/mol for the two steps, respectively.