I found it surprising that the trichloromethyl group could be chemically reduced into a methyl group quite rapidly in the presence of thiophenol, but once again a failed reaction in the lab gave us the opportunity to learn some nuances about the chemical reactivity of organic compounds. Even more surprising was the fact that this reduction occured through a mechanism in which chlorine atoms behave as electrophiles and not as nucleophiles.
We proposed the mechanism shown in figure 1 to be consistent with the 1H-NMR kinetic experiment (Figure 2) which shows the presence of the intermediary sulfides and leads to the observed phenyl-disulfide as the only isolable byproduct. The proposed mechanism invokes the presence of σ-holes on chlorine atoms to justify the attack of thiophenolate towards the chlorine atom leaving a carbanion behind during the first step. The NMR spectra were recorded at 195K which implies that the energy barriers had to be very low; the first step has a ~3kcal/mol energy barrier at this temperature.
To calculate these energy barriers we employed the BMK functional as implemented in Gaussian09. This functional came highly recommended to this purpose and I gotta say it delivered! The optimized geometries of all transition states and intermediaries were then taken to an MP2 single point upon which the maximum electrostatic potential on each atom (Vmax) was calculated with MultiWFN. In figure 3 we can observe the position and Vmax value of σ-holes on chlorine atoms as suggested by the mapping of electrostatic potential on the electron density of various compounds.
We later ran the same MP2 calculations on other CCl3 groups and found that the binding to an electron withdrawing group is necessary for a σ-hole to be present. (This fact was already present in the literature, of course, but reproducing it served us to validate our methodology.)
We are pleased to have this work published in PhysChemChemPhys. Thanks to Dr. Moisés Romero for letting us into his laboratory and to Guillermo Caballero for his hard work both in the lab and behind the computer; Guillermo is now bound to Cambridge to get his PhD, we wish him every success possible in his new job and hope to see him again in a few years, I’m sure he will make a good job at his new laboratory.
A few weeks back we wrote about using WFN(X) files with MultiWFN in order to find σ-holes in halogen atoms by calculating the maximum potential on a given surface. We later found out that using a chk file to generate a wfn(x) file using the guess=(read,only) keyword didn’t retrieve the MP2 wavefunction but rather the HF wavefunction! Luckily we realized this problem very quickly and were able to fix it. We tried to generate the wfn(x) file with the following keywords at the route section
#p guess=(read,only) density=current
but we kept retrieving the HF values, which we noticed by running the corresponding HF calculation and noticing that every value extracted from the WFN file was exactly the same.
So, if you want a WFN(X) file for post processing an MP2 (or any other post-HartreFock calculation for that matter) ask for it from the beginning of your calculation in the same job. I still don’t know how to work around this or but will be happy to report it whenever I do.
PS. A sincere apology to all subscribers for getting a notification to this post when it wasn’t still finished.
Some atomic properties such as an atomic charge are isotropic, but every now and then some derivations of them become anisotropic, for example the plotting of the Molecular Electrostatic Potential (MEP) on the electron density surface can exhibit some anisotropic behavior; quantifying it can be a bit challenging.
It is well known that halogen atoms such as Chlorine can form so-called halogen-bonds of the type R-Cl-R in crystals with a near perfect 180° angle. This finding has lead to the discovery of σ-holes in halogens. σ-Holes are electrophilic portions of the anisotropic electrostatic potential in an otherwise nucleophilic atom. Recently, Guillermo “Memo” Caballero and I calculated the MEP for a series of trichloromethyl-containing compounds at the MP2/cc-pVQZ level of theory and the mapping shows evidence of such σ-holes as seen in Figure1. Those small blue portions on an otherwise red atom indicate that some electron density is missing in that position, which by the way is located at 180° away from the carbon atom.
But having the picture is not enough. We want to quantify just how strong are those σ-holes to effectively attract a nucleophile and perhaps perform some chemistry on the C-Cl bond. That’s when we resorted to MultiWFN, a Multifunctional Wavefunction Analyzer developed by Tian Lu (卢天) at the Beijing Kein Research Center for Natural Sciences. You can check the project leader list of publications here. Among many other capabilities, MultiWFN is able to print details about properties along a surface.
In order to work with MultiWFN you need to generate a *.wfn file, if you have a previous Gaussian calculation for which you want to analyze their surface you can run a guess=only calculation in order to extract the wavefunction from the checkpoint file. Here is a dummy of the input for such calculation
%chk=oldfile.chk # output=wfn geom=check guess=(only,read) density=current --blank-- Title Card --blank-- 0 1 --blank-- filename.wfn --blank--
In our case, having a post-Hartree-Fock calculation, the use of density=current is mandatory to get the MP2 density matrix and not just the HF one. Running this calculation will generate the file filename.wfn which is now used with MultiWFN. When starting MultiWFN you get to see a terminal window like the one below in which you are asked to input the path of your wfn file:
After loading it you will get the following window with the various options available. Type 12// (these two slashes are mandatory) to get the quantitative analysis of molecular surface option.
Then you will be asked to define some elements of that surface (we used the default options 0)
The following screenshot shows the results section in which several maxima and minima of electrostatic potential were found (7 and 11 in our case); a star is placed on the side of the global maximum. The value of the MEP at those points is given in Hartrees, eV and kcal/mol which I personally hate because there isn’t such a thing as a mole of ‘potentials’ (same argument as giving an orbital’s energy in kcal/mol, moles of what? orbitals? Personally, I don’t like it even if its valid).
Their visualizer is activated through the option 0 and although it is far from pretty it is quite good enough to find the numbers corresponding to maxima and minima of the MEP on the isodensity surface. If we look for the maxima then we find for our example (CHCl3) that a maximum is located in front of each Cl atom in a straight line from the C atom. Now we get to put a number on the mapped isosurface provided by Gaussian or even import the file into Chimera.
We are still working our way around MultiWFN, I hope we can find the batch option, it would be most useful. In the mean time, Guillermo and I will continue to search for σ-holes in chlorinated reagents. Thanks to Guillermo for his ongoing work in this and other topics within the realm of organic reactivity.
Have you any suggestions or ideas to work with MultiWFN? Please share them!