We investigate the rich stereochemistry of cocaine and its diastereoisomers from a theoretical perspective using density functional theory. The relative stability of the eight considered isomers is discussed, and a comparison of the corresponding internal coordinates is given. Our results reveal that the S-pseudococaine isomer is the most stable conformation, whereas the natural occurring isomer (R-cocaine) lies higher in energy. The different isomers' chemical reactivity is discussed based on the calculation of the hardness, softness, electrophilicity and dipole moment. It was found that the dipole moment varies over a broad range from 0.65 to 4.60 D, whereas the other properties are slightly modified. The solvent effect on the energy stability of the cocaine isomers was studied by considering chloroform, dimethyl-sulfoxide (DMSO) and water as implicit solvents. Our calculations show that the different isomers' energy order and their energy gaps are slightly modified due to solvent effects. However, in all cases, the S-pseudococaine remains the most stable isomer. However, the dipole moment and the chemical reactivity of the cocaine isomers increase with the solvent polarity. Keywords: Cocaine isomers,DFT, stability, solvent effect, chemical reactivity.