Next Article in Journal
Analysis of Vibration Characteristics of Bridge Structures under Seismic Excitation
Next Article in Special Issue
Modelling Heterogeneous Anomalous Dynamics of Radiation-Induced Double-Strand Breaks in DNA during Non-Homologous End-Joining Pathway
Previous Article in Journal
High-Throughput Polar Code Decoders with Information Bottleneck Quantization
Previous Article in Special Issue
Stochastic Compartment Model with Mortality and Its Application to Epidemic Spreading in Complex Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Entropy Production in Reaction–Diffusion Systems Confined in Narrow Channels

by
Guillermo Chacón-Acosta
1,* and
Mayra Núñez-López
2
1
Applied Mathematics and Systems Department, Universidad Autónoma Metropolitana-Cuajimalpa, Vasco de Quiroga 4871, Mexico City 05348, Mexico
2
Department of Mathematics, Instituto Tecnológico Autónomo de México, Río Hondo 1, Col. Progreso Tizapán, Mexico City 01080, Mexico
*
Author to whom correspondence should be addressed.
Entropy 2024, 26(6), 463; https://doi.org/10.3390/e26060463
Submission received: 1 May 2024 / Revised: 26 May 2024 / Accepted: 27 May 2024 / Published: 29 May 2024

Abstract

:
This work analyzes the effect of wall geometry when a reaction–diffusion system is confined to a narrow channel. In particular, we study the entropy production density in the reversible Gray–Scott system. Using an effective diffusion equation that considers modifications by the channel characteristics, we find that the entropy density changes its value but not its qualitative behavior, which helps explore the structure-formation space.

1. Introduction

Many chemical and biological systems exhibit ordered structures that emerge spontaneously and remain stable over time. These structures arise and persist due to processes involving mass transport and chemical reactions under far-from-equilibrium conditions. For instance, Turing’s diffusion-driven instability mechanism describes how small fluctuations in reaction–diffusion systems can be amplified, forming spatio-temporal patterns [1,2]. These patterns arise from the interaction between chemical species, one considered inhibitory and the other activating. Diffusion induces the system’s uniform steady state, which is stable in the absence of diffusion, to become unstable, promoting the emergence of patterns. Transport originates in zones with higher inhibitor concentrations that block the activator, whereas low inhibitor concentrations do not, allowing the activator to grow [2]. This phenomenon is counterintuitive since diffusion is generally associated with a decrease in concentration gradients over time. The reaction–diffusion systems that have been used to elucidate pattern formation typically involve non-linear reactions. They often focus only on forward reactions, assuming their rates significantly surpass those of reverse reactions that have been used to elucidate these phenomena [3]. Consequently, for a complete understanding of these processes’ thermodynamics further study is needed. In classical thermodynamics, the evolution of systems within the linear regime is predictable. However, far-from-equilibrium systems may evolve into organized structures known as dissipative structures, as revealed by some analyses of irreversible thermodynamics [4,5,6].
Numerical studies have revealed some particular connections between entropy production and the emergence of patterns. Entropy production analysis in cycles of chemical oscillations or self-organizing biological systems has become possible in some reversible systems [7]. For instance, in chemical oscillations like the Belousov–Zhabotinsky reaction or the Oregonator model, entropy production averaged over a cycle may exceed or fall below the corresponding unstable equilibrium state, depending on the bifurcation parameter [8]. These processes are relevant, for instance, in self-healing surfaces involving pattern formation [9] or in enzymatic reactions of photosynthetic models that exhibit standard oscillatory and damped dynamics [10]. However, establishing thermodynamic rules for studying the evolution of non-equilibrium patterns remains challenging.
The Gray–Scott model is a powerful tool capable of replicating various Turing-like patterns. This model finds practical applications, accurately describing the autocatalytic reactions observed in chemical and enzymatic reactions in biochemistry and biology [11]. The original Gray–Scott model corresponds to two irreversible reactions for U and V, with P being an inert product [12]. This system models the chemical reaction U + 2 V 3 V , which consumes U and produces V. Therefore, the amount of both substances must be controlled to maintain the reaction. This is done by adding at the feed rate f and removing at the kill rate k. The reduction and subsequent removal of V can also be described by another chemical reaction V P . For this reaction, P is an inert product, meaning it does not react, and the parameter k controls the rate of the second reaction. The Gray–Scott system is defined by two equations that describe the behavior of two reacting substances:
U t = D u 2 U U V 2 + ( 1 U ) f ,
V t = D v 2 V + U V 2 ( f + k ) V ,
where D u and D v are the diffusion coefficients for each chemical species. In order to remove V faster than the growth of U, the rate k is added to f and multiplied by V, since it is assumed that the disappearance of V also depends on its concentration.
The reversible Gray–Scott system has been proposed to study entropy production [3], featuring a cubic autocatalytic reaction coupled with a first-order reaction. Escher, Ross [13], and Yoshida [14] also previously studied it in the context of continuous-flow well-stirred tank reactors (CSTR), where far-from-equilibrium situations can also be achieved. Its phase diagram presents a stationary-stable equilibrium point, saddle-node, and Hopf bifurcations. Heterogeneous structures typically arise near these bifurcation curves [15,16].
Entropy production in self-replicating pulses, initiated by perturbing the equilibrium point, has been investigated [3,17]. While entropy production increases with the size of each pulse, its rate gradually decreases. Thus, the entropy production value can serve as a thermodynamic metric for evaluating pattern stability. Additionally, total entropy change, incorporating outward entropy flux using a relative chemical potential, has been examined within this model [18,19]. Specific trends in entropy production evolution have been observed during the formation of heterogeneous structures. In contrast to uniform steady states, during this process a compromise occurs between maximization and minimization of entropy production, reflecting the influence of several interacting mechanisms [20].
Incorporating a thermal component into the system has revealed transitions between different morphologies of Turing patterns regulated by temperature [21]. These transitions occur with lower entropy production compared to other parameters, such as kinetic constants and reactant fluxes [22]. Formalisms such as the Boltzmann equation [23] and extended irreversible thermodynamics [24] have also drawn further insights into entropy production and entropy flows in reaction–diffusion systems.
On the other hand, the study of transport processes and reactions within confined domains holds significant interest due to its wide range of applications, including carbon nanotubes [25], zeolites [26], biosensors [27], biological membranes [28], separation processes [29,30], and microfluidics, among others [31]. Confinement plays a crucial role in influencing the transport properties by restricting the system’s accessible configuration space [32,33]. The Fick–Jacobs approximation for diffusion in narrow channels simplifies this process by assuming a fast equilibration of the transverse degrees of freedom. It focuses on the longitudinal motion by projecting along the channel axis, resulting in an effective one-dimensional equation with a position-dependent diffusion coefficient [33]. Numerous studies have proposed different functions of this diffusion coefficient under diverse conditions, such as in symmetrical channels [34,35,36,37], asymmetrical channels [38,39,40], channels with curved midlines [41,42,43], or on curved surfaces [44,45,46]. Additionally, confinement effects on chemical reactions [47] and pattern formation in reaction–diffusion systems [48], including diffusive predator–prey models [49] and long-range effects [50], were also explored. Significantly, the geometry of the channel walls introduces modifications to the Turing instability conditions, the dispersion relation, and the range of unstable modes [48].
This work studies how confinement in a narrow channel changes the entropy production density in a reaction–diffusion system. By introducing entropic flow into the Fick–Jacobs approximation, as in [51], it is clear that wall geometry affects entropy production. In the following section, we review the definitions of entropy production for reaction–diffusion systems, focusing on the reversible Gray–Scott model. Subsequently, we explore the effects of geometry on a confined reaction–diffusion system within the local Fick–Jacobs approximation and, particularly, within the mean-field effective framework. We conclude by presenting the effective entropy production calculated with both approaches and the results concerning entropy production density by diffusive flux.

2. Entropy Production in 2D Reaction–Diffusion Systems

Entropy production comprises two distinct components in reaction–diffusion systems, as discussed in [3,5,6,7,8]. The first component arises from chemical reactions’ entropy produced per unit area (or length). The second component is the rate of entropy production due to diffusion. This total entropy production is calculated by integrating the entropy production density across the entire spatial domain:
σ = σ R + σ D .
The contribution of chemical reactions is
σ R = j ( ν j + ν j ) ln ν j + ν j ,
where ν j ± are the forward and backward reaction rates of the j-th reaction. Far from equilibrium, backward rates can be considered small, and the system’s affinity tends to infinity, leading to linear entropy production. In this scenario, chemical oscillations and structure formation can occur [4]. The contribution to entropy production density due to the diffusion of an ideal solution is the product of the diffusive flux and the corresponding driving force:
σ D = k J k · F k = k ( D k P k ) · ln P k = k D k P k P k x 2 + P k y 2 ,
where D k and P k are the diffusivity and concentration of the k-th chemical species; the gas constant R, to be multiplied on the right-hand side, is chosen as 1 for simplicity of calculation [3,7]. Unlike previous studies [18,19], the contribution of entropy flow into and out of the system is not considered.
Irreversible chemical reaction systems are commonly used to model pattern formation, with the reverse reactions often omitted, as mentioned earlier. Accordingly, entropy production diverges for irreversible chemical reactions, as seen in Equation (4). This divergence prevents the calculation and use of entropy production in irreversible reactions [19]. Therefore, reversible chemical reactions, such as the reversible Gray–Scott model, are considered, to overcome this issue. It was first introduced by Mahara et al. in [3,17]. However, Escher, Ross [13], and Yoshida [14] had already considered a similar model for CSTR. The model consists of a cubic autocatalytic reaction U + 2 V 3 V coupled with a first-order reaction V P corresponding to the following system of equations:
U t = D u 2 U k 1 U V 2 + ( 1 U ) f + k 1 V 3 ,
V t = D v 2 V + k 1 U V 2 ( f + k 2 ) V k 1 V 3 + k 2 P ,
P t = D p 2 P + k 2 V ( f + k 2 ) P ,
where U, V, and P are the 2D concentrations of the involved chemical species, k i and k i are the rate constants of the forward and backward reactions, respectively, with i = 1 , 2 , and f is the feeding rate constant. When backward reaction rates vanish, the original Gray–Scott model is recovered. The second reaction rate and the feeding rate describe the system by fixing the first reaction rate, the diffusion coefficients, and the backward rates. This system can produce many different patterns, such as spots and lines, and there are even indications of regions where different patterns can coexist [15]. The stability analysis of the system and the regions of the phase diagram ( k 2 , f ) for which this system can form patterns in one and two dimensions are well-documented [15,17,20]. Studies of the phase diagram of the system indicate that it is monostable outside the saddle-node bifurcation curve but bi-stable between the saddle-node bifurcation curve and the Hopf bifurcation curve. Heterogeneous structures are often found precisely around these bifurcation curves [20].
For simplicity and to minimize the system parameters when solving the equations, the control parameters chosen were k 1 = 1 for the first direct reaction, k 1 = k 2 = 10 3 for the reverse reactions, and the diffusion coefficients D u = 2 × 10 5 , D v = 10 5 and D p = 10 6 . The initial condition was taken as ( U , V , P ) | t = 0 = ( 0.50 , 0.25 , 0.00 ) , as adopted by previous contributions [3,20]. The system size in 2D is 100 × 100 with 2500 time steps, and Neumann boundary conditions were used.
The entropy production density for this 2D system is the sum of these two terms:
σ R = ( k 1 U V 2 k 1 V 3 ) ln k 1 U k 1 V + ( k 2 V k 2 P ) ln k 2 V k 2 P ,
σ D = D u U U x 2 + U y 2 + D v V V x 2 + V y 2 + D p P P x 2 + P y 2 ,
whereas total entropy production is
σ ˜ = Ω d A σ R + σ D .
In reference [17], it was observed that entropy production from reactive terms is approximately 20 times larger than that arising from the diffusion contribution, yet the latter remains important in affecting the dynamics of pattern formation. For this reason, we will only study the diffusive contribution, since the geometry of the boundaries will primarily influence this term. Figure 1 shows a typical two-dimensional pattern of this system and the corresponding diffusive contribution to the entropy density.

3. Pattern Formation in Reaction–Diffusion on Narrow Channels

When chemical species diffuse in a confined space, their transport is driven by stochastic dynamics and the available space for the motion, i.e., entropic effects. In the ideal case, it is possible to model the effect of periodic channel wrinkles on the diffusion of a particle as an overdamped Brownian motion in a one-dimensional periodic potential [52]:
γ x ˙ ( t ) = V ( x ) + F + 2 γ k B T ξ ( t ) ,
where γ is the friction coefficient, V ( x ) is the L-periodic potential, F is the static external force, k B is the Boltzmann constant, T is the temperature, and ξ is a Gaussian noise. In the free case, i.e., V = 0 , the Einstein relation D 0 = k T / γ is recovered. For the case F = 0 , the effective Lifson–Jackson coefficient for external potential is recovered [53]:
D = D 0 L 2 0 L d x 0 L d y e [ V ( x ) V ( x y ) ] / k T .
When considering non-interacting point particles diffusing in 2D, the process can be described by the corresponding Fokker–Planck equation associated to Equation (12), where the Neumann boundary conditions model the confinement.
It is well known that the system is well-approximated by 1D Brownian diffusion in the presence of entropic barriers [33]. In narrow channels characterized by a transverse length that is significantly smaller than the length of the channel axis, the diffusion equation can be projected to one effective dimension. In this case, the marginal concentrations of the k-species are
ρ k ( x , t ) = y 1 ( x ) y 2 ( x ) d y P k ( x , y , t ) ,
where y i ( x ) , i = 1 , 2 , are the functions that determine the geometry of the boundaries and allow us to introduce the so-called width function as w ( x ) = y 2 ( x ) y 1 ( x ) . This projection assumes small amplitude boundary fluctuations w ( x ) 1 and fast equilibration of the transverse degrees of freedom or uniform transverse concentrations. In this case, the 2D equation is replaced by the 1D Fick–Jacobs–Zwanzig operator:
x D k ( x ) w ( x ) x ρ k ( x , t ) w ( x ) ,
where D k ( x ) is the position-dependent diffusion coefficient of the species k, which usually depends on variations of the channel’s width function. According to [48], the reaction kinetic function changes correspondingly, as w F ( ρ k / w ) . Also in [48], it was found that dispersion relations and the range of unstable modes will be affected by the geometry of the system and, therefore, the pattern formation. The major contribution was to explore the impact of channel geometry on reaction–diffusion systems and, consequently, on instability conditions for systems that are constrained by the effective one-dimensional projection technique.
However, the Fick–Jacobs equation is local and it is always desirable to have an equation defined in some particular domain. In [54], the authors present a theoretical link based on the macroscopic diffusion coefficient that is dependent on the porosity, tortuosity, and constriction parameters, by averaging the spatial-dependent diffusion coefficient obtained in the Fick–Jacobs (FJ) approximation. The authors compared the total fluxes predicted by the FJ reduction scheme and the effective medium theory with the numerical results arising from the two-dimensional diffusion equation, in order to show the use of the effective diffusion coefficient, given as follows:
ρ ^ k t = D e f f k W w 2 ρ ^ k x 2
where
D e f f k = D k ϕ δ τ ,
In Equation (16), ρ ^ k is the averaged marginal concentration; D e f f k is the mean-effective diffusion coefficient, in terms of constriction factor δ , tortuosity τ , and porosity ϕ = w / W , where W is the approximate width of the rectangular domain containing the pore. This coefficient is analogous to the Lifson–Jackson averaged coefficient, applicable only to periodic channels [53]. By contrast, the coefficient in Equation (16) can be used for any pore, as long as its length is much greater than its width.

4. Effect of Channel Geometry on Entropy Production

4.1. Fick–Jacobs Approximation for Reaction–Diffusion Equations and Entropy Production

Accordingly, by constraining the dynamics of the reversible Gray–Scott model to a narrow channel domain, the system of equations becomes as follows:
u t = x D u ( x ) w ( x ) x u ( x , t ) w ( x ) k 1 w 2 u v 2 + 1 u w f w + k 1 w 2 v 3 ,
v t = x D v ( x ) w ( x ) x v ( x , t ) w ( x ) + k 1 w 2 u v 2 ( f + k 2 ) v k 1 w 2 v 3 + k 2 p ,
p t = x D p ( x ) w ( x ) x p ( x , t ) w ( x ) + k 2 v ( f + k 2 ) p ,
In this case, u , v , p are the Equation (14) marginal concentrations, noting also that some reaction rates are modified by the channel geometry through the width function.
The chemical contribution to the entropy production density at the channel is
σ r = 1 w 2 k 1 u v 2 k 1 v 3 ln k 1 u k 1 v + ( k 2 v k 2 p ) ln k 2 v k 2 p .
whereas for the effective diffusive component, the Fick–Jacobs probability current and the corresponding driving force, respectively, are given as follows [51]:
J k = D k ( x ) w ( x ) x ρ k ( x , t ) w ( x ) ,
F k = x ln ρ k w ,
and so the contribution to the entropy is σ d = k J k F k ; thereby,
σ d = k D k ( x ) ρ k ( x , t ) w x ρ k w 2 = k D k ( x ) ρ k ( x , t ) ρ k x ρ k w w x 2 ,
  = D u ( x ) u w x u w 2 + D v ( x ) v w x v w 2 + D p ( x ) p w x p w 2 .
Note that if the channel does not vary with the longitudinal coordinate we recover the standard 1D case. Instead, if w 0 , some terms dependent on the variations of the channel geometry arise, i.e., entropic forces modifying the entropy production density emerge.

4.2. Effective Diffusion Coefficient in Terms of Availability, Tortuosity, and Constriction

Using the approach of effective diffusion coefficients based on the channel’s geometric features, the problem is simplified to a one-dimensional system. The averaged concentration is defined as
ρ ^ k = 1 w 0 L d y P k .
From this perspective, the reaction–diffusion system becomes a 1D system with modified diffusion coefficients:
u ^ t = D e f f u W w 2 u ^ x 2 k 1 u ^ v ^ 2 + ( 1 u ^ ) f + k 1 v ^ 3 ,
v ^ t = D e f f v W w 2 v ^ x 2 + k 1 u ^ v ^ 2 ( f + k 2 ) v ^ k 1 v ^ 3 + k 2 p ^ ,
p ^ t = D e f f p W w 2 p ^ x 2 + k 2 v ^ ( f + k 2 ) p ^ ,
where D e f f k is the effective diffusion coefficient for each species, now dependent on geometric factors such as constriction, porosity, and tortuosity, which vary according to the specific shape of each channel.
The density of entropy production from reaction will be the same as in the one-dimensional case:
σ ^ r = w k 1 u ^ v ^ 2 k 1 v ^ 3 ln k 1 u ^ k 1 v ^ + w ( k 2 v ^ k 2 p ^ ) ln k 2 v ^ k 2 p ^ ,
whereas the one corresponding to diffusion flux is as follows:
J k = D e f f k W ρ ^ k x , F k = x ln ρ ^ k ,
such that
σ ^ d = W k D e f f k ρ ^ k ρ ^ k x 2 .
It is worth noting that all these equations apply to the averaged concentration. Therefore, to compare with the one-dimensional marginal concentration, it is necessary to use the approximation ρ k ϕ W ρ ^ k , where ϕ represents availability and W is the effective width of the rectangle where the process occurs.
The reversible Gray–Scott system was integrated using the same parameters as in the first section, with f = 0.028 and k 2 = 0.055 , over a line of 200 points and 5000 steps, with periodic boundary conditions and an initial condition ( 0.5 , 0.25 , 0 ) perturbed by sinusoidal oscillations. Figure 2 shows patterns that appear in the one-dimensional system after some time and persist for a considerable period.
On the other hand, the system (25)–(27) was solved under the same conditions but in the constrained channel shown in Figure 3, with a width function w ( x ) = 1 0.8 sin 7 π x / 200 . According to [54], this channel’s geometric parameters are ( τ = 1 , δ = 0.433 , ϕ = 0.386 ) and the corresponding effective width is W = 2 . It was found that the pattern’s maximum decreases due to geometric effects, as shown in Figure 4.
Finally, the density of entropy production by diffusion at the channel’s center, as a time function, is presented. Figure 5 displays both cases for the function ln ( σ D ) . Again, in the confined case, production is reduced compared to free diffusion. It is noted that there are two moments when this function changes abruptly, corresponding to changes in pattern formation. The monotonous part relates to partial semi-stationary patterns formed in each interval, and the characteristic pattern type is shown in the inset.

5. Discussion

The study of thermodynamics in the context of pattern formation within reaction–diffusion systems offers a fascinating challenge. This is particularly true when studying systems with reversible reactions, where the formation of dissipative structures can be interpreted as the coexistence and competition of two distinct states with high- and low-entropy production. Furthermore, it has been observed that confinement modifies the characteristics of the patterns in irreversible systems, depending on the geometrical parameters of the confining walls. This work examined the reversible Gray–Scott model within a narrow channel. Using the effective Fick–Jacobs model for channel diffusion, we observed that geometric modifications to the pattern are preserved in this reversible system. Additionally, we showed that the entropy production density resulting from diffusive flux is a compelling indicator of the dynamics during structure formation. Notably, while the constrained geometry of the channel influences the value of this function, its qualitative behavior remains consistent.
Our subsequent work may delve into the parameter space to examine the behavior of entropy production across various dissipative structures, similar to what was done in [20]. However, we will also address how the system’s geometry influences this behavior, through local and effective theoretical approaches. It would also be intriguing to explore the impact of geometrical effects on entropy production in non-reversible systems, as in [7,8,10], particularly on systems where the reaction kinetics are influenced by the position and distance between species, as in [49]. A comprehensive approach could involve expanding the parameter space to include the geometrical aspects of confinement, such as slope or curvature within the local Fick–Jacobs approach, as well as tortuosity, availability, or porosity, and constriction in the context of effective theory [46,54]. Such an investigation would require more demanding numerical analyses, which we aim to undertake in forthcoming studies.

Author Contributions

Conceptualization, G.C.-A.; methodology, G.C.-A. and M.N.-L.; formal analysis, G.C.-A.; investigation, G.C.-A. and M.N.-L.; visualization, G.C.-A.; writing—original draft preparation, G.C.-A. and M.N.-L.; writing—review and editing, G.C.-A. and M.N.-L. All authors have read and agreed to the published version of the manuscript.

Funding

M.N.-L. acknowledges support from the Asociación Mexicana de Cultura, A.C.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. The original contributions presented in the study are included in the article.

Acknowledgments

We thank the Special Issue editors for their kind invitation. G.C.-A. acknowledges D.A. León-Velasco for discussions on the matter.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Turing, A.M. The chemical basis of morphogenesis. Philos. Trans. R. Soc. Lond. Ser. B 1952, 237, 37. [Google Scholar]
  2. Murray, J.D.L. Mathematical Biology II: Spatial Models and Biomedical Applications, 3rd ed.; Springer: Berlin/Heidelberg, Germany, 2003. [Google Scholar]
  3. Mahara, H.; Suematsu, N.J.; Yamaguchi, T.; Ohgane, K.; Nishiura, Y.; Shimomura, M. Three-variable reversible Gray–Scott model. J. Chem. Phys. 2004, 121, 8968. [Google Scholar] [CrossRef] [PubMed]
  4. Glansdorff, P.; Prigogine, I. Thermodynamic Theory of Structure, Stability and Fluctuations; Wiley: New York, NY, USA, 1971. [Google Scholar]
  5. Nicolis, G.; Prigogine, I. Self-Organization in Non-Equilibrium System; Wiley: New York, NY, USA, 1977. [Google Scholar]
  6. Kondepudi, D. Introduction to Modern Thermodynamics; John Wiley & Sons, Ltd.: Chichester, UK, 1998. [Google Scholar]
  7. Hanson, M.P. Spatial structures in dissipative systems. J. Chem. Phys. 1974, 60, 3210. [Google Scholar] [CrossRef]
  8. Irvin, B.R.; Ross, J. Calculation of the rate of entropy production for a model chemical reaction. J. Chem. Phys. 1988, 89, 1064. [Google Scholar] [CrossRef]
  9. Kagan, E. Turing Systems, Entropy, and Kinetic Models for Self-Healing Surfaces. Entropy 2010, 12, 554–569. [Google Scholar] [CrossRef]
  10. López-Agudelo, V.A.; Barragán, D. Entropy production in oscillatory processes during photosynthesis. Photochem. Photobiol. Sci. 2014, 13, 23–31. [Google Scholar] [CrossRef] [PubMed]
  11. Gray, P.; Scott, S.K. Chemical Oscillations and Instabilities, Nonlinear Chemical Kinetics; Clarendon: Oxford, UK, 1990. [Google Scholar]
  12. Gray, P.; Scott, S.K. Autocatalytic reactions in the isothermal, continuous stirred tank reactor: Isolas and other forms of multistability. Chem. Eng. Sci. 1983, 38, 29–43. [Google Scholar] [CrossRef]
  13. Escher, C.; Ross, J. Multiranges of flow rate with bistability and limit cycles for Schlögl’s mechanism in a CSTR. J. Chem. Phys. 1983, 79, 3773. [Google Scholar] [CrossRef]
  14. Yoshida, N. Entropy production in a chemical system involving an autocatalytic reaction in an isothermal, continuous stirred tank reactor. J. Chem. Phys. 1990, 92, 2593. [Google Scholar] [CrossRef]
  15. Mahara, H.; Susuki, K.; Jahan, R.A.; Yamaguchi, T. Coexisting stable patterns in a reaction-diffusion system with reversible Gray–Scott dynamics. Phys. Rev. E 2008, 78, 066210. [Google Scholar] [CrossRef] [PubMed]
  16. Kalantarova, J.; Uğrlu, D. Structural stability and stabilization of solutions of the reversible three-component Gray-Scott system. Math. Meth. Appl. Sci. 2019, 42, 3687–3699. [Google Scholar] [CrossRef]
  17. Mahara, H.; Yamaguchi, T.; Shimomura, M. Entropy production in a two-dimensional reversible Gray-Scott system. Chaos 2005, 15, 047508. [Google Scholar] [CrossRef] [PubMed]
  18. Mahara, H.; Yamaguchi, T. Entropy balance in distributed reversible Gray–Scott model. Phys. D 2010, 239, 729–734. [Google Scholar] [CrossRef]
  19. Mahara, H.; Yamaguchi, T. Calculation of the Entropy Balance Equation in a Non-equilibrium Reaction-Diffusion System. Entropy 2010, 12, 2436–2449. [Google Scholar] [CrossRef]
  20. Huang, W.-L.; Li, J. Compromise between minimization and maximization of entropy production in reversible Gray–Scott model. Chem. Eng. Sci. 2016, 155, 233–238. [Google Scholar] [CrossRef]
  21. Serna, H.; Muñuzuri, A.P.; Barragán, D. Thermodynamic and morphological characterization of Turing patterns in non-isothermal reaction–diffusion systems. Phys. Chem. Chem. Phys. 2017, 19, 14401. [Google Scholar] [CrossRef] [PubMed]
  22. Silva-Dias, L.; Lopez-Castillo, A. Turing patterns modulation by chemical gradient in isothermal and non-isothermal conditions. Phys. Chem. Chem. Phys. 2020, 22, 7507. [Google Scholar] [CrossRef]
  23. Mátyás, L.; Gaspard, P. Entropy production in diffusion-reaction systems: The reactive random Lorentz gas. Phys. Rev. E 2005, 71, 036147. [Google Scholar] [CrossRef]
  24. Serdyukov, S. Macroscopic Entropy of Non-Equilibrium Systems and Postulates of Extended Thermodynamics: Application to Transport Phenomena and Chemical Reactions in Nanoparticles. Entropy 2018, 20, 802. [Google Scholar] [CrossRef]
  25. Berezhkovskii, A.; Hummer, G. Single-File Transport of Water Molecules through a Carbon Nanotube. Phys. Rev. Lett. 2002, 89, 064503. [Google Scholar] [CrossRef]
  26. Chester, A.W.; Derouane, E.G. Zeolite Characterization and Catalysis; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  27. Keyser, U.; Koeleman, B.; Dorp, S.V.; Krapf, D.; Smeets, R.; Lemay, S.; Dekker, N.; Dekker, C. Direct force measurements on DNA in a solid-state nanopore. Nat. Phys. 2006, 2, 473–477. [Google Scholar] [CrossRef]
  28. Hille, B. Ion Channels of Excitable Membranes, 3rd ed.; Sinauer Associates, Inc.: Sunderland, MA, USA, 2001. [Google Scholar]
  29. Gouaux, E.; MacKinnon, R. Principles of selective ion transport in channels and pumps. Science 2005, 310, 1461. [Google Scholar] [CrossRef] [PubMed]
  30. Malgaretti, P.; IPagonabarraga, I.; Rubí, J.M. Entropic Electrokinetics: Recirculation, Particle Separation, and Negative Mobility. Phys. Rev. Lett. 2014, 113, 128301. [Google Scholar] [CrossRef]
  31. Hänggi, P.; Marchesoni, F. Artificial Brownian motors: Controlling transport on the nanoscale. Rev. Mod. Phys. 2009, 81, 387–442. [Google Scholar] [CrossRef]
  32. Malgaretti, P.; Pagonabarraga, I.; Rubí, J.M. Confined brownian ratchets. Front. Phys. 2013, 68, 1–21. [Google Scholar] [CrossRef]
  33. Dagdug, L.; Peña, J.; Pompa-García, I. Diffusion under Confinement. A Journey through Counterintuition; Springer Nature: Cham, Switzerland, 2024. [Google Scholar]
  34. Zwanzig, R. Diffusion Past an Entropy Barrier. J. Phys. Chem. 1992, 96, 3926–3930. [Google Scholar] [CrossRef]
  35. Reguera, D.; Rubi, J.M. Kinetic equations for diffusion in the presence of entropic barriers. Phys. Rev. E 2001, 64, 061106. [Google Scholar] [CrossRef] [PubMed]
  36. Kalinay, P.; Percus, J.K. Projection of two-dimensional diffusion in a narrow channel onto the longitudinal dimension. J. Chem. Phys. 2005, 122, 204701. [Google Scholar] [CrossRef] [PubMed]
  37. Kalinay, P.; Percus, J.K. Corrections to the Fick-Jacobs equation. Phys. Rev. E 2006, 74, 041203. [Google Scholar] [CrossRef] [PubMed]
  38. Bradley, R.M. Diffusion in a two-dimensional channel with curved midline and varying width: Reduction to an effective one-dimensional description. Phys. Rev. E 2009, 80, 061142. [Google Scholar] [CrossRef] [PubMed]
  39. Berezhkovskii, A.M.; Szabo, A. Time scale separation leads to position-dependent diffusion along a slow coordinate. J. Chem. Phys. 2011, 135, 074108. [Google Scholar] [CrossRef] [PubMed]
  40. Dagdug, L.; Pineda, I. Projection of two-dimensional diffusion in a curved midline and narrow varying width channel onto the longitudinal dimension. J. Chem. Phys. 2012, 137, 024107. [Google Scholar] [CrossRef] [PubMed]
  41. Ogawa, N. Diffusion in a curved tube. Phys. Lett. A 2013, 377, 2465. [Google Scholar] [CrossRef]
  42. Dagdug, L.; García-Chung, A.A.; Chacón-Acosta, G. On the description of Brownian particles in confinement on a non-Cartesian coordinates basis. J. Chem. Phys. 2016, 145, 074105. [Google Scholar] [CrossRef] [PubMed]
  43. Chávez, Y.; Chacón-Acosta, G.; Dagdug, L. Unbiased diffusion of Brownian particles in a helical tube. J. Chem. Phys. 2018, 148, 214106. [Google Scholar] [CrossRef] [PubMed]
  44. Chacón-Acosta, G.; Pineda, I.; Dagdug, L. Diffusion in narrow channels on curved manifolds. J. Chem. Phys. 2013, 139. [Google Scholar] [CrossRef] [PubMed]
  45. Pineda, I.; Chacón-Acosta, G.; Dagdug, L. Diffusion coefficients for two-dimensional narrow asymmetric channels embedded on flat and curved surfaces. Eur. Phys. J. Spec. Top. 2014, 223, 3045. [Google Scholar] [CrossRef]
  46. Ledesma-Durán, A.; León-Velasco, D.A.; Chacón-Acosta, G.; Juárez-Valencia, L.H. Surface diffusion in narrow channels on curved domains. Phys. Rev. E 2023, 107, 034801. [Google Scholar] [CrossRef]
  47. Ledesma-Durán, A.; Hernéz-Hernéz, S.I.; Santamaría-Holek, I. Generalized Fick-Jacobs Approach for Describing Adsorption-Desorption Kinetics in Irregular Pores under Nonequilibrium Conditions. J. Phys. Chem. C 2016, 120, 7810. [Google Scholar] [CrossRef]
  48. Chacón-Acosta, G.; Núñez-López, M.; Pineda, I. Turing instability conditions in confined systems with an effective position-dependent diffusion-coefficient. J. Chem. Phys. 2020, 152, 024101. [Google Scholar] [CrossRef]
  49. Núñez-López, M.; Chacón-Acosta, G. Pattern formation in a predator-prey system with finite interaction range in a channel-like region using the Fick-Jacobs diffusion approach. Phys. D 2022, 433, 133194. [Google Scholar] [CrossRef]
  50. Chacón-Acosta, G.; León-Ramírez, A.; González-Gaxiola, O. Biharmonic Fick–Jacobs diffusion in narrow channels. Phys. A 2023, 628, 129155. [Google Scholar] [CrossRef]
  51. Carusela, M.F.; Rubi, J.M. Entropy production and rectifcation efficiency in colloid transport along a pulsating channel. J. Phys. Condens. Matter 2018, 30, 244001. [Google Scholar] [CrossRef] [PubMed]
  52. Burada, P.S.; Hänggi, P.; Marchesoni, F.; Schmid, G.; Talkner, P. Diffusion in Confined Geometries. ChemPhysChem 2009, 10, 45–54. [Google Scholar] [CrossRef] [PubMed]
  53. Lifson, S.; Jackson, J.L. On the Self-Diffusion of Ions in a Polyelectrolyte Solution. J. Chem. Phys. 1962, 36, 2410. [Google Scholar] [CrossRef]
  54. Ledesma-Durán, A.; Hernández, S.I.; Santamaría-Holek, I. Relation between the porosity and tortuosity of a membrane formed by disconnected irregular pores and the spatial diffusion coefficient of the Fick-Jacobs model. Phys. Rev. E 2017, 95, 052804. [Google Scholar] [CrossRef]
Figure 1. Patterns for the U ( x , y , t ) concentration and diffusive entropy production σ D ( x , y , t ) for the reversible Gray–Scott model Equations (6)–(8), both evaluated at t = 2000 . In (a), the color gradient ranges from bright yellow tones representing a higher concentration of U to dark blue tones indicating a lower concentration. In (b), bold tones reflect lower values of entropy density, while lighter tones indicate higher values.
Figure 1. Patterns for the U ( x , y , t ) concentration and diffusive entropy production σ D ( x , y , t ) for the reversible Gray–Scott model Equations (6)–(8), both evaluated at t = 2000 . In (a), the color gradient ranges from bright yellow tones representing a higher concentration of U to dark blue tones indicating a lower concentration. In (b), bold tones reflect lower values of entropy density, while lighter tones indicate higher values.
Entropy 26 00463 g001
Figure 2. Example of spatio-temporal U concentration patterns in 1D reversible Gray–Scott model.
Figure 2. Example of spatio-temporal U concentration patterns in 1D reversible Gray–Scott model.
Entropy 26 00463 g002
Figure 3. Diagram of the geometry of the constrained channel with width function w ( x ) = 1 0.8 sin 7 π x / 200 , with aspect ratio 4:1, in the cell L × W .
Figure 3. Diagram of the geometry of the constrained channel with width function w ( x ) = 1 0.8 sin 7 π x / 200 , with aspect ratio 4:1, in the cell L × W .
Entropy 26 00463 g003
Figure 4. Comparison of the spatio-temporal patterns. Patterns for the u ( x , y , t ) concentration for both free 1D diffusion and diffusion on a constrained channel for the reversible Gray–Scott model at t = 4000 .
Figure 4. Comparison of the spatio-temporal patterns. Patterns for the u ( x , y , t ) concentration for both free 1D diffusion and diffusion on a constrained channel for the reversible Gray–Scott model at t = 4000 .
Entropy 26 00463 g004
Figure 5. Diffusive entropy production density for free 1D diffusion, and diffusion on constrained channel at the channel’s center, as a function of time.
Figure 5. Diffusive entropy production density for free 1D diffusion, and diffusion on constrained channel at the channel’s center, as a function of time.
Entropy 26 00463 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chacón-Acosta, G.; Núñez-López, M. Entropy Production in Reaction–Diffusion Systems Confined in Narrow Channels. Entropy 2024, 26, 463. https://doi.org/10.3390/e26060463

AMA Style

Chacón-Acosta G, Núñez-López M. Entropy Production in Reaction–Diffusion Systems Confined in Narrow Channels. Entropy. 2024; 26(6):463. https://doi.org/10.3390/e26060463

Chicago/Turabian Style

Chacón-Acosta, Guillermo, and Mayra Núñez-López. 2024. "Entropy Production in Reaction–Diffusion Systems Confined in Narrow Channels" Entropy 26, no. 6: 463. https://doi.org/10.3390/e26060463

APA Style

Chacón-Acosta, G., & Núñez-López, M. (2024). Entropy Production in Reaction–Diffusion Systems Confined in Narrow Channels. Entropy, 26(6), 463. https://doi.org/10.3390/e26060463

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

Article Metrics

Back to TopTop