I am not very familiar with scanning a PES for a bulk structure, and I'd be very appreciative for some suggestions.
The potential energy surface (PES) is a 3N-dimensional function for a bulk system containing N atoms (in reality 3N-3 to account for the trivial translational degrees of freedom). For a bulk structure, N typically represents the number of atoms in a simualtion cell with periodic boundary conditions, which is of the order of $10^2$-$10^3$, so the function is of very high dimension. The number $10^2$-$10^3$ comes from a typical DFT calculation; if you use cheaper methods (e.g. force fields) then that number can be larger, or if the method is more expensive (e.g. quantum chemistry) then the number is somewhat smaller. In any case, you are sampling a high dimensional function, which is an extremely hard problem.
So what can be done to scan this potential energy surface? These are a few options:
- Phonon calculations. A (meta)stable bulk structure sits at a local minimum of the PES. A phonon calculation determines the Hessian around this minimum, that is, the curvature of the PES around the minimum. Phonon calculations are relatively straight-forward to perform (say at the DFT level), and typically provide very useful information about the low-energy part of the PES. This is useful for calculating thermodynamic properties of a crystal at relatively low temperatures (when atoms do not move very far away from the minimum). One can explore a little further away by including terms beyond the second derivative in the expansion about the minimum (referred to as anharmonic terms), but these still essentially give a local view of the PES around the minimum.
- Transition states. Considering two structures associated with two nearby minima in the PES, it is possible to locate the minimum energy path between the two. This path goes through a saddle point of the PES, which is called a transition state. A well-known technique to find this is the nudged elastic band method, which requires knowledge of the two minima and calculates the saddle point. Another method that can allow one to "jump" from one minimum to another nearby minimum is molecular dynamics. These methods still provide a local view of the PES, but they explore it in a wider region than a phonon calculation would.
- Structure prediction. If you want a more uniform sampling of the PES, then you are entering the realm of structure prediction. These methods try to explore as many minima as possible of the potential energy surface. The basic idea is to generate structures at random (corresponding to arbitrary points in the PES) and then relax them to the local minimum. Repeating this many times allows you to explore basins around different minima of the PES. There are different methods that allow you to do this, which differ by how they generate the "next" structure (e.g. stochastic, genetic, particle swarm algorithms). These methods provide a global view of the potential energy surface, but lack in detail, so, for example, would have to be complemented with phonon calculations to explore the important low-energy region of a material at low temperatures.
David Wales has a book on some of this, you can find more details in his website.