载流子输运仿真方法、装置、介质及电子设备

技术领域

本发明属于量子仿真计算技术领域，特别是载流子输运仿真方法、装置、介质及电子设备。

背景技术

现代集成电路作为电子信息设备的核心取得了广泛的应用，其稳定性一直得到各界关注。随着现代集成电路工艺特征尺寸的不断缩小，量子效应对集成电路中包括的半导体器件电学性能的影响越来越显著。电学性能由半导体器件中载流子(物理学中将电子和空穴统称为载流子)的输运决定，对半导体器件中载流子输运的理解程度决定了未来集成电路的发展。目前对半导体器件中载流子输运的研究缺乏一个行之有效的方法。

发明内容

本发明的目的是提供一种载流子输运仿真方法、装置、介质及载流子设备，旨在加强对半导体器件中载流子输运的研究。

本申请的一个实施例提供了一种载流子输运仿真方法，所述方法包括：

确定半导体器件中载流子输运的物理仿真模型、初始条件和/或边界条件；

确定求解所述物理仿真模型对应的泛定方程；

基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度，以实现所述半导体器件中载流子输运的仿真。

可选的，所述物理仿真模型包括以下至少一种：半经典模型、封闭量子模型、开放量子模型；所述半经典仿真模型对应的泛定方程为泊松方程；所述封闭量子模型对应的泛定方程为所述泊松方程和薛定谔方程，所述开放量子模型对应的泛定方程为所述泊松方程和所述薛定谔方程对应的格林函数方程。

可选的，所述确定半导体器件中载流子输运的初始条件和/或边界条件，包括：

构建半导体器件的几何模型；

基于有限体积法将所述几何模型网格化，得到多个第一控制体积；

确定每个所述第一控制体积的初始条件和/或边界条件。

可选的，若所述物理仿真模型为所述半经典模型，所述基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度，包括：

基于每个所述第一控制体积的初始条件和/或边界条件、所述泊松方程确定数值矩阵形式的解析方程组；

对所述数值矩阵形式的解析方程组进行变换，得到残差形式的解析方程组；

基于牛顿迭代算法对所述残差形式的方程组进行求解，得到第一静电势；

基于所述第一静电势确定所述半导体器件中的载流子密度。

可选的，若所述物理仿真模型为所述封闭量子模型，所述基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度，包括：

基于每个所述第一控制体积的初始条件和边界条件、所述泊松方程确定数值矩阵形式的第一解析方程组；

确定每个所述第一控制体积对应的所述薛定谔方程，得到数值矩阵形式的第二解析方程组；

确定初始载流子密度；

基于所述初始载流子密度和所述第一解析方程组对所述第二解析方程组迭代求解，确定所述半导体器件中的载流子密度。

可选的，若所述物理仿真模型为所述开放量子模型，所述基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度，包括：

基于每个所述第一控制体积的初始条件和边界条件、所述泊松方程确定数值矩阵形式的第三解析方程组；

确定每个所述第一控制体积对应的所述格林函数方程，得到数值矩阵形式的第四解析方程组；

确定初始静电势；

基于所述初始静电势和所述第三解析方程组对所述第四解析方程组迭代求解，确定所述半导体器件中的载流子密度。

可选的，所述确定所述半导体器件中的载流子密度之后，所述方法还包括：

基于所述半导体器件中的载流子密度确定所述几何模型中载流子密度大于或等于预设阈值的目标区域；

基于有限体积法将所述目标区域网格化，得到多个第二控制体积，所述第二控制体积小于所述第一控制体积；

基于多个所述第二控制体积确定新的所述半导体器件中的载流子密度。

本申请的又一实施例提供了一种载流子输运仿真装置，所述装置包括：

第一确定单元，用于确定半导体器件中载流子输运的物理仿真模型、初始条件和/或边界条件；

第二确定单元，用于确定求解所述物理仿真模型对应的泛定方程；

第三确定单元，用于基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度，以实现所述半导体器件中载流子输运的仿真。

本申请的又一实施例提供了一种存储介质，所述存储介质中存储有计算机程序，其中，所述计算机程序被设置为运行时执行上述任一项中所述的方法。

本申请的又一实施例提供了一种电子装置，包括存储器和处理器，所述存储器中存储有计算机程序，所述处理器被设置为运行所述计算机程序以执行上述任一项中所述的方法。

与现有技术相比，本发明提供的一种载流子输运仿真方法，将不同的物理仿真模型对应的泛定方程与半导体器件中载流子输运的初始条件和/或边界条件结合，从而确定该半导体器件中的载流子密度，以实现半导体器件中载流子输运的仿真，进而实现了对半导体器件中载流子输运的研究。

附图说明

图1为本发明实施例提供的一种载流子输运仿真方法的计算机终端的硬件结构框图；

图2为本发明实施例提供的一种载流子输运仿真方法的流程示意图；

图3为本发明实施例提供的一种第一控制体积的结构示意图；

图4为本发明实施例提供的一种载流子输运仿真装置的结构示意图。

具体实施方式

下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能解释为对本发明的限制。

本发明实施例首先提供了一种载流子输运仿真方法，该方法可以应用于电子设备，如计算机终端，具体如普通电脑、量子计算机等。

下面以运行在计算机终端上为例对其进行详细说明。图1为本发明实施例提供的一种载流子输运仿真方法的计算机终端的硬件结构框图。如图1所示，计算机终端可以包括一个或多个(图1中仅示出一个)处理器102(处理器102可以包括但不限于微处理器MCU或可编程逻辑器件FPGA等的处理装置)和用于存储基于载流子输运仿真方法的存储器104，可选地，上述计算机终端还可以包括用于通信功能的传输装置106以及输入输出设备108。本领域普通技术人员可以理解，图1所示的结构仅为示意，其并不对上述计算机终端的结构造成限定。例如，计算机终端还可包括比图1中所示更多或者更少的组件，或者具有与图1所示不同的配置。

存储器104可用于存储应用软件的软件程序以及模块，如本申请实施例中的分子的参数配置方法对应的程序指令/模块，处理器102通过运行存储在存储器104内的软件程序以及模块，从而执行各种功能应用以及数据处理，即实现上述的方法。存储器104可包括高速随机存储器，还可包括非易失性存储器，如一个或者多个磁性存储装置、闪存、或者其他非易失性固态存储器。在一些实例中，存储器104可进一步包括相对于处理器102远程设置的存储器，这些远程存储器可以通过网络连接至计算机终端。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。

传输装置106用于经由一个网络接收或者发送数据。上述的网络具体实例可包括计算机终端的通信供应商提供的无线网络。在一个实例中，传输装置106包括一个网络适配器(Network Interface Controller，NIC)，其可通过基站与其他网络设备相连从而可与互联网进行通讯。在一个实例中，传输装置106可以为射频(Radio Frequency，RF)模块，其用于通过无线方式与互联网进行通讯。

需要说明的是，真正的量子计算机是混合结构的，它包含两大部分：一部分是经典计算机，负责执行经典计算与控制；另一部分是量子设备，负责运行量子程序进而实现量子计算。而量子程序是由量子语言如QRunes语言编写的一串能够在量子计算机上运行的指令序列，实现了对量子逻辑门操作的支持，并最终实现量子计算。具体的说，量子程序就是一系列按照一定时序操作量子逻辑门的指令序列。

在实际应用中，因受限于量子设备硬件的发展，通常需要进行量子计算模拟以验证量子算法、量子应用等等。量子计算模拟即借助普通计算机的资源搭建的虚拟架构(即量子虚拟机)实现特定问题对应的量子程序的模拟运行的过程。通常，需要构建特定问题对应的量子程序。本发明实施例所指量子程序，即是经典语言编写的表征量子比特及其演化的程序，其中与量子计算相关的量子比特、量子逻辑门等等均有相应的经典代码表示。

量子线路作为量子程序的一种体现方式，也称量子逻辑电路，是最常用的通用量子计算模型，表示在抽象概念下对于量子比特进行操作的线路，其组成包括量子比特、线路(时间线)、以及各种量子逻辑门，最后常需要通过量子测量操作将结果读取出来。

不同于传统电路是用金属线所连接以传递电压信号或电流信号，在量子线路中，线路可看成是由时间所连接，亦即量子比特的状态随着时间自然演化，在这过程中按照哈密顿运算符的指示，一直到遇上逻辑门而被操作。

一个量子程序整体上对应有一条总的量子线路，本发明所述量子程序即指该条总的量子线路，其中，该总的量子线路中的量子比特总数与量子程序的量子比特总数相同。可以理解为：一个量子程序可以由量子线路、针对量子线路中量子比特的测量操作、保存测量结果的寄存器及控制流节点(跳转指令)组成，一条量子线路可以包含几十上百个甚至成千上万个量子逻辑门操作。量子程序的执行过程，就是对所有的量子逻辑门按照一定时序执行的过程。需要说明的是，时序即单个量子逻辑门被执行的时间顺序。

需要说明的是，经典计算中，最基本的单元是比特，而最基本的控制模式是逻辑门，可以通过逻辑门的组合来达到控制电路的目的。类似地，处理量子比特的方式就是量子逻辑门。使用量子逻辑门，能够使量子态发生演化，量子逻辑门是构成量子线路的基础，量子逻辑门包括单比特量子逻辑门，如Hadamard门(H门，阿达马门)、泡利-X门(X门)、泡利-Y门(Y门)、泡利-Z门(Z门)、RX门、RY门、RZ门等等；多比特量子逻辑门，如CNOT门、CR门、iSWAP门、Toffoli门等等。量子逻辑门一般使用酉矩阵表示，而酉矩阵不仅是矩阵形式，也是一种操作和变换。一般量子逻辑门在量子态上的作用是通过酉矩阵左乘以量子态右矢对应的矩阵进行计算。

参见图2，图2为本发明实施例提供的一种载流子输运仿真方法的流程示意图，可以包括如下步骤：

步骤201：确定半导体器件中载流子输运的物理仿真模型、初始条件和/或边界条件。

其中，载流子为电流的载体；半导体器件中，载流子包括两种类型：电子和空穴。

其中，对于随着时间发展变化的物理过程，某一时刻的状态将影响该时刻以后的过程，该时刻的状态便是初始条件；边界条件是指在求解区域边界上所求解的变量或其导数随着时间和地点的变化规律。边界条件包括给定端点值的第一类边界条件，给定梯度值的第二类边界条件，给定端点值和梯度值的第三类边界条件。本发明中所述的边界条件可以为上述三种边界条件中的任意一种，在此不作限定。

步骤202：确定求解所述物理仿真模型对应的泛定方程。

其中，泛定方程用于表示物理量在空间和时间中的变化规律。

其中，所述物理仿真模型包括以下至少一种：半经典模型、封闭量子模型、开放量子模型；所述半经典仿真模型对应的泛定方程为泊松方程；所述封闭量子模型对应的泛定方程为所述泊松方程和薛定谔方程，所述开放量子模型对应的泛定方程为所述泊松方程和所述薛定谔方程对应的格林函数方程。

具体的，半经典模型对应的泊松方程为：

其中，所述

其中，

其中，所述

具体的，封闭量子模型对应的薛定谔方程为：

其中，所述

其中，

其中，所述f

需要说明的是，封闭量子模型的求解需要联立泊松方程和薛定谔方程，这里的载流子密度不在是半经典模型中的载流子密度的表达式，为上述封闭量子模型中的载流子密度的表达式。

具体的，开放量子模型对应的格林函数方程为：

其中，所述(E+iη)为能量，所述

其中，

其中，所述f

同样，需要说明的是，开放量子模型的求解也需要联立泊松方程和格林函数方程，这里的载流子密度不在是半经典模型或封闭量子模型中的载流子密度的表达式，为上述开放量子模型中的载流子密度的表达式。

其中，上述三种模型，

步骤203：基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度，以实现所述半导体器件中载流子输运的仿真。

具体的，在所述确定半导体器件中载流子输运的初始条件和/或边界条件方面，包括：构建半导体器件的几何模型；基于有限体积法将所述几何模型网格化，得到多个第一控制体积；确定每个所述第一控制体积的初始条件和/或边界条件。

示例性的，如图3所示，图3为本发明实施例提供的一种第一控制体积的结构示意图，图3左边为本发明实施例提供的一种第一控制体积的三维立体图，图3右边为本发明实施例提供的一种第一控制体积的二维平面图。如图3左边所示，第一控制体积的中心P(图中未示出)的西部、东部、背部、南部、底部和顶部的第一控制体积的中心分别对应为W、E、N、S、B和T；如图3右边所示，其中，w、e、n、s、b、t分别为第一控制体积的中心P与W、E、N、S、B、T交界面对应的中心点，P的长Δx为[b,t]，P的宽Δy为[w,e]，P的高Δz为[s,n]，P的体积D

需要说明的是，几何模型被划分为多个与P类似的第一控制体积，每个第一控制体积可以相等，也可以不等；一般来说，只有端点(或是边缘)处的第一控制体积才存在边界条件，几何模型内部其他的第一控制体积不存在边界条件。

可以看出，通过有限体积法可以将几何模型划分成多个控制体积，从而方便后续通过积分形式的守恒方程来构建泛定方程的离散方程，从而实现泛定方程的求解。

具体的，若所述物理仿真模型为所述半经典模型，在所述基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度方面，包括：基于每个所述第一控制体积的初始条件和/或边界条件、所述泊松方程确定数值矩阵形式的解析方程组；对所述数值矩阵形式的解析方程组进行变换，得到残差形式的解析方程组；基于牛顿迭代算法对所述残差形式的方程组进行求解，得到第一静电势；基于所述第一静电势确定所述半导体器件中的载流子密度。

示例性的，网格化后，包含参数(α

对于上述图3中的P，由定积分的基本性质：

取积分平均上式转化为：

运用一阶导的中心差分近似

从而可以得出

其中的[α

定义

从而上式可以简记为

对于泊松方程：

α

对于网格化的多个第一CV(控制体积,Control Volume)，则可以得到多个上述方程，将多个上述方程转换到数值矩阵的形式的解析方程组：

其中，解析函数

上述残差形式的方程组，可以通过牛顿迭代法进行求解，牛顿迭代法如下：

其中，J为雅克比矩阵；

可以看出，对于非线性的泊松方程，首先将离散化的方程组写成数值矩阵形式的解析方程组，然后对数值矩阵形式的解析方程组进行变换，得到残差形式的解析方程组，最后通过牛顿迭代算法对残差形式的方程组进行求解，实现了电势的确定。

与现有技术相比，本发明提供的一种载流子输运仿真方法，将半经典模型对应的泊松方程和半导体器件中载流子输运的初始条件和/或边界条件结合，从而确定该半导体器件中的载流子密度，以实现半导体器件中载流子输运的仿真，进而实现了对半导体器件中载流子输运的研究。

具体的，若所述物理仿真模型为所述封闭量子模型，在所述基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度方面，包括：基于每个所述第一控制体积的初始条件和边界条件、所述泊松方程确定数值矩阵形式的第一解析方程组；确定每个所述第一控制体积对应的所述薛定谔方程，得到数值矩阵形式的第二解析方程组；确定初始载流子密度；基于所述初始载流子密度和所述第一解析方程组对所述第二解析方程组迭代求解，确定所述半导体器件中的载流子密度。

其中，第一解析方程组如下：

其中，第二解析方程组如下：

其中，本征函数

需要说明的是，边界条件

需要说明的是，第一解析方程组的推导与半经典模型中泊松方程的离散化过程相同，对于P处的薛定谔方程，只需要将φ换成特征函数ψ，其余过程也与泊松方程的离散化过程相类似，其中的参数如下：

可以看出，[M]仅包含泊松方程中的拉普拉斯算符项，[M

进一步的，在所述基于所述初始载流子密度、所述第一解析方程组对所述第二解析方程组迭代求解，得到所述半导体器件中的载流子密度方面，包括：基于所述初始载流子密度求解所述第一解析方程组，得到第二静电势；基于所述第二静电势确定估计载流子密度；基于所述估计载流子密度求解所述第二解析方程组，得到第一特征函数；基于所述第一特征函数确定第一目标载流子密度；若所述第一目标载流子密度与所述初始载流子密度的差值大于第一预设误差，则将所述第一目标载流子密度作为新的所述初始载流子密度，然后执行所述基于所述初始载流子密度求解所述第一解析方程组，得到第二静电势；若所述第一目标载流子密度与所述初始载流子密度的差值小于或等于所述第一预设误差，则将所述第一目标载流子密度作为所述半导体器件中的载流子密度。

需要说明的是，初始载流子密度为预先给定的；估计载流子密度为估计值，可以更快的迭代求解出第一目标载流子密度；第一目标载流子密度与初始载流子密度的第一预设误差也是预先给定的。

与现有技术相比，本发明提供的一种载流子输运仿真方法，将封闭量子模型对应的泊松方程和薛定谔方程与半导体器件中载流子输运的初始条件和/或边界条件结合，从而确定该半导体器件中的载流子密度，以实现半导体器件中载流子输运的仿真，进而实现了对半导体器件中载流子输运的研究。

具体的，若所述物理仿真模型为所述开放量子模型，在所述基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度方面，包括：基于每个所述第一控制体积的初始条件和边界条件、所述泊松方程确定数值矩阵形式的第三解析方程组；确定每个所述第一控制体积对应的所述格林函数方程，得到数值矩阵形式的第四解析方程组；确定初始静电势；基于所述初始静电势和所述第三解析方程组对所述第四解析方程组迭代求解，确定所述半导体器件中的载流子密度。

其中，第三解析方程组如下：

其中，第四解析方程组如下：

其中，

Γ

得到的。Broadening函数的规模和

其中，g

需要说明的是，第三解析方程组的推导与半经典模型中泊松方程的离散化过程相同，对于P处的格林函数方程，只需要将φ换成格林函数G，其余过程也与泊松方程的离散化过程相类似，其中的参数如下：

具体实现过程参见上述的泊松方程离散化的推导过程，在此就不一一详述了。

可以看出，[M]仅包含泊松方程中的拉普拉斯算符项，

具体的，在所述基于所述初始静电势、所述第三解析方程组对所述第四解析方程组迭代求解，得到所述半导体器件中的载流子密度方面，包括：基于所述初始静电势求解所述第四解析方程组，得到第二特征函数；基于所述第二特征函数确定第一初始载流子密度；基于所述第一初始载流子密度求解所述第三解析方程组，得到目标静电势；基于所述目标静电势求解所述第四解析方程组，得到第三特征函数；基于所述第三特征函数确定第二目标载流子密度；若所述第二目标载流子密度与所述第一初始载流子密度的差值大于第二预设误差，则将所述目标静电势作为新的所述初始静电势，然后执行所述基于所述初始静电势求解所述第四解析方程组，得到第二特征函数；若所述第二目标载流子密度与所述第一初始载流子密度的差值小于或等于所述第二预设误差，则将所述第二目标载流子密度作为所述半导体器件中的载流子密度。

需要说明的是，初始静电势为预先给定的；第二目标载流子密度与第一初始载流子密度的第二预设误差也是预先给定的。

与现有技术相比，本发明提供的一种载流子输运仿真方法，将开放量子模型对应的泊松方程和格林函数方程与半导体器件中载流子输运的初始条件和/或边界条件结合，从而确定该半导体器件中的载流子密度，以实现半导体器件中载流子输运的仿真，进而实现了对半导体器件中载流子输运的研究。

进一步的，在所述确定所述半导体器件中的载流子密度之后，所述方法还包括：基于所述半导体器件中的载流子密度确定所述几何模型中载流子密度大于或等于预设阈值的目标区域；基于有限体积法将所述目标区域网格化，得到多个第二控制体积，所述第二控制体积小于所述第一控制体积；基于多个所述第二控制体积确定新的所述半导体器件中的载流子密度。

需要说明的是，预设阈值是预先给定的；确定第二静电势的方法可以是与第一控制体积对应的求解方式相同，也可以是通过其他方式去确定。例如，第一对应体积对应的求解方法为泊松方程，第二控制体积的求解方法为泊松方程或泊松方程与薛定谔方程的联立或泊松方程与格林函数方程的联立，在此不做限定。

与现有技术相比，本发明提供的一种载流子输运仿真方法，将开放量子模型对应的泊松方程和格林函数方程与半导体器件中载流子输运的初始条件和/或边界条件结合，从而确定该半导体器件中的载流子密度，以实现半导体器件中载流子输运的仿真，进而实现了对半导体器件中载流子输运的研究。

参见图4，图4为本发明实施例提供的一种载流子输运仿真装置的结构示意图，与图2所述的流程相对应，所述装置包括：

第一确定单元401，用于确定半导体器件中载流子输运的物理仿真模型、初始条件和/或边界条件；

第二确定单元402，用于确定求解所述物理仿真模型对应的泛定方程；

第三确定单元403，用于基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度，以实现所述半导体器件中载流子输运的仿真。

具体的，所述物理仿真模型包括以下至少一种：半经典模型、封闭量子模型、开放量子模型；所述半经典仿真模型对应的泛定方程为泊松方程；所述封闭量子模型对应的泛定方程为所述泊松方程和薛定谔方程，所述开放量子模型对应的泛定方程为所述泊松方程和所述薛定谔方程对应的格林函数方程。

具体的，在所述确定半导体器件中载流子输运的初始条件和/或边界条件发方面，第一确定单元401，具体用于：

构建半导体器件的几何模型；

基于有限体积法将所述几何模型网格化，得到多个第一控制体积；

确定每个所述第一控制体积的初始条件和/或边界条件。

具体的，若所述物理仿真模型为所述半经典模型，在所述基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度方面，第三确定单元403，具体用于：

基于每个所述第一控制体积的初始条件和/或边界条件、所述泊松方程确定数值矩阵形式的解析方程组；

对所述数值矩阵形式的解析方程组进行变换，得到残差形式的解析方程组；

基于牛顿迭代算法对所述残差形式的方程组进行求解，得到第一静电势；

基于所述第一静电势确定所述半导体器件中的载流子密度。

具体的，若所述物理仿真模型为所述封闭量子模型，在所述基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度方面，第三确定单元403，具体用于：

基于每个所述第一控制体积的初始条件和边界条件、所述泊松方程确定数值矩阵形式的第一解析方程组；

确定每个所述第一控制体积对应的所述薛定谔方程，得到数值矩阵形式的第二解析方程组；

确定初始载流子密度；

基于所述初始载流子密度和所述第一解析方程组对所述第二解析方程组迭代求解，确定所述半导体器件中的载流子密度。

具体的，若所述物理仿真模型为所述开放量子模型，在所述基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度方面，第三确定单元403，具体用于：

基于每个所述第一控制体积的初始条件和边界条件、所述泊松方程确定数值矩阵形式的第三解析方程组；

确定每个所述第一控制体积对应的所述格林函数方程，得到数值矩阵形式的第四解析方程组；

确定初始静电势；

基于所述初始静电势和所述第三解析方程组对所述第四解析方程组迭代求解，确定所述半导体器件中的载流子密度。

具体的，在所述确定所述半导体器件中的载流子密度之后，第三确定单元403，还用于：

基于所述半导体器件中的载流子密度确定所述几何模型中载流子密度大于或等于预设阈值的目标区域；

基于有限体积法将所述目标区域网格化，得到多个第二控制体积，所述第二控制体积小于所述第一控制体积；

基于多个所述第二控制体积确定新的所述半导体器件中的载流子密度。

与现有技术相比，本发明提供的一种载流子输运仿真方法，将不同的物理仿真模型对应的泛定方程与半导体器件中载流子输运的初始条件和/或边界条件结合，从而确定该半导体器件中的载流子密度，以实现半导体器件中载流子输运的仿真，进而实现了对半导体器件中载流子输运的研究。

本发明的再一实施例提供了一种存储介质，所述存储介质中存储有计算机程序，其中，所述计算机程序被设置为运行时执行上述任一项中方法实施例中的步骤。

具体的，在本实施例中，上述存储介质可以被设置为存储用于执行以下步骤的计算机程序：

确定半导体器件中载流子输运的物理仿真模型、初始条件和/或边界条件；

确定求解所述物理仿真模型对应的泛定方程；

基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度，以实现所述半导体器件中载流子输运的仿真。

本实施例中，上述存储介质可以包括但不限于：U盘、只读存储器(Read-OnlyMemory，简称为ROM)、随机存取存储器(Random Access Memory，简称为RAM)、移动硬盘、磁碟或者光盘等各种可以存储计算机程序的介质。

本发明的再一实施例还提供了一种电子装置，包括存储器和处理器，所述存储器中存储有计算机程序，所述处理器被设置为运行所述计算机程序以执行上述任一项中方法实施例中的步骤。

具体的，上述载流子装置还可以包括传输设备以及输入输出设备，其中，该传输设备和上述处理器连接，该输入输出设备和上述处理器连接。

具体的，在本实施例中，上述处理器可以被设置为通过计算机程序执行以下步骤：

确定半导体器件中载流子输运的物理仿真模型、初始条件和/或边界条件；

确定求解所述物理仿真模型对应的泛定方程；

基于所述初始条件和/或所述边界条件、所述泛定方程确定所述半导体器件中的载流子密度，以实现所述半导体器件中载流子输运的仿真。

以上依据图式所示的实施例详细说明了本发明的构造、特征及作用效果，以上所述仅为本发明的较佳实施例，但本发明不以图面所示限定实施范围，凡是依照本发明的构想所作的改变，或修改为等同变化的等效实施例，仍未超出说明书与图示所涵盖的精神时，均应在本发明的保护范围内。