專家信息：

吳忠慶，男，1974年生于浙江義烏，中國科學技術大學地球和空間科學學院教授。

教育及工作經歷：

1997年，獲四川大學水利水電工程系學士學位。

2000年，獲浙江大學物理系碩士學位。

2003年，獲清華大學物理系博士學位。

2003-2005年，獲清華大學高等研究中心博士后。

2005起先后在明尼蘇達大學化工材料系和南加州大學化工材料系博士后，用第一性原理計算研究材料在高溫高壓下的物性。 

研究方向：

1. 礦物的晶格動力學和熱力學性質。

2. 礦物的熱彈性性質和聲速。

3. 相變、雜質對礦物性質的影響。

4. 平衡同位素分餾系數。 

承擔科研項目情況：

1. Elasticity of ferropericlase at lower mantle PT condition

Elasticity of ferropericlase at 300K (red), 1000K (yellow), 2000K (green), 3000K(blue), and 4000K(purple)

Wu Z., Justo J. F., Wentzcovitch R., “Elastic Anomalies in a Spin-Crossover System: Ferropericlase at Lower Mantle Condition”, Physical Review Letters, 110, 228501 (2013).

Ferropericlase is second abundant mineral at lower mantle. Iron atoms of ferropericlase, which has magnetic moment at ambient pressure, will lose magnetic moment at lower mantle pressure. This transition is called high spin to low spin transition. This study demonstrates that spin transition significantly reduces the bulk modulus but has not such effect on shear modulus because of cubic symmetry of ferropericlase. But if cubic symmetry of ferropericase is broken by non-hydrostatic pressure, spin transition also reduces the shear modulus. This can well explain the discrepancy on shear modulus between two previous experimental measurements. The extrapolation of the experimental shear velocity of ferropericlase to higher temperature and concentration by Murakami et al [2012] introduced large error. The lower mantle composition based on their extrapolation need to be readjusted. Spin transition occurs at broad depth of lower mantle. Its dramatic effect on the elastic modulus must significantly impact our understanding to Earth’s interior. The systematical high pressure and temperature elasticity of ferropericlase obtained by this study will be very useful for us to uncover more the mystery of Earth’s interior. 

Spin effect on other properties of ferropericlase has been discussed in Wu et al (2009) and Wentzcovitch et al (2009). Spin transition of ferropericlase on mantle convection has been discussed in Shahnas et al (2011). 

2. An analytical approach for thermal elasticity of crystals

Elasticity of ferropericlase of forsterite

Wu Z., and Wentzcovitch, R. M, “Quasiharmonic thermal elasticity of crystals: An analytical approach”, Physical Review B, 83, 184115 (2011).

Thermal elastic properties of crystals are defined from second-order derivatives of free energies with respect to strain. In general, this will require calculating free energy under strain. These are exceedingly demanding computations requiring up to ∼103 parallel jobs running on tens or more processors each. Here, we developed an analytical formula to describe second-order derivatives of free energies with respect to strain. This approach decreases the computational effort, i.e., CPU time and human labor, by up to two orders of magnitude because the approach avoids calculating phonon density of states under strain. Results for the major mantle phases periclase (MgO) and forsterite (α-Mg2SiO4) show excellent agreement with previous first-principles results and experimental data. 

The approach has been used to calculated the elasticity of olivine, wadsleyite (Valdez et al GRL 2013), ringwoodite (Valdez et al EPSL 2012), and ferropericlase (Wu et al PRL 2013). 

3. A first principles-based method for anharmonic free energy

Thermal expansion of diamond with (purple) and without (green) anharmonic correction.

Wu Z., “Calculating the anharmonic free energy from first principles”, Physical Review B, 81, 172301 (2010). Wu Z., and Wentzcovitch, R. M, “Effective semiempirical ansatz for computing anharmonic free energies”, Physical Review B, 79, 104304 (2009).

We developed a method to compute anharmonic free energy neglected in quasiharmonic calculations by introducing parametrized temperature-dependent modification of the vibrational density of states. The introduced parameter is a constant and can be calculated through a first-principles molecular-dynamics simulation. The thermodynamics properties of forsterite,periclase, and diamond improve considerably after inclusion of anharmonic effects. 

The method has been used to develop the precise high temperature pressure standard of MgO (Wu et al JGR 2009).

4. Equilibrium fractionation factors of Mg among major upper mantle minerals

Application of clinopyroxene–garnet Mg isotope thermometry to estimating equilibrium temperature of eclogites. 

Huang F., Chen L., Wu Z., Wang W., “First-principles calculations of equilibrium Mg isotope fractionations between garnet, clinopyroxene, orthopyroxene, and olivine: implications for Mg isotope thermometry”, Earth and Planet Science Letters, 367 61-70,2013.

Equilibrium fractionation factors of Mg isotopes at high temperature and pressure are investigated for pyrope,majorite,clinopyroxene,orthopyroxene,andolivine using density functional theory. The results reveal significant inter- mineral Mg isotope fractionations due to different coordination environment of Mg in minerals. Specifically,pyrope,where Mg is in eight-fold coordination,is more enriched in light Mg isotopes than olivine and pyroxene where Mg is in six-fold coordination.The calculated isotope fractionation for Mg isotopes between clinopyroxene and pyrope is consistent with observations from natural eclogites (Li etal.,2011; Wang etal.,2012). Our calculation also reveals substantial pressure effect onMg isotope fractionations among garnet,pyroxene,and olivine. Equilibrium fractionation of Mg isotopes between pyrope and pyroxene(and olivine) could be used as a novel and independent thermometry with precision much better than the traditional garnet– clinopyroxene Fe Mg exchange thermometry because such fractionation is significantly greater than the current uncertainties of Mg isotope analyses. These Mg isotope thermometries could have wide applications in high-temperature mafic igneous and metamorphic rocks where garnet co-exists with pyroxene or olivine. 

科研成果：

1. 發現調制的鎬鈦酸鉛超晶格會出現具有優越鐵電性能的新相；在超薄的鐵電薄膜中偶極子會形成清晰的蝸旋結構的條文疇，這些發現幫助理解了薄膜的鐵電臨界尺寸效應、疇結構在電場下不同尋常的反應特性。

2. 利用第一性原理計算所得的MgO狀態方程與超高壓沖擊波實驗數據非常符合的特點，提出了一個MgO高溫高壓絕對壓標，并在2010年被基于絕對壓標方法的實驗全面驗證，為解決高壓實驗壓強不一致的問題邁出了關鍵一步。發展了兩個方法：第一個是第一性的計算非諧自由能方法，成倍的增加結果的適用溫度范圍；第二個是彈性常數的第一性原理計算新方法，其計算量不到常規的方法的十分之一，同時不降低計算精度，為計算眾多結構復雜的礦物的彈性奠定了基礎。這些有特色的研究，為吳忠慶教授在科研道路上的攀登鋪就了牢固的基石。

3. 鐵方鎂石是下地幔主要礦物之一，其所含鐵的原子磁距如何隨深度變化，以及進一步如何影響地幔特性是地球科學的一個基本問題。2003年首次在該礦物觀測到鐵的自旋轉變掀起了這個轉變對礦物物性特別是彈性影響的研究熱潮。吳忠慶教授首次獲得了整個地幔溫度和壓強下的鐵方鎂石彈性，從理論上指出鐵方鎂石的立方對稱性決定了鐵自旋轉變只會顯著降低體模量、對剪切模量沒有類似的效果。但在高壓實驗極可能出現的偏壓會破壞鐵方鎂石的立方對稱性，從而導致部分實驗觀測到剪切模量的反常降低，這很好地解釋了發表在《科學》雜志上不同工作間的分歧。在與地幔波速比較來限定地幔的成分時，通常需要將實驗測定的礦物波速數據往高溫高壓進行大范圍的外推，2012年Murakami等人借助于這種外推在《自然》雜志發文稱下地幔主要是鈣鈦礦組成的，與以往高溫高壓研究所得的結論完全不同，吳忠慶教授通過計算發現Murakami等人的外推顯著低估了鐵方鎂石橫波波速，他們關于下地幔的成分結論是不可靠的。研究結果發表在《物理評論快報》上，《科學》雜志在編輯特選欄目以《Spinning Iron in the Mantle》為題介紹了該工作。

4. 利用新發展的彈性方法研究各種礦物的彈性，包括上面自旋轉變方面的工作，其它關于不同鐵含量下的橄欖石及其高壓相的彈性特性已經分別發表在2012年的《地球和行星科學快報》和2013年的《地球物理研究快報》雜志上。

5. 第一性原理計算同位素平衡分餾方面的工作中，發現壓強、相變等可以顯著影響礦物間Mg、Si同位素分餾特性，揭示輝石和石榴子石之間大的Mg平衡同位素分餾方面，并據此建立的一個Mg同位素地質溫度計，上下地幔的主要礦物之間有很大的平衡Si同位素分餾，上下地幔的Si同位素組成因此很可能是不均一的，從而影響地核中Si含量的估計，這些工作分別發表在《地球和行星科學快報》和《地球化學與宇宙化學學報》上。 

論文專著：

已在PNAS，JGR，PRB等雜志發表學術論文30余篇。

代表性論文：

1. Wu, Z., Huang, F., Huang, S., "Isotope fractionation induced by phase transformation: First-principles investigation for Mg2SiO4", Earth and Planetary Science Letters, 409, 339-347, 2015.

2. Wu, Z., "themodynamic properties of wadsleyite with anharmonic effect", Earthquake science, in press.

3. Feng, C., Qin, T., Huang, S., Wu, Z., Fang Huang, "First-principles investigations of equilibrium calcium isotope fractionation between clinopyroxene and Ca-doped orthopyroxene", Geochimica et Cosmochimica Acta (2014).

4. Huang, F., Wu, Z., Shichun Huang, and Fei Wu. "First-principles calculations of equilibrium silicon isotope fractionation among mantle minerals", Geochimica et Cosmochimica Acta, 140 509-520, 2014.

5. Wu, Z., and Renata M. Wentzcovitch. "Spin crossover in ferropericlase and velocity heterogeneities in the lower mantle", Proceedings of the National Academy of Sciences 111, no.29 (2014): 10468-10472.(This work is highlighted by Science)

6. Rui Yang, Wu, Z.. "Elastic properties of stishovite and the CaCl2-type silica at the mantle temperature and pressure: An ab initio investigation", Earth and Planetary Science Letters, Volume 404, 15 October 2014, Pages 14-21, ISSN 0012-821X.

7. Wu Z., Justo J. F., Wentzcovitch R., "Elastic Anomalies in a Spin-Crossover System: Ferropericlase at Lower Mantle Condition", Physical Review Letters, 110, 228501 (2013). (This work is highlighted by Science)

8. Huang F., Chen L., Wu, Z., Wang W., "First-principles calculations of equilibrium Mg isotope fractionations between garnet, clinopyroxene, orthopyroxene, and olivine: implications for Mg isotope thermometry", Earth and Planet Science Letters, 367 61-70,

9. Valdez, M., Wu, Z., Yu, Y., Wentzcovitch, R. "Thermal elasticity of (Fex,Mg1-x)2SiO4 olivine and wadsleyite", Geophysical research letters, 40 290-294,2013. 

10. Valdez, M., Wu, Z., Yu, Y., Revenaugh, J., Wentzcovitch, R. "Thermoelastic properties of ringwoodite (Fex,Mg1-x)2SiO4: Its relationship to the 520 km seismic discontinuity", Earth and Planet Science Letters, 351 115-122, 2012. 

11. Wu, Z., and Wentzcovitch, R. M，Quasiharmonic thermal elasticity of crystals: An analytical approach", Physical Review B, 83, 184115 (2011).

12. Wu, Z., Nakano, A., Kalia, R. K., and Vashishta, P.," Vibrational and thermodynamic properties of β- HMX: a first-principles investigation", Journal of Chemical Physics, 134, 204509 (2011).

13. Wu, Z, Vashishta, P., Kalia, R., and Nakano, A., "First-principles calculations of the structural and dynamic properties, and the equation of state of crystalline iodine oxides I2O4, I2O5, and I2O6", Journal of Chemical Physics, 134, 204501 (2011).

14. Shahnas M. H.; Peltier W. R.; Wu, Z, and Wentzcovitch R M, "The high-pressure electronic spin transition in iron: Potential impacts upon mantle mixing", Journal of Geophysical Research，116, B08205,(2011). 

15. Wu, Z., "Calculating the anharmonic free energy from first principles",Physical Review B, 81, 172301 (2010).

16. Wentzcovitch, R. M., Yu, Y., and Wu, Z., (2010) "Thermodynamic Properties and Phase Relations in Mantle Minerals Investigated by First Principles Quasiharmonic Theory", Rev. Mineral. Geochem. 71,59-98.

17. Wentzcovitch, R. M., Yu, Y., and Wu, Z. "First Principles Quasiharmonic Thermoelasticity of Mantle Minerals", Reviews in Mineralogy & Geochemistry, 71 59,

18. Hsu, H., Umemoto, K., Wu, Z., Wentzcovitch, R. M., "Spin-State Crossover of Iron in Lower-Mantle Minerals: Results of DFT plus U Investigations", Reviews in Mineralogy & Geochemistry, 71 169, 2010.

19. Shimojo, F., Wu, Z., Nakano, A.,Kalia, RK., and Vashishta P.,"Density functional study of 1,3,5-trinitro-1,3,5-triazine molecular crystal with van der Waals interactions", Journal of Chemical Physics 132, 094106: 1-8, 2010 

20. Wu, Z., Justo, J. F., da Silva, C. R. S., de Gironcoli, S., and Wentzcovitch, R. M, "Anomalous thermodynamics properties of ferropericlase throughout its spin crossover transition", Physical Review B, 80, 014409 (2009).

21. Wu, Z., and Wentzcovitch, R. M, Effective semiempirical ansatz for computing anharmonic free energies", Physical Review B, 79, 104304 (2009).

22. Wentzcovitch, R. M., Justo, J. F., Wu, Z., da Silva, C. R. S., Yuen, D. A., and Kohlstedt D.," Anomalous compressibility of ferropericlase throughout the iron spin crossover", Proceedings of the National Academy of Sciences of the United States of America, 106:8447-8452 (2009). 

23. Yu, Y., Wu, Z., and Wentzcovitch, R. M.,"α to β to γ transformations in MgSiO4 in Earth's transition zone", Earth and Planet Science Letters, 273,115-122 (2008)

24. Wu, Z., Wentzcovitch, R.M., Umemoto, K., Li, B., Hirose, K., and J. Zheng, "PVT relations in MgO: an ultra-high PT scale for planetary sciences applications" Journal of Geophysical Research,113, B06204 (2008)

25. Sun, T., Umemoto, K., Wu, Z., Zheng J-C, and Wentzcovitch, R. M., "Lattice Dynamics and Thermal Equation of State of Platinum", Physical review B, 78, 024304 (2008)

26. Yu J., Wu, Z., Liu ZR, Yan QM, Wu J, Duan WH, "Phase diagram of ferroelectric BaTiO3 ultrathin films under open-circuit conditions", Journal of Physics-Condensed Matter, 20 135203(2008). 

27. Wu, Z., and Wentzcovitch, R. M., "Vibrational and thermodynamic properties of wadsleyite: a density functional study", Journal of Geophysical Research, 112, B12202 (2007)

28. Wu, Z., Duan WH, Huang ND, Wu J, and Gu BL, "Self-organization nanodomain structure in ferroelectric ultrathin films", Nanotechnology 18,325703(2007)

29. Wu, Z., Huang ND, Liu ZR, Duan WH, Wu J, and Gu BL, "Unusual vortex structure in ultrathin Pb(Zr0.5Ti0.5)O3 films", Journal of Applied Physsics 101,014112(2007). 

30. Wu, Z., Huang ND, Duan WH, Wu J, and Gu BL, "The phase diagram of ultrathin Pb(Zr0.5Ti0.5)O3 films under strain," Applied Physics Letters 86, 202903(2005).

31. Wu, Z., Duan WH, Wu J, Gu BL, and Zhang XW, "Dielectric properties of relaxor ferroelectric films," Journal of Applied Physics 98, 094105(2005).

32. Wu, Z., Huang ND, Liu ZR, Duan WH, Wu J, Gu BL and Zhang XW, "Ferroelectricity in Pb(Zr0.5Ti0.5)O3 thin films: critical thickness and 180o stripe domain," Physical Review B 70, 104108 (2004).

33. Huang ND, Liu ZR, Wu, Z., Duan WH, Gu BL and Zhang XW, "Huge enhancement of electromechanical responses in compositionally modulated Pb(Zr1-xTix)O3," Physical Review Letters 91, 067602 (2003).

34. Wu, Z., Duan WH, Wang Y, Gu BL and Zhang XW, "Effect of defect-induced internal field on the aging of relaxors," Physical Review B 67, 052101 (2003).

35. Wu, Z., Duan WH, Liu ZR, Gu BL and Zhang XW, "Effect of tunneling frequency on relaxor behavior," Microelectric Engineering 66: 676-682 (2003).

36. Wu, Z., Duan WH, Liu ZR, Gu BL and Zhang XW, "Multiple tunneling channels order-disorder ferroelectric model and field-induced phase transition in relaxors," Physical Review B 65, 174119 (2002).

37. Wu, Z. and Zhang XJ, "Charge fluctuations in cuprate superconductors," arxiv: cond-mat/0210455.

38. Wu, Z., Liu ZR, and Gu BL, "Order-Disorder phase transition and dielectric mechanism in relaxor ferroelectrics," Tsinghua Science and Technology 6, 97-108 (2001).

榮譽獎勵：

資料更新中……

Phy.Rev.Lett.發表了我們組關于“自旋轉變下鐵鐵方鎂石高溫高壓下彈性”的工作，
Science 在編輯特選欄目以"Spinning Iron in the Mantle"為題介紹了該工作。

 日前，我校地球和空間科學學院地震與地球內部物理實驗室吳忠慶教授與來自美國明尼蘇達大學的研究人員合作，利用第一性原理計算研究了高溫高壓下鐵自旋 轉變對鐵方鎂石高溫高壓彈性的影響，解釋了不同實驗在自旋轉變對橫波波速影響存在分歧的成因，指出Murakami 等人［Nature 485，90（2012）］外延的鐵方鎂石橫波波速數據 存在極大的誤差，由此得出的下地幔主要是鈣鈦礦的結論是不可靠的。這些系統完整的鐵方鎂石高溫高壓彈性數據為我們進一步利用自旋轉變效應限定地球內部成分提供了關鍵的基礎。研究結果近期發表在《物理評論快報》（Phys. Rev. Lett. 110, 228501 (2013)）上，Science雜志在編輯特選欄目以“Spinning Iron in the Mantle”為題介紹了該工作。

地幔中含有大量的鐵，這些鐵原子的磁矩在高壓下是否會發生改變，以及這種改變對礦物性質會帶來怎么樣的影響，一直是地球科學研究領域的基本問題。2003和2004年，Badro等人從實驗上觀測到下地幔兩大主要礦物——鈣鈦礦和鐵方鎂石——中的鐵都會在下地幔的壓強下經歷高自旋到低自旋的轉變，掀起了自旋轉變對礦物特性包括彈性影響的研究熱潮。吳忠慶教授和合作者用第一性原理計算研究了自旋轉變對彈性的影響。第一性原理計算是一個解量子力學方程的方法，它無須引入任何經驗參數，較易實現高溫高壓條件且所得的結果跟實驗可以媲美，已成為獲得礦物高溫高壓物性的重要方法。理論分析顯示，自旋轉變會顯著降低縱波波速，但類似的現象不會在橫波中發生，這是由鐵方鎂石的立方對稱性決定的。但靜壓環境很難完全在高壓實驗實現，偏壓會破壞鐵方鎂石的立方對稱性，從而出現橫波波速的反常降低。這很好地解釋了不同實驗觀察到橫波波速行為的分歧，也跟這些實驗所采用的傳壓媒質維持靜壓條件的能力一致。

通過比較各種礦物在地幔溫壓條件下的波速和地球深部的波速可以最直接有效的限定地幔的組分。然而，受限于高壓高溫等實驗條件，相關礦物在地幔溫度和壓強條件下的物性的測量數據非常很有限，必須將這些有限的實驗數據進行外延才能討論下地幔的成分等關鍵問題。2012年，Murakami 等人借助這種外延得出了下地幔主要是鈣鈦礦組成的結論(>93%)，與以往高溫高壓研究所得的結論（地幔巖，~80%鈣鈦礦）完全不同。避免外延的的一個有效方法是利用第一性原理計算。吳忠慶教授和合作者用第一性原理計算發現Murakami 等人的外延數據有問題。因為計算的結果跟所有測量的實驗數據都能很好符合，但Murakami 等人的外延鐵方鎂石橫波波速則普遍比計算的結果要小9%。計算結果表明地幔巖（pyrolite模型）是下地幔很好的組分模型，結果說明直接測量地幔對應的溫度壓強和成分條件下的物性以避免外延的必要性和迫切性。

本研究課題是在美國自然科學基金和國家自然科學基金面上項目及中科院外國專家局“創新團隊國際合作伙伴計劃”的資助下完成的。

來源：中國科學技術大學地球和空間科學學院地震與地球內部物理實驗室
 

中國科大研究揭示地幔中部結構成因

中國科學技術大學地球和空間科學學院地震與地球內部物理實驗室吳忠慶教授合作研究揭示，下地幔礦物中鐵的自旋轉變會導致地幔中部的縱波波速對溫度變化不敏感，解釋了多個地震層析成像觀測到的不尋常波速結構的成因，從而為利用自旋轉變效應認識地球內部結構打開了一扇門。工作發表在近期的《美國科學院院刊》上，《科學》雜志的編輯特選欄目以“Ample explanation for seismic variation”為題介紹了該工作。這是繼2013年吳忠慶教授等人在《物理評論快報》發表自旋轉變對鐵方鎂石彈性不尋常的影響的工作后的又一個重要進展。

自旋轉變下熱異常產生的P波波速結構示意圖

鐵方鎂石中的鐵從有磁矩的高自旋態轉變到沒有磁矩的低自旋態會顯著降其波速，但不同于其它礦物相變有一個明顯的波速間斷面，鐵方鎂石中的鐵是逐漸從高自旋態過渡到低自旋態，是一個平滑的相變，這導致利用地震學手段探測該轉變的努力沒有成功，該轉變被認為是地震學不可見的。吳教授仔細分析了自旋轉變下的彈性數據，發現自旋轉變可在地震層析成像圖中留下可觀測的多個特征，例如它會導致縱波波速在~1750公里深度對溫度變化不敏感，其后果就是起源于深部的地幔柱會在~1750公里深度中斷，這樣的現象已經在夏威夷等地多個熱點下的地幔柱縱波成像圖中觀測到，但一直沒有得到很好的理解。地幔柱是起源于核幔邊界的熱物質上涌流，穿越整個地幔到達巖石圈底部，是板塊運動的主要驅動力之一，能夠解釋地表熱點軌跡，大火成巖省形成等眾多地質觀測現象，因此直接證實地幔柱是地球科學非常重要的一個課題，利用熱物質波速慢的特點，地震學成像是目前觀測地幔柱主要手段。但實際觀測發現多個地幔柱在～1750公里深度波速并沒有明顯變慢，既上涌熱物質在這個深度左右突然變的不熱，這很難理解也跟地幔柱的整個概念有沖突，而吳教授的研究結果表明上涌物質其實還是熱的，只不過自旋轉變下縱波波速在這個深度對溫度變化不敏感，給人溫度不高的錯覺（圖1），吳教授的工作給地幔柱模型強有力的支持。由于找到了識別自旋轉變的特征，可以預期自旋轉變一定會象橄欖石系列相變、后鈣鈦礦相變一樣有力地促進了我們對地球內部結構的認識。

吳忠慶教授系列關于自旋轉變的工作受到中國自然科學基金、973計劃、中科院海外創新團隊等基金的支持，并應邀在AOGS2014, AGU做邀請報告。

來源：中國科學技術大學新聞網 2014-08-05

探秘地球深處

——記中國科學技術大學地球和空間科學學院吳忠慶教授 

電影《地心歷險記》講述了一位地質學家深信地底另有世界，于是與侄兒前往冰島，從火山口潛入地心，體驗了地下不為人知的幻妙世界�；氐浆F實，我們所生活的地球內部究竟是怎樣的？

科學發展到今天，人們能上“九天攬月”卻仍然“入地無門”。因為那是一個復雜的高溫高壓系統，地球核心的壓力約360Pa，溫度超過5000℃，研究物質在高溫高壓下物性因此也就成為探索地球深部物質組成、結構和變化規律的關鍵一環。從全球范圍來看，高壓科研仍處于萌芽階段，中國起步更晚。正是在這一需求下，吳忠慶教授作為計算研究礦物高溫高壓物性的海外杰出人才被引進回國。 

初嘗科研快樂 

吳忠慶教授1997年畢業于四川大學水利水電工程系獲學士學位，源于發自心底對物理的興趣，同年考入浙江大學物理系攻讀碩士學位，在這里受到高溫超導研究的耳濡目染，重新思考高溫超導體中的空穴配對的概念，提出新觀點解釋高溫超導的多個實驗結果，初嘗獲得新發現的快樂，感受到科學研究迷人的魅力。自此，也堅定了從事科學研究的決心。

2000年，吳忠慶考入清華大學物理系攻讀博士學位，在顧秉林院士和段文輝教授組從事鐵電材料特性的計算研究。在寬松、自由、活躍的學術氛圍下，他發現調制的鎬鈦酸鉛超晶格會出現具有優越鐵電性能的新相；在超薄的鐵電薄膜中偶極子會形成清晰的蝸旋結構的條文疇，這些發現幫助理解了薄膜的鐵電臨界尺寸效應、疇結構在電場下不同尋常的反應特性。

盡管當時吳忠慶教授在鐵電方面研究工作取得了可喜的成果，但他卻不安于現狀，在導師顧秉林院士和段文輝教授的鼓勵下，決定出國留學，開拓視野、拓寬研究領域。

探索高溫高壓下的礦物物性 

2005年6月，吳忠慶教授來到美國明尼蘇達大學化工材料系進行博士后研究，開始從事第一性原理計算研究礦物高溫高壓下的物性。合作導師Renata Wentzcovitch教授是第一性原理計算研究礦物高溫高壓物性領域的杰出學者，近些年來領導的課題組在礦物物理領域的兩大進展——后鈣鈦礦相變和鐵的自旋轉變中做出了重要貢獻。其中，計算礦物彈性更是她領導的研究組的一大特點。

吳忠慶教授介紹說，諸如層析成像等地震學反演是他們獲得地球內部結構的主要手段，但要解讀所得的波速結果需要他們對眾多礦物在地球內部溫壓下的性質有系統的了解，比較礦物聚合體在高溫高壓下的波速與地球內部波速是最直接限定地球內部成分的方法。礦物高溫高壓下的物性特別是彈性特性，因此是礦物物理學的核心內容。與高壓實驗相比，第一性原理計算可以容易地實現高溫高壓條件，所得結果可以跟實驗媲美，在研究中扮演越來越重要的角色。

在國外，吳忠慶教授在該領域取得多個重要成果：例如，利用第一性原理計算所得的MgO狀態方程與超高壓沖擊波實驗數據非常符合的特點，提出了一個MgO高溫高壓絕對壓標，并在2010年被基于絕對壓標方法的實驗全面驗證，為解決高壓實驗壓強不一致的問題邁出了關鍵一步。發展了兩個方法：第一個是第一性的計算非諧自由能方法，成倍的增加結果的適用溫度范圍；第二個是彈性常數的第一性原理計算新方法，其計算量不到常規的方法的十分之一，同時不降低計算精度，為計算眾多結構復雜的礦物的彈性奠定了基礎。這些有特色的研究，為吳忠慶教授在科研道路上的攀登鋪就了牢固的基石。 

回國后再次起航 

鑒于國內高溫高壓礦物物性研究是一個相對比較薄弱的領域，但同時又是地球科學非常重要的一個方向，中國科學技術大學地球和空間科學學院陳曉非院長力主在2010年引進吳忠慶教授。在學校和學院的大力支持下，他圍繞深化地球內部結構認識這個主題開展研究，在幾年時間里在該領域做出了多個有影響的成果，為迅速提升我國在該領域的國際影響力打下了一個好的開局。

鐵方鎂石是下地幔主要礦物之一，其所含鐵的原子磁距如何隨深度變化，以及進一步如何影響地幔特性是地球科學的一個基本問題。2003年首次在該礦物觀測到鐵的自旋轉變掀起了這個轉變對礦物物性特別是彈性影響的研究熱潮。吳忠慶教授首次獲得了整個地幔溫度和壓強下的鐵方鎂石彈性，從理論上指出鐵方鎂石的立方對稱性決定了鐵自旋轉變只會顯著降低體模量、對剪切模量沒有類似的效果。但在高壓實驗極可能出現的偏壓會破壞鐵方鎂石的立方對稱性，從而導致部分實驗觀測到剪切模量的反常降低，這很好地解釋了發表在《科學》雜志上不同工作間的分歧。在與地幔波速比較來限定地幔的成分時，通常需要將實驗測定的礦物波速數據往高溫高壓進行大范圍的外推，2012年Murakami等人借助于這種外推在《自然》雜志發文稱下地幔主要是鈣鈦礦組成的，與以往高溫高壓研究所得的結論完全不同，吳忠慶教授通過計算發現Murakami等人的外推顯著低估了鐵方鎂石橫波波速，他們關于下地幔的成分結論是不可靠的。研究結果發表在《物理評論快報》上，《科學》雜志在編輯特選欄目以《Spinning Iron in the Mantle》為題介紹了該工作。

鐵方鎂石的鐵自旋轉變是一個漸變的過程，不會導致波速跳變，因此似乎是地震學不可見的，吳忠慶教授通過分析鐵方鎂石彈性數據發現，鐵的自旋轉變產生的效果完全可以通過層析成像等手段觀測到，例如它會導致P波波速在～1750公里深度對溫度變化不敏感，起源于深部的地幔柱在該深度“斷開”，這樣的現象已經在多個層析成像結果中出現，但一直沒有得到很好的理解。地幔柱是起源于核幔邊界的熱物質上涌流，穿越整個地幔到達巖石圈底部，是板塊運動的主要驅動力之一，能夠解釋地表熱點軌跡，大火成巖省形成等眾多地質觀測現象，因此直接證實地幔柱是地球科學非常重要的一個課題，地震學成像是目前觀測地幔柱主要手段，利用了熱物質波速慢的特點。但實際觀測發現多個地幔柱在～1750公里深度波速并沒有明顯變慢，既上涌熱物質在這個深度左右突然變的不熱，這是很難理解的，也跟地幔柱的整個概念有沖突，而吳忠慶教授的研究結果表明上涌物質在該深度其實還是熱的，只不過由于鐵自旋轉變，P波波速在這個深度對溫度變化不敏感，在P波圖像上感受不到高溫。吳教授的工作給地幔柱模型強有力的支持。可以預期自旋轉變一定會象橄欖石系列相變、后鈣鈦礦相變一樣有力地促進了我們對地球內部結構的認識。研究結果發表在《美國科學院院刊》。

吳忠慶教授的研究工作重點主要是兩大方面。一是利用新發展的彈性方法研究各種礦物的彈性，這包括上面自旋轉變方面的工作，其它關于不同鐵含量下的橄欖石及其高壓相的彈性特性已經分別發表在2012年的《地球和行星科學快報》和2013年的《地球物理研究快報》雜志上。另一個是第一性原理計算同位素平衡分餾方面的工作，發現壓強、相變等可以顯著影響礦物間Mg、Si同位素分餾特性，揭示輝石和石榴子石之間大的Mg平衡同位素分餾方面，并據此建立的一個Mg同位素地質溫度計，上下地幔的主要礦物之間有很大的平衡Si同位素分餾，上下地幔的Si同位素組成因此很可能是不均一的，從而影響地核中Si含量的估計，這些工作分別發表在《地球和行星科學快報》和《地球化學與宇宙化學學報》上。

地球科學是一個交叉性特別強的學科，吳忠慶教授在研究中非常關注這一特點，在中科大地空學院自由的學術氛圍中積極地開展各種合作。

在認識自然界的過程中，充滿各種挑戰和困難；在探秘地球深處的過程中，充滿無限拓新和神奇，而這些探索對地球形成與演化規律、開發地球資源及減輕自然災害（如地震等）是有著重要意義的。我們期待吳忠慶教授在探索未知的路上有更多發現和進展！

來源：科學中國人 2014年第8期