Sodium nitrite as a corrosion inhibitor of copper in simulated cooling water

The corrosion inhibition behavior of sodium nitrite (NaNO2) towards pure copper (99.95%) in simulated cooling water (SCW) was investigated by means of electrochemical impedance spectroscopy (EIS) and dynamic electrochemical impedance spectroscopy (DEIS). NaNO2 interferes with metal dissolution and reduce the corrosion rate through the formation or maintenance of inhibitive film on the metal surface. Surface morphologies illustrated that the surface homogeneity increased on adding sodium nitrite. Sodium nitrite’s adsorption on copper surface followed the modified form of Langmuir, Freundlich and Frumkin isotherms. Physiosorption mode was involved in the corrosion protection. Electrochemical results revealed an corrosion resistance of copper increases on increasing the inhibitor concentration. The DEIS results indicated that copper corrosion mechanism could be hindered by 50% even after interval of 24 h by optimum concentration of sodium nitrite. The maximum inhibition was achieved with 2000 ppm of NaNO2. With this concentration, inhibition efficiency of up to 61.8% was achievable.

www.nature.com/scientificreports/ water lies in a neutral range. Although the pH of the prepared SCW was recorded as 7.07 initially, on dissolving the maximum concentration of inhibitor, it increased to 7.90. Copper is fairly inert metal, however aggressive ions in SCW like chlorides, sulfates, carbonates and bicarbonates can cause copper to degrade and corrode as shown in the following equations.
The partially soluble cuprous chloride in dilute NaCl reacts to form cuprite or cuprous oxide (Cu 2 O). This cuprous oxide, in the presence of dissolved salts in water, oxidized to cupric hydroxide (Cu(OH) 2 ). The dissolved salts in SCW not only lead to corrosion but also scaling and depositions. Normally copper in cooling systems tends to form a protective passive layer of Cu 2 O which later on evolves to CuO/Cu(OH) 2 which can be confirmed by checking the potential and the pH of electrochemical process in the E-pH diagram of aqueous corrosion of copper 25,26 . Higher temperatures lead to conversion of Cu (OH) 2 to CuO completely, but as the analysis is done at room temperature, a predominant presence of Cu(OH) 2 is expected in the passive layer.
Chloride is particularly corrosive to copper, even at basic values of pH. Chloride and sulfates are well known to antagonistically accelerate the degradation of copper. The small aggressive anions migrate through the film to the areas of high positive charge density. Substitution of monovalent chloride species for divalent charged oxygen species might occur leading to the release of Cu 2+ ions in the bulk solution (corrosion). Calcium and magnesium along with the bicarbonate ions present in SCW protect the copper by forming a deposited film, however localized corrosion may still occur beneath these deposited layers. Hence pitting corrosion is expected to occur profusely on the copper metal in such scenarios.
In order to tackle copper corrosion, corrosion inhibitors are the most convenient and economic options. The role of inhibitor to suppress the anodic or cathodic reactions or both of them, determine the type of corrosion inhibitor it might be. If the addition of corrosion inhibitor to corrosive environment reduces anodic dissolution of copper, the corrosion inhibitor provides anodic inhibition. Corrosion inhibitors, including inorganic or organic compounds like azoles, amines and amino acids effectively protect copper and copper alloys. This action is associated with the chelation effect of the functional groups in them or the formation of impermeable barrier between the inhibitor and the copper, preventing its dissolution. Also, presence of vacant d orbitals in copper atom leads to coordinate interactions with heteroatom such as nitrogen, sulfur and oxygen or π-interaction amidst aromatic rings via the heterocyclic section of the inhibitor. Many organic and inorganic substances have been utilized as corrosion inhibitors of copper in acidic, alkaline and neutral media. While the chromates which increase the corrosion rate by increasing the rate of cathodic reactions, have been rejected because of their toxicity, the tetraborates and molybedates have been observed to form instable protective films on copper surface. The benzotriazoles act by formation of a protective monolayer or multilayer film in presence of oxidants or by anodic polarization. For alkaline and neutral mediums, benzotriazoles have given impressively high corrosion protection capacities yet their non-biodegradability poses a danger of flow of non-biodegradable substances in high quantities in general waterways 27 . Sodium is a neutral species whose effect is non-discernable in cooling water systems' corrosion. In this research article sodium nitrite has been utilized to check its effect on corrosion of copper in SCW. It is hypothesized and proved that nitrite will suppress the anodic dissolution of copper 28 . The authors in the current research article have tried to explain the exact function of sodium nitrite in the protection of copper in SCW. Sodium nitrite is nontoxic and ecofriendly compound used as a food preservative. It is a well-known fact that the corrosion inhibitors with heteroatoms like 'N' , tend to adsorb on the surface to form inhibitive film as a barrier between electrolytic solution and metal to stop the corrosion reaction from occurring 29 . The inhibition mechanism of nitrite ion described by previous researchers shows that it starts by conversion of nitrite ions to nitrous oxide and oxygen ions as follows.
These O 2− ions oxidize cuprous ions to cupric. In an alkaline solution, stable passive film of CuO/Cu (OH) 2 forms on the metal surface 30 .

Results
EIS measurements. Impedance spectroscopy is a common practice carried out to estimate the resistance between the interface of a metal and a corrosive medium. Figure 1a shows the EIS data of copper samples subjected to the various concentrations of NaNO 2 in SCW. The Nyquist plots indicate the capacitive behavior of the spectra which means that the corrosion process of the copper electrode includes charge transfer resistance. The diameter of capacitive semicircle increased with addition of NaNO 2 , implying that the corrosion protective behavior of copper was improved on increasing the inhibitor concentration. The Bode phase plots showed a single constant in the impedance spectrum which was again generally associated with charge transfer process (Fig. 1b). The EIS graphs indicated that the copper surface was covered by a thin film, a passive layer which prevented the attack of the corrosive ions present in SCW. The phase angle diagrams also displayed that maximum phase angle was increased each time the concentration of NaNO 2 increased in the simulated cooling water. This proves again that addition of NaNO 2 was beneficial to the growth of inhibitive barrier.
This might be indicative of the fact that the passive layer of oxides or the thin film of Cu 2+ /inhibitor complex covering the copper surface thickens or improve with increasing inhibitor concentration. The same can be understood by the values of Y o associated with growth of double layer capacitance. Although the values of R ct increase on adding increasing its concentration, an increase in the values of R s is also noted. Usually film www.nature.com/scientificreports/ forming inhibitors get adsorbed to the metal and R s tends to decrease. In case of NaNO 2 , it is generally considered anodic inhibitor which helps in forming the thin protective layer of oxide on the copper surface 30 . This implies that due to this oxide barrier on surface, the ion transfer between copper and bulk solution decreases each time the inhibitor concentration is increased that results in increasing R s . The impedance data is enlisted in Table 1.
The impedance values were accurately fitted to R(QR) circuit and the χ 2 values of 10 −3 were obtained to ensure the correctness of the this fit (Fig. 1c). R s is the non-compensated solution resistance, R ct is the charge transfer resistance and Q is the constant phase element, CPE used in place of 'C dl ' double layer capacitance to account for the non-ideal behavior of the working electrode 31 . 'n' is the quantification of the capacitor's heterogeneity whose  www.nature.com/scientificreports/ value may be in the range of 0 < n < 1. Table 1 displays improvement in surface homogeneity with increasing inhibitor concentration. Table 1 distinctly specifies that the values of R ct increased with increasing NaNO 2 concentration. R ct values obtained for uninhibited solution are resultant of primarily adsorbed water molecules and other ions on the metal/electrolyte interface, which were later displaced by inhibitor species. A potential of 0.2 V/Ag/AgCl corresponding to pH value of 7.9 might be attributed to the initial formation of Cu 2 O, which evolved to CuO/Cu   www.nature.com/scientificreports/ (OH) 2 at later stages in uninhibited electrolytic solution. However it can also change to Cu 2+ /Inhibitor in the inhibited solution 26 . At the maximum concentration, NaNO 2 is 61.8% efficient inhibitor according to EIS.
DEIS measurements. DEIS analysis at different time intervals to get an insight for tracking the changes in the adsorbed inhibitor film of nitrite on the copper surface in SCW was carried out up to 24 h of sample immersion. Figure 2a represents the DEIS spectra of 24 h for copper specimens tested at room temperature after attaining a stable open circuit potential (OCP). In order to gain better understanding of the dissolution pattern of formation of Cu 2+ ions on the copper surface and prevention of this corrosive reaction by sodium nitrite in SCW the DEIS experiments were conducted for all the environments at different time ranges of 2 h , 6 h, 12 h and 24 h, respectively. It was observed that R ct improved with time indicating better performance of inhibitor with increase in its concentration. The equivalent circuit used in the DEIS analysis was similar to EIS analysis (Fig. 1c). As usual, the best fit was determined based on the low values of χ 2 in the range of 10 −3 . The DEIS parameters thus obtained have been given in Table 2. The influence of the addition of the maximum inhibitive concentration of sodium nitrite to SCW can be clearly understood by looking at the charge transfer versus time plot in Fig. 2b.
The values of R ct after 2 h interval with DEIS coincide with EIS result of maximum inhibitor concentration which indicated the consistency of the tests performed. In the case of the uninhibited solution it was observed that R ct dips for initial 0 to 14,400 s which could be associated with initial attack of the aggressive ions in the SCW. After 14,400 s a very gradual constant increase in the values of R ct was observed which was resistant to the aggressive solution. In the copper specimen inhibited by 2000 ppm of sodium nitrite in SCW, there was no dip in R ct values which is due to the inhibition offered by nitrite ions and prevention of dissolution of copper as well as prevention of chloride ions in destroying the passive layer of oxide which formed on the copper surface in aqueous solutions. However, the charge transfer value for the first spectra obtained, immediately after immersion in inhibited solution was lower than that of the blank SCW solution.
The previous researches conducted on the behavior of nitrite ion towards copper and its oxides in different media have concluded that the growth of resistance of initially formed Cu 2 O on the copper surface is hindered by the nitrite ion, however after some time when the Cu 2 O evolves into Cu 2 O + CuO/Cu(OH) 2 , the nitrite ions starts assisting in growth of resistance offered by oxide film by reacting with Cu(OH) 2 25,26 . In the absence of nitrite ions in the uninhibited solution the growing film of Cu 2 O initially formed on copper specimen exhibited some resistance, along with attack of aggressive ions which displaced oxygen and lead to anodic dissolution. However interestingly, when inhibitor containing nitrite ion was added to the blank SCW solution, it hindered the growth of cuprous oxide layer on the surface which led to ease of attack on SCW ions on the surface and corrosion. Hence, the R ct began to evolve from a lower value in the first spectra, than that observed in blank solution ( Table 2). However exactly after 14,400 s the development of copper (II) hydroxide along with cupric oxide favoured the corrosion protection by formation of copper nitrite on the surface. This led to the elevation of R ct values higher than in blank solution after 2 h and a steady inhibition efficiency in the range of 50-59% was observed over a period of 24 h. At 2 h interval 54% IE was achieved which might be due to acceleration of the growth of passive oxide layer by nitrite ion. At 6 h IE attains a maximum of 59.5%. At 12 h the IE depleted slightly to 58% and further reduced to 50.5% at 24 h interval. Generally it was expected that nitrite ion would keep assisting the growth of passive layer and IE would maintain an increasing trend till 24 h. This lowering of IE suggested that nitrite ion was itself involved in some reaction on the copper surface. If the process of physiosorption is involved in such association, after a set period of time the decline in IE could be explained by the mode of adsorption involved which will be discussed in the following section. Note worthily, after 14,400 s the nitrite ion interferes more with the anodic dissolution of copper and higher R ct is obtained compared to the blank SCW solution. The results of this analysis suggested that sodium nitrite was effective in prevention of copper corrosion in SCW.
Adsorption isotherms. The isotherms involved in depicting the process of adsorption of nitrite ion on the copper surface provided better understanding of the corrosion inhibition in the current study. The adsorption process tends to obey modified Langmuir adsorption isotherm or El-Awady isotherm, Freundlich adsorption isotherm and Frumkin adsorption isotherm based on the values of correlation "R 2 " which should be close to unity. The adsorption parameters thus obtained have been enlisted in Table 3.
The experimental data were best suited to modified Langmuir isotherm which is generally called El-Awady isotherm (Fig. 3a) 32 . El-Awady isotherm is described by www.nature.com/scientificreports/ where y is number of inhibitor molecules accommodated in unit active site, ' θ ' the surface coverage, 'C' the concentration and 'K ads ' the constant related to the adsorptive equilibrium constant.
When this 1/y is below 1 it suggests multilayer adsorption. Currently, the values of 1/y obtained were more than one (1/y = 1.4) which meant that each inhibiting molecule was attached to more than one active site on the (4) log www.nature.com/scientificreports/ copper surface. This might be because of the presence of two hydroxide or two nitrite ions in association with the anodic site. The adsorptive equilibrium constant 'K ads ' (obtained from the intercept of the graphs of adsorption isotherms) helped in the derivation of adsorption free energy 'ΔG ads ' as per following equation: where the water molecules concentration is 1000 gm/L, 'R' is the universal gas constant (R = 8.314 J K −1 mol −1 ) and T delineate the absolute temperature (K) which is room temperature in the current analysis. Freundlich isotherm (Fig. 3b) suggests the adsorptions and interactions according to the following equations where 1/n = 1.15 is a constant depending on the nature of the adsorbed molecule, K ads (3.2 × 10 3 ppm −1 )is the adsorption-desorption constant which displays interaction strength of the adsorbed layer. In fact, when its value is near zero, it means non homogenous surface; whereas, a value below unity means a strong adsorption (chemisorption), similarly, above unity it is indicative of moderate adsorption (comprehensive/mixed). The value of 1/n was slightly more than unity which suggested moderate physiosorption on the copper surface The correlation coefficient for Freundlich isotherm was 0.98 33 . The Frumkin adsorption isotherm (Fig. 3c) assumes copper surface as heterogeneous and interaction between the adsorbed molecules as per following equation: α (0.37) is related to the molecular interaction in adsorbed layer.where α is the lateral interaction parameter that describes the interaction in the adsorbed layer. The molecular interaction parameter can be assumed positive or negative: when α < 0, there is repulsion in the adsorbed layer, and otherwise attraction 34 .
Negative values of ΔG ads meant a spontaneous adsorption. When the values of (ΔG ads ) are up to − 20 kJ mol −1 , they suggest electrostatic interactions between the charged molecules and the charged metal (physisorption) and when these values are around − 40 kJ mol −1 , they suggest chemisorption which is as a result of sharing or transfer of electrons from inhibitor molecules to the metal surface to form a coordinate bond (chemisorption) 35 . The values of ΔG ads listed in Table 3 were all negative which suggested a spontaneous adsorption of inhibitor on copper surface. The range of the values suggested the predominant role of physiosorption, which implies a minor possibility of formation of a complex on the metal substrate. However El-Awady isotherm suggested the role of comprehensive type of adsorption which might occur in the presence of two heteroatoms N and O directly involved in the protection of copper surface. Summarizing, from the adsorption analysis it was suggested that more than one species of inhibitor was associated with each active site on the copper surface and the adsorbed species were interacting with each other. After 2 h CuO/Cu(OH) 2 was formed on the surface of copper. There is also a possibility of it undergoing a displacement reaction for nitrite ion to get adsorbed on the copper surface.
Surface analysis. Surface analysis led to important inferences about the adsorption of NaNO 2 on copper surface in SCW. The improvement of surface properties, quantitative and qualitative judgment of the molecules adsorbed on the copper specimens could only be understood through surface analysis. One of the important objectives of the study of surface composition of inhibited specimen was to know whether the adsorption is solely dependent on oxygen or nitrogen is also associated with the copper surface.
AFM. To judge the efficiency and role of NaNO 2 in corrosion prevention and maintaining the homogeneity of the copper surface when immersed in SCW, 3D and 2D images obtained by AFM were observed for each specimen. The freshly polished copper specimen had a mean roughness factor R m of 27.4 nm (Fig. 4a). When this polished specimen was immersed in SCW for 2 h, the roughness was degraded to R m = 246.8 nm (Fig. 4b). The abrasion lines apparent in the polished specimen were dissolved in the specimen immersed in SCW.
The corroded specimen in the uninhibited solution, also displayed blue colored regions in 2D view which might be attributed to the pitting caused by chloride ions (Fig. 4a). According to the scale of the images, the blue regions depicted are the areas where the metal was been removed by corrosion. The localized corrosion is very clearly evident in the corroded specimen (Fig. 4b). The roughness of the specimen was less when it was immersed in a solution of SCW inhibited by 2000 ppm of NaNO 2 , apparent reduction of blue regions was observed. The mean roughness improved to R m = 32.9 nm (Fig. 4c). Apparently few blue regions (pits) are seen in the 2D image of inhibited specimen in Fig. 4c, which suggests less pitting and low corrosion. Addition of sodium nitrite clearly improved the surface roughness of the copper specimen was very apparent in the AFM images. Figure 5a,b,c,d,e and f represented the SEM and associated EDX graphs of copper specimen having different morphologies and compositions in different environments. The changes occurring on the copper www.nature.com/scientificreports/ www.nature.com/scientificreports/ surface were apparent from the images obtained by SEM and elemental composition of the surface suggested by EDX. The SEM of polished copper in Fig. 5a displayed an un-corroded, fairly smooth surface with a few abrasion lines, which could be attributed to surface preparation by SiC paper. The relative EDX graph indicated that the surface was devoid of any corrosion products and is pure copper specimen (Fig. 5b). The SEM of the copper sample immersed in SCW in Fig. 5c displayed a considerably corroded and pitted surface which can be associated www.nature.com/scientificreports/ www.nature.com/scientificreports/ with the presence of chloride, sulphur and oxygen signals in the EDX spectra (Fig. 5d). On the surfaces of corroded specimen localized corrosion was more apparent along with some uniform corrosion. The amount of chlorides and sulfide ions on the surface was bearing a low signal because of their low millimolar concentration in SCW. However it is interesting to note that even in such concentrations these aggressive ions have posed considerable local damage to the copper surface. The oxygen signal was more apparent than the chloride signal and even its weight percentage is higher than that of chloride, which indicated the presence of oxides on the corroded surface, particularly cuprous oxide in case of uninhibited specimen. The primary feature noticed in the SEM of the copper specimen inhibited by 2000 ppm NaNO 2 in Fig. 5e was that considerable corrosion and pitting which was previously observed, now highly reduced compared to the uninhibited specimen corroded in SCW.

SEM/EDX.
To get a more clear inference, the EDX values of the elements present on the copper substrate displayed the occurrence of considerable amounts of N, O and trace amounts of S. The amount of Cl − ions was reduced in the inhibited surface according to the EDX graph (Fig. 5f) which suggests that nitrite ions have clearly prevented the association of chloride with the copper surface. The signal of oxygen was a little low in inhibited specimen because sodium nitrite is more favorable to the development of cupric oxide and hydroxides on the copper surface, which could have undergone displacement by inhibitor molecules too and hence a low signal for oxygen was obtained. The EDX spectrum also displayed an improvement in copper percentage in inhibited specimen.
Optical profilometry. The optical profiling was done to know the surface roughness, materials loss and pits statistics on the copper surface to understand the effect of SCW and inhibitor on it.
The 3D image very clearly depicted the red regions on the polished sample (Fig. 6a), the uncorroded red areas diminished in the uninhibited corroded specimen (Fig. 6b), while the inhibited sample showed corroded as well as more protected areas in the 3D image (Fig. 6c). The 2D image and the height profile clearly displayed the occurrence of localized and pitting corrosion in corroded specimen which was immersed in SCW (Fig. 6b). The inhibited specimen displayed an improvement in the depth of valley ratio and average roughness of the whole sample in the presence on 2000 ppm NaNO 2 in SCW (Fig. 6c). The quantitative data is present in the Table 4. Table 5, frontier orbital energies, energy gap, electronegativity, hardness, fraction of electrons transferred, electron donating power, electron accepting power and dipole moment values calculated for nitrite ion at various calculation levels of the theory are given as detailed (Fig. 7). Frontier orbital energies, namely HOMO/high occupied molecular orbital energy and LUMO/low unoccupied molecular orbital energy, are widely considered in the analysis of corrosion inhibition performances of molecules. High values of E HOMO facilitate the adsorption of inhibitor molecules on metal surfaces and therefore enhance the inhibition efficiency. On the other hand, the energy level of the lowest unoccupied molecular orbital represents the ability of the molecule to accept electrons from an electron donor. Energy gap, namely, the difference between HOMO and LUMO orbital energies is an important indicator of the reactivities of molecules. It should be noted that high energy gap values belong to the molecules with low reactivity. This situation can be explained in the light of Maximum Hardness Principle.

Computational analysis. In
According to hard and soft acid-base principle introduced by Pearson, hardness, softness and energy gap are closely related parameters to each other. Pearson defined the hardness as the resistance towards electron  www.nature.com/scientificreports/ cloud polarization or deformation of chemical species. The hard molecules with high energy gap values are not effective against the metal corrosion. On the other hand, soft molecules with low energy gap values are good corrosion inhibitor. In the hard and soft acid base classification of Pearson, NO 2 − ion appears between borderline bases. Here it is apparent that nitrite ion is more effective against the corrosion of copper surface compared to  www.nature.com/scientificreports/ nitrate ion that is a hard base. It could be seen from the data given in Table 5 that hardness values calculated for nitrite ion are quite low in general. According to maximum hardness principle 36 hard molecules are more stable compared to soft ones. Adsorption behavior of nitrite ion on copper surface supported the validity of this electronic structure principle. The tendency of a chemical species to attract electrons is explained with the help of electronegativity concept. The molecules having high electronegativity values are not effective against the metal corrosion. Electronegativity values calculated for studied inhibitor molecule were not that high. Another useful parameter used in corrosion inhibition studies is electrophilicity index. An electronic structure principle known as minimum electrophilicity principle correlates the electrophilicity index with chemical stability. According to this principle the natural direction of a chemical reaction is toward a state of minimum electrophilicity. From here, it can be concluded that the molecules with low electrophilicity value are more stable. In contradiction with this principle, in corrosion studies, it is widely reported that the molecules having high electrophilicity values are not good corrosion inhibitor. Nitrite ion with low electrophilicity value is a good corrosion inhibitor for copper surface. This result shows that minimum electrophilicity principle 37 is not successful in the explanation of corrosion inhibition activities of molecules. According to electron donating and electron accepting power parameters derived by Gazquez and coworkers, good corrosion inhibitors should have high electron donating power values. It is apparent from the related table that electron donating power values calculated for studied ion is quite high. Adsorption energy values obtained molecular dynamics simulation approach shows that the interaction between nitrite ion and copper surface is quite powerful. It can be understood from the data given in Table 6 that nitrite ion interacts more strongly with Cu (111) surface compared to other ions. Theoretically obtained results are in good agreement with experimental observations.

Conclusions
The effect of addition of NaNO 2 to the electrochemical and surface properties of copper in SCW was observed. The following inferences were made: 1. A maximum concentration of NaNO 2 at 2000 ppm prevented 61.8% of corrosion as suggested by EIS. The EIS parameters suggested the formation of a thin barrier film which improved with increase in inhibitor concentration. The first inference which was derived was that nitrite ions have promoted the growth of passive oxide layer on copper surface. 2. The DEIS analysis gave some important insights into role of nitrite ion in inhibition of copper corrosion in SCW. The nitrite ion diminished the growth of initial oxide layer formation of Cu 2 O but at exactly after two hours improved the inhibition of copper when Cu 2 O + Cu(OH) 2 was formed after 14,400 s. The inhibition effect was 50-59% over the period of 24 h. The diminishing of the IE values at 12 and 24 h could be associated to the loss of nitrite ions from the copper surface which were physiosorbed and gave a good protection (59.5%) at 6 h interval. This kind of behaviour of NaNO 2 could not be detected and explained solely through EIS which raised the importance of DEIS in corrosion measurements. Another important inference derived from the DEIS data was that this inhibitor continued to be effective even after 24 h. 3. The adsorption studies suggested spontaneous and predominant physiosorption of sodium nitrite on copper in SCW. The adsorption studies also suggested that each active site on the copper surface was associated with more than one adsorbing species. This could be attributed to formation of either Cu(OH) 2 or Cu(NO 2 ) 2 . The studies also suggested the lesser interaction among the adsorbed species. The comprehensive type adsorption suggested by El-Awady isotherm indicated some role of chemisorption in the process too, which might be due to the formation of [Cu/Inhibitor] complex using oxygen or nitrogen. 4. The inhibitor NaNO 2 was really effective in improving the surface homogeneity as suggested by AFM and OP analysis. 5. One of the most interesting inferences was derived from SEM/EDX analysis. The inhibition which was observed cannot be attribute solely to the formation of passive layer Cu 2 O + CuO/Cu(OH) 2. Considerable amounts of nitrogen and oxygen detected through EDX analysis in inhibited specimen suggested an association of nitrite ion to the copper surface. 6. The parameters derived by the HOMO and LUMO energy of system suggest that nitrite ion has a low electrophilicity value and high electron donating power. The molecular dynamic simulation indicated a powerful interaction between nitrite ion and copper surface i.e. sodium nitrite might be regarded as an anodic inhibitor, but through experimentations it has been observed that it may act as film forming inhibitor for copper in simulated cooling water. 7. All the analysis performed to establish the role of sodium nitrite in protecting copper in SCW suggested that it was a very efficient inhibitor for cooling water systems fabricated with copper.

Methods
Materials and media. The copper samples certified to be 99.99% pure obtained from Sofıa Med Ad, Bulgaria, were utilized as test specimens. The test specimens' working surface was abraded by SiC papers ranging from 360 to 2000 grit, degreased with alcohol and subsequently washed with distilled water. The simulated cooling water was prepared by dissolving NaCl (7.5 mmol/L), NaHCO 3 (2.0 mmol/L), Na 2 SO 4 (3.5 mmol/L), MgSO 4 (0.25 mmol/L) and CaCl 2 (0.50 mmol/L) in distilled water 38 .
Electrochemical impedance spectroscopy (EIS). The experiments were performed in a system where Ag/AgCl acted as reference electrode, a counter electrode of platinum (Pt) and working electrode of copper with exposed surface area 0.785 cm 2 formed the investigation cell. The steady state OCP was obtained at interval of two hours. EIS was conducted at 10 mV amplitude signal peak-to-peak corrosion potential (E corr ) with a 10 mHz www.nature.com/scientificreports/ to 100 kHz of frequency range 39 . The output was processed by the ZsimpWin 3.21 software. The % inhibition efficiency 'IE' values were computed according to the following equation.
where R ct (i) represents the charge transfer resistance in the presence of NaNO 2 while R ct is the value for uninhibited solution, respectively. Dynamic electrochemical impedance spectroscopy (DEIS). The effectiveness of inhibitors against scaling and corrosion of copper over longer time periods can be monitored through DEIS. DEIS was performed using current perturbation signal given by National Instrument Ltd PCI-4461 digital-analog card and acquisition of voltage response signal of copper sample. The galvanostatic conditions were established by Slepski Galvanostat which also served as current-voltage converter during the process 40 . Multisinusoidal signal, which consisted of 20 elementary current sinusoids of optimized phase shift and amplitude value, was used as perturbation source. The sampling frequency was 12.8 kHz, and the frequency range from 300 mHz to 4.5 kHz, with 8 points per decade of frequency was applied. The DEIS assures the live assessment of the corrosion process on working electrode. The impedance spectra obtained by DEIS were analyzed with the same program that interpreted the EIS data.
Surface analysis. Visual realization of the degree of corrosion occurring on the copper specimens' surface in terms of heterogeneity/roughness requires high precision techniques like AFM and SEM/EDX. The current experiments were conducted on copper specimens in inhibited and uninhibited conditions at room temperature. SEM/EDX study was done by FEG250 (FEI, Holland) with EDX attachment. The instrument for AFM was 'Park Systems XE-100E model. The scan rate was 0.2 Hz and optical resolution was 1 μm. The area of (9.09 × 9.09) μm 2 of test coupon was visualized in non-contact mode. The optical profilometry was conducted by Phase View optik profiler.
Computational methods. In the Conceptual Density Functional Theory introduced by Parr and Pearson, descriptors of chemical reactivity like chemical potential (µ), electronegativity (χ), hardness (η) and softness (σ) are derived with respect to number of electron of total electronic energy at a constant external potential. Softness is the multiplicative inverse of the hardness. Electronegativity is reported as the negative value of chemical potential. Within the framework of finite differences approach, the relation with ground state ionization energy (I), ground state electron affinity (A), total electronic energy E and the number of electrons (N) of the aforementioned descriptors are presented via the following equations 41 .
Koopmans Theorem 42 can be considered as a bridge between Molecular Orbital Theory and Conceptual Density Functional Theory. Koopmans theorem states that the negative values of HOMO and LUMO orbital energies relate to ionization energy and electron affinity of molecules is a parallel approach to predict the molecular ionization energy and electron affinity values. In the light of this theorem, the following relations I = −E HOMO and A = −E LUMO are given. According to electrophilicity index (ω) given by Parr, Szentpaly and Liu 43 , electrophilic power of chemical species can be predicted using their electronegativity (or chemical potential) and hardness values. After this well-known study, Chattaraj proposed the nucleophilicity (ε) as the multiplicative inverse of electrophilicity.
In recent years, Gazquez and coworkers 44 made detailed studies regarding to the prediction of electron donation and electron accepting abilities of molecules and they proposed two new parameters called as electron donating power (ω − ) and electron accepting power (ω + ) and formulated these parameters based on ionization energy and electron affinity concepts via following equations. www.nature.com/scientificreports/ The fraction of electrons transferred (ΔN) from inhibitor molecule to metal surface is an important parameter. The equation proposed for the calculation of this parameter has been suggested considering electronegativity equalization principle imparted to science by Sanderson. Until recently, for the calculation of the fraction of electrons transferred, electronegativity of metal atom was in use but some researchers proposed the use of work function determined for metal surface instead of metal electronegativity. The equation for the calculation of ΔN values of adsorption processes is given as: wherein, it is important to note that χ inh , η Cu and η inh stand for electronegativity of inhibitor, hardness of metal and hardness of inhibitor, respectively. In the calculations, η Cu = 0 is taken assuming that for a metallic bulk I = A. Work function value reported for Cu (111) surface is 4.80 eV 45 . Molecular dynamics simulation. Adsorption properties of nitrite ion (NO 2 − ) were estimated using Molecular Dynamics Simulation Approach on Forcite module from Accelrys Inc. As model metal surface Cu (111) was selected. To see the effect of solvent, the calculations were made in both vacuum and aqueous media. The calculations were made in a simulation box where nitrite ion is in a contact with copper surface at 303 K. Molecular Dynamic simulation approach is well established for the analysis of metal-inhibitor interaction. It is notably important that in the calculation of interaction (E interaction ) and binding (E binding ) energies regarding to adsorption process of nitrite ion on Cu (111) surface, the following equations are considered 46 .
In the given equations, E total stands for the total energy of all system. E surface+H2O represent the total energy of Cu (111) surface with H 2 O molecules. E inhibitor is the total energy of inhibitor alone.