A nonlinear energy method has been developed recently to construct global in time smooth solutions near Maxwellians for the Boltzmann equation. As the key feature of such a method, the linearized collision operator L can be shown to be positive definite for solutions to the Boltzmann equation close to Maxwellians, even in the presence of spatial dependence. The main idea and applications of such a method will be discussed.