二极管器件SPICE电路模型对锥形离子通道 I- V特性曲线的模拟

doi:10.7503/cjcu20230384

摘要/Abstract

摘要：

离子整流性是纳米离子通道的一个重要特征, 具有整流性的离子通道体系也被称为纳米流体二极管. 本文比较了离子通道的泊松-能斯特-普朗克(PNP)方程组模型和固体半导体的扩散-漂移模型, 提出可以使用二极管器件的仿真电路模拟器(SPICE)电路模型对离子通道体系的电流-电压( I- V)曲线进行模拟.以锥形离子通道的PNP数值模型的计算结果为基础, 通过对这一体系进行讨论, 给出一个锥形离子通道的SPICE电路模型, 它可以较好地模拟 I- V特性曲线. 离子通道SPICE电路模型的建立可用于研究纳米流体二极管作为一个器件在电路中的应用.

关键词: 离子通道, 整流性, 纳米流体二极管, 泊松-能斯特-普朗克方程, 仿真电路模拟器电路模型

Abstract:

Ion rectification is a main property of nanochannel. The ion-channel system with ion current rectification is so-called nanofluidic diode. In this paper， Poisson-Nernst-Planck（PNP） equations model of ion-channel and diffusion-drift model of semiconductor diode were compared， and Simulation Program with Integrated Circuit Emphasis（SPICE） model of diode was suggested， which could be used for I- V characteristic curves simulation of ion-channel system. Comparing I- V curves of several kinds of ion-channel from reported works， that is found that besides the saturated current term of ideal PN junction diode SPICE model， the reverse saturated current term should be used to build SPICE model of ion-channel systems within the ‒2—+2 V experiment range. Usually， PNP equations numerical model has a good fit to conical ion-channel I- V curve experiment. So， firstly based on PNP equations model， a conical ion-channel 2D model with a 10 μm length， a 400 nm base diameter and a 10 nm tip diameter was built in COMSOL multiphysics software， ion-channel systems with a series of surface charge densities in 1， 10 and 100 mmol/L KCl aqua solutions were simulated. After discussed I- V curve properties of conical ion-channel system based on this model， especially composition of diffusion and electrophoresis currents， the suggested SPICE model with the reverse saturated current term of diode would be built. In the SPICE model of conical ion-channel system， junction potential and series resistance of diode could be estimated from PNP model simulation， reverse saturated current and grading coefficient of diode would be set to fit PNP model simulation. As a result， SPICE model fit PNP equations model well simulating I- V curves of conical ion-channel system. SPICE model of ion-channel system would be useful to research the property of nanofluidic diode as an electrical element in electric circuit.

Key words: Ion channel, Rectification, Nanofluidic diode, Poisson-Nernst-Planck equation, Simulation program with integrated circuit emphasis（SPICE） model

中图分类号:

O646

TrendMD:

江嘉乔, 鲁冰新, 肖天亮, 翟锦. 二极管器件SPICE电路模型对锥形离子通道 I- V特性曲线的模拟. 高等学校化学学报, 2023, 44(12): 20230384.

JIANG Jiaqiao, LU Bingxin, XIAO Tianliang, ZHAI Jin. Simulation of SPICE Model of Diode for Conical Ion-channel I- V Curve. Chem. J. Chinese Universities, 2023, 44(12): 20230384.

图/表 8

Fig.1 Geometry of the conical ion channel PNP model and its boundary conditionsLeft boundary（red） is assigned as applied voltage and bulk concentration， right boundary（blue） is assigned as ground potential and bulk concentration， boundaries（yellow） of ion-channel wall are assigned with a uniform surface charge density.

Fig.2 Simulated I⁃ V curve and its composition of diffusion and electrophoresis currents(A), electric voltage and ions concentrations distributions on the axis of channel system under ±2 V applied voltages of the 10 mmol/L bulk concentration, -5 mC/m 2 surface charge density conical ion channel PNP model(B)

Fig.3 Simulated I⁃ V curves(A) and corresponding electric voltage distributions on the axis of channel system(B) under 0 V applied voltages of the 10 mmol/L bulk concentration conical ion channel PNP model with different surface charge densities

Fig.4 Scheme of SPICE model for the conical ion channel(A), simulated I⁃ V curves of the 10 mmol/L bulk concentration conical ion channel SPICE model with different surface charge densities(B)

Table 1 Parameters for SPICE model and simulated results of PNP and SPICE models of the 10 mmol/L bulk concentration conical ion channel

Surface charge density/（mC·m ^-2）	Model	I_SR/μA	R_S/kΩ	M	V_J /mV	I/μA
Surface charge density/（mC·m ^-2）	Model	I_SR/μA	R_S/kΩ	M	V_J /mV	-2 V	-1 V	-0.5 V	0.5 V	1 V	2 V
-10	PNP	—	—	—	—	-0.803	-0.510	-0.348	1.577	4.202	11.34
	SPICE	0.041	158	0.6	13	-0.812	-0.531	-0.345	2.291	5.246	11.37
-5	PNP	—	—	—	—	-0.849	-0.585	-0.419	1.243	3.165	8.405
	SPICE	0.054	215	0.5	7.4	-0.848	-0.589	-0.402	1.792	3.972	8.476
-2	PNP	—	—	—	—	-1.234	-0.890	-0.591	0.982	2.321	5.791
	SPICE	0.074	320	0.45	3.1	-1.233	-0.862	-0.588	1.305	2.784	5.818
-1	PNP	—	—	—	—	-1.773	-1.178	-0.684	0.901	1.972	4.605
	SPICE	0.088	410	0.45	1.6	-1.774	-1.183	-0.746	1.069	2.234	4.609

Fig.5 Simulated I⁃ V curves of PNP model(A), corresponding electric voltage distribution on the axis of channel system(B) under 0 V applied voltages, simulated I⁃ V curves of SPICE model(C) of the 1 mmol/L bulk concentration conical ion channel with different surface charge densities, simulated I⁃ V curves of PNP model(D), corresponding electric voltage distributions on the axis of channel system(E) under 0 V applied voltages, simulated I⁃ V curves of SPICE model(F) of the 100 mmol/L bulk concentration conical ion channel with different surface charge densities

Table 2 Parameters for SPICE model and simulated results of PNP and SPICE models of the 1 mmol/L bulk concentration conical ion channel

Surface charge density/（mC·m ^-2）	Model	I_SR/μA	R_S/kΩ	M	V_J /mV	I/μA
Surface charge density/（mC·m ^-2）	Model	I_SR/μA	R_S/kΩ	M	V_J /mV	-2 V	-1 V	-0.5 V	0.5 V	1 V	2 V
-1	PNP	—	—	—	—	-0.0778	-0.0484	-0.0325	0.1556	0.4094	1.0660
	SPICE	0.0050	1650	0.6	19.5	-0.0778	-0.0510	-0.0332	0.2187	0.5017	1.0880
-0.5	PNP	—	—	—	—	-0.0857	-0.0583	-0.0414	0.1221	0.3081	0.7820
	SPICE	0.0065	2350	0.5	10.5	-0.0853	-0.0592	-0.0404	0.1651	0.3649	0.7770
-0.2	PNP	—	—	—	—	-0.1285	-0.0912	-0.0596	0.0974	0.2266	0.5418
	SPICE	0.0091	3450	0.45	4.3	-0.1284	-0.0892	-0.0603	0.1221	0.2595	0.5411
-0.1	PNP	—	—	—	—	-0.1884	-0.1200	-0.0687	0.0882	0.1938	0.4378
	SPICE	0.0111	4350	0.45	2.2	-0.1881	-0.1233	-0.0760	0.1017	0.2116	0.4357

Table 3 Parameters for SPICE model and simulated results of PNP and SPICE models of the 100 mmol/L bulk concentration conical ion channel

Surface charge density/（mC·m ^-2）	Model	I_SR/μA	R_S/kΩ	M	V_J /mV	I/μA
Surface charge density/（mC·m ^-2）	Model	I_SR/μA	R_S/kΩ	M	V_J /mV	-2 V	-1 V	-0.5 V	0.5 V	1 V	2 V
-10	PNP	—	—	—	—	-17.80	-11.88	-6.891	8.840	19.55	45.53
	SPICE	0.177	43.0	0.5	0.125	-17.83	-11.50	-7.096	10.610	22.24	45.99
-5	PNP	—	—	—	—	-23.50	-13.70	-7.354	8.426	17.72	39.14
	SPICE	0.205	48.0	0.5	0.066	-23.54	-14.20	-7.918	9.383	19.39	39.75
-2	PNP	—	—	—	—	-28.11	-14.88	-7.643	8.066	16.52	34.63
	SPICE	0.350	55.5	0.45	0.026	-28.06	-15.55	-7.931	8.400	17.13	34.80
-1	PNP	—	—	—	—	-29.74	-15.28	-7.744	7.959	16.11	33.03
	SPICE	0.365	59.5	0.45	0.013	-29.69	-15.39	-7.708	7.952	16.14	32.66

参考文献 31

1	Jiang L.， Bioinspired Intelligent Nanostructured Materials， Science Press， Beijing， 2015， 19—35
	江雷. 仿生智能纳米材料，北京：科学出版社， 2015， 19—35
2	Sisson A. L.， Shah M. R.， Bhosale S.， Matile S.， Chem. Soc. Rev.， 2006， 35， 1269—1286
3	Vlassiouk I.， Smirnov S.， Siwy Z.， ACS Nano， 2008， 2（8）， 1589—1602
4	Siwy Z.， Fuliński A.， Am. J. Phys.， 2004， 72（5）， 567—574
5	Siwy Z. S.， Adv. Funct. Mater.， 2006， 16（6）， 735—746
6	Siwy Z.， Fulinski A.， Phys. Rev. Lett.， 2002， 89， 198103
7	Fulinski A.， Kosinska I. D.， Siwy Z.， New J. Phys.， 2005， 7， 132
8	Woermann D.， Phys. Chem. Chem. Phys.， 2003， 5， 1853—1858
9	Cervera J.， Schiedt B.， Ramirez P.， Europhys. Lett.， 2005， 71（1）， 35—41
10	Daiguji H.， Oka Y.， Shirono K.， Nano Lett.， 2005， 5（11）， 2274—2280
11	Cheng L. J.， Guo L. J.， Chem. Soc. Rev.， 2010， 39， 923—938
12	Li D.， Principle of Electrical Chemistry， Beihang University Press， Beijing， 2008
	李荻. 电化学原理. 北京：北京航空航天大学出版社， 2008
13	Bard A. J.， Faulkner L. R.， Electrochemical Methods Fundamentals and Applications， 2nd Ed.， Translated by Shao Y. H.， Zhu G. Y.， Dong X. D.， Zhang B. L.， Chemical Industry Press， Beijing， 2013
	邵元华，朱果逸，董献堆，张柏林［译］. 电化学方法原理和应用. 第2版，北京：化学工业出版社， 2013
14	Huang K.， Han R.， Solid State Physics， Higher Education Press， Beijing， 1988， 344—368
	黄昆，韩汝奇. 固体物理学. 北京：高等教育出版社， 1988， 344—368
15	Kittel C.， Introduction to Solid State Physics， 8th Ed.， Translated by Xiang J. Z.， Wu X. H.， Chemical Industry Press， Beijing， 2011， 249—267
	项金钟，吴兴惠［译］. 固体物理导论. 第8版，北京：化学工业出版社， 2011， 249—267
16	Hu A.， Zhang W. Y.， Solid State Physics， 1st Ed.， Higher Education Press， Beijing， 2005， 244—266
	胡安，章维益. 固体物理学. 第1版，北京：高等教育出版社， 2005， 244—266
17	Jing Y. K.， Pan Z. H.， Computer Simulation， 2008， 25（9）， 316—320
	荆玉珂，潘志浩. 计算机仿真， 2008， 25（9）， 316—320
18	Harikesh P. C.， Yang C. Y.， Wu H. Y.， Zhang S.， Donahue M. J.， Caravaca A. S.， Huang J. D.， Olofsson P. S.， Berggren M.， Tu D.， Fabiano S.， Nature Mater.， 2023， 22， 242—248
19	Antognetti P.， Massobrio G.， Semiconductor Device Modeling with Spice， 2nd Ed.， McGraw⁃Hill， New York， 1993， 39—41
20	Xu Y.， Jiang J.， Meng Z.， Zhang Q.， Li X.， Zhai J.， J. Phys. Chem. C， 2016， 120， 17342—17347
21	Li R.， Sui X.， Li C.， Jiang J.， Zhai J.， Gao L.， ACS Appl. Mater. Interfaces， 2018， 10， 3241—3247
22	Xie G.， Tian W.， Wen L.， Xiao K.， Zhang Z.， Liu Q.， Hou G.， Li P.， Tian Y.， Jiang L.， Chem. Commun.， 2015， 51， 3135—3138
23	Wang Y.， Chen H.， Jiang J.， Zhai J.， Small， 2021， 17， 2102880
24	Zhang H.， Hou J.， Hu Y.， Wang P.， Ou R.， Jiang L. Liu J. Z.， Freeman B. D.， Hill A. J.， Wang H.， Sci. Adv.， 2018， 4， eaaq0066
25	Fu L.， Wang Y.， Jiang J.， Lu B.， Zhai J.， ACS Appl. Mater. Interfaces， 2021， 13， 35197—35206
26	Zhou D.， Meng Z. Y.， Zhang M. H.， Zhai J.， Acta Chim. Sinica， 2015， 73， 716
	周迪，孟哲一，张明辉，翟锦. 化学学报， 2015， 73， 716
27	Xu Y. L.， Meng Z. Y.， Zhai J.， Acta Chim. Sinica， 2016， 74， 538
	许阳蕾，孟哲一，翟锦. 化学学报， 2016， 74， 538
28	Xu Y.， Zhang M.， Tian T.， Shang Y.， Meng Z.， Jiang J.， Zhai J.， Wang Y.， NPG Asia Mater.， 2015， 7， e215
29	Zhang Q.， Zhang Z.， Zhou H.， Xie Z.， Wen L.， Liu Z.， Zhai J.， Diao X.， Nano Res.， 2017， 10（11）， 3715—3725
30	Xu Y.， Sui X.， Jiang J.， Zhai J.， Gao L.， Adv. Mater.， 2016， 28， 10780—10785
31	Wang X.， Chen Y.， Meng Z.， Zhang Q.， Zhai J.， J. Phys. Chem. C， 2018， 122， 24863—24870

[1]	陈伟鹏, 孔祥玉, 闻利平. 水凝胶基仿生离子通道及其智能离子传输[J]. 高等学校化学学报, 2023, 44(6): 20220772.
[2]	曹忠, 郭佩弦, 薛书蕾, 肖忠良, 谭淑琪, 亢璇, 冯泽猛, 印遇龙, 马天骥. 肽适体修饰纳米孔道对L⁃精氨酸的灵敏检测[J]. 高等学校化学学报, 2023, 44(6): 20220758.
[3]	何芮, 王晓娜, 王歆悦, 邱娅, 杨胜勇, 尹宗宁. 钙离子通道阻滞剂自乳化给药系统的形成机理[J]. 高等学校化学学报, 2015, 36(1): 131.
[4]	刘军,那万里,辛伟红,刘丽丹,于波,麻彤辉,杨红 . 20(S)-原人参二醇促进CFTR氯离子通道开放[J]. 高等学校化学学报, 2008, 29(4): 731.
[5]	郭继东, 邢巍, 陆天虹, 杨辉. 短杆菌肽在十八烷基硫醇自组装单层膜上的通道行为[J]. 高等学校化学学报, 2002, 23(5): 958.
[6]	韩晓军, 宋永海, 吴正岩, 张柏林, 汪尔康. 钙离子对支撑磷脂膜离子通道行为的诱导作用[J]. 高等学校化学学报, 2001, 22(8): 1298.