In this work, we assess the accuracy of Generalized Hartree-Fock (GHF) trial wavefunctions in ph-AFQMC calculations of strongly correlated molecules. GHF trial wavefunctions have previously been used successfully in constrained AFQMC calculations of lattice models, including the two-dimensional Hubbard model and the Hubbard-Kanamori model. Selfconsistent AFQMC calculations using trials of GHF form were performed on the hole-doped three-band Hubbard model. Symmetry-broken multi-determinant GHF expansions have been explored as trial wavefunctions for AFQMC calculations on Hubbard models of varying geometries, and symmetry-projected GHF wavefunctions were shown to be very promising for the two-dimensional Hubbard model. The GHF basis is also more amenable to augmentation of the Hamiltonian with relativistic terms, such as the spin-orbit coupling, which has been explored in the context of both molecules 32 and Hubbard models. There are many attractive features and properties of GHF wavefunctions, which lead us to hypothesize that they may be a suitable choice as trial wavefunctions for ph-AFQMC.