How to change the simulator used by Qiskit's `QAOA` (or get around the "Maximum allowed dimension exceeded" error)?

Question

I'm running MinimumEigenOptimizer with QAOA on a large quadratic program:

qp = QuadraticProgram(....) 

meo = MinimumEigenOptimizer(QAOA(....)) 

However, I'm getting the feared ValueError: Maximum allowed dimension exceeded from qiskit.

It seems to me like there are 2 barriers to overcome for this:

1. Qiskit's maximum number of qubits for the given simulator.

2. numpy's maximum array size.

I tried to address 1) by choosing a different simulator. In the previous version of QAOA there used to be an option for choosing the quantum_instance (https://qiskit.org/documentation/stubs/qiskit.algorithms.QAOA.html). But since the recent migration, I can't find this option: (https://qiskit.org/documentation/stubs/qiskit.algorithms.minimum_eigensolvers.QAOA.html). How do I change the simulator (or another option if it exists) to allow a larger number of qubits?