DFT Computational and Experimental Study of Electroless Nickel Deposition

The kinetic mechanism of the electroless deposition of NiP from hypophosphite solutions was investigated from first principles with the intent of determining a set of elementary reactions that could effectively identify the most relevant steps of this complex process. The outcome of the computational approach is discussed in terms of experimental results obtained by electrochemical methods and surface analysis techniques [1,2]. The computational study was performed without considering explicitly the effect that complexing agents introduced in solution may have on the deposition kinetics. Simulations were performed with Density Functional Theory DFT using the B3LYP exchange and correlation functionals [3], the 6-311+G(d,p) basis set for P,C,H and O atoms [4], and the Stuttgart effective core potential basis sets for Ni [5]). The interaction of Ni++ and H2PO2- with the Ni surface was studied using Ni clusters of different sizes (from 4 to 12 atoms). Solvation in water was modeled using the implicit polarizable continuum model [6]; calculations were performed using the G09 computational suite [7]. Simulations were initially devoted to the identification of the first reaction steps that lead to the adsorption of Ni++ on the Ni surface and to its successive reduction by means of the hypophosphite anion. Though several theories were proposed to describe this key step in the literature, our search was guided by the original insight published several years ago [8,9], according to which the reduction of Ni++ requires the formation of an intermediate species containing the NiOH group. Such species would then take active part in the transfer of the OH group to H2PO2-, thus favoring its oxidation to H2PO3-. In addition, it was imposed that the searched elementary mechanism satisfies the experimental observation [10], obtained with marked deuterium species, that each molecule of hydrogen that is evolved contextually to the deposition of NiP contains two H atoms coming from two H2PO2- anions. The possibility that NiOH+ is formed in solution prior to deposition was investigated determining computationally the pKa of Ni++. Results revealed that the formation of NiOH+ is not favored in the pH ranges in which NiP is usually deposited. A similar result was obtained when considering the possibility that NiOH+, after being formed in solution, is stabilized through the formation of a complex with H2PO2-. A thermodynamically favored set of reactions was instead found when it was