Ouissam
El Bakouri
a,
Miquel
Solà
*a and
Jordi
Poater
*bcd
aInstitut de Química Computacional i Catàlisi (IQCC) and Departament de Química, Universitat de Girona, Campus Montilivi, 17071 Girona, Catalonia, Spain. E-mail: miquel.sola@udg.edu
bDepartament de Química Inorgànica i Orgànica & Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Catalonia, Spain. E-mail: jordi.poater@gmail.com
cDepartment of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Vrije Universiteit Amsterdam, De Boeleaan 1083, NL-1081HV Amsterdam, The Netherlands
dInstitució Catalana de Recerca i Estudis Avançats (ICREA), Pg. Lluís Companys 23, 08010 Barcelona, Catalonia, Spain
First published on 11th March 2016
Despite the fact that B and Al belong to the same group 13 elements, the B62− cluster prefers the planar D2h geometry, whereas Al62− favours the Oh structure. In this work, we analyse the origin of the relative stability of D2h and Oh forms in these clusters by means of energy decomposition analysis based on the turn-upside-down approach. Our results show that what causes the different trends observed is the orbital interaction term, which combined with the electrostatic component do (Al62− and Ga62−) or do not (B62−) compensate the higher Pauli repulsion of the Oh form. Analysing the orbital interaction term in more detail, we find that the preference of B62− for the planar D2h form has to be attributed to two particular molecular orbital interactions. Our results are in line with a dominant delocalisation force in Al clusters and the preference for more localised bonding in B metal clusters. For mixed clusters, we have found that those with more than two B atoms prefer the planar structure for the same reasons as for B62−.
Among clusters, those made by group 13 atoms are particularly important.5 Both B and Al belong to the same group 13, and thus present a similar electronic structure, [He]2s22p1 and [Ne]3s23p1, respectively. However, when they form small clusters, B clusters adopt a planar conformation as the equilibrium structure;6–9 whereas Al clusters present a three-dimensional (3D) closed shape.10–13 The most relevant examples are B62− and Al62− clusters, which were obtained experimentally as lithium salts in the form of LiB6− and LiAl6−.14–16 B62− adopts a planar D2h geometry in its low-lying singlet state, whereas the Al62− cluster is octahedral. Both shapes of the metal clusters are kept when lithium salts are formed.
The chemical bonding of B62− and Al62− has been widely analysed in previous studies.14,17,18 In particular, Alexandrova et al.18 highlighted the fact that B62− is able to 2s–2p hybridize and to form 2-center–2-electron (2c–2e) B–B covalent localised bonds. On the other hand, 3s–3p hybridisation in the Al62− cluster is more difficult due to larger s–p energy separation, which hampers the formation of directional covalent Al–Al bonds.19 In this case, bonding comes from the combination of radial and tangential p-orbitals that result in extensive delocalisation.20 Indeed, the Al62− cluster displays octahedral aromaticity,14,21 whereas planar D2h B62− is considered σ- and π-antiaromatic.17,18,22,23 Thus, as pointed out by Alexandrova et al.,18,24–26 covalent and delocalised bonding shows opposite effects in determining the molecular structure of many clusters. Huynh and Alexandrova analysed the whole series BnAl6−n2− (n = 0–6), from B62− till Al62− by substituting one B by Al each time, concluding that covalent bonding is a resilient effect that governs the cluster shape more than delocalisation does. Indeed, the planar structure of B62− persists until n = 5, the reason being the strong tendency to form 2c–2e B–B bonds in case the cluster contains two or more B atoms.18 Similar results were reported by Fowler and Ugalde in larger clusters of group 13. In particular, these authors found that B13+ prefers a planar conformation27 in contrast to Al13−,28 which adopts an icosahedral geometry. Interestingly, in closo boranes and substituted related species, like B6H62− or B12I122−, the delocalised 3D structure is preferred. However, successive stripping of iodine in B12I122− leads to a B12 planar structure with some localised 2c–2e B–B bonds.29,30 Similarly, for B6Hn− clusters, the clusters are planar for n ≤ 3 and become tridimensional for n ≥ 4.31
As can be seen in Scheme 1, both 2D D2h planar and 3D Oh geometries for X62− (X = B, Al) can be obtained joining the same two X3− cluster fragments.14,17 Therefore, X62− species in D2h and Oh geometries are particularly suitable for an energy decomposition analysis (EDA)32–35 based on the turn-upside-down approach.36–39 In this approach, two different isomers are formed from the same fragments and the bonding energy is decomposed into different physically meaningful components using an EDA. Differences in the energy components explain the reasons for the higher stability of the most stable isomer. For instance, using this method we provided an explanation of why the cubic isomer of Td geometry is more stable than the ring structure with D4h symmetry for (MX)4 tetramers (X = H, F, Cl, Br, and I) if M is an alkalimetal and the other way round if M belongs to group 11 transition metals.38 Therefore, the application of this type of analysis to B62− and Al62− clusters will disclose the factors that make the planar D2h structure more stable for boron and the octahedral one for aluminium. As said before, boron clusters favour localised covalent bonds whereas aluminium clusters prefer a more delocalised bonding. With the present analysis, we aim to provide a more detailed picture of the reasons for the observed differences. The analysis will be first applied to the above referred B62− and Al62− clusters, and then further complemented with Ga62−. Finally, X2Y42− and X3Y32− (X, Y = B, Al, Ga) mixed clusters in their distorted D2h planar and 3D D4h geometries will also be discussed.
Relative energies between the planar and 3D species were also calculated using the Gaussian 09 program45 at the coupled cluster level46 with single and double excitation (CCSD)47 and with triple excitation treated perturbatively (CCSD(T))48 using Dunning's correlation consistent augmented triple-ζ (aug-cc-pVTZ)49,50 at optimised BLYP-D3(BJ)/TZ2P molecular geometries.
The bonding energy corresponding to the formation of X62− for both D2h and Oh symmetries from two anionic quintet tetraradicals, fragment 1 (αααα) + fragment 2 (ββββ) (see Scheme 1), is made up of two major components (eqn (1)):
ΔE = ΔEdist + ΔEint | (1) |
ΔEint = ΔVelstat + ΔEPauli + ΔEoi + ΔEdisp | (2) |
In bond-energy decomposition,51–53 open-shell fragments were treated with spin-unrestricted formalism but, for technical reasons, spin-polarisation was not included. This error causes the studied bond to become in the order of a few kcal mol−1 too strong. To facilitate a straightforward comparison, the EDA results were scaled to match exactly the regular bond energies (the correction factor is consistently in the range 0.97–0.98 in all model systems and does therefore not affect trends). A similar scheme based on the same EDA approach was used by Frenking and coworkers54,55 and by some of us36,37,56 to estimate the strength of π-cyclic conjugation in typical (anti)aromatic organic compounds and in metallabenzenes and metallacyclopentadienes.
Let us mention here that, as already mentioned in the introduction, some of the analysed metal clusters exist experimentally as lithium salts.14–16 On the other hand, these dianionic systems are unstable against the ejection of an electron. However, their molecular and electronic structure is very similar to that of their corresponding lithium salts, which justifies the analysis of the chemical bonding of these doubly charged systems, as it is not affected by the presence of a lithium cation.
Finally, the metalloaromaticity57 of these clusters was evaluated at the BLYP/aug-cc-pVDZ level of theory with the optimized BLYP-D3(BJ)/TZ2P geometries by means of multicentre electron sharing indices (MCIs).58–60 MCIs provide a measure of electron sharing among the atoms considered,59 in the present case the six atoms that form each of the clusters studied. MCI values have been calculated using the ESI-3D program.61,62
Fig. 1 Geometries of X62− metal clusters analysed with D2h and Oh symmetries. Distances in Å and angles in degrees. |
Table 1 encloses the energy differences between Oh and D2h clusters. For B62−D2h symmetry is more stable than Oh by 67.5 kcal mol−1, the latter not being a minimum.18 Meanwhile the opposite trend is obtained in the other two metal clusters, for which Oh is lower in energy by 15.8 (Al62−) and 9.3 kcal mol−1 (Ga62−) than D2h structures. These trends are confirmed by higher level CCSD(T)/aug-cc-pVTZ single point energy calculations at the same BLYP-D3(BJ)/TZ2P geometries (values also enclosed in Table 1). The relative energies of B62−, Al62−, and Ga62− between Oh and D2h symmetries are now −38.7, +44.8 and +46.6 kcal mol−1, respectively. CCSD(T) values systematically favour Oh as compared to D2h structures by about 20–30 kcal mol−1. However, the qualitative picture remains the same.
Clusters | BLYP-D3(BJ)/TZ2Pa | CCSD(T)/aug-cc-pVTZb | MCIc | ||||
---|---|---|---|---|---|---|---|
O h | D 2h | O h | D 2h | O h | D 2h | ||
a B2Ga42− (D2h) has not been obtained because optimization breaks the symmetry; whereas B3Al32− and B3Ga32− (Oh) have not been obtained because the strength of the B3 unit causes the systems to be planar and to avoid a 3D geometry. b Single point energy calculations at BLYP-D3(BJ)/TZ2P geometries. c MCI calculated at the BLYP/aug-cc-pVDZ level of theory with the BLYP-D3(BJ)/TZ2P optimized geometries. d Local minima. e One imaginary frequency. f One small imaginary frequency due to numerical integration problems. g Two imaginary frequencies. h Three imaginary frequencies. | |||||||
X62− | B62− | 67.5e | 0.0d | 38.7 | 0.0 | 0.062 | −0.052 |
Al62− | 0.0d | 15.8e | 0.0 | 44.8 | 0.077 | 0.068 | |
Ga62− | 0.0f | 9.3f | 0.0 | 46.6 | 0.083 | 0.071 |
The aromaticity of these X62− metal clusters was evaluated by means of the MCI electronic criterion. The six-membered MCIs are enclosed in Table 1. In all cases, the Oh system is more aromatic than the D2h one, in agreement with the larger electronic delocalisation of the former, as discussed in the Introduction.21 MCI values confirm the octahedral aromaticity21 of Oh Al62− and the antiaromatic character of D2h B62−.17,18,22,23 Interestingly, MCI values point out the clear aromatic character of all 3D clusters that do not contain boron (MCI = 0.074–0.077); whereas mixed B2Al42−, Al2B42−, and Ga2B42−D4h clusters containing boron atoms are less aromatic (MCI = 0.032–0.047). For planar structures, there are basically two groups of clusters. First, the group formed by B62− and B2Al42− has eight valence electrons distributed in two π-MOs and two σ-MOs (vide infra). Therefore, having four π-electrons and four σ-electrons, they are σ- and π-antiaromatic species. Second, the group formed by Al62−, Ga62−, Al2B42−, Al2Ga42−, Ga2B42−, and Ga2B42− have eight valence electrons distributed in one π-MO and three σ-MOs (vide infra) and, therefore, they are σ- and π-aromatic species.
With the aim to obtain a deeper insight into the origin of 2D to 3D relative energies an energy decomposition analysis was performed, following the reaction presented in Scheme 1. As pointed out above, both systems can be constructed from two identical X3− anionic fragments, both in their quintet state in order to form the corresponding new bonds. Three of these bonds are of σ character, two tangential (σT) and one radial (σR), and one π character (see Fig. 2). It must be pointed out that, very recently, Mercero et al. have proven the multiconfigurational character of some of the lowest-lying electronic states of Al3−.19 In the case of the quintet state of Al3−, which is the fragment used in our calculations, the authors showed that the electronic configuration of the four valence electrons is also derived from the occupation of two σ-type tangential and one σ-type radial molecular orbitals arising from the 3px and 3py atomic orbitals, and one π-type orbital arising from the 3pz ones. This quintet state was found to be dominated by one-single configuration with a coefficient of 0.92 in the multiconfigurational wavefunction.19 Moreover, the energy difference between the ground state and the quintet state was almost the same when computed at DFT or at the MCSCF levels of theory.19 This seems to indicate that DFT methods give reasonable results for this quintet state. Finally, the T1 test64 applied to clusters collected in Table 1 was found to be always less than 0.045, thus indicating the relatively low multiconfigurational character of these species. It is commonly accepted that CCSD(T) produces acceptable results for T1 values as high as 0.055.65
The different terms of the EDA for B62−, Al62−, and Ga62− clusters are enclosed in Table 2. First we notice that the total bonding energies (ΔE) are much larger for B62− than for Al62− or Ga62−. For the former, ΔE are −100.2 (Oh) and −179.5 kcal mol−1 (D2h), whereas for the two latter are in between −19.0 and −38.1 kcal mol−1. This trend correlates with the shorter B–B bond lengths mentioned above. Table 2 also encloses the relative EDA energies between the two clusters. The B3− fragment taken from the B62− system in its D2h symmetry is the one that suffers the largest deformation, i.e. the largest change in geometry with respect to the fully relaxed B3− cluster in the quintet state (ΔEdist = 12.5 kcal mol−1), whereas the rest of the systems present small values of ΔEdist (0.0–1.7 kcal mol−1). However, differences in ΔE are not due to distortion energies (indeed ΔEdist values follow the opposite trend as ΔE), but to interaction energies (ΔEint).
B62− | Al62− | Ga62− | |||||||
---|---|---|---|---|---|---|---|---|---|
D 2h + D2h → D2h | O h + Oh → Oh | Δ(ΔE) | D2h + D2h → D2h | O h + Oh → Oh | Δ(ΔE) | D 2h + D2h → D2h | O h + Oh → Oh | Δ(ΔE) | |
ΔEint | −192.0 | −101.4 | −90.6 | −20.7 | −39.8 | 19.1 | −19.1 | −31.0 | 11.9 |
ΔEPauli | 533.5 | 735.3 | −201.8 | 225.7 | 348.0 | −122.3 | 269.6 | 384.5 | −114.9 |
ΔVelstat | −239.0 | −291.9 | 52.9 | −96.3 | −166.5 | 70.2 | −138.0 | −207.5 | 69.5 |
ΔEoi | −483.4 | −542.8 | 59.4 | −146.9 | −217.4 | 70.5 | −146.7 | −203.4 | 56.7 |
ΔEdisp | −3.2 | −2.1 | −1.1 | −3.2 | −3.9 | 0.7 | −4.0 | −4.7 | 0.6 |
ΔEdist | 12.5 | 1.3 | 11.2 | 0.0 | 1.7 | −1.7 | 0.1 | 1.4 | −1.3 |
ΔE | −179.5 | −100.2 | −79.3 | −20.7 | −38.1 | 17.4 | −19.0 | −29.6 | 10.6 |
Thus, we focus on the decomposition of ΔEint into ΔEPauli, ΔVelstat, ΔEoi, and ΔEdisp terms. As a general trend, in all three X62− clusters ΔEPauli is larger for the Oh than the D2h cluster (Δ(ΔEPauli) = −201.8, −122.3, and −114.9 kcal mol−1 for B62−, Al62−, and Ga62−, respectively), so making it less stable. The overlaps between doubly occupied MOs are larger in the more compact Oh structure that, consequently, has larger ΔEPauli. The larger difference in ΔEPauli between the Oh and D2h structures in the case of B62− as compared to Al62− and Ga62− is attributed to the particularly short B–B distances that increase the overlap between doubly occupied MOs of each B3− fragment. At the same time, the Oh form presents larger (more negative) electrostatic interactions (Δ(ΔVelstat) = 52.9, 70.2, and 69.5 kcal mol−1 for B62−, Al62−, and Ga62−, respectively). It is usually the case that higher destabilising Pauli repulsions go with larger stabilising electrostatic interactions. The reason has to be found in the fact that both interactions increase in the absolute value when electrons and nuclei are confined in a relatively small space. The electrostatic interaction together with orbital interaction (Δ(ΔEoi) = 59.4, 70.5, and 56.7 kcal mol−1 for B62−, Al62−, and Ga62−, respectively) terms favour the Oh structure. However, in the case of Oh B62−, Δ(ΔVelstat) and Δ(ΔEoi) cannot compensate Δ(ΔEPauli), which causes the D2h system to be the lowest in energy. The opposite occurs for Al62− and Ga62−. Finally, the dispersion term almost does not affect the relative energies, as the difference in dispersion is only in the order of ca. 1.0 kcal mol−1. Therefore, what causes the different trend observed for B62− on one side, and Al62− and Ga62− on the other side is basically the ΔEoi term, which combined with the ΔVelstat component does (Al62− and Ga62−) or does not (B62−) compensate the higher ΔEPauli of the Oh form.
The comparison of the MOs diagrams of B62− and Al62−, built from their X3− fragments, justify the trends of ΔEoi (see Fig. 2 and 3). Both D2h and Oh clusters are built from the same fragments; the only difference is that the two tangential fragHOMO(σT(b2)) and fragHOMO−1(σT(a1)) MOs of Al3− are degenerate when obtained from Al62− in its Oh geometry, whereas they are not when generated from the D2h system, although they still are very close in energy. As discussed from the EDA, Oh is more stable than D2h because of more stabilizing electrostatic and orbital interactions, which compensate its larger Pauli repulsion. Fig. 2 also encloses the overlaps for the interactions between the four SOMOs of the Al3− fragments to form the MOs of the metal clusters in both geometries. We take the Al3− fragments in their quintet states with three unpaired σ- and one unpaired π-electrons, all of them with spin α in one fragment and β in the other. A more negative ΔEoi in Oh Al62− is justified from the larger 〈SOMO|SOMO〉 overlaps, especially for t2gOHOMOah and OHOMObh (0.360 compared to 0.225 and 0.232 for b2uDHOMO2h and agDHOMO−12h, respectively). D2h only presents a larger overlap for the π fragment orbital (0.251 for b3uDHOMO−22h and 0.124 for t2gOHOMOch). Meanwhile both of them have almost the same overlap for the combination of the radial MO (σR) fragment (fragHOMO−2), with 〈SOMO|SOMO〉 = 0.298 and 0.301 for agDHOMO−32h and OHOMO−1h, respectively. Overall, the higher orbital interaction term of the Oh system can be explained by the larger 〈SOMO|SOMO〉 overlaps of two of the t2g delocalised molecular orbitals for this cluster (see Fig. 2). The energies of the occupied MOs of Al62− formed are higher than those of the Al3− SOMOs because we move from a monoanionic fragment to a dianionic molecule.
Now it is the turn to visualize the MOs of B62−. The fragments for B3− are the same as those for Al3− (see Fig. 3). However, the first difference appears in the MOs for B62− with D2h symmetry. In this case, it would be more reasonable to build the MOs of this molecule from two triplet (not quintet) B3− fragments. The reason is the different occupation of the MOs when compared to the D2h Al62− species. In D2h B62−, the HOMO corresponds to the antibonding π MO. To reach doubly occupied bonding (b3uDHOMO−42h) and antibonding (b2gDHOMO2h) π MOs, the π MO (fragHOMO−3) should be doubly occupied. Furthermore, the tangential σT(a1) fragHOMO does not participate in any occupied MO of this metal cluster and only generates virtual MOs. Consequently, MOs of B62− are better formed from two B3− fragments in their triplet state (see red electron in Fig. 3). On the other hand, B62− with Oh follows the same trend as Al62−, and in this case the same SOMOs in their quintet state are involved. At this point, it is worth mentioning that, as pointed out by Mercero et al., due to the strong multiconfigurational character of this species, one must be cautious with the electronic configuration, especially for the triplet state, as radial and tangential MOs are very close in energy.19
To make results comparable, Table 2 gathers the EDA of Oh and D2h B62− from two B3− fragments in their quintet states. Also in this case ΔEoi is more favourable for Oh than for D2h, however, at a lower extent when compared to Al62−. There are two main reasons for such a decrease of the strength of ΔEoi in Oh compared to D2h. First, and more importantly, because the DHOMO−22h formed presents a much larger 〈SOMO|SOMO〉 overlap than t2gOHOMO−1h (0.518 in the former vs. 0.338 in the latter). In particular, this DHOMO−22h MO contributes to the 2c–2e B–B localised bonds that are related to the larger covalent character of this structure. And second, because the π-interaction between the two π SOMO fragments is much larger in the case of D2h (0.225 vs. 0.059 for D2h and Oh, respectively). Nevertheless, these two more favourable orbital interactions are not enough to surpass the ΔEoi term of the Oh cluster. However, as compared to Al62−, for B62− the Δ(ΔEoi) term favours the Oh system to a less extent and cannot compensate the higher ΔEPauli term of the Oh form, thus making the planar geometry to be more stable in this case. This is related to the determinant force of the formed covalent bonding, involving more localised MOs than for Al62−. Such a larger covalent component in B62− is also supported by the covalent character of the interaction between the two fragments calculated as % covalency = (ΔEoi/(ΔEoi + ΔVelstat + ΔEdisp)) × 100. This formula results in B62−: 65–67% (Oh, D2h), Al62−: 56–60% (Oh, D2h), and Ga62−: 49–51% (Oh, D2h); thus confirming again the larger covalency found in B62−.
Finally, as done usually in the turn-upside-down approach,36–39,56,66,67 instead of building X62− in Oh symmetry from the corresponding two X3− fragments obtained from the Oh structure, we can build the Oh system from two X3− fragments extracted from the X62− cluster in D2h symmetry, and viceversa (see Tables S2–S4 in the ESI‡). The main conclusions remain unaltered and confirm that the D2h structures suffer a lower Pauli repulsion whereas those of Oh symmetry have more favourable electrostatic and orbital interactions. The interplay between the Pauli repulsion on the one hand and electrostatic and orbital interactions on the other determines the most favorable symmetry in each case.
Just to conclude this section, we must point out that the whole EDA and turn-upside-down analyses were performed with fragments in their quintet state. However, as we commented before this is not the most reasonable way to build B62− in D2h symmetry. Table S5 (ESI‡) contains the EDA for Oh and D2h B62− systems using B3− fragments in their triplet states. Results show that although the different terms are larger in the absolute value, the trends discussed above are not affected, and the D2h cluster is favoured mainly because of smaller Pauli repulsions.
Fig. 4 Geometries of mixed metal clusters analysed with planar and 3D geometries. Distances in Å and angles in degrees. |
The EDA was also performed for this series of six mixed metal clusters (see Table 3) with the aim to further understand the determinant force towards the most stable cluster. For the X2Y42− clusters, the EDA was carried out taken YXY− fragments in their quintet states. For Al3Ga32−, the fragments were Al3− and Ga3− in the quintet state too. For those systems for which the out-of-plane geometry is the most stable, the combination of more favourable electrostatic and orbital interactions, even though presenting larger Pauli repulsion, gives the explanation to the trend observed. This is the same behaviour already discussed above for both Al62− and Ga62− systems. On the other hand, when D2h symmetry is the cluster lower in energy, as for Al2B42− and Ga2B42− metal clusters, even though the D4h system presents more stable electrostatic interaction, now the orbital interactions in combination with less unfavourable Pauli repulsion favour the D2h symmetry. This latter behaviour differs from that of B62−, for which the orbital interactions also favour the Oh symmetry, thus making Pauli repulsion the determinant factor towards the preference for planar D2h B62−.
ΔEint | ΔEPauli | ΔVelstat | ΔEoi | ΔEdisp | ||
---|---|---|---|---|---|---|
B2Al42− | D 4h | −52.1 | 440.1 | −202.7 | −285.7 | −3.9 |
D 2h | −40.4 | 243.4 | −98.1 | −182.5 | −3.3 | |
ΔE | 11.7 | −196.7 | 104.6 | 103.2 | 0.6 | |
Al2B42− | D 4h | −75.1 | 584.0 | −251.7 | −404.1 | −3.3 |
D 2h | −139.6 | 556.6 | −238.6 | −454.5 | −3.3 | |
ΔE | −64.6 | −27.4 | 13.2 | −50.4 | 0.0 | |
Al2Ga42− | D 4h | −35.0 | 381.2 | −201.0 | −210.6 | −4.6 |
D 2h | −19.2 | 283.1 | −147.1 | −151.4 | −3.8 | |
ΔE | 15.8 | −98.1 | 53.8 | 59.3 | 0.8 | |
Ga2B42− | D 4h | −83.8 | 590.4 | −262.5 | −408.2 | −3.5 |
D 2h | −157.5 | 540.0 | −225.6 | −468.6 | −3.2 | |
ΔE | −73.7 | −50.4 | 36.8 | −60.4 | 0.3 | |
Ga2Al42− | D 4h | −38.4 | 370.1 | −188.0 | −216.1 | −4.3 |
D 2h | −20.6 | 218.0 | −90.8 | −144.4 | −3.6 | |
ΔE | 17.8 | −152.0 | 97.3 | 71.7 | 0.8 | |
Al3Ga32− | D 3h | −36.8 | 381.0 | −197.8 | −215.7 | −4.2 |
C 3v | −20.7 | 254.4 | −122.8 | −148.7 | −3.6 | |
ΔE | 16.1 | −126.6 | 75.0 | 67.0 | 0.6 |
Footnotes |
† This work is dedicated to Prof. Evert Jan Baerends as a proof of our admiration for his brilliant contributions to chemistry and of our gratitude for helping us to understand chemistry better. |
‡ Electronic supplementary information (ESI) available. See DOI: 10.1039/c6cp01109h |
This journal is © the Owner Societies 2016 |