<div dir="ltr">Hi Jen,<br><br>sure, simply read the ground state wave function and specify the appropriate excitation in the MOM section. Alternatively, you can start the calculation using the atomic guess, if you wish, but it is usually less stable. <br><br>It is also possible to restart the MOM calculation using an MOM wave function for the same excited state. In this case simply read the appropriate wave function and provide the empty MOM section (&MOM / &END MOM). The MOM procedure will then select those orbitals which overlap most with the occupied orbitals from your restart file.<br><br>Sergey<br><br>On Thursday, September 15, 2016 at 6:35:25 PM UTC+1, Jen E. wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Thanks very much! I am in the process of trying out the recommended steps. I've managed to converge the ground state and determine which orbitals I want to occupy. Is the proper next step to read back in the converged ground state wavefunction and then add the MOM tag to the input and re-run the calc? <br><br>Best, <br>Jen<br><br>On Wednesday, September 14, 2016 at 2:10:04 AM UTC-4, Matt W wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>the geo_opt bug was fixed shortly after the posts above - so a recent executable should work fine (the bug was actually for MOM in combination with the MO_CUBES keyword).</div><div><br></div><div>We don't really have a definite protocol for running the method beyond: Converge ground state; Check orbitals; Run MOM calc with appropriate excitations. You might need to set ALPHA in &MIXING quite small (< 0.1). </div><div><br></div><div>For a well separated S1 state, I'd hope it might work quite OK.</div><div><br></div><div>Matt<br><br>On Tuesday, September 13, 2016 at 9:59:13 PM UTC+1, Jen E. wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi All, <br><br>I am a new user of CP2k and I am trying to use the MOM method to converge the S1 excited state of a molecule. I have a few related questions.. 1.) I was wondering if there is a specific protocol for running a MOM calculation  and 2.) Has there been any resolution in getting the MOM method to work with a geometry optimization? <br><br>Thanks very much! <br><br>Jen <br><br>On Tuesday, April 19, 2016 at 4:58:16 AM UTC-4, Matt W wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">HI Tassem,<div><br></div><div>it should work, in principle. One think that could cause trouble is if the step size is large during the geometry optimization so the orbitals change too much.</div><div><br></div><div>If you send me some output offline, we can have a look.</div><div><br></div><div>Matt<br><br>On Monday, April 18, 2016 at 4:10:30 PM UTC+1, Al-Moatasem El-Sayed wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi all,<div><br></div><div>I am trying to run a geometry optimization while using the maximum overlap method to occupy an excited state. When I try to do this, the first SCF cycle successfully occupies the higher state I would like to be occupied. I also manage to occupy the state I would like if I just run a single point calculation. However, when the geometry optimization begins, the electron moves back into the lowest unoccupied orbital. So I was wondering if it is actually possible to use the maximum overlap method with a geometry optimization.</div><div><br></div><div>Thanks,</div><div>Tassem</div></div></blockquote></div></div></blockquote></div></blockquote></div></div></blockquote></div></blockquote></div>