1、第41卷第6 期2024年12 月J.At.Mol.Phys.,2024,41:062005(12pp)原子与分子物理学报JOURNAL OF ATOMIC AND MOLECULAR PHYSICS导电石墨烯孔道内双电层结构的研究Vol.41 No.6Dec.2024董乐,蒋更平(武汉科技大学高性能钢铁材料及其应用省部共建协同创新中心国际钢铁研究院,武汉430 0 8 1)摘要:为有效开发和利用新能源,人们迫切需要高性能的超级电容器提供能量的存储和转换,在超级电容器中双电层结构扮演着关键性的角色:本文利用分子动力学方法通过建立开放的石墨烯纳米孔道(12 nm),研究了KCl溶液在纳米孔道内的
2、双电层结构,同时也比较了恒电量模拟(Q)和恒电势模拟法(U)下双电层结构的异同:结果表明在恒电势模拟法考虑了导电石墨烯壁的镜像作用使结果更符合实验中的材料系统,而石墨烯壁的镜像作用能额外吸附离子从而增强孔道内部的阴阳离子,这可能有助于电极电容的提升:通过对不同孔道高度的研究,本文发现水分子作为介电材料在水基超级电容器中发挥着决定性的作用,它能在很大程度上抵消不同离子和不同孔道高度下双电层的变化,从而在不同情况下获得了相似的电容。关键词:分子动力学模拟;超级电容器;恒电势模拟方法;石墨烯中图分类号:0 48 5Study onelectric double layers inside the c
3、onducting graphene nanopores文献标识码:AD0I:10.19855/j.1000-0364.2024.062005DONG Le,JIANG Geng-Ping(Collaborative Innovation Center for Advanced Steels,International Research Institute for Steel Technology,Wuhan 430081,China)Abstract:In order to effectively develop and utilize new energy,people urgently
4、need high-performance super-capacitors to provide energy storage and conversion.Electric double layer structure plays a key role in superca-pacitors.In this work,the molecular dynamics simulation was taken to an open graphene nanochannel(1 2nm),and the double-layers structure of KCl solution inside
5、the nanochannel was studied.Meanwhile,the sim-ilarities and diferences of the double-layer structure under fixed charge simulation(Q)and constant potentialsimulation(U)were compared.The results show that the mirror image effect of the conducting graphene wall isconsidered in the constant potential s
6、imulation method,which makes the results more consistent with the practicalmaterials.The mirror image effect can absorb extra ions.Thus it enhances the concentratiosn of anion and cationin the pore,which may contribute to the improvement of electrode capacitance.Through the study of differentchannel
7、 heights,it is found that water molecules play a decisive role as dielectric materials in water-based su-percapacitor.It can largely offset the change of electric double layer under different ions and different channelheights,resulting in a similar capacitance under dfferent conditions.Key words:Mol
8、ecular dynamics simulation;Supercapacitor;Constant-potential simulation method;Graphenepacitor)作为超级电容器的一种,是一种新型储能1 引 言装置,其充放电过程完全没有涉及到化学变化,双电层电容器(Electrical Double-Layer Ca-超级电容器储能具有充放电速度快、功率密度大、收稿日期:2 0 2 3-0 3-2 9基金项目:国家自然科学基金青年科学基金项目(2 190 52 15)作者简介:董乐(1998 一),女,湖北武汉人,硕士研究生,研究方向为分子动力学模拟.E-mail:2
9、 6 36 6 6 2 7 13 q q.c o m通讯作者:蒋更平.E-mail:G e n g p i n g j i a n g w u s t.e d u.c n062005-1第41卷循环稳定性优异、清洁环保等优势,因此得到人们的广泛关注:石墨烯作为最早被发现的二维材料,是一种单层的蜂窝状的碳晶体,拥有许多优良的物理、化学、电学、力学性质2-4。通过对石墨烯内部结构的研究,发现表面化学和溶液成分的调控能实现比传统材料高得多的比电容和更稳定的循环性能5。Huang等人 提出了一个在小微孔中形成的模型,即electrici n c y l i n d e r 模型,有别于传统的双电层模型
10、Jiang等人7 通过分子动力学模拟研究发现较大的离子在较强的外加电场作用下会去溶剂化并穿透水层直接吸附于电极表面形成亥姆霍兹层Paek等人8 通过对碳纳米管和石墨烯表面双电层结构的研究揭示了电极外表面曲率半径对实验结果的影响Kalluri 等人9 展示了水合离子在纳米孔道约束下的脱水作用(h1.3nm时,孔道内同离子升高的幅度,与对离子降低的幅度相同这意味着在这一范围内,孔道的溶液以离子交换的形式极化2 9 即,以一个阳离子置换一个阴离子的形式来平衡表面电荷而当h1.3nm),在h=1.8nm内孔道内,得到了反常的离子分布如U=+0.5V时,孔道内的阳离子浓062005-10-61.2-4-
11、2h(A)1.41.6h(nm)02468101.82.0第41卷度基本不变(0.6 99M)溶液的极化作用仅由阴离子的注人完成的:如果进一步增大阳极的电压(U=+1.0 V),则阳离子浓度降低的幅度(1.84 M)与阴离子升高的幅度(0.2 2 M)相匹配.也就是说在1.8 nm的阳极孔道中,首先以对离子注入完了溶液的极化,然后伴随有阴阳离子的交换过程类似的反常现象也可以在1.8 nm的阴极孔道中观察得到当U=0.5V或-1.0 V时,孔道内部的CI-浓度基本不变,始终维持在比较高的浓度(0.48 9)这一反常结果可能是应为石墨烯壁上双电层结构干涉相作用和镜像作用所共同造成的.通过以上的恒电
12、势模拟,本文也计算了不同孔道内部双电层结构的电容,即C=Q/AU.Q为电极上平均的电荷数目在恒电势模拟中U为阴阳两电极间的电势差,这一计算得到的电容是阴阳两电极所组成的电容器的整体电容,表1为电压为U=0.5V与1.0 V时的电容我们发现这一结果与此前Jiang7的文章中恒电荷条件下计算得到的结果相差不大如U=0.5V或U=1.0V时,比电容大小为1.41 1.44F/g.而Jiang7的文章中当Q=0.041 e时,KCl 溶液的比电容也为1.44 F/g.表1不同孔道宽度下电容器的质量比电容Table 1 The specific capacitances of nanochannel w
13、ith thedifferent channel heights孔道宽度(nm)0.5V11441.21431.51451.81412147比较不同孔道宽度时,计算得到的双电层电容却几乎保持不变尽管在不同孔道高度下内部的双电层有着截然不同的分布结构(离子峰的分布和高度)但最终的电极电容却几乎一致相较于离子,中性的水分子的数目在溶液中占主导地位(1:55.4)因此,水分子对超级电容器的储能性能的起到关键作用7 不同离子分布导致的电容变化,反过来被作为电介质物质的水分子所屏蔽这最终导致了类似的电极电容。董乐,等:导电石墨烯孔道内双电层结构的研究质量比电容(F/g)1.0V1411441501431
14、44第6 期4 结 论本章使用分子动力学模拟研究了纳米尺度下的石墨烯孔道内的双电层结构,采用恒电量模拟(FCM)与恒电势模拟(CPM)两种不同方法分别对不同孔道高度的下的石墨烯孔道进行了充电与不充电的模拟,研究了不同孔道宽度、电量与电极电压下水分子与阴阳离子的在孔道内的分布,主要得到以下结论:1.在不充电情况下,对比了恒电量(Q=0)与恒电势(U=0)模拟下孔道内部的溶液分布,两者有着类似的水分子的分布结构但恒电势模拟下,极板具有镜像作用这增强了孔道内部的离子浓度。2随着孔道高度的降低,离子进人孔道内部需要脱去周围的水合结构这导致了内部的离子浓度低于孔道外的溶液区域(1.3 1.8 M)同时随
15、着孔道高度的降低,孔道的约束作用强烈,这导致了在入口处更高的能垒,这降低了孔道内的离子浓度,随着电极电压的上升,在大孔径内(h 1.3nm)溶液的极化以离子交换为主:而当h1.3nm,极化的主要方式是对离子的注人:3通过计算不同孔道高度下的电容,发现电极电容基本保持不变,因为水分子是影响电容的重要因素对于在不同孔道内由不同离子分布而产生的变化,会反过来被水分所屏蔽,这导致了类似的电极电容,参考文献:1Chen J,Wei H,Fu N,et al.Facile synthesis of nitro-gen-containing porous carbon as electrode materi
16、alsfor superior-performance electrical double-layer ca-pacitorsJJ.J.Mater.Sci.,2018,53:2137.2Huang G R,Chen J.Synthesis and applications of gra-phene J.Carbon Tech.,2009,28:35(in Chinese)黄桂荣,陈建石墨烯的合成与应用J炭素技术,2 0 0 9,2 8:353Huang G J,Chen Z G,Li M D,et al.Surface func-tionalization modification of gr
17、aphene and graphene ox-ide J.Acta Chim.Sin.,2016,74:789(i n Ch i-nese)【黄国家,陈志刚,李茂东,等石墨烯和氧化石墨烯的表面功能化改性J化学学报,2 0 16,74:7894Hu Y J,Jin J,Zhang H,et al.Preparation,function-alization and application of graphene in chemistry J.062005-11第41卷Acta Phys.-Chim.Sin.,2010,26:2073(i n Ch i-nese)【胡耀娟,金娟张卉,等石墨烯的制备
18、、功能化及在化学中的应用J物理化学学报,2 0 10,26:20735 Yaseen M,Khattak M A K,Humayun M,et al.AReview of supercapacitors:materials design,modifica-tion,and applicationsJ.Energies,2021,14:7779.6Huang J,Sumpter B G,Meunier V.Theoretical modelfor nanoporous carbon supercapacitors J.Angew.Chem.Int.Ed.Engl.,2008,120:530.7J
19、iang G,Cheng C,Li D,et al.Molecular dynamicssimulations of the electric double layer capacitance ofgraphene electrodes in mono-valent aqueous electro-lytesJ.Nano Res.,2016,9:174.8Paek E,Pak A J,Hwang G S.Curvature effects onthe interfacial capacitance of carbon nanotubes in an i-onic liquidJ.J.Phys.
20、Chem.C,2015,117:23539.9Kalluri R K,Biener M M,Suss M E,et al.Unravelingthe potential and pore-size dependent capacitance ofslit-shaped graphitic carbon pores in aqueous electro-lytesJ.Phys.Chem.Chem.Phys.,2013,15:2309.10 Zhang Y,Cummings P T.Effects of solvent concentra-tion on the performance of io
21、nic-liquid/carbon super-capacitors J.ACS Appl.Mater.Inter.,2019,11:42680.11 Merlet C,Pan C,Rotenberg B,et al.Highly con-fined ions store charge more efficiently in supercapaci-torsJ.Nat.Commun.,2013,4:2701.12 Yin L,Huang Y,Chen H,et al.A mean-field theo-ry on the differential capacitance of asymmetr
22、ic ionicliquid electrolytes.II.Accounts of ionic interactionsJ.Phys.Chem.Chem.Phys.,2018,20:17606.13 Wu Y C,Ye J,Jiang G,et al.Electrochemical char-acterization of single layer graphene/electrolyte inter-face:effect of solvent on the interfacial capacitanceJ.Angew.Chem.Int.Ed.Engl.,2021,60:13317.14
23、Gu C,Yin L,Li S,et al.Differential capacitance ofionic liquid and mixture with organic solvent J.Elec-trochim.Acta,2021,367:137517.15 Li S,Zhao Z,Liu X.Electric double layer structureand capacitance of imidazolium-based ionic liquidswith FSI-and Tf anions at graphite electrode by mo-lecular dynamic
24、simulations J.J.Electroanal.Chem.,2019,851:113452.16 SShalabi A,Daniels L,Scott M,et al.Differential ca-pacitance of ionic liquid interface with graphene:Theeffects of correlation and finite size of ions J.Elec-原子与分子物理学报trochim.Acta,2019,319:423.17 Bossa G V,Caetano D L Z,de Carvalho S J,et al.Diffe
25、rential capacitance of an electrical double layerwith asymmetric ion sizes in the presence of hydrationinteractionsJ.Electrochim.Acta,2019,321:134655.18 Feng G,Qiao R,Huang J,et al.Atomistic insight onthe charging energetics in subnanometer pore superca-pacitorsJ.J.Phys.Chem.C,2010,114:18012.19 Zeng
26、 L,Wu T,Ye T,et al.Modeling galvanostaticcharge-discharge of nanoporous supercapacitors J.Nat.Comput.Sci.,2021,1:725.201Neto A J P,da Silva D A C,Goncalves V A,et al.An evaluation of the capacitive behavior of supercapaci-tors as a function of the radius of cations using simula-tions with a constant
27、 potential method J.Phys.Chem.Chem.Phys.,2022,24:3280.21Demir B,Searles D J.Investigation of the ionic liquidgraphene electric double layer in supercapacitors usingconstant potential simulations J.Nanomaterials,2020,10:2181.22 Bo Z,Li C,Yang H,et al.Design of supercapacitorelectrodes using molecular
28、 dynamics simulations J.Nanomicro.Lett.,2018,10:1.23 Mackerell A D,Bashford D,Bellott M,et al.All-at-om empirical potential for molecular modeling and dy-namics studies of proteins J.J.Phys.Chem.B,1998,102:3586.24 JJorgensen W L,Chandrasekhar J,et al.Comparison ofsimple potential functions for simul
29、ating liquid waterJ.J.Chem.Phys.,1983,79:926.25Hoover W G.Canonical dynamics:Equilibrium phase-space distributions J.Phys.Rev.A,1985,31:1695.26Hiroshi,Nakano,Hirofumi,et al.A chemical poten-tial equalization approach to constant potential polariz-able electrodes for electrochemical-cell simulationsJ
30、.J.Chem.Phys.,2019,151:164123.27Deng L J.Application of electrode polarization phe-nomenon J.University Chem.,2022,37:2107049(i n Ch i n e s e)【邓玲娟电极极化现象的应用J.大学化学,2 0 2 2,37:2 10 7 0 4928 Li P,Song L F,Merz K M.Systematic parameteriza-tion of monovalent ions employing the nonbonded modelJ.J.Chem.Theory Comput.,2015,1l:1645.29 Simon P,Gogotsi Y.Perspectives for electrochemicalcapacitors and related devices J.Nat.Mater.,2020,19:1151.062005-12第6 期