Theoretical Investigation on H₂O₂-Ng (He, Ne, Ar, Kr, Xe, and Rn) Complexes Suitable for Stereodynamics: Interactions and Thermal Chiral Rate Consequences

Abstract

Although molecular collisions of noble gases (Ng) can be theoretically used to distinguish between the enantiomers of hydrogen peroxide - H₂O₂ (HP), little is known about the effects of HP-Ng interactions on the chiral rate. In this work, the chiral rate as a function of temperature (CRT) between enantiomeric conformations of HP and Ng (Ng=He, Ne, Ar, Kr, Xe, and Rn) are presented at MP2(full)/aug-cc-pVTZ level of theory through a fully basis set superposition error (BSSE) corrected potential energy surface. The results show that: (a) the CRT is highly affected even at a small decrease in the height of trans-barrier; (b) its smallest values occur with Ne for all temperatures between 100 and 4,000 K; (c) that the decrease of CRT shows an inverse correlation with respect to the average valence electron energy of the Ng and (d) Ne and He may be the noble gases more suitable for study the oriented collision dynamics of HP. In addition to binding energies, the electron density ρ and its Laplacian ∇²ρ topological analyses were also performed within the atoms in molecules (AIM) theory in order to determine the nature of the HP-Ng interactions. The results of this work provide a more complete foundation on experiments to study HP's chirality using Ng in crossed molecular beams without a light source.

1. Introduction

Hydrogen peroxide - H₂O₂ (HP) is a molecule of interest in a large and diverse number of fields in addition to its industrial uses. For example, it has emerged as a major metabolite in redox signaling and regulation (Antunes and Brito, 2017; Sies, 2017), and its presence was observed in Martian atmosphere (Encrenaz et al., 2004) and also on the surface of Jupiter's moon Europa (Carlson et al., 1999). The HP is interesting since it is simplest molecule that exhibits internal (torsional) rotation and chirality. Furthermore, this molecule can form dimers (Dobado and Molina, 1993; González et al., 1997), clusters (Yu and Yang, 2011), complexes with water (Mo et al., 1994; González et al., 1997) and with biologically important molecules such as adenine (Dobado and Molina, 1999), DNA (Piatnytskyi et al., 2016), glycine (Shi and Zhou, 2004) or nitrosamines (Roohi et al., 2010). These features indicate that HP should be a better proton donor for hydrogen bonding than water. Thus, the understanding of how the relative orientation of the O-H can lead to a weakly complex or a chemical reaction has also been paid considerable attention due to its implication in atmospheric chemistry and oxidation reactions (Lundell et al., 1998, 2001; Daza et al., 2000; Goebel et al., 2000, 2001a,b, 2002; Molina et al., 2002; Pehkonen et al., 2004; Mucha and Mielke, 2009; Grzechnik et al., 2013). Moreover, HP's properties have been investigated, such as its isolated chirality (Roncaratti and Aquilanti, 2010), stereomutation (Fehrensen et al., 2007; Bitencourt et al., 2008), size-dimensional wave packets (Wang et al., 2012), spectroscopy (Hunt et al., 1965; Małyszek and Koput, 2013; Al-Refaie et al., 2015) and rotation barriers (Song et al., 2005).

On the other hand, the hydrogen peroxide seems to be a prototypical model to be used into experiments to observe chirality in crossed molecular beam without a light source (Palazzetti et al., 2013), a frontier in research of stereodynamics which is still at early stages (Su et al., 2013; Lombardi and Palazzetti, 2018). In these kind of experiments, the molecular orientation control on the intense continuous beam is mandatory to the phenomena of chiral selectivity to be demonstrated (Aquilanti et al., 2005). For this reason, the interaction between HP and atoms, molecules and ions is so relevant to sterodynamics studies (Barreto et al., 2007, 2010; Lombardi et al., 2011; Roncaratti et al., 2014; Leal et al., 2016).

In the present paper we investigated the dynamics of the chiral molecule HP interacting by van der Waals forces with noble gases Ng (Ng=He, Ne, Ar, Kr, Xe and Rn) in order to obtain the chiral rate as a function of temperature (CRT) for these complexes, as well as its consequences for the chiral kinetic interconversion when tunneling effect is included. The aim of this work was to understand how the kinetic interconversion of the two HP's isomers is affected along collisional events through a basis set superposition error (Boys and Bernardi, 1970) (BSSE) corrected potential energy surface (PES). In addition, the electron density ρ(r) and its Laplacian ∇²ρ(r) topological analyses were also performed within the atoms in molecules (AIM) theory in order to determine the nature of the intermolecular interactions.

2. Methodology

2.1. Computational Details

All calculations were performed using the Gaussian09 package (Frisch et al., 2009). The structures were optimized without constraints at MP2(full) method in conjunction with aug-cc-pVTZ-PP for Xe and Rn (Peterson et al., 2003) and aug-cc-pVTZ for the remaining atoms (Dunning, 1989; Woon and Dunning Jr, 1993; Wilson et al., 1999). Vibrational frequencies at the same level of theory were also performed in order to ensure that each minimum has only positive frequencies and that each transition state has only a single imaginary frequency, as well as to obtain the zero point vibrational energy (ZPE). The counterpoise method of Boys and Bernardi (1970) was used to correct the BSSE for binding energy.

AIM analysis (Matta and Boyd, 2007) and graphic representations were performed with the AIMALL program (Keith, 2017) using the MP2(full) density (wavefunction) as input as described in the AIM theory (Dobado et al., 1998; Cortés-Guzmán and Bader, 2005).

2.2. Overview of the Transition State Theory

The transition state theory (TST)(Truhlar et al., 1996) was developed primarily by Henry Eyring (Eyring, 1935) and Michael Polanyi (Polanyi and Wigner, 1928) between 1928 and 1935. The TST is an improvement over the so-called theory of collisions (Lewis, 1967), and it is widely used to calculate the rate constants of chemical reactions.

The start point of TST is the existence of a transition state (TS) between the reagents and products. Located at the top of the potential energy barrier and it assumes a quasi-equilibrium between reactants and activated transition state complexes. For a bimolecular reaction given by

The TS is characterized by a single imaginary frequency along the reaction coordinate of the molecular system which is represented here by . In its turn, the reaction coordinate can be represented by angular changes in bond distances during the chemical reaction (Henkelman et al., 2002).

The equation that determines the reaction rate is known as the Eyring equation, given by where 0 < κ(T) ≤ 1 is the so-called transmission coefficient, k_B is the Boltzmann constant, h is the Planck constant, is the standard molar partition function, N_A is the Avogadro constant, R is the gas constant and is the barrier energy with zero-point energy correction. In addition, the TS, R₁ and R₂ subscripts stand for the transition state and reagents, respectively. Thus, the rate constant is determined by the parameters that characterize both reagents and the TS.

The general partition function is formed by the product of translational q^trans., rotational q^rot., vibrational q^vib. and electronic q^ele. partition functions. The translational partition function for a free particle with mass m moving along the length dimension l_x can be evaluated by considering that the separation of energy levels is small and that a large number of states are accessible at room temperatures. Therefore, the energy levels should be continuous and the sum contribution of the translational partition function becomes an integral. Which the solution for the three-dimensional case is (Atkins et al., 2013)

Although the system can be excited at normal modes, the energy levels are discrete for the rotational mode. The three degrees of freedom of spatial rotation and the three moments of inertia I_A, I_B and I_C must be taken into account for a non-linear molecule (Atkins et al., 2013), thus where σ is the so-called number of symmetry. The vibrational mode has reasonably spaced energy which must be taken into account since they are partially occupied. As a consequence, the vibrational partition function is strictly calculated as a sum over the occupied states. In the case of n vibrational degrees of freedom, the vibrational partition function is given by the product of n partition functions, where ν_i is each of the fundamental vibrational frequencies. In most cases, only the lowest energy state is occupied and the electronic energies should not contribute considerably to the total partition function (Atkins et al., 2013). A good approximation is to disregard the contributions of the nuclear and electronic spins and to vanish the fundamental energy level for the electronic partition function. Under these considerations the electronic partition function should be equal to unity (Atkins et al., 2013)

On the other hand, the coefficient κ(T) represents the tunneling effect of the reaction coordinate of the chemical system and it is usually important for light atoms or molecules at low temperatures. Thus, tunneling estimates were made using both Wigner (Polanyi and Wigner, 1928) and Eckart (Eckart, 1930) methods.

The Wigner tunneling correction proposes a parabolic potential, where E_b corresponds the energy potential barrier of MEP, is the imaginary frequency of transition state and s is the coordinate reaction. This implies in a transmission coefficient given by Bell (1959)

For very low temperatures, the Wigner tunneling effect is not very effective, and for this reason, it was also employed Eckart tunneling correction (Truhlar et al., 1985).

The Eckart tunneling correction uses a potential of the type where α is a parameter described by and μ is the reduced mass of the system. These parameters determine the barrier width. Here it is important to note that the A and B can be positive, negative or zero. The A parameter corresponds to the energy difference V(x → −∞) and V(x → +∞), and B is a parameter that measures the height of the barrier given by

So the most usual form for the Eckart's potential in the study of reaction rates is (Truhlar et al., 1985) where s is the coordinate of the reaction and s₀ is the reaction coordinate corresponding to the maximum of the barrier, which is given by Finally, the transmission probability (Bell, 1980), obtained through the solution of the Schrödinger equation with Eckart's potential, is expressed by the following equation where k, β and δ depend on , A, B and energy (E).

The quantum tunneling correction κ(T) can thus be calculated from the ratio between the quantum rate k_quan.(T) and the classical rate k_class.(T) in which the particles cross the barrier. Thus, the Eckart tunneling correction with transmission coefficient is given by where integration is performed over all possible energies.

3. Results and Discussion

3.0.1. Geometric Parameters, Interactions and AIM Analysis

The details about the generation of the potential energy surface are described in another work of our group (Roncaratti et al., 2014), so it will be commented briefly here. First, all HP geometry parameters were kept frozen at their equilibrium values of D_OO = 1.45Å, D_OH = 0.966Å and the angle HOO = 100.8°. The Ng's position is expressed in terms of the polar coordinates as represented in Figure 1, where R is the distance of Ng relative to the middle of O-O bond and α is the polar angle with respect to an axis perpendicular to the O-O bond (z axis). The two planes defined by O-O-H atoms are then rotated around the O-O axis, with steps of 1°. In addition, α was equal to 0°, 45°, 90° and R distance was varied from 2 to 5Å with steps 0.1Å.

Figure 1

Topological studies performed on this adjusted potential energy surface (PES) showed that the HP and HP-Ng complexes have two overall minimum configurations, termed cis (labeled as θ₋) and trans (labeled as θ₊), separated by two potential barriers, denoted here as cis-barrier and trans-barrier. The potential energy curves (PEC) obtained from the PES are then presented in Figure 2.

Figure 2

The PECs obtained for HP-Ng complexes are similar in shape and depth. The cis-barriers for the HP-Ng complexes are all smaller than the respective value for the free HP. The free HP has a cis-barrier of 7.5594 kcal/mol whereas the values for the complexes increase monotonically from 6.9828 kcal/mol for HP-Rn up to 7.5107 kcal/mol for HP-He. In addition the trans-barrier values are also lower than the respective value for the free HP, which is 1.0427 kcal/mol, and their values are 1.0928 kcal/mol for HP-He, 1.0817 kcal/mol for HP-Ne, 1.0651 kcal/mol for HP-Ar, 1.0676 kcal/mol for HP-Kr, 1.0736 kcal/mol for HP-Rn and 1.0749 kcal/mol for HP-Xe complexes. For the free HP, the cis-barrier and trans-barrier experimental energies (Hunt et al., 1965) are 7.0334 ± 0.0715 kcal/mol and 1.1036 ± 0.0114 kcal/mol, respectively. These values are in a good agreement with our results. However, for the HP-Ng complexes we did not find experimental data for comparison.

The results concerning geometric parameters, interactions and their characterization are summarized in Tables 1–3. The geometrical parameters obtained at MP2(full)/aug-cc-pVTZ, optimized without any constraints for two minimum structures and transition states (presented as cis and trans barriers) are given in Table 1 together with the graphical representation in Figure 3. Table 2 lists the binding energies corrected and uncorrected for the BSSE. Table 3 shows the numerical results for AIM analysis and Figure 4 depicts the ∇²ρ(r) contour plots for cis and trans barrier configurations.

Figure 3

Figure 4

The PESs yield θ₋ and θ₊ as true minima, i.e. without any imaginary frequencies in accordance with results from literature (Maciel et al., 2006; Roncaratti et al., 2014). In addition, all transition state structures displayed a well characterized imaginary frequency around 600 cm⁻¹ for cis and 400 cm⁻¹ for trans barriers (see Supplementary Information for further details). Figure 5 describes a schematic representation of vibrational modes of the isolated HP in the transition state with the actual frequencies and an imaginary frequency, which represents the frequency along the reaction coordinate.

Figure 5

For the HP-Ng complexes, the geometrical parameters are almost the same when compared with isolated HP in agreement with the weak interaction of these systems. The HP-Ng distances increase from He up Rn. On average, they are close to 2.55Å(He), 2.65Å(Ne), 2.75Å(Ar), 2.85Å(Kr), 3.00Å(Xe), and 3.06Å(Rn).

Regarding the binding energies, HP-He and HP-Ne are all repulsive, being less repulsive for the cis barrier configuration. This can be understood as a consequence of the fact that the noble gases turn out to be the hardest elements (Furtado et al., 2015) and this hardness decreases when the Ng atomic number is increased (the hardness in this context is a resistance to changes in its electronic population Furtado et al., 2015 coupled to Ng's high electronegativity Allen and Huheey, 1980). Although the BSSE increases monotonically from He to Rn, the binding energies also become more attractive as the atomic number increases.

For the four structures of each HP-Ng, the higher binding values are always observed for the cis-barrier configuration. As it will pointed latter, the decrease of the rate through the two barriers are not correlated with the binding energy, suggesting the hyperconjugation effects on HP may be important for the decrease of the interaction rate.

Regarding the AIM analysis, the existence of (3,−1) bond critical point (BCP) and its associated atomic interaction line indicates that electronic charge density is accumulated between the linked nuclei (Bader, 1991). In its turn, the values of the charge density ρ(r) in BCP are small while their corresponding ∇²ρ(r) are positive in accordance with a closed shell type of interaction. As a consequence, all configurations of all complexes show an interaction of a van der Waals type. Since higher ellipticity suggests conjugation and hyperconjugation effects of electron delocalization, these effects seem more pronounced in the HP-He and HP-Ne complexes. Another interesting feature is that all cis-barrier configurations of all complexes show a (3,+1) BCP indicating a cyclic nature.

3.0.2. Thermal Chiral Rate Analysis

The temperature dependence of the rate constant for cis to trans (i.e., through trans-barrier) and trans to cis (i.e., through cis-barrier) conformations for HP and HP-Ng complexes are presented in Figure 6. These results, in addition to conventional rate, are also exhibited with Eckart's or Wigner's tunneling corrections.

Figure 6

It was found that for the entire 100 K up to 4,000 K range the HP-Ne has the lowest rates for both barriers among all noble gas complexes, followed by HP-He. This result suggests that Ne and He are the noble gases more suitable for study the oriented collision dynamics with HP. In fact, the decrease of CRT shows an inverse correlation with respect the average valence electron energy (Allen, 1989), which follows the sequence (from higher to lower values): Ne, He, Ar, Kr; with Xe and Rn having very close values.

Nevertheless, there is a trend of rate increase as are move from Ar up to Rn. It is interesting to note that although this behavior is very similar regarding the cis-barrier for all rates (conventional, Wigner e Eckart), it seems that the tunneling is more important to describe the trans barrier's rate, where there is a significant difference for Eckart's values specially in the 100–200 K range when compared to respective Wigner and conventional results.

The final thermal rate constant can be expressed in the two familiar Arrhenius forms. In this work, the first is the Arrhenius modified form given by where A is the pre-exponential factor, T a temperature, n is a real number, R is the universal gas constant and E_a is the activation energy. The second is the d-Arrhenius form (Aquilanti et al., 2010; Silva et al., 2013; Carvalho-Silva et al., 2017) expressed by where d is a parameter that yield the degree of deformation of the exponential function.

The curve obtained by the reaction rate constant vs. the temperature can be fitted (Ramalho et al., 2011) to obtain the parameters A, n and E_a for the Arrhenius modified form, as presented in Table 4, and the parameters A, d, and E_a for d-Arrhenius form, as presented in Table 5. This feature confirms the trend of lower k(T) observed for HP-Ng complexes for both barriers when compared to isolated HP.

It can be also observed in Figure 6 that the trans to cis conformation rate of HP is lower (in the range 100–200 K) than the corresponding ones for HP-Ar, HP-Kr, HP-Xe, and HP-Rn. In the case of the chiral transition from cis to trans, the rates of all HP-Ng complexes are lower than that of the isolated HP. These results showed that the transition rate from cis to trans is greater than the corresponding trans to cis for both the isolated HP molecule and for all HP-Ng complexes. This suggests that the most important barrier that separates the chiral configurations of the isolated HP and the HP-Ng complexes is the trans-barrier, since it is the smallest. The energy of the HP's trans-barrier is relatively small (1.0427 kcal/mol) compared to its cis-barrier (7.5595 kcal/mol) as already seen in Figure 2.

An interesting result is presented in the Table 6. Although the increase in the trans-barrier of the HP-Ng complexes relative to HP is considerably small (see Figure 7), the change in the transition rate from cis to trans is relatively high. This is verified for high (4,000 K), room (298.15 K) and also for low temperatures (100 K). The most pronounced decrease in the rate corresponds to the HP-Ne complex, in which the decrease of the trans-barrier of just 0.0389 kcal/mol (see Figure 6) corresponds to a decrease of over 60% for the cis-trans transition rate, followed by HP-He. It is also interesting to note that this small change in energy barrier but with a substantial change in rate was also observed for other HP-Ng complexes. For example, HP-Ar complex showed a decrease of just 0.0224 kcal/mol but a 28.68% decrease of rate at 100 K.

Figure 7

Finally, at a temperature close to 300K, the Boltzmann distribution shows that about 16% of HP's population has higher energy than the trans-barrier with thermal fluctuations of approximately 1.7686 kcal/mol (Ball and Brindley, 2016). It has also been found that at low temperatures the chiral interconversion quantum encapsulation time of HP is very small. At a temperature of 100 K this time is <1 pico-second (Bitencourt et al., 2008), and at temperatures close to 0 K which can reach 3 pico-seconds.

4. Conclusions

The obtained results indicate that the chiral transition rate of trans to cis configuration of hydrogen peroxide in the presence of the noble gases He and Ne were the lowest over the entire temperature range of 4,000–100 K.

The AIM analysis shows that the interaction between H₂O₂ and the noble gases should be a van der Waals type. Although the H₂O₂ acts as an acid in the context of this investigation, the high hardness and high electronegativity of the nobles gases hold their electrons very tight to permit a covalence bond between H₂O₂ and Ng. On the other hand, it seems that both He and Ne are better able to affect the hyperconjugation effect and destabilizing repulsion among the lone pairs that are responsible for rotational barriers (Song et al., 2005). This may explain why the chiral transition rate decreases more for the complexes composed by Ne and He atoms, the hardest and more electronegative noble gases (Furtado et al., 2015).

Finally, the trans-barrier plays an important role because it is much smaller than the cis-barrier. The results showed that a small increase in the trans-barrier height in the complexes is responsible for a significant decrease in the rate of transition from cis to trans. Thus, these effects may contribute to the feasibility of separating one or the other enantiomer of the H₂O₂ molecule.

Acknowledgments

We gratefully acknowledge the financial support from the Brazilian Research Councils CNPq and FAPDF.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

• 1

  AllenL. C. (1989). Electronegativity is the average one-electron energy of the valence-shell electrons in ground-state free atoms. J. Am. Chem. Soc.111, 9003–9014. 10.1021/ja00207a003

  • CrossRef
  • Google Scholar

• 2

  AllenL. C.HuheeyJ. E. (1980). The definition of electronegativity and the chemistry of the noble gases. J. Inorg. Nuclear Chem.42, 1523–1524. 10.1016/0022-1902(80)80132-1

  • CrossRef
  • Google Scholar

• 3

  Al-RefaieA. F.OvsyannikovR. I.PolyanskyO. L.YurchenkoS. N.TennysonJ. (2015). A variationally calculated room temperature line-list for H₂O₂. J. Mol. Spectrosc.318, 84–90. 10.1016/j.jms.2015.10.004

  • CrossRef
  • Google Scholar

• 4

  AntunesF.BritoP. M. (2017). Quantitative biology of hydrogen peroxide signaling. Redox Biol.13, 1–7. 10.1016/j.redox.2017.04.039

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  AquilantiV.BartolomeiM.PiraniF.CappellettiD.VecchiocattiviF.ShimizuY.et al. (2005). Orienting and aligning molecules for stereochemistry and photodynamics. Phys. Chem. Chem. Phys.7, 291–300. 10.1039/B415212C

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  AquilantiV.MundimK. C.ElangoM.KleijnS.KasaiT. (2010). Temperature dependence of chemical and biophysical rate processes: Phenomenological approach to deviations from arrhenius law. Chem. Phys. Lett.498, 209–213. 10.1016/j.cplett.2010.08.035

  • CrossRef
  • Google Scholar

• 7

  AtkinsP.De PaulaJ.FriedmanR. (2013). Physical Chemistry: Quanta, Matter, and Change. New York, NY: Oxford University Press.

  • Google Scholar

• 8

  BaderR. F. (1991). A quantum theory of molecular structure and its applications. Chem. Rev.91, 893–928. 10.1021/cr00005a013

  • CrossRef
  • Google Scholar

• 9

  BallR.BrindleyJ. (2016). The life story of hydrogen peroxide iii: chirality and physical effects at the dawn of life. Orig. Life Evol. Biospheres46, 81–93. 10.1007/s11084-015-9465-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  BarretoP. R.PalazzettiF.GrossiG.LombardiA.MacielG.VilelaA. (2010). Range and strength of intermolecular forces for van der waals complexes of the type h2xn-rg, with x= o, s and n= 1, 2. Int. J. Quantum Chem.110, 777–786. 10.1002/qua.22127

  • CrossRef
  • Google Scholar

• 11

  BarretoP. R.VilelaA. F.LombardiA.MacielG. S.PalazzettiF.AquilantiV. (2007). The hydrogen peroxide- rare gas systems: Quantum chemical calculations and hyperspherical harmonic representation of the potential energy surface for atom- floppy molecule interactions. J. Phys. Chem. A111, 12754–12762. 10.1021/jp076268v

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  BellR. (1980). The Tunnel Effect in Chemistry.

  • Google Scholar

• 13

  BellR. P. (1959). The tunnel effect correction for parabolic potential barriers. Trans. Faraday Soc.55, 1–4. 10.1039/tf9595500001

  • CrossRef
  • Google Scholar

• 14

  BitencourtA. C.RagniM.MacielG. S.AquilantiV.PrudenteF. V. (2008). Level distributions, partition functions, and rates of chirality changing processes for the torsional mode around o–o bonds. J. Chem. Phys.129:154316. 10.1063/1.2992554

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  BoysS. F.BernardiF. d. (1970). The calculation of small molecular interactions by the differences of separate total energies. some procedures with reduced errors. Mol. Phys.19, 553–566. 10.1080/00268977000101561

  • CrossRef
  • Google Scholar

• 16

  CarlsonR.AndersonM.JohnsonR.SmytheW.HendrixA.BarthC.et al. (1999). Hydrogen peroxide on the surface of europa. Science283, 2062–2064. 10.1126/science.283.5410.2062

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  Carvalho-SilvaV. H.AquilantiV.de OliveiraH. C.MundimK. C. (2017). Deformed transition-state theory: Deviation from a rrhenius behavior and application to bimolecular hydrogen transfer reaction rates in the tunneling regime. J. Comput. Chem.38, 178–188. 10.1002/jcc.24529

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  Cortés-GuzmánF.BaderR. F. (2005). Complementarity of qtaim and mo theory in the study of bonding in donor–acceptor complexes. Coord. Chem. Rev.249, 633–662. 10.1016/j.ccr.2004.08.022

  • CrossRef
  • Google Scholar

• 19

  DazaM. C.DobadoJ.MolinaJ. M.VillavecesJ. L. (2000). Structure and bonding of h 2 o 2··· x complexes with (x= no+, cn-, hcn, hnc, co). Phys. Chem. Chem. Phys.2, 4089–4094. 10.1039/b001885f

  • CrossRef
  • Google Scholar

• 20

  DobadoJ.Martínez-GarcíaH.MolinaJ. M.SundbergM. R. (1998). Chemical bonding in hypervalent molecules revised. application of the atoms in molecules theory to y3 x and y3 xz (y= h or ch3; x= n, p or as; z= o or s) compounds. J. Am. Chem. Soc.120, 8461–8471. 10.1021/ja980141p

  • CrossRef
  • Google Scholar

• 21

  DobadoJ.MolinaJ. (1999). Adenine- hydrogen peroxide system: Dft and mp2 investigation. J. Phys. Chem. A103, 4755–4761. 10.1021/jp990671n

  • CrossRef
  • Google Scholar

• 22

  DobadoJ.MolinaJ. M. (1993). Ab initio molecular orbital calculation of the hydrogen peroxide dimer: study of basis set superposition error. J. Phys. Chem.97, 7499–7504. 10.1021/j100131a018

  • CrossRef
  • Google Scholar

• 23

  DunningT. H.Jr. (1989). Gaussian basis sets for use in correlated molecular calculations. i. the atoms boron through neon and hydrogen. J. Chem. Phys.90, 1007–1023. 10.1063/1.456153

  • CrossRef
  • Google Scholar

• 24

  EckartC. (1930). The penetration of a potential barrier by electrons. Phys. Rev.35:1303.

  • Google Scholar

• 25

  EncrenazT.BézardB.GreathouseT.RichterM.LacyJ.AtreyaS.et al. (2004). Hydrogen peroxide on mars: evidence for spatial and seasonal variations. Icarus170, 424–429. 10.1016/j.icarus.2004.05.008

  • CrossRef
  • Google Scholar

• 26

  EyringH. (1935). The activated complex in chemical reactions. J. Chem. Phys.3, 107–115. 10.1063/1.1749604

  • CrossRef
  • Google Scholar

• 27

  FehrensenB.LuckhausD.QuackM. (2007). Stereomutation dynamics in hydrogen peroxide. Chem. Physics, 338, 90–105. 10.1016/j.chemphys.2007.06.012

  • CrossRef
  • Google Scholar

• 28

  FrischM. J.TrucksG. W.SchlegelH. B.ScuseriaG. E.RobbM. A.CheesemanJ. R.et al. (2009) Gaussian 09 Revision e.01. Wallingford, CT: Gaussian Inc.

  • Google Scholar

• 29

  FurtadoJ.De ProftF.GeerlingsP. (2015). The noble gases: how their electronegativity and hardness determines their chemistry. J. Phys. Chem.119, 1339–1346. 10.1021/jp5098876

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  GoebelJ.AultB. S.Del BeneJ. E. (2000). Matrix isolation and ab initio study of the hydrogen-bonded complex between h2o2 and (ch3) 2o. J. Phys. Chem.104, 2033–2037. 10.1021/jp9941716

  • CrossRef
  • Google Scholar

• 31

  GoebelJ. R.AntleK. A.AultB. S.Del BeneJ. E. (2002). Matrix isolation and ab initio study of 1: 1 hydrogen-bonded complexes of h2o2 with hf, hcl, and hbr. J. Phys. Chem.106, 6406–6414. 10.1021/jp020520b

  • CrossRef
  • Google Scholar

• 32

  GoebelJ. R.AultB. S.Del BeneJ. E. (2001a). Matrix isolation and ab initio study of 1: 1 hydrogen-bonded complexes of h2o2 with nh3 and n (ch3) 3. J. Phys. Chem.105, 6430–6435. 10.1021/jp010813p

  • CrossRef
  • Google Scholar

• 33

  GoebelJ. R.AultB. S.Del BeneJ. E. (2001b). Matrix isolation and ab initio study of 1: 1 hydrogen-bonded complexes of h2o2 with phosphorus and sulfur bases. J. Phys. Chem.105, 11365–11370. 10.1021/jp013262

  • CrossRef
  • Google Scholar

• 34

  GonzálezL.MóO.YáñezM. (1997). High-level ab initio versus dft calculations on (h2o2) 2 and h2o2–h2o complexes as prototypes of multiple hydrogen bond systems. J. Comput. Chem.18, 1124–1135.

  • Google Scholar

• 35

  GrzechnikK.MierzwickiK.MielkeZ. (2013). Matrix-isolated hydrogen-bonded and van der waals complexes of hydrogen peroxide with ocs and cs2. ChemPhysChem14, 777–787. 10.1002/cphc.201200696

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  HenkelmanG.JóhannessonG.JónssonH. (2002). Methods for finding saddle points and minimum energy paths, in Theoretical Methods in Condensed Phase Chemistry, ed SchwartzS. D. (Dordrecht: Springer), 269–302.

  • Google Scholar

• 37

  HuntR. H.LeacockR. A.PetersC. W.HechtK. T. (1965). Internal-rotation in hydrogen peroxide: the far-infrared spectrum and the determination of the hindering potential. J. Chem. Phys.42, 1931–1946. 10.1063/1.1696228

  • CrossRef
  • Google Scholar

• 38

  KeithT. A. (2017). Aimall (version 17.11.04). TK gristmill software: Overland Park, KS. Available online at: aim.tkgristmill.com

  • Google Scholar

• 39

  LealL. A.da CunhaW. F.RoncarattiL. F.SilvaG. M. E.GarganoR. (2016). H₂O₂–ng dynamics predictions using an accurate potential energy surface. Mol. Phys.114, 440–445. 10.1080/00268976.2015.1078507

  • CrossRef
  • Google Scholar

• 40

  LewisW. M. (1967). Studies in catalysis. Part ix. the calculation in absolute measure of velocity constants and equilibrium constants in gaseous systems, in Selected Readings in Chemical Kinetics, eds BackM. H.LaidlerK. J. (Pergamon: Elsevier), 36–40.

  • Google Scholar

• 41

  LombardiA.PalazzettiF. (2018). Chirality in molecular collision dynamics. J. Phys. Condensed Matter30:063003. 10.1088/1361-648X/aaa1c8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  LombardiA.PalazzettiF.MacielG.AquilantiV.SevryukM. (2011). Simulation of oriented collision dynamics of simple chiral molecules. Int. J. Quantum Chem.111, 1651–1658. 10.1002/qua.22816

  • CrossRef
  • Google Scholar

• 43

  LundellJ.JolkkonenS.KhriachtchevL.PetterssonM.RäsänenM. (2001). Matrix isolation and ab initio study of the hydrogen-bonded h2o2–co complex. Chem. Eur. J. 7, 1670–1678. 10.1002/1521-3765(20010417)7:8<1670::AID-CHEM16700>3.0.CO;2-N

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  LundellJ.PehkonenS.PetterssonM.RäsänenM. (1998). Interaction between hydrogen peroxide and molecular nitrogen. Chem. Phys. Lett.286, 382–388. 10.1016/S0009-2614(98)00003-7

  • CrossRef
  • Google Scholar

• 45

  MacielG. S.BitencourtA. C. P.RagniM.AquilantiV. (2006). Studies of the dynamics around the o–o bond: orthogonal local modes of hydrogen peroxide. Chem. Phys. Lett.432, 383–390. 10.1016/j.cplett.2006.10.073

  • CrossRef
  • Google Scholar

• 46

  MałyszekP.KoputJ. (2013). Accurate ab initio potential energy surface and vibration-rotation energy levels of hydrogen peroxide. J. Comput. Chem.34, 337–345. 10.1002/jcc.23137

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  MattaC.BoydR. (2007). Quantum Theory of Atoms in Molecules: Recent Progress in Theory and Application. New York, NY: Wiley-VCH.

  • Google Scholar

• 48

  MoO.YanezM.RozasI.ElgueroJ. (1994). Structure, vibrational frequencies, and thermodynamic properties of hydrogen peroxide dimers: an ab initio molecular orbital study. J. Chem. Phys.100, 2871–2877. 10.1063/1.466429

  • CrossRef
  • Google Scholar

• 49

  MolinaJ. M.DobadoJ. A.DazaM. C.VillavecesJ. L. (2002). Structure and bonding of weak hydrogen peroxide complexes. J. Mol. Struct.580, 117–126. 10.1016/S0166-1280(01)00602-9

  • CrossRef
  • Google Scholar

• 50

  MuchaM.MielkeZ. (2009). Photochemistry of the glyoxal–hydrogen peroxide complexes in solid argon: formation of 2-hydroxy-2-hydroperoxyethanal. Chem. Phys. Lett.482, 87–92. 10.1016/j.cplett.2009.09.082

  • CrossRef
  • Google Scholar

• 51

  PalazzettiF.TsaiP.-Y.LombardiA.NakamuraM.CheD.-C.KasaiT.et al. (2013). Aligned molecules: chirality discrimination in photodissociation and in molecular dynamics. Rendiconti Lincei24, 299–308. 10.1007/s12210-013-0248-y

  • CrossRef
  • Google Scholar

• 52

  PehkonenS.LundellJ.KhriachtchevL.PetterssonM.RäsänenM. (2004). Matrix isolation and quantum chemical studies on the h 2 o 2–so 2 complex. Phys. Chem. Chem. Phys.6, 4607–4613. 10.1039/B410223A

  • CrossRef
  • Google Scholar

• 53

  PetersonK. A.FiggenD.GollE.StollH.DolgM. (2003). Systematically convergent basis sets with relativistic pseudopotentials. ii. small-core pseudopotentials and correlation consistent basis sets for the post-d group 16–18 elements. J. Chem. Phys.119, 11113–11123. 10.1063/1.1622924

  • CrossRef
  • Google Scholar

• 54

  PiatnytskyiD.ZdorevskyO.PerepelytsyaS.VolkovS. (2016). Formation of complexes of hydrogen peroxide molecules with dna. Ukr. J. Phys.61, 219–225. 10.15407/ujpe61.03.0219

  • CrossRef
  • Google Scholar

• 55

  PolanyiM.WignerE. (1928). Über die interferenz von eigenschwingungen als ursache von energieschwankungen und chemischer umsetzungen. Zeitschrift für Phys. Chem.139, 439–452. 10.1515/zpch-1928-13930

  • CrossRef
  • Google Scholar

• 56

  RamalhoS. S.da CunhaW. F.BarretoP. R.NetoP. H.RoncarattiL. F.e SilvaG. M.GarganoR. (2011). Thermal rate constant calculation of the nf+ f reactive system multiple arrangements. J. Phys. Chem.115, 8248–8254. 10.1021/jp204229d

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  RoncarattiL.LealL.PiraniF.AquilantiV.e SilvaG.GarganoR. (2014). Chirality of weakly bound complexes: the potential energy surfaces for the hydrogen-peroxide- noble-gas interactions. J. Chem. Phys.141:134309. 10.1063/1.4897136

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  RoncarattiL. F.AquilantiV. (2010). Whittaker–hill equation, ince polynomials, and molecular torsional modes. Int. J. Quantum Chem.110, 716–730. 10.1002/qua.22255

  • CrossRef
  • Google Scholar

• 59

  RoohiH.NowrooziA.BavafaS.AkbaryF.EshghiF. (2010). Interaction between nh2no and h2o2: a quantum chemistry study. Int. J. Quantum Chem.110, 1972–1981. 10.1002/qua.22364

  • CrossRef
  • Google Scholar

• 60

  ShiY.ZhouZ.-Y. (2004). Density functional theory study of the hydrogen bonding interaction complexes of hydrogen peroxide with glycine. J. Mol. Struct.674, 113–119. 10.1016/j.theochem.2003.12.018

  • CrossRef
  • Google Scholar

• 61

  SiesH. (2017). Hydrogen peroxide as a central redox signaling molecule in physiological oxidative stress: oxidative eustress. Redox Biol.11, 613–619. 10.1016/j.redox.2016.12.035

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  SilvaV. H.AquilantiV.de OliveiraH. C.MundimK. C. (2013). Uniform description of non-arrhenius temperature dependence of reaction rates, and a heuristic criterion for quantum tunneling vs classical non-extensive distribution. Chem. Phys. Lett.590, 201–207. 10.1016/j.cplett.2013.10.051

  • CrossRef
  • Google Scholar

• 63

  SongL.LiuM.WuW.ZhangQ.MoY. (2005). Origins of rotational barriers in hydrogen peroxide and hydrazine. J. Chem. Theory Comput.1, 394–402. 10.1021/ct049843x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  SuT.-M.PalazzettiF.LombardiA.GrossiG.AquilantiV. (2013). Molecular alignment and chirality in gaseous streams and vortices. Rendiconti Lincei24, 291–297. 10.1007/s12210-013-0249-x

  • CrossRef
  • Google Scholar

• 65

  TruhlarD. G.GarrettB. C.KlippensteinS. J. (1996). Current status of transition-state theory. J. Phys. Chem.100, 12771–12800. 10.1021/jp953748q

  • CrossRef
  • Google Scholar

• 66

  TruhlarD. G.IsaacsonA. D.GarrettB. C. (1985). Generalized Transition State Sheory, Volume 4.Boca Raton, FL: CRC Press.

  • Google Scholar

• 67

  WangC.ZhangD. H.SkodjeR. T. (2012). A six-dimensional wave packet study of the vibrational overtone induced decomposition of hydrogen peroxide. J. Chem. Phys.136:164314. 10.1063/1.4705755

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68

  WilsonA. K.WoonD. E.PetersonK. A.DunningT. H.Jr. (1999). Gaussian basis sets for use in correlated molecular calculations. ix. The atoms gallium through krypton. J. Chem. Phys.110, 7667–7676. 10.1063/1.478678

  • CrossRef
  • Google Scholar

• 69

  WoonD. E.Dunning JrT. H. (1993). Gaussian basis sets for use in correlated molecular calculations. iii. The atoms aluminum through argon. J. Chem. Phys.98, 1358–1371. 10.1063/1.464303

  • CrossRef
  • Google Scholar

• 70

  YuC.-Y.YangZ.-Z. (2011). A systemic investigation of hydrogen peroxide clusters (h2o2) n (n= 1- 6) and liquid-state hydrogen peroxide: based on atom-bond electronegativity equalization method fused into molecular mechanics and molecular dynamics. J. Phys. Chem.115, 2615–2626. 10.1021/jp111284t

  • Pubmed Abstract
  • CrossRef
  • Google Scholar