There are a few reasons to like gaussian, and many reasons not to. Gaussian is fast, and their whitepapers are great resources for learning...

Another update: My jobs have run long enough now that I can confirm that using your own z matrix still does not constrain the bond lengths...

Update 2: The coordinates are actually gradients, and so aren't terribly informative to a casual user like myself. See this post for how...