Having a proper input geometry for the calculations is absolutely crucial in order to obtain quality results. There is however no universal method for that, but rather one has to combine common physical and chemical sense with various methods. In some cases something similar like the Li4H2 cluster approach will be useful, otherwise relaxing an equivalent bulk system (smaller or larger than the actual central region) works better.
