CN112784472A - 循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法 - Google Patents

Abstract

本发明公开了一种循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，包括：建立循环神经网络，循环神经网络为长短时记忆网络；长短时记忆网络包括T个按时间顺序排列的LSTM细胞，每个LSTM细胞具有输入值x_t和输出值h_t，LSTM细胞内具有参数(W,b)；将根据量子条件主方程得到的电流的散粒噪声谱S(ω)，替代输入值x_t；利用量子条件主方程中的密度矩阵迹，替代输出值h_t；利用前后时刻的量子条件主方程中的密度矩阵迹之间联系，替代参数(W,b)。本发明建立了循环神经网络中的长短时记忆网络和量子条件主方程的联系，利用量子系统产生的散粒噪声谱的数据，解决求解量子条件主方程时方程无限循环闭合的难题，实现循环神经网络对量子条件主方程的模拟。

Description

循环神经网络模拟量子输运过程中的量子条件主方程的模拟 方法

技术领域

本发明涉及一种循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法。

背景技术

量子输运现象作为介观系统中的一种重要的物理现象，近年来被广泛研究。对于传统器件，信噪比可以通过抑制散粒噪声来提高，但是在量子点构成的系统中，散粒噪声不会被无限的减少。事实上，量子器件的量子输运噪声并不是一定有害的。这些和精细时间有关的散粒噪声能敏感的反应输运过程中精细的动力学，丰富的量子输运性质以及在其中的精细的能量尺度。因此，在研究低维介观纳米器件的输运特点的过程中，对量子散粒噪声系统的测试和分析会成为一个重要的理论工具和方法。

对于理论计算，我们需要去面对一个含有噪声和量子点的开放系统并研究其性质。前人提出了很多方法去研究它，包括Buttiker和Beenaker等人提出的散射矩阵方法、非平衡格林函数方法、量子主方程方法。不同于前人的研究，李兴奇等人沿着Gurvitz的方法提出了一个研究电荷比特细致输运过程的条件主方程。

尽管量子条件主方程能够描述细致的描述电荷的输运过程，但是进一步研究此过程相关的物理量却是非常困难的，因为其是一个无限递归的微分方程系统。因此，对输运过程中的量子条件主方程的求解就显得极为重要。

因此，提供循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，属于本领域亟待解决的问题。其中，模拟量子条件主方程的循环神经网络可以用于指导微纳量子器件的设计。

发明内容

本发明的目的在于克服现有技术的不足，提供循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法。

本发明的目的是通过以下技术方案来实现的：

本发明的第一方面，提供循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，包括以下步骤：

建立一个循环神经网络，所述循环神经网络为长短时记忆网络；所述长短时记忆网络包括T个按时间顺序排列的LSTM细胞，每个LSTM细胞具有输入值x_t和输出值h_t，输出值h_t会传入下一时刻的LSTM细胞中，LSTM细胞内具有参数(W,b)；

将根据量子条件主方程得到的电流的散粒噪声谱，替代输入值x_t；利用量子条件主方程中的密度矩阵迹，替代输出值h_t；利用前后时刻即t-1时刻和t时刻的量子条件主方程中的密度矩阵迹之间联系，替代参数(W,b)；

利用量子输运过程中产生的散粒噪声谱的数据来训练所述循环神经网络从而达到模拟量子条件主方程的目的；所述量子输运过程对应一个可实现的物理实际系统。

进一步地，所述量子条件主方程由二能级量子电荷比特输运系统推导；所述二能级量子电荷比特输运系统包括量子点系统S和电源V，量子点系统S的左电极L与电源V的正极连接，量子点系统S的右电极R与电源V的负极连接；所述二能级量子电荷比特输运系统的总哈密顿量为：

式中，H_s表示的是量子点系统S的哈密顿量，H_E表示的是左电极L和右电极R的哈密顿量，H′表示的是量子点系统S和电极之间相互作用的哈密顿量。

进一步地，假设量子点系统S和环境之间的相互作用不是很强，把H′当作微扰来处理，根据二阶矩累积展开和Lindblad方程，得到描述量子输运过程的量子主方程：

式中，刘维尔超算符定义为：

G(t,τ)是与量子点系统S哈密顿量H_s有关的传播子；量子点系统S的约化密度矩阵为ρ(t)＝Tr_E[ρ_T(t)]，<(…)>＝Tr_E[(…)ρ_E]，ρ_E表示电极的密度矩阵；i表示虚数单位，ρ(t)表示在t时刻的密度矩阵，τ表示小于时间t的任意时刻，

表示ρ(t)对时间的一阶导数；

对电极所处的希尔伯特空间进行划分，用E⁽ⁿ⁾表示在Δt时间内有n个电子经过量子点系统S时电极处于的空间，则电极所在的空间就可以表示为

把所述假设带到量子主方程中，得到量子条件主方程：

这里

表示在Δt时间内由n个电子穿过量子点系统S时系统的约化密度矩阵，ρ⁽ⁿ⁾(t)表示在Δt时间内由n个电子穿过量子点系统S时系统的约化密度矩阵，

表示ρ⁽ⁿ⁾(t)对时间的一阶导数。

进一步地，H_s、H_E和H′的具体形式为：

式中，s表示电子的自旋，↑,↓分别表示自旋向上和自旋向下；j表示能级，∈_j表示第j个能级的能量；

分别表示电子处于第j个能级上且自旋为s的产生/湮灭算符；ω是两个电子占据同一能级但自旋不相同的库伦作用能，n_j↑、n_j↓、n_js分别表示电子占据第j个能级且自旋为↑,↓,s时的粒子数算符；C_E是与占据能级的电子数有关的电荷能；α表示电极；k表示的是电子的动量；∈_αks表示电极上动量为k的电子的能量，考虑到电极上的电子处于热统计平衡状态，其分布函数为：

μ表示费米能量，考虑外部电压是对成的加在系统上的，这里费米能量等于μ_L＝eV/2,μ_R＝-eV/2；T表示的是温度，就是量子输运系统处于的温度，k_B表示玻尔兹曼常数；

表示α电极上自旋为s、动量为k的电子的产生算符；c_αks表示α电极上自旋为s、动量为k的电子的湮灭算符；τ_αkj表示系统和环境之间的相互作用强度，h.c.表示厄米共轭。

进一步地，假设

并在马尔可夫近似下，所述量子条件主方程改写为：

式中，

是电极中电子的谱函数；ρ⁽ⁿ⁾即ρ⁽ⁿ⁾(t)，

表示处于第j个能级上、自旋为s的电子的产生算符，

表示处于第j个能级上、自旋为s的电子湮灭算符，γ表示谱函数

中的自变量；

流过量子点系统S中的电流表示为：

式中，P(n,t)表示在Δt时间内由n个电子穿过量子点系统的概率，e表示单位电荷，n表示单位时间穿过量子点系统S的电子数目；

根据MacDonald公式，电流的散粒噪声谱表示为：

式中，ω表示散粒噪声S(ω)函数中的自变量。

进一步地，所述二能级量子电荷比特输运系统的密度矩阵在不同时刻之间的联系，利用Kraus算符即量子隐马尔可夫进行表示：即

式中m表示不同的K,Km是第m个Kraus算符；该公式与所述量子主方程等价；

而在量子条件主方程下，将ρ(t)＝∑_nρ⁽ⁿ⁾(t)带入上式就有

根据量子条件主方程的改写方程，ρ⁽ⁿ⁾(t+Δt)和ρ⁽ⁿ⁾(t),ρ^(n-1)(t),ρ⁽ⁿ⁺¹⁾(t)相关，因此联合上式可知：

即该式与量子条件主方程的改写方程有关系，目标就是构建不含时的映射

进一步地，所述LSTM细胞的输入值x_t和输出值h_t的关系，由下式方程给出：

f_t＝σ(W_f·[h_t-1,x_t]+b_f)

i_t＝σ(W_i·[h_t-1,x_t]+b_i)

O_t＝σ(W_O·[h_t-1,x_t]+b_O)

h_t＝O_t×tanh(C_t)

式中，(W_f,W_i,W_c,W_o,b_f,b_i,b_c,b_o)作为所述参数(W,b)；f_t为遗忘门限层的输出，i_t和

为输入门限层的输出，O_t和h_t为输出门限层的输出。

进一步，利用散粒噪声谱的数据来构建密度矩阵迹之间即t-1时刻和t时刻之间的联系，即构建

把Tr[ρ⁽ⁿ⁾(t)]和长短时记忆网络中的h_t参数相对应；

将长短时记忆网络中的参数(W,b)就可以充当

的作用。

进一步地，所述方法还包括以下步骤：

用Tr[ρ⁽ⁿ⁾(t)]对总电流的贡献度的方法来有效的截断n值，n表示量子条件主方程的改写方程中的粒子数；定义评估函数：

式中，M是数值实验上能取到的最大粒子数值，P_M对应为有M个电子流过量子点系统的概率值；

通过不断调整M的值，绘制出了E(M)随M的变化图像，根据所述变化图像来确定M的取值。

进一步地，所述量子输运过程中产生的散粒噪声谱的数据分为训练数据和测试数据；

利用训练数据来训练所述循环神经网络，得到训练数据的误差随迭代次数的第一关系；利用所述测试数据来测试所述循环神经网络，得到测试数据的误差随迭代次数的第二关系；

利用第一关系和第二关系确定长短时记忆网络模拟量子条件主方程的模拟效果。

本发明的有益效果是：

(1)在本发明的一示例性实施例中，建立了循环神经网络中的长短时记忆网络和量子条件主方程的联系，得到二者的等价关系。同时利用量子系统产生的散粒噪声谱的数据，解决求解量子条件主方程时方程无限循环闭合的难题，实现循环神经网络对量子条件主方程的模拟。

(2)在本发明的又一示例性实施例中，公开了量子条件主方程的推导前提，即二能级量子电荷比特输运系统的具体实现方式；同时在本发明的又一示例性实施例中，公开了长短时记忆网络的具体结构。

(3)在本发明的又一示例性实施例中，用Tr[ρ⁽ⁿ⁾(t)]对总电流的贡献度的方法来有效的截断n值，进一步解决求解量子条件主方程时方程无限循环闭合的难题，实现循环神经网络对量子条件主方程的模拟。

附图说明

图1为本发明一示例性实施例公开的方法流程图；

图2为本发明一示例性实施例公开的技术实现路线图；

图3为本发明一示例性实施例公开的二能级量子电荷比特输运系统结构示意图；

图4为本发明一示例性实施例公开的量子隐马尔可夫计算图；

图5为本发明一示例性实施例公开的长短时记忆网络计算图；

图6为本发明一示例性实施例公开的长短时记忆网络的LSTM细胞结构示意图；

图7为本发明一示例性实施例公开的截断判断中E(M)随M的变化图；

图8为本发明一示例性实施例公开的训练数据的误差随迭代次数的关系图；

图9为本发明一示例性实施例公开的测试数据的误差随迭代次数的关系图。

具体实施方式

下面结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

在本申请使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本申请。在本申请和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。

应当理解，尽管在本申请可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本申请范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。

此外，下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。

具体地，在下述示例性实施例中，通过推导出描述量子输运过程的量子条件主方程，然后发现量子隐马尔可夫过程和量子主方程有一定的联系，经过展开量子隐马尔可夫过程的计算图发现了其和循环神经网络之间的联系。然后用量子输运过程(一个可实现的物理实际系统)中产生的噪声谱的数据来训练循环神经网络从而达到模拟量子条件主方程的目的。可以用于微纳量子器件的设计。

参见图1，图1示出了本发明的一示例性实施例提供的循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，包括以下步骤：

建立一个循环神经网络，所述循环神经网络为长短时记忆网络；所述长短时记忆网络包括T个按时间顺序排列的LSTM细胞，每个LSTM细胞具有输入值x_t和输出值h_t，输出值h_t会传入下一时刻的LSTM细胞中，LSTM细胞内具有参数(W,b)；

将根据量子条件主方程得到的电流的散粒噪声谱S(ω)(其中，散粒噪声谱S(ω)表示实际计算值，在有实验条件的情况下，可以从量子输运系统中采集)，替代输入值x_t；利用量子条件主方程中的密度矩阵迹Tr[ρ⁽ⁿ⁾(t)]，替代输出值h_t；利用前后时刻即t-1时刻和t时刻的量子条件主方程中的密度矩阵迹之间联系

替代参数(W,b)；

利用量子输运过程中产生的散粒噪声谱的数据来训练所述循环神经网络从而达到模拟量子条件主方程的目的；所述量子输运过程对应一个可实现的物理实际系统。

具体地，在该示例性实施例中，建立了循环神经网络中的长短时记忆网络和量子条件主方程的联系，得到二者的等价关系。同时利用量子系统产生的散粒噪声谱的数据，解决求解量子条件主方程时方程无限循环闭合的难题，实现循环神经网络对量子条件主方程的模拟。

需要说明的是，该三个参数与长短时记忆网络的三个参数进行对应，是因为：量子条件主方程和循环神经网络是在展开计算图上是等价的，即图4和图5(下述示例性实施例将详细进行展开)；同时，由于长短时记忆网络是与量子输运系统相关的，而散粒噪声谱则是用于描述该量子输运系统/输运过程，因此利用散粒噪声谱作为长短时记忆网络的输入参数x_t。

也就是说，x_t参数就是网络的输入参数，是一个序列数据，散粒噪声谱也是一个序列数据，在实际的操作中，就把散粒噪声谱这个序列数据当作是x_t输入到网络中，也就是说第一步输入散粒噪声谱的第一个数据，第二步输入散粒噪声谱的第二个数据，以此类推。φ和网络参数的对应要从宏观的角度来看，h的前一步和这一步两个步骤之间的变化是通过φ联系起来的，然后h参数的上一步和这一步之间的联系是通过网络参数联系起来的，不管是网络参数还是φ，他们都做了同一个事情，就是联系某一个量的上一步和这一步，然后又因为h和ρ对应，我们就把φ和网络参数对应。

在其他具体示例性实施例中，该训练后的循环神经网络或长短时记忆网络可以用于指导设计微纳量子器件的技术领域。

需要说明的是，在训练过程中，输入值x_t即量子输运过程中产生的散粒噪声谱S(ω)的数据为已知参数，参数(W,b)为待训练的参数，输出值h_t即量子条件主方程中的密度矩阵迹Tr[ρ⁽ⁿ⁾(t)]是未知数据。需要说明的是，输出值h_t是未知的，是通过长短时记忆网络计算出来的，在计算过程中，只需给定一个初始值h₀，后面的时间步都是通过长短时记忆网络计算得到的。

具体地，下述示例性实施例首先通过推导出描述量子输运过程的量子条件主方程，然后发现量子隐马尔可夫过程和量子主方程有一定的联系，经过展开量子隐马尔可夫过程的计算图发现了其和循环神经网络之间的联系，如图2所示。

更优地，在一示例性实施例中，所述量子条件主方程由二能级量子电荷比特输运系统推导；如图3所示，所述二能级量子电荷比特输运系统包括量子点系统S和电源V，量子点系统S的左电极L与电源V的正极连接，量子点系统S的右电极R与电源V的负极连接；电子在外部电压的激励下从量子点中流过。

所述二能级量子电荷比特输运系统的总哈密顿量为：

式中，H_s表示的是量子点系统S的哈密顿量，H_E表示的是左电极L和右电极R的哈密顿量，H′表示的是量子点系统S和电极之间相互作用的哈密顿量。

更优地，而在又一示例性实施例中，H_s、H_E和H′的具体形式为：

式中，s表示电子的自旋，↑,↓分别表示自旋向上和自旋向下；j表示能级，∈_j表示第j个能级的能量；

分别表示电子处于第j个能级上且自旋为s的产生/湮灭算符；ω是两个电子占据同一能级但自旋不相同的库伦作用能，n_j↑、n_j↓、n_js分别表示电子占据第j个能级且自旋为↑,↓,s时的粒子数算符；C_E是与占据能级的电子数有关的电荷能；α表示电极；k表示的是电子的动量；∈_αks表示电极上动量为k的电子的能量，考虑到电极上的电子处于热统计平衡状态，其分布函数为：

μ表示费米能量，考虑外部电压是对成的加在系统上的，这里费米能量等于μ_L＝eV/2,μ_R＝-eV/2；T表示的是温度，就是量子输运系统处于的温度，k_B表示玻尔兹曼常数；

表示α电极上自旋为s、动量为k的电子的产生算符；c_αks表示α电极上自旋为s、动量为k的电子的湮灭算符；τ_αkj表示系统和环境之间的相互作用强度，h.c.表示厄米共轭。

更优地，在一示例性实施例中，假设量子点系统S和环境之间的相互作用不是很强，可以把H′当作微扰来处理，根据二阶矩累积展开和Lindblad方程，得到描述量子输运过程的量子主方程：

式中，刘维尔超算符定义为：

G(t,τ)是与量子点系统S哈密顿量H_s有关的传播子(格林函数)；量子点系统S的约化密度矩阵为ρ(t)＝Tr_E[ρ_T(t)]，<(…)＞＝Tr_E[(…)ρ_E]，ρ_E表示电极的密度矩阵；i表示虚数单位，ρ(t)表示在t时刻的密度矩阵，τ表示小于时间t的任意时刻，

表示ρ(t)对时间的一阶导数；

如果我们对电极所处的希尔伯特空间进行划分，用E⁽ⁿ⁾表示在Δt时间内有n个电子经过量子点系统S时电极处于的空间，则电极所在的空间就可以表示为

值得提的是，当n＝0时，由于没有电子经过量子点系统，此时所处的希尔伯特空间由左右两个独立的电极的子空间所张成，即

把所述假设带到量子主方程(4)中，得到量子条件主方程：

这里

表示在Δt时间内由n个电子穿过量子点系统S时系统的约化密度矩阵，，ρ⁽ⁿ⁾(t)表示在Δt时间内由n个电子穿过量子点系统S时系统的约化密度矩阵，

表示ρ⁽ⁿ⁾(t)对时间的一阶导数。

更优地，在一示例性实施例中，假设

并在马尔可夫近似下，所述量子条件主方程改写为：

式中，

是电极中电子的谱函数；ρ⁽ⁿ⁾即ρ⁽ⁿ⁾(t)，

表示处于第j个能级上、自旋为s的电子的产生算符，

表示处于第j个能级上、自旋为s的电子湮灭算符，γ表示谱函数

中的自变量。下述内容将式(6)称为量子条件主方程的改写方程。

流过量子点系统S中的电流表示为：

式中，P(n,t)表示在Δt时间内由n个电子穿过量子点系统S的概率，e表示单位电荷，n表示单位时间穿过量子点系统S的电子数目；

根据MacDonald公式，电流的散粒噪声谱表示为：

式中，ω表示散粒噪声S(ω)函数中的自变量，可以类似理解为傅里叶变换中的频率。

上述示例性实施例推导出描述量子输运过程的量子条件主方程，下述内容发现量子隐马尔可夫过程和量子主方程有一定的联系，经过展开量子隐马尔可夫过程的计算图发现了其和循环神经网络之间的联系。具体地：

更优地，在一示例性实施例中，所述二能级量子电荷比特输运系统的密度矩阵在不同时刻之间的联系，利用Kraus算符即量子隐马尔可夫进行表示：即

式中m表示不同的K,Km是第m个Kraus算符；该公式与所述量子主方程即公式(4)等价；

而在量子条件主方程下，将ρ(t)＝∑_nρ⁽ⁿ⁾(t)带入上式就有

根据量子条件主方程的改写方程即公式(6)，可以看出ρ⁽ⁿ⁾(t+Δt)和ρ⁽ⁿ⁾(t),ρ^(n-1)(t),ρ⁽ⁿ⁺¹⁾(t)相关，因此联合上式可知：

即该式与量子条件主方程的改写方程(6)有关系，而我们的目标就是构建不含时的映射

下述内容对比了量子隐马尔科夫的计算图和循环神经网络的计算图，发现两者有很高的相似性。具体地：

图4为量子隐马尔可夫计算图。对于量子隐马尔科夫，密度矩阵的演化过程实际上是一个不随时间变化的参数的循环计算的过程，即每一次变换都是利用φ参数。而图5为长短时记忆网络计算图，在长短时记忆网络训练后，每一次计算的(W,b)参数不变。因此可以看到，过展开量子隐马尔可夫过程的计算图发现了其和循环神经网络(长短时记忆网络计算图)之间的联系。

具体地，长短时记忆网络是循环神经网络的一个子类，在处理时间序列数据上有很大的优势。图6示出了长短时记忆网络LSTM细胞的具体结构图。

更优地，在一示例性实施例中，所述LSTM细胞的输入值x_t和输出值h_t的关系，由下式方程给出：

式中，(W_f,W_i,W_c,W_o,b_f,b_i,b_c,b_o)作为所述参数(W,b)；f_t为遗忘门限层的输出，i_t和

为输入门限层的输出，O_i和h_i为输出门限层的输出。该部分的内容属于现有技术的具体内容，在此不进行赘述。

x_t是当前的输入值，h_t可以作为当前的输出值进行输出并且传入下一时刻的LSTM细胞中。

更优地，在一示例性实施例中，由于电流的散粒噪声谱(即式(8))的计算只用到了密度矩阵迹

的信息，所以只能用散粒噪声谱的数据来构建密度矩阵迹之间即t-1时刻和t时刻之间的联系，即构建

(构建出来的关系不含时间)；

密度矩阵在演化的过程中始终保持迹为1，因此根据ρ(t)＝∑_nρ⁽ⁿ⁾(t)可知，在量子条件主方程下(即式(5))只能保证总体密度矩阵的迹Tr[ρ(t)]＝Tr[∑_nρ⁽ⁿ⁾(t)]＝∑_nTr[ρ⁽ⁿ⁾(t)]＝1保持不变，并不能保证Tr[ρ⁽ⁿ⁾(t)]不发生变化；这个公式(Tr[ρ⁽ⁿ⁾(t)])不会恒等于1，所以才有模拟的必要；

密度矩阵的演化过程实际上是一个不随时间变化的参数的循环计算的过程，这与循环神经网络的计算思想是一致的，即循环神经网络训练后，参数(W,b)保持不变；

把Tr[ρ⁽ⁿ⁾(t)]和长短时记忆网络中的h_t参数相对应，因为h_t参数在经过SoftMax函数作用后也会有Tr[h_t]＝Tr[∑_nh_t ⁽ⁿ⁾]＝∑_nTr[h_t ⁽ⁿ⁾]＝1，这与Tr[ρ⁽ⁿ⁾(t)]的行为相同(即整体求和为1，但是单体变化)；

长短时记忆网络中的参数(W,b)就可以充当

的作用，联系上一步和下一步的关系。

我们的目标就是用二能级量子系统产生的噪声谱的数据构建这样一个关系。

但是这里还存在一个问题，在方程(6)中粒子数n的取值可以取到任意值，但是在数值实验中我们必须对n做一个有效的截断。

更优地，在一示例性实施例中，所述方法还包括以下步骤：

用Tr[ρ⁽ⁿ⁾(t)]对总电流的贡献度的方法来有效的截断n值，n表示量子条件主方程的改写方程中的粒子数；定义评估函数：

式中，M是数值实验上能取到的最大粒子数值，P_M对应为有M个电子流过量子点系统的概率值；

通过不断调整M的值，绘制出了E(M)随M的变化图像，如图7所示。可以根据图7来确定M的取值，在其中一示例性实施例中，取M＝50，Tr[ρ⁽⁵⁰⁾(t)]对总电流的贡献只占0.4％。

更优地，所述量子输运过程中产生的散粒噪声谱的数据分为训练数据和测试数据；

利用训练数据来训练所述循环神经网络，得到训练数据的误差随迭代次数的第一关系；利用所述测试数据来测试所述循环神经网络，得到测试数据的误差随迭代次数的第二关系；

利用第一关系和第二关系确定长短时记忆网络模拟量子条件主方程的模拟效果。

具体地，图8展示了训练数据的误差随迭代次数的关系，图9展示了测试数据的误差随迭代次数的关系。从图中中可以看到，误差随着迭代次数逐渐减小直至收敛，这说明我们通过长短时记忆网络很好的构建了

也就是说，我们利用长短时记忆网络很好的模拟了量子条件主方程。

显然，上述实施例仅仅是为清楚地说明所作的举例，而并非对实施方式的限定，对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其他不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明创造的保护范围之中。

Claims (10)

1.循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，其特征在于：包括以下步骤：

建立一个循环神经网络，所述循环神经网络为长短时记忆网络；所述长短时记忆网络包括T个按时间顺序排列的LSTM细胞，每个LSTM细胞具有输入值x_t和输出值h_t，输出值h_t会传入下一时刻的LSTM细胞中，LSTM细胞内具有参数(W,b)；

将根据量子条件主方程得到的电流的散粒噪声谱，替代输入值x_t；利用量子条件主方程中的密度矩阵迹，替代输出值h_t；利用前后时刻即t-1时刻和t时刻的量子条件主方程中的密度矩阵迹之间联系，替代参数(W,b)；

利用量子输运过程中产生的散粒噪声谱的数据来训练所述循环神经网络从而达到模拟量子条件主方程的目的；所述量子输运过程对应一个可实现的物理实际系统。

2.根据权利要求1所述的循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，其特征在于：所述量子条件主方程由二能级量子电荷比特输运系统推导；所述二能级量子电荷比特输运系统包括量子点系统S和电源V，量子点系统S的左电极L与电源V的正极连接，量子点系统S的右电极R与电源V的负极连接；所述二能级量子电荷比特输运系统的总哈密顿量为：

式中，H_s表示的是量子点系统S的哈密顿量，H_E表示的是左电极L和右电极R的哈密顿量，H′表示的是量子点系统S和电极之间相互作用的哈密顿量。

3.根据权利要求2所述的循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，其特征在于：假设量子点系统S和环境之间的相互作用不是很强，把H′当作微扰来处理，根据二阶矩累积展开和Lindblad方程，得到描述量子输运过程的量子主方程：

式中，刘维尔超算符定义为：

G(t,τ)是与量子点系统S哈密顿量H_s有关的传播子；量子点系统S的约化密度矩阵为ρ(t)＝Tr_E[ρ_T(t)]，<(…)>＝Tr_E[(…)ρ_E]，ρ_E表示电极的密度矩阵；i表示虚数单位，ρ(t)表示在t时刻的密度矩阵，τ表示小于时间t的任意时刻，ρ(t)表示ρ(t)对时间的一阶导数；

对电极所处的希尔伯特空间进行划分，用E⁽ⁿ⁾表示在Δt时间内有n个电子经过量子点系统S时电极处于的空间，则电极所在的空间就可以表示为

把所述假设带到量子主方程中，得到量子条件主方程：

这里

表示在Δt时间内由n个电子穿过量子点系统S时系统的约化密度矩阵，ρ⁽ⁿ⁾(t)表示在Δt时间内由n个电子穿过量子点系统S时系统的约化密度矩阵，ρ⁽ⁿ⁾(t)表示ρ⁽ⁿ⁾(t)对时间的一阶导数。

4.根据权利要求3所述的循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，其特征在于：H_s、H_E和H′的具体形式为：

式中，s表示电子的自旋，↑,↓分别表示自旋向上和自旋向下；j表示能级，∈_j表示第j个能级的能量；

分别表示电子处于第j个能级上且自旋为s的产生/湮灭算符；ω是两个电子占据同一能级但自旋不相同的库伦作用能，n_j↑、n_j↓、n_js分别表示电子占据第j个能级且自旋为↑,↓,s时的粒子数算符；C_E是与占据能级的电子数有关的电荷能；α表示电极；k表示的是电子的动量；∈_αks表示电极上动量为k的电子的能量，考虑到电极上的电子处于热统计平衡状态，其分布函数为：

μ表示费米能量，考虑外部电压是对成的加在系统上的，这里费米能量等于μ_L＝eV/2,μ_R＝-eV/2；T表示的是温度，就是量子输运系统处于的温度，k_B表示玻尔兹曼常数；

表示α电极上自旋为s、动量为k的电子的产生算符；c_αks表示α电极上自旋为s、动量为k的电子的湮灭算符；τ_αkj表示系统和环境之间的相互作用强度，h.c.表示厄米共轭。

5.根据权利要求4所述的循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，其特征在于：假设

并在马尔可夫近似下，所述量子条件主方程改写为：

式中，

是电极中电子的谱函数；ρ⁽ⁿ⁾即ρ⁽ⁿ⁾(t)，

表示处于第j个能级上、自旋为s的电子的产生算符，

表示处于第j个能级上、自旋为s的电子湮灭算符，γ表示谱函数

中的自变量；

流过量子点系统S中的电流表示为：

式中，P(n,t)表示在Δt时间内由n个电子穿过量子点系统S的概率，e表示单位电荷，n表示单位时间穿过量子点系统S的电子数目；

根据MacDonald公式，电流的散粒噪声谱表示为：

式中，ω表示散粒噪声S(ω)函数中的自变量。

6.根据权利要求5所述的循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，其特征在于：所述二能级量子电荷比特输运系统的密度矩阵在不同时刻之间的联系，利用Kraus算符即量子隐马尔可夫进行表示：即

式中m表示不同的K,Km是第m个Kraus算符；该公式与所述量子主方程等价；

而在量子条件主方程下，将ρ(t)＝∑_nρ⁽ⁿ⁾(t)带入上式就有

根据量子条件主方程的改写方程，ρ⁽ⁿ⁾(t+Δt)和ρ⁽ⁿ⁾(t),ρ^(n-1)(t),ρ⁽ⁿ⁺¹⁾(t)相关，因此联合上式可知：

即该式与量子条件主方程的改写方程有关系，目标就是构建不含时的映射

7.根据权利要求6所述的循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，其特征在于：所述LSTM细胞的输入值x_t和输出值h_t的关系，由下式方程给出：

f_t＝σ(W_f·[h_t-1,x_t]+b_f)

i_t＝σ(W_i·[h_t-1,x_t]+b_i)

O_t＝σ(W_O·[h_t-1,x_t]+b_O)

h_t＝O_t×tanh(C_t)

式中，(W_f,W_i,W_c,W_o,b_f,b_i,b_c,b_o)作为所述参数(W,b)；f_t为遗忘门限层的输出，i_t和

为输入门限层的输出，O_i和h_i为输出门限层的输出。

8.根据权利要求7所述的循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，其特征在于：利用散粒噪声谱的数据来构建密度矩阵迹之间即t-1时刻和t时刻之间的联系，即构建

把Tr[ρ⁽ⁿ⁾(t)]和长短时记忆网络中的h_t参数相对应；

将长短时记忆网络中的参数(W,b)就可以充当

的作用。

9.根据权利要求5所述的循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，其特征在于：所述方法还包括以下步骤：

用Tr[ρ⁽ⁿ⁾(t)]对总电流的贡献度的方法来有效的截断n值，n表示量子条件主方程的改写方程中的粒子数；定义评估函数：

式中，M是数值实验上能取到的最大粒子数值，P_M对应为有M个电子流过量子点系统的概率值；

通过不断调整M的值，绘制出了E(M)随M的变化图像，根据所述变化图像来确定M的取值。

10.根据权利要求1所述的循环神经网络模拟量子输运过程中的量子条件主方程的模拟方法，其特征在于：所述量子输运过程中产生的散粒噪声谱的数据分为训练数据和测试数据；

利用训练数据来训练所述循环神经网络，得到训练数据的误差随迭代次数的第一关系；利用所述测试数据来测试所述循环神经网络，得到测试数据的误差随迭代次数的第二关系；

利用第一关系和第二关系确定长短时记忆网络模拟量子条件主方程的模拟效果。

Images