Modeling Chemical Reactivity at the Interfaces of Emulsions: Effects of Partitioning and Temperature

Bulk phase chemistry is hardly ever a reasonable approximation to interpret chemical reactivity in compartmentalized systems, because multiphasic systems may alter the course of chemical reactions by modifying the local concentrations and orientations of reactants and by modifying their physical properties (acid-base equilibria, redox potentials, etc.), making them—or inducing them—to react in a selective manner. Exploiting multiphasic systems as beneficial reaction media requires an understanding of their effects on chemical reactivity. Chemical reactions in multiphasic systems follow the same laws as in bulk solution, and the measured or observed rate constant of bimolecular reactions can be expressed, under dynamic equilibrium conditions, in terms of the product of the rate constant and of the concentrations of reactants. In emulsions, reactants distribute between the oil, water, and interfacial regions according to their polarity. However, determining the distributions of reactive components in intact emulsions is arduous because it is physically impossible to separate the interfacial region from the oil and aqueous ones without disrupting the existing equilibria and, therefore, need to be determined in the intact emulsions. The challenge is, thus, to develop models to correctly interpret chemical reactivity. Here, we will review the application of the pseudophase kinetic model to emulsions, which allows us to model chemical reactivity under a variety of experimental conditions and, by carrying out an appropriate kinetic analysis, will provide important kineticparameters.


Introduction
The three-dimensional interfaces between two immiscible liquids (oil and water) found in emulsions and other colloidal systems have been widely used as models of membrane function to mimic their behavior, as models of lipid-based foods and to model important (bio)chemical reactions such as the oxidation of lipids [1][2][3][4]. A variety of fundamental processes involved in biocatalysis, such as ion pumping, electron transport, membrane fusion, and photosynthesis, have all been investigated in such interfacial systems for further guidance in developing new biomaterials and medicines and to better understand life science [5][6][7][8][9][10][11].
Emulsions are composed of a liquid dispersed as small droplets in a second immiscible liquid. The droplet interface is encompassed by a narrow region (2-20 nm thick) surrounding each emulsion droplet, Figure 1 [12]. The interfacial region typically does not contribute significantly to the total volume of an emulsion unless the droplet sizes are very small [13], but has important consequences in the fate of chemical reactions that take place in the emulsions. For example, lipid oxidation reactions and those between antioxidants and lipid radicals are largely affected by the presence of interfacial regions, because they are highly anisotropic regions composed of a mélange of oil, water, and surfactants [2,[14][15][16]. Conceptually, emulsions can be divided into three regions with distinct solvent properties: the core of droplets, the continuous phase, and the interfacial region ( Figure  1). The interfacial layer surrounding each emulsion droplet has a crucial role in the applicability and properties of emulsions because many important chemical, biological, physical, and technological processes take place at interfaces [16][17][18][19]. Reactive components can partition between the different regions, reducing or increasing the concentrations of these components in the site of reaction. The differential solubility of reactants (and of reaction products) in the various regions of the system can alter their local concentrations, making possible to carry out reactions at rates much higher (or much lower) than those in bulk homogeneous phases, allowing the tuning of reaction rates. It is not strange, therefore, that the interfaces between two immiscible liquids can serve as convenient models for investigating the fate of chemical reactions in more complex natural systems such as cell membranes [20][21][22][23].
To interpret chemical reactivity in emulsions, a kinetic model based on the formalism of the pseudo-phase was devised [19,24,25]. In this article, we will summarize its main features and applications. We will also show some illustrative examples that highlight the importance of reactant partitioning and of the main outcomes of the model, for example, determining activation parameters of reactions taking place in the interfacial region of emulsions. The model was originally developed aiming at solving an important, but refractory, problem in food science: how to quantify the effective concentration of antioxidants added to slow lipid oxidation in the oil, water, and interfacial regions of food grade, lipid-based emulsions.
In the pseudo-phase model, the whole solution is divided into oil, interfacial, and aqueous regions (pseudo-phases), with the surfactant located in the boundary between the oil and water regions [19,24,[26][27][28][29][30]. Each region has different solvent properties and is treated as a separate phase or pseudo-phase. Reactants (and other components) partition between the regions thermodynamically, depending on the values of the Gibbs free energy of transfer. This implies that the distribution of the reagents is always described by partition constants whose value is the same throughout the course of the reaction.
However, some complications may arise when considering thick emulsions where transport of reactants between regions can be rate-limiting [13,31]. For the sake of simplicity, the review will focus exclusively on fluid emulsions where restrictions to the free movement of reactants and other components, if any, do not affect molecular transport. Conceptually, emulsions can be divided into three regions with distinct solvent properties: the core of droplets, the continuous phase, and the interfacial region ( Figure 1). The interfacial layer surrounding each emulsion droplet has a crucial role in the applicability and properties of emulsions because many important chemical, biological, physical, and technological processes take place at interfaces [16][17][18][19]. Reactive components can partition between the different regions, reducing or increasing the concentrations of these components in the site of reaction. The differential solubility of reactants (and of reaction products) in the various regions of the system can alter their local concentrations, making possible to carry out reactions at rates much higher (or much lower) than those in bulk homogeneous phases, allowing the tuning of reaction rates. It is not strange, therefore, that the interfaces between two immiscible liquids can serve as convenient models for investigating the fate of chemical reactions in more complex natural systems such as cell membranes [20][21][22][23].
To interpret chemical reactivity in emulsions, a kinetic model based on the formalism of the pseudo-phase was devised [19,24,25]. In this article, we will summarize its main features and applications. We will also show some illustrative examples that highlight the importance of reactant partitioning and of the main outcomes of the model, for example, determining activation parameters of reactions taking place in the interfacial region of emulsions. The model was originally developed aiming at solving an important, but refractory, problem in food science: how to quantify the effective concentration of antioxidants added to slow lipid oxidation in the oil, water, and interfacial regions of food grade, lipid-based emulsions.
In the pseudo-phase model, the whole solution is divided into oil, interfacial, and aqueous regions (pseudo-phases), with the surfactant located in the boundary between the oil and water regions [19,24,[26][27][28][29][30]. Each region has different solvent properties and is treated as a separate phase or pseudo-phase. Reactants (and other components) partition between the regions thermodynamically, depending on the values of the Gibbs free energy of transfer. This implies that the distribution of the reagents is always described by partition constants whose value is the same throughout the course of the reaction.
However, some complications may arise when considering thick emulsions where transport of reactants between regions can be rate-limiting [13,31]. For the sake of simplicity, the review will focus exclusively on fluid emulsions where restrictions to the free movement of reactants and other components, if any, do not affect molecular transport. Prior to considering modeling the chemical reactivity, we will briefly discuss some properties of emulsions and their dynamic aspects.

Emulsions: Classification, Physical Stabilization, and Main Properties
Emulsions are mixtures of two immiscible liquids such as oil and water, where one of the liquids is dispersed as small spherical droplets in the other one. They can be prepared, Scheme 1, by either shaking together the oil and water or by adding one of the liquids drop by drop to the other liquid with agitation. However, the surface of each droplet is an interface between hydrophilic and hydrophobic molecules. The contact between the molecules of the two immiscible liquids is highly unfavorable because of the increased enthalpy caused by the contact between hydrophilic and hydrophobic molecules, which is not compensated by the increased entropy of mixing. As a result, emulsions are thermodynamically unstable colloidal systems, which tend to separate the two liquids until the contact area between them is the lowest possible (i.e., minimal Gibbs free energy). Prior to considering modeling the chemical reactivity, we will briefly discuss some properties of emulsions and their dynamic aspects.

Emulsions: Classification, Physical Stabilization, and Main Properties
Emulsions are mixtures of two immiscible liquids such as oil and water, where one of the liquids is dispersed as small spherical droplets in the other one. They can be prepared, Scheme 1, by either shaking together the oil and water or by adding one of the liquids drop by drop to the other liquid with agitation. However, the surface of each droplet is an interface between hydrophilic and hydrophobic molecules. The contact between the molecules of the two immiscible liquids is highly unfavorable because of the increased enthalpy caused by the contact between hydrophilic and hydrophobic molecules, which is not compensated by the increased entropy of mixing. As a result, emulsions are thermodynamically unstable colloidal systems, which tend to separate the two liquids until the contact area between them is the lowest possible (i.e., minimal Gibbs free energy). Scheme 1. Formation of emulsions after mixing oil, an aqueous solution, and the appropriate surfactant with some agitation. Addition of surfactants is needed to stabilize kinetically the droplets.
To stabilize kinetically emulsions, surfactant molecules are added to the binary oilwater system. Surfactants are amphiphilic molecules having both hydrophilic, water-soluble and lipophilic, oil soluble residues that interact simultaneously with the oil and water, so that they locate at the droplet surface with their hydrophilic part of the emulsifier pointing to the water and its lipophilic part pointing to the oil, forming a three-dimensional boundary surrounding the surface of the droplets, Figure 1. This arrangement makes emulsions be stabilized kinetically during periods of time ranging from minutes to years, depending on the composition of the emulsion. Thus, an appropriate choice of the surfactant is crucial in the formation of the emulsion and its long-term stability. Surfactants can be anionic, cationic, non-ionic, and zwitterionic, Scheme 2, and they play a major role on the formation, stability, and functional characteristics of the emulsion. Scheme 1. Formation of emulsions after mixing oil, an aqueous solution, and the appropriate surfactant with some agitation. Addition of surfactants is needed to stabilize kinetically the droplets.
To stabilize kinetically emulsions, surfactant molecules are added to the binary oilwater system. Surfactants are amphiphilic molecules having both hydrophilic, watersoluble and lipophilic, oil soluble residues that interact simultaneously with the oil and water, so that they locate at the droplet surface with their hydrophilic part of the emulsifier pointing to the water and its lipophilic part pointing to the oil, forming a three-dimensional boundary surrounding the surface of the droplets, Figure 1. This arrangement makes emulsions be stabilized kinetically during periods of time ranging from minutes to years, depending on the composition of the emulsion. Thus, an appropriate choice of the surfactant is crucial in the formation of the emulsion and its long-term stability. Surfactants can be anionic, cationic, non-ionic, and zwitterionic, Scheme 2, and they play a major role on the formation, stability, and functional characteristics of the emulsion. A valuable tool for classifying surfactants and their main applications involves the so-called hydrophilic lipophilic balance (HLB) value. The HLB of an emulsifier represents a balance between the hydrophilic (polar) portion of the surfactant and the lipophilic (non-polar) groups of the emulsifier. The more lipophilic surfactants have low HLB val-Scheme 2. Chemical structures of some common surfactants employed to stabilize kinetically emulsions. A valuable tool for classifying surfactants and their main applications involves the so-called hydrophilic lipophilic balance (HLB) value. The HLB of an emulsifier represents a balance between the hydrophilic (polar) portion of the surfactant and the lipophilic (non-polar) groups of the emulsifier. The more lipophilic surfactants have low HLB values (<9.0); meanwhile, the more hydrophilic ones have a high HLB number (>11.0). When two or more surfactants are blended, the resulting HLB of the blend is easily calculated taking into consideration the HLB of each surfactant and its relative proportion. Table 1 illustrates the range of HLB values that are most suited for a particular application.  [32,33].

<3
Surface films 3-6 Water-in-oil emulsions 7-9 Wetting agents 8-15 Oil-in-water emulsions 13-15  Emulsions are thermodynamically unstable systems and tend to phase separation with time. Scheme 3 shows the main emulsion destabilization mechanisms of emulsions and nanoemulsions. Nanoemulsions are, in general, more resistant to physical destabilization via creaming, sedimentation, coalescence, and/or flocculation than those with larger droplets, in part, because of the efficient steric stabilization provided by the surfactant [35]. In addition, Brownian motion usually overcomes gravitational forces in emulsions with small droplets. On the contrary, nanoemulsions are mostly susceptible to destabilization through the Ostwald ripening mechanism [9,32,35]. Emulsions are thermodynamically unstable systems and tend to phase separation with time. Scheme 3 shows the main emulsion destabilization mechanisms of emulsions and nanoemulsions. Nanoemulsions are, in general, more resistant to physical destabilization via creaming, sedimentation, coalescence, and/or flocculation than those with larger droplets, in part, because of the efficient steric stabilization provided by the surfactant [35]. In addition, Brownian motion usually overcomes gravitational forces in emulsions with small droplets. On the contrary, nanoemulsions are mostly susceptible to destabilization through the Ostwald ripening mechanism [9,32,35].
Emulsions are thermodynamically unstable systems and tend to phase separation with time. Scheme 3 shows the main emulsion destabilization mechanisms of emulsions and nanoemulsions. Nanoemulsions are, in general, more resistant to physical destabilization via creaming, sedimentation, coalescence, and/or flocculation than those with larger droplets, in part, because of the efficient steric stabilization provided by the surfactant [35]. In addition, Brownian motion usually overcomes gravitational forces in emulsions with small droplets. On the contrary, nanoemulsions are mostly susceptible to destabilization through the Ostwald ripening mechanism [9,32,35]. Conceptually, emulsions can be divided into three regions with distinct solvent properties: the core of droplets, the continuous phase, and the interfacial region ( Figure  1). The interfacial layer surrounding each emulsion droplet has a crucial role in the applicability and properties of emulsions because many important chemical, biological, physical, and technological processes take place at interfaces [16][17][18][19]. The properties of the oil droplet surface control, to a large extent, the physical stability of emulsions. The nature and concentration of surfactants, together with the composition of the aqueous phase (pH, ionic strength, and ions) determine the electrostatic charge and the thickness of the interfacial layer, both parameters being deeply involved in the physical stability of emulsions [13]. However, the interfacial layer also has an important role in the fate of chemical reactions because reactants in emulsions are compartmentalized [19]. Polar molecules distribute preferentially, between the aqueous and interfacial regions; meanwhile, nonpolar molecules distribute preferentially between the oil and interfacial regions. Molecules of moderate hydrophobicity are distributed between the three regions and those that are both oil and water insoluble are exclusively located in the interfacial region. Conceptually, emulsions can be divided into three regions with distinct solvent properties: the core of droplets, the continuous phase, and the interfacial region ( Figure 1). The interfacial layer surrounding each emulsion droplet has a crucial role in the applicability and properties of emulsions because many important chemical, biological, physical, and technological processes take place at interfaces [16][17][18][19]. The properties of the oil droplet surface control, to a large extent, the physical stability of emulsions. The nature and concentration of surfactants, together with the composition of the aqueous phase (pH, ionic strength, and ions) determine the electrostatic charge and the thickness of the interfacial layer, both parameters being deeply involved in the physical stability of emulsions [13]. However, the interfacial layer also has an important role in the fate of chemical reactions because reactants in emulsions are compartmentalized [19]. Polar molecules distribute preferentially, between the aqueous and interfacial regions; meanwhile, nonpolar molecules distribute preferentially between the oil and interfacial regions. Molecules of moderate hydrophobicity are distributed between the three regions and those that are both oil and water insoluble are exclusively located in the interfacial region.

Dynamic Aspects of Emulsions and Mass Transfer
Emulsions are highly dynamic systems and diffusion within droplets, exchange of material between droplets, and self-diffusion of droplets in the continuous phase are common dynamic processes having different time scales, which may affect their potential use [31,[37][38][39][40]. Molecular diffusion in emulsions has attracted considerable interest to control the macroscopic physicochemical properties of the system. For example, Ostwald ripening involves the transfer of oil through the aqueous phase, while coalescence involves oil transfer achieved by fusion of the emulsifier film (a droplet-droplet direct interaction) [41]. Interfacial adsorbed surfactants are also dynamic, experiencing lateral and transversal diffusion processes, undergoing continuous reorganization, and are transferred between emulsion droplets [20,42,43]. Surfactants are able to form micelles in the water phase, and they may form swollen micelles (which are several nm in size) that contain oil molecules. These swollen micelles may carry small fractions of oil from small to big droplets and may contribute to the ripening of the emulsions.
Molecular diffusion in emulsions is also relevant to control the kinetic processes of the various chemical reactions that take place at the liquid-liquid interfacial boundaries because emulsion droplets are broken and reorganized constantly, being oil, surfactants, and other emulsion components (rapidly) transferred between droplets [39,44,45]. Diffusion of reactive molecules within or between emulsion droplets depends not only on their molecular nature but also on the particular properties of the components of the emulsions (composition, polarity and viscosity of the oil, components of the water phase, nature of the emulsifier, etc.). Diffusion can be either faster or slower than the undergoing chemical reaction. If the observed rate depends on the rate of molecular diffusion, the reactions are said to be under microscopic diffusion control; meanwhile, if the observed rate depends on the rate at which solutions mix, they are said to be macroscopic diffusion control. When the reaction is slow in comparison with the emulsion dynamics, the reaction "sees" the emulsion droplets as static entities, so that the observed rate of the reaction can be computed as the sum of the rates in all regions of the emulsion. On the contrary, when the diffusion of reactants is rate-limiting, the dynamics of the droplets needs to be taken into account to explain the fate of the chemical reactions.

Kinetics of Chemical Reactions at Interfaces: Diffusive vs. Reactive Systems
Most kinetic theories have been mainly developed to describe the fate of reactions taking place in homogeneous systems, where ideal conditions are usually assumed. However, emulsions are complex systems containing oils, aqueous solutions (where ionic and nonionic compounds can be dissolved in different extents), and an interfacial region connecting the oil and water phases created upon addition of surfactants. In emulsions, several processes occurring on a variety of time scales take place simultaneously. For example, emulsion droplets break and reform continuously, and emulsion components such as oil, water, surfactants, antioxidants, lipid radicals, or other reactive molecules exchange between droplets. Thus, the rate of the chemical reactions that takes place at the interfacial region may depend on the rate of transport of emulsion components between the three regions of the emulsion and eventually may be rate-limiting [19,24,37].
To predict whether the dynamics of the droplets affect the chemical reaction, a practical approach, based on the Smoluchowski theory for Brownian motion [46] and the application of the Acree-Curtin-Hammet principle [47], was developed by Bravo et al. [48]. The approach only considers "ideal" systems, where there are no physical barriers limiting the transport of molecules between the oil, aqueous, and interfacial regions and, therefore, where there are no significant deviations from simple diffusion laws, molecular mobility, partitioning, and volume exclusion effects, etc. [48].
The Acree-Curtin-Hammett (ACH) principle [47,49] is a well-known concept described in many physical organic and kinetics textbooks commonly used to explain the selectivity ratios for stereoselective reactions such as that depicted in Scheme 4A. Briefly, the ACH principle states that: "if the magnitude of the rates involving the equilibrium process of A and B (k1, k − 1) are much higher than those rates of the chemical reactions leading to the formation of products (kA, kB), then the system is under kinetic control, and the final product ratio is equal to the ratio of the product forming rate constants, Equation (1), where K is the equilibrium constant (K = k1/k − 1)". On the contrary, if the rates of diffusion are lower than the rates of the chemical reaction, then the reaction is diffusion-controlled and the product distribution depends, in a complex manner, on both the relative size of the diffusion rate constants and the chemical rate constants. [ The Acree-Curtin-Hammett ideas can be extended to settings where the reactants partition between regions of different solvent properties as depicted in Scheme 4B [19,48]. Only when the rates of transport of molecules (i.e., the rates of diffusion) between regions (k1, k − 1, k2, k − 2, k3, k − 3, k4, k − 4) are much faster than the rates of the chemical reactions (kW, kI, kO), the overall rate of the reaction is the sum of the rates of the reactions in every phase. Under these conditions, the system is said to be in "dynamic equilibrium" throughout the course of the chemical reaction. On the contrary, when the diffusion of reactants is rate-limiting, the dynamics of the droplets needs to be taken into account to explain the fate of the chemical reactions [50][51][52][53], and the product distribution depends, in a complex manner, on both the relative size of the diffusion rate constants and the chemical rate constants. constant (K = k1/k−1)". On the contrary, if the rates of diffusion are lower than the rates of the chemical reaction, then the reaction is diffusion-controlled and the product distribution depends, in a complex manner, on both the relative size of the diffusion rate constants and the chemical rate constants.
[P ] The Acree-Curtin-Hammett ideas can be extended to settings where the reactants partition between regions of different solvent properties as depicted in Scheme 4B [19,48]. Only when the rates of transport of molecules (i.e., the rates of diffusion) between regions (k1, k−1, k2, k−2, k3, k−3, k4, k−4) are much faster than the rates of the chemical reactions (kW, kI, kO), the overall rate of the reaction is the sum of the rates of the reactions in every phase. Under these conditions, the system is said to be in "dynamic equilibrium" throughout the course of the chemical reaction. On the contrary, when the diffusion of reactants is rate-limiting, the dynamics of the droplets needs to be taken into account to explain the fate of the chemical reactions [50][51][52][53], and the product distribution depends, in a complex manner, on both the relative size of the diffusion rate constants and the chemical rate constants. A simple mathematical model, grounded on the transition state theory (TST) and diffusion laws, permits us to analyze the conditions required to distinguish between kinetically controlled and diffusion-controlled reactive systems based on molecular and solvent properties. For the sake of clarity, only the main assumptions and conclusions of the model will be presented and the interested reader is referred to specialized publications [46,48,50,51].
The model considers two neutral reactive, spherical, molecules A and B, having radii rA and rB, respectively, in a fluid of a given viscosity (for example, in water) with an energy enough to overcome the activation barrier. The model assumes that the observed rate does not depend on the rate at which the reactants are mixed (i.e., mixing effects are unimportant), that reactants are uniformly distributed, and that the orientation and steric A simple mathematical model, grounded on the transition state theory (TST) and diffusion laws, permits us to analyze the conditions required to distinguish between kinetically controlled and diffusion-controlled reactive systems based on molecular and solvent properties. For the sake of clarity, only the main assumptions and conclusions of the model will be presented and the interested reader is referred to specialized publications [46,48,50,51].
The model considers two neutral reactive, spherical, molecules A and B, having radii r A and r B , respectively, in a fluid of a given viscosity (for example, in water) with an energy enough to overcome the activation barrier. The model assumes that the observed rate does not depend on the rate at which the reactants are mixed (i.e., mixing effects are unimportant), that reactants are uniformly distributed, and that the orientation and steric effects are unimportant. Therefore, the estimated rate constants will be an upper limit because if any of the assumptions does not apply, then the true rate constant will be a fraction of that estimated.
If the reactive molecules react immediately when their distance is equal to, or less than, the sum of their radii (i.e., upon an encounter), the reaction is under full microscopic diffusion control, and an approximate value for the rate of the reaction between neutral molecules can be found by employing Fick s law of diffusion. Since the reactive molecules are usually much larger in size than the solvent molecules, Stoke's law, Equation (2), applies, and Equation (3) can be derived [50,54]. In Equations (2) and (3), D is the diffusion coefficient, k B is the Boltzmann's constant, T is the temperature, η is the dynamic viscosity of the media, and r the hydrodynamic radius of the reactive molecules.
In water, at T = 25 • C, the viscosity is η~0.9 cP s, and for molecules of similar size (r A ≈ r B ), a value of k D ≈ 7 × 10 9 dm 3 mol −1 s −1 can be estimated for the bimolecular reaction. This k D value represents an upper limit value for diffusion-controlled reactions of neutral molecules. That is, neutral molecules cannot react faster than the rate at which they diffuse in the solvent. If reactive ions are involved, the value is expected to be somewhat higher because the electrostatic interactions increase their diffusivity. For example, values of k D ranging (3 − 15) × 10 10 dm 3 mol −1 s −1 have been estimated for the reaction between hydrogen ions and various anions [50].
If the chemical process is not fast enough to remove the reactive molecules completely as they reach the critical distance r AB , then the observed rate constant, k obs , depends on the relative values of k D and the rate constant of the chemical process, k chem , Equation (4). Note that Equation (4) reduces to Equation (3) when k chem >>> 4π(D A + D B )r AB ; meanwhile, when k chem <<< 4π(D A + D B )r AB , k obs = k chem (i.e., no diffusion control).
An analysis of the variation of k obs with k chem for a given value of k D , Equation (4), allows us to predict if the observed rate constant is, or is not, affected by the diffusion of reactants and in what extent it is. For the purpose, we can plot the variation of k obs with k chem for a given upper limit value of k D . Thus, by setting the previously estimated value of k D = 7 × 10 9 dm 3 mol −1 s −1 (T = 25 • C) for a bimolecular reaction in aqueous solution, the variations of k obs with log(k chem ), Figure 3, can be determined. As shown, diffusion-controlled reactions (k obs ≈ k D ) are those whose rate of the chemical reaction is k chem ≥ 10 12 M −1 s −1 . Therefore, for kinetically controlled reactions, k obs ≈ k chem and k chem should be <10 9 dm 3 mol −1 s −1 . For intermediate k chem values (10 9 -10 11 dm 3 mol −1 s −1 ), the reaction is controlled by both diffusion and activation in a different extent, depending on the particular k chem value. For example, when k chem~kD (=7 × 10 9 dm 3 mol −1 s −1 ), k obs ≈ 3 × 10 9 , i.e., approximately 50% of diffusion control.  The model also considers the role of the activation energy of the chemical reaction that, according to the TST theory, largely affects the chemical reaction rate. Comparison of the kchem values in aqueous solution, Table 2, with the diffusion rate values estimated in Figure 3 indicates that for any chemical reaction with Ea higher than 25 kJ mol −1 , kobs = kchem, i.e., for most thermally activated reactions that may take pace in emulsions, the reactants are distributed according to their solubilities in the different regions of the emulsion before they react, and the rates of diffusion of the reactants are much higher than those of the chemical reactions and the system will be under dynamic equilibrium.
The model also considers the role of the activation energy of the chemical reaction that, according to the TST theory, largely affects the chemical reaction rate. Comparison of the k chem values in aqueous solution, Table 2, with the diffusion rate values estimated in Figure 3 indicates that for any chemical reaction with E a higher than~25 kJ mol −1 , k obs = k chem , i.e., for most thermally activated reactions that may take pace in emulsions, the reactants are distributed according to their solubilities in the different regions of the emulsion before they react, and the rates of diffusion of the reactants are much higher than those of the chemical reactions and the system will be under dynamic equilibrium. For extremely rapid reactions, for example, radical reactions where there are shortlived radicals or photochemical reactions where excited-state molecules may play a significant role, the reactions may be diffusion-controlled, as their activation energies are <10 kJ mol −1 , and the rate of the reaction depends on the dynamics of fast processes (proton transfer, solvation of the probe, solute, or counterion entry and exit from the aggregates) that take place in the nanoseconds to microsecond scales [37,55,56]. Some of these reactions, for instance, the bimolecular quenching of the excited state of the probes by other molecules present in the medium, have been used to determine the mobilities of molecules, to estimate local microviscosities and encounter probabilities in the medium, and to investigate ultrafast proton transfer reactions [37,56].
In summary, the rates of common thermally activated reactions are much lower than diffusion, and unless there are physical restrictions to movement (for example, if polymer networks are located in the interfacial region or if electrostatic forces are relevant), diffusion is not rate-limiting. Both diffusion-controlled and activation-controlled reactions are important and commonly observed in nature. Examples of diffusion-controlled reactions include heterogeneous catalyzes, cell metabolism, gas diffusion through solid, polymer chain growth kinetics, colloid or crystal growth, precipitation, fluorescence quenching and other photochemical processes, some electrochemical reactions, combustion, etc. In general, diffusion-controlled reactions are those with low activation energies, so that the reaction time is negligible in comparison with that required for the diffusion of the molecules to collide. On the contrary, if the reaction times are much higher than those required for the diffusion of the reactants (high activation energies), the reactions are controlled by activation.
Our research on antioxidants in fluid emulsions over the past decade demonstrates that the dynamic equilibrium condition holds provided that the emulsions are either kinetically stable or gently stirred after initial mixing because the intermolecular interactions between molecules and ions in emulsions are generally weak and their diffusivities should be orders of magnitude faster than the rates of thermal reaction [3,4,19,24,25,34,[57][58][59][60][61][62][63].

Modeling Chemical Reactivity at the Interfaces of Emulsions
In bulk solution, reactants need to get close to each other to form the so-called "encounter complex" before the chemical reaction takes place, so that the observed (or measured) rate of a bimolecular reaction depends first on the relative rates of these processes. According to the transition state theory (TST), a bimolecular reaction A + B → Products is initiated when the formation of the "encounter pair" is followed by the formation of the activated complex, which must move on to the products for the reaction to proceed, otherwise the reaction does not take place [50,53]. The rate law is given by Equation (5) In multiphasic systems, reactants (and other molecules) are exchanged between the different regions, Scheme 4B, and after bulk mixing is complete, the reactants distribute according to their relative solubilities between the oil, interfacial, and aqueous regions [19,24]. Assuming that the dynamic equilibrium conditions holds, the overall rate of the reaction is given by Equation (6).
where subscripts W, I, and O indicate the oil, interfacial, and aqueous regions, respectively, k W , k I , and k O are the second-order rate constants, parenthesis ( ) denote concentrations in moles per liter of the volume of a particular region, and Φ is the volume fraction of a region, defined as Φ = V region /V Total .
To correctly interpret chemical reactivity in colloidal systems, theoretical frameworks provided by different kinetic models are needed. Quantitative interpretations of experimental results, in terms of reactivity, can only be achieved if the local reagent concentrations and the intrinsic rate constants in the various regions can be obtained from the overall, apparent, experimental kinetic data. For instance, Menger and Portnoy [64] first developed pseudophase models to interpret chemical reactivity in micellar systems. The model was later expanded by Romsted [65][66][67] to incorporate ion-exchange effects. García-Río et al. developed kinetic models to quantitatively explain the influence of microemulsions on chemical reactivity [68] and Sanchez et al. employed the Bronsted relationship to rationalize chemical reactivity in microemulsions [69]. Certainly, it is up to the imagination of researchers to design and employ different chemical probes to assess the distributions of reactants of interest.
However, unraveling reaction mechanisms such as that depicted in Scheme 4B, whose overall rate k obs is given by Equation (6), is extremely complex, and very intricate kinetic equations are frequently obtained because of the huge amount of variables needed to fully describe the partitioning of reactants and the rate of the reaction. Moreover, the rate constants in each region are distinct because the solvent properties of each region are totally different (differences may be of several orders of magnitude). To overcome the problem, judicious simplifications are made by selecting suitable chemical reactions under appropriate conditions. Most frequently, simplifications are achieved through: (1) a careful choice of the chemical structure of the reactants, for example, by employing hydrophobic or hydrophilic reactants, may favor the chemical reaction in one region with respect to the other ones [29,67,68,[70][71][72], (2) by using chemical probes with specific molecular or spectral characteristics; (3) by choosing carefully and appropriate experimental conditions; (4) by using ionic surfactants to promote the concentration of reactants or its electrostatic separation in the interfacial region [28,73,74].
For example, N-nitroso-N-methy-p-toluensulfonamide has a well-established denitrosation mechanism [75], which was exploited to get insights towards the polarity of the Stern layer in micelles [76,77], the role of alcohols in micellar systems [77], and into the solvent properties of microemulsions [27]. An intelligent choice of molecular structures of reactants undergoing nitrosation and denitrosation reactions has also been exploited by Garcia-Río et al. [26,27] and Iglesias et al. [78,79] to investigate the use of microemulsions as reaction media. Aliaga et al. employed a series of 4-alkanoyloxy-1,1,6,6-tetramethylpyperidineoxyl (TEMPO) radicals of different amphiphilic characters as probes to investigate antioxidant activities in micelles and emulsions [80][81][82].
Romsted and coworkers [19,24,62,83] made use of the rich chemistry of arenediazonium ions to investigate the interfacial composition of association colloids and to investigate the distribution of polyphenolic antioxidants in intact emulsions under a variety of experimental conditions [19,24,62]. They synthesized the hydrophobic 4-hexadecylbenzenediazonium, 16-ArN 2 + , tetrafluoroborate, which is reduced in the presence of polyphenolic antioxidants to model the observed reactivity of the probe in emulsions and to estimate the distributions of the antioxidants, which otherwise could not be determined because of the (physical) impossibility of detaching the interfacial region from the aqueous and oil regions [24,[84][85][86][87]. The chemical structures of 16-ArN 2 + and of some representative antioxidants are shown in Scheme 5.
Romsted and coworkers [19,24,62,83] made use of the rich chemistry of arenediazonium ions to investigate the interfacial composition of association colloids and to investigate the distribution of polyphenolic antioxidants in intact emulsions under a variety of experimental conditions [19,24,62]. They synthesized the hydrophobic 4-hexadecylbenzenediazonium, 16-ArN2 + , tetrafluoroborate, which is reduced in the presence of polyphenolic antioxidants to model the observed reactivity of the probe in emulsions and to estimate the distributions of the antioxidants, which otherwise could not be determined because of the (physical) impossibility of detaching the interfacial region from the aqueous and oil regions [24,[84][85][86][87]. The chemical structures of 16-ArN2 + and of some representative antioxidants are shown in Scheme 5.

Scheme 5.
Structures of the chemical probe 16-ArN2 + and of some reactive antioxidants.
The use of the 16-ArN2 + in emulsions considerably simplifies Scheme 4B. The compound 16-ArN2 + is both oil and water insoluble because it is a cation and because of its long hydrophobic tail. Therefore, it locates at the interfacial region of colloidal aggregates, so that their reactive group is anchored in the interfacial region of the emulsion where reacts with other molecules [62,83,88,89]. In this case, it is not necessary to consider the distribution of one of the reactants; the number of unknown parameters decreases significantly, and Scheme 6 now applies. The use of the 16-ArN 2 + in emulsions considerably simplifies Scheme 4B. The compound 16-ArN 2 + is both oil and water insoluble because it is a cation and because of its long hydrophobic tail. Therefore, it locates at the interfacial region of colloidal aggregates, so that their reactive group is anchored in the interfacial region of the emulsion where reacts with other molecules [62,83,88,89]. In this case, it is not necessary to consider the distribution of one of the reactants; the number of unknown parameters decreases significantly, and Scheme 6 now applies. Because the reactive group of 16-ArN2 + is located only in the interfacial region, Scheme 6, and bearing in mind that the oil and water concentrations of 16-ArN2 + are negligible, i.e., [ (7), ΦI is the surfactant volume fraction defined as ΦI = Vsurf/Vtotal. Because the reactive group of 16-ArN 2 + is located only in the interfacial region, Scheme 6, and bearing in mind that the oil and water concentrations of 16 (7), Φ I is the surfactant volume fraction defined as Φ I = V surf /V total .
Scheme 6 can be further simplified by, for instance, carefully choosing the chemical characteristics of the antioxidant so that the antioxidant can be located exclusively in the interfacial region or can only partition between two regions. Some representative examples are given in the next section. For the sake of simplicity, we will first consider non-ionic emulsions, where the droplet surface is not charged and electrostatic effects are negligible, and will consider, in a different section, the case of ionic emulsions.

Main Equations Derived from the Application of the Pseudophase Model Emulsions
Predicting "a priori" the partitioning of a given antioxidant between the oil, interfacial, and aqueous regions of emulsions is truly challenging. Values of the partition constant in binary systems cannot be employed in emulsions because of the physical impossibility of isolating the interfacial region from the oil and aqueous ones, which makes it impossible to predict the values of the partition constants in emulsions (P I O and P I W ) from the value of the partition constant P O W in binary systems. However, knowledge of P O W values may, to a first approximation, help to envisage qualitatively how the distribution of the molecule may be in emulsions. Binary oil-water systems can be considered as limit cases of oil-in-water systems when no surfactant is added, that is, when the concentration of surfactant is low enough to allow phase separation. Two limiting situations can be envisaged, when the antioxidant is hydrophilic (low P O W values), so that their concentration in the oil region is negligible, and the antioxidant only partitions between the aqueous and interfacial regions. Conversely, very hydrophobic antioxidants (high P O W values) are water insoluble and only partition between the oil and interfacial regions. Cautions must be taken, however, because P O W values for food-grade oils cannot be extrapolated from the values in octanolwater binary systems (probably the most common system employed to determine the hydrophobicity of a molecule) because the solvents are totally different [3,19,24,90].

Reactions Where the Antioxidant Partitions between the Three Regions
Antioxidants of moderate hydrophobicity usually distribute between the three regions as shown in Scheme 6. Changes in the polarity of antioxidants can be achieved, for example, by grafting inert alkyl chains in the aromatic ring of polyphenols, so that homologous series of antioxidants bearing the same reactive groups but of different hydrophobicity can be prepared [91,92].
Combination of Equation (7) with the mass balance equations for AO and 16-ArN 2 + , and bearing in mind the definition of the relevant partition constants, Equations (8)- (10) can be derived for k obs [24]. Equation (10) encloses three unknowns, k I , P I O , and P I W , and their individual values cannot be obtained in a single experiment, and a second independent mathematical relationship is needed, which can be obtained by carrying out a second set of kinetic experiments in emulsions with a different Φ W /Φ O ratio to obtain a second set of (k obs , Φ I ) pairs of data. [25,[93][94][95] After some mathematical treatment, Equation (11) can be derived, where the parameters a and b are defined by the Equations (12) and (13) [19,24,25]. The partition constants P I O and P I W , defined by Equations (11) and (12), are then obtained by employing Equation (13) from two sets of kinetic data at different Φ W /Φ O ratios and solving two equations in two unknowns. Once the partition constants are known, the k I value is obtained from Equation (12). Figure 4 illustrates the variations of k obs with Φ I at two Φ W /Φ O ratios.
Molecules 2021, 26, x FOR PEER REVIEW 14 of 26 Figure 4. Variations of kobs with I for methyl caffeate in 4:6 (squares) and 1:9 (circles) olive oil emulsions. The solid lines are the theoretical curves obtained by fitting the experimental data to Equations (16) and (17). Adapted from ref. [96].
A second, much more straightforward approach [25,97], is to measure the partition constant PO W for the antioxidant distribution in binary oil-water mixtures because the value of this partition constant equals to the ratio of the partition constants I W P / I O P , Equation (15). This approach implies to consider the binary system as a constraining condition of the emulsion when the surfactant concentration goes to zero, i.e., when no surfactant is added, Scheme 2.    (16) and (17). Adapted from ref. [96].
A second, much more straightforward approach [25,97], is to measure the partition constant P W O for the antioxidant distribution in binary oil-water mixtures because the value of this partition constant equals to the ratio of the partition constants P I W / P I O , Equation (15). This approach implies to consider the binary system as a constraining condition of the emulsion when the surfactant concentration goes to zero, i.e., when no surfactant is added, Scheme 2.
Values for P I O and P I W are obtained by combining the value for P W O (Equation (14)) with fits of k obs versus Φ I at constant [AO T ] and Φ W :Φ O (Equation (13)) and solving two equations in two unknowns, and once P I O and P I W are known, Equation (12) is used to determine k I . Further details on the calculations, advantages, and limitations of the approaches can be found elsewhere [25,[93][94][95].

Reactions Where the Antioxidant Partitions between Two Regions
We now consider those cases where the antioxidant only partitions between two regions. A number of polar antioxidants such as phenolic acids (gallic, caffeic, catechin) are oil insoluble and mainly distribute between the aqueous and interfacial region, Scheme 7. On the contrary, very hydrophobic antioxidants (such as α-tocopherol (vitamin E) and octyl and lauryl gallates) are water insoluble and only distribute between the oil and interfacial regions, Scheme 7. In both cases, only one partition constant is needed to describe their distribution: that between the water interface, P I W , and that between the oil and interfacial region, P I O , defined by Equations (8) and (9).
Molecules 2021, 26, x FOR PEER REVIEW 15 of 26 oil and interfacial regions, Scheme 7. In both cases, only one partition constant is needed to describe their distribution: that between the water interface, PW I , and that between the oil and interfacial region, PO I , defined by Equations (8) and (9). Bearing in mind Equation (7), the definition of the partition constants, Equations (8) and (9), and the corresponding mass balances, Equations (16)- (19), can be derived. Equations (15) and (16) describe the dependence of kobs on both the concentration of the antioxidant and the medium effects and were used to calculate PW I , kI and PO I , kI values, respectively. Equations (17) and (18) are the reciprocal of Equations (16) and (17), respectively, and predict that plots of 1/kobs vs. ΦI should be linear with a positive intercept. Figure 5 illustrates the variations of kobs vs. ΦI and 1/kobs vs. ΦI for gallic acid and lauryl gallate, and values for kI, PW I and PO I are obtained from the slopes and intercepts of Equations (18) and (19). Scheme 7. Distribution of a hydrophilic (oil insoluble) AO (left) and a hydrophobic (water-insoluble) antioxidant (right) in emulsions. PW I and PO I are the partition constants of the AO between the aqueous and interfacial and the oil and interfacial regions, respectively, and kI is the interfacial rate constant for the reaction between 16-ArN2 + and AO.  Bearing in mind Equation (7), the definition of the partition constants, Equations (8) and (9), and the corresponding mass balances, Equations (16)- (19), can be derived. Equations (15) and (16) describe the dependence of k obs on both the concentration of the antioxidant and the medium effects and were used to calculate P I W , k I and P I O , k I values, respectively. Equations (17) and (18) are the reciprocal of Equations (16) and (17), respectively, and predict that plots of 1/k obs vs. Φ I should be linear with a positive intercept. Figure 5 illustrates the variations of k obs vs. Φ I and 1/k obs vs. Φ I for gallic acid and lauryl gallate, and values for k I , P I W and P I O are obtained from the slopes and intercepts of Equations (18) and (19).  (16) and (18) (A) and (17) and (19)

Reactions Taking Place Exclusively at the Interfacial Region
This represents a limiting case in which both reactants (16-ArN2 + and the antioxidant) do not partition, as illustrated in Scheme 8, and therefore, are found entirely in the interfacial region.    (16) and (18) (A) and (17) and (19)

Reactions Taking Place Exclusively at the Interfacial Region
This represents a limiting case in which both reactants (16-ArN 2 + and the antioxidant) do not partition, as illustrated in Scheme 8, and therefore, are found entirely in the interfacial region.  (16) and (18) (A) and (17) and (19)

Reactions Taking Place Exclusively at the Interfacial Region
This represents a limiting case in which both reactants (16-ArN2 + and the antioxidant) do not partition, as illustrated in Scheme 8, and therefore, are found entirely in the interfacial region.  This case is formally equivalent to that of a bimolecular reaction taking place in a vessel of volume equal to that of the interfacial region, and the observed rate constant for the reaction is given by Equation (20), which predicts that k obs values should decrease upon increasing Φ I , as illustrated in Figure 6.
Scheme 8. Limiting case where both reactants do not partition and are located in the interfacial region of a non-ionic emulsion.  The reciprocal of Equation (20) predicts that plots of 1/k obs vs. Φ I must be straight lines with intercepts at 1/k obs = 0 as illustrated in Figure 6. Some representative values of k I , obtained by least squares fitting of (1/k obs , Φ I ) pairs of data, are shown in Table 3. Table 3. Some representative values of the interfacial rate constant for the reaction between 16-ArN 2 + and octyl and lauryl gallates (OG and LG, respectively) in octane-based emulsions. Note that the intercept is statistically zero as predicted by Equation (9). Data extracted from ref. [99].

AO Emulsion
Intercept Because the values of k I reflect the polarity of the reaction site, the small variations in k I upon changing the oil/water ratio suggest they have negligible effects on the polarity of the interfacial region of both emulsions. The results in Table 3 also suggest that changes in the length of the alkyl chain of the antioxidant have negligible effects on k I values.

Reactions in Ionic Emulsions
When the surfactants employed to stabilize kinetically emulsions are ionic, the surfaces droplets become either negatively or positively charged, Scheme 9, and this charge has significant kinetic consequences because: (1) Electrostatic forces become important because droplet interfaces may attract or repel ions present in the aqueous phase, decreasing their concentration substantially or producing high local concentrations [30,83].
(2) Charged interfaces significantly modify the interfacial molarities of ions such as H + or OH − , having major effects on the pK a of weak acids present in emulsions. Anionic interfaces have higher interfacial H + concentrations relative to those in the bulk, but decrease the interfacial OH − concentration with respect to that in the bulk aqueous phase. Cationic emulsions do the opposite [30,65,100,101].
The pseudophase ion exchange (PIE) model provides a suitable theoretical and operative framework to interpret chemical reactivity at charged emulsion interfaces. The PIE model was first developed for modeling chemical reactivity in micellar solutions of ionic surfactants (Scheme 2) to describe specific ion effects and ion distributions [65,67,100]. The PIE model assumes that the interfacial regions of ionic colloidal systems behave as selective ion exchangers, where counterions (ions of opposite charge to that of the surface) exchange with H + . For example, in anionic sodium dodecyl sulfate micelles, Na + ions are counterions, as well as those from added NaCl. The dynamic equilibria of counterions Figure 7 illustrates the kinetic profiles obtained for the reaction between the antioxidant t-butylhydroquinone and 16-ArN2 + in ionic cetyltrimethylammonium (CTAB) and sodium dodecyl sulfate, SDS, emulsions [85,87]. In CTAB emulsions, Figure 7A, in the absence of salt, the decrease in kobs with increasing ΦI is produced by dilution of TBHQ upon increasing the interfacial volume (i.e., increasing surfactant concentration) [19].  Figure 7 illustrates the kinetic profiles obtained for the reaction between the antioxidant t-butylhydroquinone and 16-ArN 2 + in ionic cetyltrimethylammonium (CTAB) and sodium dodecyl sulfate, SDS, emulsions [85,87]. In CTAB emulsions, Figure 7A, in the absence of salt, the decrease in k obs with increasing Φ I is produced by dilution of TBHQ upon increasing the interfacial volume (i.e., increasing surfactant concentration) [19]. The observed effects on kobs depend on the nature of the reaction considered. For the reaction between ArN2 + ions with antioxidants, ion exchange has important consequences because most antioxidants are (poly)phenolic compounds and the reactive species is the deprotonated form of the phenol, AO − [106]. Thus, at constant [AOT] and ΦO:ΦW, the main effect of increasing ΦI (i.e., addition of CTAB) on kobs is to dilute both the anionic and neutral forms of TBHQ. The effective H + concentration in the interfacial region (mole per liter of interfacial volume) should be reduced to approximately the same amount as the TBHQ and TBHQconcentrations. Therefore, the decrease in kobs with increasing ΦI is caused by dilution of the TBHQ anion within the increasing volume of the totality of the interfacial regions of the emulsions droplets as in nonionic emulsions, Figure 7A [3,19,85].
In SDS emulsions, Figure 7B, the variation of kobs with ΦI decreases initially, passing through a shallow minimum, and then increases. The shallow minimum was attributed to two competing effects, the simultaneous decrease in H + concentration in the interfacial region and the dilution of TBHQ -, which leads to an increase in kobs because of the increase in the concentration of the reactive species TBHQ -. Details can be found elsewhere [19,85,87].

Effects of Temperature on the Kinetics in Emulsions: Interfacial Activation Parameters
The activation parameters for the reaction of 16-ArN2 + with antioxidants in the interfacial region of emulsions can be obtained from the variations of the rate constant kI (see Scheme 6) at different temperatures [107], which are obtained as described in Section 4.1. This is the only method reported in the literature, up to date, that permits estimations of activation energies of reactions taking place in the interfacial region of emulsified systems [59,63,108,109].
Figures 8A illustrates the effect of temperature on kobs for the reaction of 16-ArN2 + with caffeic acid at increasing ΦI values. The straight lines shown in Figure 8B are plots of 1/kobs vs. ΦI, from where the kI values were obtained (Equations (12)-(15)). Similar plots were obtained when employing different antioxidants, [59,63,108,109], and Table 4 shows some values for interfacial rate constants kI for α-tocopherol, caffeic acid and propyl (PG), octyl (OG) and lauryl (LG) gallates (see Scheme 5) at different temperatures. The observed effects on k obs depend on the nature of the reaction considered. For the reaction between ArN 2 + ions with antioxidants, ion exchange has important consequences because most antioxidants are (poly)phenolic compounds and the reactive species is the deprotonated form of the phenol, AO − [106]. Thus, at constant [AO T ] and Φ O :Φ W , the main effect of increasing Φ I (i.e., addition of CTAB) on k obs is to dilute both the anionic and neutral forms of TBHQ. The effective H + concentration in the interfacial region (mole per liter of interfacial volume) should be reduced to approximately the same amount as the TBHQ and TBHQ − concentrations. Therefore, the decrease in k obs with increasing Φ I is caused by dilution of the TBHQ anion within the increasing volume of the totality of the interfacial regions of the emulsions droplets as in nonionic emulsions, Figure 7A [3,19,85].
In SDS emulsions, Figure 7B, the variation of k obs with Φ I decreases initially, passing through a shallow minimum, and then increases. The shallow minimum was attributed to two competing effects, the simultaneous decrease in H + concentration in the interfacial region and the dilution of TBHQ − , which leads to an increase in k obs because of the increase in the concentration of the reactive species TBHQ − . Details can be found elsewhere [19,85,87].

Effects of Temperature on the Kinetics in Emulsions: Interfacial Activation Parameters
The activation parameters for the reaction of 16-ArN 2 + with antioxidants in the interfacial region of emulsions can be obtained from the variations of the rate constant k I (see Scheme 6) at different temperatures [107], which are obtained as described in Section 4.1. This is the only method reported in the literature, up to date, that permits estimations of activation energies of reactions taking place in the interfacial region of emulsified systems [59,63,108,109]. Figure 8A illustrates the effect of temperature on k obs for the reaction of 16-ArN 2 + with caffeic acid at increasing Φ I values. The straight lines shown in Figure 8B are plots of 1/k obs vs. Φ I , from where the k I values were obtained (Equations (12)- (15)). Similar plots were obtained when employing different antioxidants, [59,63,108,109], and Table 4 shows some values for interfacial rate constants k I for α-tocopherol, caffeic acid and propyl (PG), octyl (OG) and lauryl (LG) gallates (see Scheme 5) at different temperatures.     The activation energy, Ea, and the activation parameters, ΔH ‡ and ΔS ‡ , for the reaction in the interfacial region were determined by employing the Arrhenius Equation (25) and that derived from the absolute rate theory, Equation (26), where kB and h refers to the Boltzmann and Plank constants, respectively. Illustrative examples are given in Figure 9.    The activation energy, E a , and the activation parameters, ∆H ‡ and ∆S ‡ , for the reaction in the interfacial region were determined by employing the Arrhenius Equation (25) and that derived from the absolute rate theory, Equation (26), where k B and h refers to the Boltzmann and Plank constants, respectively. Illustrative examples are given in Figure 9. ln

Conclusions and Perspective
Understanding the reactivity and behavior of solute molecules at the interface between immiscible liquids (oil and water) is very important at both the fundamental and theoretical levels because the anisotropic environment of the interfacial region is expected to affect the course of chemical reactions in a way significantly different from that in bulk fluids. Significant differences are also expected from those reactions taking place in a single bulk adjacent to a reactive interface (such as in heterogeneous catalysis). The latter requires the arrival of only one molecule at the interface; meanwhile, reactions taking place in three-dimensional interfaces such as those in emulsions require the formation of an encounter pair, i.e., need the simultaneous arrival at the same location within the interface of two reactive species.
By combining a careful choice of reactants, their properties, and experimental conditions with the theoretical framework provided by pseudophase models, new approaches to model chemical reactivity in three-dimensional interfaces can be obtained. Pseudophase models are powerful tools to model chemical reactivity and to correctly interpret experimental kinetic data, provided that the specificity of the reaction systems is incorporated into the model [71,110,111]. They can be applied to a variety of reactive systems and to reactions carried out under restricted geometry conditions. In fact, the basic equations derived from the application of the pseudophase model are, formally, similar to the mathematical relationships obtained for enzymatic reactions [112,113], and those to model heterogeneous catalysis [50]. However, the use of the pseudophase models requires considering the effects of surface potentials and partitioning effects between regions.
Unfortunately, no complete and systematic theory exists for modeling surface potentials nor solvation effects, and both theoretical and experimental studies of chemical reactions taking place at liquid interfaces are still far from being completely understood, and more experimental data and theoretical models are needed to provide a full understanding of the reaction dynamics at interfaces [114][115][116]. On the theoretical side, progress is needed to understand how the structure of the neat water/organic interface is composed, and how it influences reactant and/or ion transfers and their solvation. The molecular orientation of the reactants at the interface is also relevant because the interfacial region is anisotropic, and certain molecular orientations may be preferred at the interface. A deeper understanding of solvation at liquid interfaces is also important for the correct interpretation of the interfacial phenomena. Solvation plays a central role in un-

Conclusions and Perspective
Understanding the reactivity and behavior of solute molecules at the interface between immiscible liquids (oil and water) is very important at both the fundamental and theoretical levels because the anisotropic environment of the interfacial region is expected to affect the course of chemical reactions in a way significantly different from that in bulk fluids. Significant differences are also expected from those reactions taking place in a single bulk adjacent to a reactive interface (such as in heterogeneous catalysis). The latter requires the arrival of only one molecule at the interface; meanwhile, reactions taking place in threedimensional interfaces such as those in emulsions require the formation of an encounter pair, i.e., need the simultaneous arrival at the same location within the interface of two reactive species.
By combining a careful choice of reactants, their properties, and experimental conditions with the theoretical framework provided by pseudophase models, new approaches to model chemical reactivity in three-dimensional interfaces can be obtained. Pseudophase models are powerful tools to model chemical reactivity and to correctly interpret experimental kinetic data, provided that the specificity of the reaction systems is incorporated into the model [71,110,111]. They can be applied to a variety of reactive systems and to reactions carried out under restricted geometry conditions. In fact, the basic equations derived from the application of the pseudophase model are, formally, similar to the mathematical relationships obtained for enzymatic reactions [112,113], and those to model heterogeneous catalysis [50]. However, the use of the pseudophase models requires considering the effects of surface potentials and partitioning effects between regions.
Unfortunately, no complete and systematic theory exists for modeling surface potentials nor solvation effects, and both theoretical and experimental studies of chemical reactions taking place at liquid interfaces are still far from being completely understood, and more experimental data and theoretical models are needed to provide a full understanding of the reaction dynamics at interfaces [114][115][116]. On the theoretical side, progress is needed to understand how the structure of the neat water/organic interface is composed, and how it influences reactant and/or ion transfers and their solvation. The molecular orientation of the reactants at the interface is also relevant because the interfacial region is anisotropic, and certain molecular orientations may be preferred at the interface. A deeper understanding of solvation at liquid interfaces is also important for the correct interpretation of the interfacial phenomena. Solvation plays a central role in understanding interfacial effects on chemical reactions, because it depends strongly on the solvation free energy of the reactants, products, and of the transition state [114].