I am unable to achieve convergence in my VC Relax calculation for a doped Bismuth structure using Quantum Espresso. When I ran the calculation for only Bismuth it converged but after I doped it with Cu, It has failed to converge. With some of the input files, I get the error: "Error in routine bfgs (1): dE0s is positive which should never happen", while with other similar doped structures the calculation terminates upon reaching the maximum number of steps. What could be the reason for this? I have attached here the input file and output file of Bi that achieved relaxation/convergence as well as the doped Bi (Bi and Cu) that did not achieve convergence. Please let me know your thoughts.