Pub Date : 2024-07-17DOI: 10.1007/s10910-024-01660-2
M. Sivashankar, S. Sabarinathan, Hasib Khan, Jehad Alzabut, J. F. Gómez-Aguilar
Kinetic chemical reactions find applications across various fields. In industrial processes, they drive the production of essential materials like fertilizers and pharmaceuticals. In environmental science, they are crucial to understanding pollution dynamics. Additionally, in biochemistry, they underpin vital cellular processes, offering insights into disease mechanisms and drug development. In this work, we present a new advancement of a dynamical system for kinetically controlled chemical reactions and the dependency of its solution on the initial conditions using mathematical techniques for fractional orders. By utilizing this fixed-point approach, we can derive the existence and uniqueness theorem of the proposed model. We further show that the chemical kinetics of the fractional model are stable through the Hyers-Ulam stability condition. We finally run a numerical simulation to verify our conclusions. The manuscript concludes with demonstrative examples.
{"title":"Stability and computational results for chemical kinetics reactions in enzyme","authors":"M. Sivashankar, S. Sabarinathan, Hasib Khan, Jehad Alzabut, J. F. Gómez-Aguilar","doi":"10.1007/s10910-024-01660-2","DOIUrl":"10.1007/s10910-024-01660-2","url":null,"abstract":"<div><p>Kinetic chemical reactions find applications across various fields. In industrial processes, they drive the production of essential materials like fertilizers and pharmaceuticals. In environmental science, they are crucial to understanding pollution dynamics. Additionally, in biochemistry, they underpin vital cellular processes, offering insights into disease mechanisms and drug development. In this work, we present a new advancement of a dynamical system for kinetically controlled chemical reactions and the dependency of its solution on the initial conditions using mathematical techniques for fractional orders. By utilizing this fixed-point approach, we can derive the existence and uniqueness theorem of the proposed model. We further show that the chemical kinetics of the fractional model are stable through the Hyers-Ulam stability condition. We finally run a numerical simulation to verify our conclusions. The manuscript concludes with demonstrative examples.</p></div>","PeriodicalId":648,"journal":{"name":"Journal of Mathematical Chemistry","volume":"62 9","pages":"2346 - 2367"},"PeriodicalIF":1.7,"publicationDate":"2024-07-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141740117","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-07-16DOI: 10.1007/s10910-024-01655-z
Michael Perfileev, Vladislav Lyakishev
This work is based on the productive idea of Mulliken about the alignment of electronegativities of atoms in the process of bond formation to their geometric mean value. The paper considers in detail the case of a binary molecule and obtains formulas for the dependence of the current values of the electronegativities of the two atoms forming the molecule on time, and finds a mathematical connection between the current and initial values of electronegativities. Also, in the work the theorem on the relation between the rates of alignment of electronegativities of atoms entering into chemical bonding is formulated and proved, and a special case of this theorem is considered.
{"title":"Theorem on rates of alignment of electronegativities of atoms in the process of formation of a chemical bond in a binary molecule","authors":"Michael Perfileev, Vladislav Lyakishev","doi":"10.1007/s10910-024-01655-z","DOIUrl":"10.1007/s10910-024-01655-z","url":null,"abstract":"<div><p>This work is based on the productive idea of Mulliken about the alignment of electronegativities of atoms in the process of bond formation to their geometric mean value. The paper considers in detail the case of a binary molecule and obtains formulas for the dependence of the current values of the electronegativities of the two atoms forming the molecule on time, and finds a mathematical connection between the current and initial values of electronegativities. Also, in the work the theorem on the relation between the rates of alignment of electronegativities of atoms entering into chemical bonding is formulated and proved, and a special case of this theorem is considered.</p></div>","PeriodicalId":648,"journal":{"name":"Journal of Mathematical Chemistry","volume":"62 9","pages":"2326 - 2333"},"PeriodicalIF":1.7,"publicationDate":"2024-07-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141644020","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-07-15DOI: 10.1007/s10910-024-01651-3
Sarika Kumari, Lok Nath Kannaujiya, Narendra Kumar, Amit K. Verma, Ravi P. Agarwal
In response to the challenges posed by complex boundary conditions and singularities in molecular systems and quantum chemistry, accurately determining energy levels (eigenvalues) and corresponding wavefunctions (eigenfunctions) is crucial for understanding molecular behavior and interactions. Mathematically, eigenvalues and normalized eigenfunctions play crucial role in proving the existence and uniqueness of solutions for nonlinear boundary value problems (BVPs). In this paper, we present an iterative procedure for computing the eigenvalues ((mu )) and normalized eigenfunctions of novel fractional singular eigenvalue problems,
$$begin{aligned} D^{2alpha } y(t) + frac{k}{t^alpha } D^alpha y(t) + mu y (t) =0,~~ 0< t<1,~~ 0< alpha le 1, end{aligned}$$
with boundary condition,
$$y'(0)=0, ~~~~y(1)=0,$$
where (D^alpha , D^{2alpha }) represents the Caputo fractional derivative, (k ge 1). We propose a novel method for computing Lagrange multipliers, which enhances the variational iteration method to yield convergent solutions. Numerical findings suggest that this strategy is simple yet powerful and effective.
为了应对分子系统和量子化学中复杂边界条件和奇异性带来的挑战,准确确定能级(特征值)和相应的波函数(特征函数)对于理解分子行为和相互作用至关重要。在数学上,特征值和归一化特征函数在证明非线性边界值问题(BVP)解的存在性和唯一性方面起着至关重要的作用。在本文中,我们提出了一种迭代过程,用于计算新型分数奇异特征值问题的特征值((mu ))和归一化特征函数,$$begin{aligned}。D^{2alpha } y(t) + frac{k}{t^alpha }D^{2alpha } y(t) + mu y (t) =0,~~ 0< t<1,~~ 0< alpha le 1, end{aligned}$$边界条件为$y'(0)=0, ~~~~y(1)=0,$$ 其中 (D^alpha , D^{2alpha }) 表示 Caputo 分数导数, (k ge 1).我们提出了一种计算拉格朗日乘数的新方法,该方法增强了变分迭代法,从而产生收敛解。数值结果表明,这一策略简单而强大有效。
{"title":"A novel hybrid variation iteration method and eigenvalues of fractional order singular eigenvalue problems","authors":"Sarika Kumari, Lok Nath Kannaujiya, Narendra Kumar, Amit K. Verma, Ravi P. Agarwal","doi":"10.1007/s10910-024-01651-3","DOIUrl":"10.1007/s10910-024-01651-3","url":null,"abstract":"<div><p>In response to the challenges posed by complex boundary conditions and singularities in molecular systems and quantum chemistry, accurately determining energy levels (eigenvalues) and corresponding wavefunctions (eigenfunctions) is crucial for understanding molecular behavior and interactions. Mathematically, eigenvalues and normalized eigenfunctions play crucial role in proving the existence and uniqueness of solutions for nonlinear boundary value problems (BVPs). In this paper, we present an iterative procedure for computing the eigenvalues (<span>(mu )</span>) and normalized eigenfunctions of novel fractional singular eigenvalue problems, </p><div><div><span>$$begin{aligned} D^{2alpha } y(t) + frac{k}{t^alpha } D^alpha y(t) + mu y (t) =0,~~ 0< t<1,~~ 0< alpha le 1, end{aligned}$$</span></div></div><p>with boundary condition, </p><div><div><span>$$y'(0)=0, ~~~~y(1)=0,$$</span></div></div><p>where <span>(D^alpha , D^{2alpha })</span> represents the Caputo fractional derivative, <span>(k ge 1)</span>. We propose a novel method for computing Lagrange multipliers, which enhances the variational iteration method to yield convergent solutions. Numerical findings suggest that this strategy is simple yet powerful and effective.</p></div>","PeriodicalId":648,"journal":{"name":"Journal of Mathematical Chemistry","volume":"62 9","pages":"2304 - 2325"},"PeriodicalIF":1.7,"publicationDate":"2024-07-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141646962","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-07-14DOI: 10.1007/s10910-024-01648-y
Nikolay A. Tikhonov
The phenomenon of isothermal supersaturation of solutions in a porous medium at ion exchange is studied on the basis of mathematical modeling. The phenomenon consists in the fact that the solution with concentration significantly higher than the maximal solubility of the substance is formed in the pores of sorbent and no precipitation occurs. The question of why sediment does not appear in the pores between the grains is investigated in the article. It is shown that the phenomenon under consideration can be explained by the effect of dynamic equilibrium between the association of condensed phase particles in the inner part of the pores, their diffusion, and decomposition near the surface of the sorbent grains caused by a change in potential. The degree of possible supersaturation of the solution is estimated depending on the process parameters. The proposed hypothesis is confirmed by quantitative studies using the available experimental data.
{"title":"Modeling of isothermal supersaturation of solutions in a porous medium: estimation of it’s possible degree","authors":"Nikolay A. Tikhonov","doi":"10.1007/s10910-024-01648-y","DOIUrl":"10.1007/s10910-024-01648-y","url":null,"abstract":"<div><p>The phenomenon of isothermal supersaturation of solutions in a porous medium at ion exchange is studied on the basis of mathematical modeling. The phenomenon consists in the fact that the solution with concentration significantly higher than the maximal solubility of the substance is formed in the pores of sorbent and no precipitation occurs. The question of why sediment does not appear in the pores between the grains is investigated in the article. It is shown that the phenomenon under consideration can be explained by the effect of dynamic equilibrium between the association of condensed phase particles in the inner part of the pores, their diffusion, and decomposition near the surface of the sorbent grains caused by a change in potential. The degree of possible supersaturation of the solution is estimated depending on the process parameters. The proposed hypothesis is confirmed by quantitative studies using the available experimental data.</p></div>","PeriodicalId":648,"journal":{"name":"Journal of Mathematical Chemistry","volume":"62 9","pages":"2289 - 2303"},"PeriodicalIF":1.7,"publicationDate":"2024-07-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141649816","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-07-12DOI: 10.1007/s10910-024-01650-4
Jun Fujisawa
A fullerene graph is a 3-connected plane cubic graph in which every face is pentagonal or hexagonal. A set of hexagons (mathcal {H}) of G is called a resonant pattern if there exists a perfect matching M of G such that exactly three edges of H is contained in M for each member H of (mathcal {H}). In this paper we prove for any natural number k that almost all of the family of k disjoint hexagons are resonant patterns in sufficiently large fullerene graphs.
富勒烯图是一个三连平面立方图,其中每个面都是五边形或六边形。如果 G 中存在完美匹配的 M,使得 H 的每一个成员 H 都有三条边包含在 M 中,那么 G 的六边形集合 (mathcal {H}) 就被称为共振图案。在本文中,我们证明了对于任意自然数 k,在足够大的富勒烯图中,几乎所有 k 个不相交的六边形族都是共振图案。
{"title":"Resonant hexagons in fullerene graphs","authors":"Jun Fujisawa","doi":"10.1007/s10910-024-01650-4","DOIUrl":"10.1007/s10910-024-01650-4","url":null,"abstract":"<div><p>A fullerene graph is a 3-connected plane cubic graph in which every face is pentagonal or hexagonal. A set of hexagons <span>(mathcal {H})</span> of <i>G</i> is called a <i>resonant pattern</i> if there exists a perfect matching <i>M</i> of <i>G</i> such that exactly three edges of <i>H</i> is contained in <i>M</i> for each member <i>H</i> of <span>(mathcal {H})</span>. In this paper we prove for any natural number <i>k</i> that almost all of the family of <i>k</i> disjoint hexagons are resonant patterns in sufficiently large fullerene graphs.</p></div>","PeriodicalId":648,"journal":{"name":"Journal of Mathematical Chemistry","volume":"62 9","pages":"2280 - 2288"},"PeriodicalIF":1.7,"publicationDate":"2024-07-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141609647","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-07-11DOI: 10.1007/s10910-024-01654-0
Jasinth Sylvia, Surath Ghosh
Throughout this research paper, we represent a highly effective Haar wavelet technique to determine the solution of the complex nonlinear dynamical system with three variables chemical reaction model. The foremost objective of this study is to represent the dynamical behavior of chemical reaction model in the sense of Caputo derivative. The convergent analysis and stability analysis of the three variable chemical reaction model are discussed. The existence and uniqueness of the given model is also verified. Furthermore, the residual error analysis for this model is also presented. In addition, graphically the numerical solutions in a 2-dimensional and 3-dimensional manner are obtained by using MATLAB (2023a).
{"title":"Solution of chemical reaction model using Haar wavelet method with Caputo derivative","authors":"Jasinth Sylvia, Surath Ghosh","doi":"10.1007/s10910-024-01654-0","DOIUrl":"10.1007/s10910-024-01654-0","url":null,"abstract":"<div><p>Throughout this research paper, we represent a highly effective Haar wavelet technique to determine the solution of the complex nonlinear dynamical system with three variables chemical reaction model. The foremost objective of this study is to represent the dynamical behavior of chemical reaction model in the sense of Caputo derivative. The convergent analysis and stability analysis of the three variable chemical reaction model are discussed. The existence and uniqueness of the given model is also verified. Furthermore, the residual error analysis for this model is also presented. In addition, graphically the numerical solutions in a 2-dimensional and 3-dimensional manner are obtained by using MATLAB (2023a).</p></div>","PeriodicalId":648,"journal":{"name":"Journal of Mathematical Chemistry","volume":"62 9","pages":"2222 - 2252"},"PeriodicalIF":1.7,"publicationDate":"2024-07-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141588596","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The Gruneisen parameter offers crucial insights into the frequency distribution of the phonon spectrum in solids. In present study, we focus on the theoretical prediction of Gruneisen parameter for magnesium chalcopyrites MgSiP2, MgSiAs2, and MgSiSb2 by using three different logarithmic equation of state (EOS) viz. Poirier Tarantola EOS, Third-Order EOS, and Bardeen EOS at varying compression values (V/V0). These EOSs are subjected to rigorous testing against the fundamental thermodynamic requirements, especially at extreme compression limits. It is observed that at low compressions, all three EOSs—Poirier Tarantola, Third-Order EOS and Bardeen EOS yield identical results. However, when estimating the Gruneisen parameter at high compression, we found that after compression range V/V0 = 0.98 for MgSiAs2 the Poirier Tarantola EOS gets deviated with other two EOSs and also after compression range V/V0 = 0.99 for MgSiP2 the Poirier Tarantola EOS gets deviated with other two EOSs and after compression range V/V0 = 0.99 for MgSiSb2 the third order EOS get deviated with other two EOS.
格鲁尼森参数是了解固体声子频谱频率分布的重要依据。在本研究中,我们通过使用三种不同的对数状态方程(EOS),即不同压缩值(V/V0)下的 Poirier Tarantola EOS、Third-Order EOS 和 Bardeen EOS,重点对黄铜矿镁硅石 MgSiP2、MgSiAs2 和 MgSiSb2 的 Gruneisen 参数进行了理论预测。根据基本热力学要求对这些 EOS 进行了严格测试,特别是在极端压缩极限下。结果表明,在低压缩条件下,所有三种 EOS--波里埃-塔兰托拉 EOS、三阶 EOS 和巴丁 EOS 得出的结果完全相同。然而,在估算高压缩时的格鲁尼森参数时,我们发现在硅铝镁的压缩范围 V/V0 = 0.98 之后,普瓦里尔-塔兰托拉方程与其他两种方程产生了偏差;在硅铝镁的压缩范围 V/V0 = 0.99 之后,普瓦里尔-塔兰托拉方程与其他两种方程产生了偏差;在硅铝镁的压缩范围 V/V0 = 0.99 之后,三阶方程与其他两种方程产生了偏差。
{"title":"Theoretical prediction of Gruneisen parameter for chalcopyrites","authors":"Shipra Tripathi, Abhi Sarika Bharti, Anjani Kumar Pandey, Chandra Kumar Dixit","doi":"10.1007/s10910-024-01645-1","DOIUrl":"10.1007/s10910-024-01645-1","url":null,"abstract":"<div><p>The Gruneisen parameter offers crucial insights into the frequency distribution of the phonon spectrum in solids. In present study, we focus on the theoretical prediction of Gruneisen parameter for magnesium chalcopyrites MgSiP<sub>2</sub>, MgSiAs<sub>2</sub>, and MgSiSb<sub>2</sub> by using three different logarithmic equation of state (EOS) viz. Poirier Tarantola EOS, Third-Order EOS, and Bardeen EOS at varying compression values (V/V<sub>0</sub>). These EOSs are subjected to rigorous testing against the fundamental thermodynamic requirements, especially at extreme compression limits. It is observed that at low compressions, all three EOSs—Poirier Tarantola, Third-Order EOS and Bardeen EOS yield identical results. However, when estimating the Gruneisen parameter at high compression, we found that after compression range V/V<sub>0</sub> = 0.98 for MgSiAs<sub>2</sub> the Poirier Tarantola EOS gets deviated with other two EOSs and also after compression range V/V<sub>0</sub> = 0.99 for MgSiP<sub>2</sub> the Poirier Tarantola EOS gets deviated with other two EOSs and after compression range V/V<sub>0</sub> = 0.99 for MgSiSb<sub>2</sub> the third order EOS get deviated with other two EOS.</p></div>","PeriodicalId":648,"journal":{"name":"Journal of Mathematical Chemistry","volume":"62 9","pages":"2265 - 2279"},"PeriodicalIF":1.7,"publicationDate":"2024-07-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141609644","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-07-11DOI: 10.1007/s10910-024-01647-z
Abhay P. Srivastava, Brijesh K. Pandey, A. K. Gupta, Anjani K. Pandey
The study utilized the theory of interionic potentials and included analytical functions to account for the volume-dependent short-range force constant. Specifically, a modified version of the Shanker equation of state was used, and expressions were established for isothermal bulk modulus and its pressure derivatives. The researcher extensively analyzed the bismuth ferrite (BiFeO3) material at pressures up to 10 GPa. The result obtained by the newly derived equation of state is compared against previously obtained equations of state, including the Shanker and Vinet equation of state and experimental data. Graphical representations demonstrate the changes in pressure, bulk modulus, and pressure derivative of bulk modulus with compression. The result shows that the newly developed equation of state provides superior outcomes compared to the Shanker and Vinet equations, particularly at high compression levels, due to the inclusion of higher-order compression terms.
{"title":"Theoretical prediction of thermoelastic properties of bismuth ferrite by a new approach","authors":"Abhay P. Srivastava, Brijesh K. Pandey, A. K. Gupta, Anjani K. Pandey","doi":"10.1007/s10910-024-01647-z","DOIUrl":"10.1007/s10910-024-01647-z","url":null,"abstract":"<div><p>The study utilized the theory of interionic potentials and included analytical functions to account for the volume-dependent short-range force constant. Specifically, a modified version of the Shanker equation of state was used, and expressions were established for isothermal bulk modulus and its pressure derivatives. The researcher extensively analyzed the bismuth ferrite (BiFeO<sub>3</sub><b>)</b> material at pressures up to 10 GPa. The result obtained by the newly derived equation of state is compared against previously obtained equations of state, including the Shanker and Vinet equation of state and experimental data. Graphical representations demonstrate the changes in pressure, bulk modulus, and pressure derivative of bulk modulus with compression. The result shows that the newly developed equation of state provides superior outcomes compared to the Shanker and Vinet equations, particularly at high compression levels, due to the inclusion of higher-order compression terms.</p></div>","PeriodicalId":648,"journal":{"name":"Journal of Mathematical Chemistry","volume":"62 9","pages":"2253 - 2264"},"PeriodicalIF":1.7,"publicationDate":"2024-07-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141587098","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-07-08DOI: 10.1007/s10910-024-01646-0
Xianfa Hu, Wansheng Wang, Bin Wang, Yonglei Fang
In this paper, two novel classes of implicit exponential Runge–Kutta (ERK) methods are studied for solving highly oscillatory systems. First of all, symplectic conditions for two kinds of exponential integrators are derived, and we present a first-order symplectic method. High accurate implicit ERK methods (up to order four) are formulated by comparing the Taylor expansion of the exact solution, it is shown that the order conditions of two new kinds of exponential methods are identical to the order conditions of classical Runge–Kutta (RK) methods. Moreover, we investigate the linear stability properties of these exponential methods. Numerical examples not only present the long time energy preservation of the first-order symplectic method, but also illustrate the accuracy and efficiency of these formulated methods in comparison with standard ERK methods.
{"title":"Cost-reduction implicit exponential Runge–Kutta methods for highly oscillatory systems","authors":"Xianfa Hu, Wansheng Wang, Bin Wang, Yonglei Fang","doi":"10.1007/s10910-024-01646-0","DOIUrl":"10.1007/s10910-024-01646-0","url":null,"abstract":"<div><p>In this paper, two novel classes of implicit exponential Runge–Kutta (ERK) methods are studied for solving highly oscillatory systems. First of all, symplectic conditions for two kinds of exponential integrators are derived, and we present a first-order symplectic method. High accurate implicit ERK methods (up to order four) are formulated by comparing the Taylor expansion of the exact solution, it is shown that the order conditions of two new kinds of exponential methods are identical to the order conditions of classical Runge–Kutta (RK) methods. Moreover, we investigate the linear stability properties of these exponential methods. Numerical examples not only present the long time energy preservation of the first-order symplectic method, but also illustrate the accuracy and efficiency of these formulated methods in comparison with standard ERK methods.</p></div>","PeriodicalId":648,"journal":{"name":"Journal of Mathematical Chemistry","volume":"62 9","pages":"2191 - 2221"},"PeriodicalIF":1.7,"publicationDate":"2024-07-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141567830","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-07-04DOI: 10.1007/s10910-024-01630-8
Homero G. Díaz-Marín, José L. Sánchez-Ponce
Enzyme catalysis in reactors for industrial applications usually require an external intervention of the species involved in the chemical reactions. We analyze the most elementary open enzyme catalysis with competitive inhibition where a time-dependent inflow of substrate and inhibitor supplies is modeled by almost periodic functions. We prove global stability of an almost periodic solution for the non-autonomous dynamical system arising from the mass-law action. This predicts a well behaved situation in which the reactor oscillates with global stability. This is a first case study in the path toward broader global stability results regarding intraspecific and monotone open reaction networks.
{"title":"Intraspecific and monotone enzyme catalysis with oscillatory substrate and inhibitor supplies","authors":"Homero G. Díaz-Marín, José L. Sánchez-Ponce","doi":"10.1007/s10910-024-01630-8","DOIUrl":"10.1007/s10910-024-01630-8","url":null,"abstract":"<div><p>Enzyme catalysis in reactors for industrial applications usually require an external intervention of the species involved in the chemical reactions. We analyze the most elementary open enzyme catalysis with competitive inhibition where a time-dependent inflow of substrate and inhibitor supplies is modeled by almost periodic functions. We prove global stability of an almost periodic solution for the non-autonomous dynamical system arising from the mass-law action. This predicts a well behaved situation in which the reactor oscillates with global stability. This is a first case study in the path toward broader global stability results regarding intraspecific and monotone open reaction networks.</p></div>","PeriodicalId":648,"journal":{"name":"Journal of Mathematical Chemistry","volume":"62 9","pages":"2160 - 2190"},"PeriodicalIF":1.7,"publicationDate":"2024-07-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141550393","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}