dft/pbe0/6-31g* hess(S0,type=analytical) dftgrid(rad_npts=150,ang_npts=974)
geom="""
N      0.000000     0.000000     0.000000
H      0.000000     0.999351     0.000000
H     -0.865464    -0.499676     0.000000
H      0.865464    -0.499676     0.000000
"""