In strongly coupled field theories, perturbation theory cannot be employed to study the low-energy spectrum. Thus, non-perturbative techniques are required. We employ the variational method, a rigorous, non-perturbative approach which provides variational upper bounds on the energy eigenstates of the theory. An essential step in the variational method is the choice of trial wave function. In this work, we study the viability of employing a neural network as our variational ansatz. As a first step towards phenomenologically interesting strongly coupled theories like quantum chromodynamics, we study scalar field theories with quartic couplings.