# Electron Charge Density Distribution from X-Ray Diffraction Study of the 4-Methoxybenzenecarbothioamide Compound

Authors: Mokhtaria Drissi; Abdelkader Chouaih; Youcef Megrouss; Fodil Hamzaoui
Journal: Journal of Crystallography
 (2013)

---

## Abstract

The molecular electron charge density distribution of the title compound is described accurately using the multipolar model of Hansen and Coppens. The net atomic charge and the in-crystal molecular dipole moment have been determined in order to understand the nature of inter- and intramolecular charge transfer.
The study reveals the nature of intermolecular interactions including charge transfer and hydrogen bonds in the title compound. In this crystal, the molecules form dimers via N–H⋯S intermolecular hydrogen bonds. The dimers are further linked by C–H⋯O hydrogen bonds into chains along the c crystallographic axis. This study has also allowed us to determine the electrostatic potential and therefore locate the electropositive part and the electronegative part in molecular scale of the title compound.


---

## Body

## 1. Introduction

The chemistry of thioamides and their derivatives continues to be of particular interest due to their interesting structural features and also due to their biological importance. Heterocyclic compounds of these types exhibit a wide spectrum of biological activity and have extensive practical applications [1–3]. In addition, thioamides are important ligands in the field of coordination chemistry [4]. Due to these positive traits, thioamides have been under study for a long time, but much of their basic chemistry remains unexplored.The heterocyclic 4-methoxybenzenecarbothioamide compound, of chemical formula C8H9NOS, crystallizes in the P212121 space group, with eight molecules in the unit cell. The general features of the structure have been described previously [5]. It contains two independent molecules with the methoxy groups oriented in opposite conformations. The mean planes of the carbothioamide groups are tilted by 7.88(15) and 11.16°(9) from the mean planes of the benzene rings.In this work, the crystallographic data were taken from another article [5] and used for the current study. The structure has been redetermined using the same X-ray diffraction data.So this work focuses on the thermal analysis and electronic charge density study of the 4-methoxybenzenecarbothioamide molecule. The scope of this study is to model aspherical features of the atomic electron density arising from chemical bonding at low temperature of the molecule by transferring the electron density multipole parameters obtained from accurate X-ray diffraction experiments and calculate the molecular dipole moment. In order to understand the nature of inter- and intramolecular charge transfer, the electrostatic potential distribution of the molecule in crystal has been analyzed.


## 2. Materials and Methods

### 2.1. High-Order Spherical Atom Refinement

A high-order spherical atom refinement (HO) on the non-H atoms is performed to get least-biased positional and thermal parameters. Hence, due to the properties of the Fourier transform, this HO procedure refines the positional and anisotropic thermal motion parameters only on core electrons. Because valence electron scattering factors diffuse only at low resolution, HO refinement gives precise positions and displacements of the core electrons only.In order to extract the aspherical features of the electron density, the following strategy was designed. First, a full-anisotropic spherical atom refinement using starting parameters from a Shelxl [6] refinement is performed. Then, the multipole refinement which has been developed by Hansen and Coppens [7] is selected to perform the electron density analysis using the MoPro package software [8].In parallel, the net atomic charge on the different atoms of the title compound and the molecular dipole moment were also estimated by theoretical calculations. Theoretical calculations were performed with the Gaussian package [9] at the HF/6-31G* level of theory and semiempirical methods.A higher-order refinement of X-ray data becomes essential to obtain accurate core positions and the associated thermal parameters of the different atoms of the molecule.The X-ray scattering factors for S, C, N, and O atoms come from the international tables for X-ray crystallography [10] while for H atom the bonded H-atom scattering curve of  Stewart et al. (1965) was used [11]. The H atoms were found by difference Fourier synthesis and their coordinates were adjusted by extending along the C–H and N–H bond directions to bond lengths of 1.105 Å and 1.020 Å, respectively. Any position error of the hydrogen atom will bring about errors strongly correlated with the dipolar population of hydrogen. Refinement details are displayed in Table 1.Table 1
Refinement details.
N
R
w
R
Spherical refinement
13
0.0274
0.03997
Multipole refinement
68
0.02358
0.03165
Wavelength,λ (Å)
0.71073
Completeness (%)
99.4
Resolution, Sinθ/λ (Å−1)
0.8
N is the number of refined parameters. R=∑|Fo|-|Fc|/∑|Fo|; wR=[∑w(|Fo|-|Fc|)2/∑w|Fo|2]1/2.


### 2.2. Charge Density Multipole Model

The Hansen-Coppens multipole formalism, as implemented in the Mopro least-squares program for multipole refinement, was used for both observed and theoretical structure factor fitting. It describes the crystal electron density as a superposition of aspherical pseudo-atoms each modeled on a multipole expansion:(1)ρatom(r→)=ρc(r→)+PVK′3ρv(k′r→)+∑llmax∑m=-1+1K′′Rl(K′′r→)PlmYlm(r→r).Hereρc and ρv are spherically averaged Hartree-Fock core and valence densities, with ρv being normalized to one electron. The Slater-type radial functions Rl(r→)=Nl·rlnlexp(-k′ξlr→) modulate the spherical harmonic angular functions Ylm± and N1 is a normalization factor. The values for parameters n=n1 and ξ were chosen according to rules provided by Coppens (1997) [12].In practice, two charge density variables, the population parameter and a kappa parameter which allows expansion and contraction of the valence shell, are added to the conventional parameters of structure analysis [13]. The population parameters Pv and Plm are floated along with κ′ and κ′′ during the refinement.The sum overm in (1) includes ±l, so that for each one, 2l +1 functions are included. The aspheric charge density was described at the octupole level (l=3) for the atoms C, N, and O and the hexadecapole for the atoms S (l=4) and at the dipole levels (l=1 or 2) for hydrogen atoms not involved and involved in strong H bonds, respectively. Charge densities of all hydrogen were considered to have cylindrical symmetry along the corresponding hydrogen-heavy atom bond.The quality of a refined model can be monitored based on the residuals and the goodness of fit apart from closely inspecting the deformation density maps. This model provides the advantage of using constraints based on assumed local symmetry on the atomic sites to limit the number of least-squares variables. Figure1 gives the local site symmetry of the different atoms of the molecule.Figure 1
Labeling of the atoms and definition of local orthogonal reference axes for the atom-centered multipolar functions.


## 2.1. High-Order Spherical Atom Refinement

A high-order spherical atom refinement (HO) on the non-H atoms is performed to get least-biased positional and thermal parameters. Hence, due to the properties of the Fourier transform, this HO procedure refines the positional and anisotropic thermal motion parameters only on core electrons. Because valence electron scattering factors diffuse only at low resolution, HO refinement gives precise positions and displacements of the core electrons only.In order to extract the aspherical features of the electron density, the following strategy was designed. First, a full-anisotropic spherical atom refinement using starting parameters from a Shelxl [6] refinement is performed. Then, the multipole refinement which has been developed by Hansen and Coppens [7] is selected to perform the electron density analysis using the MoPro package software [8].In parallel, the net atomic charge on the different atoms of the title compound and the molecular dipole moment were also estimated by theoretical calculations. Theoretical calculations were performed with the Gaussian package [9] at the HF/6-31G* level of theory and semiempirical methods.A higher-order refinement of X-ray data becomes essential to obtain accurate core positions and the associated thermal parameters of the different atoms of the molecule.The X-ray scattering factors for S, C, N, and O atoms come from the international tables for X-ray crystallography [10] while for H atom the bonded H-atom scattering curve of  Stewart et al. (1965) was used [11]. The H atoms were found by difference Fourier synthesis and their coordinates were adjusted by extending along the C–H and N–H bond directions to bond lengths of 1.105 Å and 1.020 Å, respectively. Any position error of the hydrogen atom will bring about errors strongly correlated with the dipolar population of hydrogen. Refinement details are displayed in Table 1.Table 1
Refinement details.
N
R
w
R
Spherical refinement
13
0.0274
0.03997
Multipole refinement
68
0.02358
0.03165
Wavelength,λ (Å)
0.71073
Completeness (%)
99.4
Resolution, Sinθ/λ (Å−1)
0.8
N is the number of refined parameters. R=∑|Fo|-|Fc|/∑|Fo|; wR=[∑w(|Fo|-|Fc|)2/∑w|Fo|2]1/2.


## 2.2. Charge Density Multipole Model

The Hansen-Coppens multipole formalism, as implemented in the Mopro least-squares program for multipole refinement, was used for both observed and theoretical structure factor fitting. It describes the crystal electron density as a superposition of aspherical pseudo-atoms each modeled on a multipole expansion:(1)ρatom(r→)=ρc(r→)+PVK′3ρv(k′r→)+∑llmax∑m=-1+1K′′Rl(K′′r→)PlmYlm(r→r).Hereρc and ρv are spherically averaged Hartree-Fock core and valence densities, with ρv being normalized to one electron. The Slater-type radial functions Rl(r→)=Nl·rlnlexp(-k′ξlr→) modulate the spherical harmonic angular functions Ylm± and N1 is a normalization factor. The values for parameters n=n1 and ξ were chosen according to rules provided by Coppens (1997) [12].In practice, two charge density variables, the population parameter and a kappa parameter which allows expansion and contraction of the valence shell, are added to the conventional parameters of structure analysis [13]. The population parameters Pv and Plm are floated along with κ′ and κ′′ during the refinement.The sum overm in (1) includes ±l, so that for each one, 2l +1 functions are included. The aspheric charge density was described at the octupole level (l=3) for the atoms C, N, and O and the hexadecapole for the atoms S (l=4) and at the dipole levels (l=1 or 2) for hydrogen atoms not involved and involved in strong H bonds, respectively. Charge densities of all hydrogen were considered to have cylindrical symmetry along the corresponding hydrogen-heavy atom bond.The quality of a refined model can be monitored based on the residuals and the goodness of fit apart from closely inspecting the deformation density maps. This model provides the advantage of using constraints based on assumed local symmetry on the atomic sites to limit the number of least-squares variables. Figure1 gives the local site symmetry of the different atoms of the molecule.Figure 1
Labeling of the atoms and definition of local orthogonal reference axes for the atom-centered multipolar functions.


## 3. Results and Discussion

### 3.1. Structural and Thermal Vibration Analysis

The main characteristic of this structure is that the eight molecules in the unit cell are engaged in seven chains formed by infinite chains of hydrogen-bonded coplanar molecules. The different chains link translationally equivalent molecules along thec crystallographic axis through C–S and N–H bonds.The rigid bond test indicated that the differences between the mean square displacements amplitudes (MSDAs) along interatomic directions have magnitudesΔ≤10×10-4Å2 for nearly all the bonded pairs of the non-H atoms [14]. In contrast, along intramolecular interatomic directions some Δ values were considerably larger, indicating that the molecule is not entirely rigid and that low-frequency soft modes of internal molecular vibration contribute significant1y to the observed mean square displacements [15]. The thermal motion analysis of the title compound has been performed using the THMA11 program [16].  Selected results are presented in Table 2.Table 2
Matrix for differences in MSDAs (mean square displacements of atoms)  (values listed are 104 MSDAs for column atom minus that for row atom; values correspond to chemical bonds).
Molecule A
C5
C1
C8
C2
C7
C4
S1
36
C3
−21
−19
N1
−4
C6
18
C4
7
C7
1
C2
16
Molecule B
C14
C13
C9
C15
C12
S2
42
N2
31
C11
−3
C12
7
C10
28
35
C15
1The rigid-body motion is described by three tensorsT, L, and S taking into account translation, liberation, and the correlation between translation and liberation, of the rigid group, respectively. These tensors are obtained by a least-squares fit refinement using the observed atomic thermal motion parameters obtained by the refinement.In terms of rigid-bodyTLS, the calculated anisotropic thermal parameters are given in the Trueblood notation as such:
(2)Uij=Tij+GijklLkl+HijklSkl+D2Ω2ninj,whereG and H are geometrical parameters. The analysis of the non-H atoms in terms of rigid-body motion yields the T, L, and S tensor values summarized in Table 3. The ellipsoids of the different atoms representing their thermal motion described before are shown with an ORTEPIII diagram [17] in Figure 2.Table 3
Rigid-body motion parameters.
L (rad2)
T (Å2)
S (rad Å)
(
0.00415
0.00036
-
0.00003
0.00011
-
0.00009
0.00029
)
(
0.02595
0.00051
0.00165
0.05841
-
0.00143
0.03275
)
(
-
0.00001
-
0.00055
-
0.00002
-
0.00011
-
0.00010
0.00050
-
0.00016
-
0.00011
0.00009
)Figure 2
ORTEPIII diagram of the title compound. Thermal ellipsoids are shown for the 30% probability level.


### 3.2. Electron Density Maps

The aspherical atom model used in multipole refinement gives structure factor phases closer to the true phases for crystals than the spherical or independent atom model does. This enables the mapping of the electron density by the Fourier synthesis in various ways using the WinGX program [18] to represent the experimental density deformation and the VMoPro program, implemented in the MoPro package, to observe the dynamic density deformation of the molecule.The residual map using a spherical atom model is shown in Figure3, from which we can notice the absence of the density on the atomic sites and the appearance of all the bond density peaks. This map confirms the high quality of the data sets and the efficiency of the formalism used for the data processing as proposed by Blessing [19].Figure 3
Experimental density map from high-order refinement. Contour map is 0.06 e·Å-3.The final electron density deformation was calculated using the following equation:(3)Δρdyn(r→)=1V∑H[|Fmul(H→)|eiφmul(H→)-|Fsph(H→)|eiφsph(H→)]e-i2πH→·r→.This map is obtained from the calculated multipole structure factorsFmul(H), where Fsph(H) is computed with atomic positions and thermal parameters obtained from the multipole refinement. Figure 4 shows this electron density deformation on the different atoms in the benzene ring section of the title compound, where one can easily observe the obvious increase of the density and the good localization peaks in the middle of covalent bonds of the cycle of the molecule. In Figure 5, the density map contains two independent molecules A and B. The density maps are given in the two representations of the methoxy groups. Figure 6(a) shows the density deformation at the plane yz and Figure 6(b) the deformation in the plane xy.Dynamic density maps of the two molecules A and B with the methoxy groups oriented in opposite conformations of the title compound. Contour map is 0.05 e·Å-3.
(a)
(b)Figure 5
Dynamic density map in the plane (S1, C1, and N1) contains two independent molecules A and B. Contour map is 0.05 e·Å-3.Dynamic density map of the methoxy groups: (a) deformation in the planexy of the molecule A; (b) deformation in the plane xy of the molecule B. Contour map is 0.05 e·Å-3.
(a)
(b)A residual density map in the molecular plane obtained in the final cycle of multipole refinement (see Figure7) shows the adequacy of the multipolar model to describe the electron experimental density of the molecule. The quasi-totality absence of the density peaks again confirms the high quality of the recorded data and the precision of the used equipment.Figure 7
Residual density map:ρres=ρ0-ρmult. Contour map is 0.05 e·Å-3.


### 3.3. Hydrogen Bond

Hydrogen bonds can be classified on the basis of charge density. The crystal structure of the title compound rests on chains of molecule joined by hydrogen bonds. The almost linear hydrogen bonding links by translation of equivalent molecules along thec crystallographic axis in the unit cell. In this crystal, the molecules form dimers via intermolecular N–H⋯S intermolecular hydrogen bonds. The dimers are further linked by C–H⋯O hydrogen bonds into chains along the c-axis.Figure8 shows the electron density in the plane of the hydrogen bonds formed by the atoms (N2, H2B, and S1), existing in the molecule. The electron density in the intermolecular hydrogen bonds is deficient in charge between the donor and acceptor of the molecule.Figure 8
Deformation density map in the plane of the hydrogen bond formed by the atoms (N2, H2B, and S1). Contour map is 0.05 e·Å-3.


### 3.4. Net Atomic Charges

Knowledge of the charge density function can lead to some important properties of the molecule such as the net charges on the different atoms, the molecular dipole moment, and the electrostatic potential around the molecule. The valence population coefficientsPvi were used to estimate the partial charges on the different atoms according to the following equation:
(4)qi=ni-Pvi,whereni is the total number of electrons of atoms i. The population coefficients are presented in Table 4. The refined atomic charges are generally of the sign expected from chemical knowledge: H atoms are positively charged, whereas O, S, and N atoms are negatively charged. The net atomic charges obtained by the X-ray experimental method using the multipolar model are compared to those derived from the theoretical calculations. Semiempirical calculations were carried out using the AM1 (Austin’s method) and ab initio calculations were carried out at the HF/6-31G* levels of theory using the Gaussian package. The results obtained by both theory and experiment methods are summarized in Table 4. The experimental results and those derived from the semiempirical and ab initio calculations are in good agreement concerning the signs of the net charges of the H atoms and those of the N, S, and O atoms.Table 4
Net charges (q) in the different atoms of the title compound.
Atoms
q
Experimental
Ab initio
Semiempirical
S1
−0.332
−0.250
−0.223
O1
−0.124
−0.243
−0.207
N1
−0.004
−0.026
−0.035
C6
0.026
−0.038
−0.166
C3
−0.020
−0.079
−0.077
C5
0.022
0.025
−0.146
C2
0.084
−0.098
0.102
C1
−0.002
−0.026
−0.212
C8
0.138
−0.037
−0.078
C7
−0.019
−0.053
0.077
C4
0.026
0.134
−0.036
H1A
0.017
0.207
0.221
H1B
0.017
0.213
0.243
H3
0.019
0.070
0.135
H4
0.047
0.060
0.133
H6
0.053
0.076
0.163
H7
0.019
0.063
0.144
H8A
0.104
0.075
0.109
H8B
0.104
0.066
0.077
H8C
0.104
0.062
0.075


### 3.5. Molecular Dipole Moment

The molecular dipole moment was first estimated using the population coefficients of the multipolar density model from the following equation:(5)μ=∑iqiri+∑i4n+33ki′′(Pxii+Pyij+Pzik),whereri is the atomic position and n is the number of electrons in the molecule. Pi are the dipolar electron populations. Dipole moment orientation in the molecule is shown in Figure 9. The resulting molecular dipole moment direction is in accordance with the evaluation of the positive sign of the net charges on the H atoms and the negative sign of the net charges on the O, S, and N atoms.Figure 9
Orientation of the molecular dipole moment of the title compound; the origin is at the center of mass of the molecule;μM: molecular dipole moment from the multipolar model; μS: molecular dipole moment from the semiempirical calculation; μA: molecular dipole moment from the ab initio calculation.In a parallel study, the net charges on the different atoms and the molecular dipole moment were also estimated by ab initio (B3LYP/6-31G*) and semiempirical (AM1) calculations. In this procedure the treatment is confined to the valence-shell electrons. The results of these calculations are presented in Table5.Table 5
Components of the dipolar moment of the charge distribution from the point charge models (AM1 and B3LYP/6-31G*) and comparison with the components derived from experimental electron density (multipolar refinement).
Methods
Models
μ
x
μ
y
μ
z
μ (Debye)
X-Ray experiment
Multipolar refinement
−2.1064
1.5995
1.5281
3.0546
Semiempirical
AM1
−1.2175
1.3216
1.3220
2.2308
Ab initio
B3LYP/6-31G*
−2.4646
2.2520
1.8834
3.8332From the results, we can see a good agreement between the value obtained with X-ray diffraction (3.0546 D) and the other values obtained by theoretical calculation (AM1 gives 2.2308 D and DFT gives 3.8332 D).Analysis of obtained dipole moments shows that the size and direction of the dipole moment vector depend mainly on the position of substituents. Furthermore, the orientation of the dipole moment vector is toward the sulfur atom.


### 3.6. Electrostatic Potential

The distribution of positive and negative charges in a crystal fully defines physical properties like the electrostatic potential and its derivatives. The electrostatic potential is important in the study of intermolecular interactions of molecular systems.For a continuous charge distribution study, the potential is obtained by integration over the space containing the distribution. At a point defined byr′, the potential in given by
(6)ϕ(r′)=∫ρtotal(r)|r-r′|dr.In whichρtotal represents both the nuclear and the electronic charge, where the integration is over the molecular volume, and r represents the atomic position relative to common origin.The integration includes the atoms of only one molecule and therefore does not include directly the effects of charge distribution of the neighboring molecules.Figure10 shows the electrostatic potential distribution in the plane of the base ring. The extension of the positive electrostatic potential around the carbothioamide and C–H groups and the regions of negative electrostatic potential around the nitro and hydroxyl groups confirm the nature of the intramolecular charge transfer as found by the orientation of the molecular dipole moment.Figure 10
The electrostatic potential maps around the molecule. The section is in the plane of the benzene ring.The obtained point charges corresponding to S atom are −0.332, −0.250, and −0.223 e from experimental, ab initio, and semiempirical methods, respectively. But the region very close to S appears positive. This is because of the fact that S atom of the title compound is surrounded by the most electropositive hydrogen atoms (H1A = 0.017, H1B = 0.017, H3 = 0.019, and H7 = 0.019 e).


## 3.1. Structural and Thermal Vibration Analysis

The main characteristic of this structure is that the eight molecules in the unit cell are engaged in seven chains formed by infinite chains of hydrogen-bonded coplanar molecules. The different chains link translationally equivalent molecules along thec crystallographic axis through C–S and N–H bonds.The rigid bond test indicated that the differences between the mean square displacements amplitudes (MSDAs) along interatomic directions have magnitudesΔ≤10×10-4Å2 for nearly all the bonded pairs of the non-H atoms [14]. In contrast, along intramolecular interatomic directions some Δ values were considerably larger, indicating that the molecule is not entirely rigid and that low-frequency soft modes of internal molecular vibration contribute significant1y to the observed mean square displacements [15]. The thermal motion analysis of the title compound has been performed using the THMA11 program [16].  Selected results are presented in Table 2.Table 2
Matrix for differences in MSDAs (mean square displacements of atoms)  (values listed are 104 MSDAs for column atom minus that for row atom; values correspond to chemical bonds).
Molecule A
C5
C1
C8
C2
C7
C4
S1
36
C3
−21
−19
N1
−4
C6
18
C4
7
C7
1
C2
16
Molecule B
C14
C13
C9
C15
C12
S2
42
N2
31
C11
−3
C12
7
C10
28
35
C15
1The rigid-body motion is described by three tensorsT, L, and S taking into account translation, liberation, and the correlation between translation and liberation, of the rigid group, respectively. These tensors are obtained by a least-squares fit refinement using the observed atomic thermal motion parameters obtained by the refinement.In terms of rigid-bodyTLS, the calculated anisotropic thermal parameters are given in the Trueblood notation as such:
(2)Uij=Tij+GijklLkl+HijklSkl+D2Ω2ninj,whereG and H are geometrical parameters. The analysis of the non-H atoms in terms of rigid-body motion yields the T, L, and S tensor values summarized in Table 3. The ellipsoids of the different atoms representing their thermal motion described before are shown with an ORTEPIII diagram [17] in Figure 2.Table 3
Rigid-body motion parameters.
L (rad2)
T (Å2)
S (rad Å)
(
0.00415
0.00036
-
0.00003
0.00011
-
0.00009
0.00029
)
(
0.02595
0.00051
0.00165
0.05841
-
0.00143
0.03275
)
(
-
0.00001
-
0.00055
-
0.00002
-
0.00011
-
0.00010
0.00050
-
0.00016
-
0.00011
0.00009
)Figure 2
ORTEPIII diagram of the title compound. Thermal ellipsoids are shown for the 30% probability level.


## 3.2. Electron Density Maps

The aspherical atom model used in multipole refinement gives structure factor phases closer to the true phases for crystals than the spherical or independent atom model does. This enables the mapping of the electron density by the Fourier synthesis in various ways using the WinGX program [18] to represent the experimental density deformation and the VMoPro program, implemented in the MoPro package, to observe the dynamic density deformation of the molecule.The residual map using a spherical atom model is shown in Figure3, from which we can notice the absence of the density on the atomic sites and the appearance of all the bond density peaks. This map confirms the high quality of the data sets and the efficiency of the formalism used for the data processing as proposed by Blessing [19].Figure 3
Experimental density map from high-order refinement. Contour map is 0.06 e·Å-3.The final electron density deformation was calculated using the following equation:(3)Δρdyn(r→)=1V∑H[|Fmul(H→)|eiφmul(H→)-|Fsph(H→)|eiφsph(H→)]e-i2πH→·r→.This map is obtained from the calculated multipole structure factorsFmul(H), where Fsph(H) is computed with atomic positions and thermal parameters obtained from the multipole refinement. Figure 4 shows this electron density deformation on the different atoms in the benzene ring section of the title compound, where one can easily observe the obvious increase of the density and the good localization peaks in the middle of covalent bonds of the cycle of the molecule. In Figure 5, the density map contains two independent molecules A and B. The density maps are given in the two representations of the methoxy groups. Figure 6(a) shows the density deformation at the plane yz and Figure 6(b) the deformation in the plane xy.Dynamic density maps of the two molecules A and B with the methoxy groups oriented in opposite conformations of the title compound. Contour map is 0.05 e·Å-3.
(a)
(b)Figure 5
Dynamic density map in the plane (S1, C1, and N1) contains two independent molecules A and B. Contour map is 0.05 e·Å-3.Dynamic density map of the methoxy groups: (a) deformation in the planexy of the molecule A; (b) deformation in the plane xy of the molecule B. Contour map is 0.05 e·Å-3.
(a)
(b)A residual density map in the molecular plane obtained in the final cycle of multipole refinement (see Figure7) shows the adequacy of the multipolar model to describe the electron experimental density of the molecule. The quasi-totality absence of the density peaks again confirms the high quality of the recorded data and the precision of the used equipment.Figure 7
Residual density map:ρres=ρ0-ρmult. Contour map is 0.05 e·Å-3.


## 3.3. Hydrogen Bond

Hydrogen bonds can be classified on the basis of charge density. The crystal structure of the title compound rests on chains of molecule joined by hydrogen bonds. The almost linear hydrogen bonding links by translation of equivalent molecules along thec crystallographic axis in the unit cell. In this crystal, the molecules form dimers via intermolecular N–H⋯S intermolecular hydrogen bonds. The dimers are further linked by C–H⋯O hydrogen bonds into chains along the c-axis.Figure8 shows the electron density in the plane of the hydrogen bonds formed by the atoms (N2, H2B, and S1), existing in the molecule. The electron density in the intermolecular hydrogen bonds is deficient in charge between the donor and acceptor of the molecule.Figure 8
Deformation density map in the plane of the hydrogen bond formed by the atoms (N2, H2B, and S1). Contour map is 0.05 e·Å-3.


## 3.4. Net Atomic Charges

Knowledge of the charge density function can lead to some important properties of the molecule such as the net charges on the different atoms, the molecular dipole moment, and the electrostatic potential around the molecule. The valence population coefficientsPvi were used to estimate the partial charges on the different atoms according to the following equation:
(4)qi=ni-Pvi,whereni is the total number of electrons of atoms i. The population coefficients are presented in Table 4. The refined atomic charges are generally of the sign expected from chemical knowledge: H atoms are positively charged, whereas O, S, and N atoms are negatively charged. The net atomic charges obtained by the X-ray experimental method using the multipolar model are compared to those derived from the theoretical calculations. Semiempirical calculations were carried out using the AM1 (Austin’s method) and ab initio calculations were carried out at the HF/6-31G* levels of theory using the Gaussian package. The results obtained by both theory and experiment methods are summarized in Table 4. The experimental results and those derived from the semiempirical and ab initio calculations are in good agreement concerning the signs of the net charges of the H atoms and those of the N, S, and O atoms.Table 4
Net charges (q) in the different atoms of the title compound.
Atoms
q
Experimental
Ab initio
Semiempirical
S1
−0.332
−0.250
−0.223
O1
−0.124
−0.243
−0.207
N1
−0.004
−0.026
−0.035
C6
0.026
−0.038
−0.166
C3
−0.020
−0.079
−0.077
C5
0.022
0.025
−0.146
C2
0.084
−0.098
0.102
C1
−0.002
−0.026
−0.212
C8
0.138
−0.037
−0.078
C7
−0.019
−0.053
0.077
C4
0.026
0.134
−0.036
H1A
0.017
0.207
0.221
H1B
0.017
0.213
0.243
H3
0.019
0.070
0.135
H4
0.047
0.060
0.133
H6
0.053
0.076
0.163
H7
0.019
0.063
0.144
H8A
0.104
0.075
0.109
H8B
0.104
0.066
0.077
H8C
0.104
0.062
0.075


## 3.5. Molecular Dipole Moment

The molecular dipole moment was first estimated using the population coefficients of the multipolar density model from the following equation:(5)μ=∑iqiri+∑i4n+33ki′′(Pxii+Pyij+Pzik),whereri is the atomic position and n is the number of electrons in the molecule. Pi are the dipolar electron populations. Dipole moment orientation in the molecule is shown in Figure 9. The resulting molecular dipole moment direction is in accordance with the evaluation of the positive sign of the net charges on the H atoms and the negative sign of the net charges on the O, S, and N atoms.Figure 9
Orientation of the molecular dipole moment of the title compound; the origin is at the center of mass of the molecule;μM: molecular dipole moment from the multipolar model; μS: molecular dipole moment from the semiempirical calculation; μA: molecular dipole moment from the ab initio calculation.In a parallel study, the net charges on the different atoms and the molecular dipole moment were also estimated by ab initio (B3LYP/6-31G*) and semiempirical (AM1) calculations. In this procedure the treatment is confined to the valence-shell electrons. The results of these calculations are presented in Table5.Table 5
Components of the dipolar moment of the charge distribution from the point charge models (AM1 and B3LYP/6-31G*) and comparison with the components derived from experimental electron density (multipolar refinement).
Methods
Models
μ
x
μ
y
μ
z
μ (Debye)
X-Ray experiment
Multipolar refinement
−2.1064
1.5995
1.5281
3.0546
Semiempirical
AM1
−1.2175
1.3216
1.3220
2.2308
Ab initio
B3LYP/6-31G*
−2.4646
2.2520
1.8834
3.8332From the results, we can see a good agreement between the value obtained with X-ray diffraction (3.0546 D) and the other values obtained by theoretical calculation (AM1 gives 2.2308 D and DFT gives 3.8332 D).Analysis of obtained dipole moments shows that the size and direction of the dipole moment vector depend mainly on the position of substituents. Furthermore, the orientation of the dipole moment vector is toward the sulfur atom.


## 3.6. Electrostatic Potential

The distribution of positive and negative charges in a crystal fully defines physical properties like the electrostatic potential and its derivatives. The electrostatic potential is important in the study of intermolecular interactions of molecular systems.For a continuous charge distribution study, the potential is obtained by integration over the space containing the distribution. At a point defined byr′, the potential in given by
(6)ϕ(r′)=∫ρtotal(r)|r-r′|dr.In whichρtotal represents both the nuclear and the electronic charge, where the integration is over the molecular volume, and r represents the atomic position relative to common origin.The integration includes the atoms of only one molecule and therefore does not include directly the effects of charge distribution of the neighboring molecules.Figure10 shows the electrostatic potential distribution in the plane of the base ring. The extension of the positive electrostatic potential around the carbothioamide and C–H groups and the regions of negative electrostatic potential around the nitro and hydroxyl groups confirm the nature of the intramolecular charge transfer as found by the orientation of the molecular dipole moment.Figure 10
The electrostatic potential maps around the molecule. The section is in the plane of the benzene ring.The obtained point charges corresponding to S atom are −0.332, −0.250, and −0.223 e from experimental, ab initio, and semiempirical methods, respectively. But the region very close to S appears positive. This is because of the fact that S atom of the title compound is surrounded by the most electropositive hydrogen atoms (H1A = 0.017, H1B = 0.017, H3 = 0.019, and H7 = 0.019 e).


## 4. Conclusion

This study has obtained good accurate results on the structure which gives high-quality descriptive model for the electron charge density distribution from X-ray diffraction experiment. It also shows that electron density can yield electronic proprieties such as dipole moment. Its determination lets us locate the region of the donor and the acceptor groups and also the direction of charge transfer within the molecule. The deformation of the electrostatic potential of the title compound shows that the electronegative potential is located on the side of methoxy group and electropositive potential is on the side of carbothioamide group.In the multipole refinement, all nonhydrogen atoms were refined anisotropically and subject to a rigid-bond constraint. H atoms anisotropic parameters are obtained fromTLS analysis. The final R factor (2.2%) attests to the quality of the data.To study the effect of electron charge transfer in biological activity more, we should perform vibrational analysis using FT-Raman and FTIR spectroscopic techniques.



---

*Source: 326457-2013-08-14.xml*