一种基于单分子和二分子化学反应网络实现M/N占空比时钟信号的方法与流程

本发明属于dna计算领域，是dna复杂设计的基础性前端。

背景技术：

2010年，davidsoloveichik等人提出“dnaasauniversalsubstrateforchemicalkinetics”，一定程度上从理论层面证明：对于任何一个形式化学反应网络(chemicalreactionnetworks，以下简称crns)，总能方便地找到其对应的dna的物理实现。也就是说，任意一个我们所设计的crns，可以映射到dna分子链置换反应上；而dna分子链置换反应，在“忽略反应物、生成物具体是什么dna分子”的前提下，可以抽象为一个由一系列形如a+b→c+d的基元反应所构成的crns，从而使设计简单化。

2013年，jianghua等人提出“digitallogicwithmolecularreactions”，对于单个比特如x，可以用式(a)这一组3个化学反应来表示，并由此提出了基于分子反应双稳态的数字逻辑，以及简单组合逻辑的设计。

另外，jianghua等人也在这篇文章中首次提出了一种化学反应时钟的设计方法，但是这种方法是建立在单分子，双分子和三分子化学反应基础上的，就目前的crns到dna映射方法来说三分子反应体系是较难实现其动力学特性的映射的，而且该时钟设计没有给出系统性的时钟参数修改方法和严格的数学论证。因此，我们现在针对这些问题和未来可能的复杂设计需求提出了一种新的设计方法。

在传统电子学中，时钟信号可以由晶体振荡器以及其他方式成功产生。本质上，所有的硅基设计方法都有相同的一个物理实质——时钟生成器必有共振电路和放大器这两个重要组成部分。具体来说，共振电路实际上是一个振荡器，意图能够不断地重复同样的操作。而放大器则通常被用来作为控制器在反馈环路中维持振荡效果。

最初由h.jiang等人于2010年提出，而后被l.ge,c.zhang等人详细阐述，一个四相化学振荡器，即rgby振荡器，主要用于在分子系统中生成化学时钟信号。其对应的四个时钟相信号则由分子类型r(ed)、g(reen)、b(lue)、y(ellow)表示。每个相时钟信号在一个周期中持续相同时间。整个振荡器共需16个形式化学反应构成。本质上，该rgby振荡器实现了相时钟信号的重复转化。

具体来说，在一个rgby振荡器中，我们称r、g、b、y为相时钟信号，而称r、g、b、y为不存在标识物种。以由r到g的相时钟信号转换为例，对应的化学反应对应于表1中的第二行。假设系统最开始的相时钟状态为r，随着外界不断地低速注入r，化学反应r+r→r意味着相信号吸收消耗着不存在标识物种，从而维持着时钟现态。一旦分子r被消耗且转换为g，时钟状态则发生改变。阈值反应r+y→g+y触发而主动力反应r+2g→3g则加速了从r到g的转变。

表1rgby振荡器所需的化学反应

在h.jiang,m.riedel,andk.parhi等人发表于2011年的文献"synchronoussequentialcomputationwithmolecularreactions,"中，符号φ表明“没有反应物”，意味着生成物从一个很大的可充能的源中生成。为便于ode(ordinarydifferentialequations，常微分方程组)分析，在下文中使用物种emp或者empi来代替φ，其中i是一个正整数。由于φ和emp或empi的浓度足够高而可被视为一个常数，这种替代意味着两者得到的不存在标识物种等效。

davidsoloveichik等人指出，只有由单分子或双分子构成的crns可以被“转译”到dna链置换反应当中。因此，在真实的dna场景中，形如表1中r+2g→3g这样的反应物有3个分子的反应，即三分子反应将会变得较难实现。而且，尽管这些反应可以被分解成反应物不超过两个的级联反应，其复杂度，尤其是反应速率的代价将会变得很大。

技术实现要素：

发明目的：针对三分子反应难以实现的问题，本发明公开了一种基于单分子和二分子化学反应网络实现m/n占空比时钟信号的方法，该方法只考虑单分子和双分子反应，具有易于实现的优点。

技术方案：一种基于单分子和二分子化学反应网络实现m/n占空比时钟信号的方法，包括如下步骤：

设m/n占空比时钟信号的状态为re0和re1，分别代表“低”和“高”逻辑状态，0时刻状态为re0；

如果n为2，构造化学反应网络生成1/2占空比时钟信号；

如果n为不等于2的偶数，循环执行步骤1，得到时钟状态物种为re0和re1的m/n占空比时钟信号：

步骤1、记录n个连续的1/2占空比时钟信号的状态物种转移次数，在所述n个连续的1/2占空比时钟信号中第n-m次和第n次由物种e0转移到e1时，m/n占空比时钟信号的状态进行re0和re1之间的转移；由此得到一个周期的m/n时钟信号，其中e0和e1为所述1/2占空比时钟信号的时钟状态物种；

如果n为奇数，循环执行步骤2，得到时钟状态物种为re0和re1的m/n占空比时钟信号：

步骤2、记录2n个连续的1/2占空比时钟信号的状态物种转移次数，在所述2n个连续的1/2占空比时钟信号中第n-m次、第n次、第2n-m次和第2n次由物种e0转移到e1时，m/n占空比时钟信号的状态进行re0和re1之间的转移；由此得到2个周期的m/n时钟信号。

所述构造化学反应网络生成1/2占空比时钟信号包括如下步骤：

设1/2占空比时钟信号的时钟状态物种为e0和e1；初始时分子系统中有状态物种为e0；

外部环境连续低速注入不存在标识物种empi，i＝1,2,3，所述不存在标识物种在分子系统中不同物种的作用下被消耗，化学反应方程如下：

在分子系统中有物种e0存在时：

在分子系统中有物种e1存在时：

在分子系统中有物种m0存在时：

在分子系统中有物种m1存在时：

在分子系统中有物种f0存在时：

在分子系统中有物种f1存在时：

不同物种浓度发生变化，生成了中间物种e0s和e1s，中间物种与时钟状态物种循环往复转移，化学反应方程如下：

中间物种与时钟状态物种消耗中间反应分子，化学反应方程如下：

强制信号物种和时钟信号物种之间进行转移，化学反应方程如下：

其中，e0，e1，f0，f1，m0，m1为中间反应分子；e0t，e1t，f0t，f1t，m0t，m1t为转移指示物种；e0s和e1s为中间物种；f0，f1，m0和m1为起延时作用的强制信号物种；

由式(1)-(9)中的32个化学反应组成的化学反应网络，使分子系统在物种e0和e1之间转移，得到1/2占空比时钟信号。

所述实现1/2占空比时钟信号的化学反应网络可以借助方图实现，方图中圆形节点代表一个物种，每个圆形节点产生或接收各种箭头；直线或环形箭头表示物种数量传递流，多边形折线型箭头表示外部方形循环和内部菱形循环之间的信息传递。

为了控制各物种的维持时间，可以给每个化学反应方程设置不同的速率常数。

优选地，式(1)-(9)中，下列化学反应的速率常数设置为10：

下列化学反应的速率常数为20：

下列化学反应的速率常数设置为2.09：

其他化学反应速率常数则设置为1000。

有益效果：与现有技术相比，本发明公开的方法具有以下优点：1、本发明公开的方法提供了一个形式化学反应网络，它实现了1/2占空比的时钟信号的功能，该方法能够更方便的用dna链置换反应实现，而无需进行局部再设计；2、本发明公开的方法只考虑单分子和双分子反应，不涉及实现难度较大的三分子反应，由此，本发明具有易于实现的优点；3、现有的利用化学反应生成m/n占空比时钟信号的方法是将rgby这四相振荡器推广到n相振荡器，用生成的1/2占空比的时钟信号分别控制n相振荡器的第m相和第n相时钟信号，从而实现m/n占空比；而本发明的方法则无需借助n相振荡器，仅仅是用到了1/2占空比时钟的两个相信号e0和e1，而没有再去关注e0，e1是如何生成的。因此，只要保证在反应环境中m/n占空比时钟设计引入的物质不影响1/2占空比的时钟，那么就可以使用其他方法生成的1/2占空比时钟作为基础。

附图说明

图1为本发明实施例一实现1/2占空比时钟信号的方图；

图2为本发明实施例一实现1/2占空比时钟信号的仿真结果；

图3为当n为偶数且n≠2时实现m/n占空比时钟信号的一般性方法；

图4为实施例二中实现3/4占空比时钟信号的crns和仿真结果；

图5为实施例三中实现5/7占空比时钟信号的crns和仿真结果。

具体实施方式

下面结合附图和具体实施方式，进一步阐明本发明。

实施例一：

本实施例为实现1/2占空比时钟信号的具体步骤。

为实现1/2占空比的时钟信号，总共需要32个化学反应，化学反应方程如式(1)-(9)。设1/2占空比时钟信号的时钟状态物种为e0和e1，e0s和e1s为中间物种，整个分子系统状态转移是四种物种按照e0→e0s→e1→e1s→e0的形式循环往复。强制信号物种f0，f1，m0和m1的产生可以被视为延迟单元，而随后触发的化学反应显示了时钟信号状态物种的转化。

实现1/2占空比时钟信号的化学反应网络可以借助方图实现，如图1所示。方图实际上揭示了在循环中物种转移的轨迹。其中，每个圆形节点代表一个物种，且该节点产生或接收各种箭头。直线或环形箭头暗示着物种数量传递流，而多边形折线型箭头表示外部方形循环和内部菱形循环之间的信息传递。外部环境连续低速注入不存在标识物种empi，i＝1,2,3；内部菱形循环表示相时钟状态的循环变化，外部方形循环来控制相时钟状态变化的延迟，同时外部方形循环也是受到内部菱形循环相时钟的控制，两者是交互控制的。

以e0节点为例，共生成总共3个箭头：1)箭头e0-1：朝向e0s的直线；2)箭头e0-2：从e0开始并返回到e0的环形曲线；3)箭头e0-3：朝向m1的箭头1的多边形折线。

有了箭头e0-2，存在和不存在的逻辑关系得以确立：外部环境连续低速注入不存在标识物种emp1，该不存在标识物种只要时钟信号物种e0存在就被快速消耗，否则就会缓慢积累到分子体系当中，对应的化学反应方程为：

e0所产生的箭头e0-3发送消息以便可低速转化m1为f0；化学反应方程为：

e0+m1→e0+f0

m1节点生成2个箭头：m1-1和m1-2。m1的箭头m1-2表示m1和m1t的消耗，则将从f0到m0的转化延长了τd时间；化学反应方程如下：

然后，f0快速吸收f0，迫使e0保持在峰值长达τd，化学反应方程如下：

如上所述，在几乎所有m1转化为f0后，m1缓慢增多，使f0存在一小段时间后急剧转化为m0；化学反应方程如下：

一旦f0消失，e0同样也在一个很小的τt1时间段内转化为e0s；化学反应方程如下：

e0s在e1出现之前快速消耗了e1的转移指示物种f1t；化学反应方程如下：

e0s+f1t→e0s

e0s不会持续存在很长时间。当e0增长时，e0s在很短的时间τt2内转化为e1；化学反应方程如下：

m0则为e1等待了τd的时间。其他物种的转化以类似的方式运作。

数值仿真结果如图2所示，图2-(a)是只含e0和e1的波形图；图2-(b)是包括e0、e1、m0、m1和f0的波形图。

上述时长τd、τt1和τt2可以通过给每个化学反应方程设置不同的速率常数来控制。在本实施例中，所有类似empi→empi+resultant的化学反应，其速率常数都设置为10；另外四个反应：2f0→2f0+f0t、2f1→2f1+f1t、2m0→2m0+m0t和2m1→2m1+m1t，其速率常数也设置为10。而2e0→2e0+e0t和2e1→2e1+e1t的速率常数为20。此外，化学反应e0+m1→e0+f0和e1+m0→e1+f1的速率常数设置为2.09。其他化学反应速率常数则设置为1000。

上述速率常数也可以根据实际情况来人为调控。

实施例二：

对于n≠2的m/n占空比的时钟信号，其实现是基于1/2占空比时钟信号。设时钟状态物种为re0和re1，分别代表“低”和“高”逻辑状态。m/n占空比的时钟信号实现方法的核心思想是“记录”1/2占空比相时钟信号中e0和e1转化的次数，即用一组反应物ctri来记录每次时钟信号的转变，其中i是从1到n的整数。即整个1/2占空比时钟信号将被“打包”成段。每一段都包含了e0和e1之间的n次转化，其指示占空比的分母n。分子m的实现是通过时钟状态物种re0和re1构建的，即将n次转化拆分成两部分，一个含有n-m次转化，而另一个则含m次。因此，目标m/n占空比时钟信号可以是由时钟状态物种re0和re1来实现。值得注意的是，n是偶数或奇数时情况会有所不同。

表2实现n为偶数(n≠2)时的m/n占空比时钟信号的crns步骤

本实施例为n为偶数时的m/n占空比时钟信号。如前面所述，分母n可以通过反应物ctri,(i＝1,2,..,n)来实现。当设置分子m时，则使用这两个反应：一个是在第n次转移时，起记录作用的反应物ctrn将时钟状态从re1改变为re0，而另一反应则是在第n-m次转移中由ctrn-m触发，将时钟状态物种由re0变为re1，整个方案如图3所示。获得相应的crns的具体步骤见表2。

表2中，步骤11需要总共n-1个化学反应，步骤12只需一个反应，而步骤13则需两个。因此当n为偶数时，除了1/2个占空比时钟信号所需的化学反应之外，总共还需要n+2个反应来实现m/n占空比时钟信号。

以3/4占空比时钟信号为例，需要记录1/2占空比时钟信号状态物种的4次转移。因此可以用公式(11)所示的化学反应方程来刻画图4(a)所示的从ctr1→ctr2→ctr3→ctr4的循环转移走向。由于公式(12)中的化学反应，目标时钟状态re0和re1在反应物ctr1和ctr4的指导下分别进行第一和第四次状态转化。仿真结果如图4(b)所示，图中虚线11和21为e0和e1的浓度上下界。

实施例三

根据n为偶数的方法，很自然地会对n是奇数时记录n次原始1/2占空比时钟信号状态的转移。但是，两者在选择“记录反应物”来改变目标时钟状态这一方面很不一样。其主要原因在于，如果仍沿用n为偶数的方法，那么下一个循环(从ctr1到ctrn)将不会进行。针对这个问题，对n是奇数的情况，本发明公开了一个可行的解决办法，即开展一个新的循环，该循环包含两次从ctr1到ctrn的基本循环。因此这新循环可以“记录”1/2时钟信号的2n次转移。触发时钟状态发生变化的反应物包括ctrn-m、ctrn、ctr2n-m和ctr2n。获得相应的crns的具体步骤见表3。

表3当n为奇数时的m/n占空比时钟信号的设计

表3中，除了1/2占空比时钟信号所需的化学反应之外，该方法还需2n+4个反应。

以5/7占空比的时钟信号为例。根据如图5(a)所示的14态环路，可以得到如公式(13)所列的化学反应方程。相应的仿真结果如图5(b)所示，图中虚线11和21为e0和e1的浓度上下界。

在实施例二和实施例三中，仅仅是用到了1/2占空比时钟的两个时钟信号状态物种e0和e1，而没有再去关注e0和e1是如何生成的。因此，只要保证在反应环境中m/n占空比时钟信号的实现过程引入的物质不影响1/2占空比的时钟信号的实现，那么就可以使用其他方法生成的1/2占空比时钟作为基础。