Recently, there has been a worldwide substantial increase in the consumption of new psychoactive substances (NPS), compounds that mimic the structure of illicit drugs, such as amphetamines or ecstasy. The producers try to avoid the law by a slight modification of illicit structures, thereby developing dozens of temporarily legal NPS every year. The current trends in the detection and monitoring of such substances demand a fast and reliable analysis. Molecular spectroscopy represents a highly effective tool for the identification of NPS and chiroptical methods can provide further information on their 3D structure, which is the key for the determination of their biological activity. We present the first systematic study of NPS, specifically butylone, combining chiroptical and vibrational spectroscopies with ab initio calculations. According to density functional theory calculations, 6 stable lowest energy conformers of butylone were found and their molecular structure was described. For each conformer, the relative abundance based on the Boltzmann distribution was estimated, their population weighted spectra predicted and compared to the experimental results. Very good agreement between the experimental and the simulated spectra was achieved, which allowed not only the assignment of the absolute configuration, but also a precise description of the molecular structure.