J. Chem. Sci. Vol. 129, No. 8, August 2017, pp. 1301–1317. © Indian Academy of Sciences DOI 10.1007/s12039-017-1322-2 REGULAR ARTICLE
Solvolysis of organophosphorus pesticide parathion with simple and α nucleophiles: a theoretical study
CHANDAN SAHU and ABHIJIT K DAS∗ Department of Spectroscopy, Indian Association for the Cultivation of Science, Jadavpur, Kolkata, West Bengal 700 032, India E-mail: [email protected]
MS received 10 April 2017; revised 2 June 2017; accepted 4 June 2017
Abstract. Density functional theory (DFT) has been used to study the solvolysis process of the organophosphorus compound, O,O-diethyl p-nitrophenyl thiophosphate (Parathion, PTH) with α-nucleophiles − − − [hydroxylamine anion (NH2O ), hydroperoxide (HOO ) and simple nucleophile methylthiolate (CH3S )in both gas and aqueous phases. Formation of a trigonal bipyramidal intermediate at the phosphorus center followed by elimination of leaving group is considered to be the major solvolyzed pathway through addition-elimination scheme. In this study, although there are two possible orientations for incoming nucleophiles with respect to the substrate, the present reaction mechanism is found to be independent of this relative orientation. The proposed addition-elimination mechanism has been first explored here. The results indicate that the α-effect is observed − − in presence of solvent. Free energy barriers for NH2O and HOO are comparable and lower than that for the − simple nulcleophile, CH3S . An important physical insight of this study is that there is a significant influence of the reaction medium on the nucleophilic reaction for solvolysis of PTH irrespective of the relative orientation of incoming nucleophile group.
Keywords. Parathion; α-nucleophiles; Solvolysis; DFT; potential energy surface; pesticide; softness; Fukui function.
1. Introduction environment friendly due to the formation of prod- ucts that have mild or acute toxicity. 8,12 Some recent Parathion (PTH), an organophosphate (OP), is widely investigations pointed out that α-nucleophiles play an used as a pesticide to protect crops from insects. 1 OPs effective role in the degradation of the above type of show significant toxicity towards mammalian organ- toxic compounds. 13,14 Several theoretical studies were isms; 2,3 they inhibit the activity of the enzyme acetyl- also carried out for the thermal unimolecular decompo- cholinesterase, which is involved in the transmission of sition, 15 hydrolysis, 16,17 α-nucleophilic destruction 18,19 nerve impulses through phosphorylation. On the other of such lethal agents and reported that α-nucleophiles hand, increasing use of these insecticides contaminates are highly effective for decomposition of toxic OP pes- the environment, particularly soil and ground water. ticides under mild condition. 20–24 As α-Nucleophiles Therefore, due to the high toxicity and bio-accumulation bear non-bonding pairs of electrons at the position R shown by these organophosphate compounds, many to the nucleophilic center, their reactivity is greater than methods have been developed for their degradation. 4,5 that would be expected on the basis of the pKa val- The methods available presently for the degradation OPs ues. 25 Moreover, the high reactivity of α-Nucleophiles are homogeneous and heterogeneous hydrolysis, 6 pho- towards phosphorus makes them reagents of special tolysis, 7 photochemical degradation, 8 biodegradation, 9 interest for the destruction of nerve gases and other the chemical treatments based on the use of nucle- organophosphorus poisons. 26 But the origin of α-effect ophiles 4,5,9 or α-nucleophiles. 10,11 However, some of remains controversial till today. 20,25,27–44 Earlier theoret- these treatments may not be very efficient or are not ical studies suggested that the origin of α-effect includes ground–state (GS) destabilization, transition-state (TS) *For correspondence stabilization, thermodynamic stabilization of products Electronic supplementary material: The online version of this article (doi:10.1007/s12039-017-1322-2) contains supplementary material, which is available to authorized users.
1301 1302 Chandan Sahu and Abhijit K Das and solvent effects. 20,25,27–44 Interestingly, solvent effect the phosphorus center of phosphate triesters having an − is also contentious for the α-effect, whether it is intrin- aryl-derived leaving group by CH3S nucleophile in gas sic properties of the nucleophiles or solvent induced phase. phenomenon. 25,28–31,33,35,42–50 The gas phase reactions In order to investigate the α-effect of the nucleophiles, between various anions and OP compounds have been we have performed a systematic computational study investigated using mass spectrometry 51,52 in absence of of the solvolysis reaction of PTH with α-nucleophiles − − − solvents and these findings reported that α-nucleophiles (NH2O , HOO ) and simple nucleophile, CH3S in show greater reactivity than normal nucleophiles. Biber- both gas and aqueous phases. As the biological activ- baum et al. 28 in their study of gas phase reaction ity of molecules depends on the electronic structure of between α-nucleophiles and alkyl chlorides using a tan- the active part of a compound and its conformation, dem flowing afterglow-selected ion flow tube instrument we first performed a conformational analysis of PTH. concluded that the α-effect is not due to an intrinsic As the nucleophilic attack at the phosphorous center property of the anion, instead, it is due to the solvent (SN2@P) is the important pathway for solvolysis of OP effect. But this is in contradiction to other experi- pesticides, 19,29,54,55 a detailed study of these pathways mental and theoretical studies of McAnoy et al., 29,53 has been discussed for these three nucleophiles for the with gas phase α-effect. In an experimental study, solvolysis of PTH. In addition, a comprehensive analy- DePuy et al. 28 reported that nucleophilicity of perox- sis has been performed to know the reactivity of the three − − − ide anion is similar to hydroxide in gas phase and nucleophiles, NH2O , HOO and CH3S in the SN2@P the α-effect is not significant in the absence of sol- reaction using conceptual density functional theory. vent. These results motivated us to investigate solvent effects on α-nucleophiles as well as on simple nucle- ophiles (no α-effect). However, depending primarily 2. Computational details on the alkylation or arylation state of OPs and the 19 To find out stable conformers of PTH, a molecular dynam- nature of nucleophile, different pathways of degrada- ics (MD) conformational search has been performed with tion were reported. It has been reported that the attack an unconstrained MD trajectory using the Verlet velocity at the phosphorus center, SN2(P) is the sole reaction algorithm and NVE thermostat along with other default 54 pathway for the reactions of Paraxon and its sulfur parameters in Gabedit V.2.3.8. 61 The PM6 semi-empirical analogue, Parathion. Another study of alkaline perhy- method, as implemented in MOPAC 2009, 62 has been used drolysis of a model VX compound reported that all the to find the minimum conformational geometries. Twenty five hydrolysis reactions proceed through the phosphorus- representative minimum structures, which are selected from centered pentavalent intermediate. 29,55 There are two the rigorous conformational search for final validation of different types of α-nucleophiles and one simple nucle- the conformational analysis, are studied by the density func- tional theory (DFT). The DFT, with the hybrid functional ophile used for nucleophilic destruction. Among various 63 α-nucleophiles, the (-N-O- type) α-nucleophiles, i.e., (M06-2X) of Truhlar and Zhao, has been employed to fully optimize the geometries of all the molecular species hydroxybenzotriozoles have been used recently for the involved in this study. The standard 6-31++G(d,p) 64 basis solvolysis of organophosphate esters, which show faster 56 set is used for all the atoms. The theoretical level for this rate of solvolysis for such toxic esters. This finding combination of method and basis set is denoted as M062X/6- motivated us to explore the solvolysis of PTH with 31++G(d,p). The normal-mode analyses have been performed − 57 NH2O nucleophiles. Previous theoretical study of at the same level of theory for reactants and products as well solvolysis of VX using hydroxide and α-nucleophile as TS geometries, and the minima are characterized with no hydroperoxide inferred that hydroperoxide is a better imaginary frequency, whereas the presence of one imaginary nucleophile compared to hydroxide as hydroperox- frequency is the characteristic of TS. To evaluate the zero- ide solely produces non-toxic products. Therefore, we point vibration energy (ZPVE) and thermal corrections to have taken HOO− as a α-nuclophile for detoxifica- the Gibbs free energy at T = 298.15 K, harmonic vibra- − tional frequencies are calculated at the M062X/6-31++G(d,p) tion. Other type of simple nucleophiles, namely CH3S and CH O− were reported 58 theoretically for solvoly- level. The first-order saddle points, which are the transition 3 states that connect the equilibrium geometries, are obtained sis of different type of phosphate triesters. Xia et al. 59 using the synchronous transit-guided quasi-Newton (STQN) reported a detail theoretical study of the methanoly- method. A parallel intrinsic reaction coordinate (IRC) calcu- sis of paraoxon and analogous reactions with sulfur lation with all transition states has been performed to confirm substituted at the key oxygen position. However, no whether these transition states connect the right minima or theoretical work is available so far in the literature not. 65,66 Single-point energy calculation are performed on − for the thiolysis of PTH by CH3S . Earlier theoreti- the M062X/6-31++G(d,p) optimized geometries at MP2/6- cal study 60 reported that there is a significant attack at 311++G(2d,2p) level. 67–72 Unless stated otherwise, energy Solvolysis of organophosphorus pesticide with nucleophiles 1303 values reported herein include zero-point vibrational energy. soft acid base (HSAB) principle of Gazquez85 and Nguyen The aqueous-phase calculations have been performed at the et al., 86 and after a generalization by Ponti, 87 it emerges that same level of theory similar to the gas-phase ones. The sol- the interaction between reaction partners is favored when it vent effect has been taken into account by the conductor-like occurs through minimal s(r), where, 73 screening solvation model (COSMO), as implemented in + − s(r) =|s (r) − s (r)| Gaussian 09. 74 All COSMO calculations in this study have been performed by using default choice of the Gaussian 09 program with the recommended standard parameters. In our 3. Results and Discussion previous study 19 of nucleophilic degradation of fenitroth- ion by anionic nucleophiles, it has been reported that the 3.1 Conformational search of PTH COSMO model provides more appropriate information on energetics than the integral equation formalism polarized Low-energy conformations of molecular systems are continuum solvation model (IEF-PCM) 75–79 when compared found through conformational search by varying geo- with the experimental findings. Therefore, we have used the metric parameters. Among the 25 selected conformers aqueous-phase energetics as well as geometric parameters of of PTH, the minima are finally optimized at M062X/6- the COSMO model in successive discussions. All electronic 31++G(d,p) level of theory. After DFT optimization, we structure calculations have been performed using the Gaus- get five lower-energy structures. The optimized geome- sian 09 suite of quantum chemistry programs. 74 tries are shown in Figure 1, with the relative energies The reactivity descriptor, which is used in this work, is related to the local softness, s(r), the global softness, S of given in the parentheses. The lowest energy conformer the system 80 and the Fukui function, 81 f(r) by the following found here is M1. The most stable ground state struc- relation: tures obtained by us are in good agreement with that reported previously by Ford-Green et al., 88 who carried ( ) = ( ) s r f r S out the conformational search for PTH using DFT. Using the finite difference approximation, as proposed by Yang and Mortier, 82 the compact form of Fukui functions 3.2 Reaction mechanisms: solvolysis of PTH of an atom k in a molecule with N electrons are, + The nucleophilic substitution reactions for solvolysis of f =[q (N + 1) − q (N)] for nucleophilic attack k k k phosphotriesters compounds illustrate that substitution − =[ ( ) − ( − )] fk qk N qk N 1 for electrophilic attack reactions mainly occur at phosphorus center through associative pathway. 89,90 There are two possibilities in where qk(N), qk(N − 1), qk(N + 1), are the electron popula- tion for the atom k in the neutral, cationic, anionic states of a associative pathway: one is addition-elimination mecha- system, respectively. The electron populations are evaluated nism and other is direct displacement pathway. Accord- at the neutral molecular geometry in both gas and aqueous ing to earlier results, 16,17,55,91,92 the multistep pathway phases. The Mulliken 83 charges calculated using the opti- occurs in sarin, acephate, O,S-dimethyl methylphos- mized geometries at M062X/6-31++G(d,p) level is utilized phonothiolate (VX type model compound) and tabun, to compute fk. The finite difference formula is used to calcu- but paraoxon, parathion, fenitrothion and demeton-S late the global softness (S), i.e., appear to be hydrolyzed through direct displacement S = 1/(IE − EA) pathway. In contrast to these results, we have obtained only the multistep pathway for our systems, although where, IE and EA are the first vertical ionization energy and several attempts have been made to find out the sin- electron affinity, respectively. IE and EA can be estimated gle step pathway. Previous theoretical work 16 based on from the Koopman’s theorem 84 as the energy of the highest occupied molecular orbital (HOMO) and the lowest unoccu- direct displacement type mechanisms was reported for pied molecular orbital (LUMO), respectively. Thus, S can be the hydrolysis of a number of phosphotriesterase (PTE) calculated from the following the expression, substrates, but the solvolysis of parathion based on addition-elimination type mechanisms by α-nuclophiles = /( − ) S 1 ELUMO EHOMO has not yet been investigated theoretically. In the present The local softness of nucleophile and electrophile can be investigation, the solvolysis of PTH is mainly due to defined as, the addition-elimination mechanism, involving the for- − − mation of trigonal bipyramidal intermediate and its s = Sf k k further decomposition. The schematic representation of + = + sk Sfk the aforesaid mechanism is as follows (Scheme 1). In our study, s+ represents the local softness of the P atom in Nucleophilic attack on the P atom of PTH followed PTH and s− represents the local softness of the O/S atom in by the release of p-nitrophenol through the SN2-type nucleophiles. From the local version of the well-known hard reaction is consistent with the addition-elimination 1304 Chandan Sahu and Abhijit K Das
Figure 1. Gas-phase-optimized geometries at M062X/6-31++G(d,p) level [(Relative energies in kcal mol−1 are given the parenthesis)] for the most stable conformers of PTH with important geometrical parameters. mechanism. Depending on the relative orientation of free energies and enthalpies calculated at M062X/6- incoming nucleophilic group, there are two possible 31++G(d,p) level are given in Table 1. pathways. In the first pathway, proton of hydroxylamine and hydroperoxide and methyl group of CH S− are in 3 3.3 Reaction of PTH with hydroxylamine anion the opposite direction to the phosphoryl sulfur atom (C1 − path), whereas in second pathway, it is positioned in the (NH2O ) same way (C2 path). All possible pathways have been First we have studied solvolysis mechanism of PTH with explored in both gas and aqueous phases. Interestingly, − NH2O in both gas and aqueous phases. The optimized second pathway is not established for gas phase. To con- structures are given in Figure 2 and the free energy pro- firm this finding, we have scanned the potential energy files obtained from M062X/6-31++G(d,p) calculation surface (PES) (given in the Supplementary Information are presented in Figure 3. All the coordinates associated [SI]) in the gas phase for the P-Nu bond formation and with the stationary points and E values calculated at it is observed that there is a steady increase in energy, M062X/6-31++G(d,p) are given in the SI. having no local minimum or maximum that might sug- In the first step of gas phase solvolysis, a reaction gest absence of any transition state. However, we have complex (RCa) is formed, having 38.0 kcal mol−1 lower bestowed light on only one pathway where orientation in free energy than the infinitely separated complex − of nucleophile is opposite to phosphoryl sulfur atom (PTH + NH2O ). Here, orientation of nucleophile is in the gas phase and established both pathways in the opposite to the phosphoryl sulfur atom. The RCa is aqueous phase. To facilitate the identification of each stabilized through the formation of a hydrogen bond stationary point, we adopt a nomenclature to character- between the O atom of the nucleophile and H atom of ize them in the following sections. In Figures 2 to 7, both –OC2H5 moiety with distances of 1.93 and 1.94 Å. a, b, c in the entire species name indicate the nucle- Free energy profile (Figure 3) shows that nucleophile − − − ophiles NH2O , HOO ,CH3S , respectively, and there attacks at the P center of PTH to form an intermedi- are two pathways associated with two relative orien- ate, IMa, through a unique trigonal bipyramidal (TBP) tations in aqueous phase, one is named as C1aq and transition state, TS1a, having activation energy ( E) −1 other is C2aq. The relative energies calculated at MP2/6- of 6.9 kcal mol and activation free energy ( G) 311++G(2d,2p) // M062X/6-31++G(d,p) level, relative of 6.8 kcal mol−1. Our calculated values of both E Solvolysis of organophosphorus pesticide with nucleophiles 1305
S S S OC2H5 OC2H5 OC2H5 ArO P ArO P Nu + Nu ArO P Nu
OC2H5 OC2H5 OC2H5 Reactant Transition State Intermediate
S S OC2H5 OC2H5 OAr + P ArO P Nu Nu C2H5O OC2H5
Product Transition State
Scheme 1. Addition-elimination reaction mechanism for the nucleophilic solvolysis of PTH. reaction of solvolysis of PTH.
and G are comparable with that reported 92 for the G values for both TS are almost similar but lower alkaline hydrolysis of paraoxon ( E = 8.7, G = than the gas phase values. This is probably due to more 8.6 kcal mol−1), methyl parathion ( E = 8.7, G = α-effect of the nucleophile in the solvent. In both path- 8.8 kcal mol−1) and fenitrothion ( E = 9.0, G = ways for two relative orientations, the product formation 9.3 kcal mol−1) at MP2/6-311++G(2d,2p) // B3LYP/6- is exergoinc in nature ( G of Pa-C1aq = −23.3 and G 31+G(d) level of theory. In the next step, -OAr group is of Pa-C2aq = −22.9 kcal mol−1). The normal-mode eliminated from IMa, which is a penta-coordinated inter- analysis and the IRC calculation confirmed that both mediate with a unique trigonal bipyramidal structure, to transition states are well connected by the intermediate form the product, Pa, situated 73.9 kcal mol−1 below in and product. The aforementioned reaction mechanisms free energy profile, via the transition state, TS2a, having are in good agreement with the addition-elimination activation free energy of 0.7 kcal mol−1. The vibrational mechanism in spite of different orientation of incom- mode of the single imaginary frequency of TS2a corre- ing nucleophilic group. As the activation energy barrier sponds to the complete rupture of the P-O bond. This for elimination of leaving group is lower compared to value is also comparable to the earlier reported results 92 addition of nucleophile group, the formation of TBP for alkaline hydrolysis of paraoxon, methyl parathion, intermediate is the rate-limiting step and this outcome 92 fenitrothion. The reason for this lower energy barrier is agrees well with the previous results. due to utmost similarity in structure between IMa and TS2a and it is consistent with the previous studies of 3.4 Reaction of PTH with hydroperoxide anion 92 paraoxon hydrolysis. In addition, from the enthalpy (HOO−) values, it is observed that formation of Pa is also exother- mic in nature. Hydroperoxide anion is also an efficient α-nucleophile Figure 3 represents the PES of the reaction mech- for solvolysis of PTH and it follows the similar type − − anisms of PTH with NH2O at M062X/6-31++G(d,p) mechanism as described above for NH2O . The opti- level in the aqueous phase. Here, two alternative path- mized stationary points and free energy profile for ways are discussed. Both pathways are started from the both phases are presented in Figure 4 and Figure 5, same substrate but orientation of incoming nucleophile respectively. Figure 5 shows reactant complex (RCb), is different. Figure 2 shows a significant structural differ- which is located 36.9 kcal mol−1 below the reactants ence between RCa-C1aq and RCa-C2aq complexes. But (PTH + HOO−), transforms further into intermediate, nucleophilic attack occurs in both the stationary points IMb, via a TBP geometric transition state, TS1b, having via TS1a-C1aq and TS1a-C2aq to form an intermedi- activation energy ( E) of 7.6 kcal mol−1 and activation ate IMa-C1aq and IMa-C2aq, respectively. The Eand free energy ( G) of 7.4 kcal mol−1. These values are 1306 Chandan Sahu and Abhijit K Das
Table 1. Relative energies ( E) calculated at MP2/6-311++G(2d,2p) // M062X/6- 31++G(d,p) level of theory and H and G calculated at M062X/6-31++G(d,p) level for all the species involved in the reactions.
Gas Aqueous Species E H G Species E H G
− − PTH + NH2O 0.0 0.0 0.0 PTHaq +NH2O aq 0.0 0.0 0.0 RCa −37.6 −49.6 −38.0 RCa-C1aq −2.8 −1.5 8.1 TS1a −30.7 −45.8 −31.2 TS1a-C1aq 1.9 0.2 11.7 IMa −55.1 −72.8 −58.7 IMa-C1aq −16.7 −22.7 −9.1 TS2a −55.3 −71.7 −58.0 TS2a-C1aq −17.0 −22.8 −8.4 Pa −71.5 −86.3 −73.9 Pa-C1aq −31.8 −34.2 −23.3 RCa-C2aq −2.5 −1.1 8.6 − PTH + HOO 0.0 0.0 0.0 TS1a-C2aq 2.2 0.3 13.3 RCb −36.3 −49.1 −36.9 IMa-C2aq −17.4 −23.4 −8.4 TS1b −28.7 −43.4 −29.5 TS2a-C2aq −17.6 −23.5 −8.2 IMb −53.2 −70.5 −55.9 Pa-C2aq −32.0 −34.6 −22.9 − − − + − TS2b 52.1 67.5 54.3 PTHaq HOOaq 0.0 0.0 0.0 Pb −59.1 −73.7 −61.6 RCb-C1aq −3.8 −2.3 7.9 TS1b-C1aq 1.3 −0.3 11.5 − PTH + CH3S 0.0 0.0 0.0 IMb-C1aq −11.3 −16.5 −2.2 RCc −25.4 −34.6 −23.6 TS2b-C1aq −10.8 −15.2 −1.4 TS1c −18.9 −29.8 −16.6 Pb-C1aq −24.1 −25.1 −13.8 IMc −26.2 −35.9 −24.0 RCb-C2aq −3.3 −1.8 8.1 TS2c −26.3 −35.9 −23.3 TS1b-C2aq 1.0 −0.9 12.2 Pc −41.4 −49.3 −39.0 IMb-C2aq −14.3 −19.1 −4.4 TS2b-C2aq −11.6 −16.0 −1.3 Pb-C2aq −21.1 −23.4 −11.6 + − PTHaq CH3Saq 0.0 0.0 0.0 RCc-C1aq −4.4 −1.5 6.7 TS1c-C1aq 6.8 7.1 18.7 IMc-C1aq 1.8 3.6 15.2 TS2c-C1aq 3.5 4.6 17.1 Pc-C1aq −15.2 −11.6 −1.9 RCc-C2aq −3.4 −1.3 9.4 TS1c-C2aq 8.1 8.0 20.2 IMc-C2aq 2.6 4.5 17.5 TS2c-C2aq 3.0 4.6 18.3 Pc-C2aq −12.3 −9.2 −0.5
All energies are in kcal mol−1. comparable with previously reported results 92 of alka- group, -OAr is eliminated to form the product com- line hydrolysis of organophosphorus pesticides. In the plexes, Pb-C1aq and Pb-C2aq via the transition state, next step, a bond cleavage is occurred in the TS, TS2b, TS2b-C1aq and TS2b-C2aq. The IRC calculation for both (activation free energy ( G), 1.6 kcal mol−1)toform the TS corroborates with the leaving of –OAr group the product, Pb. IRC calculation for the corresponding through the intermediate and the product complex. Here, TS confirms the formation of final product via complete formation of the intermediate is rate-determining step breaking of the P-O bond. and E values for corresponding TS is lower in aque- In case of aqueous phase, two possible pathways are ous phase compared to gas phase due to more α-effect explored. Both the reactant complexes, RCb-C1aq and in solvent. RCb-C2aq are lower in free energy in the energy profile, and nucleophilic attack occurs irrespective of orienta- 3.5 Reaction of PTH with methanethiolate anion tions of the PTH and forms unique trigonal bipyramidal − (CH3S ) intermediate structures through the TS1b-C1aq ( G = . −1 = . −1 3 6 kcal mol ) and TS1b-C2aq ( G 4 1 kcal mol ), In contrast to the above-mentioned α-nucleophiles, respectively. Next in the elimination process of leaving − CH3S is also an efficient simple nucleophile for the Solvolysis of organophosphorus pesticide with nucleophiles 1307
Figure 2. Optimized geometries with important parameters in the gas and aqueous phases obtained at − M062X/6-31++G(d,p) level for all the species involved in the reaction of parathion with NH2O . All bond lengths and angels are in Å and (o), respectively. detoxification of PTH and other organophosphate esters. An adduct (RCc), with 23.6 kcal mol−1 lower in free The optimized geometries in gas and aqueous phases energy with respect to isolated reactants, is formed are given in Figure 6 and the potential energy pro- during the initial nucleophilic approach towards PTH files obtained from M062X/6-31++G(d,p) calculation in gas phase. The RCc once formed transforms into are presented in Figure 7. another intermediate, IMc, through a TBP geometric 1308 Chandan Sahu and Abhijit K Das
− Figure 3. Free energy profile for the reaction of parathion with NH2O at M062X/6-31++G(d,p) level of theory in gas (black color line) and aqueous (violet and cyan color lines) phases. C1 (violet line) and C2 (cyan line) correspond to the relative orientation in aqueous phase.
transition state, TS1c, having activation energy ( E) tive orientations of reactant complex, RCc-C1aq and −1 of 6.5 kcal mol and activation free energy ( Gof RCc-C2aq. But nucleophilic attack at P center occurs in −1 7.0 kcal mol . After surmounting the transition state, both the stationary points via TS1c-C1aq and TS1c-C2aq TS1c, the intermediate, IMc, a penta-coordinated inter- to form the intermediates, IMc-C1aq and IMc-C2aq, mediate with a TBP structure, is formed. In the very respectively. The G values for both TS are 12.0 and next step, p-nitrophenol can be eliminated from IMc 10.8 kcal mol−1 but higher than the gas phase values. to form the product, Pc, via the transition state, TS2c, This is probably due to the fact that the ionic nucle- having activation free energy of 0.7 kcal mol−1.The ophile is more stabilized in solvent than in gas phase. vibrational mode of the single imaginary frequency of In the elimination step, the leaving group –OAr is TS2c corresponds to the breaking of P-O bond. The eliminated from both intermediates having orientations greater similarity between intermediate and transition through TS2c-C1aq and TS2c-C2aq to form solvolyzed state structures is responsible for lowering of the free products. The normal-mode analysis and the IRC cal- energy barrier. culation confirmed that both transition states are well − The aqueous phase reaction of PTH with CH3S connected between the intermediate and product struc- passes through the formation of two different rela- tures. The afore-mentioned reaction mechanisms are in Solvolysis of organophosphorus pesticide with nucleophiles 1309
Figure 4. Optimized geometries with important parameters in the gas and aqueous phases obtained at M062X/6-31++G(d,p) level for all the species involved in the reaction of parathion with HOO−. All bond lengths are in Å.
− good agreement with addition-elimination mechanism attack is the rate-determining step. As CH3S is a sim- in spite of the orientation of different incoming nucle- ple nucleophile, there is no possibility of any enhanced ophilic group. Moreover, the first step of nucleophilic nucleophilic effect like other α-nucleophiles. Though 1310 Chandan Sahu and Abhijit K Das
Figure 5. Free energy profile for the reaction of parathion with HOO− at M062X/6-31++G(d,p) level of theory in gas (black color line) and aqueous (violet and cyan color lines) phases. C1 (violet line) and C2 (cyan line) correspond to the relative orientation in aqueous phase.
− CH3S is a simple nucleophile, it also shows multistep in both gas and aqueous phases. The calculated G addition-elimination pathways. value for TS1c is 0.4 kcal mol−1 lower compared to the corresponding transition state TS1b and 0.2 kcal mol−1 higher than the TS1a in gas phase. But from these 3.6 Reactivity trend of three nucleophiles relative energetic values ( G), the trend in reactivity cannot be predicted in gas phase. In case of aqueous We now discuss about the nucleophilic reactivity based phase, G is similar for TS1a-C1 and TS1b-C1 but energetic and conceptual parameters values. Nucle- aq aq ophilicity can be correlated by comparing the rate both are lower in G than TS1c-C1aq. The order of determining transition state structure for different nucle- reactivity for nucleophilic attack for C1aq can be pre- ophiles. It is clear from Table 1 that the first transition dicted from the aqueous phase energetics data and it is − ≈ − > − state of nucleophilic attack is the rate determining state NH2O HOO CH3S , whereas in case of C2aq Solvolysis of organophosphorus pesticide with nucleophiles 1311
Figure 6. Optimized geometries with important parameters in the gas and aqueous phases obtained at − M062X/6-31++G(d,p) level for all the species involved in the reaction of parathion with CH3S . All bond lengths are in Å.
pathway, G value for TS1b-C2aq is lower than that Hand G suggest that the proposed mechanism for of TS1a-C2aq and TS1c-C2aq. Hence, reactivity trend the solvolysis of PTH in both phases is realistic. − − − should be HOO > NH2O > CH3S . So, the reac- Conceptual DFT analysis (by means of electron den- tivity trend is different in aqueous phase and it may be sity as the fundamental property) has been performed to concluded that there is significant enhancement of α- find the trend of reactivity of α-nucleophiles and simple effect for α-nucleophiles in aqueous medium compared nucleophiles with PTH in both phases and try to corre- to simple nucleophiles. But the reasonable values of E, late it with the reactivity trend obtained from Hand 1312 Chandan Sahu and Abhijit K Das
− Figure 7. Free energy profile for the reaction of parathion with CH3S at M062X/6-31++G(d,p) level of theory in gas (black color line) and aqueous (violet and cyan color lines) phases. C1 (violet line) and C2 (cyan line) correspond to the relative orientation in aqueous phase.