Abstract

Optimized geometrical structures of alkali metal phenoxyacetates were obtained using B3LYP/6-311++G** method. Geometric and magnetic aromaticity indices, dipole moments, and energies were calculated. Atomic charges on the atoms of phenoxyacetic acid molecule and its alkali metal salts were calculated by Mulliken, APT (atomic polar tensor), NPA (natural population analysis), MK (Merz-Singh-Kollman method), and ChelpG (charges from electrostatic potentials using grid-based method) methods. The theoretical wavenumbers and intensities of IR as well as chemical shifts in NMR spectra were obtained and compared with experimental data. The effect of alkali metals on molecular structure of phenoxyacetic acid appears in the shift of selected bands along the series of alkali metal salts. The correlations between chosen bands and some metal parameters, such as electronegativity, ionization energy, and atomic, and ionic radius, have been noticed.

1. Introduction

Antimicrobial activity, anticancer, analgesic and anti-inflammatory are some of the important properties of phenoxyacetic acids, which are widely used as plant growth regulators and herbicides. Due to their solubility in water they can move in agricultural ecosystems and cause the pollution of surface and ground waters [1]. It seems to be interesting to study the formation of phenoxyacetate salts with common ions in soil. This paper is a continuation of our reports on the metal effect on the electronic system of carboxylic acids [2, 3]. The influence of lithium, sodium, potassium, rubidium, and cesium on the electronic system of the salt molecules was analyzed in this paper.

2. Experimental

Lithium, sodium, potassium, rubidium, and cesium phenoxyacetates were prepared by dissolving the powder of phenoxyacetic acid in the water solution of the appropriate metal hydroxide in a stoichiometric ratio (1 : 1). All reagents were obtained from Aldrich Chemical Company. The solution was left at the room temperature for 24 h until the sample crystallized in the solid state. Precipitants were filtered, washed by water, and dried under reduced pressure at 110°C. Obtained complexes were anhydrous—in the IR spectra of solid state samples the lack of bands characterized for crystallizing water was observed.

The IR spectra were recorded with the Equinox 55, Bruker FT-IR spectrometer within the range 4000–400 cm−1. Samples in the solid state were measured in KBr pellets. The resolution of spectrometer was 1 cm−1. UV/Vis spectra were recorded with the DR/4000 U, HACH spectrophotometer. The NMR spectra of DMSO solution were recorded with the NMR AC 200 F Bruker unit. TMS was used as an internal reference.

The density functional (DFT) hybrid method B3LYP/6-311++G** was used to calculate optimized geometrical structures of studied compounds (Figure 1). All theoretical calculations were performed using the GAUSSIAN 09 [4] program running on a PC computer.

3. Results and Discussion

3.1. Calculated Data

Optimized geometrical structures of lithium (LiPA), sodium (NaPA), and potassium (KPA) phenoxyacetates were calculated using B3LYP/6-311++G** method. The bond lengths and the angles between bonds in studied molecules are presented in Table 1. The changes of almost all bond lengths in aromatic ring are noticed in the order from lithium to sodium phenoxyacetate molecules. The significant increasement of C1–C2 and C6–C1 bond lengths is observed in the series LiPA < NaPA < KPA (atomic numbers are shown in Figure 1). The bond lengths of C3–C4 and C4–C5 insignificantly increase whereas C2–C3 and C5–C6 insignificantly decrease in the above order. The significant decrease of C1–O7, C9–O10, and C9–O11 bond lengths is also observed in the following order: LiPA > NaPA > KPA. But the significant increase is noticed for O7–C8, C8–C9, O10–M, and O11–M bond lengths in the same series. As they might have been expected, the differences between O10–M and O11–M significantly decrease. In the case of angles, the increase in the order LiPA < NaPA < KPA was observed for C1–C2–C3, C2–C3–C4, and C5–C6–O1 angles in the aromatic ring, but the decrease in this order was noticed for C3–C4–C5 and C6–C1–C2. The angles between carbon, oxygen, and metal atoms in alkil chain of molecules increase with exception of C8–C9–O10 and C8–C9–O11 angles.

Geometric and magnetic aromaticity indices [57], dipole moments, and energies were calculated and also shown in Table 1. Almost all aromaticity indices decrease in the following order: LiPA > NaPA > KPA. Values of energy also decrease in the same series. It might indicate the increase of stability of salt molecules in this order. Mulliken, APT, NPA, MK, and ChelpG atomic charges on the atoms of phenoxyacetic acid molecule and its alkali metal salts are gathered in Table 2. The increase of total charge of carboxylate group (negative values) are observed in the order LiPA < NaPA < KPA when Mulliken and APT methods were used for calculation. In NPA method there are no changes whereas the decrease is observed for MK and ChelpG methods. The total charge in aromatic ring is not changed in above series irrespective of used method.

The theoretical wavenumbers of IR and Raman spectra as well as chemical shifts in NMR spectra were obtained and compared with experimental spectra.

3.2. FT-IR, FT-Raman, and UV Spectra

The vibrational (in KBr) and electronic spectra of synthesized alkali metal phenoxyacetates were recorded and presented in Figure 2. In the UV spectra there were no changes in the wavelengths. The influence of studied metal cations on the vibrational structure of phenoxyacetates appearsin the shift of selected bands along the metal series. The wavenumbers of asymmetric stretching vibration bands of carboxylate group in IR spectra decrease in the series LiPA→CsPA; however, for rubidium phenoxyacetate the insignificant increase is noted in comparison to potassium salt. In Raman spectra there are no changes in this series. For symmetric stretching vibration bands in Raman spectra the increase of the wavenumbers is noticed from lithium to cesium phenoxyacetate, but in IR spectra the wavenumbers increase to potassium phenoxyacetate and from KPA to CsPA the decrease is observed. The magnitudes of separation between wavenumbers due to asymmetric and symmetric stretching vibration bands Δ 𝜈 a s s = 𝜈 a s ( C O O ) 𝜈 s ( C O O ) in Raman spectra decrease in the series LiPA > NaPA = KPA > RbPA > CsPA. The corresponding values for lithium, sodium, potassium, rubidium, and cesium phenoxyacetates in IR spectra are equal to 382, 361, 353, 357, and 361, respectively. The wavenumbers of asymmetric in-plane deformations 𝛽 a s (COO) decreased in the following series LiPA > NaPA > KPA > RbPA > CsPA, whereas the magnitudes of separation between wavenumbers due to asymmetric and symmetric in plane deformations of C O O group increase in this order. The values of Δ 𝛽 = 𝛽 s ( C O O ) 𝛽 a s ( C O O ) are equal to 518, 524, 528, 530, and 536 cm−1. The wavenumbers as well as intensity of the aromatic ring bands did not significantly change.

The correlations between chosen bands and some metal parameters, such as electronegativity, ionization energy, atomic and ionic radius, atomic mass, and inverse of atomic mass [8] have been studied. The best correlation was obtained for ionic radius and ionization energy for band of 𝜈 (C–O–C) etheric group vibration ( 𝑅 2 = 0 . 9 7 8 ).

3.3. NMR Spectra

Values of the experimental as well as calculated (B3LYP/6-311++G**) chemical shifts in 1H and 13C NMR spectra of lithium, sodium, potassium, rubidium, and cesium phenoxyacetates are presented in Table 3. Data for the last two salts are only experimentally obtained. Experimental chemical shifts show decreasing tendency in the series Li > Na > K > Rb > Cs salts for H8a and H8b protons, but for protons connected with the aromatic ring the decrease was observed only for lithium, sodium, and potassium phenoxyacetates, and for rubidium and cesium salts these values increase in comparison to potassium one. The changes in 13C NMR spectra indicated decreasing tendency from lithium to potassium phenoxyacetate, and then the increase from potassium to rubidium, and finally the insignificant decrease for cesium salt were observed for almost all carbon atoms in the aromatic ring (except of C1 and C4) and C9 atom. For C1 and C8 carbon atoms the opposite tendency was observed, whereas an insignificant decreasing tendency was observed for C4 carbon atom in the series LiPA = NaPA > KPA > RbPA > CsPA.

4. Conclusions

Some regular changes were observed in the calculated values of studied alkali metal phenoxyacetates. Almost all bond lengths decreased in the series LiPA > NaPA > KPA. All geometric aromaticity indices decreased in the same order. Values of energy also decreased in the studied series. On the other hand the shift of selected bands in IR, Raman, and NMR spectra along the phenoxyacetate series was also observed.

The calculated parameters were compared to the experimental characteristics of studied compounds. Good correlation between experimental and theoretical IR and Raman spectra was noted. The correlation coefficients ( 𝑅 2 ) for IR spectra for lithium, sodium, and potassium phenoxyacetates amount to 0.9986, 0.9990, and 0.9984, respectively. The corresponding values for Raman spectra are 0.9970, 0.9975, and 0.9973. The linear correlation between calculated and experimental data of NMR spectra was also observed. Correlation coefficients for 13C NMR are in the range 0.9002–0.9063 and for 1H NMR the range is 0.9845–0.9878.

It was interesting to note the linear correlation between wavenumbers of 𝜈 s (COO) and 𝜈 a s (COO) bands and distances between C9–O11 and O10–M. The corresponding correlation coefficients for 𝜈 s (COO) wavenumbers amount to 0.9951 and 0.9767 and for 𝜈 a s (COO) 0.9999 and 0.9469. Good correlation between chemical shifts in 13C NMR spectra and Mulliken atomic charges on C1 and C6 atoms was also found ( 𝑅 2 = 0 . 9 7 1 8 and 0.9932, resp.).

Acknowledgment

The presented work was supported by Białystok University of Technology (Theme no. S/WBiIŚ/1/2012).