I guess what you're looking for is to run a geometry optimization. This is quite straightforward in ATK, using the calculateOptimizedGeometry() function, although ATK might not be the most accurate code for molecules. You can also set up the relaxation from the Scripter in VNL.