I was trying to optimize a molecule using DFT/b3lyp in Gaussian. But after the optimization, it was found that some bonds are missing from the molecule in the optimized structure. And also in certain case, the optimized mlolecule was entirely different from the input structure. Can anyone help me why it happens and how can I solve this problem?