Next Article in Journal
A Low-Spin CoII/Nitroxide Complex for Distance Measurements at Q-Band Frequencies
Previous Article in Journal
Magnetic Force Microscopy in Physics and Biomedical Applications
Previous Article in Special Issue
The Underexplored Field of Lanthanide Complexes with Helicene Ligands: Towards Chiral Lanthanide Single Molecule Magnets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Magnetic Field Effect on the Oxidation of Unsaturated Compounds by Molecular Oxygen

by
Evgenii M. Pliss
1,* and
Mikhail E. Soloviev
2
1
Director Institute of Fundamental and Applied Chemistry, P.G. Demidov Yaroslavl State University, Sovetskaya 14, 150003 Yaroslavl, Russia
2
Institute Digital Systems, Yaroslavl State Technical University, Moskovskiy 88, 150023 Yaroslavl, Russia
*
Author to whom correspondence should be addressed.
Magnetochemistry 2022, 8(4), 44; https://doi.org/10.3390/magnetochemistry8040044
Submission received: 17 February 2022 / Revised: 30 March 2022 / Accepted: 4 April 2022 / Published: 11 April 2022
(This article belongs to the Special Issue Reviews on Slow-Relaxation Molecules)

Abstract

:
A quantum-chemical analysis of the effect of a constant magnetic field on radical formation in the processes of chain oxidation of organic compounds by molecular oxygen is presented. The calculation of the total electronic energies and thermodynamic functions of the compounds involved in the reactions was performed by the density functional method with the hybrid exchange-correlation functional of Becke, Lee, Yang and Parr DFT B3LYP/6-311G** using the NWChem software package. The effect of the magnetic field on the individual stages of chain oxidation is associated with the evolution of radical pairs. It is assumed that the dipole–dipole interaction in a radical pair is not averaged by the diffusion of radicals and should be taken into account. To a large extent, the magnetic field effect (MFE) value is influenced by the ratio between the relaxation time of the oscillatory-excited state in the radical pair (tvib) and the relaxation time of the inter-combination transitions (tst). Although the developed technique refers to liquid-phase reactions, it can be used to study the MFE for oxidation of biologically significant compounds in multiphase systems, such as micelles, liposomes and membranes.

1. Introduction

1.1. Background and Motivation of the Research

The magnetic dependence of biosystems usually shows itself at the level of elementary biochemical processes in which spin particles originate or participate. The unpaired electrons in them happen to be spin magnetism carriers; it is these unpaired electrons that interact with magnetic fields [1,2,3,4,5,6,7,8,9]. The magnetic effect manifests itself in the competition of various channels of transformation in the elementary stages of a reaction and is caused by the dependence of the efficiency of the chemical process on the spin state of a pair of reacting particles (the theory of radical pairs) [2,5,10,11,12,13,14,15], as well as the magnetosensitivity of transitions between spin states. The effect of stationary magnetic fields on the recombination of radical pairs in molecular oxygen oxidation is one of the well-known mechanisms of interaction of the magnetic field with various biological systems. It can increase the activity, concentration and lifetime of free radicals, which causes oxidative stress, genetic mutation and/or apoptosis [1,2,4,9,10,16]. At the same time, the main role of singlet–triplet transitions in radical pairs is played by Δg, HFC and relaxation mechanisms [2,5,10,14,15,17,18,19,20,21,22,23]. A number of studies have revealed the effects of the influence of weak magnetic fields (1 Tl) on individual stages of chain oxidation of unsaturated compounds [10,17,18,19,20]. The analysis of these effects can shed some light on the characteristics of the mechanisms of oxidation of biologically significant polyunsaturated fatty acids and their esters, oxidation in micelles, membranes, etc.

1.2. The Basic Mechanism

A theoretical analysis of the effect of a stationary magnetic field (MF) on the detailed mechanisms of radical chain oxidation of organic compounds and experimental data obtained indicates that the stages in generation (initiation), propagation by the reaction of oxygen addition to alkyl radicals and chain termination are magnetically selective [2,8,9,10,17,18,19,20].
A convenient model in such studies can be created from unsaturated compounds of different classes, of which the oxidation mechanism is described by a well-known kinetic Scheme 1 [24,25].
Reaction (2) can also occur as an addition to the π-bond: RO2 + CH2=C< → ROOCH2-C < (≡ R), and the primary oxidation product is not hydroperoxide but, rather, peroxide. In accordance with the Scheme 1, the following equation is valid for the oxidation rate (W), where Wi is the chain initiation rate [24]:
W = k 2 [ RH ] W i 0.5 k 5 0.5 [ 1 + k 4 k 5 k 2 [ RH ] k 1 [ O 2 ] + k 3 k 5 ( k 2 [ RH ] k 1 [ O 2 ] ) 2 ] 0.5
At moderate temperatures (≤100 °C) and [O2] ≤ 0.3 mM at low transformation depths (≤1%), reactions (3) and (4) can be excluded and the oxidation rate (W) obeys the well-known equation [24]:
d[O2]/dt = W = k2 · k5−0.5 · [RH] · Wi0.
(if initiation chains are generated by the reaction of molecular oxygen with a substrate molecule (O2 + RH), the value of chain nucleation rate W0 is used instead of Wi [24,25]).
Obtaining kinetic information about the influence of the magnetic field on each individual stage is experimentally difficult since we deal with a multi-stage radical chain process. Therefore, to find the magnetic field effect (MFE), a methodology of the research of the detailed mechanism of chain oxidation reactions (a detailed description can be found in monographs) [24,25], and a set of experimental and theoretical methods, including quantum calculations and kinetic computer modeling, are used.
Detailed kinetic analyses of the chain propagation and termination reactions were carried out in [10,17,18,19,20]. Therefore, the present work is only devoted to the quantum-chemical analysis of the effect of a constant magnetic field on chain oxidation initiation reactions in the reactions of initiator decay and chain generation in the processes of molecular oxygen addition to the π-bond of the oxidized substrate.

2. Methods

The calculation of total electronic energies and thermodynamic functions of the compounds involved in the reactions were performed using the quantum chemical density functional method [26,27] with the hybrid exchange-correlation functional of Bechke [28], Lee, Yang and Parr [29] DFT B3LYP/6-311G** [30] using the NWChem software package [31]. The calculation of ∆g-tensors and A-tensors of radicals was carried out using the DFT method in the approximation of a two-component relativistic regular approximation of zero-order (approximate two-component relativistic zeroth-order regular approximation) (ZORA) [32,33], implemented in the NWChem package [34].
In the process of molecular dynamic modeling, the energy and gradient at each step of the integration of the classical equation of motion were calculated using the same DFT B3LYP/6-311G** quantum chemical method. The initial rate of distribution corresponded to the Maxwell–Boltzmann distribution. A thermostat was used to maintain the temperature [35], and the relaxation parameter was equal to 100 integration steps in time.

3. Results and Discussions

3.1. Reactions Chain Initiation of Chain Oxidation

The AIBN was used as an initiator decay by a known mechanism ˃(CN)CN=NC(CN)˂ → ˃ 2(CN)C + N2 through the formation of geminal pairs [r● ●r] [2,10,12,13]. During the thermal initiator decay, a radical pair arises under the action of an oscillatory quantum of energy, that is, before the decay, the molecule is in an excited state of oscillation. At the same time, its initial spin state is a singlet.
In the future, depending on the relaxation time of the vibrationally excited state, the pair can recombine by emitting an oscillatory quantum or the radicals can escape into the volume. In this case, the release of radicals into the volume can occur in the initial singlet state or a singlet–triplet conversion can occur in a radical pair [10,11,12,13,14]. Henceforth, depending on the relaxation time of the oscillational excited state, the pair can recombine by emitting an oscillatory quantum, or the radicals can exit into the volume. In this case, the release of radicals into the volume can occur either in the initial singlet state or a singlet–triplet conversion can occur in a radical pair [10,11,12,13,14].
When analyzing the MFE on the initiation reaction, it is necessary to take into account the general laws of the spin dynamics of radical pairs which are also considered in detail in a number of review papers [10,11,12,13,14].
Among all the effects that determine the MFE on the probability of a singlet–triplet transition in a radical pair, we must first distinguish the interaction of unpaired electrons with an external magnetic field, characterized by the g-tensor (the Zeeman effect) and the interaction of electrons with magnetic nuclei of radicals that are characterized by the hyperfine interaction tensor A (the Fermi effect). In qualitative analysis of the MFE on the oxidation kinetics, it is sufficient to compare the orders of these effects for interacting radicals [17,18,19,20]. However, for a detailed analysis, a computer modeling of the effects under consideration is necessary since the mismatch of the oscillatory dynamics of radicals in a pair in the pre-emission state can change both their g-tensors and hyperfine interaction tensors, which can affect the probability of singlet–triplet conversion.
In ordinary organic molecules, the excited triplet level has a significantly higher energy than the main singlet and spin-forbidden singlet–triplet transitions are unlikely. At the same time, for initiators with lower bond decay energy, the number of oscillatory levels increases rapidly with increasing temperature and the probability of singlet–triplet transitions should also increase simultaneously. Such effects are observed, for example, during thermo- and photoinitiated radical polymerization [36,37]. In this regard, quantum chemical calculations and molecular dynamic modeling of the processes of radical decay of AIBN in singlet and triplet states were carried out in order to assess the possible influence of oscillatory excited states on the probability of singlet–triplet transitions. The integration of the equations of motion was carried out in accordance with the Verlet algorithm; the integration step in time was 0.2419 fs. In the course of quantum chemical calculations, cross sections of the potential energy surface were studied along the oscillatory coordinate, which is the distance R between carbon atoms C1…C5 (Figure 1), thus characterizing the distance between the fragments of the radical pair.
For this purpose, conditional optimization of the geometry of the system was performed at fixed values of the coordinates of atoms C1 and C5 with different values of R. Figure 2 shows the dependence of the difference in the energy of the ΔE system in the singlet spin state with respect to the global minimum on the value of R.
Molecular structures are shown to be near their local extremes (local maximum and local minimum). It can be seen that dependence can be approximated by a smooth curve. After reaching the maximum, the value of which is about 135 kJ/mol, a local minimum corresponding to the formation of a singlet radical pair is observed. The radical pair structure, with a symmetrical placement of radicals relative to the nitrogen molecule located in the center, corresponds to the concert decay mechanism characteristic of azo initiators [38]. The nitrogen molecule located between the radicals will prevent their recombination in the singlet state and promote the release of radicals into the volume. Figure 3 shows a similar dependence for a system in a triplet state.
The minimum energy of the triplet state lies significantly higher than of the main singlet state, which corresponds to the prohibition of singlet–triplet transitions in the main oscillatory state. However, with the increase in R, the situation changes. When the R value reaches about 0.42 nm, there is a break in the cross section of the potential surface with a sharp drop in the energy value. At the same time, as can be seen from the figure, there is a significant restructuring of the molecular structure with the release of a nitrogen molecule into the volume and the formation of a triplet radical pair. According to the spin prohibition, such a pair should not recombine, and the radicals should also be released into the volume. This conclusion is confirmed by the results of molecular dynamic modeling in the quantum-classical approximation.
Figure 4 shows the time dependences of R in molecular dynamic modeling for AIBN molecules in singlet and triplet states.
For a molecule in a singlet state the R distance does not undergo significant changes, and instead fluctuates relative to the average value corresponding to the minimum potential energy. For a molecule in the triplet state, fluctuations are also initially observed relative to approximately the same value of R; however, after a time of about 250 fs, the distance between the radicals begins to systematically increase, which indicates the decay of the radical pair and the release of radicals into the volume. At this moment a rearrangement of the structure takes place accordant to a break at the intersection of the potential surface (Figure 3). This fact is proved by energy distribution types for the molecule in different spin states obtained in molecular dynamics simulation (Figure 5).
At this moment a rearrangement of the structure takes place accordant to a break at the intersection of the potential surface (Figure 3). This fact is proved by energy distribution types for the molecule in different spin states obtained in molecular dynamics simulation (Figure 5).
Curve 1 corresponds to the energy distribution of the molecule in the singlet state. Curve 2 corresponds to the distribution of the molecule in the triplet state at simulation time up to 250 fs. Curve 3 corresponds to the distribution of the molecule in the triplet state at simulation time after 250 fs. There are characteristic symmetric distributions for all three curves; however, for a molecule in a triplet state, the dispersion of distributions is slightly higher. It can be seen that curves 2 and 3 correspond to different states. A sharp decrease in energy after 250 fs for the triplet molecule corresponds to a gap in the section of the potential surface in Figure 3 with the exit of the nitrogen molecule from the radical pair. At the same time, as expected, the triplet radical pair is not stable, and the radicals leave into volume. The time required for the relaxation transition from state 2 to state 3 can be characterized as the relaxation time of the oscillatory-excited state—tvib.
The question of whether singlet–triplet transitions occur during the decay of the initiator, and accordingly, whether an external magnetic field will affect the probability of decay, is not obvious. Triplet–singlet transitions are observed during the photoinduced decay of initiators [36,37], while an external magnetic field can have a significant effect. That is, the lifetime of triplet radical pairs is sufficient for an inter-combination transition, but a singlet radical pair must recombine quickly. An analysis of the decomposition products of azo initiators suggests that such transitions do take place. Particularly in the case of the thermally initiated decomposition, all geminal radical pairs in the absence of inter-combination transitions should be singlets; however, the yield of the recombination products of geminal radical pairs seems to be underestimated, given the significant cellular effect [24,38].
On the one hand, the identical structure of the radicals in the AIBN decomposition products, which results in the equality of their g-factors in the main oscillatory state, should lead to the absence of the influence of the magnetic field on the probability of conversion by the Δg mechanism. On the other hand, when the process occurs via the cleavage of two bonds with the formation of two radicals and a nitrogen molecule, as in the singlet radical pair in Figure 2, the correlation between the formed radicals will be weakened so that they can be at different oscillatory levels with different values of the oscillatory excited state time (tst). In this case, their g- and A-tensors will differ, which will contribute to the inter-combination transition. In order to evaluate this effect, g- and A-tensors of the (CH3)2C∙(CN) radicals of singlet and triplet AIBN molecules were calculated at the minimum energy and after 1000 fs of molecular dynamic modeling, which is when the correlation between them is practically absent, and their configurations are very likely to correspond to different oscillating levels. The Table 1 shows the eigenvalues of g-tensors calculated at these times, and as an invariant of A-tensors, the Aeff values calculated by the formula [29] are given
A e f f = ( A A T ) 1 2
where the A symbol denotes the A-tensor in diagonal form, and summation is performed over all the protons of the radical.
From the above data, the eigenvalues of the radical g-tensors change during the simulation for both singlet and triplet radical pairs. Although these changes are not very significant, but taking into account the fact that, as can be seen from Figure 5, the energy distributions of the singlet (curve 1) and triplet states (curve 3) do intersect, at certain times the singlet-triplet transition may turn out to be favorable, so such transitions can occur in a strong magnetic field.
The Aeff values indicate that, in strong magnetic fields (MF), S-T0 inter-combination transitions are also possible via the mechanism of hyperfine interaction with characteristic relaxation times tst.
Thus, the calculation performed indicates the possible effect of an external constant (MF) on the kinetics of the release of radicals into the volume during the AIBN thermal decomposition, which in turn should lead to an increase in the rate of oxidation initiation; this corresponds to the observed experimental data [19]. The degree of influence of an external MF on the rate of initiation will obviously depend in each specific case on the ratio of the relaxation times tvib and tst.

3.2. Chain Generation by Reaction with Dioxygen

Let us consider the reaction of oxygen addition to a double bond in the example of styrene. Here, as in the previous case, the MFE mechanism is realized through the influence of a magnetic field on the probability of a triplet–singlet conversion of a radical pair during the reaction.
In [17,19], based on the experimental data obtained, it was shown that the effect of MF is determined by the mechanism of interaction of O2 with the π-bond. The reaction proceeds through the formation of a triplet biradical, which initiates oxidation chains: O2 + CH2=C< → O2 + C=C → [OOC−C]TOOCH2-C< via reactions with O2 and π-bonds. The Zeeman interaction stimulates triplet–singlet conversion through the T–S transition leading to the formation of non-radical compounds of molecular products (dioxetane and its decay products [24,38]): Magnetochemistry 08 00044 i001
As a result, the share of radical initiation decreases, and the sign of MFE will be negative [19]. At the same time, not all details of this mechanism have been described due to the lack of appropriate quantum chemical calculations.
In the previous case, the magnetic field leads to an increase in S–T transitions; then, the opposite picture takes place in the nucleation reaction: the magnetic field changes the probability of T–S transitions. Since the stable spin state of molecular oxygen is a triplet, the pre-reaction complex, including styrene and oxygen molecules bound by intermolecular interaction, is also a triplet. As a result of the reaction, primary radical pairs may arise, which are also in the triplet state. When unpaired electrons of radical pairs are attached to olefin double bonds, a chain oxidation process begins. At the same time, as a result of triplet–singlet transitions, molecular compounds can be formed in radical pairs. Since in this case, unlike with the azo initiator initiations, singlet states are lower in energy than the initial triplet, such processes will remove the resulting compounds from the chain process and the MFE sign will be negative. Figure 6 shows the calculated variants of the structures of compounds in the minimum potential energy, reflecting the various stages of formation and evolution of radical pairs that occur when oxygen is attached to the double bond of a styrene molecule.
For each structure, the values of the change in its energy (ΔE) with respect to the sum of the total electronic energies of isolated styrene and oxygen molecules are indicated. The T symbol at the Roman numeral corresponding to the structure number indicates a triplet state; the S symbol is a singlet state. The IT structure is an initial pre-reaction complex. Structure III is an internal radical pair that occurs when a diatomic oxygen molecule is attached to one of the carbon bonds of a styrene molecule. The IVT structure is a radical pair that occurs when an oxygen molecule is attached to a styrene double bond by two atoms simultaneously. In [19], it was assumed that the IVT structure emerges during the transformation of the IIT structure. However, as can be seen from the energy changes values, this process is unlikely and the process leading to the IIIT structure is more likely. Both radical pairs IVT and IIIT can initiate two kinetic oxidation chains; however, the difference between them is that in the case of the IVT structure, these kinetic chains will start on the same molecule, whereas in the case of the IIIT structure, they will start on two different molecules. The second molecule initiating the chain, in this case, will be a hydroxyl radical. When a constant magnetic field is applied, T–S transitions will be initiated, which in this reaction can be both T+ → S and T → S depending on the magnetic field strength, and with increasing field strength, the proportion of T → S transitions should decrease. Since these transitions will lead to the formation of molecular compounds, the share of radical pairs involved in the initiation of oxidation chains will decrease and the MFE sign will be negative.
The reagents evolution scheme in the studied initiation process is shown in Figure 7.
Since, for all possible types of radical pairs, the final singlet states of the reaction products are lower in energy than triplet ones, the course of these reactions will lead to a decrease in the share of the chain oxidation mechanism as well as to the rate of the entire oxidation process (Figure 7).

4. Conclusions

It can be concluded that the influence of the magnetic field on individual stages of chain oxidation is associated with the evolution of radical pairs. The detected effects may be due to the kinetic features of the detailed mechanism of the chain process. The answer to this question should be sought in a combination of theoretical analysis with experimental study of elementary reactions. The calculations presented in this paper suggest a strong dependence of the MFE on the dipole–dipole interaction in the radical pair. This interaction is not averaged by the diffusion of radicals and must be taken into account. To a large extent, the MFE value is influenced by the ratio of the relaxation time of the oscillatory-excited state in the radical pair, tvib, and the relaxation time of the inter-combination transitions, tst. Although the developed technique refers to liquid-phase reactions, it can nevertheless be used to study the MFE for the oxidation of biologically significant compounds in multiphase systems, such as micelles, liposomes and membranes.

Author Contributions

All authors (E.M.P. and M.E.S.) made an equal contribution to the article. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian National Scientific Fund, grant number 14-23-00018.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No any data supporting results have been published anywhere before.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hayashi, H. Introduction to Dynamic Spin Chemistry: Magnetic Field Effects on Chemical and Biochemical Reactions; World Scientific Printers (S) Pte Ltd.: Singapore, 2004; 236p. [Google Scholar] [CrossRef]
  2. Buchachenko, L. Magneto-Biology and Medicine; Nova Science Publishers: New York, NY, USA, 2014; Volume 144, p. 3. [Google Scholar]
  3. Kokhan, V.; Matveeva, M.; Mukhametov, A.; Shtemberg, A. Risk of defeats in the central nervous system during deep space missions. Neurosci. Biobehav. Rev. 2016, 71, 621–632. [Google Scholar] [CrossRef] [PubMed]
  4. Binhi, V.; Prato, F. Biological effects of the hypomagnetic field: An analytical review of experiments and theories. PLoS ONE 2017, 12, e0179340. [Google Scholar] [CrossRef] [PubMed]
  5. Salikhov, K.; Sagdeev, R.M.Y.; Buchachenko, A.A. Spin Polarisation and Magnetic Field Effects in Radical Reaction; Salikhov, K.M., Molin, Y.N., Sagdeev, R.A., Buchachenko, R.L.A., Eds.; Elsevier: Amsterdam, The Netherlands; Akademia Klado: Budapest, Hungary, 1984; 419p. [Google Scholar]
  6. Kazin, V.; Guzov, E.; Pliss, E.; Moshareva, V.; Makaryin, V.; Levshin, N.; Baranov, A. The effect of a constant magnetic field on components of protein structures in human blood. Biophysics 2017, 62, 998–1007. [Google Scholar] [CrossRef]
  7. Zhang, B.; Tian, L. Reactive oxygen species: Potential regulatory molecules in response to hypomagnetic field exposure. Bioelectromagnetics 2020, 41, 573–580. [Google Scholar] [CrossRef]
  8. Buchachenko, A.; Pliss, E. Isotope Effects of Hydrogen and Atom Tunnelling//Russian Chemical Reviews. 2016. V. 85. I. 6. p. 557–564. Available online: https://www.elibrary.ru/item.asp?id=27066411 (accessed on 10 February 2022).
  9. Sampson, C.; Keens, R.; Kattnig, D. On the magnetosensitivity of lipid peroxidation: Two- versus three-radical dynamics. Phys. Chem. Chem. Phys. 2019, 21, 13526. [Google Scholar] [CrossRef] [Green Version]
  10. Lukzen, N.; Ivanov, K.; Sadovsky, V.; Sagdeev, R. Magnetic field effect on recombination of radicals diffusing on a two-dimensional plane. J. Chem. Phys. 2020, 152, 34103. [Google Scholar] [CrossRef]
  11. Eichwald, C.J.; Walleczek, J. Low-frequency-dependent effects of oscillating magnetic fields on radical pair recombination in enzyme kinetics. J. Chem. Phys. 1997, 107, 4943. [Google Scholar] [CrossRef]
  12. Woodward, J. Radical pairs in solution. Prog. React. Kinet. Mech. 2002, 27, 165–207. [Google Scholar] [CrossRef]
  13. Buchachenko, A. Experimental testing of molecular dynamic function of pairs by the isotope selectivity of radical recombination. Russ. Chem. Bull. 1995, 44, e1571–e1577. [Google Scholar] [CrossRef]
  14. Adam, B.; Sinayskiy, I.; Petruccione, F. An open quantum system approach to the radical pair mechanism. Sci. Rep. 2018, 8, 15719. [Google Scholar] [CrossRef]
  15. Buchachenko, A.; Bukhvostov, A.; Ermakov, K.; Kuznetsov, D. A specific role of magnetic isotopes in biological and ecological systems. Physics and biophysics beyond. Prog. Biophys. Mol. Biol. 2020, 155, 1–19. [Google Scholar] [CrossRef] [PubMed]
  16. Usselman, R.; Hill, I.; Singel, D.; Martino, C. Spin biochemistry modulates Reactive Oxygen Species (ROS) production by radio frequency magnetic fields. PLoS ONE 2014, 9, 93065. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Pliss, E.; Grobov, A.; Kuzaev, A.; Buchachenko, A. Magnetic field effect on the oxidation of hydrocarbons by molecular oxygen. Mendeleev Commun. 2017, 27, 246–247. [Google Scholar] [CrossRef]
  18. Pliss, E.; Grobov, A.; Kuzaev, A.; Buchachenko, A.L. Magnetic field effect on the oxidation of organic substances by molecular oxygen. J. Phys. Org. Chem. 2019, 32, e3915. [Google Scholar] [CrossRef]
  19. Pliss, E.; Grobov, A.; Kuzaev, A.; Buchachenko, A. Magnetic field as a means to identify initiating reaction in oxidation of organic substances by molecular oxygen. Mendeleev Comm. 2020, 30, 433–435. [Google Scholar] [CrossRef]
  20. Pliss, E.; Grobov, A.; Kuzaev, A.; Buchachenko, A. Magnetic field effects on the initiation of chain oxidation. Mendeleev Commun. 2021, 31, 341–342. [Google Scholar] [CrossRef]
  21. Fedin, M.; Purtov, P.; Bagryanskaya, E. Spin relaxation of radicals in low and zero magnetic field. J. Chem. Phys. 2003, 118, 192–201. [Google Scholar] [CrossRef]
  22. Gohdo, M.; Wakasa, M.; Kitahama, Y.; Sakaguchi, Y.J. Magnetic field effects due to the relaxation mechanism observed for the photo-induced electron transfer reaction of zinc (II) tetraphenylporphyrin and 2-methyl-1, 4-naphthoquinone. Photochem. Photobiol. A Chem. 2008, 199, 130–135. [Google Scholar] [CrossRef]
  23. McKenzie, I. Spin relaxation of a short-lived radical in zero magnetic field. Phys. Chem. Chem. Phys. 2011, 13, 1168–1173. [Google Scholar] [CrossRef]
  24. Denisov, E.; Afanas’ev, I. Oxidation and Antioxidants in Organic Chemistry and Biology; CRC Press: Boca Raton, FL, USA, 2005; p. 992. [Google Scholar] [CrossRef]
  25. Pliss, E.; Safiulin, R.; Zlotsky, S. Inhibited Oxidation of Unsaturated Compounds. In Kinetics, Mechanism, Correlation of Structure with Reactionary Ability; LAP LAMBERT Academic Publishing: Saarbruchen, Germany, 2012; p. 130. [Google Scholar]
  26. Hohenberg, P.; Kohn, W. Inhomogeneous electron gas. Phys. Rev. 1964, 136, B864–B871. [Google Scholar] [CrossRef] [Green Version]
  27. Kohn, W.; Sham, L. Self-consistent equations including exchange and correlation effects. Phys. Rev. 1965, 140, A1133. [Google Scholar] [CrossRef] [Green Version]
  28. Becke, D.A. The performance of a family of density functional methods. J. Chem. Phys. 1993, 98, 5648. [Google Scholar] [CrossRef] [Green Version]
  29. Miehlich, B.; Savin, S.; Preis, N. Results obtained with the correlation energy density functionals of Becke and Lee, Yang and Parr. Chem. Phys. Lett. 1989, 157, 1477–1489. [Google Scholar] [CrossRef]
  30. Mueller, M. Fundamentals of quantum chemistry. Molecular spectroscopy and modern electronic structure computing. Science 2001, 2, 291–292. [Google Scholar] [CrossRef] [Green Version]
  31. Apra, E.; Bylaska, E.J.; De Jong, W.A.; Govind, N.; Kowalski, K.; Straatsma, T.P.; Valiev, M.; van Dam, H.J.; Alexeev, Y.; Anchell, J.; et al. NWChem: Past, present, and future. J. Chem. Phys. 2020, 152, 184102. [Google Scholar] [CrossRef] [PubMed]
  32. Van Lenthe, E.; Baerends, E.; Snijders, J. Relativistic total energy using regular approximations. J. Chem. Phys. 1993, 101, 9783–9792. [Google Scholar] [CrossRef]
  33. Van Lenthe, E.; Ehlers, A.; Baerends, E. Geometry optimization in the Zero Order Regular Approximation for relativistic effects. J. Chem. Phys. 1999, 110, 8943–8953. [Google Scholar] [CrossRef] [Green Version]
  34. Nichols, P.; Govind, N.; Bylaska, E.; de Jong, W. Gaussian basis set and planewave relativistic spin−orbit methods in NWChem. J. Chem. Theory Comput. 2009, 5, 491–549. [Google Scholar] [CrossRef]
  35. Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007, 126, 014101. [Google Scholar] [CrossRef] [Green Version]
  36. Khudyakov, I.; Arsu, N.; Jockusch, S.; Turro, J. Magnetic and spin effects in the photoinitiation of polymerization. Desig. Monomers Polym. 2003, 6, 91–101. [Google Scholar] [CrossRef] [Green Version]
  37. Rintoul, I. Kinetic control of aqueous polymerization using radicals generated in different spin states. Processes. Open Access Chem. Biolog. Proc. Eng. J. 2017, 5, 15. [Google Scholar] [CrossRef] [Green Version]
  38. Denisov, E.; Denisova, T.; Pokidova, T. Handbook of Radicals Initiators; Wiley: New York, NY, USA, 2003; p. 886. [Google Scholar] [CrossRef] [Green Version]
Scheme 1. The mechanism of radical chain oxidation of unsaturated compounds, where I—initiator, RH—unsaturated compound or its saturated analog.
Scheme 1. The mechanism of radical chain oxidation of unsaturated compounds, where I—initiator, RH—unsaturated compound or its saturated analog.
Magnetochemistry 08 00044 sch001
Figure 1. The AIBN structure at the minimum energy.
Figure 1. The AIBN structure at the minimum energy.
Magnetochemistry 08 00044 g001
Figure 2. The dependence of the energy changes relative to the global minimum of the AIBN molecule in the singlet state on the distance between carbon atoms C1 and C5 (Figure 1).
Figure 2. The dependence of the energy changes relative to the global minimum of the AIBN molecule in the singlet state on the distance between carbon atoms C1 and C5 (Figure 1).
Magnetochemistry 08 00044 g002
Figure 3. The dependence of the energy change relative to the global minimum of the AIBN molecule in the triplet state on the distance between carbon atoms C1 and C5 (Figure 1).
Figure 3. The dependence of the energy change relative to the global minimum of the AIBN molecule in the triplet state on the distance between carbon atoms C1 and C5 (Figure 1).
Magnetochemistry 08 00044 g003
Figure 4. The dependence of the distance between C1 and C5 atoms on the time of molecular dynamic modeling: 1—a molecule in a singlet state, 2—in a triplet state.
Figure 4. The dependence of the distance between C1 and C5 atoms on the time of molecular dynamic modeling: 1—a molecule in a singlet state, 2—in a triplet state.
Magnetochemistry 08 00044 g004
Figure 5. The energy distributions of the AIBN molecule in molecular dynamic modeling: 1—in the singlet state, 2—in the triplet state before the time of 250 fs, 3—in the triplet state after the simulation time of 250 fs.
Figure 5. The energy distributions of the AIBN molecule in molecular dynamic modeling: 1—in the singlet state, 2—in the triplet state before the time of 250 fs, 3—in the triplet state after the simulation time of 250 fs.
Magnetochemistry 08 00044 g005
Figure 6. Structures emerging during the reaction of formation and evolution of radical pairs during auto initiation of styrene oxidation: a—IT, ∆E = −20.8 kJ/mol; b—IIT, ∆E = 74.6 kJ/mol; c—IIT, ∆E = −33.0 kJ/mol; d—IVT, ∆E = 289.5 kJ/mol; e—IS, ∆E = −80.9 kJ/mol; f—IIS, ∆E = −46.3 kJ/mol.
Figure 6. Structures emerging during the reaction of formation and evolution of radical pairs during auto initiation of styrene oxidation: a—IT, ∆E = −20.8 kJ/mol; b—IIT, ∆E = 74.6 kJ/mol; c—IIT, ∆E = −33.0 kJ/mol; d—IVT, ∆E = 289.5 kJ/mol; e—IS, ∆E = −80.9 kJ/mol; f—IIS, ∆E = −46.3 kJ/mol.
Magnetochemistry 08 00044 g006
Figure 7. Scheme of the evolution of radical pairs during auto initiation of styrene oxidation.
Figure 7. Scheme of the evolution of radical pairs during auto initiation of styrene oxidation.
Magnetochemistry 08 00044 g007
Table 1. C∙(CN) radicals of singlet and triplet AIBN molecules at minimum energy and after 1000 fs of molecular dynamic modeling.
Table 1. C∙(CN) radicals of singlet and triplet AIBN molecules at minimum energy and after 1000 fs of molecular dynamic modeling.
Parametert = 0t = 1 ps
SingletTriplet
gxx2.0022262.0023882.002331
gyy2.0031322.0034122.003585
gzz2.0038512.0040302.003647
Aeff ·10−8, Hz2.391.842.80
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pliss, E.M.; Soloviev, M.E. Magnetic Field Effect on the Oxidation of Unsaturated Compounds by Molecular Oxygen. Magnetochemistry 2022, 8, 44. https://doi.org/10.3390/magnetochemistry8040044

AMA Style

Pliss EM, Soloviev ME. Magnetic Field Effect on the Oxidation of Unsaturated Compounds by Molecular Oxygen. Magnetochemistry. 2022; 8(4):44. https://doi.org/10.3390/magnetochemistry8040044

Chicago/Turabian Style

Pliss, Evgenii M., and Mikhail E. Soloviev. 2022. "Magnetic Field Effect on the Oxidation of Unsaturated Compounds by Molecular Oxygen" Magnetochemistry 8, no. 4: 44. https://doi.org/10.3390/magnetochemistry8040044

APA Style

Pliss, E. M., & Soloviev, M. E. (2022). Magnetic Field Effect on the Oxidation of Unsaturated Compounds by Molecular Oxygen. Magnetochemistry, 8(4), 44. https://doi.org/10.3390/magnetochemistry8040044

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop