Fix: correct the way to use two-center-intergral and support more hydrogen-like orbitals strategy

[kirk0830]

This cannot solve the bad condition of QO overlap.

Reminder

• Have you linked an issue with this pull request?
• Have you noticed possible changes of behavior below or in the linked issue?
• Have you explained the changes of codes in core modules of ESolver, HSolver, ElecState, Hamilt, Operator or Psi? (ignore if not applicable)

Linked Issue

Fix #3607

What's changed?

• correct the last param in initializing two-center-integral.
• add one more function to generate hydrogen-like orbital, for the most realistic, and the smallest space

Any changes of core modules? (ignore if not applicable)

No.

[dyzheng]

LGTM