Comparison of Model Systems (M+)n·[CrX63−] and M3CrX6 + 18MX Based on Quantum-Chemical Calculations (X: F, Cl)

Kremenetsky, Vyacheslav; Kuznetsov, Sergey

doi:https://doi.org/10.1155/2016/8324636

Journal of Chemistry

On this page

Abstract Introduction Results and Discussion Conclusions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2016 | Article ID 8324636 | https://doi.org/10.1155/2016/8324636

Comparison of Model Systems (M⁺)_n·[CrX₆³⁻] and M₃CrX₆ + 18MX Based on Quantum-Chemical Calculations (X: F, Cl)

Vyacheslav Kremenetsky¹and Sergey Kuznetsov¹

Academic Editor: Stefano Caporali

Received27 Sept 2016

Accepted28 Nov 2016

Published21 Dec 2016

Abstract

On the basis of quantum-chemical calculations the most stable particle compositions are estimated in such model systems as (M⁺·[CrCl₆] and M₃CrCl₆ + 18MCl (M = Na, K, and Cs). In all systems these particles are positively charged. For systems (M⁺·[CrCl₆], (M⁺·[CrF₆], M₃CrF₆ + 18MCl, M₃CrF₆ + 18MF, and M₃CrCl₆ + 18MCl (M = Na, K, and Cs) a number of energy parameters characterizing the state of the system before and after electron transfer are calculated. The results indicate the possibility of electron transfer from the cathode to the melt system, which is in the initial state. However, this possibility cannot be realized in systems where LUMOs (lowest unoccupied molecular orbitals) have purely ligand character. In this case, the preliminary deformation of a cationic shell of electroactive species is required; it transforms the initial system to the transition state. However, in all considered systems the search of the transition state should be carried close to the initial state . This greatly simplifies a problem and transforms it from a purely theoretical sphere to the field of practical tasks that do not require exceptional cost of computer time.

1. Introduction

Previously [1–4], quantum-chemical calculations of such model systems as (M⁺·[NbF₇], (M⁺·[NbCl₆], and (M⁺·[CrCl₆] (M = Na, K, Cs; ) were carried out. These systems contain a complex anion with a cation outer-sphere (OS) shell. It was found that compositions with an intermediate number of OS cations , where is the limiting number of OS cations bound to a given complex, are thermodynamically most stable. Such systems (designating them as the type I system) are often used for modeling of the charge transfer process in melts [5–9]. The activation energy of electron transfer is estimated here based on Marcus theory [10]. These systems are suitable for express-evaluation of the charge transfer parameters, because it does not require large expenditures of computer time. However, it is clear that such simple systems can have only limited application. In particular, the effect of the anion composition of the electrolyte cannot be incorporated in the type I systems.

To verify the main findings for model systems I quantum-chemical computations of some parameters have been carried out in extended model systems of type II M₃CrX₆ + 18MX (M = Na, K, and Cs; X: F, Cl).

Some results for the CrF₆-containing systems I and II were given previously [11–13]; for CrCl₆-containing systems I and II calculations were made in the nonrelativistic basis sets for M: Na and K in the main [2, 4, 14, 15]. Here we provide data for the systems (M⁺·[CrCl₆] and M₃CrCl₆ + 18MCl in a quasi-relativistic ECP basis for M: Na, K, and Cs and an additional analysis of all these systems.

2. Computational Methods

The quantum-chemical calculations were performed with the Firefly program package [16], partially based on the codes of the GAMESS (US) program [17], by the density functional theory methods (DFT/UHF). The spin-polarized version of the Kohn-Sham equations and the B3LYP hybrid functional were used. All calculations are made with the quasi-relativistic ECP basis set of Stuttgart/Dresden groups [18–20], more exactly, with Stuttgart RSC 1997 ECP for Cr, K, and Cs and Stuttgart RLC ECP for F, Cl, and Na. In addition Na-systems in some cases used CRENBL ECP basis set [20]. The energy values correspond to true minima of the potential energy surface (the imaginary frequencies are absent). In all complexes chromium is in a high-spin state: for Cr(III) and for Cr(II). Spin contamination can be ignored: for the Cr(III) and Cr(II) particles, the values were 3.76–3.82 and 6.01–6.03, respectively (3.75 and 6.00 in the ideal case).

3. Results and Discussion

3.1. Composition of the Most Stable Particles

In systems II calculation of the interaction energy of the second coordination sphere fragments (M⁺ with both the complex and the external environment of a given fragment can be made. Figures 1(a) and 1(b) show the plots of the energy of cationic OS shell formation in systems I (M⁺·[CrCl₆]³⁻ (Figure 1(a)) and II M₃CrCl₆ + 18MCl (Figure 1(b)) calculated by where , , and are the energies of systems I (or an equivalent fragment of systems II), free [CrCl₆]³⁻ complex, and free M⁺ cation, respectively. The energy value reflects the total effect of the interaction of OS cations M⁺ with the complex [CrCl₆]³⁻ and with each other.

(a)

(b)

Dependence of this type always has a minimum at some intermediate value. The existence of the minimum is mainly caused by an increase in the repulsion between OS cations as their number increases. The composition of the system at point is the most stable. In the case of systems I for all M (Figure 1(a)). The maximal number of OS cations retained by the chromium complex in systems I is 6.

Figure 1(b) shows the energy of the OS shell formation for analogous fragments of systems II (M⁺·[CrCl₆]³⁻. Note, in systems II, the OS cations are ranked by increasing the (Cr-M) distances (); that is, with increase of the number on Figure 1(b) the distance (Cr-) increases, too. The second coordination sphere of chromium in these systems contains 8-9 cations; however the minimum corresponds to equal to 4 (Na, K) or 5 (Cs) only. Figure 2 shows the examples of model structures for systems I and II.

(a)

(b)

Thus, no significant displacement of the minimum was observed during the transition from systems I to systems II. This allows considering data based on results for systems I as a valid initial estimate.

Unlike the systems I, for the systems II the interaction energy of the second coordination sphere (M with the complex and the rest of the system can be calculated directly. If the interaction energy of this shell with the complex exceeds the energy of its interaction with the outer environment, we can state that there is a dynamic equilibrium in the system responsible for the existence of rather stable complex species of definite composition [11]. Relevant data are presented below.

Figure 3 shows the energy of interaction of the (M⁺·[CrCl₆]³⁻ fragment with the outer environment (the rest of the systems II) as a function of the OS cations number. The values were calculated from the following equation:

Here the symbols , , and denote energies of the entire model system, its fragment (M_n)·[CrCl₆]³⁻, and the rest of systems II, respectively.

The maximal value of this energy is observed at . The compositions closest in energy to the composition with the minimum energy are given in parentheses. Thus, the composition of the fragment (M⁺·[CrCl₆]³⁻ at which its interaction with the environment is minimal () is close to the most stable compositions mentioned in Figure 3 (, ).

Finally, let us consider the function ΔE versus in Figure 4. It was obtained in the following way. First, two types of the interaction energy of the cation OS shell (M⁺ were calculated: (i) the energy of interaction with the complex and (ii) the energy of interaction with the rest of systems II . The ΔE value in Figure 4 is equal to the difference:

If we denote such fragments as [CrCl₆], (M_n), and the rest of system by characters , , and , respectively, then the first energy corresponds to equation . The second energy, that is, , characterizes the equation , while the energy ΔE corresponds to equilibrium .

As follows from Figure 4, in the range –6, the energy of OS shell interaction (M⁺· with the complex is larger than with the external rest of system in all model systems. Maximum of this function corresponds to the value .

Thus, the most stable particles in systems II are (M⁺·[CrCl₆]³⁻ with (Na, K) or 5 (Cs), as follows from Figure 1(b). According to Figure 4, such fragments can be considered as a single particle relatively weakly bound to the surroundings; they have a positive charge for all M.

3.2. Comparison of the System Energies before and after Electron Transfer

Let us introduce the following notation. is initial state of the system (before electron transfer); is final state of the system (after electron transfer); is state of the system with the geometric structure , but after the electron transfer (ET); is state of the system with the geometric structure , but before ET; is the energy of system in state ; ; ; . Electron transfer here means the transfer of one electron in the cathode process to the complex of chromium: [CrX₆]³⁻ + ē → [CrX₆]⁴⁻ (X – F, Cl).

The energy values ΔE are given in the Table 1 for systems I (M⁺·[CrF₆]³⁻ and (M⁺·[CrCl₆]³⁻ (–4) and systems II M₃CrF₆ + 18MCl, M₃CrF₆ + 18MF, and M₃CrCl₆ + 18MCl (M – Na, K, Cs).

Here we should pay attention to the following. In the range of –4, the energy values ΔE2 and ΔE3 change the sign in the system (M⁺·[CrF₆]³⁻. Therefore, in this range of the n value mechanism of ET will depend on the composition of the electroactive particle. In particular for Na-system at there is intermediate structure state in which the intersection of the potential energy surfaces of states, before and after ET, occurs. This fact is the basis for estimation of the ET activation energy according to the Marcus theory [10]. If composition of the electroactive species in the Na-system is close to the one with equal 4, then ET may occur directly to the complex in the initial state (in case when there are no other restrictions). These considerations are also applicable to K-, Cs-systems.

In the system (M⁺·[CrCl₆]³⁻ only functions ΔE2 change the sign. However, this is quite enough to use for states with and 4, the conclusions obtained above in relation to the initial state . Low stability of the state in comparison with the state is due to the strong influence of the Jahn-Teller effect on the structure of the [CrCl₆]⁴⁻ complex after ET. For this reason, the state is strongly destabilized and the ΔE3 energy is less than zero for all values of .

In the extended systems II the ΔE values are less than zero for all M. In this case, formally electron transfer can occur at the initial state . Some comments on this conclusion are made in the next section.

4. Conclusions

Results of this work indicate the possibility of ET directly to the initial particle , rather than through a transition state . This finding needs to be clarified.

Earlier [21] analysis of the nature of frontier molecular orbitals was held in systems M₃CrX₆ + 18MCl (X: F, Cl; M: Na, K, and Cs) and it was shown that ET through the OS cations is only possible in system Cs₃CrF₆ + 18CsCl. In other systems, the nature of the LUMO (lowest unoccupied molecular orbital) requires preliminary adjustment of the OS shell to provide direct contact ligands of the first coordination sphere with the electrode surface. That is, in these systems, ET has to take place in a transition state .

In work [21] it was also suggested that for the system Cs₃CrF₆ + 18CsCl charge transfer takes place in a state with geometric structure which is close to the structure of the initial system. In this work we obtain confirmation of this assumption.

In the Cs₃CrF₆ + 18CsCl system ET occurs with abnormally high rate and apparently is limited by diffusion. The approach of the active particles to the electrode surface requires the overcoming of the activation barrier, which depends on the particle charge. Data related to the composition and the possible charge of such particles can be obtained from data such as in Figures 1, 3, and 4. An example of use of such approach is contained in the work [22] where the abnormal rate of ET in the system Na₂NbCl₆ + 18MCl in comparison with systems M₂NbCl₆ + 18MCl (M – K, Cs) has been explained by distinction of charge of the electroactive species.

Furthermore, according to calculations, in some cases a higher rate of ET can occur for particles with less stability [13].

Thus, analysis of the mechanism of electrochemical ET in melts requires taking into account a number of factors. Nevertheless, the main conclusion of this work remains true: the search of the transition state in the considered systems should be carried close to the initial state . This greatly simplifies a problem and transforms it from a purely theoretical sphere to the field of practical tasks that do not require exceptional cost of computer time.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by the Russian Foundation for Basic Research (Project no. 15-03-02290-a).

References

S. A. Kuznetsov, V. G. Kremenetskii, A. V. Popova, O. V. Kremenetskaya, and A. V. T. Kalinnikov, “Unusual effect of the second coordination sphere on the standard charge transfer rate constants for the Nb(V)/Nb(IV) redox couple in chloride–fluoride melts,” Doklady Physical Chemistry, vol. 428, no. 2, pp. 770–773, 2009.
View at: Publisher Site | Google Scholar
V. G. Kremenetsky, O. V. Kremenetskaya, S. A. Kuznetsov, and V. T. Kalinnikov, “Quantum-chemical approach to the estimation of the composition of stable complex species in alkali metal halide melts,” Doklady Physical Chemistry, vol. 437, part 2, pp. 75–77, 2011.
View at: Publisher Site | Google Scholar
V. G. Kremenetsky, O. V. Kremenetskaya, and A. D. Fofanov, “Influence of the outer-sphere shell on the parameters of niobium complexes: quantum-chemical calculation,” Russian Metallurgy (Metally), vol. 2011, no. 8, pp. 748–750, 2011.
View at: Publisher Site | Google Scholar
Y. V. Stulov, V. G. Kremenetsky, and S. A. Kuznetsov, “Effect of the second coordination sphere on the standard rate constants of charge transfer for the Cr(III)/Cr(II) redox couple in chloride melts,” ECS Transactions, vol. 50, no. 11, pp. 135–152, 2012.
View at: Google Scholar
F. Bouyer and G. Picard, “Evidence of cations' influence for structural rearrangements in cryolitic melts,” in Proceedings of the 54th International Meeting of Physical Chemistry: Fast Elementary Processes in Chemical and Biological Systems, vol. 364 of AIP Conference Proceedings, pp. 532–539, 1996.
View at: Google Scholar
L. Joubert, G. Picard, B. Silvi, and F. Fuster, “Topological analysis of the electron localization function: a help for understanding the complex structure of cryolitic melts,” Journal of the Electrochemical Society, vol. 146, no. 6, pp. 2180–2183, 1999.
View at: Publisher Site | Google Scholar
V. I. Shapoval, V. V. Solov'ev, and V. V. Malyshev, “Electrochemically active species and multielectron processes in ionic melts,” Russian Chemical Reviews, vol. 70, no. 2, pp. 161–176, 2001.
View at: Publisher Site | Google Scholar
V. V. Soloviov and L. A. Chernenko, “Influence modeling of cationic composition melt on structural features of the electrochemically active particles during simultaneous and sequential electron transfer,” Research & Reviews in Electrochemistry, vol. 4, no. 5, pp. 168–173, 2013.
View at: Google Scholar
S. V. Volkov, I. A. Novoselova, V. V. Solovyev, and L. A. Chernenko, “Cation catalysis in high-temperature electrochemical synthesis of tungsten carbides in molten salts,” Research & Reviews in Electrochemistry, vol. 3, no. 3, pp. 81–89, 2012.
View at: Google Scholar
R. A. Marcus, “On the theory of oxidation-reduction reactions involving electron transfer. I,” The Journal of Chemical Physics, vol. 24, no. 5, pp. 966–978, 1956.
View at: Publisher Site | Google Scholar
S. A. Kuznetsov and V. G. Kremenetsky, “Comparison of some structural parameters of the model systems nM⁺[CrX₆]^3- and ₃M⁺[CrX₆]^3- + 18MCl (M–Na, K, Cs; X–F, Cl; n = 1–6),” ECS Transactions, vol. 64, no. 4, pp. 183–188, 2014.
View at: Publisher Site | Google Scholar
V. G. Kremenetsky and O. V. Kremenetskaya, “Quantum-chemical validation of the relative stability of complex species with different outer-sphere composition in the M₃CrF₆ + 18MCl systems (M = Na, K),” Russian Journal of Inorganic Chemistry, vol. 58, no. 12, pp. 1523–1526, 2013.
View at: Publisher Site | Google Scholar
Y. V. Stulov, V. G. Kremenetsky, and S. A. Kuznetsov, “Electrochemical and quantum-chemical studies of chromium(III,II) fluoride complexes in alkali chloride melts,” Russian Journal of Electrochemistry, vol. 50, no. 9, pp. 815–823, 2014.
View at: Publisher Site | Google Scholar
Y. V. Stulov, V. G. Kremenetsky, and S. A. Kuznetsov, “Influence of the second coordination sphere on the electrochemical behavior of chromium complexes in chloride melts: experimental and calculation methods,” International Journal of Electrochemical Science, vol. 8, no. 5, pp. 7327–7344, 2013.
View at: Google Scholar
Yu. V. Stulov, V. G. Kremenetsky, O. V. Kremenetskaya, A. D. Fofanov, and S. A. Kuznetsov, “Experimental and calculation methods of studying the effect of second coordination sphere on standard rate constants of charge transfer for Cr(III)/Cr(II) redox couple in chloride melts,” Russian Journal of Electrochemistry, vol. 47, no. 8, pp. 948–958, 2011.
View at: Publisher Site | Google Scholar
A. A. Granovsky, http://classic.chem.msu.su/gran/gamess/index.html.
M. W. Schmidt, K. K. Baldridge, J. A. Boatz et al., “The general atomic and molecular electronic structure system,” Journal of Computational Chemistry, vol. 14, pp. 1347–1363, 1993.
View at: Google Scholar
D. J. Feller, “The role of databases in support of computational chemistry calculations,” Journal of Computational Chemistry, vol. 17, no. 13, pp. 1571–1586, 1996.
View at: Publisher Site | Google Scholar
K. L. Schuchardt, B. T. Didier, T. Elsethagen et al., “Basis set exchange: a community database for computational sciences,” Journal of Chemical Information and Modeling, vol. 47, no. 3, pp. 1045–1052, 2007.
View at: Publisher Site | Google Scholar
EMSL Basis Set Library, https://bse.pnl.gov/bse/portal.
V. G. Kremenetsky and S. A. Kuznetsov, “Anomaly of the charge transfer rate in the CsCl-K₃CrF₆ melt: quantum-chemical analysis,” International Journal of Electrochemical Science, vol. 10, no. 8, pp. 6164–6174, 2015.
View at: Google Scholar
V. G. Kremenetsky and S. A. Kuznetsov, “Quantum-chemical estimation of the temperature effect on the parameters of the model electrochemical systems containing chloride complexes of niobium,” in Proceedings of the 25th Abstracts of EUCHEM 2014 Conference on Molten Salts and Ionic Liquids, p. 142, Tallinn, Estonia, July 2014.
View at: Google Scholar

Copyright

Copyright © 2016 Vyacheslav Kremenetsky and Sergey Kuznetsov. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

679

Downloads

814

Citations