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–HS intermolecular hydrogen bonds. The dimers are further linked by C–HO 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.

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 [13]. 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.

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:

Here and are spherically averaged Hartree-Fock core and valence densities, with being normalized to one electron. The Slater-type radial functions modulate the spherical harmonic angular functions and is a normalization factor. The values for parameters 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 and are floated along with and during the refinement.

The sum over in (1) includes , so that for each one, +1 functions are included. The aspheric charge density was described at the octupole level () for the atoms C, N, and O and the hexadecapole for the atoms S () and at the dipole levels ( 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. Figure 1 gives the local site symmetry of the different atoms of the molecule.

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 the 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 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.

The rigid-body motion is described by three tensors , , and 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-body , the calculated anisotropic thermal parameters are given in the Trueblood notation as such:

where and are geometrical parameters. The analysis of the non-H atoms in terms of rigid-body motion yields the , , and 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.

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 Figure 3, 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].

The final electron density deformation was calculated using the following equation:

This map is obtained from the calculated multipole structure factors (), where () 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 and Figure 6(b) the deformation in the plane .

A residual density map in the molecular plane obtained in the final cycle of multipole refinement (see Figure 7) 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.

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 the crystallographic axis in the unit cell. In this crystal, the molecules form dimers via intermolecular N–HS intermolecular hydrogen bonds. The dimers are further linked by C–HO hydrogen bonds into chains along the -axis.

Figure 8 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.

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 coefficients were used to estimate the partial charges on the different atoms according to the following equation:

where is the total number of electrons of atoms . 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.

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:

where is the atomic position and is the number of electrons in the molecule. 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.

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 Table 5.

From 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 by , the potential in given by

In which represents both the nuclear and the electronic charge, where the integration is over the molecular volume, and 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.

Figure 10 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.

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 from analysis. The final 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.

Supplementary Materials

The multipolar coefficients (molecule A) described as the model of Hansen & Coppens.

  1. Supplementary Table