This study examines a fully parabolic predator-prey chemo-alarm-taxis system under homogeneous Neumann boundary conditions in a bounded domain $Ω\subset \mathbb{R}^n$ with a smooth boundary $\partialΩ$. Under specific parameter conditions, it is shown that the system admits a unique, globally bounded classical solution. The convergence of the solution is established through the construction of an appropriate Lyapunov functional. In addition, numerical simulations are presented to validate the asymptotic behaviour of the solution. The results highlight the significant role of chemotaxis and alarm-taxis coefficients in determining the existence and stability of predator-prey models, as discussed in the literature.
পত্রের ID : 2510.14728শিরোনাম : রাসায়নিক আকর্ষণ এবং সতর্কতা-আকর্ষণের সমন্বিত প্রভাব শিকারী-শিকার ব্যবস্থার বৈশ্বিক অস্তিত্ব এবং স্থিতিশীলতার উপরলেখক : Gnanasekaran Shanmugasundaram (NIT Tiruchirappalli), Jitraj Saha (NIT Tiruchirappalli), Rafael Díaz Fuentes (University of Cagliari)শ্রেণীবিভাগ : math.AP (আংশিক অবকল সমীকরণের বিশ্লেষণ)জমা দেওয়ার সময় : ২০২৫ সালের অক্টোবর ১৬ তারিখপত্রের লিঙ্ক : https://arxiv.org/abs/2510.14728 এই গবেষণা সীমাবদ্ধ ডোমেইন Ω ⊂ R n \Omega \subset \mathbb{R}^n Ω ⊂ R n এ মসৃণ সীমানা ∂ Ω \partial\Omega ∂ Ω সহ সমজাতীয় Neumann সীমানা শর্তের অধীনে সম্পূর্ণ প্যারাবোলিক শিকারী-শিকার রাসায়নিক-সতর্কতা আকর্ষণ ব্যবস্থা পরীক্ষা করে। নির্দিষ্ট প্যারামিটার শর্তের অধীনে, সিস্টেমের একটি অনন্য বৈশ্বিক সীমাবদ্ধ ধ্রুবক সমাধান বিদ্যমান প্রমাণ করা হয়েছে। উপযুক্ত Lyapunov ফাংশনাল নির্মাণের মাধ্যমে সমাধানের সংগতি প্রতিষ্ঠা করা হয়েছে। অতিরিক্তভাবে, সংখ্যাসূচক অনুকরণের মাধ্যমে সমাধানের অসিম্পটোটিক আচরণ যাচাই করা হয়েছে। ফলাফলগুলি রাসায়নিক আকর্ষণ এবং সতর্কতা-আকর্ষণ সহগের শিকারী-শিকার মডেলের অস্তিত্ব এবং স্থিতিশীলতা নির্ধারণে গুরুত্বপূর্ণ ভূমিকা তুলে ধরে।
সতর্কতা আহ্বান প্রক্রিয়া : গবেষণা সামুদ্রিক ইকোসিস্টেমে "চোর সতর্কতা" অনুমানের উপর ভিত্তি করে, যেখানে শিকার সংকেত প্রেরণ করে মাধ্যমিক শিকারীদের আকর্ষণ করে, যা প্রধান শিকারীদের হুমকি দেয়। সাধারণ উদাহরণ হল ডাইনোফ্ল্যাজেলেট যা কপেপড হস্তক্ষেপের প্রতিক্রিয়ায় জৈব আলোকসজ্জা সংকেত প্রেরণ করে, যা মাছকে মাধ্যমিক শিকারী হিসাবে আকর্ষণ করে।রাসায়নিক আকর্ষণ ঘটনা : সূক্ষ্মজীব রাসায়নিক গ্রেডিয়েন্টের প্রতিক্রিয়ায় গতিশীলতা, যা ইকোলজি, ইমিউনোলজি, জীববিজ্ঞান এবং চিকিৎসায় ব্যাপক প্রয়োগ রয়েছে।এই গবেষণা যে মূল সমস্যা সমাধান করতে চায় তা হল:
রাসায়নিক আকর্ষণ এবং সতর্কতা-আকর্ষণ উভয়ই অন্তর্ভুক্ত করে এমন তিন-প্রজাতির পারস্পরিক ক্রিয়া সিস্টেম কীভাবে গাণিতিকভাবে মডেল করা যায় কোন শর্তে সিস্টেমের সমাধান বৈশ্বিকভাবে বিদ্যমান এবং সীমাবদ্ধ থাকে সিস্টেমের দীর্ঘমেয়াদী অসিম্পটোটিক আচরণ কীভাবে হয় বিদ্যমান সাহিত্য প্রধানত নিম্নলিখিতগুলিতে ফোকাস করে:
একক আকর্ষণ প্রক্রিয়ার গবেষণা (বিশুদ্ধ রাসায়নিক আকর্ষণ বা বিশুদ্ধ সতর্কতা-আকর্ষণ) দ্বিমাত্রিক স্থানের স্থানীয় ফলাফল দুটি আকর্ষণ প্রক্রিয়া একসাথে বিবেচনা করার কঠোর গাণিতিক বিশ্লেষণের অভাব প্রথমবারের জন্য প্রতিষ্ঠা করা হয়েছে রাসায়নিক আকর্ষণ এবং সতর্কতা-আকর্ষণ উভয়ই অন্তর্ভুক্ত করে এমন শিকারী-শিকার সিস্টেমের বৈশ্বিক অস্তিত্ব, সীমাবদ্ধতা এবং সংগতির কঠোর গাণিতিক তত্ত্বউদ্ভাবনী Lyapunov ফাংশনাল নির্মাণ করা হয়েছে সিস্টেমের বৈশ্বিক অসিম্পটোটিক স্থিতিশীলতা বিশ্লেষণের জন্যবৈশ্বিক ধ্রুবক সমাধান প্রমাণ করা হয়েছে n ≥ 1 n \geq 1 n ≥ 1 মাত্রার স্থানে অস্তিত্ব এবং সমান সীমাবদ্ধতাসম্পূর্ণ স্থিতিশীলতা বিশ্লেষণ কাঠামো প্রতিষ্ঠা করা হয়েছে , যার মধ্যে সহ-অস্তিত্ব অবস্থা, বিলুপ্তি অবস্থা এবং আধা-সহ-অস্তিত্ব অবস্থার অসিম্পটোটিক আচরণ রয়েছেসংখ্যাসূচক যাচাইকরণ প্রদান করা হয়েছে , দ্বিমাত্রিক এবং ত্রিমাত্রিক স্থানে তাত্ত্বিক ফলাফল যাচাই করা হয়েছেগবেষণা নিম্নলিখিত সংযুক্ত অরৈখিক প্যারাবোলিক আংশিক অবকল সমীকরণ বিবেচনা করে:
{ u t = d 1 Δ u − χ 1 ∇ ⋅ ( u ∇ ( v w ) ) + μ 1 u ( 1 − u + a 1 v + a 2 w ) v t = d 2 Δ v − χ 2 ∇ ⋅ ( v ∇ z ) + μ 2 v ( 1 − v − a 3 u + a 4 w ) w t = d 3 Δ w + ξ ∇ ⋅ ( w ∇ z ) + μ 3 w ( 1 − w − a 5 u − a 6 v ) z t = d 4 Δ z + α v + β w − γ z \begin{cases}
u_t = d_1\Delta u - \chi_1\nabla \cdot (u\nabla(vw)) + \mu_1u(1-u+a_1v+a_2w) \\
v_t = d_2\Delta v - \chi_2\nabla \cdot (v\nabla z) + \mu_2v(1-v-a_3u+a_4w) \\
w_t = d_3\Delta w + \xi\nabla \cdot (w\nabla z) + \mu_3w(1-w-a_5u-a_6v) \\
z_t = d_4\Delta z + \alpha v + \beta w - \gamma z
\end{cases} ⎩ ⎨ ⎧ u t = d 1 Δ u − χ 1 ∇ ⋅ ( u ∇ ( v w )) + μ 1 u ( 1 − u + a 1 v + a 2 w ) v t = d 2 Δ v − χ 2 ∇ ⋅ ( v ∇ z ) + μ 2 v ( 1 − v − a 3 u + a 4 w ) w t = d 3 Δ w + ξ ∇ ⋅ ( w ∇ z ) + μ 3 w ( 1 − w − a 5 u − a 6 v ) z t = d 4 Δ z + αv + βw − γ z
যেখানে:
u ( x , t ) u(x,t) u ( x , t ) : মাধ্যমিক শিকারীর ঘনত্বv ( x , t ) v(x,t) v ( x , t ) : প্রধান শিকারীর ঘনত্বw ( x , t ) w(x,t) w ( x , t ) : শিকারের ঘনত্বz ( x , t ) z(x,t) z ( x , t ) : রাসায়নিক আকর্ষক ঘনত্বχ 1 \chi_1 χ 1 : সতর্কতা-আকর্ষণ সহগ (মাধ্যমিক শিকারী প্রধান শিকারী-শিকার ঘনত্ব গ্রেডিয়েন্টের দিকে গতিশীল)χ 2 \chi_2 χ 2 : রাসায়নিক আকর্ষণ সহগ (প্রধান শিকারী রাসায়নিক সংকেতের দিকে)ξ \xi ξ : রাসায়নিক বিকর্ষণ সহগ (শিকার রাসায়নিক সংকেত থেকে দূরে)d i d_i d i : বিস্তার সহগμ i \mu_i μ i : লজিস্টিক বৃদ্ধি সহগa i a_i a i : প্রজাতি-মধ্যস্থ পারস্পরিক ক্রিয়া সহগবহু-পদক্ষেপ পদ্ধতি ব্যবহার করা হয়:
প্রথম পদক্ষেপ : L p ( Ω ) L^p(\Omega) L p ( Ω ) নর্ম অনুমান প্রতিষ্ঠা করা (লেমা 3.1)দ্বিতীয় পদক্ষেপ : Neumann তাপ অর্ধ-গ্রুপ কৌশল ব্যবহার করে L ∞ ( Ω ) L^{\infty}(\Omega) L ∞ ( Ω ) সীমাবদ্ধতা অর্জন করা (লেমা 3.2)তৃতীয় পদক্ষেপ : গ্রেডিয়েন্ট অনুমান W 1 , ∞ ( Ω ) W^{1,\infty}(\Omega) W 1 , ∞ ( Ω ) প্রতিষ্ঠা করা (লেমা 3.4-3.6)বিভিন্ন সমতা অবস্থার জন্য সংশ্লিষ্ট শক্তি ফাংশনাল নির্মাণ করা হয়:
সহ-অস্তিত্ব অবস্থা :
E 1 ( t ) = ∫ Ω ( u − u ∗ − u ∗ ln u u ∗ ) + Γ 1 ∫ Ω ( v − v ∗ − v ∗ ln v v ∗ ) + Γ 2 ∫ Ω ( w − w ∗ − w ∗ ln w w ∗ ) + 1 2 ∫ Ω ( z − z ∗ ) 2 E_1(t) = \int_\Omega \left(u-u^*-u^*\ln\frac{u}{u^*}\right) + \Gamma_1\int_\Omega \left(v-v^*-v^*\ln\frac{v}{v^*}\right) + \Gamma_2\int_\Omega \left(w-w^*-w^*\ln\frac{w}{w^*}\right) + \frac{1}{2}\int_\Omega (z-z^*)^2 E 1 ( t ) = ∫ Ω ( u − u ∗ − u ∗ ln u ∗ u ) + Γ 1 ∫ Ω ( v − v ∗ − v ∗ ln v ∗ v ) + Γ 2 ∫ Ω ( w − w ∗ − w ∗ ln w ∗ w ) + 2 1 ∫ Ω ( z − z ∗ ) 2
যেখানে Γ 1 = μ 1 a 1 μ 2 a 3 \Gamma_1 = \frac{\mu_1a_1}{\mu_2a_3} Γ 1 = μ 2 a 3 μ 1 a 1 , Γ 2 = μ 1 a 2 μ 3 a 5 \Gamma_2 = \frac{\mu_1a_2}{\mu_3a_5} Γ 2 = μ 3 a 5 μ 1 a 2
μ > 0 \mu > 0 μ > 0 বিদ্যমান রয়েছে, যেমন যখন min { μ 2 , μ 3 } > μ \min\{\mu_2,\mu_3\} > \mu min { μ 2 , μ 3 } > μ , শর্ত সন্তুষ্ট অ-ঋণাত্মক প্রাথমিক ডেটার জন্য, সিস্টেম (1.1) একটি অনন্য ধ্রুবক সমাধান ( u , v , w , z ) (u,v,w,z) ( u , v , w , z ) বিদ্যমান রয়েছে এবং সমানভাবে সীমাবদ্ধ:
∥ u ( ⋅ , t ) ∥ L ∞ ( Ω ) + ∥ v ( ⋅ , t ) ∥ W 1 , ∞ ( Ω ) + ∥ w ( ⋅ , t ) ∥ W 1 , ∞ ( Ω ) + ∥ z ( ⋅ , t ) ∥ W 1 , q ( Ω ) ≤ C \|u(\cdot,t)\|_{L^{\infty}(\Omega)} + \|v(\cdot,t)\|_{W^{1,\infty}(\Omega)} + \|w(\cdot,t)\|_{W^{1,\infty}(\Omega)} + \|z(\cdot,t)\|_{W^{1,q}(\Omega)} \leq C ∥ u ( ⋅ , t ) ∥ L ∞ ( Ω ) + ∥ v ( ⋅ , t ) ∥ W 1 , ∞ ( Ω ) + ∥ w ( ⋅ , t ) ∥ W 1 , ∞ ( Ω ) + ∥ z ( ⋅ , t ) ∥ W 1 , q ( Ω ) ≤ C
গবেষণা চারটি শ্রেণীর অসিম্পটোটিক আচরণ প্রতিষ্ঠা করেছে:
সহ-অস্তিত্ব অবস্থা (উপপাদ্য 1.2): প্যারামিটার নির্দিষ্ট শর্ত সন্তুষ্ট করলে, সমাধান সহ-অস্তিত্ব সমতা অবস্থা ( u ∗ , v ∗ , w ∗ , z ∗ ) (u^*,v^*,w^*,z^*) ( u ∗ , v ∗ , w ∗ , z ∗ ) এ সূচকীয়ভাবে সংগত হয়মাধ্যমিক শিকারী একক অস্তিত্ব অবস্থা (উপপাদ্য 1.3): সমাধান ( 1 , 0 , 0 , 0 ) (1,0,0,0) ( 1 , 0 , 0 , 0 ) এ সংগত হয়আধা-সহ-অস্তিত্ব অবস্থা (উপপাদ্য 1.4):শিকার বিলুপ্তি অবস্থা: ( u ˉ , v ˉ , 0 , z ˉ ) (ū,v̄,0,z̄) ( u ˉ , v ˉ , 0 , z ˉ ) প্রধান শিকারী বিলুপ্তি অবস্থা: ( u ^ , 0 , w ^ , z ^ ) (û,0,ŵ,ẑ) ( u ^ , 0 , w ^ , z ^ ) স্থানিক ডোমেইন : [ − 0.5 , 0.5 ] n [-0.5,0.5]^n [ − 0.5 , 0.5 ] n (n = 2 , 3 n=2,3 n = 2 , 3 )গ্রিড : 2D এর জন্য 101 × 101 101×101 101 × 101 , 3D এর জন্য 50 × 50 × 50 50×50×50 50 × 50 × 50 প্রাথমিক শর্ত : রেডিয়াল সমরূপ ঘণ্টা আকৃতির ফাংশনপ্যারামিটার : d i = χ 1 = χ 2 = ξ = μ i = α = β = 1 d_i=\chi_1=\chi_2=\xi=\mu_i=\alpha=\beta=1 d i = χ 1 = χ 2 = ξ = μ i = α = β = 1 , γ = 2 \gamma=2 γ = 2 চারটি সাধারণ উদাহরণের মাধ্যমে তাত্ত্বিক পূর্বাভাস যাচাই করা হয়েছে:
উদাহরণ 5.1 : a i = 0.5 a_i=0.5 a i = 0.5 হলে, সমাধান সহ-অস্তিত্ব অবস্থা ( 1.286877 , 0.428973 , 0.144368 , 0.286793 ) (1.286877, 0.428973, 0.144368, 0.286793) ( 1.286877 , 0.428973 , 0.144368 , 0.286793 ) এ সংগত হয়উদাহরণ 5.2 : নির্দিষ্ট প্যারামিটারে সমাধান ( 1 , 0 , 0 , 0 ) (1,0,0,0) ( 1 , 0 , 0 , 0 ) এ সংগত হয়উদাহরণ 5.3-5.4 : দুটি আধা-সহ-অস্তিত্ব অবস্থার সংগতি যাচাই করা হয়েছেঅরৈখিক পদ প্রক্রিয়াকরণ : সতর্কতা-আকর্ষণ পদ χ 1 ∇ ⋅ ( u ∇ ( v w ) ) \chi_1\nabla \cdot (u\nabla(vw)) χ 1 ∇ ⋅ ( u ∇ ( v w )) এর জটিলতাসংযুক্ত সিস্টেম বিশ্লেষণ : চারটি সমীকরণের শক্তিশালী সংযোগ প্রকৃতিগ্রেডিয়েন্ট অনুমান : উচ্চ-মাত্রিক স্থানে প্রযুক্তিগত জটিলতাGagliardo-Nirenberg অসমতা Neumann তাপ অর্ধ-গ্রুপ তত্ত্ব সর্বোচ্চ Sobolev নিয়মিততা Young অসমতা এবং Hölder অসমতা Keller-Segel মডেল (1970): রাসায়নিক আকর্ষণের ভিত্তি কাজHaskell-Bell মডেল (2021): সতর্কতা-আকর্ষণের প্রথম গাণিতিক মডেলিংবিদ্যমান সম্প্রসারণ : প্রধানত একক প্রক্রিয়া বা নিম্ন-মাত্রিক ক্ষেত্রে কেন্দ্রীভূতবিদ্যমান কাজের তুলনায়, এই পত্র প্রথমবারের জন্য:
দুটি আকর্ষণ প্রক্রিয়ার সমন্বিত প্রভাব কঠোরভাবে প্রক্রিয়া করেছে উচ্চ-মাত্রিক স্থানের বৈশ্বিক তত্ত্ব প্রতিষ্ঠা করেছে সম্পূর্ণ স্থিতিশীলতা বিশ্লেষণ কাঠামো প্রদান করেছে উপযুক্ত প্যারামিটার শর্তের অধীনে, chemo-alarm-taxis সিস্টেম বৈশ্বিক সীমাবদ্ধ ধ্রুবক সমাধান বিদ্যমান রয়েছে সিস্টেমের দীর্ঘমেয়াদী আচরণ প্যারামিটার সম্পর্ক দ্বারা নির্ধারিত হয়, চারটি ভিন্ন সমতা অবস্থায় সংগত হতে পারে রাসায়নিক আকর্ষণ এবং সতর্কতা-আকর্ষণ সহগ সিস্টেম স্থিতিশীলতা নির্ধারণে মূল ভূমিকা পালন করে "চোর সতর্কতা" অনুমানের গাণিতিক সম্ভাব্যতা যাচাই করেছে প্রাকৃতিক জগতে জটিল তিন-প্রজাতির পারস্পরিক ক্রিয়া গতিশীলতা ব্যাখ্যা করেছে ইকোলজি সংরক্ষণ কৌশলের জন্য তাত্ত্বিক ভিত্তি প্রদান করেছে প্যারামিটার শর্ত : স্থিতিশীলতা ফলাফল শক্তিশালী প্যারামিটার সীমাবদ্ধতা প্রয়োজন করেরৈখিকীকরণ বিশ্লেষণ : স্থানীয় স্থিতিশীলতা বিশ্লেষণ রৈখিকীকরণের উপর ভিত্তি করেসংখ্যাসূচক যাচাইকরণ : শুধুমাত্র নির্দিষ্ট প্যারামিটার পরিসরে যাচাইকরণ করা হয়েছেপ্যারামিটার শর্তের সীমাবদ্ধতা শিথিল করা আরও জটিল প্রতিক্রিয়া পদ বিবেচনা করা স্থানিক প্যাটার্ন গঠন গবেষণা করা র্যান্ডম পরিবেশে সম্প্রসারণ করা তাত্ত্বিক কঠোরতা : সম্পূর্ণ গাণিতিক প্রমাণ কাঠামো প্রদান করেছেপদ্ধতি উদ্ভাবনী : একাধিক বিশ্লেষণ কৌশল চতুরভাবে সংযুক্ত করেছেফলাফল সম্পূর্ণতা : অস্তিত্ব, সীমাবদ্ধতা এবং স্থিতিশীলতা অন্তর্ভুক্ত করেছেসংখ্যাসূচক যাচাইকরণ : 2D এবং 3D সংখ্যাসূচক পরীক্ষা তাত্ত্বিক ফলাফল সমর্থন করেছেপ্যারামিটার শর্ত জটিল : বাস্তব প্রয়োগে প্যারামিটার যাচাইকরণ কঠিনজৈব বাস্তবতা : কিছু গাণিতিক অনুমান অত্যধিক আদর্শায়িত হতে পারেগণনা জটিলতা : উচ্চ-মাত্রিক ক্ষেত্রে সংখ্যাসূচক গণনা চ্যালেঞ্জতাত্ত্বিক অবদান : আকর্ষণ সিস্টেম গবেষণার জন্য নতুন বিশ্লেষণ কাঠামো প্রদান করেছেপ্রয়োগ মূল্য : ইকোলজি গাণিতিক মডেলিংয়ের জন্য গুরুত্বপূর্ণ সরঞ্জাম প্রদান করেছেপদ্ধতি মূল্য : Lyapunov ফাংশনাল নির্মাণ পদ্ধতি সম্প্রসারণযোগ্যসামুদ্রিক ইকোসিস্টেম মডেলিং কোষ জীববিজ্ঞানে সংকেত প্রবাহ ইমিউন সিস্টেম গতিশীলতা বিশ্লেষণ ইকোলজি সংরক্ষণ কৌশল প্রণয়ন পত্রটি 46টি সম্পর্কিত সংদর্ভ উদ্ধৃত করেছে, যা রাসায়নিক আকর্ষণ তত্ত্ব, সতর্কতা-আকর্ষণ মডেলিং, আংশিক অবকল সমীকরণ তত্ত্ব এবং অন্যান্য একাধিক ক্ষেত্রের গুরুত্বপূর্ণ কাজ অন্তর্ভুক্ত করে, যা গবেষণার জন্য দৃঢ় তাত্ত্বিক ভিত্তি প্রদান করে।
সামগ্রিক মূল্যায়ন : এটি গাণিতিক জীববিজ্ঞান ক্ষেত্রে গুরুত্বপূর্ণ তাত্ত্বিক তাৎপর্য সহ একটি উচ্চ-মানের পত্র, যা প্রথমবারের জন্য chemo-alarm-taxis সিস্টেমের সম্পূর্ণ গাণিতিক তত্ত্ব প্রতিষ্ঠা করেছে, পদ্ধতি কঠোর, ফলাফল গভীর, এবং সম্পর্কিত ক্ষেত্রের উন্নয়নে গুরুত্বপূর্ণ প্রচারমূলক ভূমিকা রয়েছে।