# How to run crystal simulation on GAMESS?

I know GAMESS is not best suited for crystal simulation but I am not ready to start learning another package as VASP. As I heard GAMESS is able to do some calculations in a "terminated" crystal (I mean using an ion in the periphery instead of repeating it multiple times). Mostly I need to calculate HOMO-LUMO gap / hessian matrix and transitional state of a surface reaction .

I looked in the web but couldn't find a good worked example of using GAMESS for crystals. Is there anywhere I can find such an example ?

Since you mention wanting to treat a surface reaction, you'd use a builder (e.g., Avogadro) to create a large cluster (e.g., Build->Supercell). So rather than one primitive unit cell, you create maybe a 2x2x2 or 3x3x3 slab. If needed, you'll terminate dangling bonds (e.g., adding hydrogens to cap oxygens in $\ce{TiO2}$ for example.. notice the terminal "uncapped" atoms on the back and right of the unit cell).
Now, this isn't as strange as it sounds. In some cases ($\ce{MgO}$ and $\ce{TiO2}$ are good examples) the properties require more advanced calculation methods like multi-configuration QM, and this may be easier to do in a quantum chemical package with a cluster, rather than a solid-state crystal code.