Journal of Chemistry
Volume 2017 (2017), Article ID 6537402, 9 pages
https://doi.org/10.1155/2017/6537402

Research Article

Synthesis, Crystal Structure, and DFT Study of Ethyl 1-(2-(Hydroxyimino)-2-phenylethyl)-3-phenyl-1H-pyrazole-5-carboxylate

Jin-Hui Zhou,1,2 Liang-Wen Zheng,² Mao-Cai Yan,¹ Mao-Jian Shi,¹ Jing Liu,¹ and Guo-Qiang Shangguan¹

¹College of Pharmacy, Jining Medical University, Rizhao 276826, China
²School of Chemistry and Chemical Engineering, University of Shandong, Jinan 250100, China

Correspondence should be addressed to Guo-Qiang Shangguan; gqshangguan@mail.jnmc.edu.cn

Received 6 July 2017; Revised 7 August 2017; Accepted 30 August 2017; Published 18 October 2017

Academic Editor: Josefina Pons

Copyright © 2017 Jin-Hui Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The crystal structure of ethyl 1-(2-(hydroxyimino)-2-phenylethyl)-3-phenyl-1H-pyrazole-5-carboxylate has been determined by X-ray single crystal diffraction. The crystal of the title compound is the monoclinic space group P2/c with unit cell parameters of  Å,  Å,  Å, ,  Å³, and . The dihedral angles formed by the planes of the central pyrazole ring and the adjacent benzene rings are 73.60(7)° and 3.55(7)°, respectively. The combination of the weak intermolecular C-HO and N-HO hydrogen-bonding interactions stabilizes the crystal packing. The geometries of its Z and E isomers and the corresponding transition state (TS), as well as the dimer of its Z isomer, are optimized using the B3LYP hybrid functional coupled with def-TZVP triple-zeta polarized basis set. The bond angles and bond lengths of the optimized structure of Z dimer are very consistent with those of its single crystal parameters. Double-hybrid functional PWPB95-D3 in combination with very highly accurate basis set def2-QZVP is employed to evaluate accurate energy of each isomer and TS. The calculated equilibrium constant between Z and E isomers corresponds to the [Z]/[E] ratio of 4.29. Mulliken atomic charges and electrostatic potential (ESP) on molecular van der Waals (vdW) surface are calculated in order to study and predict the intermolecular interactions. The molecular total energies and frontier orbital analysis are also discussed.

1. Introduction

Pyrazole oximes are one class of heterocyclic compounds which play a key role in the field of chemistry and pharmacology. It has been found that they exhibit various biological activities that are often related to low toxicity and increase a great amount of beneficial influences on human health. Therefore, they have attracted wide interest for decades owing to their diverse biological properties such as insecticidal [1], acaricidal [2], antitumor [3–5], fungicidal [6], antiviral [7, 8], antibacterial [9–11], COX inhibitory [12], and fluorescent [13, 14] properties.

In our paper, the structure of ethyl 1-(2-(hydroxyimino)-2-phenylethyl)-3-phenyl-1H-pyrazole-5-carboxylate is determined by single crystal X-ray diffraction [15, 16]. Its configuration of oxime bond and intermolecular interactions is confirmed. We have conducted density functional theory (DFT) calculations to prove that selected hybrid functionals and basis sets are particularly suitable for the calculations of the structural parameters, intermolecular hydrogen bonds, and equilibrium constant between Z and E isomers for oximes. Besides, this work provides a possibility to predict these parameters when a certain compound lacks its single crystal. Based on the reliable structure model, the DFT calculations are applied to investigate the optimized structural parameters, which is contrasted with the experimental X-ray diffraction structure. The geometries of Z and E isomers and the corresponding transition state (TS), as well as the dimer of its Z isomer, are optimized using B3LYP hybrid functional coupled with def-TZVP triple-zeta polarized basis set. The bond angles and bond lengths of optimized structure of Z dimer are very consistent with those of its single crystal parameters. Double-hybrid functional PWPB95-D3 in combination with very highly accurate basis set def2-QZVP is employed to evaluate the accurate energy of each isomer and TS. The calculated equilibrium constant between Z and E isomers is calculated to be 4.29 ([Z]/[E]), which is close to the value estimated according to the experimental NMR spectra. The single crystal structure is Z configuration. Mulliken atomic charges and electrostatic potential (ESP) on molecular van der Waals (vdW) surface are calculated to study and predict intermolecular interactions. It is expected that the intermolecular hydrogen bonds must exist in title compound’s molecular crystal. Indeed, the actual crystal structure nicely supports this point. In addition, the frontier orbital analysis and molecular total energies are discussed.

2. Experimental

2.1. Materials and Instrumentation

Chemical shifts of ¹³C NMR and ¹H NMR spectra are recorded in ppm and NMR spectra are reported on a Bruker AVANCE 400 MHz spectrometer (Bruker, Billerica, MA, USA) with the application of solvent shifts as internal criterion for ¹³C and ¹H (CDCl₃: ¹³C = 77.0, ¹H = 7.26). MS spectra are determined by LTQ Orbitrap Hybrid Mass Spectrometer (Thermo Fisher Scientific, Waltham, MA, USA). X-ray single crystal experiments are performed on a Bruker SMART CCD single crystal diffractometer (Bruker, Billerica, MA, USA).

2.2. Preparation of Title Compound

Synthesis of this compound (Scheme 1) refers to a previous literature [3]. Ethyl 1-(2-oxo-2-phenylethyl)-3-phenyl-1H-pyrazole-5-carboxylate (349.38 mg, 1.0 mmol) and hydroxylamine hydrochloride (347.45 mg, 5.0 mmol) in EtOH (40 ml) are added to a 100 ml flask with a magnetic stirrer. The flask is stirred and heated slowly. After reflux for about 1 h, sodium acetate anhydrous (410.15 mg, 5.0 mmol) is added to the flask, which maintains continuous reflux for 1.5–4 h. The solvent is removed; the resulting residue is partitioned with water and ethyl acetate layers afterwards. Brine and water are used to wash the organic layer successively, and then it is dried over MgSO₄, evaporating under reduced pressure to give a residue. Recrystallization from absolute ethanol or silica gel column chromatography is used to purify the residue. The obtained target compound is a white solid. Single crystal is collected from the solution, which is obtained after the evaporation of solvent at low temperature for days.

Scheme 1: Synthesis of ethyl 1-(2-(hydroxyimino)-2-phenylethyl)-3-phenyl-1H-pyrazole-5-carboxylate.

2.3. Crystal Structure Determination

Crystals suitable for diffraction experiment are collected from the solutions. In addition to the data collection, Table 1 also summarizes the structural refinement details. The structure is solved and then refined using the OLEX2 program suite [17] that is equipped with SHELXS and SHELXL program [18]. All nonhydrogen atoms are placed from Fourier map directly by SHELXS and refined anisotropically. Hydrogen atoms on all nonhydrogen atoms are put in calculated positions, and their displacement parameters and coordinates are forced to depend on the carrier atoms. The CIF file could be downloaded from http://www.ccdc.cam.ac.uk.

Table 1: Crystal data of title compound.

CCDC number
1060837

Empirical formula
C₂₀H₁₉N₃O₃

Formula weight
349.38

Crystal system, space group
Monoclinic, P2/c

Temperature (K)
293.15

, , (Å)
8.634(2), 9.616(2), 22.190(3)

(Å³)
1818.3(4)

(°)
99.265(2)

4

Radiation type
MoKa

(mm⁻¹)
0.088

Crystal size (mm)
0.14 × 0.19 × 0.23

Data collection

Diffractometer
Bruker SMART CCD area-detector diffractometer

Absorption correction
Multiscan (SADABS in SAINT; Bruker, 2002)

min, max
0.97, 0.99

Number of measured, independent, and observed [I > 2σ(I)] reflections
15059, 4196, 2690

0.0307

Refinement [F² > 2σ(F²)], (F²), S
0.0573, 0.1751

Number of reflections
4196

Number of parameters
239

H-Atom treatment
Mixed

, (Å⁻³)
0.71, −0.26

2.4. Computational Details

The geometries of Z and E isomers as well as the corresponding transition state (TS) are optimized using B3LYP hybrid functional [19] coupled with def-TZVP triple-zeta polarized basis set [20] within Gaussian 03 package [21]. Then the optimized structures are subjected to frequency analysis at the same level to obtain thermodynamic correction values. Vibrational frequency analysis shows that the optimized structures correspond to the lowest points on the potential energy surfaces without virtual frequencies, which means that the optimized structure is stable. Double-hybrid functional PWPB95-D3 [22] in combination with very highly accurate basis set def2-QZVP [23] is employed to evaluate the accurate energy of each isomer and TS. In this process, the resolution of identity (RI) technique [24] is applied to accelerate electron integral calculations. The molecular electrostatic potential at the B3LYP/TZVP level is analyzed by Multiwfn 3.3.7 software [25] and then plotted using VMD 1.9.1 program [26]. The composition (%) of atomic orbitals in molecular frontier orbitals is analyzed by Multiwfn software.

3. Results and Discussion

3.1. Crystal Structure

As shown in Figure 1, the asymmetrical unit contains an independent molecule of title compound. The detailed analysis of the molecule shows that all the bond lengths and bond angles are within the normal ranges. The molecule is nonplanar. The dihedral angles that are made by the plane of the central pyrazole ring (ring-2) and its adjacent benzene rings (ring-1 and ring-3) are 73.60(7)° and 3.55(7)°, respectively. The bond lengths of N1-C1 and N1-O1 are 1.298(3) Å and 1.382(3) Å, respectively, which confirms that the bond between N1 and O1 atoms is oxime bond. The torsion angle of O1-N1-C1-C2 is 0.1(3)°, which reveals that the configuration of the double bond is Z.

Figure 1: The asymmetrical unit of title compound. Color codes: C (gray), O (red), N (blue), and H (green). Displacement ellipsoids are drawn at the 30% probability level.

The supramolecular structure of title compound is formed by various weak interactions, including hydrogen-bonding interactions, aromatic π-π stacking interactions, and C-Hπ interactions. As shown in Figure 2, a dimer is formed through intermolecular hydrogen bond O1-H1. The 3D supramolecular structure is obtained through the weak hydrogen bonds (C14-H14, C17-H17, C19-H19). The structure expanded along an axis through C17-H17 as a C(7) motif. The hydrogen-bonding data are shown in Table 2. In addition to the hydrogen-bonding interactions, there are several aromatic π-π stacking interactions. The aspectant π-π stacking can be found between ring-2 and ring-3 at position and between ring-1 and ring-3 at position. The structure expanded along direction through O1-H1 and Cg3. The π-π stacking interactions are characterized by a centroid to centroid distance of 3.794(3) Å with ring slippage of 1.102 Å for the former and a centroid to centroid distance of 4.044(2) Å for the latter. The π-π stacking data are presented in Table 2.

Table 2: Hydrogen-bond and C-Hπ geometry (Å and °) of title .

D-HA
D-H
HA
DA
D-HA

O1-H1N3
0.849(10)
1.976(10)
2.825(2)
177(4)

C14-H14O1
0.93
2.58
3.471(3)
162.1

C17-H17O1
0.93
2.71
3.517(3)
146.3

C19-H19O3
0.93
2.61
3.287(3)
129.7

C11-H11Cg1
0.93
2.95
3.733(6)
143

Cg3Cg3
—
—
3.794(3)
0.0(3)

Cg2Cg3
—
—
4.044(2)
3.38(9)

code: (i) , , ; (ii) , , ; (iii) , , ; (iv) , , ; (v) , , . Cg is the centroid of the ring.

Figure 2: The intermolecular hydrogen-bonding interactions (green-dashed lines), π-π stacking interactions (carmine-dashed lines), and C-Hπ interactions (orange-dashed lines) in title compound. Uncorrelated hydrogen atoms are omitted for clarity. H-Atoms on C11, C14, C17, and C19 are H11, H14, H17, and H19, respectively. Color codes: C (black), O (red), N (blue), and H (white).

3.2. Optimized Geometries

It should be mentioned that the optimized hydrogen bonds are 13%–17% longer than those in the corresponding crystal structure of Z isomer monomer. This phenomenon is common, as the position of the hydrogen atoms always cannot be measured correctly. The distances between nonhydrogen atoms are in consistency with the experimental values, whose deviation does not exceed 0.02 Å. The bond angles are also in consistency with the corresponding experimental values with a deviation smaller than 2° (Table 3). The differences between some calculated and experimental dihedral angles are much large, which can be attributed to the different treated methods. The optimization is performed in vacuo, and the dihedral angles are quite flexible. However, in the crystal structures, all the bonds and bond angles would be restricted obviously. Figure 3 shows the plots of Z and E configurations of title compound. The positions of their hydroxyl groups are different, and other geometries remain almost unchanged except for the dihedral angles between benzene rings. The bond angles and bond lengths of optimized structure of Z dimer are very consistent with those of the structural parameters of single crystal.

Table 3: Selected experimental and theoretical geometrical parameters of title compound.

Bond type
Calculated data
Experimental data
Atomic number

Monomer
Dimer

O-H
0.96366
0.82034
0.85(1)
O(1)-H(1)

N-O
1.40710
1.38079
1.382(3)
N(1)-O(1)

C=N
1.27849
1.29849
1.298(3)
C(1)-N(1)

C-C
1.51621
1.50113
1.502(3)
C(1)-C(2)

C-C
1.48329
1.47365
1.472(3)
C(1)-C(15)

N-N
1.33374
1.34626
1.347(2)
N(2)-N(3)

C=O
1.21574
1.19728
1.198(3)
C(6)-O(3)

C-O
1.34447
1.32638
1.326(3)
C(6)-O(2)

C-N
1.37160
1.35496
1.355(2)
C(5)-N(2)

C=N
1.34089
1.33854
1.338(2)
C(3)-N(3)

C-C
1.51153
1.48621
1.486(5)
C(7)-C(8)

C=C
1.38405
1.37172
1.371(3)
C(4)-C(5)

C=C
1.40663
1.39875
1.398(3)
C(3)-C(4)

N-O-H
102.132
109.472
104(2)
N(1)-O(1)-H(1)

O-N=C
113.709
111.984
111.8(2)
C(1)-N(1)-O(1)

N=C-C
116.326
114.806
114.7(2)
N(1)-C(1)-C(15)

C-C-C
120.278
120.387
120.5(2)
C(16)-C(15)-C(1)

N=C-C
123.507
123.628
123.7(2)
N(1)-C(1)-C(2)

C-C-N
111.805
112.751
112.8(2)
N(2)-C(2)-C(1)

C-N-N
118.218
118.666
118.6(2)
N(3)-N(2)-C(2)

N-C-C
124.600
123.927
123.9(2)
N(2)-C(5)-C(6)

C-C-O
125.689
126.451
126.4(2)
O(3)-C(6)-C(5)

C-C-O
110.964
109.733
109.8(2)
O(2)-C(6)-C(5)

N-C-C
120.850
121.304
121.3(2)
N(3)-C(3)-C(9)

Figure 3: The plots of Z and E configurations of title compound.

The selected experimental and theoretical geometrical parameters of monomer and dimer of Z isomer are shown in Table 3. Compared with the crystal structure values, the bond lengths N(1)-O(1) and O(1)-H(1) of Z isomer dimer have a deviation of about 0.02 Å; the band angle N(1)-O(1)-H(1) of Z isomer dimer has a deviation of about 5°. Except for the above deviations, the other bond lengths of title compound are almost the same. The causes of deviations may be ascribed to the intermolecular hydrogen bonding in crystal.

3.3. Isomers’ Energies

The single-point energies of Z and E isomers are −1163.366685225933 Hartree and −1163.365234699443 Hartree at the high-accurate PWPB95-D3/def2-QZVP level of theory, respectively. Therefore, the conformational energy of E isomer is 3.81 kJ/mol higher than that of Z isomer. Based on the Gibbs energy correction by the above method, the Gibbs energy difference between Z and E isomers is 3.61 kJ/mol. Therefore, the equilibrium constant between Z and E isomers can be calculated by employing the following equation: . At 298.15 K, the equilibrium constant is calculated to be 0.2331, which corresponds to the [Z]/[E] ratio of 4.29. It is close to the value of 3.7 that is estimated according to the experimental NMR spectra [3].

3.4. Transition State

The transition state (TS) (Figure 4) of the Z-E tautomerism is located. Obviously, the hydroxyl C-N-O is almost in a straight line with a bond angle of 177.85°. The TS energy is 224.88 kJ/mol and 228.69 kJ/mol higher than that of E and Z isomers, respectively.

Figure 4: The plot of TS configuration of title compound.

3.5. Frontier Orbital Energy Analysis and Molecular Total Energies

The frontier orbital energy levels and molecular total energies are demonstrated in Table 4, and their corresponding frontier orbital plots are shown in Figure 6, which are calculated at the B3LYP/def-TZVP level of theory. The most important orbitals in title compound are the frontier molecular orbitals, which are called LUMO and HOMO. The LUMO stands for an ability to get an electron, and HOMO expresses the ability to donate an electron [27] as an electron acceptor. Therefore, they are the most significant factors that influence the bioactivity, and the study on the frontier orbital energy can offer helpful information about the biological mechanism. The energy gap between HOMO and LUMO is 0.17039 Hartree, which indicates that this molecule is chemically inactive, as it cannot be easily excited by small energy. In Figure 6, HOMO is basically on phenyl diazine part of molecule, while LUMO is almost localized on the whole molecule except for the phenyl ring of phenyl diazine part. Therefore, the electrons transit from the phenyl diazine part to almost whole molecule. The HOMO-LUMO gaps of Z and E isomers are obtained at the B3LYP/def-TZVP level of theory, where their HOMO-LUMO gaps are 4.637 eV and 4.665 eV, respectively.

Table 4: Total energy and frontier orbital .

Energy

/Hartree
−1163.83

/Hartree
−0.26884

/Hartree
−0.25915

/Hartree
−0.25683

/Hartree
−0.24092

/Hartree
−0.22923

0.17039

/Hartree
−0.05884

/Hartree
−0.04354

/Hartree
−0.03342

/Hartree
−0.01384

/Hartree
−0.01049

Hartree = 4.35974417 × 10⁻¹⁸; 1 J = 27.2113845 eV;

Table 5 shows the composition (%) of atomic orbitals in HOMO and LUMO orbitals, respectively, where the percentage over 0.5% is shown. Obviously, the atomic orbitals of N1, C2, N3, C11, C14, C25, C29, and C31 occupy the most composition (%) of HOMO and N1, N3, O7, C10, C14, O16, C17, and C18 occupy almost the whole composition (%) of LUMO. The atomic number in molecule is shown in Figure 5. The understanding of the composition (%) of atomic orbitals in molecule would be helpful for the studies on the electron- and/or charge-transfer character.

Table 5: The composition (%) of atomic orbitals in HOMO or LUMO of molecule.

Atom
HOMO percentage
LUMO percentage

N1
7.00%
13.94%

N3
11.41%
11.12%

N8
0.88%
0.75%

C9
0.52%
0.96%

C10
2.26%
11.37%

C29
8.81%
0.92%

C11
11.02%
0.66%

C14
8.30%
12.51%

O16
0.61%
13.50%

H22
0.88%
1.07%

C33
3.25%
0.76%

Figure 5: Optimized geometry of title compound and its atomic numbers.

Figure 6: Representative frontier molecular orbital plots of molecule.

3.6. Mulliken Atomic Charges and Molecular Electrostatic Potential

The computed Mulliken atomic charges of all atoms of title compound are shown in Table 6. There are small charges in terms of all the nitrogen atoms. All the oxygen atoms (O4, O7, and O16) and some of carbon atoms (C2, C20, C23, and C42) have large negatively charged ones, which can have a facile interaction with the positively charged part of the receptor. The atom with the most negative charge is O16, followed by C42 and O4. Therefore, C42 is ready for electrophilic attack. All carbon atoms have negative charges except for the cases of C17 and C11, which would be responsible for nucleophilic attack sites.

Table 6: Mulliken atomic charges of title compound.

Atom
Charge

N1
−0.061300

C2
−0.147918

N3
0.057386

O4
−0.302941

H5
0.288510

C6
−0.002504

O7
−0.186096

N8
0.013476

C9
−0.072248

C10
−0.047429

C11
0.016290

C12
−0.111760

C25
−0.096077

H26
0.106782

C27
−0.093525

H28
0.109823

C29
−0.138813

H30
0.130015

C31
−0.105148

H32
0.106597

C33
−0.087616

H34
0.102026

C35
−0.070456

H13
0.174212

C14
−0.014417

H15
0.160320

O16
−0.355410

C17
0.105030

C18
−0.092596

H19
0.120824

C20
−0.227730

H21
0.199381

H22
0.183246

C23
−0.130862

H24
0.107364

C37
−0.054152

H38
0.129783

H39
0.127008

C40
−0.097878

H41
0.110371

C42
−0.327594

H43
0.126005

H44
0.119205

H45
0.124974

H36
0.105841

Electrostatic potential (ESP) on molecular van der Waals (vdW) surface is important to study and predict the intermolecular interactions [28–30]. It is well known that intermolecular recognition and electrostatic dominated noncovalent interactions tend to occur in terms of electrostatically complementary manner. The investigation of ESP of title compound is supposed to be useful to understand its intermolecular interactions. The ESP mapped vdW surface coped with surface extrema of Z isomer is shown in Figure 7(a).

(a)

(b)

Figure 7: (a) ESP mapped molecular van der Waals surface of Z isomer: the unit is in kcal/mol. Maxima of ESP and surface local minima are marked as orange and green spheres, respectively. The values of some prominent surface extremes are labeled on the graph, and the italic font labels global maximum and minimum. (b) An illustration of intermolecular H-bond in the crystal of Z isomer: only a pair of dimers is shown for clarity, and the red-dashed lines correspond to the H-bonds.

From Figure 7(a), it can be seen that the global maximum of ESP on van der Waals surface is present in the orange sphere with +44.3 kcal/mol, which is the direct consequence of the significantly positively charged hydroxyl hydrogen. Such a high value of ESP indicates that this hydroxyl group must be able to behave as a strong H-bond donor. In Figure 7(a), there are two notable regions where the color is deep blue, and both of them have the corresponding ESP minimum, which reflect the remarkable negative contribution to ESP due to the lone pair of the nitrogens. As highlighted by italic font in the figure, the most negative point is the green sphere containing nitrogen atom in the five-membered ring, implying that it should be the most favorable site to form H-bonds or interact with Lewis acids. Since this compound simultaneously possesses a strong H-bond donor and various H-bond acceptors, it is expected that the hydrogen bonds exist in its molecular crystal. Indeed, the actual crystal structure nicely validated this point. As shown in Figure 7(b), two pairs of equivalent hydrogen bonds can be clearly observed between two adjacent molecules in the crystal structure. Furthermore, we note that each of the hydrogen bonds can be approximately viewed as the result of fusion of global maximum and minimum of ESP, exhibiting that the formation of the hydrogen bonds completely follows the maximal electrostatic complementarity principle. The hydrogen bond length and bond angle are 2.82429 Å and 170.793°, respectively, which are very consistent with the corresponding structural parameters of single crystal in Table 3.

4. Conclusion

In summary, ethyl 1-(2-(hydroxyimino)-2-phenylethyl)-3-phenyl-1H-pyrazole-5-carboxylate is synthesized and featured by NMR, IR, MS, and X-ray single crystal diffraction. The geometries of Z and E isomers and the corresponding transition state (TS), as well as the dimer of Z isomer, are optimized using B3LYP hybrid functional coupled with def-TZVP triple-zeta polarized basis set. The bond angles and bond lengths of optimized structure of Z dimer are almost in full accord with crystal structural parameters. The combination of the weak intermolecular C-HO and N-HO hydrogen-bonding interactions stabilizes the crystal packing. The calculated hydrogen bond length and bond angle are 2.82429 Å and 170.793°, respectively, which is very consistent with the corresponding structural parameters of single crystal. According to the crystal structure, the torsion angle of H1-O1-N1-C1 is 161(2)°, which reveals that the dominant configuration of double bond is Z configuration. Double-hybrid functional PWPB95-D3 in combination with very highly accurate basis set def2-QZVP is employed to evaluate the accurate energy of each isomer and TS. The calculated equilibrium constant between Z and E isomers is 4.29 ([Z]/[E]), which is close to the value of 3.7 that is estimated according to the experimental NMR spectra.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this work.

Acknowledgments

The theoretical calculations are performed in the laboratory of Professor Yuxiang Bu. In addition, the authors thank Professor Baoxiang Zhao very much for his helpful discussions. Their institution is School of Chemistry and Chemical Engineering, University of Shandong.

References

H. Dai, J. Chen, H. Li et al., “Synthesis and bioactivities of novel pyrazole oxime derivatives containing a 5-trifluoromethylpyridyl moiety,” Molecules, vol. 21, no. 3, article no. 276, 2016.

Y. Shi, S. Wang, H. He et al., “Synthesis and bioactivity of novel pyrazole oxime ester derivatives containing furan moiety,” Chinese Journal of Organic Chemistry, vol. 35, no. 8, pp. 1785–1791, 2015.

L.-W. Zheng, Y. Li, D. Ge et al., “Synthesis of novel oxime-containing pyrazole derivatives and discovery of regulators for apoptosis and autophagy in A549 lung cancer cells,” Bioorganic & Medicinal Chemistry Letters, vol. 20, no. 16, pp. 4766–4770, 2010.

H. Dai, S. Ge, G. Li et al., “Synthesis and bioactivities of novel pyrazole oxime derivatives containing a 1,2,3-thiadiazole moiety,” Bioorganic & Medicinal Chemistry Letters, vol. 26, no. 18, pp. 4504–4507, 2016.

H.-J. Park, K. Lee, S.-J. Park et al., “Identification of antitumor activity of pyrazole oxime ethers,” Bioorganic & Medicinal Chemistry Letters, vol. 15, no. 13, pp. 3307–3312, 2005.

S.-L. Wang, Y.-J. Shi, H.-B. He, Y. Li, Y. Li, and H. Dai, “Synthesis and bioactivity of novel pyrazole oxime derivatives containing oxazole ring,” Chinese Chemical Letters, vol. 26, no. 6, article no. 3284, pp. 672–674, 2015.

G. Ouyang, X.-J. Cai, Z. Chen et al., “Synthesis and antiviral activities of pyrazole derivatives containing an oxime moiety,” Journal of Agricultural and Food Chemistry, vol. 56, no. 21, pp. 10160–10167, 2008.

G. Ouyang, Z. Chen, X.-J. Cai et al., “Synthesis and antiviral activity of novel pyrazole derivatives containing oxime esters group,” Bioorganic and Medicinal Chemistry, vol. 16, no. 22, pp. 9699–9707, 2008.

X. H. Liu, L. P. Zhi, B. A. Song, and H. L. Xu, “Synthesis, characterization and antibacterial activity of new 5-aryl pyrazole oxime ester derivatives,” Chemical Research in Chinese Universities, vol. 24, no. 4, pp. 454–458, 2008.

B. Sangepu, B. Gandu, G. Anupoju, and V. Jetti, “Synthesis of Isoxazole, 1, 2, 4-Oxadiazole and (1H-Pyrazol-4-yl)-methanone Oxime Derivatives from N-Hydroxy-1H-pyrazole-4-carbimidoyl Chloride and their Biological Activity,” Journal of Heterocyclic Chemistry, vol. 53, no. 3, pp. 754–761, 2016.

S. N. Shelke, C. H. Gill, M. S. More, S. B. Kale, S. M. Sonawane, and B. K. Karale, “Synthesis of novel benzo[d]oxazole derivatives as antimicrobial agents,” Organic Chemistry: An Indian Journal, vol. 3, no. 1, pp. 40–44, 2007.

R. R. Ranatunge, M. Augustyniak, U. K. Bandarage et al., “Synthesis and Selective Cyclooxygenase-2 Inhibitory Activity of a Series of Novel, Nitric Oxide Donor-Containing Pyrazoles,” Journal of Medicinal Chemistry, vol. 47, no. 9, pp. 2180–2193, 2004.

L.-W. Zheng, Z.-L. Gong, W.-L. Liu, Y.-R. Liu, and B.-X. Zhao, “Synthesis, X-ray crystal structure and fluorescent spectra of novel pyrazolo[1,5-a]pyrazin-4(5H)-one derivatives,” Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, vol. 81, no. 1, pp. 372–379, 2011.

D. S. Bolotin, M. V. Il'in, I. E. Kolesnikov et al., “Fluorescent (pyrazolyl acetoxime)ZnIIcomplexes: Synthetic, structural, and photophysical studies,” Inorganica Chimica Acta, vol. 455, pp. 9–14, 2017.

S. Viveka, G. Vasantha, Dinesha, S. Naveen, N. K. Lokanath, and G. K. Nagaraja, “Structural, spectral, and theoretical investigations of 5-methyl-1-phenyl-1H-pyrazole-4-carboxylic acid,” Research on Chemical Intermediates, vol. 42, no. 5, pp. 4497–4511, 2016.

A. Barakat, N. Shivalingegowda, N. Renuka et al., “Crystal structure of 3-(thiophen-2-yl)-5-(p-tolyl)-4,5-dihydro-1H-pyrazole-1-carboxamide,” Zeitschrift fur Kristallographie - New Crystal Structures, vol. 231, no. 1, pp. 267–269, 2016.

O. V. Dolomanov, L. J. Bourhis, R. J. Gildea, J. A. K. Howard, and H. Puschmann, “OLEX2: a complete structure solution, refinement and analysis program,” Journal of Applied Crystallography, vol. 42, no. 2, pp. 339–341, 2009.

G. M. Sheldrick, “Crystal structure refinement with SHELXL,” Acta Crystallographica Section C: Structural. Chemistry, vol. 71, pp. 3–8, 2015.

P. J. Stephens, F. J. Devlin, C. F. Chabalowski, and M. J. Frisch, “Ab Initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields,” The Journal of Physical Chemistry C, vol. 98, no. 45, pp. 11623–11627, 1994.

A. Schäfer, C. Huber, and R. Ahlrichs, “Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr,” The Journal of Chemical Physics, vol. 100, no. 8, pp. 5829–5835, 1994.

M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian 03, Gaussian, Pittsburgh, PA, USA, 2003.

L. Goerigk and S. Grimme, “Efficient and accurate double-hybrid-meta-GGA density functionals- evaluation with the extended GMTKN30 database for general main group thermochemistry, kinetics, and noncovalent interactions,” Journal of Chemical Theory and Computation, vol. 7, no. 2, pp. 291–309, 2011.

D. Rappoport and F. Furche, “Property-optimized Gaussian basis sets for molecular response calculations,” The Journal of Chemical Physics, vol. 133, no. 13, Article ID 134105, 2010.

F. Weigend, “A fully direct RI-HF algorithm: Implementation, optimised auxiliary basis sets, demonstration of accuracy and efficiency,” Physical Chemistry Chemical Physics, vol. 4, no. 18, pp. 4285–4291, 2002.

T. Lu and F. Chen, “Multiwfn: a multifunctional wavefunction analyzer,” Journal of Computational Chemistry, vol. 33, no. 5, pp. 580–592, 2012.

W. Humphrey, A. Dalke, and K. Schulten, “VMD: visual molecular dynamics,” Journal of Molecular Graphics, vol. 14, no. 1, pp. 33–38, 1996.

H. Becker and J. Fleming, “Frontier Orbitals and Organic Chemical Reactions,” Journal Für Praktische Chemie, vol. 320, no. 5, pp. 879-880, 1978.

J. S. Murray and P. Politzer, “Electrostatic Potentials: Chemical Applications,” in Encyclopedia of Computational Chemistry, vol. 2, pp. 912–920, John Wiley and Sons, New York, USA, 1998.

J. S. Murray and P. Politzer, “The electrostatic potential: an overview,” Wiley Interdisciplinary Reviews: Computational Molecular Science, vol. 1, no. 2, pp. 153–163, 2011.

P. Politzer and J. S. Murray, “Molecular Electrostatic Potentials and Chemical Reactivity,” in In Reviews in Computational Chemistry, vol. 2, pp. 273–312, John Wiley and Sons, New York, USA, 1991.