Molecular simulation techniques are used to find the basal spacing of organoclay on the basis of the energy minimum, using the canonical NVT ensemble. Then, the interlayer structure and mobility of alkyl chains are explored based on the interlayer atomic density profiles. Besides the basic lateral-monolayer arrangement, lateral-bilayer accompanied by partial a pseudo-trilayer and a transition structure between the two basic lateral models are observed. The later provides an excellent explanation about the reflection at 16 angstrom on XRD patterns in the literature. The atomic density profiles reveal that nitrogen atoms show stronger layering behavior than carbon atoms do. Our simulation demonstrates that the molecular mobility of the confined alkyl chains decreases from lateral-monolayer to lateral-bilayer with the increase of the intercalated surfactant. This is in accordance with the suggestion deduced from experiments. Furthermore, our simulation indicates that the mobility of the alkyl chains strongly depends on the surfactant arrangement rather than the surfactant packing density.