<div>Hi all,<br></div><div><br></div><div>I have a system that I've run and have good wavefunctions for. I need to try and calculate STM images for this system. I want to add additional ghost basis functions into the vacuum space to improve the description of the wavefunctions here. </div><div><br></div><div>So the question: is  there a hack to let me do this but keep the converged wavefunctions? </div><div><br></div><div>At the moment it reverts to atomic guess if the number of atoms changes. What I'd want is to read in the mo_coeffs from the file and zero the coeffs of the ghost atoms and reortho...is this doable?</div><div><br></div><div>Cheers,</div><div><br></div><div>Matt</div>