2025-11-15T19:55:11.527544

Solid-state dewetting of axisymmetric thin film on axisymmetric curved-surface substrates: modeling and simulation

Duan, Li, Zhou

In this work, we consider the solid-state dewetting of an axisymmetric thin film on a curved-surface substrate, with the assumption that the substrate morphology is also axisymmetric. Under the assumptions of axisymmetry, the surface evolution problem on a curved-surface substrate can be reduced to a curve evolution problem on a static curved substrate. Based on the thermodynamic variation of the anisotropic surface energy, we thoroughly derive a sharp-interface model that is governed by anisotropic surface diffusion, along with appropriate boundary conditions. The continuum system satisfies the laws of energy decay and volume conservation, which motivates the design of a structure-preserving numerical algorithm for simulating the mathematical model. By introducing a symmetrized surface energy matrix, we derive a novel symmetrized variational formulation. Then, by carefully discretizing the boundary terms of the variational formulation, we establish an unconditionally energy-stable parametric finite element approximation of the axisymmetric system. By applying an ingenious correction method, we further develop another structure-preserving method that can preserve both the energy stability and volume conservation properties. Finally, we present extensive numerical examples to demonstrate the convergence and structure-preserving properties of our proposed numerical scheme. Additionally, several interesting phenomena are explored, including the migration of 'small' particles on a curved-surface substrate generated by curves with positive or negative curvature, pinch-off events, and edge retraction.

academic

अक्षसममित पतली फिल्म का ठोस-अवस्था विआर्द्रीकरण अक्षसममित वक्र-सतह सब्सट्रेट पर: मॉडलिंग और सिमुलेशन

मूल जानकारी

पेपर ID: 2501.00783
शीर्षक: अक्षसममित पतली फिल्म का ठोस-अवस्था विआर्द्रीकरण अक्षसममित वक्र-सतह सब्सट्रेट पर: मॉडलिंग और सिमुलेशन
लेखक: Zhenghua Duan, Meng Li*, Chunjie Zhou (झेंगझोऊ विश्वविद्यालय गणित और सांख्यिकी संस्थान)
वर्गीकरण: math.NA cs.NA
प्रकाशन समय: 3 जनवरी 2025 (arXiv प्रीप्रिंट)
पेपर लिंक: https://arxiv.org/abs/2501.00783

सारांश

यह पेपर अक्षसममित वक्र-सतह सब्सट्रेट पर अक्षसममित पतली फिल्म की ठोस-अवस्था विआर्द्रीकरण घटना का अध्ययन करता है। अक्षसममित मान्यता के तहत, वक्र-सतह सब्सट्रेट पर सतह विकास समस्या को स्थिर वक्र सब्सट्रेट पर वक्र विकास समस्या में सरल किया जा सकता है। विषमदिश सतह ऊर्जा के थर्मोडायनामिक भिन्नता के आधार पर, लेखकों ने विषमदिश सतह विसरण द्वारा नियंत्रित तीव्र इंटरफेस मॉडल और संबंधित सीमा शर्तों को कठोरता से प्राप्त किया है। सतत प्रणाली ऊर्जा क्षय और आयतन संरक्षण नियमों को संतुष्ट करती है, जो संरचना-संरक्षण संख्यात्मक एल्गोरिदम डिजाइन करने के लिए प्रेरणा प्रदान करता है। सममित सतह ऊर्जा मैट्रिक्स का परिचय देकर, एक नवीन सममित भिन्नता सूत्र प्राप्त किया गया है। भिन्नता सूत्र की सीमा शर्तों के सावधानीपूर्वक विवेकीकरण के माध्यम से, अक्षसममित प्रणाली के लिए बिना शर्त ऊर्जा-स्थिर पैरामीट्रिक परिमित तत्व सन्निकटन स्थापित किया गया है। चतुर संशोधन विधि लागू करके, ऊर्जा स्थिरता और आयतन संरक्षण दोनों गुणों को बनाए रखने वाली संरचना-संरक्षण विधि विकसित की गई है।

अनुसंधान पृष्ठभूमि और प्रेरणा

समस्या की महत्ता

ठोस-अवस्था विआर्द्रीकरण (Solid-State Dewetting, SSD) एक ऐसी घटना है जिसमें ठोस पतली फिल्म पिघलने बिंदु से बहुत कम तापमान पर अस्थिर हो जाती है, फिल्म विआर्द्र हो जाती है या एकत्रित हो जाती है, जिससे जटिल आकृति परिवर्तन होता है और सब्सट्रेट पर छोटे कणों का निर्माण होता है। यह घटना प्रकाशीय और चुंबकीय उपकरणों, पतली फिल्मों, सेंसर और उत्प्रेरक निर्माण जैसे क्षेत्रों में व्यापक अनुप्रयोग रखती है।

मौजूदा अनुसंधान की सीमाएं

आयामी प्रतिबंध: मौजूदा सैद्धांतिक अनुसंधान मुख्य रूप से समतल सब्सट्रेट पर केंद्रित है, टोपोलॉजिकल पैटर्न वाले सब्सट्रेट पर कम ध्यान दिया गया है
द्विआयामी सीमा: मौजूदा अनुसंधान मुख्य रूप से द्विआयामी परिदृश्य पर केंद्रित है, त्रिआयामी वक्र-सतह सब्सट्रेट पर SSD का अनुसंधान सीमित है
संख्यात्मक विधि की कमी: ऊर्जा स्थिरता और आयतन संरक्षण दोनों को बनाए रखने वाली कुशल संख्यात्मक विधियों की कमी है

अनुसंधान प्रेरणा

त्रिआयामी पतली फिल्मों और वक्र-सतह सब्सट्रेट में आमतौर पर घूर्णन समरूपता होती है, जो जटिल त्रिआयामी SSD को वक्र सब्सट्रेट पर वक्र-आधारित प्रणाली में सरल कर सकती है
मजबूत विषमदिश प्रभावों को संभालने में सक्षम संख्यात्मक विधियों का विकास आवश्यक है
वक्र सब्सट्रेट की ज्यामितीय विशेषताएं पतली फिल्म के विकास व्यवहार पर महत्वपूर्ण प्रभाव डालती हैं

मुख्य योगदान

सैद्धांतिक मॉडलिंग: थर्मोडायनामिक भिन्नता के माध्यम से अक्षसममित वक्र-सतह सब्सट्रेट पर त्रिआयामी अक्षसममित पतली फिल्म के SSD का तीव्र इंटरफेस मॉडल प्राप्त किया गया, जो कमजोर और मजबूत विषमदिश दोनों स्थितियों को शामिल करता है
संख्यात्मक विधि नवाचार:
- सममित सतह ऊर्जा मैट्रिक्स का परिचय, नवीन सममित भिन्नता सूत्र प्राप्त किया
- सीमा शर्तों के सावधानीपूर्वक विवेकीकरण के माध्यम से, बिना शर्त ऊर्जा-स्थिर पैरामीट्रिक परिमित तत्व सन्निकटन डिजाइन किया
- आयतन संरक्षण और ऊर्जा स्थिरता दोनों को बनाए रखने वाली संरचना-संरक्षण विधि विकसित की
एल्गोरिदम सत्यापन: व्यापक संख्यात्मक प्रयोगों के माध्यम से प्रस्तावित विधि के अभिसरण और संरचना-संरक्षण गुणों को सत्यापित किया
भौतिक घटना अन्वेषण: सकारात्मक/नकारात्मक वक्रता सब्सट्रेट पर छोटे कणों का प्रवास, जकड़न घटनाएं और किनारे संकुचन सहित विभिन्न दिलचस्प घटनाओं का अध्ययन किया

विधि विवरण

कार्य परिभाषा

अक्षसममित वक्र-सतह सब्सट्रेट पर अक्षसममित पतली फिल्म के ठोस-अवस्था विआर्द्रीकरण विकास प्रक्रिया का अध्ययन, जिसमें शामिल है:

इनपुट: प्रारंभिक पतली फिल्म आकृति, सब्सट्रेट ज्यामिति, सामग्री पैरामीटर
आउटपुट: समय के साथ पतली फिल्म का विकास प्रक्षेपवक्र
बाधाएं: आयतन संरक्षण, ऊर्जा क्षय, संपर्क रेखा सीमा शर्तें

गणितीय मॉडल

कुल मुक्त ऊर्जा

प्रणाली की कुल मुक्त ऊर्जा को इस प्रकार व्यक्त किया जाता है: $W = \int\int_S \gamma_{FV}(N) dS + (\gamma_{FS} - \gamma_{VS}) A(\Gamma_o/\Gamma_i)$

जहां $S$ पतली फिल्म/वाष्प इंटरफेस है, $\Gamma_i$ और $\Gamma_o$ आंतरिक और बाहरी संपर्क रेखाएं हैं।

अक्षसममित सरलीकरण

अक्षसममित मान्यता के तहत, सतह को इस प्रकार पैरामीट्रिज किया जा सकता है: $S(s,\phi) := (r(s)\cos\phi, r(s)\sin\phi, z(s))$

कुल ऊर्जा को सरल किया जाता है: $W = 2\pi\int_0^L \gamma(\theta)r ds + 2\pi(\gamma_{FS} - \gamma_{VS})\int_{c_l}^{c_r} \hat{x}(c)dc$

तीव्र इंटरफेस मॉडल

थर्मोडायनामिक भिन्नता के माध्यम से नियंत्रण समीकरण प्राप्त किया जाता है: $X_t \cdot n = \frac{1}{r}(r\mu_s)_s, \quad 0 < s < L(t), t > 0$

रासायनिक क्षमता को इस प्रकार परिभाषित किया जाता है: $\mu = (\gamma(\theta) + \gamma''(\theta))\kappa - \frac{\gamma(\theta)z_s + \gamma'(\theta)r_s}{r}$

सीमा शर्तों में शामिल हैं:

संपर्क रेखा शर्त: सुनिश्चित करता है कि संपर्क रेखा वक्र सब्सट्रेट पर रहती है
शिथिल संपर्क कोण शर्त: संपर्क कोण के गतिशील विकास का वर्णन करता है
शून्य द्रव्यमान प्रवाह शर्त: आयतन संरक्षण सुनिश्चित करता है

सममित भिन्नता सूत्र

सममित मैट्रिक्स का परिचय: $B(\theta) = \begin{pmatrix} \gamma(\theta) & -\gamma'(\theta) \\ \gamma'(\theta) & \gamma(\theta) \end{pmatrix} \begin{pmatrix} \cos 2\theta & \sin 2\theta \\ \sin 2\theta & -\cos 2\theta \end{pmatrix} + S(\theta)\left[\frac{1}{2}I - \frac{1}{2}\begin{pmatrix} \cos 2\theta & \sin 2\theta \\ \sin 2\theta & -\cos 2\theta \end{pmatrix}\right]$

इस मैट्रिक्स के माध्यम से, रासायनिक क्षमता समीकरण को इस प्रकार लिखा जा सकता है: $r\mu n = \partial_s[rB(\theta)\partial_s X] - \gamma(\theta)e_1$

संख्यात्मक विधि

पैरामीट्रिक परिमित तत्व विधि

भिन्नता सूत्र को हल करने के लिए पैरामीट्रिक परिमित तत्व विधि (PFEM) का उपयोग किया जाता है:

समय विवेकीकरण: अर्ध-निहित प्रारूप
स्थान विवेकीकरण: रैखिक परिमित तत्व
सीमा प्रसंस्करण: सीमा शर्तों का सावधानीपूर्वक विवेकीकरण

ऊर्जा-स्थिर विधि

डिजाइन की गई संख्यात्मक योजना असतत ऊर्जा स्थिरता को संतुष्ट करती है: $W(X^{m+1}) - W(X^m) \leq 0$

संरचना-संरक्षण विधि

संशोधन पद $\delta f^{m+\frac{1}{2}}$ के माध्यम से सटीक आयतन संरक्षण प्राप्त किया जाता है: $V(X^{m+1}) - V(X^m) = 0$

प्रयोगात्मक सेटअप

परीक्षण मामले

पेपर तीन विशिष्ट प्रारंभिक विन्यास पर विचार करता है:

गोलाकार पतली फिल्म: अर्धगोलाकार सब्सट्रेट पर त्रिज्या 1.5 की गोलाकार पतली फिल्म
वलयाकार पतली फिल्म: अर्धगोलाकार सब्सट्रेट पर मोटाई 0.5 की अक्षसममित वलयाकार पतली फिल्म
गोलाकार पतली फिल्म: गोलाकार सब्सट्रेट पर त्रिज्या 1.8 की गोलाकार पतली फिल्म

सामग्री पैरामीटर

संपर्क रेखा प्रवासन दर: $\eta = 100$
सतह ऊर्जा पैरामीटर: $\sigma = -\frac{\sqrt{3}}{2}$
विषमदिश सतह ऊर्जा: $\gamma(\theta) = 1 + \beta\cos(4\theta)$ $γ (θ) = 1 + β cos (4 θ)$
- $\beta = 0$ : समदिश
- $0 \leq \beta \leq \frac{1}{15}$ : कमजोर विषमदिश
- $\beta > \frac{1}{15}$ : मजबूत विषमदिश

संख्यात्मक पैरामीटर

जाल आकार: $h = 2^{-7}$ से $2^{-8}$
समय चरण: $\Delta t = 2^{-9}$ से $2^{-10}$
अभिसरण सहनशीलता: $10^{-8}$

प्रयोगात्मक परिणाम

अभिसरण सत्यापन

मैनिफोल्ड दूरी मीट्रिक के माध्यम से संख्यात्मक त्रुटि को मापा जाता है: $e^{h,\Delta t}(t) := |(\Omega^{h,\Delta t}(t)\backslash\Omega^r(t)) \cup (\Omega^r(t)\backslash\Omega^{h,\Delta t}(t))|$

परिणाम दिखाते हैं:

स्थान अभिसरण क्रम: द्वितीय-क्रम अभिसरण, सैद्धांतिक अपेक्षा के अनुरूप
समय अभिसरण क्रम: प्रथम-क्रम अभिसरण
समदिश, कमजोर विषमदिश और मजबूत विषमदिश सभी स्थितियों में अच्छी अभिसरण बनाए रखता है

संरचना-संरक्षण गुणों का सत्यापन

ऊर्जा स्थिरता

संख्यात्मक प्रयोग असतत ऊर्जा के एकदिष्ट ह्रास की पुष्टि करते हैं:

विभिन्न समय चरणों के तहत ऊर्जा वक्र सभी घटते हुए प्रवृत्ति दिखाते हैं
विषमदिश तीव्रता ऊर्जा स्थिरता को प्रभावित नहीं करती है

आयतन संरक्षण

ऊर्जा-स्थिर विधि: सापेक्ष आयतन त्रुटि लगभग $10^{-3}$ क्रम की है
संरचना-संरक्षण विधि: मशीन परिशुद्धता पर आयतन संरक्षण प्राप्त करता है (त्रुटि $\sim 10^{-15}$ )

भौतिक घटना अनुसंधान

कण प्रवास घटना

वक्र सब्सट्रेट पर देखी गई:

सकारात्मक वक्रता सब्सट्रेट: छोटे कण कम वक्रता वाले क्षेत्रों में प्रवास करते हैं
नकारात्मक वक्रता सब्सट्रेट: कण व्यवहार विषमदिश तीव्रता से प्रभावित होता है
विषमदिश प्रभाव: मजबूत विषमदिश प्रवास दर को महत्वपूर्ण रूप से प्रभावित करता है

जकड़न घटना

लंबी पतली फिल्म विकास प्रक्रिया में जकड़न से गुजरती है:

जकड़न के बाद दो स्वतंत्र भागों में पुनः आरंभ किया जाता है
अंततः छोटे द्वीप और वलयाकार पतली फिल्मों का निर्माण होता है

किनारे संकुचन

अर्ध-अनंत सीढ़ीदार पतली फिल्म का किनारे संकुचन:

पतली फिल्म धीरे-धीरे सब्सट्रेट के कोने पर चढ़ती है
विषमदिश तीव्रता जितनी अधिक होती है, संकुचन दर उतनी ही धीमी होती है

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

अक्षसममित वक्र सब्सट्रेट पर SSD का संपूर्ण गणितीय मॉडल सफलतापूर्वक स्थापित किया गया
ऊर्जा स्थिरता और आयतन संरक्षण दोनों को बनाए रखने वाली कुशल संख्यात्मक विधि विकसित की गई
वक्र सब्सट्रेट ज्यामिति के पतली फिल्म विकास व्यवहार पर महत्वपूर्ण प्रभाव का पता चला

सीमाएं

अक्षसममित मान्यता: संभाले जा सकने वाली ज्यामितीय आकृतियों को सीमित करता है
सतह विसरण प्रभुत्व: लोचदार प्रभाव और रासायनिक प्रतिक्रिया को नजरअंदाज करता है
कम्प्यूटेशनल जटिलता: त्रिआयामी समस्या की कम्प्यूटेशनल लागत अभी भी अधिक है

भविष्य की दिशाएं

गैर-अक्षसममित ज्यामिति तक विस्तार
बहु-भौतिकी युग्मन प्रभावों पर विचार
स्व-अनुकूली जाल विधि विकसित करना
वास्तविक औद्योगिक समस्याओं में अनुप्रयोग

गहन मूल्यांकन

शक्तियां

सैद्धांतिक कठोरता: थर्मोडायनामिक्स के प्रथम सिद्धांत से शुरू करके, संपूर्ण और कठोर व्युत्पत्ति
विधि नवाचार: सममित भिन्नता सूत्र और संरचना-संरक्षण प्रारूप में नवाचार है
संख्यात्मक स्थिरता: बिना शर्त ऊर्जा स्थिरता, समय चरण प्रतिबंध से बचा जाता है
भौतिक महत्व: संख्यात्मक परिणाम समृद्ध भौतिक घटनाओं को प्रकट करते हैं
प्रयोगात्मक पूर्णता: अभिसरण, स्थिरता और अनुप्रयोग के कई पहलुओं को शामिल करता है

कमियां

अनुप्रयोग सीमा: अक्षसममित मान्यता विधि की सार्वभौमिकता को सीमित करती है
कम्प्यूटेशनल दक्षता: बड़े पैमाने की त्रिआयामी समस्याओं के लिए कम्प्यूटेशनल लागत अधिक हो सकती है
पैरामीटर संवेदनशीलता: कुछ सामग्री पैरामीटर के चयन में विस्तृत चर्चा की कमी है
प्रायोगिक सत्यापन: वास्तविक प्रयोगात्मक डेटा के साथ तुलना की कमी है

प्रभाव

शैक्षणिक योगदान: वक्र सब्सट्रेट पर SSD अनुसंधान के लिए महत्वपूर्ण सैद्धांतिक और संख्यात्मक उपकरण प्रदान करता है
व्यावहारिक मूल्य: नैनो-उपकरण डिजाइन और पतली फिल्म इंजीनियरिंग में अनुप्रयोग किया जा सकता है
पद्धति महत्व: संरचना-संरक्षण संख्यात्मक विधि अन्य ज्यामितीय विकास समस्याओं तक विस्तारित की जा सकती है

लागू परिदृश्य

अक्षसममित विशेषताओं वाले पतली फिल्म उपकरण डिजाइन
वक्र सतह पर नैनो संरचना तैयारी
उत्प्रेरक कण आकृति नियंत्रण
प्रकाशीय और चुंबकीय उपकरण सतह इंजीनियरिंग

संदर्भ

पेपर में 55 उच्च-गुणवत्ता संदर्भ हैं, जो SSD सैद्धांतिक मॉडलिंग, संख्यात्मक विधि और प्रायोगिक अनुसंधान के मुख्य कार्यों को शामिल करते हैं, साहित्य सर्वेक्षण व्यापक और प्राधिकृत है।

समग्र मूल्यांकन: यह एक उच्च-गुणवत्ता की संख्यात्मक विश्लेषण पेपर है, जिसमें सैद्धांतिक मॉडलिंग, संख्यात्मक विधि और अनुप्रयोग सत्यापन के सभी पहलुओं में महत्वपूर्ण योगदान है। विधि में अच्छा गणितीय आधार और संख्यात्मक गुण हैं, जो त्रिआयामी वक्र सब्सट्रेट पर ठोस-अवस्था विआर्द्रीकरण अनुसंधान के लिए शक्तिशाली उपकरण प्रदान करता है।