This article describes an approach in determination of equilibrium geometries and harmonic frequencies of molecules by the Ornstein-Uhlenbeck diffusion quantum Monte Carlo method based on the floating spherical Gaussians. In conjunction with a projected and renormalized Hellmann-Feynman gradient and an electronic energy at variational Monte Carlo and diffusion quantum Monte Carlo, respectively, the quasi-Newton algorithm implemented with the Broyden-Fletcher-Goldfarb-Shanno updated Hessian was used to find the optimized molecular geometry. We applied this approach to N-2 and H2O molecules. The geometry and harmonic frequencies calculated were consistent with some sophisticated ab initio calculated values within reasonable statistical uncertainty. (C) 2004 American Institute of Physics.