I am studying adsorption of organic molecules on a surface using a slab model. To find low energy configurations of molecule on surface, I am using quench MD.
fix 10 all nvt temp 300.0 300.0 0.1
run 1000000
unfix 10
fix 10 all nvt temp 300.0 0.1 0.1
run 1000000
unfix 10
minimize 1.0e-9 1.0e-9 10000 100000
And repeated cycles of this procedure.
Is there any other better method to get the low energy configuration?
Please understand that this is not really a topic for the LAMMPS category, since this has nothing to do with specific LAMMPS features but is a question about MD methodology overall. I thus suggest you look through text books and review publications that discuss searching for low (free) energy configurations or rare events (those are related). There are several possible methods that can be used to more efficiently sample the available phase space. Once you have settled on some suitable approaches you can come back here and discuss how to realize those in LAMMPS.