Summary Grating spacing dependence of nonvolatile holographic recording in doubly doped lithium niobate crystals is theoretically investigated allowing arbitrary charge transport lengths. It is shown that the nonvolatile refractive index modulation initially increases with increasing grating spacing, then a saturation behavior arises because of the dominant bulk photovoltaic effect. Although different charge transport length results in different nonvolatile refractive index modulation, the grating spacing dependence of nonvolatile holographic recording obeys almost the same rules for arbitrary charge transport lengths. The experimental results obtained by recording nonvolatile holograms in LiNbO 3 :Cu:Ce crystals with different grating spacing are consistent with the theoretical analyses.