I know that we can easily transform Molecular Hamiltonians to qubit operators using qiskit nature (FermionicOp Mappers). But what about an arbitrary Hamiltonian? Suppose that I have a first-quantized Hamiltonian as follows:
$$ H(r_1,r_2)=-\frac{1}{2m_1}\nabla^{2}_1-\frac{1}{2m_2}\nabla^{2}_2-\frac{1}{\lvert r_1 -r_2 \rvert}+\frac{1}{2}\omega^2 (m_1 r_1^2+m_2 r_2^2) $$ where $m_1$ and $m_2$ are the masses of the particles and $\omega$ is the frequency of a harmonic oscillator. What is the most straightforward way (including qiskit or qiskit nature) to transform this Hamiltonian to qubit operators and ready for the next quantum computations like VQE algorithm?
