Abstract

Unsaturated fatty acids have great interest by their activities as industrial materials in novel applications. In the present work, the cis-3-hexenoic acid (HA) adsorbed on the Ni(111) surface was studied by first-principles calculations using the Vienna Ab Initio Simulation Package (VASP). The most stable location for HA is presented on top site of Ni(111), although the energies are very similar for all the adsorption sites. The surface-molecule interaction takes place between the carboxyl group of HA and surrounding Ni atoms. The adsorption is weak and consequently the metal-molecule length is enhanced. The carboxyl group is elongated and weakened after adsorption giving rise to a shift in stretching frequencies. There are notable changes on the magnetic moments values of Ni surface atoms neighboring to the molecule that mainly induced magnetic moments on O and H atoms. Noticeable charge transfer occurs in 3d 4s, p Ni orbitals and 2s C, 2s p O, 1s H orbitals of carboxyl group. The surface presents positive work function changes after adsorption as a consequence of an electron back-donation. During interaction, the significance of the magnetic effects over the structural effects is evidenced. This sets the stage for a future adsorption process improvement based on the modification of the surface magnetic properties.

1. Introduction

The development of thin organic films on various substrates has attracted strong interest due to the use of organic films in novel electronic applications, such as optoelectronic devices, sensors, and field effect transistors [1]. Acids usually form close-packed monolayer films with chains highly oriented on substrates. Previous studies have reported the adsorption of stearic acid on gold [2], aluminum [3], and copper and silver [4]. However, there are no detailed reports of organic acids on magnetic metals.

Unsaturated fatty acids have attracted big interest by their behavior as medical and industrial materials. The knowledge of dipole moment of molecules to nanotechnology is related within the functionality of these compounds producing designed products for specific applications [5]. The design of materials in nanoscale can be rearranged through interactions such as electrostatic dipole, van der Waals interactions, and hydrophobic and hydrophilic interactions. An example of this arrangement is the induced self-aggregation where the arrangement of the molecules is driven under control through an external magnetic field or electrical field [6]. On the other hand, unsaturated fatty acids are widely distributed in biological tissues as a main constituent of biomembranes. These acids promote fluidity and control the phase transition behavior depending on the environmental condition [7]. In addition, in technology for new materials, hydrophobic nanoparticles of copper capped with oleic acid were prepared to prevent oxidation [8]. In spite of such importance, only a little is known about the molecular-order structures and the physicochemical properties of unsaturated acids on magnetic metals. The development of modern theoretical surface science provides an opportunity to investigate surfaces and adsorbate structures on the atomic scale with useful applications in industrial technologies. A detailed study of the interaction of organic films adsorbed on the metal substrates is usually focused in the surface geometry and the electronic properties. In addition, when the substrate is magnetic, many other effects can also be studied such as magnetic proximity effect between the molecule and the magnetic substrates [9]. Focusing on magnetic substrates, there is a vast chemical industry based on the catalytic properties of nickel. Nickel is one of the most versatile heterogeneous metal catalysis. Its chemical stability makes it a very convenient election in a lot of applications. The adsorption of molecular hydrogen on an ultrathin layer of Ni(111) hydride has been researched by Shan et al. [10]. Carrasco and colleagues have theoretically studied the adsorption of CO and C on Ni(111) surface [11]. The characterization of stearic acid adsorption on Ni(111) surface by experimental and first-principles study approach has been investigated by Liang et al. [12]. Focusing on the adsorbed organic acid, specifically, the cis-3-hexenoic acid (C5H9COOH) is an unsaturated acid that presents a double bond in the center of the chain and, due to their regular and simple organic structure, the obtained results can be used in the prediction of the adsorption and catalytic behavior to other industrial acids of longer chain (e.g., oleic acid) [13] with a significant reduction in computational requirements (e.g., system size and computational time).

In this work, the cis-3-hexenoic acid (HA) adsorbed on Ni(111) surface is studied by first-principles calculation. The physical-chemical, magnetic, and structural effects taking place during adsorption are analyzed. The theory and model is described in the next section.

2. Theory and Model

The calculations were performed in the framework of the Density Functional Theory (DFT) using the Vienna Ab Initio Simulation Package [14] applying Grimme’s D2 correction [15]. The electron projector augmented wave method was used and the generalized gradient approximation with the Perdew-Burke-Ernzerhof functional was employed. Cut-off energy of 400 eV was employed for the fixed convergence of the plane-wave expansion. This value is in agreement with that reported in [16, 17] and it is based on a previous test where the calculation error was lower than 0.01 eV. A set of 3 × 3 × 1 Monkhorst-Pack -points was used to sample the Brillouin Zone while the ground state was obtained using a Methfessel-Paxton smearing of 0.2 eV. The computations were performed at the spin-polarized level in order to evaluate the magnetic properties of the systems. The Ni(111) surface of FCC stacking layered structure was represented with a periodically repeated slab containing five layers of atoms separated in the normal direction by a vacuum region. The width of this gap was optimized to avoid the interaction between slabs. The lattice constants for Ni bulk were optimized obtaining a value of 3.53 Å that is in agreement with the experimental value [18]. A box of 20 × 20 × 20 Å3 was used to obtain the isolated molecule energy.

The molecule was placed on one side of the slab and its geometry was permitted to relax completely together with the uppermost three uppermost layers of the metallic slab. The adsorption energy () was calculated as the difference between the energy of the HA molecule adsorbed on Ni surface () and the sum of the free Ni surface () and the isolated HA molecule () energies:

The adsorption energies were calculated in order to compare the stability of the molecule adsorbed on the different surface sites. A negative value indicates an exothermic adsorption process.

In the present work, the binding preference site and the adsorption energies of HA molecule on the Ni(111) surface model are investigated. Some electronic characteristics as the density of orbital states (DOS), the changes in the work function, and the magnetic moment are evaluated. The vibrational properties are also calculated. The influence of spin polarization on adsorption with the hope of future applicability in the optimization of adsorption and catalytic behavior is analyzed.

3. Results and Discussion

3.1. Adsorption Energy and Site Preference

An accurate calculation of the adsorption energy for a given system is the starting point to understanding of the mechanism of catalytic activity. The adsorption energy describes the strength of the chemical bond between adsorbate and substrate. Different adsorption sites for cis-3-hexenoic acid (HA) on Ni(111) surface were considered; specifically, they are top, hollow, and bridge sites (see Figure 1). In Table 1, the equilibrium interatomic distances and the adsorption energies for HA adsorbed on the different sites of the Ni(111) surface are summarized. The internal geometry of the molecule does not change significantly in presence of the substrate, independently of the adsorption site. The energies are similar for all the sites; the more stable location for HA is presented on top Ni(111) site. Figure 2 shows lateral and top views of HA molecule adsorbed on top Ni(111) site. The cis-3-hexenoic acid is adsorbed on Ni(111) surface via carboxyl group. The acid group (COOH) of the molecule interacts with surrounding Ni atoms while the long chain extends away from the surface. The molecule-surface distances depend slightly on the adsorption site (see Table 1). In general, the height of the molecule on the surface is established by the expansion of the metal d orbitals. For the higher coordinated sites, the height of the molecule on the surface is also influenced by the nearest neighbor distance of the substrate atoms. The distance between the molecule and the nearest metal atoms decreases if the extension of the interacting metal orbitals decreases. At the same time, the molecule will adsorb at higher heights if the metal has a small lattice constant. In our studio, the HA molecule adsorbs at the distances shown in Table 1 as compromise between the small lattice constant of nickel (compared with other transition metals) and the influence of Ni d orbitals. The local density of states (LDOS) graphics of the s, p, and d band obtained from a spin-polarized calculation indicates that Ni surface exhibits dispersion of s and p bands much larger than d band, indicating the much more contracted nature of Ni d orbitals. The computed electronic configuration for the surface atoms of bare Ni(111) is d8.3 s0.5 p0.4. The electron occupation of the valence s orbitals seems to be slightly lower and that of the p orbitals slightly greater. Similar characteristics have been reported by Wong and Hoffmann [19], where the role of s, p Ni orbitals and the contracted nature of d orbitals were also emphasized in the discussion. The authors indicated that there is little doubt that the d orbitals are involved in adsorbate-surface bonding, but the combined effect of the s and p bands is much more important [19]. The changes in the DOS of these orbitals after adsorption are thus much less obvious than those of the d bands; the s and p bands of a transition metal are very diffuse and their band widths are tens of eV. The authors remark that besides the d bands, the s and p bands also play a significant role in the bonding [19]. As a result, the authors expose that the role of the s and p orbitals on chemisorption involving transition metals is seldom discussed, and they claim that these orbitals should interact with the adsorbate [19]. The contribution of Ni s, p and d orbitals has also been reported by Germán et al. for the cyclopentene (c-C5H8) dehydrogenation to cyclopentadienyl anion (c-C5H5) on Ni(111) [20], while Liang et al. have reported that Ni 4s orbital plays an important role in the adsorption of the stearic acid on Ni(111) [12]. Gupta remarks the importance of p orbital during H hydride formation in Ni metals [21]. A wavefunction analysis indicates the large delocalization of the states of s-p character, the narrow d band weak interacting atomic like states, and the importance of Ni p orbital in the bonding [21].

At this stage, it is useful to compare the adsorption geometry of noble metals. The noble metals present differences from other metals; this reveals the completely filled d band and consequently a different bonding mechanism. Adsorption on the noble metals is notably weaker and hence the metal-molecule lengths are increased. Even though nickel is not a noble metal, its behavior corresponds in some way because it presents a quasi-filled d band. According to calculations, is weak, decreases for on-top site, and increases for sites with higher coordination (see Table 1).

Adsorption may also correlate with the electronegativity of the adatoms, that is, with how strongly an atom will attract electrons from the surrounding atoms. In relationship, Ni atoms have more affinity to interact with H atom than O atoms of carboxyl group and, therefore, the H-surface distance is presented as the shortest. The electronegativity differences and, in consequence, the electron exchange contribute to the adsorption behavior (the charge density redistribution is analyzed in Section 3).

In next sections, in absence of a preferential adsorption site, the HA molecule located on top site of Ni(111) was selected to study the electronic, magnetic characteristics and the vibrational properties of the system.

3.2. Magnetic Moment

The dipole moment change was analyzed considering different contributions. As postulated above, the surface-molecule interaction takes place between the COOH group of the acid of the molecule and the surrounding Ni atoms. The local magnetic moments for the different atoms of our slab model were calculated. The atom number indicates the position of the atom in the Ni(111) slab: 1 (surface), 2 (interface), 3 (medium), and 4 (bulk) (see Figure 1). Table 2 shows the magnetic moment of Ni atoms and the OH atoms of carboxyl group. The results show that magnetic moments were mainly induced on the O and H atoms by proximity to the Ni atoms. The magnetic moment is enhanced for HA atoms neighboring to the Ni surface where the proximity effect is stronger. Small magnetic moment for O and H atoms of HA molecule is obtained. Conversely, the magnetic moments of interacting Ni() atoms increase at least by 8.27% with respect to the clean surface. Greater changes of the magnetic moments values of Ni atoms were obtained; it seems that the physical cause of this electron enrichment is inherent to the chemisorption bonding with the HA molecule.

In general, the magnetic moment for surface Ni() and interface Ni(2) atoms increases while the magnetic moment for medium Ni(3) and bulk Ni(4) atoms decreases after adsorption (see Table 2). On the other hand, the magnetic moment for surface Ni() atom of pure Ni differs from that of Ni(4) atom in the bulk of pure Ni: 0.643 μB in comparison with 0.654 μB for bulk Ni. This is a property that occurs for the presence of surface states and contracted bands. The surface Ni() atom in the Ni(111) slab has the largest Ni magnetic moment (0.701 μB) after HA adsorption, which is enhanced by 9.42% with respect to that of bulk Ni(4) (0.636 μB). It confirms the importance of band hybridization on the Ni magnetization. The magnetic moment of the C atom (of COOH group) also changes, in a fewer percent comparing with O and H atoms; the reasons are that C 2s, 2p population was changed by adsorption inducing a charge transfer and hybridization interaction. The HA adsorption on top site has the surface Ni atoms with lowest magnetic moments value (comparing adsorption on all the Ni(111) sites); thus, the system presents the better stability.

3.3. Electronic Structure: d Band Centre and Charge Density Redistribution

The electronic structure supplied important information about the interaction between adsorbate and substrate and their effect on the binding energy and the geometrical structure. The position of the d band center of the isolated surface is another important quality which is associated with the strength of the molecule-surface interaction. As the width of the d band increases, the energy is reduced and the distance of the molecule on the surface should increase. The d band center takes part in the bonding for many molecule-surface interactions where the major contact is due to hybridization of the adsorbate frontiers orbitals and the substrate d orbitals.

The density of states (DOS) for the HA adsorbed on top Ni(111) surface is shown in Figure 3(a). For comparison, the densities of states of both the clean Ni(111) surface (Figure 3(b)) and the isolated HA molecule (Figure 3(c)) were also calculated. In the spin-polarized calculations, the average between the Ni(111) d band centers obtained for up and down cases was taken from the Ni d band local density of states (LDOS) [22]. The Ni d band expands from −5 to −0.5 eV and, therefore, the center is located to about 4.5 eV. These values are in agreement with ab initio and semiempirical results reported in [2226]. The system’s total DOS nearby the Fermi energy is dominated by the Ni surface, but, in a deeper energy level, the system DOS is affected by the adsorbed molecule (Figure 3(a)). The location of d-states with respect to the molecule frontiers orbitals is an important factor to predict the facility to the metallic surface to form a chemical bond with the molecule. In this work, the molecule-surface overlapping occurs within −8 and −0.5 eV; Ni d-states and also s, p Ni states are involved in the interaction with HA. As the molecule approaches the surface, the localized HA orbitals are shifted to lower energies. The HA orbitals have a lesser task in the interaction with the surface because they are completely occupied and characterized by Pauli repulsion. In general, after adsorption, the changes on the DOS peaks increase with coordination: the higher the coordination is the greater the DOS peaks change. In our study, there are no significant changes because Figure 3 shows the molecule located on the lower coordination site.

As known, the difference in spin-up and spin-down states is the basis of the Ni magnetism. The spin unbalance decreases considerably upon HA adsorption. The formation of molecule-surface interactions increases the DOS with respect to Ni pure, showing a coupling between Ni orbitals and COOH orbitals, which induces a larger magnetic moment per Ni atom. Accordingly, there are electronic charge density changes in the substrate and adsorbate orbitals. The formation of the surface-HA bond exhibits electron rearrangement between the bonding species. The HA charge density is redistributed and some of it is transferred into antibonding orbitals, which weakens the COOH bond. The Ni surface had influence on the adsorbed molecular electronic properties; consequently, we have calculated the charge transfer during adsorption. The results show that the charge transfer mainly happens on the interface between carboxyl and Ni surface (see Table 3). For HA adsorption on top Ni(111), the charge transfer percentages are as follows: 4s Ni (between −0.61 and −1.20%), 4p Ni (between +1.86% and +3.45%), and 3d Ni (between +0.24% and 0.49%) for Ni nearest neighbors and 2s C (−0.47%), 2p C (−0.12%), 2s (+0.54%), 2p (−0.93%), 2s (−0.12%), 2p (−0.72%), and 1s H (−11.72%) of carboxyl group. The population of the second nearest Ni atom modifies a very little amount; this indicates a majority participation of the Ni first neighbor orbitals; the influence is mainly limited to the system interface.

3.4. Work Function

When the molecule is adsorbed on the surface the work function change is immediately associated with the change of the electric dipole in the surface and it can be calculated as the difference between the work function of the metal surface with and without the adsorbed molecule. The electron donation from molecule to surface is presumed to reduce the work function, while the back-donation from surface to molecule increases this property. Accordingly, metallic surface with positive work function change after adsorption presents an electron back-donation that conduces to a positive dipole moment. Evaluating the Ni surface work function upon the HA chemisorption with respect to the clean Ni surface, it is observed that for all the adsorption sites positive work function values are obtained. The increase of the work function is due to the downward dipole direction, induced by the HA molecule. The dipole moment change, produced by the HA on the Ni(111) surface, generates an extraction of the electron from the Ni surface. The work function values present a variation in a range of 0.6–0.9 eV. Small variation is evidenced due the absence of a preferential adsorption site.

3.5. Vibrational Frequencies

In general the adsorption energy for the molecule located on a surface stable site corresponds to the absolute minimum while the adsorption energies calculated for other high-symmetry sites correspond to a local minimum energy on the potential energy surface. In this work, we have corroborated that the obtained geometry for HA on Ni(111) surface is in true minima and not on the transition state because calculated HA vibrational frequencies were all real.

Table 4 shows the calculated stretching frequencies of the HA molecule on top Ni(111). The stretching vibration mode for the HA molecule changes after adsorption. The HA adsorption on Ni(111) surface produces a shift to lower frequency values for OH peak and higher frequency values for C=O peak of COOH. These changes in the stretching frequencies are due to the formation of the interaction between the molecule and the metallic surface. The corresponding molecule geometry shows small changes: elongation of the C=O and O–H bonds with respect to the isolated molecule while the C–C, C–O, and C=C bonds shrink after adsorption (see Table 5). The carboxyl group is elongated and weakened giving rise to a shift in stretching frequencies. Weaker and therefore longer C=O and O–H bonds imply molecule-metal bonding and result in the shift in C=O and O–H stretching frequencies. In general, the relationship of bond distance changes on the coordination is smallest for the 3d metals. The tendency is presented as follows: the higher surface reactivity, the stronger bonding and population of antibonding orbitals, reflecting more elongated bond lengths and shifting stretching frequencies.

4. Conclusions

In the present work, the cis-3-hexenoic acid (HA) adsorbed on Ni(111) surface was studied by first-principles calculations using the VASP code. The more stable location for HA is presented on top site of Ni(111), although the energies are very similar for all the adsorption sites. The surface-molecule interaction takes place between carboxyl group of HA molecule and the surrounding Ni atoms of surface. Adsorption is weak and the metal-molecule length is enhanced. The carboxyl group is elongated and weakened after adsorption giving rise to a shift in stretching frequencies related to the chemisorption bonding with HA molecule.

There are notable changes of the magnetic moments values of Ni surface atoms neighboring to the molecule that mainly induced magnetic moments on O and H atoms of carboxyl group by proximity. The physical cause of this electron enrichment is inherent to the bonding with the HA molecule. The Ni(111) surface depicts positive work function change after adsorption as a consequence of an electron back-donation. Noticeable charge transfer occurs in 4s p 3d Ni orbitals and 2s C, 2s p O, 1s H orbitals of carboxyl group.

During interaction, the significance of the magnetic effects over the structural effects is evidenced. This study sets the stage for a future adsorption process improvement based on the modification of the surface magnetic properties.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Our work was supported by SCyT UTN, SCyT UNS, and PIP-CONICET 0341. A. Juan, G. Brizuela, and S. Simonetti are members of CONICET.