Author Topic: An old question-Pulay mixing inversion failed. Using only last step  (Read 2355 times)

0 Members and 1 Guest are viewing this topic.

Offline songll0808

  • Regular QuantumATK user
  • **
  • Posts: 6
  • Reputation: 0
    • View Profile
Hi everyone, I have set up a two probe system. When I make scf calculation, the calculation would finish and give the message:Pulay mixing inversion failed. Using only last step.

I have check my system, and I can make sure that the structure is all right. I also improve the accuracy of the energy contour integral parameters, but the results is the same.

Thank for any replies.

Offline Nordland

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 812
  • Reputation: 18
    • View Profile
Have you removed anything in the output?

What are you mixing factor?

Something is strange. It appears to converged in the first step for the two-probe geometry - so perfectly that the pulay mixer warns. I have never seen this before,
 and therefore I suspect that something is wrong.

Offline fanjiaping

  • Heavy QuantumATK user
  • ***
  • Posts: 67
  • Reputation: 0
    • View Profile
 hi nordland:
       how can you know the system converged from the result?

Offline Anders Blom

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 5416
  • Country: dk
  • Reputation: 89
    • View Profile
    • QuantumATK at Synopsys
It converged numerically in 4 steps, that's what the log file shows. In step 4 you see
Code
# sc  4 : q =  412.90491 e  Etot = -4308.82930 Ry  dRho =  0.0000E+00  dEtot =  0.0000E+00 Ry
which is below the tolerance for dEtot. However, having a convergence to exactly dE=0 and dRho=0 is very suspicious and shows something is wrong. Indeed, as has been discussed separately offline, the geometry was not correctly set up. With the correct geometry, all works fine.