孫文祥，李志文，劉林，梁偉華，王英龍

(1.河北民族師范學(xué)院 物理系，河北 承德 067000;2.河北民族師范學(xué)院 管理工程系，河北 承德 067000;

3.北華航天工業(yè)學(xué)院 經(jīng)濟(jì)管理系，河北 廊坊 065000;4.河北大學(xué) 物理科學(xué)與技術(shù)學(xué)院，河北 保定 071002)

應(yīng)力對(duì)熱電體電熱效應(yīng)的影響

孫文祥1，李志文2，劉林3,4，梁偉華4，王英龍4

(1.河北民族師范學(xué)院 物理系，河北 承德 067000;2.河北民族師范學(xué)院 管理工程系，河北 承德 067000;

3.北華航天工業(yè)學(xué)院 經(jīng)濟(jì)管理系，河北 廊坊 065000;4.河北大學(xué) 物理科學(xué)與技術(shù)學(xué)院，河北 保定 071002)

在一定溫度范圍內(nèi)，采用Landau-Devonshire自由能理論，討論了應(yīng)變對(duì)電場(chǎng)導(dǎo)致的PbZr0.4Ti0.6O3(PZT)以及[P(VDF-TrFE)](65/35)偏氟乙烯和三氟乙烯共聚物薄膜和塊材的系統(tǒng)熵的影響，以及熵極大時(shí)的溫度變化.結(jié)果表明：隨著應(yīng)變的增加，聚合物和PZT的場(chǎng)致溫變和等溫熵變最大值時(shí)的溫度分別線性增大和減小.為了增大薄膜在較低工作溫度下的熵變，要求聚合物的應(yīng)變大于0，PZT的應(yīng)變小于0.

熱電體；電熱效應(yīng)；應(yīng)力；PZT

在2006年以前，Thacher[1]對(duì)鐵電和反鐵電Pb(Zr, Ti)O3的實(shí)驗(yàn)研究發(fā)現(xiàn)由于熱電體的電熱效應(yīng)引起的溫度變化為0.15～0.3 K，Annaorazov[2]等人的實(shí)驗(yàn)發(fā)現(xiàn)，電熱效應(yīng)引起的Fe49Rh51材料鐵磁反鐵磁相變溫度的變化在外加應(yīng)變時(shí)可以到達(dá)5 K.Xiao[3]等人研究(1-x)Pb(Mg1/3Nb2/3)O3-xPbTiO3材料，發(fā)現(xiàn)其溫度變化最大只有1 K，可見(jiàn)相變溫度的變化較小，因此在實(shí)際應(yīng)用中沒(méi)有引起人們的足夠重視[4-8].直到Mischenko等人[9]在350 nm厚的PbZr0.95Ti0.05O3薄膜在外加高電場(chǎng)(90 MV/m)作用下出現(xiàn)‘巨電熱效應(yīng)’，226 ℃時(shí)絕熱溫度改變達(dá)到12 K之后，才使熱電體的電熱效應(yīng)實(shí)際應(yīng)用于汽車、航天航空、食品等行業(yè)成為可能.熱電體還可應(yīng)用在超微規(guī)模上的芯片制冷，從而解決集成電路的高溫問(wèn)題.并且傳統(tǒng)蒸汽壓縮制冷系統(tǒng)對(duì)環(huán)境有害，而電熱效應(yīng)制冷是環(huán)保的固態(tài)制冷技術(shù).最近，Neese等人[10]在聚合物(偏二氟乙烯-三氟乙烯)[P(VDF-TrFE)](55/45)(物質(zhì)的量比)也獲得了大的絕熱溫度變化，其物理成因尚不明確.鐵電張弛振蕩器以其電熱現(xiàn)象顯著、相變溫度范圍廣是制冷設(shè)備的理想材料，但這尚未受到足夠的關(guān)注.

分別對(duì)典型的PZT鐵電材料以及偏氟乙烯和三氟乙烯共聚物[P(VDF-TrFE)](物質(zhì)的量比為65/35)，在特定電場(chǎng)變化下對(duì)電場(chǎng)誘發(fā)的薄膜熵變進(jìn)行了分析和討論，考察薄膜應(yīng)變對(duì)熱電體的電熱效應(yīng)的影響，以求獲得應(yīng)變對(duì)場(chǎng)致熵變的影響規(guī)律，能夠?qū)﹄姛嵝?yīng)的實(shí)際應(yīng)用提供一定的參考價(jià)值.

1 理論基礎(chǔ)

對(duì)于無(wú)約束、剛性自由的鐵電體塊材樣品，其Landau-Devonshire理論給出的吉布斯自由能體密度可以表示為[11-14]

(1)

式中α,β和γ是電介質(zhì)剛度系數(shù)，E,P和G0分別表示電場(chǎng)強(qiáng)度、極化強(qiáng)度和無(wú)外加電場(chǎng)情況下順電相自由能體密度,C是Curie常數(shù),ε0為真空介電常數(shù),T0是居里-外斯溫度.

對(duì)沉積在襯底上的薄膜材料，由于二者晶格常數(shù)的不同引起晶格失配以及熱力學(xué)膨脹，存在著內(nèi)應(yīng)力，自由能密度表達(dá)式中的電介質(zhì)剛度系數(shù)需要進(jìn)行修正.考慮力學(xué)邊界條件和面束縛，薄膜的自由能體密度改寫(xiě)為[15]

(2)

其中α1和β1是修正的電介質(zhì)剛度系數(shù)，其數(shù)學(xué)表達(dá)式見(jiàn)參考文獻(xiàn)[15].

假設(shè)薄膜和襯底的晶格常數(shù)分別為a0,as，由于襯底和薄膜的晶格常數(shù)不同,就會(huì)存在晶格失配，從而引起應(yīng)變.定義u=(as-a0)/as為由于晶格失配所引起的應(yīng)變，若as>a0，則應(yīng)變u>0，薄膜受到拉應(yīng)力；否則，若as

鐵電材料在穩(wěn)態(tài)下的自發(fā)極化強(qiáng)度P，滿足?Gbulk/?P=0或?Gfilm/?P=0.可以得到

E=αP+βP3+γP5,

(3)

E=α1P+β1P3+γP5.

(4)

由方程(3)、(4)可以得到不同外加電場(chǎng)的P-T曲線及不同溫度下的P-E曲線(電滯回線)；替換式(1)或(2)中的P(T,E)可以得到自由能體密度曲線方程.依據(jù)S=-(?G/?P)E，可以得到單位體積的熵SEX.假如外加電場(chǎng)從E1變到E2時(shí)，電場(chǎng)誘發(fā)的單位體積等溫熵變通過(guò)下式確定：

ΔS(T,ΔE)=SEX(T,E2)-SEX(T,E1).

(5)

場(chǎng)致熵變使得鐵電體系統(tǒng)的溫度變化ΔT，絕熱溫度變化用下式計(jì)算：

(6)

式中|ΔSm|為等溫熵變最大值，Tmax是|ΔSm|出現(xiàn)時(shí)的溫度，C1是材料的比熱容，ρ是其密度.

2 結(jié)果與分析

本模擬中使用的材料參數(shù)匯于表1中.

表1 鐵電薄膜的參數(shù)[16-17]

對(duì)于鐵電體材料偏氟乙烯和三氟乙烯共聚物[P(VDF-TrFE)](物質(zhì)的量比65/35)放置在任意的襯底上，使得鐵電薄膜與襯底之間出現(xiàn)任意可能的晶格失配應(yīng)變，取應(yīng)變u為連續(xù)變化，以覆蓋薄膜與任意材料之間的可能的應(yīng)變.假定晶體熱釋電軸沿著薄膜厚度方向，外加電場(chǎng)E、極化強(qiáng)度P與此方向平行，薄膜的下表面固定在基片上而上表面除了局部加緊，其他是自由的.而下端和襯底之間存在相互作用，屬于部分夾持狀態(tài)[15].

選擇外加電場(chǎng)E1=300 MV/m與E2=400 MV/m，在溫度為450～650 K時(shí)對(duì)外加電場(chǎng)導(dǎo)致的附加熵變進(jìn)行模擬，如圖1所示.

圖1 聚合物在不同應(yīng)變下熵變隨工作溫度的變化

從圖2中可以明顯看到，無(wú)論是塊材還是附著在某一襯底上的薄膜，在圖示的整個(gè)溫度范圍內(nèi)，熵變都是先隨溫度的增加而增大，在某一特定溫度(記為Tmax)時(shí)，達(dá)到極大值(記為-ΔSm)，之后又隨溫度的升高而降低，存在著熵變極值.

薄膜的熵變最大值(為18.21 J·kg-1·K-1)與應(yīng)變無(wú)關(guān)，基本不隨u值的變化發(fā)生變化.而熵變最大值對(duì)應(yīng)的特定溫度Tmax卻與應(yīng)變存在關(guān)聯(lián)，隨著應(yīng)變的增加而逐漸升高.從圖2可以看出，場(chǎng)致熵變導(dǎo)致的溫度變化從7.35 K到9.58 K均勻增長(zhǎng)，隨應(yīng)變線性變化；與熵變極值對(duì)應(yīng)的Tmax隨應(yīng)變線性增加.

圖3是聚合物以間隔溫變?yōu)?0 K在工作溫度從450 K到650 K時(shí)場(chǎng)致熵變的變化.在一個(gè)確定的工作溫度下，熵變隨著應(yīng)變的增加呈現(xiàn)拋物線形式的變化，先是逐漸升高在某個(gè)應(yīng)變值的地方達(dá)到極大獲得最大熵變，電熱效應(yīng)最明顯，再隨著應(yīng)變的增加熵變轉(zhuǎn)而降低，電熱效應(yīng)變?nèi)?這表明鐵電晶體薄膜在一定的工作溫度下工作時(shí)，只有合適的應(yīng)變才能使電熱效應(yīng)明顯，其他的應(yīng)變所獲得的溫度改變都是較小的，沒(méi)能充分發(fā)揮出電熱效應(yīng)的最大效用.與-ΔSm對(duì)應(yīng)的應(yīng)變?yōu)?0.056,-0.030,-0.002,0.026,0.053.

圖2 聚合物熵變極值溫度和場(chǎng)致溫變與應(yīng)變的關(guān)系

圖3 不同工作溫度下聚合物熵變與應(yīng)變的關(guān)系

同樣對(duì)PZT進(jìn)行對(duì)應(yīng)模擬，選擇外加電場(chǎng)E1=15 MV/m與E2=53 MV/m，對(duì)外加電場(chǎng)導(dǎo)致的附加熵變與應(yīng)變的關(guān)系進(jìn)行模擬，得到圖4.從圖4看出，最大熵變出現(xiàn)時(shí)對(duì)應(yīng)的應(yīng)變值為0.003 2,0.001 4,-0.000 4,-0.002 2,-0.004 0，呈線性減小.圖5的虛線表明了這一變化趨勢(shì).

圖4 不同工作溫度下PZT熵變與應(yīng)變的關(guān)系

圖5 PZT熵變極值溫度和場(chǎng)致溫變與應(yīng)變的關(guān)系

比較圖3和圖4中的曲線，發(fā)現(xiàn)2種材料對(duì)應(yīng)變的響應(yīng)是不同的.圖3顯示u<0時(shí)的壓應(yīng)力使得薄膜在較低的工作溫度，就獲得了比較明顯的熵變，電熱效應(yīng)效果較好.圖4卻表明，對(duì)于PZT來(lái)說(shuō)，u>0時(shí)的拉應(yīng)力才能使薄膜材料工作在較低的溫度，獲得較好的效果.這是因?yàn)椴煌N類的材料對(duì)應(yīng)變的響應(yīng)不同，拉應(yīng)力時(shí)應(yīng)選擇鐵電體材料來(lái)獲得較好的電熱效應(yīng).

將薄膜附著在襯底上出現(xiàn)晶格失配應(yīng)變是必然的，在應(yīng)用熱電體的電熱性質(zhì)時(shí)就要考慮相應(yīng)應(yīng)變和電場(chǎng)的共同作用下，薄膜熵變極大值所對(duì)應(yīng)的溫度Tmax，來(lái)確定晶體薄膜的工作溫度，以及熵變最大所引起薄膜溫度的變化也將是最大的，可以獲得比較好的實(shí)用效果.

3 結(jié)論

本文討論了偏氟乙烯和三氟乙烯共聚物以及PZT在一定的溫度范圍內(nèi)連續(xù)應(yīng)變對(duì)電場(chǎng)導(dǎo)致的系統(tǒng)熵的影響，系統(tǒng)溫度變化趨勢(shì)，熵極大時(shí)對(duì)應(yīng)的溫度變化.結(jié)果表明，隨著應(yīng)變的增加，聚合物的場(chǎng)致溫變和Tmax均線性增大，PZT的場(chǎng)致溫變和Tmax均線性減小.為了使得薄膜在較低的工作溫度下熵變比較明顯，要求聚合物的u<0，PZT的u>0.聚合物與PZT具備成為電制冷工質(zhì)材料的潛力,合理的替代或者摻雜改進(jìn)其工作溫度與電熱效應(yīng),將對(duì)電制冷材料的發(fā)展有重要作用.

[1]THACHER P D. Electrocaloric effects in some ferroelectric and antiferroelectric Pb(Zr, Ti)O3compounds[J]. J Appl Phys, 1968, 39:1996-2002.

[2]ANNAORAZOV M P, NIKITIN S A, TYURIN A L, et al. Anomalously high entropy change in FeRh alloy[J]. J Appl Phys, 1996, 79:1689-1695.

[3]XIAO D Q, WANG Y C, ZHANG R L, et al. Electrocaloric properties of (1-x)Pb(Mg1/3Nb2/3)O3-xPbTiO3ferroelectric ceramics near room temperature[J]. Mat Chem and Phys, 1998, 57:182-185.

[4]LAWLESS W N. Specific heat and electrocaloric properties of KTaO3at low temperatures[J]. Phys Rev B, 1977, 16:433-439.

[5]OLSEN R B, BUTLER W F, PAYNE D A, et al. Observation of a polarocaloric (electrocaloric) effect of 2 ℃ in lead zirconate modified with Sn4+and Ti4+[J]. Phys Rev Lett, 1980, 45:1436-1438.

[6]RADEBAUGH R. Behavior of the pyroelectric coefficient at low temperatures[J]. Phys Rev Lett, 1978, 40:572-574.

[7]LOMBARDO G, POHL R. O. Electrocaloric effect and a new type of impurity mode[J]. Phys Rev Lett, 1965, 15:291-293.

[8]TUTTLEA B A, PAYNEA D A. Ferroelectrics[M]. London: Gordon and Breach Science Publisher, 1981:603-606.

[9]MISCHENKO A S, ZHANG Q, SCOTT J F, et al. Giant Electrocaloric Effect in Thin-Film PbZr0.95Ti0.05O3[J]. Science, 2006, 311:1270-1271.

[10]NEESE B, CHU B, LU S G, et al. Large electrocaloric effect in ferroelectric polymers near room temperature[J]. Science, 2008, 321:821-823.

[11]ALPAY S P, MISIRLIOGLU I B, SHARMA A, et al. Structural characteristics of ferroelectric phase transformations in single-domain epitaxial films[J]. J Appl Phys, 2004, 95:8118-8123.

[12]WANG Xingyuan, WANG Yinglong, YANG Rongjia. Lattice model for strained nanoscale ferroelectric capacitors: Investigation on fundamental size limits in ferroelectricity[J]. Appl Phys Lett, 2009, 95:142910-1-142910-3.

[13]WANG B, XIA R. Dynamic process of domain switching in ferroelectric films[J]. J Appl Phys, 2003, 94:3384-3389.

[14]YANG Wei, CHEN Longqing. Computer simulation of the dynamics of 180° ferroelectric domains[J]. J Am Ceram Soc, 1995, 78:2554-2556.

[15]LIU P F, MENG X J, CHU J H, et al. Effect of in-plane misfit strains on dielectric and pyroelectric response of poly(vinylidene fluoride-trifluoroethylene) ferroelectric polymer[J]. J Appl Phys, 2009, 105:114105-1-114105-4.

[16]CHEN L, NAGARAJAN V, RAMESH R, et al. Nonlinear electric field dependence of piezoresponse in epitaxial ferroelectric lead zirconate titanate thin films[J]. J Appl Phys, 2003, 94:5147-5152.

[17]PERTSEV N A, KUKHAR V G, KOHLSTEDT H, et al. Phase diagrams and physical properties of single domain epitaxial Pb(Zr1-xTix) O3thin films[J]. Phys Rev B, 2003, 67:054107-1-054107-10.

Effectsofstressesonelectrocaloriceffectsofthethermoelectricbody

SUNWenxiang1，LIZhiwen2，LIULin3，4，LIANGWeihua4，WANGYinglong4

(1.Department of Physics,Hebei National Normal University, Chengde 067000, China;

2.Department of Management Engineering, Hebei National Normal University, Chengde 067000, China;

3.Department of Economy and Management, North China Institute of Aerospace Engineering; Langfang065000, China;4.College of Physics Science and Technology, Hebei University, Baoding 071002, China)

In a certain temperature range, using Landau-Devonshire free energy theory, the effects of strain on system entropy of both bulk and thin film of PbZr0.4Ti0.6O3(PZT) and [P(VDF-TrFE)](65/35) vinylidene fluoride-trifluoroethylene polymers and temperature change when entropy is maximum have been discussed. The results showed that as strain increases, temperature change caused by electric field and temperature when isothermal entropy change was maximum of polymers and PZT increased linearly and decreased linearly separately. To increase the entropy change of the films when they works at a low temperature, the strain of the polymers should be above zero, and the strain of PZT should be below zero.

thermoelectric body; electrocaloric effects; strain; PZT

10.3969/j.issn.1000-1565.2013.04.005

2012-12-08

河北省自然科學(xué)基金資助項(xiàng)目(E2012201035)

孫文祥(1964-)，男，河北滄州人，河北民族師范學(xué)院副教授，主要從事功能材料理論方向研究.

王英龍(1965-)，男，河北定州人，河北大學(xué)教授，主要從事功能材料方向研究.E-mail:hdwangyl@hbu.cn

O469

A

1000-1565(2013)04-0360-05

(責(zé)任編輯孟素蘭)