A numerical method of solving time-dependent Hartree-Fock-Bogoliubov equation with Gogny interaction

A numerical method to solve time-dependent Hartree-Fock-Bogoliubov equation is proposed to treat the case with Gogny effective interaction. To check the feasibility of the method, it is applied to oxygen isotope O20 and small amplitude oscillations are calculated. The conservation of the nucleon numbers as well as energy expectation value is demonstrated. The strength distributions of the small amplitude quadrupole oscillations are also shown.