Having this https://github.com/aromanro/HartreeFock/issues/16 implemented, it's just an optimization problem. The code could compute the molecule structure.
