  We extend the recently introduced phaseless auxiliary-field quantum Monte
Carlo (QMC) approach to any single-particle basis, and apply it to molecular
systems with Gaussian basis sets. QMC methods in general scale favorably with
system size, as a low power. A QMC approach with auxiliary fields in principle
allows an exact solution of the Schrodinger equation in the chosen basis.
However, the well-known sign/phase problem causes the statistical noise to
increase exponentially. The phaseless method controls this problem by
constraining the paths in the auxiliary-field path integrals with an
approximate phase condition that depends on a trial wave function. In the
present calculations, the trial wave function is a single Slater determinant
from a Hartree-Fock calculation. The calculated all-electron total energies
show typical systematic errors of no more than a few milli-Hartrees compared to
exact results. At equilibrium geometries in the molecules we studied, this
accuracy is roughly comparable to that of coupled-cluster with single and
double excitations and with non-iterative triples, CCSD(T). For stretched bonds
in H$_2$O, our method exhibits better overall accuracy and a more uniform
behavior than CCSD(T).
