Synthesis, Characterization, and DFT Studies of N -(3,5- Bis (triﬂuoromethyl)benzyl)stearamide

: The novel N -(3,5- bis (triﬂuoromethyl)benzyl)stearamide 3 was prepared in moderate yield by a solventless direct amidation reaction of stearic acid 1 with 3,5- bis (triﬂuoromethyl)benzylamine 2 at 140 ◦ C for 24 h under metal- and catalyst-free conditions. This practical method was conducted in air without any special treatment or activation. The fatty acid amide 3 was fully characterized by IR, UV–Vis, 1D and 2D NMR spectroscopy, mass spectrometry, and elemental analysis. More-over, molecular electrostatic potential studies, determination of quantum descriptors, fundamental vibrational frequencies, and intensity of vibrational bands were computed by density functional theory (DFT) using the B3LYP method with 6-311+G(d,p) basis set in gas phase. Simulation of the infrared spectrum using the results of these calculations led to good agreement with the observed spectral patterns. metal- and catalyst-free direct amidation.


Introduction
The amide bond is one of the most important and fascinating functional groups in organic chemistry and biochemistry owing to its widespread occurrence in natural products, peptides, proteins, and a plethora of other biomolecules [1,2]. Medicinally, the amide-forming reactions are pivotal in the pharmaceutical industry and medicinal chemistry for the preparation of active pharmaceutical ingredients [3,4], insecticides [5,6], polymers [7], and a vast number of bioactive molecules [8][9][10][11][12]. Most notably, the amide functional group is present in approximately a quarter of clinically approved drugs and two-thirds of all drug candidates [13]. In particular, the fatty acid amides (FAAs) have received augmented interest due to their wide range of physiological and pharmacological activities. Members of FAAs include the N-acyltaurines such as N-arachidonoyltaurine that activates TRPV1 and TRPV4 calcium channels of the kidney [14]; the N-acylethanolamines (NAEs) such as N-arachidonoylethanolamine (anandamide), an endogenous ligand for the cannabinoid CB1 and CB2 receptors in the mammalian brain [15,16]; the N-acylamino acids (NAAs) such as N-arachidonylglycine that suppresses pain via a peripheral action [17]; and the primary fatty acid amides (PFAMs) such as oleamide that shows an important therapeutic potential in the management of sleep disorders and pain [18,19], as illustrated in Figure 1.
The procedures most used in the pharmaceutical industry for amide synthesis involve the acylation reaction of activated carboxylic acid derivatives such as acyl chlorides, anhydrides, or esters with amines, as well as the direct amidation of carboxylic acids with amines employing stoichiometric quantities of diverse coupling reagents, which generate large quantities of waste leading to poor atom-economy and tricky purification procedures [20][21][22]. Hence, in the past two decades, there was a significant increase in Molbank 2021, 2021, M1215 2 of 9 the development of sustainable catalytic methods for amide bond formation under mild conditions with a broad synthetic scope [23][24][25][26]. In particular, we recently reported a solventless direct amidation of non-activated carboxylic acids with amines employing a biogenic CuO−CaCO 3 catalyst for amide synthesis in moderate to high yields under normal atmospheric conditions [27]. Despite the elegant progress that has been accomplished, most catalytic direct amidation methods are unattractive for large-scale production owing to the limited substrate scope, low energy efficiency, poor sustainability, and not industrially relevant levels of recyclability of most current catalysts [23,24]. Importantly, the metal-and catalyst-free direct amidation of carboxylic acids with amines has been neglected for many years owing to its high temperatures and long reaction times; nonetheless, this thermal approach can eliminate some of the previously reported synthetic drawbacks. In connection with the ongoing development of efficient and simple protocols for the construction of carbon-nitrogen bonds [27][28][29][30][31] and our current studies on the synthetic utility of benzylamine derivatives [32,33], we describe a solventless synthesis of N-(3,5bis(trifluoromethyl)benzyl)stearamide 3 through an uncatalyzed process under normal atmospheric conditions. Moreover, molecular electrostatic potential studies, determination of quantum descriptors, fundamental vibrational frequencies, and intensity of vibrational bands were computed by the B3LYP method using 6-311+G(d,p) basis set in gas phase. The procedures most used in the pharmaceutical industry for amide synthesis involve the acylation reaction of activated carboxylic acid derivatives such as acyl chlorides, anhydrides, or esters with amines, as well as the direct amidation of carboxylic acids with amines employing stoichiometric quantities of diverse coupling reagents, which generate large quantities of waste leading to poor atom-economy and tricky purification procedures [20][21][22]. Hence, in the past two decades, there was a significant increase in the development of sustainable catalytic methods for amide bond formation under mild conditions with a broad synthetic scope [23][24][25][26]. In particular, we recently reported a solventless direct amidation of non-activated carboxylic acids with amines employing a biogenic CuO−CaCO3 catalyst for amide synthesis in moderate to high yields under normal atmospheric conditions [27]. Despite the elegant progress that has been accomplished, most catalytic direct amidation methods are unattractive for large-scale production owing to the limited substrate scope, low energy efficiency, poor sustainability, and not industrially relevant levels of recyclability of most current catalysts [23,24]. Importantly, the metal-and catalyst-free direct amidation of carboxylic acids with amines has been neglected for many years owing to its high temperatures and long reaction times; nonetheless, this thermal approach can eliminate some of the previously reported synthetic drawbacks. In connection with the ongoing development of efficient and simple protocols for the construction of carbon-nitrogen bonds [27][28][29][30][31] and our current studies on the synthetic utility of benzylamine derivatives [32,33], we describe a solventless synthesis of N-(3,5-bis(trifluoromethyl)benzyl)stearamide 3 through an uncatalyzed process under normal atmospheric conditions. Moreover, molecular electrostatic potential studies, determination of quantum descriptors, fundamental

Synthesis
We report a metal-and catalyst-free synthesis of N-(3,5-bis(trifluoromethyl)benzyl) stearamide 3 by a direct amidation reaction between equimolar amounts of stearic acid 1 with 3,5-bis(trifluoromethyl)benzylamine 2 at 140 • C for 24 h under solvent-free conditions, as depicted in Scheme 1. It should be noted that the water vapor condensed on the walls of the open-topped tube was removed with a small piece of cotton attached to a spatula to displace the equilibrium towards amide formation. After the specified reaction time, the mixture was allowed to cool to ambient temperature, and the resulting crude product was purified by flash chromatography on silica gel using a mixture of DCM/MeOH (100:1, v/v) as eluent to furnish the fatty acid amide 3 in 41% yield. Conceptually, the condensation reaction between a carboxylic acid and an amine is the most attractive and practical synthetic approach to form the amide bond via the expulsion of a molecule of water. Nonetheless, the amide 3 was isolated in moderate yield (41%) because the innate acidity of the stearic acid and the basicity of the 3,5-bis(trifluoromethyl)benzylamine can generate ionized forms that limit its thermal condensation. Despite the moderate yield, the collateral formation of substituted ammonium carboxylate salt was not detected by thin-layer chromatography (TLC) and 1 H NMR spectroscopy. In consequence, the releasing of water as the unique by-product gives to this approach a remarkable eco-friendly quality. It should be noted that the synthesis and full characterization of amide 3 have not been reported on Reaxys and SciFinder databases. For that reason, a full spectroscopic and analytical characterization was performed in this work (see Section 3 and Figures S1-S5). The IR, 1D NMR, and MS spectra and elemental analysis suggest that the structure of the isolated white solid corresponds to the amide 3. Furthermore, 2D HSQC, HMBC, and COSY experiments permit the assignment of all proton and carbon atoms, confirming the proposed structure for amide 3 without ambiguity (see Supplementary Material).
the most attractive and practical synthetic approach to form the amide bond via the expulsion of a molecule of water. Nonetheless, the amide 3 was isolated in moderate yield (41%) because the innate acidity of the stearic acid and the basicity of the 3,5bis(trifluoromethyl)benzylamine can generate ionized forms that limit its thermal condensation. Despite the moderate yield, the collateral formation of substituted ammonium carboxylate salt was not detected by thin-layer chromatography (TLC) and 1 H NMR spectroscopy. In consequence, the releasing of water as the unique by-product gives to this approach a remarkable eco-friendly quality. It should be noted that the synthesis and full characterization of amide 3 have not been reported on Reaxys and SciFinder databases. For that reason, a full spectroscopic and analytical characterization was performed in this work (see Section 3 and Figure S1-S5). The IR, 1D NMR, and MS spectra and elemental analysis suggest that the structure of the isolated white solid corresponds to the amide 3. Furthermore, 2D HSQC, HMBC, and COSY experiments permit the assignment of all proton and carbon atoms, confirming the proposed structure for amide 3 without ambiguity (see Supplementary Material). As expected, some aromatic carbons appeared as quartets with intensities 1:3:3:1 due to the splitting caused by the fluorine atoms on trifluoromethyl groups attached to the benzene ring. A molecular ion with m/z 509 in the mass spectrum is also consistent with the structure 3. Ultimately, two well-defined absorptions are noted in the UV-Vis spectrum recorded in acetonitrile at 205.0 and 265.7 nm corresponding to π→π*(C=O)) and n→π*(C=O)) transitions, respectively (see Supplementary Material).

Molecular Electrostatic Potential (MEP) Mapping
In order to study the electronic behavior of amide 3, the geometric optimization was computed using density functional theory (DFT), employing Becke's three-parameter hybrid function with the nonlocal correlation of Lee-Yang-Parr (B3LYP) method at 6-311+G(d,p) basis set in gas phase [34,35]. The atomic coordinates of the optimized structure of 3 are given in Table S1. Afterward, the molecular electrostatic potential (MEP) was calculated to find the charge distribution on the surface of 3, as well as to determine the sites with higher and lower electrostatic potential, as shown in Figure 2. It is important to mention that electrostatic potentials at the surface are represented by different colors, such as red, electron-rich; blue, electron-deficient; light blue, slightly electron-deficient; and yellow, slightly electron-rich. The MEP analysis for amide 3 showed an electron-rich region around the carbonyl oxygen, whereas an electron-deficient region is located around the amide nitrogen; these findings are in good agreement with the resonance As expected, some aromatic carbons appeared as quartets with intensities 1:3:3:1 due to the splitting caused by the fluorine atoms on trifluoromethyl groups attached to the benzene ring. A molecular ion with m/z 509 in the mass spectrum is also consistent with the structure 3. Ultimately, two well-defined absorptions are noted in the UV-Vis spectrum recorded in acetonitrile at 205.0 and 265.7 nm corresponding to π→π*(C=O)) and n→π*(C=O)) transitions, respectively (see Supplementary Material).

Molecular Electrostatic Potential (MEP) Mapping
In order to study the electronic behavior of amide 3, the geometric optimization was computed using density functional theory (DFT), employing Becke's three-parameter hybrid function with the nonlocal correlation of Lee-Yang-Parr (B3LYP) method at 6-311+G(d,p) basis set in gas phase [34,35]. The atomic coordinates of the optimized structure of 3 are given in Table S1. Afterward, the molecular electrostatic potential (MEP) was calculated to find the charge distribution on the surface of 3, as well as to determine the sites with higher and lower electrostatic potential, as shown in Figure 2. It is important to mention that electrostatic potentials at the surface are represented by different colors, such as red, electron-rich; blue, electron-deficient; light blue, slightly electron-deficient; and yellow, slightly electronrich. The MEP analysis for amide 3 showed an electron-rich region around the carbonyl oxygen, whereas an electron-deficient region is located around the amide nitrogen; these findings are in good agreement with the resonance hybrid contributions in amides [36]. In addition, two green regions located around both 3,5-bis(trifluoromethyl)phenyl moiety and long aliphatic chain showed potential sites for lipophilic interactions. hybrid contributions in amides [36]. In addition, two green regions located around both 3,5-bis(trifluoromethyl)phenyl moiety and long aliphatic chain showed potential sites for lipophilic interactions.

Frontier Orbitals and Global Reactivity Descriptors
The frontier molecular orbitals (FMO) energy levels highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) for amide 3 were computed using the B3LYP method with 6-311+G(d,p) basis set in gas phase. The HOMO and LUMO energies represent the ability to donate and gain an electron, respectively. The HOMO-LUMO energy gap and the pictorial illustration of HOMO and LUMO are shown in Figure 3. For the optimized structural models, the HOMO is located over the amide nitrogen (N2pz), which indicates that the majority of charge transfer occurs between amide nitrogen and carbonyl carbon, in line with the resonance model in amides [36]. Furthermore, the LUMO is mainly located over the 3,5-bis(trifluoromethyl)phenyl moiety with a minor contribution of the π*C=O orbital, which indicates the small contribution of carbonyl oxygen to the LUMO. As a result, the basic amide nitrogen and the acidic carbonyl group in amide 3 become less reactive, which is evidenced by the high HOMO-LUMO gap (5.54 eV). The calculated ionization potential (IP), electron affinity (EA), electrophilicity index (ω), chemical potential (μ), electronegativity (χ), and hardness (η) for amide 3 are illustrated in Table 1. The Koopmans' theorem in density functional theory has shown that

Frontier Orbitals and Global Reactivity Descriptors
The frontier molecular orbitals (FMO) energy levels highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) for amide 3 were computed using the B3LYP method with 6-311+G(d,p) basis set in gas phase. The HOMO and LUMO energies represent the ability to donate and gain an electron, respectively. The HOMO-LUMO energy gap and the pictorial illustration of HOMO and LUMO are shown in Figure 3. For the optimized structural models, the HOMO is located over the amide nitrogen (N2p z ), which indicates that the majority of charge transfer occurs between amide nitrogen and carbonyl carbon, in line with the resonance model in amides [36]. Furthermore, the LUMO is mainly located over the 3,5-bis(trifluoromethyl)phenyl moiety with a minor contribution of the π* C=O orbital, which indicates the small contribution of carbonyl oxygen to the LUMO. As a result, the basic amide nitrogen and the acidic carbonyl group in amide 3 become less reactive, which is evidenced by the high HOMO-LUMO gap (5.54 eV).
hybrid contributions in amides [36]. In addition, two green regions located around both 3,5-bis(trifluoromethyl)phenyl moiety and long aliphatic chain showed potential sites for lipophilic interactions.

Frontier Orbitals and Global Reactivity Descriptors
The frontier molecular orbitals (FMO) energy levels highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) for amide 3 were computed using the B3LYP method with 6-311+G(d,p) basis set in gas phase. The HOMO and LUMO energies represent the ability to donate and gain an electron, respectively. The HOMO-LUMO energy gap and the pictorial illustration of HOMO and LUMO are shown in Figure 3. For the optimized structural models, the HOMO is located over the amide nitrogen (N2pz), which indicates that the majority of charge transfer occurs between amide nitrogen and carbonyl carbon, in line with the resonance model in amides [36]. Furthermore, the LUMO is mainly located over the 3,5-bis(trifluoromethyl)phenyl moiety with a minor contribution of the π*C=O orbital, which indicates the small contribution of carbonyl oxygen to the LUMO. As a result, the basic amide nitrogen and the acidic carbonyl group in amide 3 become less reactive, which is evidenced by the high HOMO-LUMO gap (5.54 eV). The calculated ionization potential (IP), electron affinity (EA), electrophilicity index (ω), chemical potential (μ), electronegativity (χ), and hardness (η) for amide 3 are illustrated in Table 1. The Koopmans' theorem in density functional theory has shown that The calculated ionization potential (IP), electron affinity (EA), electrophilicity index (ω), chemical potential (µ), electronegativity (χ), and hardness (η) for amide 3 are illustrated in Table 1. The Koopmans' theorem in density functional theory has shown that the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) energies of typical molecules are compared with the ionization potential (IP) and electron affinity (EA), respectively [37]. The average value of HOMO and LUMO energies is related to the electronegativity (χ), as χ = (IP + EA)/2 [38]. The HOMO−LUMO energy gap is correlated to the hardness (η), which is a good indicator of the chemical stability [39]. Parr and coworkers incorporated the term electrophilicity index (ω), as ω = µ 2 /2η where µ is the chemical potential taking the average value µ = −(IP + EA)/2 [40]. The negative value of the chemical potential (µ = −4.61 eV) for amide 3 is associated with a good chemical stability resulting in a high melting point. The magnitude of the hardness (η) supported by the HOMO-LUMO energy gap is 5.54 eV, indicating a low chemical reactivity. Besides, the low HOMO energy (−7.38 eV) indicates the difficulty of electron donation of the amide 3 due to its high electronegativity value (χ = 4.61 eV). Electrophilicity is an important descriptor of reactivity that allows the quantitative classification of the global electrophilic nature of a molecule within a relative scale. Consequently, Domingo and coworkers established an electrophilicity (ω) scale for the classification of organic molecules as strong electrophiles with ω > 1.5 eV, moderate electrophiles with 0.8 < ω < 1.5 eV, and marginal electrophiles with ω < 0.8 eV [38]. As shown in Table 1, the amide 3 can be regarded as a strong electrophile (ω = 1.92 eV).

Vibrational Analysis
For visual comparison, a superposition of experimental and calculated [B3LYP/6-311+G(d,p)] FT-IR spectra is illustrated in Figure 4. The vibrational frequencies from DFT calculations are often overestimated and commonly scaled by empirical factors. A scaling factor of 1.0 was used in this study to fit the calculated frequencies to the experimental ones.

General Information
Stearic acid and 3,5-bis(trifluoromethyl)benzylamine were purchased from Sigma-Aldrich and used without further purification. Starting materials were weighed and handled in air at ambient temperature. Silica gel aluminum plates (Merck 60 F 254 , Darmstadt, Germany) were used for analytical TLC. The IR absorption spectrum was recorded at room temperature from 4000 to 500 cm −1 employing a Nicolet iS50 FTIR spectrometer (Thermo Fischer Scientific Inc., Madison, WI, USA) equipped with an attenuated reflectance accessory. 1

Computational Study
The geometry optimization of the N-(3,5-bis(trifluoromethyl)benzyl)stearamide 3 was performed using density functional theory (DFT), employing Becke's three-parameter hybrid function with the nonlocal correlation of Lee-Yang-Parr (B3LYP) method with 6-311+G(d,p) basis set in gas phase. Theoretical calculations were performed by DFT using Gaussian 16 [34,35,41]. The molecular visualization, optimized geometry, and molecular electrostatic potential were obtained using GaussView 6.0.16. The multiwfn algorithm was used for these calculations. It should be noted that DFT calculations were performed for the prediction of chemical stability and reactivity of amide 3 with the help of quantum chemical descriptors such as frontier molecular orbitals and HOMO-LUMO energy gap, as well as global reactivity descriptors such as potential (IP), electron affinity (EA), electrophilicity index (ω), chemical potential (µ), electronegativity (χ), and hardness (η). The molecular electrostatic potential (MEP) analysis is presented as a powerful tool for the knowledge of charge distribution, and its results can be useful in determining how amide 3 molecules could interact with each other.

Conclusions
In summary, we described an operationally simple and sustainable method for the thermal direct amidation of a nonactivated carboxylic acid with amine in the absence of a solvent, catalyst, coupling reagent, and drying agent under normal atmospheric conditions. The releasing of water as the unique by-product gives this approach a remarkable ecofriendly quality. Accomplishing the synthesis of the amide 3 contributed to expanding our previous amide library [27]. The MEP total density showed an electron-rich region around the carbonyl oxygen, whereas an electron-deficient region is located around the amide nitrogen. The frontier molecular orbital energies and global reactivity descriptors showed that the molecule 3 would have a low reactivity owing to its high HOMO-LUMO energy gap (5.54 eV). Besides, the low HOMO energy (−7.38 eV) would indicate the difficulty of electron donation due to its high electronegativity value (χ = 4.61 eV). Simulation of the infrared spectrum led to good agreement with the observed spectral patterns. Ultimately, the fatty acid amide 3 could be used as an analog of the endogenous cannabinoid receptor ligand arachidonylethanolamide (anandamide) with potential applications in medicinal chemistry and drug discovery.

Data Availability Statement:
The data presented in this study are available in this article.