CN114662424A - Thin film deposition simulation method, device, electronic device and storage medium - Google Patents
Thin film deposition simulation method, device, electronic device and storage medium Download PDFInfo
- Publication number
- CN114662424A CN114662424A CN202210548370.3A CN202210548370A CN114662424A CN 114662424 A CN114662424 A CN 114662424A CN 202210548370 A CN202210548370 A CN 202210548370A CN 114662424 A CN114662424 A CN 114662424A
- Authority
- CN
- China
- Prior art keywords
- thin film
- deposition
- phase field
- model
- film deposition
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Laminated Bodies (AREA)
Abstract
本发明涉及数据模拟分析技术领域,具体公开了一种薄膜沉积模拟方法、装置、电子设备及存储介质,其中,方法包括以下步骤:设定代表薄膜凝固情况的第一相场参量及代表输入蒸汽局部密度的第二相场参量;根据第一相场参量和第二相场参量建立用于模拟物理气相沉积的弹道沉积相场模型;根据第一相场参量耦合弹道沉积相场模型和固态薄膜的杨氏模量,以使弹道沉积相场模型与力学平衡方程耦合而建立薄膜沉积模型;根据薄膜沉积模型进行薄膜沉积模拟;该方法实现薄膜沉积过程的预测模拟,且能根据模拟结果反映薄膜沉积形态和应力分布特点。
The invention relates to the technical field of data simulation analysis, and specifically discloses a thin film deposition simulation method, device, electronic device and storage medium, wherein the method comprises the following steps: setting a first phase field parameter representing the solidification condition of the thin film and a representative input steam The second phase field parameters of local density; the ballistic deposition phase field model for simulating physical vapor deposition is established according to the first phase field parameters and the second phase field parameters; the ballistic deposition phase field model and the solid thin film are coupled according to the first phase field parameters The Young's modulus of the film is obtained by coupling the ballistic deposition phase field model with the mechanical equilibrium equation to establish the film deposition model; the film deposition simulation is carried out according to the film deposition model; the method realizes the prediction and simulation of the film deposition process, and can reflect the film deposition process according to the simulation results. Deposit morphology and stress distribution characteristics.
Description
技术领域technical field
本申请涉及数据模拟分析技术领域,具体而言,涉及一种薄膜沉积模拟方法、装置、电子设备及存储介质。The present application relates to the technical field of data simulation analysis, and in particular, to a thin film deposition simulation method, device, electronic device and storage medium.
背景技术Background technique
对器件进行镀膜可以大大提升器件的性能表现,镀膜技术在各种工业领域、高精器件制造领域中都有着重要的应用。Coating the device can greatly improve the performance of the device, and the coating technology has important applications in various industrial fields and high-precision device manufacturing.
物理气相沉积(physical vapor deposition,PVD)技术是沉积薄膜的重要技术手段。Physical vapor deposition (PVD) technology is an important technical means for depositing thin films.
为了正确研究和理解薄膜特性与气相沉积条件之间的联系,现有技术一般通过实验手段来了解不同材料在不同沉积条件下的薄膜形成过程,但实验需依赖昂贵的实验仪器并消耗大量实验耗材,若采用数字化模拟方式则需要充分考虑及相、晶粒和晶界的形成等互相关联的变化因素,现有技术缺少提供一种能准确描述关于材料特性和薄膜演变关系的模拟模型。In order to correctly study and understand the relationship between film properties and vapor deposition conditions, the existing technology generally uses experimental means to understand the film formation process of different materials under different deposition conditions, but the experiments need to rely on expensive experimental instruments and consume a lot of experimental consumables If the digital simulation method is adopted, the interrelated change factors such as the formation of phases, grains and grain boundaries need to be fully considered. The prior art lacks a simulation model that can accurately describe the relationship between material properties and film evolution.
针对上述问题,目前尚未有有效的技术解决方案。There is currently no effective technical solution for the above-mentioned problems.
发明内容SUMMARY OF THE INVENTION
本申请的目的在于提供一种薄膜沉积模拟方法、装置、电子设备及存储介质,实现薄膜沉积过程的预测模拟,且能根据模拟结果反映薄膜沉积形态和应力分布特点。The purpose of this application is to provide a film deposition simulation method, device, electronic device and storage medium, which can realize the prediction simulation of the film deposition process, and can reflect the film deposition morphology and stress distribution characteristics according to the simulation results.
第一方面,本申请提供了一种薄膜沉积模拟方法,用于建立数字化薄膜模型以模拟薄膜沉积过程,所述方法包括以下步骤:In a first aspect, the present application provides a thin film deposition simulation method for establishing a digital thin film model to simulate a thin film deposition process, the method comprising the following steps:
设定代表薄膜凝固情况的第一相场参量及代表输入蒸汽局部密度的第二相场参量;Set the first phase field parameter representing the solidification of the film and the second phase field parameter representing the local density of the input steam;
根据所述第一相场参量和第二相场参量建立用于模拟物理气相沉积的弹道沉积相场模型;establishing a ballistic deposition phase field model for simulating physical vapor deposition according to the first phase field parameter and the second phase field parameter;
根据所述第一相场参量耦合所述弹道沉积相场模型和固态薄膜的杨氏模量,以使所述弹道沉积相场模型与力学平衡方程耦合而建立薄膜沉积模型;Coupling the ballistic deposition phase field model and the Young's modulus of the solid thin film according to the first phase field parameter, so as to couple the ballistic deposition phase field model with the mechanical equilibrium equation to establish a thin film deposition model;
根据所述薄膜沉积模型进行薄膜沉积模拟。Thin film deposition simulations are performed according to the thin film deposition model.
本申请的薄膜沉积模拟方法耦合弹道沉积相场模型和固态薄膜的杨氏模量,使最终获取的薄膜沉积模型能反映薄膜的力学性能变化,能将薄膜模拟沉积过程中的相场变化和薄膜力学性能结合模拟,实现薄膜沉积过程的预测模拟。The thin film deposition simulation method of the present application couples the ballistic deposition phase field model and the Young's modulus of the solid thin film, so that the finally obtained thin film deposition model can reflect the change of the mechanical properties of the thin film, and can simulate the phase field change and thin film during the deposition process of the thin film. The mechanical properties are combined with the simulation to realize the prediction simulation of the thin film deposition process.
所述的薄膜沉积模拟方法,其中,所述根据所述第一相场参量耦合所述弹道沉积相场模型和固态薄膜的杨氏模量,以使所述弹道沉积相场模型与力学平衡方程耦合而建立薄膜沉积模型的步骤包括:The thin film deposition simulation method, wherein the ballistic deposition phase field model and the Young's modulus of the solid thin film are coupled according to the first phase field parameters, so that the ballistic deposition phase field model and the mechanical equilibrium equation are The steps to build a thin film deposition model by coupling include:
根据所述第一相场参量获取所述弹道沉积相场模型的局部固体密度参量,耦合所述局部固体密度参量和所述固态薄膜的杨氏模量,获取关于所述局部固体密度参量的杨氏模量;Obtain the local solid density parameter of the ballistic deposition phase field model according to the first phase field parameter, couple the local solid density parameter and the Young's modulus of the solid thin film, and obtain the Young's modulus of the local solid density parameter modulus;
根据所述关于所述局部固体密度参量的杨氏模量和预设的位移场获取应力方程;obtaining a stress equation according to the Young's modulus with respect to the local solid density parameter and a preset displacement field;
根据所述应力方程和所述力学平衡方程生成所述薄膜沉积模型。The thin film deposition model is generated from the stress equation and the mechanical equilibrium equation.
所述的薄膜沉积模拟方法,其中,所述根据所述第一相场参量获取所述弹道沉积相场模型的局部固体密度参量的步骤包括:The method for simulating thin film deposition, wherein the step of obtaining the local solid density parameter of the ballistic deposition phase field model according to the first phase field parameter includes:
定义所述局部固体密度参量为f d ,所述局部固体密度参量满足:The local solid density parameter is defined as f d , and the local solid density parameter satisfies:
f d =(f+f max )/2; f d = ( f + f max )/2;
f为所述第一相场参量,f max 为所述第一相场参量的最大值。 f is the first phase field parameter, and f max is the maximum value of the first phase field parameter.
在该示例的薄膜沉积模拟方法中,f d 能表征对应区域中的固态薄膜的占比,即f d 能充当相场中固相的体积分数。In the thin film deposition simulation method of this example, f d can characterize the proportion of solid thin film in the corresponding region, that is, f d can act as the volume fraction of solid phase in the phase field.
所述的薄膜沉积模拟方法,其中,所述局部固体密度参量的杨氏模量具有残余模量系数。The thin film deposition simulation method, wherein the Young's modulus of the local solid density parameter has a residual modulus coefficient.
在气相模量可以忽略的情况下,本申请的方法设置残余模量系数能使获取的薄膜沉积模型避免数值奇异性。In the case where the gas-phase modulus is negligible, the method of the present application sets the residual modulus coefficient to avoid numerical singularities in the obtained thin film deposition model.
所述的薄膜沉积模拟方法,其中,所述根据所述薄膜沉积模型进行薄膜沉积模拟的步骤包括:The method for simulating thin film deposition, wherein the step of performing thin film deposition simulation according to the thin film deposition model includes:
根据反应条件及薄膜材质设定所述弹道沉积相场模型的现象学参数及所述固态薄膜的杨氏模量和泊松比;Setting the phenomenological parameters of the ballistic deposition phase field model and the Young's modulus and Poisson's ratio of the solid film according to the reaction conditions and the film material;
将所述弹道沉积相场模型的现象学参数及所述固态薄膜的杨氏模量和泊松比输入至所述薄膜沉积模型中,并利用所述薄膜沉积模型模拟薄膜沉积生成薄膜图像信息。The phenomenological parameters of the ballistic deposition phase field model and the Young's modulus and Poisson's ratio of the solid thin film are input into the thin film deposition model, and the thin film deposition model is used to simulate thin film deposition to generate thin film image information.
该示例的薄膜沉积模拟方法利用该沉积模型生成薄膜图像信息能直观地显示薄膜演变过程,便于用户了解材料特性与薄膜形成过程的关系,以进行实验或生产的控制调整。The film deposition simulation method of this example uses the deposition model to generate film image information, which can visually display the film evolution process, which is convenient for users to understand the relationship between the material properties and the film formation process, so as to control and adjust the experiment or production.
所述的薄膜沉积模拟方法,其中,所述根据所述薄膜沉积模型进行薄膜沉积模拟的步骤还包括:The method for simulating thin film deposition, wherein the step of performing thin film deposition simulation according to the thin film deposition model further comprises:
利用所述薄膜沉积模型模拟薄膜沉积生成应力分布图像信息。Using the thin film deposition model to simulate thin film deposition to generate stress distribution image information.
所述的薄膜沉积模拟方法,其中,所述弹道沉积相场模型包括关于所述第一相场参量的薄膜相控制方程和关于所述第二相场参量的蒸汽控制方程。The method for simulating thin film deposition, wherein the ballistic deposition phase field model includes a thin film phase governing equation for the first phase field parameter and a vapor governing equation for the second phase field parameter.
第二方面,本申请还提供了一种薄膜沉积模拟装置,用于建立数字化薄膜模型以模拟薄膜沉积过程,所述装置包括:In a second aspect, the present application also provides a thin film deposition simulation device for establishing a digital thin film model to simulate a thin film deposition process, the device comprising:
设定模块,用于设定代表薄膜凝固情况的第一相场参量及代表输入蒸汽局部密度的第二相场参量;a setting module for setting a first phase field parameter representing the solidification condition of the thin film and a second phase field parameter representing the local density of the input steam;
建模模块, 用于根据所述第一相场参量和第二相场参量建立用于模拟物理气相沉积的弹道沉积相场模型;a modeling module, configured to establish a ballistic deposition phase field model for simulating physical vapor deposition according to the first phase field parameter and the second phase field parameter;
耦合模块,用于根据所述第一相场参量耦合所述弹道沉积相场模型和固态薄膜的杨氏模量,以使所述弹道沉积相场模型与力学平衡方程耦合而建立薄膜沉积模型;a coupling module, configured to couple the ballistic deposition phase field model and the Young's modulus of the solid thin film according to the first phase field parameter, so as to couple the ballistic deposition phase field model with the mechanical equilibrium equation to establish a thin film deposition model;
模拟模块,用于根据所述薄膜沉积模型进行薄膜沉积模拟。The simulation module is used for performing thin film deposition simulation according to the thin film deposition model.
本申请的薄膜沉积模拟装置耦合弹道沉积相场模型和固态薄膜的杨氏模量,使最终获取的薄膜沉积模型能反映薄膜的力学性能变化,能将薄膜模拟沉积过程中的相场变化和薄膜力学性能结合模拟,实现薄膜沉积过程的预测模拟。The thin film deposition simulation device of the present application couples the ballistic deposition phase field model and the Young's modulus of the solid thin film, so that the finally obtained thin film deposition model can reflect the mechanical property change of the thin film, and can simulate the phase field change and thin film during the thin film deposition process. The mechanical properties are combined with the simulation to realize the prediction simulation of the thin film deposition process.
第三方面,本申请还提供了一种电子设备,包括处理器以及存储器,所述存储器存储有计算机可读取指令,当所述计算机可读取指令由所述处理器执行时,运行如上述第一方面提供的所述方法中的步骤。In a third aspect, the present application also provides an electronic device, including a processor and a memory, where the memory stores computer-readable instructions, when the computer-readable instructions are executed by the processor, the operation is as described above The steps in the method provided by the first aspect.
第四方面,本申请还提供了一种存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时运行如上述第一方面提供的所述方法中的步骤。In a fourth aspect, the present application further provides a storage medium on which a computer program is stored, and when the computer program is executed by a processor, the steps in the method provided in the above-mentioned first aspect are executed.
由上可知,本申请提供了一种薄膜沉积模拟方法、装置、电子设备及存储介质,其中,方法基于代表薄膜凝固情况的第一相场参量和固态薄膜的杨氏模量的关系耦合弹道沉积相场模型和固态薄膜的杨氏模量,使最终获取的薄膜沉积模型能反映薄膜的力学性能变化,能将薄膜模拟沉积过程中的相场变化和薄膜力学性能结合模拟,实现薄膜沉积过程的预测模拟,且能根据模拟结果反映薄膜沉积形态和应力分布特点。As can be seen from the above, the present application provides a thin film deposition simulation method, device, electronic device and storage medium, wherein the method is based on the relationship between the first phase field parameter representing the solidification condition of the thin film and the Young's modulus of the solid thin film coupled with ballistic deposition The phase field model and the Young's modulus of the solid-state thin film enable the final film deposition model to reflect the mechanical properties of the thin film. Predictive simulation, and can reflect the film deposition morphology and stress distribution characteristics according to the simulation results.
附图说明Description of drawings
图1为本申请实施例提供的薄膜沉积模拟方法的流程图。FIG. 1 is a flowchart of a thin film deposition simulation method provided by an embodiment of the present application.
图2(a)为根据本申请实施例提供的薄膜沉积模拟方法获取的薄膜沉积模型的薄膜沉积过程中起始阶段的薄膜图像示意图。FIG. 2( a ) is a schematic diagram of a thin film image at an initial stage of a thin film deposition process of a thin film deposition model obtained according to the thin film deposition simulation method provided in an embodiment of the present application.
图2(b)为根据本申请实施例提供的薄膜沉积模拟方法获取的薄膜沉积模型的薄膜沉积过程中第一中间阶段的薄膜图像示意图。FIG. 2( b ) is a schematic diagram of a thin film image in the first intermediate stage of the thin film deposition process of the thin film deposition model obtained according to the thin film deposition simulation method provided in the embodiment of the present application.
图2(c)为根据本申请实施例提供的薄膜沉积模拟方法获取的薄膜沉积模型的薄膜沉积过程中第二中间阶段的薄膜图像示意图。FIG. 2( c ) is a schematic diagram of a thin film image of the second intermediate stage of the thin film deposition process of the thin film deposition model obtained according to the thin film deposition simulation method provided in the embodiment of the present application.
图2(d)为根据本申请实施例提供的薄膜沉积模拟方法获取的薄膜沉积模型的薄膜沉积过程中结束阶段的薄膜图像示意图。FIG. 2(d) is a schematic diagram of a thin film image at the end stage of the thin film deposition process of the thin film deposition model obtained according to the thin film deposition simulation method provided in the embodiment of the present application.
图3为与图2(d)的薄膜图像对应的应力分布示意图。Fig. 3 is a schematic diagram of the stress distribution corresponding to the thin film image of Fig. 2(d).
图4为本申请实施例提供的薄膜沉积模拟装置的结构示意图。FIG. 4 is a schematic structural diagram of a thin film deposition simulation apparatus provided by an embodiment of the present application.
图5为本申请实施例提供的电子设备的结构示意图。FIG. 5 is a schematic structural diagram of an electronic device provided by an embodiment of the present application.
附图标记:201、设定模块;202、建模模块;203、耦合模块;204、模拟模块;301、处理器;302、存储器;303、通信总线。Reference numerals: 201, setting module; 202, modeling module; 203, coupling module; 204, simulation module; 301, processor; 302, memory; 303, communication bus.
具体实施方式Detailed ways
下面将结合本申请实施例中附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围,而是仅仅表示本申请的选定实施例。基于本申请的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本申请保护的范围。The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. The components of the embodiments of the present application generally described and illustrated in the drawings herein may be arranged and designed in a variety of different configurations. Thus, the following detailed description of the embodiments of the application provided in the accompanying drawings is not intended to limit the scope of the application as claimed, but is merely representative of selected embodiments of the application. Based on the embodiments of the present application, all other embodiments obtained by those skilled in the art without creative work fall within the protection scope of the present application.
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。同时,在本申请的描述中,术语“第一”、“第二”等仅用于区分描述,而不能理解为指示或暗示相对重要性。It should be noted that like numerals and letters refer to like items in the following figures, so once an item is defined in one figure, it does not require further definition and explanation in subsequent figures. Meanwhile, in the description of the present application, the terms "first", "second", etc. are only used to distinguish the description, and cannot be understood as indicating or implying relative importance.
第一方面,请参照图1-图3,图1是本申请一些实施例中的一种薄膜沉积模拟方法,用于建立数字化薄膜模型以模拟薄膜沉积过程,方法包括以下步骤:In the first aspect, please refer to FIG. 1 to FIG. 3 . FIG. 1 is a thin film deposition simulation method in some embodiments of the present application, which is used to establish a digital thin film model to simulate a thin film deposition process. The method includes the following steps:
S1、设定代表薄膜凝固情况的第一相场参量及代表输入蒸汽局部密度的第二相场参量;S1. Set the first phase field parameter representing the solidification of the film and the second phase field parameter representing the local density of the input steam;
具体地,由于本申请实施例的薄膜沉积模拟方法主要用于模拟PVD(物理气相沉积)形成薄膜的过程,故其需要模拟气相的目标材料沉积形成固相的目标材料的过程,因此,需要设计表征薄膜关于不同时间状态及空间位置的凝固情况的相场参量,即设计能反映薄膜沉积模拟过程中不同位置在不同时间上是否已形成薄膜的相场参量,将该关于是否已形成薄膜的相场参量定义为第一相场参量,同时,还需要设计表征关于不同空间位置及时间状态的气体密度情况的相场参量,该气体密度的相场参量能反映薄膜沉积的趋势,将该关于气体密度情况的相场参量定义为第二相场参量。Specifically, since the thin film deposition simulation method of the embodiment of the present application is mainly used to simulate the process of PVD (Physical Vapor Deposition) to form a thin film, it needs to simulate the process of depositing a target material in a gas phase to form a target material in a solid phase. Therefore, it is necessary to design The phase field parameters that characterize the solidification of the film with respect to different time states and spatial positions, that is, the phase field parameters designed to reflect whether a film has been formed at different positions at different times during the film deposition simulation process. The field parameter is defined as the first phase field parameter. At the same time, it is also necessary to design the phase field parameter that characterizes the gas density in different spatial positions and temporal states. The phase field parameter of the gas density can reflect the trend of film deposition. The phase field parameter of the density case is defined as the second phase field parameter.
更具体地,输入蒸汽局部密度为考虑对应区域中是否包含固相薄膜的气相的目标材料密度,因此,当第二相场参量为0时能表征对应的区域无气相。More specifically, the input vapor local density is the target material density considering whether the corresponding region contains the gas phase of the solid phase thin film. Therefore, when the second phase field parameter is 0, it can indicate that the corresponding region has no gas phase.
更具体地,第一相场参量能反映出某空间位置在某时间中属于凝固的薄膜区域、固气界面、真空区域、非凝固区域中的何种状态。More specifically, the first phase field parameter can reflect which state a certain spatial position belongs to in a solidified thin film region, a solid-gas interface, a vacuum region, and a non-solidified region at a certain time.
更具体地,第二相场参量反映出某空间位置在某时间中气相的目标材料的密度。More specifically, the second phase field parameter reflects the density of the target material in the gas phase at a certain spatial location at a certain time.
S2、根据第一相场参量和第二相场参量建立用于模拟物理气相沉积的弹道沉积相场模型;S2, establishing a ballistic deposition phase field model for simulating physical vapor deposition according to the first phase field parameter and the second phase field parameter;
具体地,物理气相沉积过程为气相的目标材料沉积形成固相的目标材料的过程,设定两个相场参量分别表征固相和气相情况,可有效简化整个模型的建立条件,即结合能表征薄膜凝固情况的第一相场参量及代表输入蒸汽局部密度的第二相场参量能便捷地建立关于气相目标材料沉积转换为固相目标材料过程的弹道沉积相场模型,在本申请实施例中,根据第一相场参量和第二相场参量建立对应的控制方程便能快速建立弹道沉积相场模型。Specifically, the physical vapor deposition process is the process of depositing the target material in the gas phase to form the target material in the solid phase. Setting two phase field parameters to represent the solid phase and gas phase conditions can effectively simplify the establishment conditions of the entire model, that is, the binding energy characterization The first phase field parameter of the film solidification condition and the second phase field parameter representing the local density of the input steam can easily establish a ballistic deposition phase field model for the process of converting the deposition of a vapor-phase target material into a solid-phase target material. In the embodiments of this application , the ballistic sedimentary phase field model can be quickly established by establishing corresponding control equations according to the first phase field parameters and the second phase field parameters.
S3、根据第一相场参量耦合弹道沉积相场模型和固态薄膜的杨氏模量,以使弹道沉积相场模型与力学平衡方程耦合而建立薄膜沉积模型;S3, coupling the ballistic deposition phase field model and the Young's modulus of the solid thin film according to the first phase field parameter, so as to couple the ballistic deposition phase field model and the mechanical equilibrium equation to establish a thin film deposition model;
具体地,在实际薄膜沉积过程中,沉积获取的薄膜自身具备弹性,单纯考虑薄膜的相变换过程无法精确有效地模拟薄膜沉积过程,因此,在薄膜沉积模拟过程中,需将考虑沉积行程薄膜的力学性能考虑进去以更准确地模拟薄膜沉积过程。Specifically, in the actual film deposition process, the deposited film itself has elasticity, and the film deposition process cannot be accurately and effectively simulated by simply considering the phase transformation process of the film. Therefore, in the film deposition simulation process, it is necessary to consider the deposition process of the film. Mechanical properties are taken into account to more accurately simulate the thin film deposition process.
更具体地,薄膜的力学性能主要体现在弹性形变范围内的应力变化上,因此,在薄膜沉积过程中引入固态薄膜的杨氏模量便能清楚反映薄膜的力学性能。More specifically, the mechanical properties of thin films are mainly reflected in the stress changes in the elastic deformation range. Therefore, the introduction of the Young's modulus of solid thin films during the thin film deposition process can clearly reflect the mechanical properties of thin films.
更具体地,在薄膜沉积过程中,气相目标材料沉积转换成固相目标材料时才具备应力,弹道沉积相场模型为根据薄膜凝固情况的第一相场参量进行建立的,故在本申请实施例中,可根据第一相场参量将固态薄膜的杨氏模量与弹道沉积相场模型进行耦合,使得弹道沉积相场模型关联于薄膜的力学性能,以使弹道沉积相场模型与力学平衡方程耦合而建立薄膜沉积模型。More specifically, in the thin film deposition process, the gas phase target material has stress when it is converted into a solid phase target material. The ballistic deposition phase field model is established according to the first phase field parameters of the film solidification situation, so it is implemented in this application. In an example, the Young's modulus of the solid film can be coupled with the ballistic deposition phase field model according to the first phase field parameter, so that the ballistic deposition phase field model is related to the mechanical properties of the film, so that the ballistic deposition phase field model and the mechanical balance are The equations are coupled to establish the thin film deposition model.
更具体地,薄膜沉积过程形成的薄膜处于稳定状态,即固相的目标材料处于力学平衡状态,故在弹道沉积相场模型与固态薄膜的杨氏模量耦合后,利用力学平衡方程对弹道沉积相场模型进行约束耦合便能获取薄膜沉积模型。More specifically, the thin film formed during the thin film deposition process is in a stable state, that is, the target material in the solid phase is in a mechanical equilibrium state. Therefore, after the ballistic deposition phase field model is coupled with the Young's modulus of the solid thin film, the mechanical equilibrium equation is used for the ballistic deposition. Constrained coupling of the phase field model can obtain a thin film deposition model.
S4、根据薄膜沉积模型进行薄膜沉积模拟。S4, performing thin film deposition simulation according to the thin film deposition model.
具体地,该薄膜沉积模型为基于弹道沉积相场模型与力学平衡方程耦合建立,能在考虑沉积形成的薄膜的应力平衡的情况下模拟物理气相沉积过程,能准确地实现薄膜沉积的预测模拟;由于薄膜沉积模型耦合了固态薄膜的杨氏模量,沉积模拟生成的模拟结果还能反映薄膜沉积形态和应力分布。Specifically, the thin film deposition model is established based on the coupling of the ballistic deposition phase field model and the mechanical equilibrium equation, which can simulate the physical vapor deposition process under the condition of considering the stress balance of the deposited thin film, and can accurately realize the prediction and simulation of thin film deposition; Since the film deposition model is coupled with the Young's modulus of the solid film, the simulation results generated by the deposition simulation can also reflect the deposition morphology and stress distribution of the film.
更具体地,目标材料并不限定为特定现有材料,本申请实施例的薄膜沉积模拟方法旨在通过数字化手段模拟物理气相沉积生产薄膜的过程,以反映相、晶粒和晶界的形成,故可根据已有的目标材料或假设的目标材料的特性调节薄膜沉积模型的设定参数进行薄膜沉积模拟。More specifically, the target material is not limited to a specific existing material, and the thin film deposition simulation method of the embodiment of the present application aims to simulate the process of producing a thin film by physical vapor deposition through digital means, so as to reflect the formation of phases, grains and grain boundaries, Therefore, the thin film deposition simulation can be performed by adjusting the setting parameters of the thin film deposition model according to the properties of the existing target material or the assumed target material.
本申请实施例的薄膜沉积模拟方法基于代表薄膜凝固情况的第一相场参量和固态薄膜的杨氏模量的关系耦合弹道沉积相场模型和固态薄膜的杨氏模量,使最终获取的薄膜沉积模型能反映薄膜的力学性能变化,能将薄膜模拟沉积过程中的相场变化和薄膜力学性能结合模拟,实现薄膜沉积过程的预测模拟,且能根据模拟结果反映薄膜沉积形态和应力分布特点。The thin film deposition simulation method of the embodiments of the present application couples the ballistic deposition phase field model and the Young's modulus of the solid thin film based on the relationship between the first phase field parameter representing the solidification of the thin film and the Young's modulus of the solid thin film, so that the final obtained thin film The deposition model can reflect the change of the mechanical properties of the thin film, and can combine the phase field change and the mechanical properties of the thin film during the simulated deposition process to realize the prediction and simulation of the thin film deposition process, and can reflect the deposition morphology and stress distribution characteristics of the thin film according to the simulation results.
在一些优选的实施方式中,根据第一相场参量耦合弹道沉积相场模型和固态薄膜的杨氏模量,以使弹道沉积相场模型与力学平衡方程耦合而建立薄膜沉积模型的步骤包括:In some preferred embodiments, the step of coupling the ballistic deposition phase field model and the Young's modulus of the solid thin film according to the first phase field parameter, so as to couple the ballistic deposition phase field model with the mechanical equilibrium equation, to establish the thin film deposition model includes:
S31、根据第一相场参量获取弹道沉积相场模型的局部固体密度参量,耦合局部固体密度参量和固态薄膜的杨氏模量,获取关于局部固体密度参量的杨氏模量;S31, obtaining the local solid density parameter of the ballistic deposition phase field model according to the first phase field parameter, coupling the local solid density parameter and the Young's modulus of the solid thin film, and obtaining the Young's modulus about the local solid density parameter;
S32、根据关于局部固体密度参量的杨氏模量和预设的位移场获取应力方程;S32, obtaining the stress equation according to the Young's modulus of the local solid density parameter and the preset displacement field;
具体地,获取局部固体密度参量的杨氏模量能使弹道沉积相场模型与薄膜应力产生关联,为了更准确地反映整个薄膜沉积模拟过程的应变情况,需要引入预设的位移场来评价水平和竖直方向的位移,在本申请实施例中,预设的位移场为,其中u、v分别水平和竖直方向的位移分量。Specifically, obtaining the Young's modulus of the local solid density parameter can correlate the ballistic deposition phase field model with the film stress. In order to more accurately reflect the strain of the entire film deposition simulation process, a preset displacement field needs to be introduced to evaluate the level and displacement in the vertical direction, in this embodiment of the present application, the preset displacement field is , where u and v are the displacement components in the horizontal and vertical directions, respectively.
更具体地,定义模拟过程中薄膜的应变为ε ij ,则有,i和j为水平方向上两个相互垂直的方向坐标标记,u为水平位移分量,故u i,j 和u j,i 为两个相互垂直方向上的水平位移分量,结合该应变和局部固体密度参量的杨氏模量便能获取弹道沉积相场模型生成的薄膜的应力方程。More specifically, defining the strain of the film during the simulation as ε ij , we have , i and j are two mutually perpendicular direction coordinate marks in the horizontal direction, u is the horizontal displacement component, so u i,j and u j,i are the two horizontal displacement components in the two mutually perpendicular directions, combining the strain and local The Young's modulus of the solid density parameter can be used to obtain the stress equation of the film generated by the ballistic deposition phase field model.
S33、根据应力方程和力学平衡方程生成薄膜沉积模型。S33, generating a thin film deposition model according to the stress equation and the mechanical equilibrium equation.
具体地,在薄膜沉积过程中,沉积形成的薄膜应为稳定的固相,故薄膜的应力应该满足力学平衡需求,因此,需要结合应力方程和力学平衡方程生成薄膜沉积模型,使薄膜沉积模型模拟沉积的薄膜满足应力平衡。Specifically, in the film deposition process, the deposited film should be a stable solid phase, so the stress of the film should meet the requirements of mechanical balance. Therefore, it is necessary to combine the stress equation and the mechanical balance equation to generate a film deposition model to simulate the film deposition model. The deposited films satisfy the stress balance.
在一些优选的实施方式中,根据第一相场参量获取弹道沉积相场模型的局部固体密度参量的步骤包括:In some preferred embodiments, the step of obtaining the local solid density parameter of the ballistic deposition phase field model according to the first phase field parameter includes:
定义局部固体密度参量为f d ,局部固体密度参量满足:Define the local solid density parameter as f d , and the local solid density parameter satisfies:
f d =(f+f max )/2 (1) f d = ( f + f max )/2 (1)
f为第一相场参量,f max 为第一相场参量的最大值。 f is the first phase field parameter, and f max is the maximum value of the first phase field parameter.
在一些优选的实施方式中,为了更清楚地描述薄膜的凝固情况,本申请实施例的方法对第一相场参量进行归一化处理,在实际模拟过程中,相场状态具有一定偏差值,因此,f取近似值来表征对应的薄膜的凝固情况,使得f≈1代表凝固的薄膜,f≈0代表固气界面,f≈-1代表真空或非凝固区域,因此有f max =1,f d ∈[0,1],使得f d 能表征对应区域中的固态薄膜的占比,即f d 能充当相场中固相的体积分数。In some preferred embodiments, in order to describe the solidification of the thin film more clearly, the method in the embodiments of the present application normalizes the first phase field parameter. In the actual simulation process, the phase field state has a certain deviation value, Therefore, f is approximated to characterize the solidification of the corresponding film, such that f ≈ 1 represents the solidified film, f ≈ 0 represents the solid-gas interface, and f ≈-1 represents the vacuum or non-solidified region, so there is f max =1, f d ∈ [0,1], so that f d can characterize the proportion of solid thin films in the corresponding region, that is, f d can act as the volume fraction of solid phase in the phase field.
在一些优选的实施方式中,局部固体密度参量的杨氏模量具有残余模量系数。In some preferred embodiments, the Young's modulus of the local solid density parameter has a residual modulus coefficient.
具体地,若直接关联局部固体密度参量和固态薄膜的杨氏模量获取局部固体密度参量的杨氏模量,局部固体密度参量为0(即f d =0)时,会引起局部固体密度参量的杨氏模量为0,从而引起数值奇异性,故本申请实施例的方法,在关联局部固体密度参量和固态薄膜的杨氏模量时设置残余模量系数,以使f d =0时,局部固体密度参量的杨氏模量不为0,在气相模量可以忽略的情况下,本申请实施例的方法获取的薄膜沉积模型可以避免数值奇异性。Specifically, if the Young's modulus of the local solid density parameter is obtained by directly correlating the local solid density parameter and the Young's modulus of the solid thin film, when the local solid density parameter is 0 (that is, f d = 0), the local solid density parameter The Young's modulus of , is 0, which causes numerical singularity. Therefore, in the method of the embodiment of the present application, the residual modulus coefficient is set when correlating the local solid density parameter and the Young's modulus of the solid thin film, so that when f d =0 , the Young's modulus of the local solid density parameter is not 0, and when the gas-phase modulus can be ignored, the thin film deposition model obtained by the method of the embodiment of the present application can avoid numerical singularity.
在一些优选的实施方式中,具有残余模量系数的局部固体密度参量的杨氏模量表示为E,满足:In some preferred embodiments, the Young's modulus of the local solid density parameter with the residual modulus coefficient, denoted E, satisfies:
(2) (2)
其中,h(f d )为相场转换函数,满足h(0)=0及h(1)=1,E 0为固态薄膜的杨氏模量,为残余模量系数,为一个不为零的小数,当f d =0时,有E=E 0 ,能使得E不等于零,从而使得气相的弹性模量不为零,从而避免数值奇异性,确保本申请实施例方法获取的薄膜沉积模型能稳定运行模拟;同理,f d =1时,有E=E 0,即在固相中,局部固体密度参量的杨氏模量与固态薄膜的杨氏模量相等。where h ( f d ) is the phase field transfer function, satisfying h (0)=0 and h (1)=1, E 0 is the Young’s modulus of the solid thin film, is the residual modulus coefficient, is a non-zero decimal, when f d =0, there is E=E 0 , can make E not equal to zero, so that the elastic modulus of the gas phase is not zero, thereby avoiding numerical singularity, and ensuring that the thin film deposition model obtained by the method of the embodiment of the present application can run stably for simulation; Similarly, when f d =1, there are E = E 0 , that is, in the solid phase, the Young's modulus of the local solid density parameter is equal to that of the solid thin film.
更优选地,在本申请实施例中,h(f d )满足:More preferably, in this embodiment of the present application, h ( f d ) satisfies:
(3) (3)
具体地,式(3)满足h(0)=0及h(1)=1,能将局部固体密度参量转换为与固态薄膜的杨氏模量关联的约束变量。Specifically, Equation (3) satisfies h (0)=0 and h (1)=1, which can convert the local solid density parameter into a constraint variable associated with the Young's modulus of the solid thin film.
基于上述对杨氏模量的处理,能使气相中的杨氏模量可以忽略,而固态薄膜中的杨氏模量为E 0,固气界面上的杨氏模量介于二者之间。值得一提的是,气相中的杨氏模量虽不为零,但没有实际物理意义,并且对固相中应力分布的影响可以忽略不计,增加残余模量系数目的仅为避免数值奇异性。Based on the above treatment of Young's modulus, the Young's modulus in the gas phase can be ignored, while the Young's modulus in the solid film is E 0 , and the Young's modulus at the solid-gas interface is between the two. . It is worth mentioning that although the Young's modulus in the gas phase is not zero, it has no practical physical meaning and has a negligible effect on the stress distribution in the solid phase. The purpose of increasing the residual modulus coefficient is only to avoid numerical singularities.
更具体地,上述公式还忽略了固相到气相过渡时的泊松比区别,有效简化模型模拟逻辑。More specifically, the above formula also ignores the difference of Poisson's ratio in the transition from solid phase to gas phase, which effectively simplifies the model simulation logic.
在一些优选的实施方式中,步骤S32获取的应力方程为:In some preferred embodiments, the stress equation obtained in step S32 is:
(4) (4)
其中,σ ij 为坐标[i,j]上的应力,υ为泊松比,δ ij 为克罗内克尔(kronecker)符号,ε kk 为关联于δ ij 的对称应变,因此,步骤S33中的力学平衡方程为,为薄膜整体应力,即σ ij 属于薄膜整体中[i,j]处的应力,利用该力学平衡方程约束应力方程应力分布能使弹道沉积相场模型转变为薄膜沉积模型。Among them, σ ij is the stress on the coordinates [ i, j ], υ is the Poisson's ratio, δ ij is the Kronecker symbol, ε kk is the symmetric strain associated with δ ij , therefore, in step S33 The mechanical equilibrium equation of is , is the overall stress of the film, that is, σ ij belongs to the stress at [ i, j ] in the whole film. Using the mechanical equilibrium equation to constrain the stress distribution of the stress equation can transform the ballistic deposition phase field model into a film deposition model.
在一些优选的实施方式中,根据薄膜沉积模型进行薄膜沉积模拟的步骤包括:In some preferred embodiments, the step of simulating thin film deposition according to the thin film deposition model includes:
S41、根据反应条件及薄膜材质设定弹道沉积相场模型的现象学参数及固态薄膜的杨氏模量和泊松比;S41. Set the phenomenological parameters of the ballistic deposition phase field model and the Young's modulus and Poisson's ratio of the solid film according to the reaction conditions and the film material;
具体地,通过调节现象学参数能改变薄膜沉积模型的弹道沉积工况和效率,通过改变固态薄膜的杨氏模量和泊松比能改变薄膜沉积模型沉积形成的薄膜的力学性能,通过调整上述数据能使薄膜沉积模型模拟出不同材质在不同状态环境下的沉积过程,使得本申请实施例的方法能根据模拟需求构建不同材料特性的薄膜沉积模型。Specifically, by adjusting the phenomenological parameters, the ballistic deposition conditions and efficiency of the thin film deposition model can be changed, and by changing the Young's modulus and Poisson's ratio of the solid thin film, the mechanical properties of the thin film deposited by the thin film deposition model can be changed. The thin film deposition model can simulate the deposition process of different materials in different state environments, so that the method of the embodiments of the present application can construct thin film deposition models with different material properties according to simulation requirements.
S42、将弹道沉积相场模型的现象学参数及固态薄膜的杨氏模量和泊松比输入至薄膜沉积模型中,并利用薄膜沉积模型模拟薄膜沉积生成薄膜图像信息。S42 , input the phenomenological parameters of the ballistic deposition phase field model and the Young's modulus and Poisson's ratio of the solid thin film into the thin film deposition model, and use the thin film deposition model to simulate thin film deposition to generate thin film image information.
具体地,利用该沉积模型生成薄膜图像信息能直观地显示薄膜演变过程,便于用户了解材料特性与薄膜形成过程的关系,以进行实验或生产的控制调整。Specifically, using the deposition model to generate film image information can intuitively display the film evolution process, which is convenient for users to understand the relationship between the material properties and the film formation process, so as to control and adjust the experiment or production.
在一些优选的实施方式中,根据薄膜沉积模型进行薄膜沉积模拟的步骤还包括:In some preferred embodiments, the step of performing the thin film deposition simulation according to the thin film deposition model further comprises:
利用薄膜沉积模型模拟薄膜沉积生成应力分布图像信息。The stress distribution image information is generated by simulating thin film deposition using the thin film deposition model.
具体地,由于本申请实施例的方法中的弹道沉积相场模型耦合了固态薄膜的杨氏模量,使得薄膜沉积模型模拟沉积生成的薄膜包含了内部应力数据,故能生成应力分布图像信息,以供用户了解薄膜形成过程中的应力分布特点,以进行实验或生产的控制调整。Specifically, since the ballistic deposition phase field model in the method of the embodiment of the present application is coupled with the Young's modulus of the solid thin film, the thin film generated by the simulated deposition of the thin film deposition model contains internal stress data, so the stress distribution image information can be generated, For users to understand the characteristics of stress distribution in the process of film formation, in order to carry out experiments or production control adjustment.
在一些优选的实施方式中,弹道沉积相场模型包括关于第一相场参量的薄膜相控制方程和关于第二相场参量的蒸汽控制方程。In some preferred embodiments, the ballistic deposition phase field model includes a thin film phase governing equation for a first phase field parameter and a vapor governing equation for a second phase field parameter.
具体地,设计薄膜相控制方程和蒸汽控制方程时应该参考弹道沉积形态进行设计对应相场的无量纲控制方程,在本申请实施例中,薄膜相控制方程为:Specifically, when designing the film phase control equation and the steam control equation, the dimensionless control equation corresponding to the phase field should be designed with reference to the ballistic deposition form. In the embodiment of the present application, the film phase control equation is:
(5) (5)
蒸汽控制方程为:The steam governing equation is:
(6) (6)
其中,为梯度算符,为二阶偏导算符,g为第二相场参量,满足g≥0,g=0表示该区域无气相,为空间位置,t为时间,a、B、D和为现象学参数,a代表生成薄膜的表面张力,B代表输入的粒子粘在薄膜表面的可能性,D为输入粒子在气相中的扩散常数,为输入粒子的垂直偏置通量强度,垂直弹道沉积时有,A为标量,代表通量强度的值,为坐标轴单位微量。in, is the gradient operator, is the second-order partial derivative operator, g is the second phase field parameter, satisfying g ≥ 0, g = 0 means that there is no gas phase in the region, is the spatial position, t is the time, a, B, D and is the phenomenological parameter, a represents the surface tension of the generated film, B represents the possibility of the input particles sticking to the surface of the film, D is the diffusion constant of the input particles in the gas phase, is the vertical bias flux intensity of the input particles, which has , A is a scalar, representing the value of flux intensity, It is a micro amount for the coordinate axis unit.
具体地,在本申请实施例中,由式(2)、(4)和力学平衡方程约束的式(5)和(6)构成了薄膜沉积模型,即本申请实施例的方法为在力学平衡方程约束下利用式(5)和(6)进行薄膜沉积模拟,求解力学平衡方程能获取薄膜的位移和应力分布。Specifically, in the embodiments of the present application, the equations (5) and (6) constrained by the equations (2), (4) and the mechanical equilibrium equation constitute the thin film deposition model, that is, the method in the embodiments of the present application is a mechanical equilibrium Under the constraints of equations, equations (5) and (6) are used to simulate the deposition of thin films, and the displacement and stress distribution of the thin films can be obtained by solving the mechanical equilibrium equation.
更具体地,根据步骤S41可知,在利用该薄膜沉积模型进行沉积模拟前,需根据材料特性设定现象学参数及固态薄膜的杨氏模量和泊松比,以实现对应目标材料的薄膜沉积模拟,如图2(a)-图2(d)所示为某目标材料的薄膜沉积模拟过程生成的薄膜图像信息,其中,图中横坐标为薄膜沉积的水平距离,纵坐标为薄膜沉积的竖直距离,右侧颜色条为局部固体密度参量f d 与颜色的关系,a=0.3、B=10、D=0.01和A=0.6、E=100GPa、υ=0.3,并采用垂直弹道进行沉积,基底位置的f初始值为1,气相中f初始值为-1,左右边界上施加周期性边界条件,而在上下边界上施加反周期性边界条件,由此,利用公式(5)和(6)进行有限元计算模拟获取如图2(a)-图2(d)所示的薄膜图像信息,能清楚反映出沉积薄膜从基底逐渐生长的过程,并具有柱状晶特征;如图3所示为图2(d)薄膜状态对应的应力分布图像信息,其中,图中横坐标为薄膜沉积的水平距离,纵坐标为薄膜沉积的竖直距离,右侧颜色条为应力大小与颜色的关系,可见图2(d)中的薄膜表面的不均匀导致图3中水平正应力分布的不均匀,可见本申请实施例的方法实现了薄膜形态和应力分布的直观演示。More specifically, according to step S41, before using the film deposition model for deposition simulation, it is necessary to set phenomenological parameters, Young's modulus and Poisson's ratio of the solid film according to the material properties, so as to realize the film deposition simulation corresponding to the target material. , as shown in Figure 2(a)-Figure 2(d), the film image information generated by the film deposition simulation process of a target material, where the abscissa in the figure is the horizontal distance of the film deposition, and the ordinate is the vertical film deposition. Straight distance, the color bar on the right is the relationship between the local solid density parameter f d and the color, a= 0.3 , B= 10 , D= 0.01 and A= 0.6, E =100GPa, υ= 0.3, and the vertical ballistic trajectory is used for deposition, The initial value of f at the substrate position is 1, and the initial value of f in the gas phase is -1. Periodic boundary conditions are imposed on the left and right boundaries, and anti-periodic boundary conditions are imposed on the upper and lower boundaries. Therefore, using formulas (5) and (6) ) to perform finite element calculation and simulation to obtain the film image information shown in Figure 2(a)-Figure 2(d), which can clearly reflect the gradual growth process of the deposited film from the substrate, and has columnar crystal characteristics; as shown in Figure 3 is the stress distribution image information corresponding to the thin film state in Fig. 2(d), in which the abscissa in the figure is the horizontal distance of film deposition, the ordinate is the vertical distance of thin film deposition, and the color bar on the right is the relationship between stress and color, It can be seen that the unevenness of the film surface in Fig. 2(d) leads to the uneven distribution of the horizontal normal stress in Fig. 3, and it can be seen that the method of the embodiment of the present application realizes an intuitive demonstration of the film morphology and stress distribution.
第二方面,请参照图4,图4是本申请一些实施例中提供的一种薄膜沉积模拟装置,用于建立数字化薄膜模型以模拟薄膜沉积过程,装置包括:In the second aspect, please refer to FIG. 4. FIG. 4 is a thin film deposition simulation device provided in some embodiments of the present application for establishing a digital thin film model to simulate a thin film deposition process. The device includes:
设定模块201,用于设定代表薄膜凝固情况的第一相场参量及代表输入蒸汽局部密度的第二相场参量;The
建模模块202, 用于根据第一相场参量和第二相场参量建立用于模拟物理气相沉积的弹道沉积相场模型;a modeling module 202, configured to establish a ballistic deposition phase field model for simulating physical vapor deposition according to the first phase field parameter and the second phase field parameter;
耦合模块203,用于根据第一相场参量耦合弹道沉积相场模型和固态薄膜的杨氏模量,以使弹道沉积相场模型与力学平衡方程耦合而建立薄膜沉积模型;The
模拟模块204,用于根据薄膜沉积模型进行薄膜沉积模拟。The simulation module 204 is used for performing thin film deposition simulation according to the thin film deposition model.
本申请实施例的薄膜沉积模拟装置基于代表薄膜凝固情况的第一相场参量和固态薄膜的杨氏模量的关系耦合弹道沉积相场模型和固态薄膜的杨氏模量,使最终获取的薄膜沉积模型能反映薄膜的力学性能变化,能将薄膜模拟沉积过程中的相场变化和薄膜力学性能结合模拟,实现薄膜沉积过程的预测模拟,且能根据模拟结果反映薄膜沉积形态和应力分布特点。The thin film deposition simulation device of the embodiment of the present application couples the ballistic deposition phase field model and the Young's modulus of the solid thin film based on the relationship between the first phase field parameter representing the solidification of the thin film and the Young's modulus of the solid thin film, so that the final obtained thin film The deposition model can reflect the change of the mechanical properties of the thin film, and can combine the phase field change and the mechanical properties of the thin film during the simulated deposition process to realize the prediction and simulation of the thin film deposition process, and can reflect the deposition morphology and stress distribution characteristics of the thin film according to the simulation results.
在一些优选的实施方式中,本申请实施例的薄膜沉积模拟装置用于执行上述第一方面提供的薄膜沉积模拟方法。In some preferred embodiments, the thin film deposition simulation apparatus of the embodiments of the present application is used to execute the thin film deposition simulation method provided in the first aspect above.
第三方面,请参照图5,图5为本申请实施例提供的一种电子设备的结构示意图,本申请提供一种电子设备,包括:处理器301和存储器302,处理器301和存储器302通过通信总线303和/或其他形式的连接机构(未标出)互连并相互通讯,存储器302存储有处理器301可执行的计算机程序,当计算设备运行时,处理器301执行该计算机程序,以执行时执行上述实施例的任一可选的实现方式中的方法。In the third aspect, please refer to FIG. 5 , which is a schematic structural diagram of an electronic device provided by an embodiment of the present application. The present application provides an electronic device including: a
第四方面,本申请实施例提供一种存储介质,其上存储有计算机程序,该计算机程序被处理器执行时,执行上述实施例的任一可选的实现方式中的方法。其中,存储介质可以由任何类型的易失性或非易失性存储设备或者它们的组合实现,如静态随机存取存储器(Static Random Access Memory, 简称SRAM),电可擦除可编程只读存储器(ElectricallyErasable Programmable Read-Only Memory, 简称EEPROM),可擦除可编程只读存储器(Erasable Programmable Read Only Memory, 简称EPROM),可编程只读存储器(Programmable Red-Only Memory, 简称PROM),只读存储器(Read-Only Memory, 简称ROM),磁存储器,快闪存储器,磁盘或光盘。In a fourth aspect, an embodiment of the present application provides a storage medium on which a computer program is stored, and when the computer program is executed by a processor, executes the method in any optional implementation manner of the foregoing embodiment. Among them, the storage medium can be realized by any type of volatile or non-volatile storage device or their combination, such as Static Random Access Memory (SRAM for short), Electrically Erasable Programmable Read-Only Memory (Electrically Erasable Programmable Read-Only Memory, EEPROM for short), Erasable Programmable Read Only Memory (EPROM), Programmable Red-Only Memory (PROM), Read Only Memory (Read-Only Memory, referred to as ROM), magnetic memory, flash memory, magnetic disk or optical disk.
综上,本申请实施例提供了一种薄膜沉积模拟方法、装置、电子设备及存储介质,其中,方法基于代表薄膜凝固情况的第一相场参量和固态薄膜的杨氏模量的关系耦合弹道沉积相场模型和固态薄膜的杨氏模量,使最终获取的薄膜沉积模型能反映薄膜的力学性能变化,能将薄膜模拟沉积过程中的相场变化和薄膜力学性能结合模拟,实现薄膜沉积过程的预测模拟,且能根据模拟结果反映薄膜沉积形态和应力分布特点。To sum up, the embodiments of the present application provide a thin film deposition simulation method, device, electronic device, and storage medium, wherein the method is based on the relationship between the first phase field parameter representing the solidification of the thin film and the Young's modulus of the solid thin film coupled with ballistic trajectory The deposition phase field model and the Young's modulus of the solid-state film enable the final film deposition model to reflect the mechanical properties of the film. It can combine the phase field changes during the film deposition process with the mechanical properties of the film to simulate the film deposition process. The prediction and simulation of the film can reflect the deposition morphology and stress distribution characteristics of the film according to the simulation results.
在本申请所提供的实施例中,应该理解到,所揭露装置和方法,可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,又例如,多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。In the embodiments provided in this application, it should be understood that the disclosed apparatus and method may be implemented in other manners. The apparatus embodiments described above are only illustrative. For example, the division of the units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components may be combined or Can be integrated into another system, or some features can be ignored, or not implemented. On the other hand, the shown or discussed mutual coupling or direct coupling or communication connection may be through some communication interfaces, indirect coupling or communication connection of devices or units, which may be in electrical, mechanical or other forms.
另外,作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。In addition, units described as separate components may or may not be physically separated, and components shown as units may or may not be physical units, that is, may be located in one place, or may be distributed to multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution in this embodiment.
再者,在本申请各个实施例中的各功能模块可以集成在一起形成一个独立的部分,也可以是各个模块单独存在,也可以两个或两个以上模块集成形成一个独立的部分。Furthermore, each functional module in each embodiment of the present application may be integrated together to form an independent part, or each module may exist alone, or two or more modules may be integrated to form an independent part.
在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。In this document, relational terms such as first and second, etc. are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such existence between these entities or operations. The actual relationship or sequence.
以上所述仅为本申请的实施例而已,并不用于限制本申请的保护范围,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。The above descriptions are merely examples of the present application, and are not intended to limit the protection scope of the present application. For those skilled in the art, various modifications and changes may be made to the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of this application shall be included within the protection scope of this application.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210548370.3A CN114662424B (en) | 2022-05-20 | 2022-05-20 | Thin film deposition simulation method, device, electronic device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210548370.3A CN114662424B (en) | 2022-05-20 | 2022-05-20 | Thin film deposition simulation method, device, electronic device and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114662424A true CN114662424A (en) | 2022-06-24 |
| CN114662424B CN114662424B (en) | 2022-09-02 |
Family
ID=82037854
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210548370.3A Active CN114662424B (en) | 2022-05-20 | 2022-05-20 | Thin film deposition simulation method, device, electronic device and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114662424B (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070017296A1 (en) * | 2005-05-10 | 2007-01-25 | Instrument Technology Research Center, National Applied Research Laboratories | Method for whole field thin film stress evaluation |
| US20140002203A1 (en) * | 2012-06-27 | 2014-01-02 | US Gov't Represented by Secretary of the Navy Chief of Naval Research Office of Counsel ONR/NRL | Method and system of an ultra high q silicon cantilever resonator for thin film internal friction and young's modulus measurements |
| US9053929B1 (en) * | 2010-05-19 | 2015-06-09 | Corporation For National Research Initiatives | Method and system for integrated MEMS and NEMS using deposited thin films having pre-determined stress states |
| US20160162616A1 (en) * | 2014-03-27 | 2016-06-09 | King Fahd University Of Petroleum And Minerals | Performance and life prediction model for photovoltaic module: effect of encapsulant constitutive behavior |
| CN110866349A (en) * | 2019-12-03 | 2020-03-06 | 湘潭大学 | Phase field analysis method of hafnium oxide based ferroelectric film based on polymorphic coexistence |
| US20200140993A1 (en) * | 2018-11-02 | 2020-05-07 | Ivworks Co., Ltd. | Apparatus, method, and recording medium storing command for controlling thin-film deposition process |
| CN112307661A (en) * | 2020-10-30 | 2021-02-02 | 湘潭大学 | Analysis of Domain and Domain Wall Motion in Hafnium Oxide-Based Ferroelectric Thin Films under Ultrafast Electrical Pulses |
| US20210118530A1 (en) * | 2019-05-27 | 2021-04-22 | Beijing University Of Technology | Multi-scale method for simulating mechanical behaviors of multiphase composite materials |
| CN114112145A (en) * | 2021-11-10 | 2022-03-01 | 中国科学院半导体研究所 | Thin film intrinsic stress measurement method, electronic device and medium |
-
2022
- 2022-05-20 CN CN202210548370.3A patent/CN114662424B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070017296A1 (en) * | 2005-05-10 | 2007-01-25 | Instrument Technology Research Center, National Applied Research Laboratories | Method for whole field thin film stress evaluation |
| US9053929B1 (en) * | 2010-05-19 | 2015-06-09 | Corporation For National Research Initiatives | Method and system for integrated MEMS and NEMS using deposited thin films having pre-determined stress states |
| US20140002203A1 (en) * | 2012-06-27 | 2014-01-02 | US Gov't Represented by Secretary of the Navy Chief of Naval Research Office of Counsel ONR/NRL | Method and system of an ultra high q silicon cantilever resonator for thin film internal friction and young's modulus measurements |
| US20160162616A1 (en) * | 2014-03-27 | 2016-06-09 | King Fahd University Of Petroleum And Minerals | Performance and life prediction model for photovoltaic module: effect of encapsulant constitutive behavior |
| US20200140993A1 (en) * | 2018-11-02 | 2020-05-07 | Ivworks Co., Ltd. | Apparatus, method, and recording medium storing command for controlling thin-film deposition process |
| US20210118530A1 (en) * | 2019-05-27 | 2021-04-22 | Beijing University Of Technology | Multi-scale method for simulating mechanical behaviors of multiphase composite materials |
| CN110866349A (en) * | 2019-12-03 | 2020-03-06 | 湘潭大学 | Phase field analysis method of hafnium oxide based ferroelectric film based on polymorphic coexistence |
| CN112307661A (en) * | 2020-10-30 | 2021-02-02 | 湘潭大学 | Analysis of Domain and Domain Wall Motion in Hafnium Oxide-Based Ferroelectric Thin Films under Ultrafast Electrical Pulses |
| CN114112145A (en) * | 2021-11-10 | 2022-03-01 | 中国科学院半导体研究所 | Thin film intrinsic stress measurement method, electronic device and medium |
Non-Patent Citations (1)
| Title |
|---|
| 袁训锋等: "强制对流作用下多晶粒生长的相场模拟", 《材料导报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114662424B (en) | 2022-09-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Seghir et al. | A hybrid approach using genetic and fruit fly optimization algorithms for QoS-aware cloud service composition | |
| Kwon et al. | Run-to-run-based model predictive control of protein crystal shape in batch crystallization | |
| Chen et al. | Spectra and Mach number scaling in compressible homogeneous shear turbulence | |
| CN110032756B (en) | Numerical Analysis Method of Flow Boundary Layer Based on Flow Function Fractional Coordinate System Transformation | |
| CN108764577A (en) | Online time series prediction technique based on dynamic fuzzy Cognitive Map | |
| WO2021135719A1 (en) | Method and apparatus for deducing physical property parameter | |
| Wang | On extended progressive and iterative approximation for least squares fitting | |
| IE20080002A1 (en) | System and method for dynamic generation of environmental operational models | |
| CN111273553A (en) | Coupling identification and decoupling control design method of two-degree-of-freedom servo holder system | |
| Gao et al. | Numerical solution of nonlinear Burgers’ equation using high accuracy multi-quadric quasi-interpolation | |
| CN108846187A (en) | Integrated circuit global wiring optimization method based on Generalized Extended Lagrange | |
| Wu et al. | An analysis of the convection–diffusion problems using meshless and meshbased methods | |
| CN118484994A (en) | A method for predicting production of refractured wells based on convolutional self-attention Transformer | |
| CN111563311A (en) | ZEMAX-based arbitrary-shape grating simulation system | |
| CN111104753A (en) | Viscous incompressible fluid simulation method based on SPH | |
| CN114700953A (en) | A particle swarm hand-eye calibration method and system based on joint zero error | |
| CN114662424B (en) | Thin film deposition simulation method, device, electronic device and storage medium | |
| Dong et al. | Intelligent injection molding: parameters self‐learning optimization using iterative gradient‐approximation adaptive method | |
| CN108958029A (en) | For parameter tuning method and system | |
| Li et al. | Coupled MLPG–FVM simulation of steady state heat conduction in irregular geometry | |
| Petersson et al. | Optimization based LPV-approximation of multi-model systems | |
| CN114580627B (en) | FPGA implementation method based on variable parameter recurrent neural network | |
| CN115906467A (en) | Data processing method and device based on battery swapping station, electronic equipment and storage medium | |
| Yu et al. | Connecting factorization and distance metric learning for social recommendations | |
| CN115495854A (en) | Parameter calibration method, device, equipment and medium of computer aided engineering model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |