Fast preparation of quantum many-body states is essential for myriad quantum algorithms and metrological applications. Here, we develop a new pathway for fast, nonadiabatic preparation of quantum many-body states that combines quantum optimal control with a variational Ansatz based on non-Gaussian states. We demonstrate our approach on the spin-boson model, a single spin interacting with the bath. We use a multipolaron Ansatz to prepare near-critical ground states. For one mode, we achieve a reduction in infidelity of up to ≈60 (≈10) times compared to linear (optimized local adiabatic) ramps; for many modes, we achieve a reduction in infidelity of up to ≈5 times compared to nonadiabatic linear ramps. Further, we show that the typical control quantity, the leakage from the variational manifold, provides only a loose bound on the state's fidelity. Instead, in analogy to the bond dimension of matrix product states, we suggest a controlled convergence criterion based on the number of polarons. Finally, motivated by the possibility of realizations in trapped ions, we study the dynamics of a system with bath properties going beyond the paradigm of (sub- and/or super-) Ohmic couplings.