I want to derive the Hamiltonian of the ion molecule H_2^+ by making use of Jordan-Wigner map. 
Thus, I used qiskit nature. I wrote down same command of given in qiskit nature. The only difference is that I change the charge number from 0 to 1. Then it generates following error message 
"Traceback (most recent call last):
  Input In [6] in <cell line: 2>
    second_q_op = es_problem.second_q_ops()
  File /opt/conda/lib/python3.8/site-packages/qiskit_nature/problems/second_quantization/electronic/electronic_structure_problem.py:94 in second_q_ops
    driver_result = self.driver.run()
  File /opt/conda/lib/python3.8/site-packages/qiskit_nature/drivers/second_quantization/electronic_structure_molecule_driver.py:179 in run
    return driver.run()
  File /opt/conda/lib/python3.8/site-packages/qiskit_nature/drivers/second_quantization/pyscfd/pyscfdriver.py:410 in run
    self._build_molecule()
  File /opt/conda/lib/python3.8/site-packages/qiskit_nature/drivers/second_quantization/pyscfd/pyscfdriver.py:460 in _build_molecule
    raise QiskitNatureError("Failed to build the PySCF Molecule object.") from exc
QiskitNatureError: 'Failed to build the PySCF Molecule object.'

Use %tb to get the full traceback."

If you know any method to fix it, please teach me. Thank you for your time.