1. Introduction
Aluminum and aluminum alloys are commonly used in the fields of national defense industry, aviation industry, automobile manufacturing, electronics and daily necessities manufacturing with their excellent properties [
1,
2]. The commonly used method of industrial aluminum production is Hall–Héroult method, which has the problems of high energy consumption, serious corrosion of equipment, environmental pollution etc. [
3,
4]. It is urgent to seek a low-temperature green production method of aluminum, and it is also a research hotspot in this field. In recent years, significant progress has been made in low temperature aluminum electrolysis. The main work focuses on the replacement of cryolite molten salt system by molten salt system with lower melting point such as chloride or fluoride. These methods can reduce the electrolysis temperature to a certain extent, but the reduction is not enough [
2,
3,
4]. The appearance of ionic liquids provides a possibility for the development of low-temperature aluminum electrolysis technology.
Ionic liquid is a kind of ionic system which is composed of some inorganic or organic anions and organic cations and is liquid at or near room temperature [
5,
6,
7,
8]. It has the advantages of non-volatile, wide electrochemical window, wide liquid temperature range, good thermal stability, conductivity and variable combination of properties, and is widely and successfully used in various fields such as electrochemistry (electrodeposition, battery), catalysis, organic synthesis, extraction and separation [
5,
6,
7,
8,
9,
10,
11,
12,
13,
14,
15,
16]. In recent years, ionic liquids, as a new type of green electrolyte, have been used for aluminum electrodeposition and electrolysis refining at or near room temperature, and some progress has been made [
9,
10,
11,
12,
13,
14]. However, there are still some problems to be solved, such as high electrolyte viscosity and dendrite phenomenon in the deposition layer [
11,
12,
13,
14]. The addition of additives is one of the common methods to address these problems [
15,
16].
The effects of additives such as chloride, acetonitrile, aromatic compounds, carbonyl small molecules, haloalkanes, nicotinic acid and their derivatives on electrodeposition of aluminum in ionic liquids have been studied, and some progress has been made [
17,
18,
19,
20,
21,
22,
23,
24,
25,
26,
27,
28,
29,
30]. Leng et al. [
24] studied the effect of methyl carbonate on electrodeposition of aluminum in [EMIM]Cl/AlCl
3 system by cyclic voltammetry. Endres et al. [
27,
28] found that Nano Al could be prepared by adding nicotinic acid into [EMIM]Cl/AlCl
3. The effect of aromatic hydrocarbons on electrodeposited aluminum in [Bmim]Cl/AlCl
3 was by cyclic voltammetry [
29]. It was found that the morphology of Al coating became smoother, more uniform and denser with the increase of additive concentration. It was found that adding 45 mmol/L of methyl carbonate can significantly improve the brightness of Al and obtain a fine, uniform and bright Al deposition layer. It can be seen that additives can reduce dendrite and energy consumption in the process of aluminum electrolytic deposition, and change the size of the aluminum film from micron scale to nano-meter scale [
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26,
27,
28,
29,
30]. However, the current researches cannot reveal the detailed information of the mechanism of the additives, such as the bonding, electron transfer, interaction between the molecules and the metal surface, and their influence on the quality of the aluminum deposition layer are not clear.
In recent years, with the development of computational techniques and methods, First Principles method has become an effective means to study the structure, properties and reaction mechanism of compounds and materials. By calculating the frontier orbital distribution, Fukui index, charge distribution of heavy atoms and other quantum chemical parameters, people can study the reaction activity, charge transfer, structure-activity relationship and so on. The microscopic interaction and reaction mechanism between molecules can be deeply revealed from these quantum chemical parameters, which provide a useful way to understand adsorption behavior between molecules and interfaces at the atomic and molecular levels. At present, it has been successfully applied to the study of various metal surface properties [
31,
32,
33]. Valencia et al. [
31,
32] studied the adsorption mechanism of [EMIM]BF
4 on the Li (100) surface with density functional theory (DFT), and found that adding a small amount of [EMIM]BF
4 to Li ion battery can prevent dendrite formation and improve battery performance. Clarke et al. [
33] studied the effect of [EtNH
3]BF
4 on the Li-ion battery system by DFT and molecular dynamics simulation (MD). It is found that the Li-F bond between BF
4 anion and lithium surface forms Li
xBF
4 clusters, which makes the interface between [EtNH
3]BF
4 and the electrode more stable and reduces the dendrite formation. These works provide a new way to explore the interaction mechanism between additives and metal surface in ionic liquids.
Aromatic compounds, as a kind of ring compounds with delocalized bonds, are widely used as additives in electrolysis and electroplating because of their stable structure and difficult decomposition. Pyridine derivatives have been widely used because of their strong polarization and leveling ability in solution and ionic liquids [
9,
10,
11,
12]. It is found that Nicotinamide, Pyridine-2-formamide and Pyridine-4-formamide can be used as effective leveling agents for electrodeposition of aluminum in [Bmim]Cl/AlCl
3 [
34]. The results showed that the order of the deposition was Nicotinamide > Pyridine-2-formamide > Pyridine-4-formamide. However, the mechanism is still unclear, especially the interaction between additives and electrode, charge transfer and adsorption behavior and characteristics. Due to its stable chemical properties (not easy to be oxidized, corroded, etc.) and low hydrogen evolution overpotential, platinum has very important applications in electrochemistry, such as auxiliary electrode, modified electrode substrate, electrocatalysis, etc. Platinum electrode is commonly used in the measurement of electrochemical properties of aqueous solution, molten salt, organic electrolyte and ionic liquid, such as cyclic voltammetry, polarization curve and electrochemical impedance. Therefore, it is very important to study the adsorption of molecules on the surface of Pt electrode. At the same time, platinum has good catalytic performance, which is used as a catalyst in fuel cell, petrochemical industry, automobile exhaust purification and other fields. In this paper, the interaction mechanism of Nicotinamide, Pyridine-2-formamide and Pyridine-4-formamide with Pt (111) surface was studied by First Principles method. The reactivity and adsorption behavior of the three additives on the Pt (111) surface were analyzed, and the interaction mechanism of Pyridine amide derivatives on the Pt surface was revealed, to provide a theoretical reference and guidance for additives’ design to electrodeposition of aluminum in ionic liquids.
2. Materials and Methods
The stable configuration, global and local reactivity parameters of Nicotinamide, Pyridine-2-formamide and Pyridine-4-formamide were calculated by commonly used B3LYP hybrid density functional in the Gaussian 09 software package [
35]. The structure was fully optimized, and the related properties were calculated at the 6–311++G (d, p) basis set. After optimization, the vibration frequency analysis was carried out to confirm the stable geometry without imaginary frequency. The results can be compared with the mp2 results. Then the physical quantity, frontier orbital distribution and Fukui index of the stable geometry are calculated and plotted by Multiwfn [
36].
The Vienna ab initio simulation package (VASP) code has been used for the periodic surface adsorption calculations [
37,
38]. A cutoff energy of 400 eV was used. The electron-ion interaction was described with the projector augmented wave (PAW) method. Exchange and correlation effects are treated in the generalized gradient approximation (GGA) in the formulation of PBE [
39]. The most popular method of DFT-D3 with the BJ-damping method proposed by Stefan Grimme, was used to consider the VDW for present system [
40]. The lattice parameters obtained by optimizing the Pt cell structure were consistent with the previous experimental value of 3.92 Å [
41]. A three-dimensional periodic structure is used to replace the infinite Pt surface. Pt (111) is selected as the research object because it is the most stable low Miller index surface of Pt surface, and used to study the absorption behavior of pyridine and other aromatic chemical compounds on Pt surface [
41,
42]. Four-layer-thick slabs (16 atoms in each layer) were used, and the additive molecule was chemisorption on the upper face only. A vacuum layer of 20 Å is constructed on the upper layer of each cell to ensure that there is no interaction between periodic structures. In the geometry optimizations, the coordinates of the molecule and two uppermost metallic layers were relaxed. Brillouin-zone integrations have been performed using a 4 × 4 × 1 Monkhorst–Pack grid and a Methfessel–Paxton smearing of 0.25 eV. The electron convergence criterion is set to 0.00001 eV by default, and the ion relaxation is set to 0.050 eV/Å. In the initial configuration, the additive molecule is placed on the Pt (111) in parallel with the previous study of imidazole system [
42], and the adsorption energy (
Ead) is defined as follows,
Among them, Ead is adsorption energy of a molecule adsorbs on the Pt (111) surface, Etotal is the total energy of the system when the additive molecule adsorbs on the surface, and Esurface and Emolecule represent the total energy of the Pt surface without adsorbed molecule and the total energy of the additive molecule respectively. Depending on the above definition, Ead of stable adsorption is negative. The density of states (DOS) and partial density of states (PDOS) was analyzed to obtain the electronic structure and properties of the relaxation structure before and after absorption.
3. Results and Discussion
3.1. Molecular Frontier Orbital Distributions
The frontier orbitals include the highest occupied orbitals (HOMO) and the lowest unoccupied orbitals (LUMO). The highest occupied orbit is the outermost orbit occupied by electrons in a molecule, which indicates the ability of the molecule to provide electrons. The lowest unoccupied orbital is the innermost orbital not occupied by electrons in a molecule, which shows the ability of a molecule to obtain electrons. Depending on the frontier orbital theory [
35,
43], the interaction between reactants in a chemical reaction only occurs between the frontier orbitals of molecules, and HOMO and LUMO determine the reaction activity of molecules. Therefore, the frontier orbitals of molecules play a major role in the chemical reactions. The reactivity of molecules can be calculated by analyzing the distribution of HOMO and LUMO. Optimized stable geometry structures of three additives are shown in
Figure 1. The HOMO and LUMO isosurface maps of these three molecules are illustrated in
Figure 2. It can be seen from
Figure 2 that HOMO of Nicotinamide and Pyridine-2-carboxamide is distributed on the pyridine ring except C3 atoms, and Pyridine-4-carboxamide is distributed on the whole pyridine ring. They are also allocated on the C, N and O atoms of substituents, and these atoms have electrophilic activity. For LUMO, Nicotinamide and Pyridine-2-carboxamide are distributed on the pyridine ring and C, N and O atoms of substituents, while Pyridine-2-carboxamide is distributed on the pyridine ring and C, N and O atoms of substituents except C13 atom. The HOMO and LUMO of the three molecules are delocalized on the pyridine ring, N atoms of amino and O atoms of carbonyl. It indicates that the pyridine rings, N atoms of amino and O atoms of carbonyl are the active regions of the molecules. Due to the atoms in the pyridine ring, N atoms and O oxygen of the three derivative molecules are basically coplanar, and it is easy to form conjugated large π bonds, which enhances the overall activity of the three molecules and makes the molecules parallel to the surface in a flat way. In addition, HOMO and LUMO cover almost all the heavy atoms in the molecule.
3.2. Molecular Local Reactivity
3.2.1. Charge Distribution
Through the analysis of the molecular frontier orbitals of the three additives, the active regions of the molecules were obtained. Taking into account the charge distribution and Fukui index of each atom in the molecules, the local reactivity of the molecules was also analyzed [
44]. The Mulliken charge distributions of atoms in the three derivatives are listed in
Table 1. It can be seen from
Table 1 that C1, C5, C10, N12 and O14 atoms in Nicotinamide molecule have more negative charges, which are consistent with the contribution of HOMO in
Figure 1a. It can be concluded that these atoms are electrophilic active centers of Nicotinamide molecule. C4 and N13 with a small amount of negative charge have a certain electrophilic activity. At the same time, these atoms contribute a lot to LUMO, which indicates that these atoms also have some nucleophilic activity. C3, N10 and O11 atoms in Pyridine-2-carboxamide molecule have more negative charges, which are consistent with the contribution of HOMO in
Figure 1b. It can be concluded that these atoms are electrophilic active centers of Pyridine-4-carboxamide molecule. C1, C8, C13 and N15 atoms have a small amount of negative charge and have a certain electrophilic activity. At the same time, these atoms contribute a lot to LUMO, which indicates that these atoms also have some nucleophilic activity. C3, C11, N8, O9 and C6 atoms in Pyridine-4-carboxamide have more negative charges, which is consistent with the contribution of HOMO in
Figure 1c. It can conclude that these atoms are electrophilic active centers of Pyridine-4-carboxamide. N13 with a small amount of negative charge has a certain electrophilic activity. At the same time, these atoms contribute a lot to LUMO, which indicates that these atoms also have some nucleophilic activity.
3.2.2. Fukui Index Analysis
Fukui index can accurately reflect the activity of each atom in the additive molecule, and can locate the nucleophilic or electrophilic reaction center of the molecule, which is an effective method to determine the molecular reaction center and its strength [
43].
f+i(
r) and
f−i(
r) denote the nucleophilic Fukui index and the electrophilic Fukui index respectively, indicating the strength of the ability of an atom
i to accept or donate electrons. The larger the value is, the stronger the ability to accept or donate electrons is [
44]. The
f+i(
r) and
f−i(
r) of all atoms in the three additives are provided in
Figure 3. It can be seen from
Figure 3 that
f+i(
r) values of C5, O14, C10, C2, C3 and N12 atoms in Nicotinamide are larger, indicating that these atoms have nucleophilic activity. The
f−i(
r) values of O14, N15, N12 and C2 atoms are bigger, indicating that these atoms have electrophilic activity. O14, N12 and C2 atoms are electrophilic active centers and nucleophilic active centers, so O14, N12 and C2 atoms should be preferentially adsorbed on the Pt (111) surface. The larger
f+i(
r) values of O11, C1, N15, C4, C3 and N10 atoms in Pyridine-2-carboxamide indicate that these atoms have nucleophilic activity. The
f−i(
r) values of O11, N15, C1, N10 and C13 atoms are larger, indicating that these atoms have electrophilic activities. O11, N15, C1 and N10 atoms are both electrophilic active centers and nucleophile active centers, so O11, N15, C1 and N10 atoms should be preferentially adsorbed on the Pt (111) surface. The larger
f+i(
r) values of O9, C1, N13, C6 and N8 atoms in Pyridine-4-carboxamide indicate that these atoms have nucleophile activity. The
f−i(
r) values of O9, N8, N13, C1, C11 and C14 atoms are larger, indicating that these atoms have electrophilic activities. O9, N8, N13, C1, C11 and C14 atoms are both electrophilic and nucleophile active centers, so O9, N8, N13 and C1 atoms should be preferentially adsorbed on the Pt (111) surface. The results are in agreement with
Section 3.2.1.
3.3. Adsorption of Three Derivative Molecules on the Pt (111) Surface
According to previous studies and literature [
41,
42], the bridge adsorption site is the most stable adsorption site for pyridine molecules on the Pt surface. Therefore, we consider that three molecules are adsorbed in the bridge adsorption site on the Pt (111) surface in a parallel manner.
Figure 4 shows the top view (a–c) and side view (d–f) of the optimized structure of Nicotinamide, Pyridine-2-formamide and Pyridine-4-formyl adsorbed on the Pt (111) surface, respectively. In the initial structure, the carbon rings are parallel to the surface at an average distance of about 2.1 Å. After optimization, the average distance between the carbon ring and the metal surface is about 2.2 Å. As can be seen from
Figure 4d–f, after optimization, it is found that the length of the C-C bond in the ring group is extended, while the C-H bond is bent upward, indicating that the hybrid state of carbon changes from planar type to non-planar type, thus forming a strong chemisorption bond.
It can be seen from
Figure 4 that the three additives are adsorbed on the Pt surface in a horizontal form, and the N atom of the amino and one Pt atom forms the N-Pt bond. However, there are some differences in the type and quantity of the bond between the pyridine ring and the Pt surface in the three additives.
Figure 4a shows the stable configuration of Nicotinamide adsorbed on the Pt (111) surface. Each C atom in the pyridine ring is bonded with Pt surface atom, and the Pt atom is stretched upward to form two π bonds and two σ bonds. C2-C3 and C4-C6 form π-bonds with one Pt atom respectively. C1 atom forms an σ-bond with a Pt atom, and N1 atom forms an N-Pt bond with a Pt atom.
Figure 4b shows the stable configuration of Pyridine-2-formamide adsorbed on the Pt (111) surface. Each C atom on the ring group is bonded with Pt surface, forming one π-bond and three σ-bonds. The N atom of the ring is not bonded with Pt atom. C3-C5 forms a π-bond with a Pt atom. C1, C2, and C6 atoms form σ-bonds with one Pt atom, respectively.
Figure 4c shows the stable configuration of Pyridine-4-formyl adsorbed on the Pt (111) surface. Each C atom on the ring group forms adsorption bonds with Pt surface, forming two π-bonds and two σ-bonds. C2-C3 is a π-bond formed by two C atoms and one Pt atom, N1-C6 is a π-bond formed by two Pt atoms and one Pt atom. C1 and C4 atoms form σ-bonds with one Pt atom respectively.
Table 2 displays the calculated adsorption energy of the additive on the Pt (111) and the distance between the corresponding atoms after adsorption. It can be seen from
Table 2 that the adsorption energy varies with the substituent position on the pyridine ring. The absolute value of the adsorption energy of meta-substituted Nicotinamide is the largest, which is −2.04 eV, indicating that it has the strongest interaction with Pt (111). The C and N atoms of Nicotinamide are bonded to the nearest Pt atom. C-Pt bond growth is about 2.2 Å, and two N-Pt bond lengths are 2.16 Å and 2.22 Å, respectively, which close to the experimental value of 2.09 Å. The adsorption energy of ortho-substituted Pyridine-2-carboxamide is −1.79 eV. The C and N atoms of the Pyridine-2-carboxamide form adsorption bonds with the nearest Pt atom. Two C atoms are adsorbed on one Pt atom, and the C-Pt bond growth is about 2.2 Å. The N atom of the substituent group is bonded with Pt atom, and the bond length is 2.23 Å. The adsorption energy of para-substituted Pyridine-4-carboxamide is −1.78 eV, and its adsorption bonding is similar to that of Nicotinamide, except that the N atom of the ring group and the adjacent C atom are adsorbed on a Pt atom. The adsorption of three additives on the Pt (111) surface is chemical adsorption [
38]. The adsorption order of the three additives was Nicotinamide > Pyridine-2-carboxamide > Pyridine-4-carboxamide.
Three pyridine derivatives are adsorbed on the Pt (111) surface by chemisorption. They are strongly adsorbed on the Pt (111) surface to form a compact organic adsorption layer, which has a certain blocking effect on the passage of electric current, hinders the reduction of metal ions, significantly reduces the discharge rate of coordination ions, increases the overvoltage of electrode reaction, and slows down the electrode reaction rate, so that the electric current is reduced As a result, the growth rate at different positions on the surface tends to be the same, which prevents the appearance of coarse-grained and dendritic coatings formed by unidirectional growth. It can significantly change the quality of cathodic electro-deposition and obtain a flat and smooth deposition layer. Nicotinamide has the strongest adsorption, so it has the best leveling effect, which is consistent with the experimental observation.
3.4. Electronic Structure Analysis
In order to further analyze the adsorption mechanism of three pyridine derivatives on the Pt (111), the density of states (DOS) and partial density of the states (PDOS) diagrams before and after adsorption was calculated and shown in
Figure 5,
Figure 6 and
Figure 7. It can be seen from
Figure 5,
Figure 6 and
Figure 7 that the contribution of the d electron of the atom on the Pt (111) plane to the total density of states is the largest, followed by the s and p orbital, but the contribution is relatively small. The p orbital electron of the additive molecule atom has the greatest contribution to the total density of states.
DOS of Nicotinamide/Pt (111) is shown in
Figure 5a, and PDOS is shown in
Figure 5b–d. After adsorption of Nicotinamide, the total DOS of Pt (111) changed to a certain extent, which showed that DOS increased in the range of −7 eV to −9 eV. Depending on the PDOS of characteristic atoms in
Figure 5b–d, the increase of DOS is mainly composed of 2p orbitals of C and N atoms and 4s and 3d orbitals of Pt atoms. In
Figure 5b, the 4s orbital of Pt and the 2p orbital of C obviously overlap at −6 eV, and the 3d orbital of Pt and the 2p orbital of C obviously overlap at −7 eV. The results show that the interaction between two atoms is the strongest, and the π-bond formed by C-C atoms bonding to a Pt atom plays a major role in the adsorption process. In the adsorption structure, the 3d and 4s orbitals of Pt and the 2p orbitals of C are relatively active, and the hybridization between these orbitals occurs in the interaction process. In
Figure 5c, there is a slight overlap between the 3d orbital of the Pt atom and the 2p orbital of the N atom near −7.2 eV. Therefore, it shows a weak interaction between Pt and N atoms. The results show that the C atom and N atom bond to one Pt atom to form σ-bond, which plays a secondary role in the adsorption process.
DOS of adsorbed Pyridine-2-carboxamide/Pt (111) is shown in
Figure 6a, and PDOS is shown in
Figure 6b–d. After adsorption, the total DOS of Pt (111) increased from −9 eV to −6.9 eV, similar to that of Nicotinamide. The change of PDOS of the C atom is similar to that of Nicotinamide system. There is a slight overlap between the 3d and 4s orbitals of Pt atom and the 2p orbitals of the N atom of amino at −8.3 eV to −7 eV. Due to the N atom of the ring not bonding with Pt atom, the π-bond formed by two C-C atoms and one Pt atom plays the main role in the adsorption process of Pyridine-2-carboxamide system, and the interaction between N atom of amino and Pt atom plays the secondary role.
DOS of Pyridine-4-carboxamide/Pt (111) is shown in
Figure 7a, and PDOS is shown in
Figure 7b–d. The total DOS of Pt (111) increased in the range of −9 eV to −7 eV after Pyridine-4-carboxamide adsorption. In
Figure 7a, the 3d, 4s and 2p orbitals of Pt and C overlap at −7.5 eV, and the 4s and 2p orbitals of Pt and C overlap at −6.3 eV. In
Figure 7c, the peak value of 2p orbital −7 eV of the N atom is the largest, which indicates that the interaction between the N atom and the metal surface is the largest. In this adsorption structure, the p orbital of N1-C6 on the ring forms π bond with the d orbital of the Pt atom, and the C atom forms σ bond with one Pt atom, which plays a secondary role.
Therefore, the π-bond between C and N atoms and Pt atom plays the main role, the σ-bond between C and N atoms and Pt atom plays the secondary role, and finally, the N atom of formamide substituent plays the role. Moreover, due to the substituents at different positions of the three additives, the bonding orbitals are different, which indicates that the interaction between the additives and Pt atoms is mainly caused by the hybridization of 2p orbitals of C and N atoms, and 4s and 3d orbitals of Pt atoms, forming a very stable internal orbital complex. It is very difficult for the ligands to separate from these complexes and form the activated complexes for the electrode reduction reaction. Therefore, the activation energy of cathodic reduction is very high and the rate of electrode reaction is very slow. This is why the additive molecules can form stable complexes with Pt and adsorb on the surface of Pt.
According to the above analysis, the additive has high surface activity and is strongly adsorbed on the metal surface. One or two Pt atoms are greatly attracted by the additive molecules, which is stabilized by the formation of C-Pt and N-Pt bonds between C, N atoms and Pt (111) surface, resulting in the electron transfer of the attracted Pt atoms to the additive molecules, and part of the electrons are also obtained on the metal surface. As a result, the additive molecules are strongly adsorbed on the Pt surface to form a film. With the increase of cathode current density, the additive may be preferentially adsorbed on the high current density part of the electrode, which may reduce the reaction speed of electrodeposition of metal ions, makes the precipitated grains finer, and makes the current distribution at the low current density part of the electrode tend to be uniform, thus inhibiting the formation of dendrites, resulting in a smooth aluminum deposition layer.