彎曲的冰如何帶電，從撓曲電量測追到表面相變與鹽水流動

一片冰受到彎曲，拉伸的一側與壓縮的一側可以產生電極化。這使原本看似缺乏壓電性的普通冰，成為機械作用與電訊號之間的轉換材料，也替雷雨雲裡冰粒碰撞帶電提供一個可量化的候選機制。Wen 等人的冰研究把這個問題從定性猜想推進到實驗量測。

雷雨雲持續把空氣、水滴與冰粒搬動。微小冰晶和較大的霰粒相撞，再被氣流與重力分開，局部接觸帶電逐漸累積成雲中的電荷分布。理解閃電，需要連起「碰撞時如何帶電」與「整朵雲如何維持電荷分離」兩個尺度。冰的撓曲電性處理前一段，雲微物理與氣流處理後一段。2025 年雲微物理研究顯示，冰粒的生成途徑也會改變模擬中的電化與閃電活動。

下面從圖的上半部看外力如何進入材料，再看下半部如何形成可量測的電訊號。重點在於：同樣「彎曲後帶電」的外觀，可以由固體的電極化，也可以由介面與液體傳輸共同造成。

水分子有極性，為何整塊冰仍缺乏壓電性

水分子的氧端與氫端具有不同電荷分布，形成電偶極。若大量偶極方向一致，材料會出現宏觀電極化；若方向彼此抵消，單個分子的極性仍然存在，整塊材料的平均極化卻可以為零。

日常大氣壓下常見的是六方冰 Ih。它的氧原子建立有規律的骨架，氫的位置則受「冰規則」約束而保有無序：每個氧附近有兩個較接近的氫，每條氧與氧之間的氫鍵只有一個氫。這些局部約束允許大量不同排列，分子極性無法自然全部沿同一方向對齊。冰的多種結構與質子有序問題，可參考 Bartels-Rausch 等人的冰物理綜述。

壓電效應描述材料受均勻應變時，能否出現與應變成正比的極化。對中心對稱晶體而言，把空間方向反轉，極化應反向，均勻應變卻保持原樣；同一材料若同時滿足這兩種描述，線性的壓電係數只能為零。這是普通冰缺乏宏觀壓電性的對稱性理由。分子是否有偶極，只回答微觀電荷分布；材料是否壓電，還要看整體結構。撓曲電理論綜述

彎曲多出另一個條件：上、下各層的應變不同。這種空間差異會在原本對稱的材料裡建立方向，使它具備另一條機械—電耦合途徑。

一個梯度讓對稱的材料也能極化

想像彎曲一把尺。外弧被拉長，內弧被壓短，中間有一層長度幾乎維持原樣。應變 ε 表示局部伸長量除以原長度；應變梯度則表示走過一段距離後，應變變了多少。

撓曲電效應的最小模型是：

P = μ × ∂ε/∂z

P 是極化，單位 C/m²；μ 是撓曲電係數，單位 C/m；z 是厚度方向。應變沒有單位，因此梯度的單位為 1/m。兩者相乘，剛好得到電荷面密度。方向與晶體取向的完整描述使用張量，例如 Pᵢ = μᵢⱼₖₗ ∂εⱼₖ/∂xₗ；這裡先固定量測方向，保留讀懂實驗所需的關係。一般理論、晶體量測

對薄梁、小變形與近似純彎曲而言，ε 隨厚度位置近似線性改變，∂ε/∂z 的量級就是曲率 κ = 1/R，R 為彎曲半徑。從這裡可得到一個重要設計判斷：若比較條件是相同表面應變，較薄的材料允許較大的曲率，撓曲電訊號也可能增強；若比較條件是相同曲率，單靠變薄就沒有同樣的增益。談尺度效應，必須先固定比較條件。

現象	主要驅動	觀察到的反應	適合回答的問題
壓電性	均勻應變	誘發極化	壓縮或拉伸能否直接轉成電訊號
撓曲電性	應變的空間梯度	誘發極化	彎曲、裂尖或接觸變形如何帶電
鐵電性	自發極化及電場切換	可切換的極化狀態	材料是否具有可保留的極性排列
流動電現象	液體在帶電介面附近流動	輸運電荷形成電流	多孔固體中的液體如何參與電訊號

前三者的關係見撓曲電綜述；最後一項在鹽冰研究中成為放大反應的來源。

原論文建立了哪些新知識

Xin Wen（ICN2／西安交通大學）等人的〈Flexoelectricity and surface ferroelectricity of water ice〉於 2025 年 8 月 27 日正式發表於 Nature Physics，第 21 卷，1587–1593 頁。2022 年已有預印本；正式發表與預印本上傳是同一研究的不同階段。正式論文、公開手稿

其實驗以超純水製成夾在電極間的多晶冰，使用動態三點彎曲，同步記錄位移與電荷，再換算成應變梯度與極化。極化對梯度近似線性，斜率給出有效撓曲電係數。公開手稿在較穩定的 203–248 K 區間報告平均 1.14 ± 0.13 nC/m；這是特定橫向量測幾何的有效係數。公開手稿，Fig. 1–2

研究另以金、鉑及鋁電極對照低溫反應，並觀察預極化電場造成的遲滯。結合介面計算，結果支持約 160 K（−113°C）的近表面鐵電相變。第三部分把已量測係數代入冰晶—霰粒碰撞模型，所得電荷量級與先前接觸帶電實驗相容。正式論文，Fig. 3–5

這三個成果各自處理不同層次：彎曲量測確立電機械反應；介面研究解釋低溫異常；碰撞模型連接自然現象。讀原文時，沿著這三段推進，比把每張圖都當成同一種「發電證明」更容易掌握論點。

從彎曲位移到極化係數，量測電路做了什麼

電荷放大器讀到的是電極上的電荷變化 Q，位移感測器讀到的是梁的彎曲。要把這兩個不同單位的訊號接起來，先用機械幾何把位移換成梯度，再以有效電極面積 A 把電荷換成 P ≈ Q/A。對照多個彎曲幅度，才能判斷是否具有穩定的比例關係。

下面用假設算例建立訊號大小的直覺。取 μ = 1 nC/m、κ = 10 m⁻¹，以及 A = 1 cm²。在上述簡化梁模型下：

• P ≈ μκ = 10⁻⁸ C/m²。
• Q ≈ PA = 10⁻¹² C，也就是 1 pC。
• 若電荷正弦振盪，頻率為 1 Hz，電流振幅 I₀ ≈ 2πfQ₀，約 6.3 pA。

這是把指定假設代入模型的試算，並非原論文的裝置輸出。它解釋了為何這類研究需要低雜訊電荷量測、穩定接線，以及同步辨認機械與電訊號。

「量到交流電流」還只是開始。電纜移動、接觸摩擦、漏電或溫度變動，也可能留下訊號。一般實驗設計可藉由反轉彎曲、改變幅度、掃描頻率與電極，檢查各候選機制的預測。這些是讀量測結果的診斷工具；不同研究完成的對照項目仍以各自方法為準。2024 年量測與材料綜述

例如，理想線性介電反應在一定頻帶內具有近似固定的 P/κ；液體流動、離子擴散或電荷陷阱則有時間尺度，改變頻率會改變振幅與相位。把頻率當成區分機制的工具，比單看一次最大的電壓更有解釋力。

薄薄的表面，為何能影響厚厚的冰

材料內部的水分子四周都有鄰居，靠近金屬的分子卻處在不同化學與電學環境。介面可以影響偶極方向、電荷補償與自由能，讓局部極性排列比體材料更有利。Stony Brook 的研究說明介紹了金屬—冰介面的量子計算如何與低溫實驗互相補充。

鐵電相變的直覺，是材料由平均缺乏有序極性的狀態，轉成具有可切換極性排列的狀態。電極若能補償排列造成的表面電荷，就可能降低極性狀態的能量代價。不同金屬與氧化表面的電學條件不同，因此改換電極同時也在改變被測系統。

原研究的金電極試片在約 165 K 出現約 7.6 nC/m 的峰值；鉑電極的峰值約 15 nC/m，出現在約 159 K；帶有原生氧化層的鋁電極則沒有相同峰值。研究者再改變預極化電場，關閉偏壓後才量測撓曲電係數，得到蝶形遲滯。電極依賴性與可保留的切換反應一起支持介面鐵電解釋。公開手稿，Fig. 3

計算端使用密度泛函理論，比較金屬介面附近的質子有序冰 XI 與無序冰 Ih，再把介面能量放入自由能模型。有序排列在金、鉑介面得到能量上的穩定化，使模型的相變溫度能與實驗對照。這把電極、分子排列與低溫峰值接成一條可檢驗的解釋鏈。公開手稿，Fig. 4

為何表面層很薄，訊號仍可能明顯？考慮一個簡化的雙表面模型。上、下表面初始極性方向相反，平直時貢獻可以抵消；彎曲時一側受拉、另一側受壓，各自的極化變化反而可能相加。表面層占材料體積雖小，它對電極訊號的貢獻還受應變與幾何放大影響。這是有限尺寸材料量測中的一般問題。表面與體材料耦合綜述

另一個需要分開理解的現象是預融。冰表面在整塊材料熔化前，就可能出現較無序、具有液體特性的準液態層。其存在範圍與性質受溫度、晶面、雜質及環境影響。這會改變表面的機械滑動與電荷移動能力。Slater 與 Michaelides，2019

因此，掃描冰的溫度同時在探索固體極化、介面排序與液體傳輸三種不同物理。約 160 K 的近表面鐵電現象，和接近融點的液體活動，應放在各自的溫度與結構條件下討論。

在這批試片中，約 248 K（−25°C）以上的電訊號相位開始明顯改變，機械反應也逐漸軟化。研究將兩者連到預融層中的可移動離子與晶界滑動。溫度掃描的價值正在這裡：它讓同一個電極訊號中的不同來源，透過相位與力學變化逐步顯現。公開手稿，Fig. 2

鹽冰把係數提高千倍，增加的是一條傳輸路徑

2025 年 9 月 15 日，Xin Wen（ICN2／西安交通大學）等人在 Nature Materials 發表〈Streaming flexoelectricity in saline ice〉。研究報告 NaCl 摻入冰後，有效撓曲電係數由純冰常見的約 1–10 nC/m 提高至 1–10 μC/m，約三個數量級；設計的彎曲裝置具有約 4,000 pC/N 的有效壓電係數。原論文

這條放大途徑發生在晶界。冰晶之間保留含鹽的液體通道，彎曲造成壓縮區與拉伸區的差異，驅動液體沿通道移動。液體攜帶介面附近可移動的電荷，形成流動電流，再被外部電路讀取。研究機構也指出，循環後反應衰減與實際功率仍是應用課題。ICN2 說明

這讓「有效材料係數」有了更具體的含義：外部看到的極化或電荷，可能同時包含固體的瞬時反應與內部通道的電荷輸運。摻鹽改變了微結構、液體與離子的活動；只看冰晶格本身，就會漏掉增強的主因。

讀數字時也要保留單位。μ 的 C/m 表示極化與梯度之比，d 的 C/N 表示裝置電荷與力的關係；從前者得到後者，需要指定裝置形狀、彎曲方式與力學條件。若再問瓦特，還必須加入頻率、負載、漏電與機械輸入功率。係數、電荷與功率分別處在材料、裝置與能量收支三個層次。

冰粒碰撞如何走到雲中的電荷分離

接觸時的變形集中在很小範圍，應變梯度可以比整片材料的緩慢彎曲更大。這種局部放大在裂尖研究中已受到重視：Wang（Max Planck Institute for Solid State Research）等人於 2019 年線上發表、2020 年卷期刊出的氧化物研究，直接觀察裂尖附近強烈的撓曲電極化。Nano Letters 原論文

若碰撞中的極化吸引可移動電荷到接觸區，顆粒分開時又有一部分電荷留在不同顆粒上，就得到可持續的帶電結果。接著，上升氣流與沉降速度不同，將小冰晶與較大的霰搬到不同位置，局部帶電才逐步變成大尺度電荷分離。

這裡可分三段思考：第一段是機械變形產生極化；第二段是可移動電荷重新分布並在分離後保留；第三段是雲中輸運累積電場。每一段都有自己的時間尺度。變形很快，離子來不及移動，結果會不同；碰撞後電荷迅速回流，長期累積也會改變。這是從局部材料反應推到氣象現象時，自然浮現的模型需求。

同樣地，2025 年的雲解析模擬研究比較不同冰生成途徑，發現次生冰過程會改變粒子數量與電化活動。它處理的是粒子族群和整朵雲的演化；與冰的局部量測放在一起，才能看見從接觸機制到雲尺度的研究拼圖。Atmospheric Chemistry and Physics

從早期陶瓷到 2026 年，研究脈絡如何展開

近代進展可整理成三條路：建立可計算的材料耦合、辨認表面與可移動電荷的貢獻，以及把梯度設計成可用的功能。下表挑選改變理解的代表研究，區分直接研究冰與其他材料的平行進展。

年份	代表研究	推進的理解	與冰研究的關係
2007	SrTiO₃ 晶體的應變梯度極化	在晶體中量測撓曲電反應並估計張量分量	建立材料量測參考
2012	機械寫入鐵電極化	局部機械作用可控制薄膜極化狀態	展示梯度可以操控極性
2013	固體撓曲電綜述	整理對稱性、尺度效應與表面貢獻	提供概念骨架
2016	氧化物半導體的增強類撓曲電反應	可移動電荷參與量測中的放大反應	提醒有效係數包含不同來源
2018	Flexo-photovoltaic effect	應變梯度打開原本受對稱性限制的光電途徑	跨領域應用，非冰的實驗
2019–2020	裂尖極化直接觀察	奈米尺度的局部梯度可造成強烈極化	支持研究接觸與破裂的動機
2025	冰的撓曲電性與近表面鐵電性	普通冰的量測、介面相變及碰撞模型	本文主論文
2025	鹽冰的流動撓曲電性	晶界液體與離子傳輸提供強烈有效反應	直接延伸至含鹽冰
2026	冰奈米電容的光學介電常數計算	改善電極與介電材料電荷的區分	同領域量測解讀，非撓曲電複驗
2026	梯度驅動單晶自極化	材料不均勻性可驅動熱循環後的極化排列	平行的功能材料進展
2026	結構化木材的高撓曲電性	化學處理與階層結構共同增強反應	平行的天然材料進展

2012 年 Lu（University of Nebraska–Lincoln）等人的工作，用探針局部施力控制鐵電薄膜的極化。它顯示機械梯度也可作為「寫入」手段，將撓曲電性由被動量測推向主動控制。論文與公開研究典藏

2018 年 Yang（University of Warwick）等人的光電研究，則利用梯度打破對稱性，拓展體光伏效應能出現的材料條件。這裡的能量輸入包含光，與外力驅動的冰裝置有不同收支；相通之處在於利用空間不均勻性創造新的電反應。論文與公開手稿

2026 年 1 月 16 日的冰奈米電容研究特別適合補充介面判讀。Mannino（Stony Brook University）等人指出，極薄電容裡如何區分電極電荷與冰的極化電荷，會影響推得的介電常數；他們提出更明確的電荷分離計算方法。研究處理的是光學頻率所對應的電子反應，讀者應與低頻分子取向反應分開。PRL 原論文、公開手稿

2026 年 4 月 9 日，Choi（UNIST）等人研究 Mn 摻雜的 PMN–PT 單晶，連結晶格與熱性質的不均勻分布、撓曲電驅動及單向自極化。外部電路的電荷補償也影響結果。這讓梯度設計與熱循環後的功能保留開始接在一起。npj Flexible Electronics

4 月 22 日，Gao（Lanzhou University）等人以去木質素與壓縮處理木材，報告 36.72 nC/m 的撓曲電係數，並展示動作感測。材料的分子結構與階層幾何一起參與反應放大。Nature Communications

這些進展共同提示：選材料時，可以同時設計化學組成、內部微結構、介面與受力方式。對冰而言，水的廉價與可就地成形只提供起點；鹽分、晶界和電極，才決定電反應在指定環境下如何運作。

發電、感測與材料研究需要不同的判斷尺度

感測器的任務是把事件轉成可辨認訊號。輸出能量很小仍可能有用，前提是訊號高於雜訊、校準可維持，且電路能讀到它。能量擷取裝置則要把電荷送進負載或儲能元件，必須計算完整週期的可交付能量。

取一個簡化的電荷源與電容模型，開路電壓 V ≈ Q/C，單次充電所存能量 U ≈ Q²/(2C)。把先前假設的 Q = 1 pC，配上 C = 10 pF，可得到 V ≈ 0.1 V、U ≈ 5 × 10⁻¹⁴ J。若每秒只有一次同等事件，這個簡化能量尺度約為 50 fW；實際週期擷取還取決於波形與負載。

這個假設算例揭示電壓與能量的差別：容易看到的電壓，可以對應很小的電荷與能量。增大面積、提高曲率、設計更強反應或增加事件頻率，會同時改變機械強度、電容、電荷損耗與壽命。能量擷取的比較應指定機械輸入條件，再量負載中的電壓與電流。

就地結冰製作寒地感測器，與主動冷卻製作發電設備，也有不同成本。前者可以利用既有低溫環境；後者要把製冷、維持冰態及機械驅動的能耗一起納入。應用判斷應回到使用情境，再選擇電壓、雜訊、能量或壽命作為主要指標。

由此可提出一個可檢驗的機制判斷：若鹽冰的放大主要由液體通道輸運驅動，改變溫度、頻率或晶界連通性，應系統性改變增益與相位。若通道活動受到抑制後，巨大增益仍維持且呈現純介電式反應，就需要重新檢查目前的機制分配。這是從傳輸模型推得的實驗方向。

這種跨層思考也可對照本站的運算架構與數值預算文章：材料係數與器件峰值之後，還有系統如何使用它的問題。冰研究在此留下的長期理解是，材料的電反應由固體、介面與內部傳輸共同決定；可利用的功能則由裝置幾何、時間尺度與外部電路接手決定。

證據範圍

原論文確立特定多晶冰試片、電極及彎曲條件下的撓曲電反應。1.14 ± 0.13 nC/m 與 203–248 K 的細節引自公開手稿；完整張量、不同冰相、天然冰與所有受力條件不能用此單一有效係數概括。圖中的路徑、符號與尺寸為說明用途，沒有重製原始量測曲線。

約 160 K 的結果支持金屬—冰介面附近的鐵電相變，不等於所有自由冰表面、整塊冰或一般冰箱溫度都具有同樣鐵電性。表面層的微觀解釋結合模型與計算；它不是將整塊冰直接量成均勻鐵電體。

碰撞模型提供與已知帶電實驗相容的量級及電荷上限，沒有完整解出分離時電荷如何保留。捕獲、回流、塑性、破裂、液體與物質轉移仍會影響自然結果。因此，論文不能被概括為「已解開閃電的全部成因」。低溫介面相變也不應直接套到溫度高得多的雷雨雲。

鹽冰的約千倍增強是特定條件下的有效係數變化，不等於純冰晶格本徵係數提高千倍，更不等於發電功率提高千倍。液體與離子傳輸帶來頻率、溫度、微結構與壽命依賴。文中電荷、電流與能量試算均為指定假設下的教學模型，並非論文裝置的重現數據。

2026 年列出的其他材料研究及冰奈米電容計算，提供平行機制與相關量測知識，不能當成獨立複驗 2025 年全部冰結論。研究脈絡以截至 2026 年 10 月 3 日能核實的代表性工作為範圍，並非所有論文的系統性文獻回顧。

截圖的敘述式標題與「No Simulation」不是正式論文題名或方法的準確摘要。原研究包含直接量測、介面量子計算及碰撞模型，三者的角色應分開閱讀。

References

1. Wen, X. et al. Flexoelectricity and surface ferroelectricity of water ice. Nature Physics 21, 1587–1593 (2025). 公開手稿。
2. Wen, X. et al. Streaming flexoelectricity in saline ice. Nature Materials 24, 1533–1537 (2025)。
3. Zubko, P., Catalan, G. & Tagantsev, A. K. Flexoelectric Effect in Solids. Annual Review of Materials Research 43, 387–421 (2013)。
4. Zubko, P. et al. Strain-Gradient-Induced Polarization in Single Crystals. Physical Review Letters 99, 167601 (2007)。
5. Bartels-Rausch, T. et al. Ice structures, patterns, and processes. Reviews of Modern Physics 84, 885–944 (2012)。
6. Lu, H. et al. Mechanical Writing of Ferroelectric Polarization. Science 336, 59–61 (2012)。
7. Narvaez, J., Vasquez-Sancho, F. & Catalan, G. Enhanced flexoelectric-like response in oxide semiconductors. Nature 538, 219–221 (2016)。
8. Yang, M.-M., Kim, D. J. & Alexe, M. Flexo-photovoltaic effect. Science 360, 904–907 (2018)。
9. Wang, H. et al. Direct Observation of Huge Flexoelectric Polarization around Crack Tips. Nano Letters 20, 88–94 (2020)。
10. Slater, B. & Michaelides, A. Surface premelting of water ice. Nature Reviews Chemistry 3, 172–188 (2019)。
11. Kołodziej, M. et al. Fundamentals of Flexoelectricity, Materials and Emerging Opportunities Toward Strain-Driven Nanocatalysts. Small 20, 2406726 (2024)。
12. Distinct effects of several ice production processes on thunderstorm electrification and lightning activity. Atmospheric Chemistry and Physics 25 (2025)。
13. Mannino, A. et al. First-Principles Nanocapacitor Simulations of the Optical Dielectric Constant in Water Ice. Physical Review Letters 136, 026202 (2026)。
14. Choi, W.-J. et al. Flexoelectricity enables piezoelectric single crystals to be self-poled. npj Flexible Electronics 10, 71 (2026)。
15. Gao, Y. et al. Enhanced strain gradient in structural wood for high flexoelectricity. Nature Communications 17, 5599 (2026)。