include procedures geometry={O; !Z-matrix for water H1,O,R; H2,O,R,H1,THETA} R=0.96 Ang !start bond distance Theta=104 !start bond angle basis=6-31g** !Pople basis set $functional=b3lyp !define fucntional (optional, b3lyp is default) freqdft !run frequency calculation