2.1. Models
To model the hydrolysis reaction, we have considered a small network formed by two B atoms, connected by a bridging oxygen atom as the reagent (see reagent in
Figure 1). The overall starting coordination number of each boron atom is three, and the non-bridging oxygen atoms are saturated with hydrogens to complete their coordination and to produce a neutral compound. We carried out attempts to simulate it without saturating the compound with hydrogens but the excessive negative charge on oxygen atoms led to unrealistic proton transfer from the water molecules. The products of the hydrolysis reaction were, instead, two B(OH)
n compounds (with
n equal to either 3 or 4).
In addition to these compounds, we introduced explicit water molecules to model a pH 7 environment, and H3O+ and OH− ions to model acidic and basic conditions. For the neutral, acidic, and basic hydrolysis pathways, we also considered systems with two additional water molecules, as they are enough to create a ring-like transition state that can alter the energy barriers.
All the possible reactions pathways involved at least one transition state with a tetra-coordinated boron atom.
2.2. Assessment of the Basis Set and Ab Initio Levels of Theory
We considered the stationary points of the reaction pathway connecting a B-O-B containing unit, namely the B
2O
5H
4 compound, and two separated orthoboric acid molecules, in an environment without a single water molecule, as a reference to check the reliability of the basis set employed as well as the MP2 level of theory. In fact, this was the simplest model studied here. The geometry of the stationary point is illustrated in
Figure 1 and the model itself is discussed in more detail in
Section 2.3.
Figure 1.
Scheme of reaction pathway for the hydrolysis mechanism at neutral pH with a single water molecule. ξ represents the reaction coordinate. The energies are not to scale; see
Table 1 for a complete report on the energy values.
Figure 1.
Scheme of reaction pathway for the hydrolysis mechanism at neutral pH with a single water molecule. ξ represents the reaction coordinate. The energies are not to scale; see
Table 1 for a complete report on the energy values.
The stationary point was first found at the MP2/6-311++G(d,p) level of theory combined with PCM.
We performed single-point CCSD(T) [
21]/Apr-cc-pVTZ calculations on each of these geometries, combined with PCM, to check that the MP2-computed energy barriers were reliable in comparison with the CCSD(T) level of theory. This was performed as Gaussian 16 currently does not allow geometry optimization at the CCSD(T) level of theory.
The computed energy differences are reported in
Table 1, and suggest that CCSD(T) single points reproduce the results of MP2 optimizations with errors < 0.01 eV for the transition state barrier, and just slightly larger for the energy of the products.
Table 1.
Energies for the stationary points of the reaction pathway at neutral pH with 1 water molecule computed with MP2 and CCSD(T) methods and different basis sets. The energies have been rescaled to have the energy of the reagents put to zero. The energy of the TS coincides with .
Table 1.
Energies for the stationary points of the reaction pathway at neutral pH with 1 water molecule computed with MP2 and CCSD(T) methods and different basis sets. The energies have been rescaled to have the energy of the reagents put to zero. The energy of the TS coincides with .
Level of Theory | Stationary Point | Energy/eV |
---|
Optimized MP2/6-311++G(d,p) + PCM | Reagents | 0 |
TS | 2.030 |
Product | −0.015 |
Single-point CCSD(T)/Apr-cc-pVTZ + PCM | Reagents | 0 |
TS | 2.038 |
Products | −0.026 |
Optimized MP2/6-311++G(2d,2p) + PCM | Reagents | 0 |
TS | 2.026 |
Products | −0.022 |
Optimized MP2/6-31+G(d) + PCM | Reagents | 0 |
TS | 2.042 |
Products | −0.034 |
Optimized MP2/6-21G + PCM | Reagents | 0 |
TS | 2.373 |
Products | −0.009 |
Optimized MP2/LanL2DZ + PCM | Reagents | 0 |
TS | 2.118 |
Products | −0.011 |
Also, we checked the effect of the basis set, in particular, using the 6-21G, 6-31+G*, 6-311++G(d,p), and 6-311++G(2d,2p) basis sets and the combined basis-set-pseudopotential LanL2DZ [
22], optimizing the stationary points at the MP2 + PCM level of theory, as also reported in
Table 1.
The energy differences remained well below the 0.01 eV threshold for the 6-311++G(d,p) and 6-311++G(2d,2p) cases, and the B-Obridge distances of the transition state changed from 1.613 Å (tetrahedral B) and 1.414 to 1.616 Å (tetrahedral B) and 1.399 Å, respectively. Overall, expanding the basis set seemed to only marginally affect the energies and very negligibly affect the structure of the models we adopted here. On the contrary, using smaller basis sets led to increased differences in the energy barrier, which were particularly evident in the case of the 6-21G basis set (>0.1 eV).
Overall, we considered the 6-311++G** as a good compromise between accuracy, computational burden, and availability in the many possible quantum chemistry codes.
2.3. Neutral pH System
To study the hydrolysis mechanism in a neutral aqueous environment, we found and considered two possible different pathways, one including a single water molecule (1 W) and another one (3 W) exploiting the presence of three water molecules.
The structures associated with the pathway 1 W are reported in
Figure 1 of
Section 2.2, as we used this model (the simplest one) to also assess the validity of the ab initio and basis set levels of theory.
Following pathway 1 W, we found a transition state (TS) resulting from the nucleophilic attack of the O atom of water on a boron atom, and at the same time, the formation of a bond between a hydrogen atom of water and the bridging oxygen atom of the B-O-B network. This mechanism thus resulted in a very unstable TS (2.03 eV at the MP2 level of theory) constituted by a 4-term ring (B-O
water-H
water-O
bridge). Indeed, the
energies, i.e., the difference between the computed energy of the transition state and the energy of the reagents (see
Figure 1), were the highest that we found for any mechanism, confirming the high unlikely and unstable nature of this TS and, consequently, of this reaction pathway. The actual energies calculated with different functionals are reported in
Table 2.
Instead, models using three water molecules gave rise to a TS with a 7-term ring (see
Figure 2), resulting in much lower transition energies (0.478 eV at the MP2 level), as shown in
Table 3.
It must be pointed out that while the reaction model with just one water molecule, the DFT methods yielded reaction barriers about 0.5 eV lower than the MP2 result, in the model using three water molecules, this difference was significantly quenched. Indeed, ωB97XD approximated the MP2 result up to 0.01 eV, and B3LYP and PBE0 showed differences of about 0.1 and 0.2 eV, respectively. Also, some functional might yield an energy for products higher than that of the reactants, about hundredths of an eV.
2.4. Acidic pH System
In presence of an H3O+ cation, to mimic acidic conditions, the energy landscape of the reaction changes significantly. In particular, in the neutral environment, we were able to locate a single transition state, involving both the protonation of the bridging oxygen and the hydroxylation of a B atom; in acidic conditions, we observed these two processes occurring in sequence, in two different stationary points.
Namely, we consistently observed the protonation of the Obridge atom in the transition state TS1, which was followed by an intermediate state INT where a water molecule bound to a boron atom, thus making the coordination of the latter change from trigonal to tetrahedral. This was followed by a second transition state that we indicated as TS2, where the bond between the bridging oxygen and the tetrahedrally coordinated boron atom was weakened and lengthened, until it eventually broke to give rise to the products.
It is important to point out that we found this mechanism occurring both in the model with just one H
3O
+ cation (named 1 W, see
Figure 3), and in that with one H
3O
+ cation and two water molecules (named 3 W, see
Figure 4), at all levels of theory.
The presence of three stationary points between the reagents and products (two transition states and one intermediate) in the case of the acidic mechanism, instead of just one (a single, concerted transition state) as in the case of the mechanism at neutral conditions, dramatically lowers the energy barrier for the 1 W model, as reported in
Table 4.
Also, it must be highlighted that for the 1 W model, the first energy barrier,
, defined as the energy difference between the reagents and the TS
1, was significantly higher than the second one,
, defined as the energy difference between the intermediate (INT) and the TS
2. A graphical representation of both
and
is shown in
Figure 3. This means that the rate-determining step in this case is the protonation of the bridging oxygen. The actual energies found at the various levels of theory for this model are reported in
Table 4.
The situation is more complicated in the 3 W model due to the increased complexity of the system. Indeed, at the DFT level of theory, the specific choice of the functional can even yield a
negative energy for either
(M06HF functional) or
(CAM-B3LYP, BLYP, HSE06, PBE, and M06 functionals); this occurs since the TS2 has a negative vibrational frequency, but it still had a lower energy with respect to the intermediate INT. A graphical representation of both
and
is shown in
Figure 4, while the actual energies are reported in
Table 5.
We investigated the reasons behind this behavior, and found them to be dependent on these functional giving a higher hydrogen bond stabilization energy, which significantly lowered the energy of the TS2.
This also explains why the 1 W model, which employs just a single H3O+ cation, did not exhibit this behavior, since in that case, the number of hydrogen bonds potentially skewing the energy budget was smaller.
In fact, using these functionals to perform single-point calculations on the stationary points found and optimized at the MP2 level produced a behavior that is largely correct but, as expected, they cannot recognize transition states or intermediates based on their vibrational frequencies. We want to underline here that B3LYP, ωB97XD, and PBE0 all gave energies and energy differences close to those found at the MP2 level, with differences of less than 0.1 eV.
Overall, the energy barriers at the MP2 level and using B3LYP, ωB97XD, and PBE0 functionals yielded similar results in the 1 W and 3 W models; this can be explained by the lack of a concerted transition state in acidic conditions that (energetically) benefits from having a larger amount of water molecules explicitly included into the model. Therefore, different from the neutral simulations, in this case, the 1 W model seemed to be reliable for obtaining accurate energy barriers. However, for the functionals that poorly reproduce H-bonds, the 3 W model could indeed lead to erratic results.
2.5. Basic pH System
We simulated alkaline conditions by using a single hydroxyl anion in our system, combined with 0 (model 1 W) or 2 (model 3 W) neutral water molecules. The hydroxyl anion attacks a boron atom, turning its coordination from trigonal to tetrahedral, giving rise to a single TS and, thus, to a rather simple energy pathway (see
Figure 5 and
Figure 6).
In particular, it must be noted that the reaction barrier became extremely small even for the 1 W model, on the order of about 0.1 eV at the MP2, B3LYP, PBE0, and ωB97XD levels of theory, as reported in
Table 6. In fact, the reaction barrier was so small that we had issues optimizing the reagents, as they very easily evolved into the transition state, which was very close in energy. Also, the reaction evolved without the need for the protonation of the bridging oxygen, which is different from the mechanisms in the neutral and acidic conditions.
In the more complex 3 W model, we observed that the reaction barriers decreased (see
Table 7) while maintaining a rather similar transition geometry for the B
2O
5H
4 and OH
− fragments but with a ring-like disposition for the other two additional water molecules. This decrease in the reaction barriers occurred at all levels of theory.