A theoretical study on the reaction of aluminum with water in the gas phase was performed using the hybrid density functional B3LYP and QCISD(T) methods with the 6-311+G(d,p) and the 6-311++G(d,p) basis sets. The results show that there are three possible reaction pathways that involve four isomers, seven transition structures, and two possible products for the reaction of aluminum with water. The two most favorable reaction pathways were found, whose intermediates and products agreed quite well with experimental results. The enthalpy and Gibbs free energy change of the reaction between A1 and H2O at 298 and 2000 K were calculated. Some results are also in good agreement with the previous calculations or experimental results.