# Modifying Qiskit Hamiltonian

I'm fairly new to the Qiskit API. I was wondering if I could get some assistance with trying to implement our technique of projecting out a new Hamiltonian from the ground-state Hamiltonian. The technique is discussed in our preprint.

There are two main challenges that I'm confronted with:

1. The qiskit constructed Hamiltonian by default consist of only the non-zero Pauli strings. How can this be changed to allow for adding and removing terms that may be initial zero.
2. How and where would one dynamically modify each Pauli coefficient in the Hamiltonian? For example it appears something like qubitop.aer_pauli[i,:] = 1.0 + 0.0j should work, but I'm not sure this is complete or an appropriate approach.

This would enable us to then use the VQE routines with Qiskit to further verify this method.

I'm looking to make these modification to the Qiskit source code just need some guidance.