Theoretical study on the interactions of halogen‐bonds and pnicogen‐bonds in phosphine derivatives with Br2, BrCl,and BrF

The MP2 ab initio quantum chemistry methods were utilized to study the halogen‐bond and pnicogen‐bond system formed between PH₂X (X = Br, CH₃, OH, CN, NO₂, CF₃) and BrY (Y = Br, Cl, F). Calculated results show that all substituent can form halogen‐bond complexes while part substituent can form pnicogen‐bond complexes. Traditional, chlorine‐shared and ion‐pair halogen‐bonds complexes have been found with the different substituent X and Y. The halogen‐bonds are stronger than the related pnicogen‐bonds. For halogen‐bonds, strongly electronegative substituents which are connected to the Lewis acid can strengthen the bonds and significantly influenced the structures and properties of the compounds. In contrast, the substituents which connected to the Lewis bases can produce opposite effects. The interaction energies of halogen‐bonds are 2.56 to 32.06 kcal·mol^?1; The strongest halogen‐bond was found in the complex of PH₂OH???BrF. The interaction energies of pnicogen‐bonds are in the range 1.20 to 2.28 kcal·mol^?1; the strongest pnicogen‐bond was found in PH₂Br???Br₂ complex. The charge transfer of lp(P) ? σ*(Br? Y), lp(F) ? σ*(Br? P), and lp(Br) ? σ*(X? P) play important roles in the formation of the halogen‐bonds and pnicogen‐bonds, which lead to polarization of the monomers. The polarization caused by the halogen‐bond is more obvious than that by the pnicogen‐bond, resulting in that some halogen‐bonds having little covalent character. The symmetry adapted perturbation theory (SAPT) energy decomposition analysis showes that the halogen‐bond and pnicogen‐bond interactions are predominantly electrostatic and dispersion, respectively.     

Coupling between Substituents as a Function of Cage Structure: Synthesis and Valence Ionized States of Bridgehead Disubstituted Parent and Hexafluorinated Bicyclo[1.1.1]pentane Derivatives C5X6Y2

He(I) photoelectron spectroscopy was used to examine the valence‐shell electronic structure of three new and seven previously known bicyclo[1.1.1]pentane derivatives, 1,3‐Y₂‐C₅X₆ (for X=H, Y=H, Cl, Br, I, CN; for X=F, Y=H, Br, I, CN). A larger series (X=H or F, Y=H, F, Cl, Br, I, At, CN) has been studied computationally with the SAC‐CI (symmetry‐adapted cluster configuration interaction) method. The outer‐valence ionization spectra calculated by the SAC‐CI method, including spin–orbit interaction, reproduced the experimental photoelectron spectra well, and quantitative assignments are given. When the extent of effective through‐cage interaction between the bridgehead halogen lone‐pair orbitals was defined in the usual way by orbital‐energy splitting, it was found to be larger than that mediated by other cages such as cubane, and was further enhanced by hexafluorination. The origin of the orbital‐energy splitting is analyzed in terms of cage structure, and it is pointed out that its relation to the degree of interaction between the bridgehead substituents is not as simple as is often assumed.     

Organometallated Halonium Salts of the Type [{Cp(CO)2Fe}2X]SbY6 (X = Cl,Br, I; Y = F,Cl)

The reaction of CpFe(CO)₂X (X = Cl, Br, I) with SbY₅ (Y = F, Cl) in toluene leads to the cationic, halogen‐bridged compounds [{Cp(CO)₂Fe}₂X]SbY₆ ( 1 – 6 ). The halide of CpFe(CO)₂X is eliminated by the Lewis acid SbY₅, and the fragment “CpFe(CO)₂⁺” reacts with further CpFe(CO)₂X to form the halogen bridge between both the organometallic substituents. The exclusive formation of the counter anion SbY₆^– is caused by the oxidizing action of the antimony pentahalides, by which SbY₃ and the interhalogens XY are always obtained. The compounds have been characterized by their NMR‐, IR‐ and Mass spectra, the compounds 1 – 3 and 6 additionally by single crystal structure analyses. They show decreasing bond angles Fe–X–Fe following the range Cl → Br → I and the VSEPR concept; the two CpFe(CO)₂ groups are staggered with the dihedral angle Cp(centre)–Fe–Fe–Cp(centre) of about 160°.     

Some measures for mediating the strengths of halogen bonds with the B–B bond in diborane(4) as an unconventional halogen acceptor

A series of halogen‐bonded complexes with diborane(4) 1 and its derivatives (Li 2 , methyl 3 , CN 4 ) as the halogen acceptors as well as with XCN, XCCH, XCF₃, XF (X = Cl, Br, I) as the halogen donors have been investigated by means of quantum chemical calculations at the MP2/aug‐cc‐pVTZ level. The result shows that the B?B bond is a good electron donor in halogen bonding, particularly for the halogen donor XF. Interestingly, for the halogen donor XF, the halogen bond becomes stronger in order of IF < BrF < ClF. It is found that the addition of electron‐donating substituents greatly strengthens the halogen bonding interaction to the point where it exceeds that of the majority of H‐bonds. When the N atom in 2 ‐BrCN is combined with another interaction, its strength has a further increase due to the cooperative effect. This study combines the boron compounds with halogen bonds and would be significant for expanding their applied fields. © 2013 Wiley Periodicals, Inc.     

Weak hydrogen and halogen bonding in 4‐[(2,2‐difluoroethoxy)methyl]pyridinium iodide and 4‐[(3‐chloro‐2,2,3,3‐tetrafluoropropoxy)methyl]pyridinium iodide

To enable a comparison between a C—H…X hydrogen bond and a halogen bond, the structures of two fluorous‐substituted pyridinium iodide salts have been determined. 4‐[(2,2‐Difluoroethoxy)methyl]pyridinium iodide, C₈H₁₀F₂NO⁺·I⁻, (1), has a –CH₂OCH₂CF₂H substituent at the para position of the pyridinium ring and 4‐[(3‐chloro‐2,2,3,3‐tetrafluoropropoxy)methyl]pyridinium iodide, C₉H₉ClF₄NO⁺·I⁻, (2), has a –CH₂OCH₂CF₂CF₂Cl substituent at the para position of the pyridinium ring. In salt (1), the iodide anion is involved in one N—H…I and three C—H…I hydrogen bonds, which, together with C—H…F hydrogen bonds, link the cations and anions into a three‐dimensional network. For salt (2), the iodide anion is involved in one N—H…I hydrogen bond, two C—H…I hydrogen bonds and one C—Cl…I halogen bond; additional C—H…F and C—F…F interactions link the cations and anions into a three‐dimensional arrangement.     

Is a MH (M = Be and Mg) radical a better electron donor in halogen‐hydride interaction?: A theoretical comparison with HMH

Quantum chemical calculations have been performed to investigate the complexes of HMH? XCCH, HMH? XCF₃, MH? XCCH, and MH? XCF₃ (M = Be and Mg; X = Cl, Br, and I) at the MP2/aug‐cc‐pVTZ level. The geometrical, energetic, and spectroscopic parameters were analyzed for these complexes. The results show that the MH is a better electron donor in the halogen‐hydride interaction than HMH. The enhancement of halogen‐hydride interaction increases in the order of Cl < Br < I, Be < Mg, and XCCH < XCF₃. The halogen‐hydride interaction was understood with natural bond orbital, atoms in molecules, and electrostatic potentials. © 2012 Wiley Periodicals, Inc.     

Magnitude and origin of the attraction and directionality of the halogen bonds of the complexes of C6F5X and C6H5X (X = I, Br, Cl and F) with pyridine

The geometries and interaction energies of complexes of pyridine with C₆F₅X, C₆H₅X (X=I, Br, Cl, F and H) and R_FI (R_F=CF₃, C₂F₅ and C₃F₇) have been studied by ab initio molecular orbital calculations. The CCSD(T) interaction energies (E_int) for the C₆F₅X–pyridine (X=I, Br, Cl, F and H) complexes at the basis set limit were estimated to be ?5.59, ?4.06, ?2.78, ?0.19 and ?4.37 kcal mol^?1, respectively, whereas the E_int values for the C₆H₅X–pyridine (X=I, Br, Cl and H) complexes were estimated to be ?3.27, ?2.17, ?1.23 and ?1.78 kcal mol^?1, respectively. Electrostatic interactions are the cause of the halogen dependence of the interaction energies and the enhancement of the attraction by the fluorine atoms in C₆F₅X. The values of E_int estimated for the R_FI–pyridine (R_F=CF₃, C₂F₅ and C₃F₇) complexes (?5.14, ?5.38 and ?5.44 kcal mol^?1, respectively) are close to that for the C₆F₅I–pyridine complex. Electrostatic interactions are the major source of the attraction in the strong halogen bond although induction and dispersion interactions also contribute to the attraction. Short‐range (charge‐transfer) interactions do not contribute significantly to the attraction. The magnitude of the directionality of the halogen bond correlates with the magnitude of the attraction. Electrostatic interactions are mainly responsible for the directionality of the halogen bond. The directionality of halogen bonds involving iodine and bromine is high, whereas that of chlorine is low and that of fluorine is negligible. The directionality of the halogen bonds in the C₆F₅I– and C₂F₅I–pyridine complexes is higher than that in the hydrogen bonds in the water dimer and water–formaldehyde complex. The calculations suggest that the C? I and C? Br halogen bonds play an important role in controlling the structures of molecular assemblies, that the C? Cl bonds play a less important role and that C? F bonds have a negligible impact.     

Comparative Computational Study of Hydrogen Abstraction Reactions of CY3H + XO− (X,Y = F,Cl, and Br)

The effects of halogen substituents on the reactivity are characterized by the hybrid B3LYP and BHandHLYP functionals of density functional theory using the aug‐cc‐pVDZ basis set. The species XO⁻ and CY₃H, where X, Y = F, Cl, and Br, have been chosen as model reactants in this work. Also, the mechanism of the hydrogen abstraction (HAT) reaction has been used to study the chemical reactivity of these anionic reactions. Our theoretical findings suggest that the relative reactivity of the CY₃H + XO⁻ reactions increases as Y goes from F to Br and decreases as X goes from F to Br. Moreover, among all reactions investigated in this study, the special role of the Y has very dominant effect on activation of the C–H bond in CY₃H when XO⁻ attacks the CY₃H. Again, through the transition state theory the rate constants at 298–1000 K are also evaluated for the HAT reactions, indicating the lower the temperature the faster is the chemical reaction.     

Free radical addition of cyclopentane and cyclohexane to halogeno derivatives of 1,2-difluoroethene

Free radical addition reactions between cyclopentane and cyclohexane and a range of difluoroalkenes, CF₂CXY (X, Y = H, F, Cl, Br) gave a series of adducts bearing difluoromethylene substituents, R-CF₂-CXYH (R = c-C₅H₉ or c-C₆H₁₁), in reasonable yield even though telomerisation and halogen transfer (when X, Y = Cl, Br) can compete. Dehydrofluorination of the adducts gave several new polyhalogenated alkenes.     

Rational Modification of the Hierarchy of Intermolecular Interactions in Molecular Crystal Structures by Using Tunable Halogen Bonds

A family of 16 isomolecular salts (3‐XpyH)₂[MX′₄] (3‐XpyH=3‐halopyridinium; M=Co, Zn; X=(F), Cl, Br, (I); X′=Cl, Br, I) each containing rigid organic cations and tetrahedral halometallate anions has been prepared and characterized by X‐ray single crystal and/or powder diffraction. Their crystal structures reflect the competition and cooperation between non‐covalent interactions: N? H???X′? M hydrogen bonds, C? X???X′? M halogen bonds and π–π stacking. The latter are essentially unchanged in strength across the series, but both halogen bonds and hydrogen bonds are modified in strength upon changing the halogens involved. Changing the organic halogen (X) from F to I strengthens the C? X???X′? M halogen bonds, whereas an analogous change of the inorganic halogen (X′) weakens both halogen bonds and N? H???X′? M hydrogen bonds. By so tuning the strength of the putative halogen bonds from repulsive to weak to moderately strong attractive interactions, the hierarchy of the interactions has been modified rationally leading to systematic changes in crystal packing. Three classes of crystal structure are obtained. In type A (C? F???X′? M) halogen bonds are absent. The structure is directed by N? H???X′? M hydrogen bonds and π‐stacking interactions. In type B structures, involving small organic halogens (X) and large inorganic halogens (X′), long (weak) C? X???X′? M interactions are observed with type I halogen–halogen interaction geometries (C? X???X′ ≈ X???X′? M ≈155°), but hydrogen bonds still dominate. Thus, minor but quite significant perturbations from the type A structure arise. In type C, involving larger organic halogens (X) and smaller inorganic halogens (X′), stronger halogen bonds are formed with a type II halogen–halogen interaction geometry (C? X???X′ ≈180°; X???X′? M ≈110°) that is electrostatically attractive. The halogen bonds play a major role alongside hydrogen bonds in directing the type C structures, which as a result are quite different from type A and B.     

Coordination Polyhedra OsXn(X = F, Cl, Br, I) in Crystal Structures

The Voronoi–Dirichlet polyhedra (VDP) and the method of intersecting spheres were used to perform crystal-chemical analysis of compounds containing complexes [Os _a X _b ] ^z–(X = F, Cl, Br, I). Atoms of Os(V) at X = F and Cl, of Os(IV) at X = Cl, Br, and of Os(III) at X = Br were found to exhibit a coordination number of 6 with respect to the halogen atoms and to form OsX₆octahedra. The coordination polyhedra of Os(III) for X = Cl, I are square pyramids OsX₄. Each Os(III) atom forms one Os–Os bond; as a consequence, the OsBr₆octahedra share a face in forming Os₂Br^3– ₉complexes, while the OsX₄pyramids (X = Cl, I) dimerize to produce [X₄Os–OsX₄]^2–ions. The influence of the valence state of the Os atoms and of the nature of the halogen atoms on the composition and structure of the complexes formed and some characteristics of the coordination sphere of Os were considered.     

The Cyclohexasilanes Si6H11X and Si6Me11X with X=F,Cl, Br and I: A Quantum Chemical and Raman Spectroscopic Investigation of a Multiple Conformer Problem

The conformers of the monohalocyclohexasilanes, Si₆H₁₁X (X=F, Cl, Br or I) and the haloundecamethylcyclohexasilanes, Si₆Me₁₁X (X=F, Cl, Br or I) are investigated by DFT calculations employing the B3LYP density functional and 6‐31+G* basis sets for elements up to the third row, and SDD basis sets for heavier elements. Five minima are found for Si₆H₁₁X—the axial and equatorial chair conformers, with the substituent X either in an axial or equatorial position—and another three twisted structures. The equatorial chair conformer is the global minimum for the X=Cl, Br and I, the axial chair for X=F. The barrier for the ring inversion is ～13 kJ mol^?1 for all four compounds. Five minima closely related to those of Si₆H₁₁X are found for Si₆Me₁₁X. Again, the equatorial chair is the global minimum for X=Cl, Br and I, and the axial chair for X=F. Additionally, two symmetrical boat conformers are found as local minima on the potential energy surfaces for X=F, Cl and Br, but not for X=I. The barrier for the ring inversion is ～14–16 kJ mol^?1 for all compounds. The conformational equilibria for Si₆Me₁₁X in toluene solution are investigated using temperature dependent Raman spectroscopy. The wavenumber range of the stretching vibrations of the heavy atoms X and Si from 270–370 cm^?1 is analyzed. Using the van′t Hoff relationship, the enthalpy differences between axial and equatorial chair conformers (H_ax?H_eq.) are 1.1 kJ mol^?1 for X=F, and 1.8 to 2.8 kJ mol^?1 for X=Cl, Br and I. Due to rapid interconversion, only a single Raman band originating from the “averaged” twist and boat conformers could be observed. Generally, reasonable agreement between the calculated relative energies and the experimentally determined values is found.     

Tuning the Photophysical Properties of Symmetric Squarylium Dyes: Investigation on the Halogen Modulation Effects

A series of symmetric squarylium dyes (SQDPA-X) with different halogen (X=F, Cl, Br, I) substituents have been developed. The photophysical properties could be facilely tuned by the halogen modulation effects. The strategy of incorporating different halogen substitutions into AIE active luminogens enables a facile approach for exploring new intriguing organic fluorescent dyes.     

Competition between halogen bond and hydrogen bond in complexes of superalkali Li3S and halogenated acetylene XCCH (X = F,Cl, Br,and I)

Complexes of superalkali Li₃S and XCCH (X = F, Cl, Br, and I) have been studied with theoretical calculations at the MP2/aug‐cc‐pVTZ level. Three types of structures are found: (A) the X atom combines with the S atom through a halogen bond; (B) the X atom interacts with the π electron of Li₃S by a π halogen bond; (C) the H atom combines with the S atom through a hydrogen bond. For A and B, a heavier halogen atom makes the interaction stronger, while for C, the change of interaction energy is not obvious, showing a small dependence on the nature of the X atom in HCCX. A is more stable than B and their difference in stability decreases as X varies from Cl to I. For the F and Cl complexes, A is weaker than C, however, the former is stronger than the latter in the Br and I complexes. The above three types of interactions have been analyzed by means of electron localization function, electron density difference, and energy decomposition, and the results show that they have similar nature and features with conventional interactions. © 2014 Wiley Periodicals, Inc.     

Donorstabilisierte Galliumdihalogenide Ga2X4·2D (X = Cl,Br; D = Donormolekül): Ein experimenteller Beitrag zu den Variationsmöglichkeiten der Ga‐Ga‐Einfachbindung

Donor‐stabilized Galliumdihalides Ga₂X₄·2D (X = Cl, Br; D = Donor): An Experimental Contribution on the Variation of the Gallium‐Gallium Single Bond During the disproportionation of metastable GaX‐solutions (X= Cl, Br) donor‐stabilized galliumdihalides are formed as oxidized products. According to X‐ray structure analyses they all exhibit dimeric entities DX₂Ga‐GAX₂D (D= THF, NHEt2, NEt3, 4‐^tButylpyridin or Br^‐), which means these compounds are isoelectronic with ethane and could schematically be regarded as representatives of catenulate or alkane‐like gallium subhalides: Ga_n(X, D)_2n+2. The gallium‐gallium bond in these compounds is shorter than in the organometallic compounds such as R₂Ga‐GaR₂. The comparison of the bonding situation in the galliumdihalides, particularly of the gallium‐gallium bond, shows clearly the influence by donor molecules as well as by halogen ligands.     

Molecular quadrupole moments and magnetic anisotropies of non-dipolar 1,3,5-trisubstituted benzenes

The molecular quadrupole moments and magnetic anisotropies of a series of non-dipolar 1,3,5-trisubstituted benzenes (C₆H₃X₃; X=H, Me, t-Bu, F, Cl, Br) are analysed to obtain information concerning the charge distributions. Whereas alkyl groups have relatively little effect on benzene. halogen substituents are strongly interactive and the moments become progressively less negative from 1.3.5-tribromo-to 1.3.5-trifluoro-benzene.     

On the Oxidation of the Three‐Dimensional Aromatics [B12X12]2− (X=F,Cl, Br,I)

The perhalogenated closo‐dodecaborate dianions [B₁₂X₁₂]^2? (X=H, F, Cl, Br, I) are three‐dimensional counterparts to the two‐dimensional aromatics C₆X₆ (X=H, F, Cl, Br, I). Whereas oxidation of the parent compounds [B₁₂H₁₂]^2? and benzene does not lead to isolable radicals, the perhalogenated analogues can be oxidized by chemical or electrochemical methods to give stable radicals. The chemical oxidation of the closo‐dodecaborate dianions [B₁₂X₁₂]^2? with the strong oxidizer AsF₅ in liquid sulfur dioxide (lSO₂) yielded the corresponding radical anions [B₁₂X₁₂] ^{? ?} (X=F, Cl, Br). The presence of radical ions was proven by EPR and UV/Vis spectroscopy and supported by quantum chemical calculations. Use of an excess amount of the oxidizing agent allowed the synthesis of the neutral perhalogenated hypercloso‐boranes B₁₂X₁₂ (X=Cl, Br). These compounds were characterized by single‐crystal X‐ray diffraction of dark blue B₁₂Cl₁₂ and [Na(SO₂)₆][B₁₂Br₁₂] ? B₁₂Br₁₂. Sublimation of the crude reaction products that contained B₁₂X₁₂ (X=Cl, Br) resulted in pure dark blue B₁₂Cl₁₂ or decomposition to red B₉Br₉, respectively. The energetics of the oxidation processes in the gas phase were calculated by DFT methods at the PBE0/def2‐TZVPP level of theory. They revealed the trend of increasing ionization potentials of the [B₁₂X₁₂]^2? dianions by going from fluorine to bromine as halogen substituent. The oxidation of all [B₁₂X₁₂]^2? dianions was also studied in the gas phase by mass spectrometry in an ion trap. The electrochemical oxidation of the closo‐dodecaborate dianions [B₁₂X₁₂]^2? (X=F, Cl, Br, I) by cyclic and Osteryoung square‐wave voltammetry in liquid sulfur dioxide or acetonitrile showed very good agreement with quantum chemical calculations in the gas phase. For [B₁₂X₁₂]^2? (X=F, Cl, Br) the first and second oxidation processes are detected. Whereas the first process is quasi‐reversible (with oxidation potentials in the range between +1.68 and +2.29 V (lSO₂, versus ferrocene/ferrocenium (Fc^0/+))), the second process is irreversible (with oxidation potentials ranging from +2.63 to +2.71 V (lSO₂, versus Fc^0/+)). [B₁₂I₁₂]^2? showed a complex oxidation behavior in cyclic voltammetry experiments, presumably owing to decomposition of the cluster anion under release of iodide, which also explains the failure to isolate the respective radical by chemical oxidation.     

Acylation Reactions of Dibenzo-7-phosphanorbornadiene: DFT Mechanistic Insights

Extensive DFT calculations provide deep mechanistic insights into the acylation reactions of tert-butyl dibenzo-7-phosphanobornadiene with PhCOX (X=Cl, Br, I, OTf) in CH₂Cl₂ solution. Such reactions are initialized by the nucleophilic P⋅⋅⋅C attack to the carbonyl group to form the acylphosphonium intermediate A⁺ together with X⁻ anion, followed either by nucleophilic X⁻⋅⋅⋅P attack (X=Cl, Br, and I) toward A⁺ to eliminate anthracene or by slow rearrangement or decomposition of A⁺ (X=OTf). In contrast to the first case (X=Cl) that is rate-limited by the initial P⋅⋅⋅C attack, other reactions are rate-limited by the second X⁻⋅⋅⋅P attack for X=Br and I and even thermodynamically prevented for X=OTf, leading to isolable phosphonium salts. The rearrangement of phosphonium A⁺ is initiated by a P-C bond cleavage, followed either by sequential proton-shifts to form anthracenyl acylphosphonium or by deprotonation with additional base Et₃N to form neutral anthracenyl acylphosphine. Our DFT results strongly support the separated acylphosphonium A⁺ as the key reaction intermediate that may be useful for the transfer of acylphosphenium in general.     

Theoretical investigation on identical anionic halide‐exchange SN2 reaction processes on N‐haloammonium cation NH3X+ (X = F,Cl, Br,and I)

CCSD(T) calculations have been used for identically nucleophilic substitution reactions on N‐haloammonium cation, X^? + NH₃X⁺ (X = F, Cl, Br, and I), with comparison of classic anionic S_N2 reactions, X^? + CH₃X. The described S_N2 reactions are characterized to a double curve potential, and separated charged reactants proceed to form transition state through a stronger complexation and a charge neutralization process. For title reactions X^? + NH₃X⁺, charge distributions, geometries, energy barriers, and their correlations have been investigated. Central barriers ΔE for X^? + NH₃X⁺ are found to be lower and lie within a relatively narrow range, decreasing in the following order: Cl (21.1 kJ/mol) > F (19.7 kJ/mol) > Br (10.9 kJ/mol) > I (9.1 kJ/mol). The overall barriers ΔE relative to the reactants are negative for all halogens: ?626.0 kJ/mol (F), ?494.1 kJ/mol (Cl), ?484.9 kJ/mol (Br), and ?458.5 kJ/mol (I). Stability energies of the ion–ion complexes ΔE_comp decrease in the order F (645.6 kJ/mol) > Cl (515.2 kJ/mol) > Br (495.8 kJ/mol) > I (467.6 kJ/mol), and are found to correlate well with halogen Mulliken electronegativities (R² = 0.972) and proton affinity of halogen anions X^? (R² = 0.996). Based on polarizable continuum model, solvent effects have investigated, which indicates solvents, especially polar and protic solvents lower the complexation energy dramatically, due to dually solvated reactant ions, and even character of double well potential in reactions X^? + CH₃X has disappeared. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Density functional theory study on the reactions of X− with CH3SY (X,Y = F,Cl, Br,I)

Gas‐phase anionic reactions X^? + CH₃SY (X, Y = F, Cl, Br, I) have been investigated at the level of B3LYP/6‐311+G (2df,p). Results show that the potential energy surface (PES) of gas‐phase reactions X^? + CH₃SY (X, Y = Cl, Br, I) has a quadruple‐well structure, indicating an addition–elimination (A–E) pathway. The fluorine behaves differently in many respects from the other halogens and the reactions F^? + CH₃SY (Y = F, Cl, Br, I) correspond to deprotonation instead of substitution. The gas‐phase reactions X^? + CH₃SF (X = Cl, Br, I), however, follow an A–E pathway other than the last two out going steps (COM2 and PR) that proceeds via a deprotonation. The polarizable continuum model (PCM) has been used to evaluate the solvent effects on the energetics of the reactions X^? + CH₃SY (X, Y = Cl, Br, I). The PES is predicted to be unimodal in the solvents of high polarity. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007