একটি সমজাতীয় এবং প্রতিসম কোষীয় ভ্রূণ কীভাবে নিজের প্রাকৃতিক প্রতিসাম্য ভেঙে একটি জেব্রার গায়ের ডোরাকাটা দাগ, চিতাবাঘের চামড়ার ছোপ অথবা মেরুদণ্ডী প্রাণীর হাতের আঙুলের সুনির্দিষ্ট বিন্যাস তৈরি করে? ১৯৫২ সালে তাঁর যুগান্তকারী গবেষণা প্রবন্ধ “দ্য কেমিক্যাল বেসিস অব মরফোজেনেসিস”-এ প্রখ্যাত গণিতবিদ অ্যালান টুরিং একটি বিস্ময়কর সত্য প্রমাণ করেন: সাধারণত ব্যাপন (diffusion) কোনো ব্যবস্থাকে সমসত্ত্ব এবং বিশৃঙ্খল করে তুললেও, রাসায়নিক বিক্রিয়া ও ভিন্ন ভিন্ন হারের ব্যাপনের যুগলবন্দী স্বয়ংক্রিয়ভাবে স্থানিক প্রতিসাম্য ভেঙে দিয়ে জটিল ম্যাক্রোস্কোপিক প্যাটার্ন বা নকশা তৈরি করতে পারে।
১. রিঅ্যাকশন-ডিফিউশন গাণিতিক কাঠামো
ধরা যাক দুটি রাসায়নিক পদার্থ বা মরফোজেনের মধ্যে মিথস্ক্রিয়া ঘটছে: একটি অ্যাক্টিভেটর $u(mathbf{r}, t)$ এবং একটি ইনহিবিটর $v(mathbf{r}, t)$। এদের স্থানিক ও কালিক বিবর্তন নির্ধারিত হয় নন-লিনিয়ার আংশিক ব্যবকলনীয় সমীকরণ দ্বারা:
[
begin{aligned}
frac{partial u}{partial t} &= f(u, v) + D_u nabla^2 u \
frac{partial v}{partial t} &= g(u, v) + D_v nabla^2 v
end{aligned}
]
এখানে $f(u, v)$ এবং $g(u, v)$ হলো নন-লিনিয়ার রাসায়নিক গতিবিদ্যা, $nabla^2$ হলো ল্যাপ্লাসিয়ান অপারেটর, এবং $D_u, D_v$ হলো যথাক্রমে অ্যাক্টিভেটর ও ইনহিবিটরের ব্যাপন গুণাঙ্ক (diffusion coefficients)।
২. টুরিং ইনস্ট্যাবিলিটির শর্তাবলী
ধরা যাক সিস্টেমটির একটি সমসত্ত্ব স্থিতাবস্থা $(u_0, v_0)$ আছে যেখানে $f(u_0, v_0) = 0$ এবং $g(u_0, v_0) = 0$। টুরিং ইনস্ট্যাবিলিটি ঘটার জন্য দুটি মৌলিক শর্ত পূরণ হতে হয়:
- ব্যাপনহীন অবস্থায় স্থানীয় স্থায়িত্ব ($D_u = D_v = 0$): জ্যাকবিয়ান ম্যাট্রিক্স $J_0$-এর ট্রেস ঋণাত্মক এবং ডিটারমিন্যান্ট ধনাত্মক হতে হবে:
[
operatorname{Tr}(J_0) = f_u + g_v 0
] - ব্যাপনজনিত অস্থিতিশীলতা ($D_u ne D_v$): স্থানিক তরঙ্গদৈর্ঘ্যের扰 perturbing mode $e^{lambda t + i mathbf{k} cdot mathbf{r}}$ যোগ করলে, একটি নির্দিষ্ট ওয়েভনাম্বার $k > 0$-এর জন্য গ্রোথ রেট $lambda(k^2) > 0$ হয়ে যায়। এর জন্য প্রয়োজনীয় শর্ত হলো:
[
D_v f_u + D_u g_v > 2 sqrt{D_u D_v det(J_0)} > 0
]
যেহেতু $f_u + g_v 0$) হতে হবে, ইনহিবিটরকে স্ব-বাধা প্রদানকারী ($g_v < 0$) হতে হবে এবং ইনহিবিটরের ব্যাপন হার অ্যাক্টিভেটরের তুলনায় অনেক বেশি হতে হবে:
[
frac{D_v}{D_u} gg 1 quad (text{স্বল্প-পাল্লার অ্যাক্টিভেশন এবং দীর্ঘ-পাল্লার ইনহিবিশন})
]
৩. তরঙ্গদৈর্ঘ্য নির্বাচন ও প্যাটার্ন রূপতত্ত্ব
যখন ব্যাপন অনুপাত একটি সংকটকালীন মান অতিক্রম করে, তখন সবচেয়ে দ্রুত বর্ধনশীল তরঙ্গসংখ্যা $k_{max}$ প্যাটার্নের বৈশিষ্ট্যসূচক তরঙ্গদৈর্ঘ্য $Lambda = 2pi / k_{max}$ নির্ধারণ করে। ডোমেইনের জ্যামিতিক সীমানা এবং নন-লিনিয়ার টার্মের ওপর নির্ভর করে সিস্টেমটি স্বতঃস্ফূর্তভাবে গোলাকার ছোপ (spots), ডোরাকাটা দাগ (stripes) অথবা গোলকধাঁধার মতো জটিল ল্যাবাইরিন্থাইন কাঠামো তৈরি করে।
