Optimization of LPE Growth Process of YIG Films Based on Multi-Physics Field Simulation

doi:10.16553/j.cnki.issn1000-985x.2025.0150

Abstract

Abstract: With the development of optical communication technology and photonic chips， yttrium iron garnet （YIG） crystals have shown broad application prospects in optical communication systems， magneto-optical isolators and other fields due to excellent magneto-optical properties. Liquid phase epitaxy （LPE） method， as a primary technique for preparing YIG films， has attracted significant attention. The key to growing YIG films via LPE lies in the precise control of solute transport， in which uniform flow field distribution and uniform temperature field environment are critical conditions to achieve consistent mass transfer. Therefore， it is of great research importance to regulate process parameters to stabilize the temperature field and flow field， thereby facilitating the growth of high-quality YIG films. In this paper， based on the method of multi-physics field simulation， a temperature-flow coupling model of YIG crystal grown by LPE method was established. The influence of crucible rotation on melt flow， axial temperature gradient and crystal growth rate was revealed by numerical simulation. Through systematic simulation and optimization， a set of optimal process parameters are obtained： the crucible rotation speed is controlled at 55~60 r/min， which could significantly suppress the disturbance of flow field and temperature field， thus maintaining the stability of the solid-liquid interface and providing favorable conditions for the high-quality epitaxial growth of YIG crystals. This study offers valuable theoretical insights and novel strategies for process optimization of controllable growth of complex oxide crystals.

Key words: crystal growth; liquid phase epitaxy; yttrium iron garnet; temperature-flow coupling; temperature field; flow field

CLC Number:

O734

LYU Bowen, WU Jiayu, ZHANG Hanxu, ZHU Senyin, ZHANG Lingli, ZHANG Yumin, WANG Xianjie, SONG Bo. Optimization of LPE Growth Process of YIG Films Based on Multi-Physics Field Simulation[J]. Journal of Synthetic Crystals, 2026, 55(1): 29-36.

Figures/Tables 6

References 23

[1]	朱森寅，张晗旭，王先杰，等. RIG型磁光薄膜的研究进展及应用［J］. 人工晶体学报， 2022， 51（9/10）： 1659-1672.
	ZHU S Y， ZHANG H X， WANG X J， et al. Research progress and application of magneto-optical thin films based on RIG［J］. Journal of Synthetic Crystals， 2022， 51（9/10）： 1659-1672 （in Chinese）.
[2]	GU S F， LI M， WU Z Z， et al. High-sensitivity and directional-identification fiber magnetic field sensor via polarimetric fiber ring laser with a fiber Faraday rotator［J］. Optical Fiber Technology， 2024， 82： 103643.
[3]	NADINOV I， KOVALENKO O， REHSPRINGER J L， et al. Functionalization of micro-size garnet at the end of optical fiber for magneto-optical applications［J］. Optics Express， 2024， 32（5）： 7651-7658.
[4]	SHIRAKASHI Z， GOTO T， TAKAGI H， et al. Reconstruction of non-error magnetic hologram data by magnetic assist recording［J］. Scientific Reports， 2017， 7： 12835.
[5]	WU Y J， YANG Q H， ZHANG D， et al. The submicron garnet film with perpendicular magnetic anisotropy prepared by liquid phase epitaxy method［J］. Journal of Magnetism and Magnetic Materials， 2020， 506： 166689.
[6]	WANG P F， LI M W， ZHOU C， et al. Numerical simulations of flow and mass transfer during large-scale potassium dihydrogen phosphate crystal growth via three-dimensional motion growth method［J］. International Journal of Heat and Mass Transfer， 2018， 127： 901-907.
[7]	TENG R， LI Y， CUI B， et al. Experiment and numerical simulation of melt convection and oxygen distribution in 400-mm Czochralski silicon crystal growth［J］. Rare Metals， 2017， 36（2）： 134-141.
[8]	林海鑫，高德东，王　珊，等. 大尺寸直拉硅单晶生长的多物理场建模与优化［J］. 人工晶体学报， 2025， 54（1）： 17-33.
	LIN H X， GAO D D， WANG S， et al. Multi-physics field modeling and optimization of large-size Czochralski silicon single crystal growth［J］. Journal of Synthetic Crystals， 2025， 54（1）： 17-33 （in Chinese）.
[9]	LIU B T， YU Y， TANG X， et al. Optimization of crucible and heating model for large-sized silicon carbide ingot growth in top-seeded solution growth［J］. Journal of Crystal Growth， 2020， 533： 125406.
[10]	GAŁĄZKA Z， WILKE H. Influence of Marangoni convection on the flow pattern in the melt during growth of Y₃Al₅O₁₂ single crystals by the Czochralski method［J］. Journal of Crystal Growth， 2000， 216（1/2/3/4）： 389-398.
[11]	WÜNSCHER M， MENZEL R， RIEMANN H， et al. Combined 3D and 2.5D modeling of the floating zone process with Comsol multiphysics［J］. Journal of Crystal Growth， 2014， 385： 100-105.
[12]	BRAESCU L. Influence of the dewetting phenomenon on the zinc distribution in CdZnTe crystals grown by dewetted Bridgman technique［C］//Photonic Fiber and Crystal Devices： Advances in Materials and Innovations in Device Applications IV. San Diego， California. SPIE， 2010.
[13]	CHOI D Y， CHUNG S J. Annealing behaviors of lattice misfit in YIG and La-doped YIG films grown on GGG substrates by LPE method［J］. Journal of Crystal Growth， 1998， 191（4）： 754-759.
[14]	YANG H D， ZHANG S. Numerical simulation of temperature field and stress field in fused deposition modeling［J］. Journal of Mechanical Science and Technology， 2018， 32（7）： 3337-3344.
[15]	LIU J， LU W Q. Preliminary study of non-isothermal phase change phenomena in vertical Bridgman crystal growth［J］. Chinese Science Bulletin， 2007， 52（5）： 701-710.
[16]	ZHANG N， LIU D. Numerical simulation of MHD oscillatory mixed convection in CZ crystal growth by Lattice Boltzmann method［J］. Results in Physics， 2018， 10： 882-890.
[17]	YAO M W， MATTHIESEN D H， CHAIT A. Numerical simulation of heat transport and fluid flow in directional crystal growth of GaAs［J］. Numerical Heat Transfer， Part A： Applications， 1996， 30（7）： 685-701.
[18]	施宇峰，王鹏飞，穆宏赫，等. 尺寸效应对坩埚下降法生长氟化钙晶体影响机制的数值模拟分析［J］. 人工晶体学报， 2024， 53（6）： 973-981.
	SHI Y F， WANG P F， MU H H， et al. Numerical simulation investigation of size effect on calcium fluoride crystals grown by vertical Bridgman method［J］. Journal of Synthetic Crystals， 2024， 53（6）： 973-981 （in Chinese）.
[19]	LIU L J， LIU X， WANG Y. Large-eddy simulation of melt turbulence in a 300-mm Cz-Si crystal growth［J］. International Journal of Heat and Mass Transfer， 2012， 55（1/2/3）： 53-60.
[20]	ROOK R， DOST S. The use of smoothed particle hydrodynamics for simulating crystal growth from solution［J］. International Journal of Engineering Science， 2007， 45（1）： 75-93.
[21]	TAUCHER L， RAMADAN Z， HAMMER R， et al. Machine learning assisted calibration of PVT simulations for SiC crystal growth［J］. CrystEngComm， 2024， 26（44）： 6322-6335.
[22]	曹　聪，刘江高，范叶霞，等. 碲锌镉晶体生长温度梯度与界面形状稳定性关系的研究［J］. 人工晶体学报， 2024， 53（4）： 641-648.
	CAO C， LIU J G， FAN Y X， et al. Relationship between temperature gradient and interfacial shape stability of CZT crystal growth［J］. Journal of Synthetic Crystals， 2024， 53（4）： 641-648 （in Chinese）.
[23]	PENG J Z， LIU X L， AUBRY N， et al. Data-driven modeling of geometry-adaptive steady heat conduction based on convolutional neural networks［J］. Case Studies in Thermal Engineering， 2021， 28： 101651.