Outline
Comptes Rendus

Theoretical calculations applied to the study of the energetics of reactions of methamphetamine synthesis and to the characterization of reactants, products and by-products
Comptes Rendus. Chimie, Volume 14 (2011) no. 5, pp. 503-510.

Abstract

Density Functional Theory (DFT) was used to perform quantum chemical calculations at B3LYP level, with the 6-31G** basis set, as implemented by the Spartan computer program. The entities studied are the reaction of synthesis of MDMA (3,4, methylene dioxymethamphetamine), and a number of analogs of methamphetamine. The analogs have also been geometry-optimized and their infrared vibrational frequencies were calculated.

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/j.crci.2010.09.001
Keywords: 4-methylenedioxyphenyl-2-propanone, MDMA (3, 4 methylene dioxymethamphetamine), Density functional calculations, Infrared spectroscopy, Ab initio calculations, Heats of formation, Energy conversion

Danielle Sapse 1; Elise Champeil 2; Anne-Marie Sapse 2

1 John Jay College of Criminal Justice, Department of Law and Police Science, 445 W 59th St, New York, NY, 10019, USA
2 John Jay College of Criminal Justice, Department of Science, 445 W 59th St, New York, NY, 10019, USA
@article{CRCHIM_2011__14_5_503_0,
     author = {Danielle Sapse and Elise Champeil and Anne-Marie Sapse},
     title = {Theoretical calculations applied to the study of the energetics of reactions of methamphetamine synthesis and to the characterization of reactants, products and by-products},
     journal = {Comptes Rendus. Chimie},
     pages = {503--510},
     publisher = {Elsevier},
     volume = {14},
     number = {5},
     year = {2011},
     doi = {10.1016/j.crci.2010.09.001},
     language = {en},
}
TY  - JOUR
AU  - Danielle Sapse
AU  - Elise Champeil
AU  - Anne-Marie Sapse
TI  - Theoretical calculations applied to the study of the energetics of reactions of methamphetamine synthesis and to the characterization of reactants, products and by-products
JO  - Comptes Rendus. Chimie
PY  - 2011
SP  - 503
EP  - 510
VL  - 14
IS  - 5
PB  - Elsevier
DO  - 10.1016/j.crci.2010.09.001
LA  - en
ID  - CRCHIM_2011__14_5_503_0
ER  - 
%0 Journal Article
%A Danielle Sapse
%A Elise Champeil
%A Anne-Marie Sapse
%T Theoretical calculations applied to the study of the energetics of reactions of methamphetamine synthesis and to the characterization of reactants, products and by-products
%J Comptes Rendus. Chimie
%D 2011
%P 503-510
%V 14
%N 5
%I Elsevier
%R 10.1016/j.crci.2010.09.001
%G en
%F CRCHIM_2011__14_5_503_0
Danielle Sapse; Elise Champeil; Anne-Marie Sapse. Theoretical calculations applied to the study of the energetics of reactions of methamphetamine synthesis and to the characterization of reactants, products and by-products. Comptes Rendus. Chimie, Volume 14 (2011) no. 5, pp. 503-510. doi : 10.1016/j.crci.2010.09.001. https://comptes-rendus.academie-sciences.fr/chimie/articles/10.1016/j.crci.2010.09.001/

Version originale du texte intégral

1 Introduction

Methamphetamines are psychostimulant and sympathomimetic drugs. They are used recreationally, but also as medications for medical conditions and for weight loss. From 1983 to 2005, many laws were passed that made it a crime to possess equipment for methamphetamine production and chemicals that can be used to synthesize them. In order to circumvent these laws, illicit or clandestine laboratories were used, where cold medications serve as precursors to methamphetamine. The Patriot Act includes a section called the Combat Methamphetamine Epidemic Act of 2005, a law which regulates the sale of these medications and limits how much non-prescription cold medicine can be bought at a time.

There are several methods of producing methamphetamines. One of the many dangers involved in this production, is the release of harmful chemicals such as phosphine and methylamine, which are dangerous to those in the surrounding area, by causing harm to the skin, eyes, and respiratory tract. These chemicals are also very harmful to the environment, as for every pound of methamphetamine produced, there are more than five pounds of toxic waste that go into rivers or the sewer system.

The chemicals used in methamphetamine production can also lead to dangerous explosions. About 15% of methamphetamine labs are found out because of a fire or explosion.

There have been 7347 incidents of methamphetamine clandestine laboratories in the United States in 2006. The greatest numbers were found in Missouri, Illinois, and Indiana.

Methamphetamines are also produced in other countries, for instance, the Netherlands, Belgium, and Australia [1]. One of the clandestine laboratories in Australia utilized Wacker oxidation of safrole to 3,4-methylenedioxyphenyl-2-propanone (MDP-2P) using palladium chloride, p-benzoquinone and methanol.

A number of studies discuss methods of synthesizing MDMA (a designer drug known also as “ecstasy” or “XTC”) from the point of view of identifying intermediates and by-products which would facilitate the recognition of illegal MDMA production in the suspected laboratory, as well as provide information about the synthetic route employed [2].

Salouros et al. [2] describe the synthesis reaction, with the obtaining of the by-product N-cyanomethyl-N-methyl–1-(3,4-methylenedioxyphenyl)-2-propylamine. This work applies quantum chemical methods to the study of the energetics of the synthesis reactions, as well as of the properties of possible intermediates and by-products. Those properties include infrared frequencies, knowledge which might help in the identifications of different structures. Indeed, Salouros et al. [2] isolated the compound identified tentatively by Swist et al. [3] as 2-dimethylamino-2-methyl-3-(3′,4′-methylenedioxyphenyl)-2-propanenitrile. Salouros et al. [2] elucidated the structure of this compound, using NMR, as N-cyanomethyl-N-methyl-1-(3′,4′-methylenedioxyphenyl)-2-propylamine. Experiments performed by Salouros et al. [2] support Swist's conclusion that this compound is a marker for the cyanoborohydride route of synthesis.

This work investigates both possible structures of the marker via the calculations of infrared frequencies. The present work also investigates 3′,4′–methylenedioxy-methamphetamines and related compounds with quantum chemical calculations.

Other methods of identification for 3′,4′-methylenedioxy-methamphetamines and related compounds include: thin layer chromatography/fluorescence [4], gas chromatography with mass spectrometry (GC-MS) [5], liquid chromatography with fluorescence detection [6], liquid chromatography with mass spectrometry (LC-MS) [7], and immunoassay [8]. Recently, MDMA has also been identified in the urine of users via NMR spectroscopy [9].

2 Methods and results

Density Functional Theory (DFT) was used to perform quantum chemical calculations at B3LYP level, with the 6-31G** basis set, as implemented by the Spartan computer program [10]. In order to investigate the necessity of diffuse functions, the reaction (1), was also investigated with the addition of diffuse functions. The results were found to be within 2.5 Kcals/mol from the results without diffuse functions. Therefore, it was deemed that it is not necessary to use diffuse functions. The entities studied are: the reaction of synthesis of MDMA, and a number of analogs of methamphetamine. These analogs have been chosen from the experimental report of Kato et al. [4].

The energetics of the reaction:

MDP-2P+CH3NH2-H2O=N-methyl-1-MDP-2-propylimineH2=MDMA(1)
was investigated by calculating the energies of the reactants, of the intermediate and of the products. Table 1 shows these energies while Figs. 1–3 show the optimized structures, with 1 being the reactant, 3,4-methylenedioxyphenyl-2-propanone, 2 being the intermediate imine and 3 being the product, MDMA (3,4-methylenedioxymethamphetamine).

Table 1

Synthesis of MDMA: energies of the reactant, 3,4-methylenedioxyphenyl-2-propanone (1); the intermediate, N-methyl-1-(3,4-methylenedioxyphenyl)-2-propylimine (2); the product, 3,4-methylenedioxymethamphetamine MDMA (3) and by products (4) and (5).

CompoundEnergy (au)Zero-point energy (Kcal/mol)Energy of solvation (Kcal/mol)
1−612.74721113.9−8.8
2−632.17683139.1−9.0
3−633.39527154.7−6.0
4−764.93053171.5−8.6
5−764.93036171.3−7.9
CH3NH2−95.8613740.3−6.3
H2O−76.4225713.1
H2−1.178506.4
ClCH2CN−592.3392623.18−4.5
Fig. 1
Fig. 2
Fig. 3

Some by-products of the reaction are shown in Figs. 4–7.

Fig. 4
Fig. 5
Fig. 6
Fig. 7

The energies of the reactions of formation of compounds 4 and 7 were also investigated and their energies were calculated:

MDMA+ClCH2CN=N-cyanomethyl-N-methyl-1-MDP-2-propylamine+HCl(2)
N-methyl-1-phenyl-2-propylamine +ClCH2CN=N-cyanomethyl-N-methyl-1-phenyl-2-propylamine+HCl(3)

The methamphetamine analogs shown in Figs. 8–11 have also been geometry-optimized and their energies are displayed in Table 2. Tables 1 and 2 display also zero-point energies and solvation energies, the latter being calculated by the semi-empirical methods [11]. Infrared vibrational frequencies were calculated for all the structures mentioned above. In order to examine the necessity of scaling the frequencies, the experimental spectrum of dimethylmethamphetamine was compared with the spectra calculated using both HF/6-31G** and B3LYP/6-31G** methods. While the HF results showed clearly the necessity of scaling, the frequencies being larger than the experimental ones, that was not the case with the B3LYP results, where the experimental frequencies and the calculated ones were very close (within the first decimal). Table 3 shows the frequencies and the intensities of the infrared peaks of Compound 2. In addition, the calculated drawn spectrum is shown. Table 4 shows the frequencies and intensities of the by-product 4, while those of compound 5 are shown in Table 5. Both spectra of the compounds shown are drawn. Table 6 shows the frequencies and intensities of the methamphetamine analogs depicted in Figs. 8 and 9, and Table 7 shows the infrared frequencies and intensities of compounds 10 and 11. Table 8 shows the frequencies and intensities of compound 7.

Fig. 8
Fig. 9
Fig. 10
Fig. 11
Table 2

Energies of compounds 6, 7, 8, 9, 10 and 11.

CompoundEnergy (au)Zero-point EnergySolvation Energy
6−444.86572145.1−2.8
7−576.40626161.7−3.4
8−672.70090172.2−4.5
9−672.71270172.7−5.1
10−672.70880165.3−6.5
11−672.70792172.7−6.5
HCl−460.800804.2−0.7
Table 3

Frequencies and intensities of the infrared peaks of compound 2.

Frequencies (cm−1)Intensities
3106.254.3
3160.812.7
3088.221.5
3070.852.5
3008.2139.8
2998.854.2
1754.692.7
1537.3276.3
1487.140.0
1325.322.6
1294.0235.3
1265.158.6
1193.023.3
1070.0118.6
960.223.8
947.824.8
Table 4

Frequencies and intensities of the infrared peaks of by-product 4.

Frequencies (cm−1)Intensities
3114.834.4
3097.559.8
3010.8136.1
2957.491.7
2945.843.9
2371.72.0
1538.7229.0
1506.123.0
1487.555.6
1304.464.2
1296.7207.1
1223.429.4
1154.924.7
1072.1131.6
978.829.2
821.926.9
Table 5

Frequencies and intensities of the infrared peaks of by-product 5.

Frequencies (cm−1)Intensities
3152.332.8
3012.0132.9
2983.858.6
2975.136.3
2335.02.7
1560.010.6
1519.815.2
1487.457.0
1261.526.4
1224.322.7
1138.842.9
1070.6122.6
1000.435.4
956.531.6
Table 6

Frequencies and intensities of the infrared peaks of by-products 8 and 9.

Compound 8Compound 9
Frequencies (cm−1)IntensitiesFrequencies (cm−1)Intensities
3141.830.63136.230.0
312932.13123.732.1
304929.53108.643.8
3009141.63100.253.4
2936.3116.83008.4151.3
2923.855.92925.5112.4
2898.1421539.4236.5
1536253.51487.167.0
149230.71304.470.3
1372.232.51210.036.8
1292.3130.61164.240.9
1291.7155.51086.143.2
1073.182.81069.9101.8
1067.137.8971.433.5
974.636.5
Table 7

Frequencies and intensities of the infrared peaks of by-products 10 and 11.

Compound 10Compound 11
Frequencies (cm−1)IntensitiesFrequencies (cm−1)Intensities
3128.431.03134.027.8
3110.138.23117.935.2
3096.657.13105.441.9
3080.250.23102.057.5
3009.7145.13095.537.2
2942.8100.13077.938.7
2900.956.63040.836.6
1539.3244.13004.9142
1486.543.02945.397.3
1346.826.52896.462.5
1310.167.11537.8245.8
1282.3157.91486.449.5
1145.031.01402.828.0
1071.6126.11298.859.8
971.039.31295.5142.7
1284.264.5
1220.641.0
1169.239.4
1071.270.0
1069.152.0
973.6.031.2
774.268.1
Table 8

Frequencies and intensities of the infrared peaks of by-product 7.

Frequencies (cm−1)Intensities
3047.735.9
2955.289.7
2929.749.4
1403.620.4
1172.430.0
1055.228.7
736.923.0

Calculated drawn spectrum for compound 2.

Calculated drawn spectrum for compound 4.

Calculated drawn spectrum for compound 5.

Calculated drawn spectrum for compound 7.

3 Discussion

The reaction of synthesis of metamphetamine, compound shown in Fig. 3, features two steps:

  • 1. the reaction of MDP-2P (Fig. 1) with a reagent such as methylamine, leading to the formation of the intermediate shown in Fig. 2 and water;
  • 2. the reduction of the intermediate by a hydrogen molecule, leading to the formation of methamphetamine.

The energetics of these reactions, calculated as shown in Methods with the DFT method (B3LYP/6-31G**) shows the first reaction to be endothermic, with an energy of 5.8 Kcal/mol and the second reaction to be exothermic, with an energy of −25.1 Kcal/mole. As such, the overall energy of the synthesis of metamphetamine is of −19.3 Kcal/mole. If the solvation energy is taken into account, since the solvation energy of MDMA is smaller by 2 Kcal/mol than the one of the reactant (compound 1), the total energy of the reaction is of −17 Kcals/mol. The differences in zero-point energies between reactant and products are small, therefore negligible.

Salouros et al. [2] propose as a by-product of these reactions the formation of N-cyanomethyl-N-methyl-1-(3′,4’methylenedioxyphenyl)-2-propylamine, as shown in Fig. 4. According to their NMR studies, the structure shown in Fig. 4 is the right one for this by-product, but this might be reinforced by examining the IR results. In order to confirm or infirm this possibility, this work calculated the infrared vibrational modes for this compound, after subjecting it to geometry optimization. Table 4 shows these values.

Swist et al. [3] propose as a by-product 2-dimethylamino-2-methyl-3-(3,4-methylenedioxyphenyl)-propanenitrile. This compound, as shown in Fig. 5 was also geometry-optimized and its infrared vibrational frequencies were calculated. These results are shown in Table 5. Examining Tables 4 and 5, one notices that in both compounds the nitril vibration shows a very low intensity. Both compounds show strong absoptions in the 3000 cm−1 region, as well as a strong peak at 1070 cm−1. Compound 4 shows a very strong absorption at 1538 cm−1, which seems to be absent in Compound 5. Another significant difference is the strong band at 1297 cm−1, present in 4 and absent in 5.

The reaction of synthesis of the compound 4 is found to be slightly endothermic, with a reaction energy of 2 Kcal/mol. Since the solvation energy difference between product and reactant is of −2 Kcal/mol, this reaction does not occur with a change in energy.

The reaction of synthesis of compound 7 is slightly exothermic, with a reaction energy of −1.32 Kcal/mol. The solvation energy in this reaction is of only 0.6 Kcals/mol.

Table 1 shows the energies of the investigated compounds and it is clear that all are minima in energy, that is, they are stable compounds. That can be seen from the fact that all the infrared frequencies are positive.

Tables 3–8 display the infrared frequencies of all the compounds discussed. The experimentalists could take the infrared spectrum of an unknown compound suspected to be a marker, a by-product or an intermediate in the synthetic route and compare the data with the quantum chemical obtained frequencies of the suspected structure in order to obtain some information about its structure.

In conclusion, the calculations can serve to help identify the structure of the investigated compounds.


References

[1] United Nations Office on Drugs and Crime. World drug report 2006. Available at: http://www.unodc.org/unodc/en/data-and-analysis/WDR-2006.html.

[2] H. Salouros; M. Collins; G. Tarrant; A.V. George J. Forensic Sci., 53 (2008), p. 1556

[3] M. Swist; J. Wilamowski; A. Parczewski Forensic Sci. Int., 152 (2005), p. 175

[4] N. Kato; S. Fujita; H. Ohta; M. Fukuba; A. Toriba; K. Hayakawa J. Forensic Sci., 48 (2008), p. 1367

[5] W. Chen; N. Poon; M. Chan; S.D. Brown; D.J. Rhodes; B.J. Pritchard Forensic Sci. Int., 48 (2003), p. 1249

[6] F. Tagliaro; Z.D. Battisti; A. Groppi; Y. Nakahara; D. Scarcella; R. Valentini; J.L. Da Costa; A.A. Da Matta Chasin; M. Concheiro; A. De Castro; O. Quintela; M. Lopez-Rivadulla; A. Cruz Forensic Sci. Int., 723 (1999), p. 195

[7] H.K. Nordgren; O. Beck J. Anal. Toxicol., 27 (2003), p. 15

[8] D. Felscher; K. Schulz; A.G. Verstraete; F.V. Heyden J. Anal. Toxicol., 45 (2000), p. 1327

[9] J. Liu; J. Decatur; G. Proni; E. Champeil Forensic Sci. Int., 194 (2010), p. 103

[10] Spartan 06, Wavefunction, Inc., Irvine, California, 92612, USA.

[11] C.C. Chambers; G.D. Hawkins; C.J. Cramer; D.G. Truhlar J. Phys. Chem., 100 (1996), p. 16385


Comments - Policy