symbols = ["H", "H"] coordinates = np.array([0.0, 0.0, -0.6614, 0.0, 0.0, 0.6614]) basis_set = "sto-3g" electrons = 2 H, qubits = qchem.molecular_hamiltonian( symbols, coordinates, basis=basis_set, )