ABSTRACT

The influence of low-frequency sonolysis on the kinetics bimolecular liquid-phase reactions was studied with due regard for the association (dimers and trimers) of starting reagents. The mathematical modeling of chemical reactions that were described by nonlinear differential equations is performed. The steady states, the singular points characteristics, the nature of concentration oscillations in the reaction system are described. With increasing frequency and amplitude of low-frequency sonic waves (up to some critical value), we observed the cessation of the reaction. This observation offers an additional tool for controlling reaction rate by the external action of low-frequency vibrations. The conclusions of the work are obtained under certain assumptions. The exact determination of the critical conditions for changes in dynamics is beyond the scope of the problem.

Keywords:

Chemical Reaction, Kinetics, Liquid Medium, Reagent Association, Low-Frequency Action, Bistability Regime, Oscillation, Resonant Change

1. Introduction

When studying the kinetics of liquid-phase reactions, the dependence of the rate constants on the concentration of the reagent, the extreme nature of the Arrhenius dependence of the rate constant on temperature were found [1] [2]. These kinetic anomalies are not described by the kinetics of a homogeneous medium. The character of these anomalies is most strictly explained by the concept of the structural organization of the liquid during the formation of associates due to intermolecular interactions [2] [3] [4]. Theoretical and experimental data show that the association of molecules affects the reactivity of the molecules that form the associates [3] [4] [5] [6].

The associative structure of a liquid reaction system opens up new possibilities for controlling the kinetics of reactions by various kinds of external influences. Although many works are devoted to the mechanical effect on liquid-phase reactions, they usually study the phenomenological aspects of the effect of high-frequency acoustic fields on the reaction system [7] [8].

Data on the effect of the association of reagents on the kinetics of reactions are still ignored in the literature, and the theory of phenomena induced by periodic mechanical action, taking into account the supramolecular structure of the reaction medium, is being developed relatively recently.

Previously, we have proposed [2] [3] [4] that some anomalies in the kinetics of bimolecular liquid-phase reactions can be explained by the formation of the associates of starting reagents. Earlier, the effect of sonolysis on the reaction kinetics in liquids was explored mostly for high-frequency acoustic waves [7] [8]. Meanwhile, noticeable chemical effects were also reported for low-frequency (10 - 100 Hz) acoustic waves [9] [10] [11]. Deeper insight into the mechanism of this phenomenon can be reached via mathematical modeling of such reactions [12] [13]. The review aims at discussing possible approaches to regulate the kinetics of chemical reactions under the action low-frequency external action.

2. Mathematical Model

In this work, we modeled the influence of low-frequency sonolysis on the stability of intermediates in the А + В → С reactions in associated liquids [2], by using the reaction of alcohol A and isocyanate B yielding urethane С as an example [1] with the formation of associates of reactant A, namely, dimer A₂ and trimer A₃.

The association of the reactants even a formally one stage bimolecular reaction turns into a multistage process:

$\text{A}+\text{A}\leftrightarrow {\text{A}}_{\text{2}}\left(k_{\text{1}},k_{-1}\right)$ (1)

$\text{A}+{\text{A}}_{\text{2}}\leftrightarrow {\text{A}}_{\text{3}}\left(k_{\text{2}},k_{-2}\right)$ (2)

$\text{A}+{\text{A}}_{\text{2}}+{\text{A}}_{\text{3}}\to {\text{2A}}_{\text{3}}\left(k_{\text{3}}\right)$ (3)

$\text{В}+\text{B}\leftrightarrow {\text{B}}_{\text{2}}\left(k_{\text{4}},k_{-4}\right)$ (4)

${\text{B}}_{\text{2}}+{\text{A}}_{\text{3}}\to \text{С}+\text{B}+\text{2A}\left(k_{\text{5}}\right)$. (5)

Here, [A], [A₂] and [A₃] are the concentrations of monomer, dimer and trimer forms of alcoyol, respectively, [B] is the concentration of isocyanate, k_i and k_–i are the rate constants for the forward (i) and reverse (–i) reactions, respectively.

The kinetic equations for the intermediates of the system under consideration have the form:

$\begin{array}{c}\text{d}\left[{\text{A}}_{\text{2}}\right]/\text{d}t=k_{\text{1}}{\left[\text{A}\right]}^{\text{2}}-k{\left(v\right)}_{-\text{1}}\left[{\text{A}}_{\text{2}}\right]-k_{\text{2}}\left[\text{A}\right]\left[{\text{A}}_{\text{2}}\right]+k{\left(v\right)}_{-\text{2}}\left[{\text{A}}_{\text{3}}\right]\\ -k_{\text{3}}\left[\text{A}\right]\left[{\text{A}}_{\text{2}}\right]\left[{\text{A}}_{\text{3}}\right]\end{array}$ (6)

$\begin{array}{c}\text{d}\left[{\text{A}}_{\text{3}}\right]/\text{d}t=k_{\text{2}}\left[\text{A}\right]\left[{\text{A}}_{\text{2}}\right]+k_{\text{3}}\left[\text{A}\right]\left[{\text{A}}_{\text{2}}\right]\left[{\text{A}}_{\text{3}}\right]-k{\left(v\right)}_{-\text{2}}\left[{\text{A}}_{\text{3}}\right]\\ –k_{\text{4}}{k}_{\text{5}}{\left[\text{B}\right]}^{\text{2}}\left[{\text{A}}_{\text{3}}\right]{\left\{k_{-\text{4}}+k_{\text{5}}\left[{\text{A}}_{\text{3}}\right]\right\}}^{-\text{1}}\end{array}$.

According to the ideas of A. Tobolsky and G. Eyring, associates are fragile formations that are easily destroyed by a mechanical field [14].

The external action with ω = 2πv generates tensile stress σ that affects reaction rate constant

$k{\left(\omega \right)}_{-i}=k_{-0i}\exp \left[\left(-E_{-i}+{\gamma }_i\sigma \right)/RT\right]$,

where γ_i is a structure-sensitive coefficient and $\sigma ={Р}_0\sin \left(\omega t\right)$. Under steady-state conditions, the concentration of dimer B₂ and intermediates is described by the following system of parametric equations:

$\text{d}x/\text{d}t=\mu -a_{1}x\exp \left[{\gamma }_{1}P\sin \left(\omega t\right)\right]-a_{2}x+by\exp \left[{\gamma }_{2}P\sin \left(\omega t\right)\right]-uxy$

$\text{d}y/\text{d}t=ax-by\exp \left[{\gamma }_{2}P\sin \left(\omega t\right)\right]+uxy-my{\left(g+y\right)}^{-1}$, (7)

where x = [A₂], y = [A₃], μ = k₁ [A]², a₁ = k₋₁, a₂ = k₂ [A], b = k₋₂, u = k₃ [A], m = k₄ [B]², $g=k_{-4}{k}_5^{-1}$.

We also assume that a change in k (ν) during compression can be neglected in view of low compressibility of liquids. Moreover, the associates can be formed in case of some favorable mutual disposition of the reagents, which is infringed due to diffusion and formation of voids under the action of tension strain. With increasing separation between the reactive species, the reaction may transfer to a mode of diffusion control. Slow association at the stage of compression is known in the literature [5].

In order to escape cumbersome calculations, we solved the problem for a situation when [A] = [B] = const, In this case, the expression for rate w for formation of product С can be represented [3] [4] as:

$w=k_{\text{4}}{k}_{\text{5}}{\left[\text{B}\right]}^{\text{2}}\left[{\text{A}}_{\text{3}}\right]{\left\{k_{-\text{4}}+k_{\text{5}}\left[{\text{A}}_{\text{3}}\right]\right\}}^{-\text{1}}=my{\left(g+y\right)}^{-\text{1}}$. (8)

Accordingly, the expression for effective rate constant (k_eff) acquires the form:

$k_{\text{eff}}=k_{\text{4}}{k}_{\text{5}}{\left\{k_{-\text{4}}+k_{\text{5}}\left[{\text{A}}_{\text{3}}\right]\right\}}^{-\text{1}}\left[{\text{A}}_{\text{3}}\right]{\left[\text{A}\right]}^{-\text{1}}\left[\text{B}\right]$.

It follows that the dependence of w on ω manifests itself in an inexplicit way.

Equations (7) describes the dynamic behavior of a system with two degrees of freedom. The following situations are possible in terms of the concentrations of intermediates A₂ and A₃: the absence of a steady-state solution, the existence of only one steady state, and, finally, the existence of two steady states. Solving system of Equations (7) makes it possible to obtain the kinetic concentration dependence for x (t) and y (t) and to study the process of reaching a stationary (equilibrium) mode.

Let us introduce the following notation:

$\begin{array}{l}\text{d}x/\text{d}t=\Phi \left(x,y\right)\\ \text{d}y/\text{d}t=\Psi \left(x,y\right)\end{array}$. (9)

Under the condition that $\Phi \left(x,y\right)=\Psi \left(x,y\right)=0$, let us find a singular point (an equilibrium point) (x₀, y₀) in the x, y phase plane from the equation

$\text{d}y/\text{d}t=\Psi \left(x,y\right)/\Phi \left(x,y\right)$, (10)

and consider changes in the integral curves x(t) and y(t) as a motion along the phase plane. In this case, the problem of finding the properties of the singular points reduced to studying the stability of the motion at the point (x₀, y₀) using a known technique [1] [15].

For this purpose, consider Equation (9) with respect to the variables (ξ, η) by representing the variables x and y in the following forms:

$\begin{array}{l}x=x_0+\eta \\ y=y_0+\xi \end{array}$. (11)

By substituting Equations (11) into Equations (9) and expanding the functions Φ (x, y) and Ψ (x, y) into series, we obtain, in the first order of smallness, a linearized system:

$\begin{array}{l}\text{d}\eta /\text{d}t={\Phi }_x\eta +{\Phi }_y\xi \\ \text{d}\xi /\text{d}t={\Psi }_x\eta +{\Psi }_y\xi \end{array}$, (12)

in which ${\Phi }_x,{\Phi }_y,{\Psi }_x,{\Psi }_y$ are values of the partial derivatives of the functions Φ (x, y) and Ψ (x, y) by x and y, respectively, at the singular point (x₀, y₀). The behavior of the solutions of system Equations (12) near the considered singular point is determined by the properties of the system of linear differential equations with constant coefficients.

The standard solution (12) can be represented as a superposition of exponentials:

$\eta ={\Phi }_x{\text{e}}^{{\delta }_{1}t}+{\Phi }_y{\text{e}}^{{\delta }_{2}t}$,

$\xi ={\Psi }_x{\text{e}}^{{\delta }_{1}t}+{\Psi }_y{\text{e}}^{{\delta }_{2}t}$,

where δ_i (i = 1, 2) are the roots of the characteristic equation

${\delta }^2+p\delta +q=0$, (13)

in which $p={\Phi }_x+{\Psi }_y$,$q={\Phi }_x{\Psi }_y+{\Phi }_y{\Psi }_x$ and. In this case, the stability of the steady-state solution is determined by the coefficients p and q of Equations (13).

3. Results and Discussion

The problem was numerically solved for the case [A] = [B] = const (open system) by varying all the parameters of the system of Equations (7). Calculation below was carried out for model systems with the following parameters: [A₂]₀ = x (0) =1.8 M, and [A₃]₀ = y (0) = 0.23 M, μ = 35 M∙s⁻¹ and m = 200 M∙s⁻¹.

For the case without external action (γ₁ = γ₂ = 0) the phase portrait of the system was determined by solving equation (10) and the linearized system of Equations (12) and (13). On the phase plane, the integral curves x (t) and y (t) (Figure 1) have the form of concentric circles emerging from a singular point at х₀ = 2.51081 and у₀ = 0.33924, which is a stable center.

The concentrations x (t) and y (t) were calculated with the same parameters from the system of Equation (7). Unusual oscillations are observed in the system the maximum of the amplitudes changes harmonically (as is typical for a stable central mode) with a rather long oscillation period (Figure 2). However,

within each individual period of harmonic oscillations, there is a fine structure that is formed by oscillations with a smaller period.

The concentrations x (t) and y (t) were calculated from (7) for different amplitudes P and frequencies ω and exhibited similar behavior. Essential change in the structure of oscillations happens at ω₀ = 0.4 s⁻¹ and for P > 2. With increasing w, functions x (t) and y (t) acquire rather complicate profiles (Figure 3) or exceedingly large spread (Figure 4), which can be explained by the presence [A₃] and high formation rate for associates B₂. Increase in w leads to the formation of unstable focuses.

For [A] = [B] = const, we obtained the kinetics of change in the concentrations of associates А₂ and А₃ as a function of amplitude Pand frequency ω at starting stable-node stability. Two critical ω values were found to occur. For low P and ω, we observed oscillations in the concentration of the associates without change in the type of equilibrium. An increase in P and/or ω results in hybrid oscillations. For low ω, P_min is close to zero. With increasing, we come to a second critical phenomenon: bifurcation of the system state. At the bifurcation point, the oscillations are insignificant but, with increasing ω, two kinds of harmonic oscillations with close periods but different amplitudes and concomitant pulsations come into being (see Figure 3).

A reasonable question arises on the applicability of the results obtained for open systems to the closed ones which are most frequent in experimental. Strictly speaking, reagent concentrations in closed systems remain unsteady only during a certain period of reaction during which the concentrations of intermediates may be considered as quasi-stationary. It has been shown experimentally that the dissipation structures can also form in closed systems as well [14].

4. Conclusions

As far as we know, this work made almost the first attempt to generalize data on the kinetics of a complex process in a liquid medium during the association of reagents.

The above results suggest that sonolysis can be used as an additional tool for controlling the liquid-phase reactions involving the associates of starting reagents. The development of ideas about the structural organization of reagents in the liquid phase, about its effect on chemical reactions, about changes in the structural organization of reagents due to external influences of a different nature, allows a new approach to solve applied problems of chemistry, biochemistry, and medicine.

Acknowledgements

This work was financially supported by FASO (grant No. 0089-2019-0002).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Kulagina, T.P., Smirnov, L.P. and Andrianova, Z.S. (2019) Influence of Low-Frequency Sonolysis on the Stability of Reactions Involving Associates. Journal of Applied Mathematics and Physics, 7, 2950-2957. https://doi.org/10.4236/jamp.2019.712205

References