The alpha-, beta-, and gamma-cyclodextrin (CyDs) dimers were studied by molecular dynamics (MD) simulations in water as an explicit solvent. The relative stability of dimers and the involved molecular interactions were determined. Three possible starting orientations were considered for the dimers: head-to-head, head-to-tail, and tail-to-tail. MD simulations were performed over a period of 5 ns to ensure the stability of the system for both the CyD dimers and monomers. The MM-PBSA methodology was used to obtain the free binding energy of the dimers and to determine the most stable arrangement for each solvated CyD. In a vacuum, MD simulations provided the head-to-head orientation as the most stable orientation for the three CyDs, while in aqueous solution the, the head-to-tail orientation was found to be the most stable for the alpha-CyD dimer and the tail-to-tail orientation the most stable for the beta- and gamma-CyD dimers.