Complexation of Nickel Ions by Boric Acid or
(Poly)borates

Anais Graff, Etienne Barrez, Philippe Baranek, Martin Bachet, Pascale

Bénézeth

► To cite this version:

Anais Graff, Etienne Barrez, Philippe Baranek, Martin Bachet, Pascale Bénézeth. Complexation of Nickel Ions by Boric Acid or (Poly)borates. Journal of Solution Chemistry, 2017, 46 (1), pp.25-43. 10.1007/s10953-016-0555-x. hal-02123959

HAL Id: hal-02123959

https://hal.science/hal-02123959v1

Submitted on 21 Jul 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Complexation of Nickel Ions by Boric Acid
or (Poly)borates

Anais Graff<fnref n="1" /> · Etienne Barrez<fnref n="1" />,<fnref n="2" /> · Philippe Baranek<fnref n="1" />,<fnref n="3" /> · Martin Bachet<fnref n="1" /> · Pascale Bénézeth<fnref n="4" />

Received: 4 August 2016/Accepted: 3 October 2016/Published online: 17 December 2016
© The Author(s) 2016. This article is published with open access at Springerlink.com

**Abstract** An experiment based on electrochemical reactions and pH monitoring was performed in which nickel ions were gradually formed by oxidation of a nickel metal electrode in a solution of boric acid. Based on the experimental results and aqueous speciation modeling, the evolution of pH showed the existence of significant nickel–boron complexation. A triborate nickel complex was postulated at high boric acid concentrations when polyborates are present, and the equilibrium constants were determined at 25, 50 and 70 °C. The calculated enthalpy and entropy at 25 °C for the formation of the complex from boric acid and $Ni^{2+}$ ions are respectively equal to (65.6 ± 3.1) kJ·mol$^{-1}$ and (0.5 ± 11.1) J·K$^{-1}$·mol$^{-1}$. The results of this study suggest that complexation of nickel ions by borates can significantly enhance the solubility of nickel metal and nickel oxide depending on the concentration of boric acid and pH. First principles calculations were investigated and tend to show that the complex is thermodynamically stable and the nickel cation in solution should interact more strongly with the $B_3O_3(OH)_4^-$ than with boric acid.

**Keywords** Aqueous · Complexation · Boric acid · Nickel · Equilibrium constant · First principles · DFT

<figure><img src="image_3.png" /></figure>

<footnote>Anais Graff
anais.graff@edf.fr</footnote>

<footnote n="1">¹ EDF R&D – Department of Material and Mechanic of Component, EDF Lab Les Renardières, Avenue des Renardières – Ecuelles, 77818 Moret-Sur-Loing Cedex, France</footnote>

<footnote n="2">² Present Address: PPSM, ENS Cachan, CNRS, Université Paris – Saclay, 94235 Cachan, France</footnote>

<footnote n="3">³ Present Address: EDF R&D, Department Economic and Technical Analysis of Energy Systems (EFESE) – Institut of Research and Development on Photovoltaic Energy (IRDEP), EDF Lab – Chatou, 6 quai Wattier, 78400 Chatou Cedex, France</footnote>

<footnote n="4">⁴ Géosciences Environnement Toulouse, CNRS, Observatoire Midi-Pyrénées, Université de Toulouse, 14 avenue Edouard Belin, 31400 Toulouse, France</footnote>
# 1 Introduction

The primary circuit of a pressurized water reactor (PWR) is chemically conditioned with boric acid as neutron adsorber [1]. Its hydrolysis giving borate ions has a major effect on the pH of the routine power plant cycles. Nevertheless, in certain conditions of temperature and concentration, other species, such as polyborates, can also be formed and need to be taken into account [2, 3]. In addition, the PWR primary circuit is subject to corrosion issues linked to the formation of different oxides, called corrosion products. Some particles of these oxides are released in the primary fluid and may be activated when they pass through the core of the reactor. Even if the quantity of these particles is very low, the radiation dose rates can be strongly influenced by their deposition on the whole primary circuit. Nickel comes from steam generator tubing alloys and is a corrosion product of major importance because of its activation into ⁵⁸Co and its impact on the radioactive contamination of the circuit. Precise knowledge of the solubility of corrosion products containing nickel and how it changes with temperature and chemistry would be valuable to understand the transport of nickel from steam generator tubes to the core. Previous experimental studies have focused on the solubility of nickel oxide and nickel metal at high temperatures and pressures in different media including PWR operating conditions but the results are subject to discrepancies [4–8]. One explanation could be the complexation of nickel ions by (poly)borates, which can increase the solubility of the solid phases.

Two principal studies have focused on the nickel–boron complex formation. Palmer et al. have studied the solubility of nickel hydroxide in different media, including low concentrations of boric acid [6]. The nickel concentrations were found to be higher in boric acid medium than all the remaining data collected in the study. However, no hypothesis was offered to explain that solubility enhancement, especially because results were in the reverse order with respect to boric acid concentration if a reaction of complexation occurred. The same researchers mentioned earlier in an EPRI report [9] that a bidendate complex could be the cause of these results by analogy with the study of the aluminate complexation by the bis–tris,2,2-bis(hydroxymethyl)-2,2',2''-nitrilotriethanol used as pH buffer [10]. They suggested the reaction:

$$Ni(OH)_2 + B(OH)_4^- \rightleftharpoons Ni(O)_2B(OH)_2^- + 2H_2O \quad (1)$$

but this explanation would be correct only if the complex is very strong, in order to raise the solubility of nickel hydroxide by orders of magnitude.

On the other hand, Shchigol has studied the solubility of the solid nickel orthoborate Ni(BO₂)₂·4H₂O in boric acid medium [11]. The measured solubility was found to increase upon addition of boric acid. To explain that enhancement, the formation of a nickel borate complex in solution was believed to occur as Ni(BO₂)₂⁻. However, the study was con- troverted [12] because the cation balance was not respected. In addition, the solid was not well characterized, the equilibrium state may not have been reached (lack of data), nickel ion complexation by chlorides or hydroxide ions had not been taken into account and the activity coefficient model was not defined. Nevertheless, a recent review of nickel chemical thermodynamics [12] suggests that a neutral complex could be formed based on Shchigol’s experimental data and according to the following reactions:

$$Ni(BO_2)_2 \cdot 4H_2O(s) \rightleftharpoons Ni^{2+} + 2B(OH)_4^- \quad (2)$$
$$ \mathrm{Ni}^{2+} + 2\mathrm{B(OH)}_{4}^{-} + \mathrm{B(OH)}_{3(\mathrm{aq})} \rightleftharpoons \mathrm{NiH}(\mathrm{BO}_{2})_{3(\mathrm{aq})} + 5\mathrm{H}_{2}\mathrm{O} \quad (3) $$

Other papers [13, 14] discussed the mechanism of a nickel–borate complexation but they did not use experimental results. Consequently, valuable experimental data are lacking to confirm the existence of a nickel–boron complex. The aim of this study is, therefore, to experimentally study complexation in solution and to determine the associated equilibrium constants at 25, 50 and 70 °C by using a pH-monitoring method and modeling. The second part of the paper is devoted to the characterization of the complex by first principles calculations, supporting its stability in water.

## 2 Proposition of Nickel–Boron Complex Based on Shchigol’s Experimental Data

Before proposing a nickel–boron complex, the speciation of boron needs to be taken into account to highlight which boron species could be relevant as a ligand. The proposed complex must also fit the experimental data obtained by Shchigol [11]. For that purpose, the software CHESS was used as a modelling tool [15]. Data used for all calculations are available in Table 1. Finally, since the complex strength derived from Shchigol’s work depends on the solubility of the solid that was present, it is necessary to experimentally confirm the presence of the complex in solution, independent from the data of Shchigol.

### 2.1 Boron Speciation

The speciation of boron strongly depends on the chemical medium and the experimental conditions such as pH, temperature, boron concentration, and counter ions. At concentrations of boron higher than 0.01 mol·kg⁻¹, different polymerized species, called polyborates, are formed and can even become dominant in solution. Various studies [2, 3, 16–19] have been conducted to understand the mechanisms of formation and the geometric configuration of polyborate species but their diversities make this study complicated and most of the results are controversial. The studies from Palmer et al. [3] and Mesmer et al. [2] about the determination of polyborate formation constants are considered as reliable because of their relatively good agreement and the existence of experimental data. Palmer et al. propose a dataset that includes either a tetraborate or a pentaborate. In this study we did not take into account the pentaborate since there are discrepancies even on its formula in the literature. Therefore, the following reactions of polyborates formation were used:

$$ \mathrm{B(OH)}_{3(\mathrm{aq})} + \mathrm{H}_{2}\mathrm{O} \rightleftharpoons \mathrm{B(OH)}_{4}^{-} + \mathrm{H}^{+} \quad (4) $$

$$ 2\mathrm{B(OH)}_{3(\mathrm{aq})} \rightleftharpoons \mathrm{B}_{2}\mathrm{O}(\mathrm{OH})_{5}^{-} + \mathrm{H}^{+} \quad (5) $$

$$ 3\mathrm{B(OH)}_{3(\mathrm{aq})} \rightleftharpoons \mathrm{B}_{3}\mathrm{O}_{3}(\mathrm{OH})_{4}^{-} + \mathrm{H}^{+} + 2\mathrm{H}_{2}\mathrm{O} \quad (6) $$

$$ 4\mathrm{B(OH)}_{3(\mathrm{aq})} \rightleftharpoons \mathrm{B}_{4}\mathrm{O}_{5}(\mathrm{OH})_{4}^{2-} + 2\mathrm{H}^{+} + 3\mathrm{H}_{2}\mathrm{O} \quad (7) $$

The results from Palmer et al. [3] are given in terms of base hydrolysis reactions, while the CHESS thermodynamic database requires reactions involving only “basic/elementary species” (in this case B(OH)₃, H⁺ and H₂O). Furthermore, Palmer et al. provided the equations of the equilibrium quotients, which involve additional parameters linked to the
<table><caption>**Table 1 CHESS parameters used for simulation**</caption><thead><tr><th>Reactions</th><th>log<sub>10</sub> K</th></tr></thead><tbody><tr><td>B(OH)<sub>3</sub>(aq) + H<sub>2</sub>O ⇌ B(OH)<sub>4</sub><sup>-</sup> + H<sup>+</sup></td><td>log<sub>10</sub> K = -36.2605 + <sup>365.18</sup>/<sub>T</sub> + 11.6402 log<sub>10</sub> T + (16.4914 - 0.023917T) log<sub>10</sub>(ρw) + log<sub>10</sub>(Kw) with pw the density of water in cm<sup>3</sup>·ml<sup>-1</sup></td></tr><tr><td>2B(OH)<sub>3</sub>(aq) ⇌ B<sub>2</sub>O(OH)<sub>5</sub><sup>-</sup> + H<sup>+</sup></td><td>log<sub>10</sub> K = -3.935 + <sup>1780.5</sup>/<sub>T</sub> + 0.95183 log(T) + log<sub>10</sub>(Kw)</td></tr><tr><td>3B(OH)<sub>3</sub>(aq) ⇌ B<sub>3</sub>O<sub>3</sub>(OH)<sub>4</sub><sup>-</sup> + H<sup>+</sup> + 2H<sub>2</sub>O<sup>a</sup></td><td>log<sub>10</sub> K = -6.495 + <sup>3219.1</sup>/<sub>T</sub> + 0.95186 log(T) + log<sub>10</sub>(Kw)</td></tr><tr><td>4B(OH)<sub>3</sub>(aq) ⇌ B<sub>4</sub>O<sub>5</sub>(OH)<sub>4</sub><sup>2-</sup> + 2H<sup>+</sup> + 3H<sub>2</sub>O<sup>a</sup></td><td>log<sub>10</sub> K = -5.031 + <sup>600.3</sup>/<sub>T</sub> - 1.3572 log<sub>10</sub>(T) + 2log<sub>10</sub>(Kw)</td></tr><tr><td>Ni<sup>2+</sup> + 2H<sub>2</sub>O ⇌ Ni(OH)<sub>2</sub>(aq) + 2H<sup>+</sup></td><td>log<sub>10</sub> K<sub>80</sub> = -0.335 - <sup>534.26</sup>/<sub>T</sub> - 4.2434 × 10<sup>-2</sup> × T</td></tr><tr><td>Ni(OH)<sub>2</sub>(S) + 2H<sup>+</sup> ⇌ Ni<sup>2+</sup> + 2H<sub>2</sub>O<sup>b</sup></td><td>log<sub>10</sub> K'<sub>80</sub> = -2.829 + <sup>4320.17</sup>/<sub>T</sub></td></tr><tr><td>Ni(B(OH)<sub>4</sub>)<sub>2</sub>(S) + 2H<sup>+</sup> ⇌ Ni<sup>2+</sup> + 2B(OH)<sub>3</sub>(aq) + 2H<sub>2</sub>O<sup>c,d</sup></td><td>log<sub>10</sub> K<sub>8</sub>(25°C) = 9.2</td></tr><tr><td>3B(OH)<sub>3</sub>(aq) + Ni<sup>2+</sup> ⇌ NiB<sub>3</sub>O<sub>4</sub>(OH)<sub>3</sub>(aq) + 2H<sup>+</sup> + 2H<sub>2</sub>O<sup>d</sup></td><td>log<sub>10</sub> K<sub>1</sub> = (-3424.6 ± 180.2)/T + (0.028 ± 0.557)</td></tr></tbody></table>

The reactions and the associated equations of the equilibrium constant at infinite dilution were taken from a Palmer et al. [3]; b the EPRI report [9]; c this study, fit to Shchigol et al. [11]; d this study, fit to experimental data given in Table 3
ionic strength for a given medium. At infinite dilution, these parameters can be removed and the activity of water is equal to 1. The associated equations of the equilibrium constants were therefore combined with the dissociation of water in order to agree with Eqs. 4–7. Since the ionic strength effect is calculated by CHESS with a model that is not the same as Palmer’s, a small deviation results on the final calculations, but because of its low value in our chemical conditions, we did not observe any significant effect of this inconsistency. The equilibrium constant equations for reactions 4–7 are available in Table 1.

These species are important for understanding the behavior of boron as a function of our chemical conditions. Figure 1 presents the speciation of boron at two different concentrations (A) 0.5 mol·kg⁻¹ and (B) 0.1 mol·kg⁻¹ and at 25 °C using the equilibrium constants calculated from the fitting equation adapted from Palmer et al. [3] and available in Table 1. It appears that polyborate species exist in solution for pH ranging from 6 to 10. The triborate species is the most abundant polyborate compared to the others and can even become a major species at 0.5 mol·kg⁻¹ of boron (Fig. 1a), reaching 35% in solution. In addition, the existence of polyborate species depends also on the total concentration of

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 1** Boron speciation diagram as a function of pH at 25 °C: boron concentration equal to **a** 0.5 mol·kg⁻¹, **b** 0.1 mol·kg⁻¹</figcaption></figure>
boron, becoming dominant when the concentration of boron increases to >0.5 mol·kg⁻¹. According to Fig. 1b, at 0.1 mol·kg⁻¹ of boron, the concentration of the triborate is reduced to less than 10%. Consequently, depending on the chemical conditions, three species are able to lead to the formation of a nickel–boron complex with a quantifiable amount in solution, which are the borate, triborate and boric acid species. Other boron species, such as diborate or tetraborate could nevertheless complex with nickel ions, but regarding to their small concentration (>2% of total boron), it could be difficult to monitor them experimentally.

## 2.2 Proposition of a Complex

Considering the speciation of boron, the use of high boric acid concentrations, up to 0.7 mol·kg⁻¹ and the pH range from 6 to 7.4 in Shchigol’s experiments [11], we can assume that polyborates are abundant under those chemical conditions. Borate ions are minor because of the acidity of the medium. As a consequence, the triborate species is chosen to be a relevant ligand. A neutral complex could be formed, as proposed by Gamsjäger et al. [12] in their review and described in Eq. 8 below:

$$B_3O_3(OH)_4^- + Ni^{2+} \rightleftharpoons NiB_3O_4(OH)_{3(aq)} + H^+ \quad (8)$$

The equation can be written differently to fit CHESS thermodynamic database requirements:

$$3H_3BO_3(aq) + Ni^{2+} \rightleftharpoons NiB_3O_4(OH)_{3(aq)} + 2H^+ + 2H_2O \quad (9)$$

with:

$$K_1 = \frac{[NiB_3O_4(OH)_{3(aq)}] \cdot \gamma_{NiB_3O_4(OH)_{3(aq)}} \times [H^+]^2 \cdot \gamma_{H^+}^2 \times a_w^2}{[Ni^{2+}] \cdot \gamma_{Ni^{2+}} \times [H_3BO_3]^3 \cdot \gamma_{H_3BO_3(aq)}^3} \quad (10)$$

where $K_1$ is the equilibrium constant, $a_w$ the activity of water, and $\gamma$ is the activity coefficient of the species in solution, calculated by the truncated Davies model, which has been preferred over the Pitzer model as, to our knowledge, no coefficients representing the interaction between the nickel cation and the boric acid species are reported in the literature. The most important contribution to the ionic strength value is from the formation of nickel ions, $Ni^{2+}$ ($Z_i^2 = 4$), and not from the acid nor from a salt (such as NaCl) as we did not want to use an electrolyte to fix the ionic strength to a constant value in our experiments. The main reason is that we did not want to introduce some species ($Na^+$, $Cl^-$) known to form complexes with boron or nickel species. This results in an increase of the ionic strength from $10^{-4}$ to $10^{-2}$ mol·kg⁻¹ upon the addition of nickel ions to the solution during a trial.

The solubility of the solid nickel orthoborate with increasing boric acid concentration was simulated with CHESS in order to determine whether the proposed triborate complex would be consistent with Shchigol’s experimental data [11] (the solubility constant of the solid used for these calculations was previously determined by fitting to Shchigol’s experimental data in dilute HCl, being $log_{10} K_S$ (25 °C) = 9.2 for the reaction reported Table 1). Results of the computation are represented in Fig. 2. The modeling agrees relatively well with experimental data at 0.5 mol·kg⁻¹ of boron when $log_{10} K_1$ is equal to -11.1, but the agreement could be better since a part of the curve is outside of the
<figure><img src="image_6.png" /><figcaption>**Fig. 2** Variation of nickel ions concentration as a function of boric acid concentration, including the formation of a triborate–nickel complex (CHESS speciation model, $log_{10} K_1 = -11.1$)</figcaption></figure>

experimental error bars. Nevertheless, those modeling results were the best and the closest obtained after trying several values of the constant $K_1$ for the proposed complexation reaction. For a better comparison, further informations is needed on the experimental investigations of Shchigol [11]. To our knowledge, the author did not specify the ionic strength model used and the synthesized solid nickel orthoborate was not characterized. In addition, other potential neutral complexes were tested by following the same simulation protocol, with different ligands such as diborate or borate species, but we did not succeed in obtaining satisfactory agreement between the experimental data and the modeling. Table 2 describes the residual errors associated to the best $log_{10} K$ determination obtained between Shchigol’s experimental data [11] and CHESS simulations for the complexation reactions involving these different ligands. These results showed that the smallest residual error, meaning the closest simulation curve to Shchigol’s data [11], was obtained with the triborate simulation. As an example, computation result for a diborate complex is also shown in Fig. 2. It did not fit the experimental data at low concentrations of nickel and the slopes are significantly different. We did not model a complex with a different charge, to be consistent with the assumption of a neutral complex. According to Fig. 3, this hypothesis is also confirmed by modeling the pH of the solution as a function of boric acid concentration, where experimental data and simulation are in quite good agreement, with the assumption of a triborate complex. Consequently, the proposed reaction of the Ni–B complexation respects the previous enunciated criteria: it is a tri-boron neutral complex, which takes into account the existence of polyborates and also matches reasonably with the experimental data of Shchigol, but this hypothesis must be experimentally confirmed independently from the literature data. According to Eq. 8, a proton is released by the complexation reaction, with an expected impact on pH. The CHESS model presented in

<table><caption>**Table 2** Residual errors obtained between the experimental data of Shchigol [11] and CHESS simulations of the solubility of nickel orthoborate with increasing boric acid concentration, involving different Ni–B complexation reactions</caption><thead><tr><th>Complexation reactions</th><th>Best simulated log<sub>10</sub> K</th><th>Residual error</th></tr></thead><tbody><tr><td>Ni<sup>2+</sup> + B(OH)<sub>3(aq)</sub> ⇌ NiBO<sub>2</sub>(OH)<sub>3(aq)</sub> + 2H<sup>+</sup></td><td>-16.6</td><td>1.6 × 10<sup>-4</sup></td></tr><tr><td>Ni<sup>2+</sup> + 2B(OH)<sub>3(aq)</sub> ⇌ NiB<sub>2</sub>O<sub>3</sub>(OH)<sub>2(aq)</sub> + 2H<sup>+</sup> + H<sub>2</sub>O</td><td>-11.5</td><td>3.8 × 10<sup>-5</sup></td></tr><tr><td>Ni<sup>2+</sup> + 3B(OH)<sub>3(aq)</sub> ⇌ NiB<sub>3</sub>O<sub>4</sub>(OH)<sub>3(aq)</sub> + 2H<sup>+</sup> + 2H<sub>2</sub>O</td><td>-11.1</td><td>1.4 × 10<sup>-5</sup></td></tr></tbody></table>
<figure><img src="image_7.png" /><figcaption>**Fig. 4** pH variation as a function of $Ni^{2+}$ concentration in boric acid media calculated with CHESS speciation model for a boron solution at 0.5 $mol·kg^{-1}$, initial pH = 4.2, t = 25 °C</figcaption></figure>

<figure><img src="image_8.png" /></figure>

Fig. 4 shows pH variations as a function of the concentration of nickel ions in boric acid medium. It has been obtained for a $0.5 mol·kg^{-1}$ boric acid solution by using the formation constant determined above ($log_{10} K_1 = -11.1$) as well as the formation constant of polyborates determined by Palmer et al. [3]. It can be seen in Fig. 4 that there is a decrease in pH when the complexation of nickel ions by the triborate occurs. Therefore, it is possible to confirm this model experimentally by pH-monitoring.

# 3 Experimental and First Principles Approaches to Study the Nickel–Boron Complexation

## 3.1 Experimental Approach

An experiment based on electrochemical reactions and pH monitoring was performed in a double wall reactor. Nickel ions were gradually formed by oxidation of a nickel metal electrode in a solution of boric acid. The apparatus, schematically represented in Fig. 5, was constituted of four different electrodes. Nickel ions come from a nickel metal electrode (WE) used as anode (Goodfellow, purity 99.98%, 1 mm in diameter, surface of 25 $cm^2$). The oxidation reaction is:

$$Ni \rightarrow Ni^{2+} + 2e^{-} \qquad (11)$$

The counter electrode (CE, cathode), where the reduction reaction occurred, is a platinum wire from Radiometer Analytical. The reaction of reduction is:
<figure><img src="image_9.png" /><figcaption>Fig. 5 Schematic of the double wall reactor with potentiostat and pH-monitoring, WE is the nickel metal working electrode, CE is the platinum counter electrode and RE is the saturated calomel reference electrode</figcaption></figure>

$$2\text{H}^+ + 2\text{e}^- \rightarrow \text{H}_2 \quad (12)$$

A saturated calomel electrode from Radiometer Analytical was used as reference (RE) and the pH was measured by a combined Ag/AgCl glass electrode from VWR as a function of time and nickel concentration using a Metrohm 780 pH meter. The electrode was calibrated with buffers at pH values of 4, 7 and 10 (from Metrohm) before each experiment. The electrode was also checked after the experiment to ensure that no potential deviation higher than the pH uncertainty had occurred. The nickel electrode was polished with a 1000 grid paper and endured a reduction cycle at –1.3 V versus SCE during 5 min before each test, to avoid the presence of NiO at the surface. The solution was continuously bubbled with nitrogen and the temperature was controlled within ±1 °C by a thermostat. According to the nickel Pourbaix diagram, nickel ions are formed in the potential range –0.2 to 0.8 V versus the hydrogen electrode (1 < pH < 7). The current was fixed at 2 mA by galvanostatic control during each experiment. This value was chosen after running several tests: the higher the current the faster the nickel ions are formed in solution. For values higher than 2 mA, the potential was out of the limits of the Ni²⁺ formation domain, whereas lower values lead to long run times (>48 h). Furthermore, the ionic strength was not fixed by salt addition, to avoid complexation with any cations other than nickel. For that reason, the impedance of the electrochemical circuit was high but this did not disturb the process of nickel ion formation in solution. All the collected data were calculated from the species activities determined by the CHESS modeling with the truncated Davie’s model as mentioned earlier. Results were obtained at 25, 50 and 70 °C.

For each temperature, runs were performed at concentrations 0.5 mol·kg⁻¹ of boric acid. Lower concentrations of boron increase the working pH range during a run. This leads to precipitation of nickel hydroxide according to its solubility diagrams. Boron solutions were obtained from boric acid powder (Alfa Aesar puratronic® 99.9995% purity) and concentrations were checked by titration with NaOH. All solutions were prepared from MilliQ water. During each run, 2 mL of the solution were periodically sampled, then
filtered (0.2 µm cellulose acetate filter) to analyze the concentration of nickel by ICP-MS (Varian 820 MS).

## 3.2 First Principles Approach

First principles calculations have been performed in order to determine different properties of the nickel–boron complexes using the NWCHEM code [20]. Their equilibrium stoichiometries and their formation's energies with respect to the reactions of complexation have been investigated at the B3LYP level [21, 22]. The “solvation model based on density” (SMD) method [23] was used to take into account the solvent effect. Regarding the basis set used to describe the electronic structure of the various molecules and complexes, Gaussian type functions (GTF) have been employed. The H, B [24] and O [25] atoms have been described with all electron GTF basis sets. Standard basis with 5s-111sp-1p*, 6s-311sp-1d* and 8s-611sp-1d* have been adopted for H, B and O, respectively. For Ni, the 3s, 3p, 3d and 4s have been treated as valence electrons combined with a Hay–Wadt small core pseudopotential as described in the literature [26–28] and a 3111sp-311d basis set. A full description of the different basis is given in the Appendix.

For the evaluation of the exchange correlation contribution to the density functional, the “xfine” grid, as defined in the manual of NWCHEM [20], was used. The convergence criterion on total energies was $10^{-8}$ au. Atomic displacements and force thresholds were $1.8 \times 10^{-3}$ and $4.5 \times 10^{-4}$ au, respectively. The precision of the obtained energies of reaction and bond length are $1 \times 10^{-6}$ au and less than 0.01 Å, respectively. With these computational conditions, the obtained data can be considered to be fully converged.

## 4 Results and Discussion

The pH variations as functions of nickel ion concentration were obtained at three different temperatures (25, 50 and 70 °C). Experimental data are available in Table 3. Nickel concentrations were an average of five replicates. The associated uncertainty was obtained

<table><caption>**Table 3** Experimental results of the pH variation obtained as a function of the molality of nickel ions</caption><thead><tr><th colspan="3">25 °C</th><th colspan="3">50 °C</th><th colspan="3">70 °C</th></tr><tr><th>Time (min)</th><th>pH</th><th>log<sub>10</sub> [Ni] (mol·kg<sup>-1</sup>)</th><th>Time (min)</th><th>pH</th><th>log<sub>10</sub> [Ni] (mol·kg<sup>-1</sup>)</th><th>Time (min)</th><th>pH</th><th>log<sub>10</sub> [Ni] (mol·kg<sup>-1</sup>)</th></tr></thead><tbody><tr><td>0</td><td>4.11</td><td>-4.15 ± 0.01</td><td>0</td><td>4.18</td><td>-8.35 ± 0.08</td><td>0</td><td>3.70</td><td>-8.26 ± 0.07</td></tr><tr><td>8</td><td>4.33</td><td>-4.32 ± 0.03</td><td>5</td><td>4.38</td><td>-4.50 ± 0.06</td><td>7</td><td>4.20</td><td>-4.21 ± 0.04</td></tr><tr><td>20</td><td>4.52</td><td>-4.05 ± 0.02</td><td>20</td><td>4.59</td><td>-4.08 ± 0.04</td><td>28</td><td>4.45</td><td>-3.87 ± 0.05</td></tr><tr><td>35</td><td>4.73</td><td>-3.82 ± 0.04</td><td>60</td><td>4.91</td><td>-3.70 ± 0.04</td><td>57</td><td>4.75</td><td>-3.65 ± 0.04</td></tr><tr><td>55</td><td>4.85</td><td>-3.65 ± 0.02</td><td>269</td><td>5.47</td><td>-3.07 ± 0.01</td><td>223</td><td>5.32</td><td>-3.20 ± 0.02</td></tr><tr><td>126</td><td>5.18</td><td>-3.30 ± 0.03</td><td>373</td><td>5.58</td><td>-2.85 ± 0.03</td><td>372</td><td>5.44</td><td>-2.91 ± 0.04</td></tr><tr><td>230</td><td>5.44</td><td>-3.04 ± 0.04</td><td>1385</td><td>5.78</td><td>-2.46 ± 0.02</td><td>1440</td><td>5.68</td><td>-2.40 ± 0.01</td></tr><tr><td>358</td><td>5.61</td><td>-2.86 ± 0.05</td><td>1828</td><td>5.80</td><td>-2.36 ± 0.01</td><td>1694</td><td>5.71</td><td>-2.32 ± 0.01</td></tr><tr><td>1370</td><td>6.03</td><td>-2.39 ± 0.02</td><td></td><td></td><td></td><td>1748</td><td>5.73</td><td>-2.30 ± 0.01</td></tr><tr><td>1846</td><td>6.12</td><td>-2.29 ± 0.01</td><td></td><td></td><td></td><td>1779</td><td>5.73</td><td>-2.29 ± 0.01</td></tr><tr><td>2743</td><td>6.10</td><td>-2.27 ± 0.01</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
for a level of confidence at 95%. During each experiment, a small amount of a black deposit was gradually formed on the counter electrode. Scanning electron microscopy identified the solid as metallic nickel. It is probably obtained by the following reaction:

$$ \mathrm{Ni}^{2+} + 2\mathrm{H}^{+} \rightleftharpoons \mathrm{H}_{2} + \mathrm{Ni}_{(\mathrm{s})} \quad (13) $$

but the reaction can be neglected and does not affect the pH measurements.

## 4.1 Determination of the Equilibrium Constants

As predicted, the pH increases with the addition of nickel ions to the solution whatever the chemical conditions. However, to demonstrate that complexation occurs, the pH must increase slower than it does in a non-complexing medium. Figure 6 represents the results for an aqueous concentration of boron equal to 0.5 mol·kg⁻¹. The “no complexation” curves obtained as a function of the temperature and represented by the solid lines are shown for a better understanding. As we can see, pH variations obtained experimentally are different from those predicted by modeling. The higher the temperature, the less the pH increases as a function of nickel ion concentration. The complexation of nickel ions by the triborate is apparent even at 25 °C. From those results, experimental data where fitted by determining the values of the equilibrium constants. Modeling results are represented by the dotted lines and the values of the calculated equilibrium constants are given in Table 4. Uncertainties were calculated by the partial differential equation method. The most significant contribution to the experimental uncertainties is the pH measurement (pH ±0.05).

At the end of the experiment, the pH was raised to 8 by addition of lithium hydroxide (Alfa Aesar puratronic®, 99.99% of purity) in order to observe if a solid precipitates consistent with CHESS calculations that indicate supersaturation for Ni(BO₂)₂·4H₂O(s). After 5 days standing, a green pale precipitate was apparent and was separated from the solution. This precipitate has also been observed by Shchigol, where he postulated the formation of fine crystals of the hexaborate NiB₆O₁₀ from the complex nickel borate solution containing an excess of orthoboric acid and according to:

$$ \mathrm{NiB}_{3}\mathrm{O}_{4}(\mathrm{OH})_{3(\mathrm{aq})}+3\mathrm{H}_{3}\mathrm{BO}_{3(\mathrm{aq})} \rightleftharpoons \mathrm{NiB}_{6}\mathrm{O}_{10(\mathrm{s})}+6\mathrm{H}_{2}\mathrm{O} \quad (14) $$

The deposit was observed by scanning electron microscopy (SEM) and X-ray diffraction (XRD). Agglomerates of thin nanoparticles less than 10 nm in diameter are present. The XRD spectrum shows an amorphous structure partially crystallized. By comparison with solids in the XRD database, several spikes could match with a nickel borate, but nickel

<figure><img src="image_10.png" /><figcaption>**Fig. 6** pH variations with increasing nickel ions concentration in boric acid 0.5 mol·kg⁻¹. The solid lines represent a CHESS speciation model, simulating a media without nickel–boron complexation for the three studied temperatures. The dashed lines represent a simulation of the media where nickel–boron complexation occurred</figcaption></figure>
<table><caption>**Table 4** Equilibrium constant of the nickel–boron complex calculated by CHESS simulation at 25, 50 and 70 °C</caption><thead><tr><th>Complex</th><th>log<sub>10</sub> K<sub>1</sub> (25 °C)</th><th>log<sub>10</sub> K<sub>1</sub> (50 °C)</th><th>log<sub>10</sub> K<sub>1</sub> (70 °C)</th></tr></thead><tbody><tr><td>NiB<sub>3</sub>O<sub>4</sub>(OH)<sub>3</sub>(aq)</td><td>-11.50 ± 0.05</td><td>-10.50 ± 0.06</td><td>-10.00 ± 0.05</td></tr></tbody></table>

hydroxide or bunsenite (NiO) were excluded. Moreover, the atomic composition obtained by energy dispersive X-ray spectrometry shows a quantity of boron 5–7 times higher compared to nickel. As a consequence, the characterization of the precipitate leads to the potential existence of a nickel hexaborate, but other technical approaches must be investigated in order to definitively confirm this hypothesis.

## 4.2 Influence of Temperature

As can be seen from Table 4, the value of the equilibrium constant increases slightly with temperature. Furthermore, as the complexation reaction, according to Eq. 9, is pseudo-isocoulombic (not strictly isocoulombic as equal numbers of like charged species are not present on either side of the reaction with $\Delta z^2 = -2$), the heat capacity change for the reaction could be small. If we consider the heat capacity change small enough to be neglected, which is often the case for strictly isocoulombic reactions, this assumption allows us to presume that the equilibrium constant is insensitive to pressure changes at constant temperature (<250 °C) [29]. Moreover, this assumption has been used or mentioned in previous studies as for example in the papers of Mesmer et al. [30] and Gu et al. [31], in particular when the pseudo-isocoulombic reaction contains only positive charges, a very regular temperature dependence can be observed and therefore the heat capacity of the reaction can be assumed to be zero. Thus, a linear equation of $log_{10} K$ as a function of temperature is often adequate to describe those reactions over wide ranges of temperature and pressure, meaning that the enthalpy $\Delta_r H^0$, and the entropy $\Delta_r S^0$, are constant. The equation of $log_{10} K_1$ was calculated from the modeling results and leads to:

$$ \log_{10} K_1 = (-3424.6 \pm 180.2)/T + (0.028 \pm 0.557) \quad (15) $$

Calculated values of the Gibbs energy change, enthalpy and entropy for the reaction (9) are respectively $\Delta_r G^0_{298K} = (65.4 \pm 0.2) \text{ kJ} \cdot \text{mol}^{-1}$, $\Delta_r H^0_{298K} = (65.6 \pm 3.1) \text{ kJ} \cdot \text{mol}^{-1}$, and $\Delta_r S^0_{298K} = (0.5 \pm 11.1) \text{ J} \cdot \text{K}^{-1} \cdot \text{mol}^{-1}$. It should be pointed out that the extrapolation of Eq. 15 to temperatures higher than 70 °C must introduce some non-negligible inconsistencies, since the reaction is only pseudo-isocoulombic. However, as this work was carried out in order to better understand the behavior of the solid phases of nickel in the PWR primary circuit chemical conditions, an approximation of the thermodynamic data at 300 °C is crucial for further investigations. For example, the results of this study suggest that the complexation reaction significantly enhances the solubility of nickel metal and nickel oxide at pH = 7 and 0.5 mol·kg⁻¹ of boron at 300 °C. Furthermore, the complex formation becomes stronger with increasing temperature; meanwhile the concentration of the triborate ligand decreases by depolymerization. The amount of soluble complex could then be lower at high temperature depending on the chemical medium.
The equilibrium constant and Gibbs energy change were also calculated for reaction (8) at 25 °C according to:

$$ \log_{10} K_{\text{Eq. 8}} = \log_{10} K_{1} - \log_{10} K_{\text{Eq. 6}} \quad (16) $$

and are equal to $\log_{10} K_{\text{Eq. 8}} = -(4.16 \pm 0.05)$ and $\Delta_{\text{r}}G^{0}_{298\text{K}} = (23.7 \pm 0.2) \text{ kJ} \cdot \text{mol}^{-1}$, respectively. Since the reaction is not isocoulombic (or pseudo-isocoulombic), the enthalpy and the entropy cannot be determined without knowing the heat capacity change.

Other experimental approaches, such as $^{11}$B NMR or Raman spectroscopy studies, were tested to characterize the complex’s structure, but we did not obtain conclusive results. One explanation is that the concentration of nickel ions is greatly limited by the solubility of nickel hydroxide when the pH is higher than 6, whereas lower pH lead to the disappearance of the triborate species according to the boron speciation. The amount of complex is consequently too low to be efficiently detected by a conventional apparatus. Although further experimental investigations are needed (EXAFS, XANES or UV spectroscopy), a first principle approach has been used in this study to confirm the complex stability.

## 5 First Principles Approach

The simplest way to test the nickel reactivity with boric acid and triborate is to use the first principles approaches to determine the equilibrium geometries and the internal energies of each molecules. Recently, this approach was used by Tossel et al. [32] to calculate absolute $pK_a$ values for weak acids in aqueous solutions and especially for the boric acid hydrolysis. By carefully using this method, the author showed that it is possible to compare the calculated equilibrium constant with the available experimental data. In this work and for that purpose, the internal energies of the system with the solvent are calculated and simulated with the SMD method [23], taking into account hydration of the molecules in the aqueous system. The internal energy change of the reaction is consequently assumed to be close enough to the Gibbs energy change $\Delta_{\text{r}}G^{0}_{298\text{K}}$, usually presented in experimental work, regarding the uncertainties and some criteria enunciated in [32]. In the present work, the internal energy change of the following complexation reactions will be approximated as the Gibbs energy change and will be called $\Delta G$. They were simulated first in the gas phase but only the results obtained by simulation in the aqueous phase are comparable to our experimental data at 25 °C.

Two different formation reactions of nickel complexes, described by the Eqs. 8 and 9, have been explored: for Eq. 8, the interaction of the nickel with existing triborate reaction $R_1$; for Eq. 9, the interaction of nickel with boric acid reaction $R_2$.

Since at this level the $H^+$ total energy is null (no electron), two ways of treating the energies of reactions 8 and 9 have been explored. The first one called C1, where the experimental value of the $H^+$ energy of solvation ($-1125 \text{ kJ} \cdot \text{mol}^{-1}$ as defined in references [32, 33]) has been used, gives the reactions Gibbs energy changes $\Delta G_{R1}$ and $\Delta G_{R2}$ for Eqs. 8 and 9, respectively:

$$ \Delta G_{R1} = E_{\text{NiB}_3\text{O}_4(\text{OH})_4} - E_{\text{B}_3\text{O}_3(\text{OH})_4} - E_{\text{Ni}^{2+}} - 1125. \quad (17) $$

$$ \Delta G_{R2} = E_{\text{NiB}_3\text{O}_4(\text{OH})_4} + 2E_{\text{H}_2\text{O}} - 3E_{\text{H}_3\text{BO}_3} - E_{\text{Ni}^{2+}} - 1125. \quad (18) $$

For the second one, C2, where $H^+$ is replaced by $\text{H}_3\text{O}^+$ in order to determine $\Delta G_{R1}$ and
$\Delta G_{R2}$ fully theoretically, the Eqs. 8 and 9 have then to be adapted. The resulting Eqs. 19 and 20 are given below and correspond to the reactions $R_1$ and $R_2$, respectively:

$$B_3O_3(OH)_4^- + Ni^{2+} + H_2O \overset{R_1}{\rightleftharpoons} NiB_3O_4(OH)_{3(aq)} + H_3O^+ \quad (19)$$

$$3H_3BO_{3(aq)} + Ni^{2+} \overset{R_1}{\rightleftharpoons} NiB_3O_4(OH)_{3(aq)} + 2H_3O^+ \quad (20)$$

The reactions Gibbs energy changes are then:

$$\Delta G_{R1} = E_{NiB_3O_4(OH)_{3(aq)}} + E_{H_3O^+} - E_{B_3O_3(OH)_4^-} - E_{Ni^{2+}} - E_{H_2O} \quad (21)$$

$$\Delta G_{R2} = E_{NiB_3O_4(OH)_{3(aq)}} + 2E_{H_3O^+} - 3E_{H_3BO_3} - E_{Ni^{2+}} \quad (22)$$

where $E_X$ represents the total energies of molecule X. In this work, these formation energies do not take into account the vibrational and solvent entropies. In order to establish the influence of hydration of the nickel cation, each reaction energy has been determined with $Ni^{2+}$ with or without its first sphere of hydration, and at two levels of approximation: (1) without the solvent (gaseous approximation); (2) with the solvent taken into account via the SMD method. For each of the species, the geometry has been optimized in order to minimize the total energy (their vibrational frequencies have been calculated to verify that the obtained configurations do not correspond to unstable points). For $B_3O_3(OH)_4^-$, the geometry obtained by Zhou et al. [18] was taken as initial geometry for the optimization. For the $NiB_3O_4(OH)_3$ complex, the electrostatic potentials around $B_3O_3(OH)_4^-$ and $Ni^{2+}$·6H$_2$O were used to find the best configuration to start the optimization. Figure 7 illustrates the obtained structures. For the different compounds, the geometries are in good agreement with experiments. For instance, the average distances B–O, O–H (in $H_3BO_3$), O–H (in $H_2O$), Ni–O (O in the first sphere of hydration) and Ni–O in the complexes are 1.373, 0.965, 0.975, 2.054 and 1.845 Å, respectively, coinciding with the various bond length of aqueous borate and nickel solution, and crystal structures of borate and nickel hydroxides that are reported in the literature [18, 34–36].

The results for $\Delta G_{R1}$ and $\Delta G_{R2}$ in the gaseous and aqueous phase are given in Table 5. These data show that, for both reactions C1 and C2, the formation of the complex $NiB_3O_4(OH)_{3(aq)}$ is thermodynamically possible. As expected, due to the electrostatic interaction, $Ni^{2+}$ is more reactant with $B_3O_3(OH)_4^-$ than with boric acid. $\Delta G_{R1}$ is higher than $\Delta G_{R2}$. The data illustrate also the solvent effect. Thus the energies in the aqueous phase are approximately five times lower than in the gaseous phase; this is explained by the fact that the interactions between the cation, $H_3BO_3$ and $B_3O_3(OH)_4^-$ are screened by the dielectric effects and the short range interactions between the solute and the solvent molecules taken into account formally in the SMD method.

The changes are more drastic when the nickel’s first sphere of hydration is treated explicitly in the simulations; then, $\Delta G_{R1}$ is three times lower than the value obtained without the hydration sphere and the reaction becomes exothermic. The trend in the results is in agreement with the one shown by the experimental data; the $R_1$ reaction’s energy is lower than that of $R_2$. However, as in references [32, 33], the results show the strong dependence of the reaction energies on the chosen thermodynamic cycle, though
<figure><img src="image_11.png" /><figcaption>**Fig. 7** Optimized molecules and complexes at the B3LYP level. **a** Ni²⁺·6H₂O, **b** H₃BO₃, **c** B₃O₃(OH)₄⁻, **d** NiB₃O₄(OH)₃ without Ni hydration sphere and **e** NiB₃O₄(OH)₃ with the Ni hydration sphere. The blue, gray, red and green atoms are H, B, O and Ni, respectively (Color figure online)</figcaption></figure>

<table><caption>**Table 5** Gibbs energy changes (in kJ·mol⁻¹) of the complex NiB₃O₄(OH)₃, ΔG<sub>R1</sub> and ΔG<sub>R2</sub> for the reaction C1 according to the relations 17 and 18 and for the reaction C2 according to the relations 21 and 22, respectively</caption><thead><tr><th></th><th>Approx.</th><th></th><th>ΔG<sub>R1</sub></th><th>ΔG<sub>R2</sub></th></tr></thead><tbody><tr><td rowspan="4">Ni²⁺</td><td rowspan="2">Gaseous phase</td><td>C1</td><td>-2074</td><td>-1741</td></tr><tr><td>C2</td><td>-1677</td><td>-936</td></tr><tr><td rowspan="2">Aqueous phase</td><td>C1</td><td>-376</td><td>-296</td></tr><tr><td>C2</td><td>-339</td><td>-222</td></tr><tr><td rowspan="4">Ni²⁺·6H₂O</td><td rowspan="2">Gaseous phase</td><td>C1</td><td>-997</td><td>-664</td></tr><tr><td>C2</td><td>-595</td><td>141</td></tr><tr><td rowspan="2">Aqueous phase</td><td>C1</td><td>25</td><td>106</td></tr><tr><td>C2</td><td>61</td><td>179</td></tr><tr><td>Experimental (this work)</td><td></td><td></td><td>24</td><td>65</td></tr></tbody></table>

Two types of nickel's environment are explored: without (Ni²⁺) and with its first sphere of hydration (Ni²⁺·6H₂O). Both reaction are treated in the gaseous and in the aqueous (with the solvent modeled with the SMD method) approximations. The experimental data are given for comparison
the semi-empirical treatment of H⁺ in C1 gives nearest results or results in best agreement with experiment, while the discrepancies with the experiment may be attributed to the complex nature of the reaction. Tossell showed that the precision on the pKₐ determination of the boric acid dissociation depends on the water molecules interacting with B(OH)₃(aq) and the formation of B(OH)₄⁻ [32]. Our approach may need to take into account these types of contributions. Anyway, taking into account the solvent effects, the most thermodynamically favorable complex formation reaction remains R₁. This result tends to show that the nickel cation in a solute should be more attracted by the existing B₃O₃(OH)₄⁻ than with boric acid.

## 6 Conclusion

By modeling previous experimental data of Shchigol [11], we were able to highlight which boron species could be a relevant ligand to complex nickel aqueous species, taking into account the speciation of boron and the chemical conditions. A complexation reaction was proposed where nickel cations react with the triborate, which is the most abundant polyborate species in the pH range between 6 and 10. Furthermore, the experimental pH monitoring data are consistent with the formation of this complex. It appears that the pH is lowered when complexation occurs, starting at 25 °C and to higher extent when temperature is increased to 70 °C. The equilibrium constants were determined at the three temperatures investigated in this study, fitted as a function of temperature and a set of thermodynamics data were deduced. Because we did not succeed in experimentally characterizing the complex by RMN or Raman methods, first principle calculations were used to study its stability. Results of the computations, taking into account the solvent, show that the complexation reaction involving the triborate species is thermodynamically favorable. For all those reasons, we consider the existence of this complex as valid. Nevertheless, it is possible that other complexes between Ni and borate may exist, in particular at lower concentration of boron when polyborates are minor, or with other polyborates and should be further investigated by using, for instance, EXAFS, XANES or UV spectroscopy.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## Appendix

### Atomic Gaussian Basis Sets Used for H, B, O and Ni

All-electron and pseudopotential basis sets have been used for H, B, O and Ni. The B and O all-electron basis sets are the same as used in references [24] for B and [25] for O; for H, the basis set is original. They are contractions of—5s-111sp-1p*, 6s-311sp-1d* and 8s-411sp-1d*—GTFs for H, B, and O, respectively. For pseudopotential, the Hay–Wadt small-core pseudopotentials [26–28] have been adopted for Ni; the associated basis set are a contraction of—3111sp-311d—GTFs. The exponents and contraction's coefficients of the pseudopotential GTFs, as well as for the all-electron basis sets, of the full set of the
<table><caption>**Table 6** Exponents and coefficients of the contracted Gaussian basis sets adopted in the present study for H, B, O and Ni. For Ni, the basis set is used in conjunction with the Hay–Wadt pseudopotential [26–28]</caption><thead><tr><th rowspan="2">Atom</th><th rowspan="2">Shell</th><th rowspan="2">Expt.</th><th colspan="2">Coeff.</th></tr><tr><th>s(d)</th><th>p</th></tr></thead><tbody><tr><td rowspan="10">H</td><td rowspan="5">s</td><td>157.418362</td><td>0.001441</td><td></td></tr><tr><td>92.942487</td><td>0.006237</td><td></td></tr><tr><td>17.727431</td><td>0.027257</td><td></td></tr><tr><td>9.844505</td><td>0.027295</td><td></td></tr><tr><td>3.516237</td><td>0.31327</td><td></td></tr><tr><td>s</td><td>0.904270</td><td>1.</td><td></td></tr><tr><td>s</td><td>0.295470</td><td>1.</td><td></td></tr><tr><td>s</td><td>0.105480</td><td>1.</td><td></td></tr><tr><td>p</td><td>1.074718</td><td>1.</td><td></td></tr><tr><td rowspan="5">B</td><td rowspan="2">sp</td><td>2.749437</td><td>-0.546313</td><td>0.171714</td></tr><tr><td>0.589916</td><td>1.133410</td><td>0.889293</td></tr><tr><td>sp</td><td>0.360657</td><td>1.</td><td>1.</td></tr><tr><td>sp</td><td>0.261831</td><td>1.</td><td>1.</td></tr><tr><td>d</td><td>0.684011</td><td>1.</td><td></td></tr><tr><td rowspan="3">O</td><td rowspan="2">sp</td><td>0.475975</td><td>1.</td><td>1.</td></tr><tr><td>0.161450</td><td>1.</td><td>1.</td></tr><tr><td>d</td><td>0.876952</td><td>1.</td><td></td></tr><tr><td rowspan="10">Ni</td><td rowspan="3">sp</td><td>25.405052</td><td>0.004128</td><td>-0.043494</td></tr><tr><td>7.273189</td><td>-0.608148</td><td>-0.115638</td></tr><tr><td>4.162815</td><td>0.392273</td><td>0.491603</td></tr><tr><td>sp</td><td>1.633718</td><td>1.</td><td>1.</td></tr><tr><td>sp</td><td>0.660162</td><td>1.</td><td>1.</td></tr><tr><td>sp</td><td>0.119693</td><td>1.</td><td>1.</td></tr><tr><td>d</td><td>50.098640</td><td>0.046157</td><td></td></tr><tr><td></td><td>14.382472</td><td>0.248248</td><td></td></tr><tr><td></td><td>4.905609</td><td>0.630232</td><td></td></tr><tr><td rowspan="2">For the all-electron basis sets of B and O, only the most diffuse GTFs are given (see Ref. [39] for a complete set of data)</td><td>d</td><td>1.672448</td><td>1.</td><td></td></tr><tr><td>d</td><td>0.503344</td><td>1.</td><td></td></tr></tbody></table>

H basis, of the B-four and O-three outer *sp* and *d* polarization's shell have been optimized using an energy criterium and are reported in Table 6. The basis set optimization was carried out using the LoptCGscript [37], which performs numerical gradient optimizations based on the conjugate gradient method [38].

## References

1. Pressurized Water Reactor Primary Water Chemistry Guidelines. vol. 1, revision 6, EPRI, Palo Alto, CA: 1014986 (2007)

2. Mesmer, R.E., Baes Jr., C.F., Sweeton, F.H.: Acidity measurements at elevated temperatures VI. Boric acid equilibriums. Inorg. Chem. **11**, 537–543 (1972)

3. Palmer, D.A., Bénézeth, P., Wesolowski, D.J.: Boric acid hydrolysis: a new look at the available data. Powerpl. Chem. **2**, 261–264 (2000)
4. Tremaine, P.R., Leblanc, J.C.: The solubility of nickel oxide and hydrolysis of Ni²⁺ in water to 573 K. J. Chem. Thermodyn. **12**, 521–538 (1980)

5. Palmer, D.A., Bénézeth, P., Xiao, C., Wesolowski, D.J., Anovitz, L.M.: Solubility measurements of crystalline NiO in aqueous solution as a function of temperature and pH. J. Solution Chem. **40**, 680–702 (2011)

6. Palmer, D.A., Gamsjäger, H.: Solubility measurements of crystalline beta-Ni(OH)₂ in aqueous solution as a function of temperature and pH. J. Coord. Chem. **63**, 2888–2908 (2010)

7. Ziemniak, S.E., Goyette, M.A.: Nickel(II) oxide solubility and phase stability in high temperature aqueous solutions. J. Solution Chem. **33**, 1135–1159 (2004)

8. You, D., Lovera, P., Plancque, G.: New Data for Thermodynamic and Kinetic Behaviour of Nickel Phases in PWR Physicochemical Conditions. Presented at the Nuclear Plant Chemistry, Sapporo (2014)

9. Palmer, D.A., Bénézeth, P., Wesolowski, D.J., Anovitz, L.M.: Impact of Nickel Oxide Solubility on PWR Fuel Deposit Chemistry. EPRI, Palo Alto (2002)

10. Wesolowski, D.J., Palmer, D.A., Begun, G.M.: Complexation of aluminate anion by bis–tris in aqueous media at 25–50 °C. J. Solution Chem. **19**, 159–173 (1990)

11. Shchigol, M.B.: Properties of cobalt and nickel borates. Russ. J. Inorg. Chem. **6**, 1361–1366 (1961)

12. Gamsjäger, H., Mompean, F.J.: Chemical Thermodynamics of Nickel. Elsevier Science and Technology, Amsterdam (2005)

13. Mukherjee, G.N.: Borate complexes in solution. Part-II. Mixed ligand complex formation of cobalt-, nickel- and zinc(II) with boric acid and some (N, N) bidentate ligands. J. Indian Chem. Soc. **79**, 45–47 (2002)

14. Bousher, A.: Unidentate complexes involving borate. J. Coord. Chem. **34**, 1–11 (1995)

15. Van der Lee, J.: Thermodynamic and Mathematical Concepts of CHESS (1998)

16. Ingri, N., Dahlen, J., Buchardt, O., Kvande, P.C., Meisingseth, E.: Equilibrium studies of polyanions. 11. Polyborates in 3.0 M Na(Br), 3.0 M Li(Br), and 3.0 M K(Br), a comparison with data obtained in 3.0 M Na(ClO₄). Acta Chem. Scand. **17**, 581–589 (1963)

17. Ingri, N.: Equilibrium studies of polyanions. II. Polyborates in NaClO₄ medium. Acta Chem. Scand. **11**, 1034–1058 (1957)

18. Zhou, Y., Fang, C., Fang, Y., Zhu, F.: Polyborates in aqueous borate solution: a Raman and DFT theory investigation. Spectrochim. Acta. **A** **83**, 82–87 (2011)

19. Momii, R.K., Nachtrieb, N.H.: Nuclear magnetic resonance study of borate-polyborate equilibria in aqueous solution. Inorg. Chem. **6**, 1189–1192 (1967)

20. Valiev, M., Bylaska, E.J., Govind, N., Kowalski, K., Straatsma, T.P., van Dam, H.J.J., Wang, D., Nieplocha, J., Apra, E., Windus, T.L., de Jong, W.A.: NWChem: a comprehensive and scalable open-source solution for large scale molecular simulations. Comput. Phys. Commun. **181**, 1477–1489 (2010)

21. Becke, A.D.: Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. **98**, 5648–5652 (1993)

22. Lee, C., Yang, W., Parr, R.G.: Development of the Colle–Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B **37**, 785–789 (1988)

23. Marenich, A.V., Cramer, C.J., Truhlar, D.G.: Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J. Phys. Chem. B **113**, 6378–6396 (2009)

24. Orlando, R., Dovesi, R., Roetti, C.: Ab-initio Hartree–Fock calculations for periodic compounds—Application to semiconductors. J. Phys. Cond. Matter **2**, 7769–7789 (1990)

25. Sophia, G., Baranek, Ph, Sarrazin, C., Rérat, M., Dovesi, R.: First principles study of the mechanisms of the pressure-induced dielectric anomalies in ferroelectric Perovskites. Ph. Trans. **86**, 1069–1084 (2013)

26. Hay, P.J., Wadt, W.R.: Ab initio effective core potentials for molecular calculations. Potentials for K to Au including the outermost orbitals. J. Chem. Phys. **82**, 299–310 (1985)

27. Hay, P.J., Wadt, W.R.: Ab initio effective core potentials for molecular calculations. Potentials for main group elements Na to Bi. J. Chem. Phys. **82**, 284–298 (1985)

28. Hay, P.J., Wadt, W.R.: Ab initio effective core potentials for molecular calculations. Potentials for transition metal atoms Sc to Hg. J. Chem. Phys. **82**, 270–283 (1985)

29. Wesolowski, D.J., Ziemniak, S.E., Anovitz, L.M., Machesky, M.L., Bénézeth, P., Palmer, D.A.: Solubility and surface adsorption characteristics of metal oxides. In: Aqueous Systems at Elevated Temperatures and Pressures: Physical Chemistry in Water, Steam and Hydrothermal Solutions. Elsevier Ltd, Amsterdam (2004)

30. Gu, Y., Gammons, C.H., Bloom, M.S.: A one-term extrapolation method for estimating equilibrium constants of aqueous reactions at elevated temperatures. Geochim. Cosmochim. Acta **58**, 3545–3560 (1994)
31. Mesmer, R.E., Marshall, W.L., Palmer, D.A., Simonson, J.M., Holmes, H.F.: Thermodynamics of aqueous association and ionization reactions at high temperatures and pressures. J. Solution Chem. **17**, 699–718 (1988)

32. Tossell, J.A.: Boric acid, “carbonic” acid, and N-containing oxyacids in aqueous solution: Ab initio studies of structure, pKₐ, NMR shifts, and isotopic fractionations. Geochim. Cosmochim. Acta **69**, 5647–5658 (2005)

33. Liptak, M.D., Shields, G.C.: Experimentation with different thermodynamic cycles on carboxylic acids using complete basis set and Gaussian-n models combined with CPCM continuum solvation methods. Int. J. Quantum Chem. **85**, 727–741 (2001)

34. Palmer, A.D., Fernandez-Prini, R., Harvey, A.H.: Aqueous Systems at Elevated Temperatures and Pressures—Physical Chemistry in Water, Steam and Hydrothermal Solutions. Elsevier Academic Press, Amsterdam (2004)

35. Hermet, P., Gourrier, L., Bantignies, J.-L., Ravot, D., Michel, T., Deabate, S., Boulet, P., Henn, F.: Dielectric, magnetic, and phonon properties of nickel hydroxide. Phys. Rev. B. **84**, 235211 (2011)

36. de P.R. Moreira, I., Illas, F., Martin, R.L.: Effect of Fock exchange on the electronic structure and magnetic coupling in NiO. Phys. Rev. B. **65**, 152 (2002)

37. Zicovich-Wilson, C.: LoptCG (Shell Procedure for Numerical Gradient Optimization). Instituto de Technologia Quimica Valencia, Valencia (2006)

38. Hestenes, M., Stiefel, E.: Methods of conjugate gradients for solving linear systems. J. Res. Natl. Bur. Stand. **49**, 409–436 (1952)

39. Variational basis sets are available at the following web site. http://www.crystal.unito.it/basis-sets.php.
Accessed 19 Nov 2016