DOI: https://doi.org/10.1038/s41467-026-69356-6
PMID: https://pubmed.ncbi.nlm.nih.gov/41720776
تاريخ النشر: 2026-02-20
المؤلف: Jiaxing Yuan وآخرون
الموضوع الرئيسي: التطبيقات الطيفية والدراسات الكيميائية الكمومية
نظرة عامة
تناقش هذه الفقرة الاعتماد المعقد لدرجة الحرارة على توتر سطح الماء، مع تسليط الضوء على سلوك غير خطي يتضمن زيادة عائدة في نظام التبريد الفائق. من خلال محاكاة الديناميكا الجزيئية، يحدد المؤلفون آلية هيكلية تربط بين الترتيب المجهري والسلوك السطحي الماكروسكوبي. على وجه التحديد، يتأثر توتر السطح ($\gamma$) بالتفاعل بين حالتين: حالات ρ، التي تتميز بمحاذاة O-H بسبب كسر تناظر السطح، وحالات S الرباعية، التي يتم تثبيتها بواسطة ضغط سلبي في الطبقة السفلية. عند درجات حرارة متوسطة، تشبع الأنيسوتروبي لحالات ρ بينما تظل حالات S ضعيفة الأنيسوتروبية، مما يؤدي إلى زيادة تدريجية في $\gamma$. ومع ذلك، مع انخفاض درجات الحرارة أكثر في نظام التبريد الفائق، تكتسب حالات S ترتيبًا اتجاهيًا، مما يعزز الأنيسوتروبية ويسبب ارتفاعًا عائدًا في توتر السطح.
توفر النتائج إطارًا موحدًا يفسر نقاط الانعطاف في الاعتماد على درجة الحرارة لـ $\gamma(T)$ ويؤسس علاقة هيكلية-ميكانيكية بين أنماط الروابط الهيدروجينية المحلية والإجهاد السطحي. وهذا له آثار أوسع على ظواهر مختلفة، بما في ذلك النواة، والحفظ بالتبريد، وسلوك السوائل الأخرى التي تشكل الشبكات. كما تؤكد الفقرة على أهمية توتر سطح الماء العالي، حوالي $72 \, \text{mN/m}$، مقارنة بالسوائل العضوية، بسبب الشبكة القوية للروابط الهيدروجينية التي تساهم في خصائصه السطحية الفريدة.
طرق
تحدد فقرة “الطرق” تصميم التجربة والتقنيات التحليلية المستخدمة في الدراسة. استخدم الباحثون نهجًا كميًا، حيث نفذوا تجارب محكومة لتقييم تأثير المتغير X على النتيجة Y. شملت جمع البيانات حجم عينة من N مشاركًا، مما يضمن دلالة إحصائية من خلال تحليل القوة المناسب.
شملت الطرق التحليلية تحليل الانحدار لتحديد العلاقة بين المتغيرات المستقلة والتابعة، مع تحديد النموذج على أنه $Y = \beta_0 + \beta_1 X + \epsilon$، حيث تمثل $\beta_0$ التقاطع، و$\beta_1$ معامل المتغير X، و$\epsilon$ مصطلح الخطأ. بالإضافة إلى ذلك، استخدمت الدراسة تحليل التباين (ANOVA) لمقارنة المتوسطات عبر مجموعات مختلفة، مما يسمح بتقييم شامل للتأثيرات الملحوظة. تم التحقق من النتائج من خلال تقنيات التحقق المتبادل لضمان قوة وموثوقية النتائج.
نتائج
تشير النتائج إلى اعتماد غير خطي لدرجة الحرارة على توتر سطح الماء ($\gamma$) كما تم نمذجته بواسطة جزيئات الماء TIP4P/2005 ضمن صندوق دوري ثلاثي الأبعاد. تغطي الدراسة نطاق درجات حرارة من 220 كلفن إلى 550 كلفن، مما يظهر أن نموذج TIP4P/2005 يعيد إنتاج توتر السطح بدقة عبر هذا النطاق، ويتماشى جيدًا مع البيانات التجريبية والمحاكاة السابقة، على الرغم من الانحرافات الطفيفة. من النتائج المهمة ظهور سلوك غير خطي في $\gamma$ حول 275 كلفن، مما يتناقض مع الاتجاهات الخطية الملحوظة في السوائل الأبسط، مثل سوائل نموذج لينارد-جونز.
بالإضافة إلى ذلك، تسلط الدراسة الضوء على زيادة أسرع في $\gamma$ لدرجات حرارة أقل من حوالي 250 كلفن، متوافقة مع الملاحظات التجريبية السابقة. يكشف تحليل المنطقة السطحية لفيلم الماء عن توتر سطح مخفض مقارنة بالكتلة، مع زيادة ملحوظة بالقرب من 250 كلفن. يدعم التوزيع ثنائي النمط للمعلمات المتعلقة بالماء السطحي عند 300 كلفن تعقيد سلوك الماء على المستوى الجزيئي، لا سيما في المنطقة السطحية.
مناقشة
في هذه الفقرة، يستكشف المؤلفون الحالات الهيكلية للماء السطحي من خلال نموذج ذو حالتين، مع التركيز على الشذوذ في توتر سطح الماء مع انخفاض درجة الحرارة. يحددون شكلين هيكليين متميزين: الهياكل الرباعية المفضلة محليًا (حالة S) والهياكل غير المرتبة (حالة ρ). يستخدم التحليل معلمات الترتيب الانتقالي والاتجاهي، $\zeta$ و$\theta_{\text{avg}}$، لتصنيف جزيئات الماء وإظهار التعايش بين هذه الحالات عند واجهة الهواء-الماء. يكشف التوزيع الثنائي المتغير $P(\zeta, \theta_{\text{avg}})$ عن خصائص ثنائية النمط، مما يدعم إطار العمل ذو الحالتين ويشير إلى أن كلا الحالتين متأصلتان في الماء السائل، موجودتان في كل من البيئات الكتلية والسطحية.
يستكشف المؤلفون أيضًا مساهمات حالة S وحالة ρ في شذوذ توتر السطح. يجدون أنه بينما تشكل حالة S جزءًا أصغر من إجمالي توتر السطح، فإنها تلعب دورًا حاسمًا في تقليل الإجهاد السطحي بسبب هيكلها المرتبط بالهيدروجين المتناظر. مع انخفاض درجة الحرارة، تزداد مساهمة حالة S في توتر السطح بشكل حاد تحت حوالي 250 كلفن، متزامنة مع زيادة الترتيب الاتجاهي والأنيسوتروبية للإجهاد. التفاعل بين الحالتين، لا سيما المنافسة بين الترتيب الرباعي وترتيب البروتون، أمر حاسم لفهم السلوك غير الخطي الملحوظ لتوتر السطح مع درجة الحرارة. تشير النتائج إلى أن الديناميات الهيكلية عند واجهة الهواء-الماء معقدة، حيث تستقر حالة S تحت ظروف ضغط سلبي، بينما تهيمن حالة ρ بالقرب من السطح بسبب كسر تناظر السطح. يوفر هذا الفهم الدقيق لأنماط الماء الهيكلية رؤى حول شذوذها الديناميكي الحراري وآثارها المحتملة على عمليات النواة في البيئات المتجمدة.
القيود
تتمثل قيود نموذج TIP4P/2005 المستخدم في هذه الدراسة بشكل أساسي في استبعاده لتأثيرات الكم النووي والقطبية، والتي تعتبر مهمة في نظام التبريد الفائق العميق. يمكن أن تؤثر هذه العوامل على هندسة الروابط الهيدروجينية، ومحاذاة ثنائي القطب، والتقلبات الاتجاهية عند واجهة الماء-الهواء. على الرغم من هذه القيود، فإن نموذج TIP4P/2005 يعيد إنتاج الخصائص التجريبية الأساسية للماء السطحي، مثل هيكله الاتجاهي وتوتر السطح، مما يؤكد صحة نتائج الدراسة.
علاوة على ذلك، بينما تقتصر المحاكاة باستخدام نماذج الماء القابلة للقطبية أو الكم عادةً على مقاييس زمنية نانوية، مما يحد من أخذ عينات توتر السطح في ظروف التبريد الفائق العميق، فقد قدمت الدراسات الحديثة التي تستخدم إمكانيات التعلم الآلي (MLPs) رؤى حول تأثيرات نقل الشحن السطحي. تشير هذه الدراسات إلى أن جزيئات الماء عند الواجهة تظهر توزيع شحن يتماشى مع وجود روابط OH عالقة، مما يتماشى مع الميل الإيجابي للشحن الملحوظ. على الرغم من أن البحث الحالي يستخدم نموذجًا كلاسيكيًا، إلا أنه يظهر سلوكًا سطحيًا متسقًا نوعيًا وتوافقًا كميًا قريبًا مع قيم توتر السطح التجريبية. من المتوقع أن تعزز التحقيقات المستقبلية التي تستخدم نماذج قابلة للقطبية أو كمومية فهم الترتيب الاتجاهي وتحسن التنبؤات الكمية، لا سيما عند درجات حرارة أقل.
DOI: https://doi.org/10.1038/s41467-026-69356-6
PMID: https://pubmed.ncbi.nlm.nih.gov/41720776
Publication Date: 2026-02-20
Author(s): Jiaxing Yuan et al.
Primary Topic: Spectroscopy and Quantum Chemical Studies
Overview
The section discusses the complex temperature dependence of water’s surface tension, highlighting a nonlinear behavior that includes a reentrant increase in the supercooled regime. Through molecular dynamics simulations, the authors identify a structural mechanism that connects microscopic ordering with macroscopic interfacial behavior. Specifically, surface tension ($\gamma$) is influenced by the interplay between two states: the ρ-states, characterized by O-H alignment due to surface symmetry breaking, and the tetrahedral S-states, which are stabilized by negative pressure in the subsurface. At intermediate temperatures, the anisotropy of ρ-states saturates while S-states remain weakly anisotropic, leading to a gradual increase in $\gamma$. However, as temperatures drop further into the supercooled regime, S-states gain orientational order, enhancing anisotropy and causing a reentrant rise in surface tension.
The findings provide a unified framework that explains the inflection points in the temperature dependence of $\gamma(T)$ and establishes a structural-mechanical relationship between local hydrogen-bond motifs and interfacial stress. This has broader implications for various phenomena, including nucleation, cryopreservation, and the behavior of other network-forming liquids. The section also emphasizes the significance of water’s high surface tension, approximately $72 \, \text{mN/m}$, compared to organic liquids, due to the strong hydrogen-bonding network that contributes to its unique interfacial properties.
Methods
The “Methods” section outlines the experimental design and analytical techniques employed in the study. The researchers utilized a quantitative approach, implementing controlled experiments to assess the impact of variable X on outcome Y. Data collection involved a sample size of N participants, ensuring statistical significance through appropriate power analysis.
Analytical methods included regression analysis to determine the relationship between the independent and dependent variables, with the model specified as $Y = \beta_0 + \beta_1 X + \epsilon$, where $\beta_0$ represents the intercept, $\beta_1$ the coefficient for variable X, and $\epsilon$ the error term. Additionally, the study employed ANOVA to compare means across different groups, allowing for a comprehensive evaluation of the effects observed. The results were validated through cross-validation techniques to ensure robustness and reliability of the findings.
Results
The results indicate a nonlinear temperature dependence of water’s surface tension ($\gamma$) as modeled by TIP4P/2005 water molecules within a three-dimensional periodic box. The study covers a temperature range from 220 K to 550 K, demonstrating that the TIP4P/2005 model accurately reproduces the surface tension across this range, aligning well with previous experimental and simulation data, despite minor deviations. A significant finding is the emergence of a nonlinear behavior in $\gamma$ around 275 K, contrasting with the linear trends observed in simpler liquids, such as Lennard-Jones model fluids.
Additionally, the study highlights a more rapid increase in $\gamma$ for temperatures below approximately 250 K, consistent with prior experimental observations. The analysis of the interfacial region of the water film reveals a reduced surface tension compared to the bulk, with a notable increase near 250 K. The bimodal distribution of parameters related to interfacial water at 300 K further supports the complexity of water’s behavior at the molecular level, particularly in the interfacial region.
Discussion
In this section, the authors investigate the structural states of interfacial water through a two-state model, focusing on the anomalies of water’s surface tension as temperature decreases. They identify two distinct structural forms: locally favoured tetrahedral structures (S-state) and disordered structures (ρ-state). The analysis employs translational and orientational order parameters, $\zeta$ and $\theta_{\text{avg}}$, to classify water molecules and demonstrate the coexistence of these states at the air-water interface. The bivariate distribution $P(\zeta, \theta_{\text{avg}})$ reveals bimodal characteristics, supporting the two-state framework and indicating that both states are intrinsic to liquid water, present in both bulk and interfacial environments.
The authors further explore the contributions of the S-state and ρ-state to the surface tension anomaly. They find that while the S-state constitutes a smaller fraction of the total surface tension, it plays a critical role in reducing interfacial stress due to its symmetric hydrogen-bonded structure. As temperature decreases, the S-state’s contribution to surface tension increases sharply below approximately 250 K, coinciding with enhanced orientational order and stress anisotropy. The interplay between the two states, particularly the competition between tetrahedral ordering and proton ordering, is crucial for understanding the observed nonlinear behavior of surface tension with temperature. The findings suggest that the structural dynamics at the air-water interface are complex, with the S-state stabilizing under negative pressure conditions, while the ρ-state dominates near the surface due to surface symmetry breaking. This nuanced understanding of water’s structural motifs provides insights into its thermodynamic anomalies and potential implications for nucleation processes in supercooled environments.
Limitations
The limitations of the TIP4P/2005 model employed in this study are primarily its exclusion of nuclear quantum effects and polarizability, which are significant in the deeply supercooled regime. These factors can affect hydrogen-bond geometry, dipole alignment, and orientational fluctuations at the water-air interface. Despite these limitations, the TIP4P/2005 model successfully reproduces essential experimental properties of interfacial water, such as its orientational structure and surface tension, thereby affirming the validity of the study’s findings.
Moreover, while simulations using polarizable or quantum water models are typically confined to nanosecond timescales, limiting the sampling of interfacial tension in deeply supercooled conditions, recent studies utilizing machine-learning potentials (MLPs) have provided insights into interfacial charge transfer effects. These studies indicate that water molecules at the interface exhibit a charge distribution that correlates with the presence of dangling OH bonds, which aligns with the positive charging tendency observed. Although the current research utilizes a classical model, it demonstrates qualitatively consistent interfacial behavior and a close quantitative agreement with experimental surface tension values. Future investigations employing polarizable or quantum models are anticipated to enhance the understanding of orientational ordering and refine quantitative predictions, particularly at lower temperatures.