JP6341922B2 - アニオン−カチオンペアを脱溶媒和しそれを電気二重層キャパシタの電極に挿入するためのトータルエネルギーを決定する方法 - Google Patents
アニオン−カチオンペアを脱溶媒和しそれを電気二重層キャパシタの電極に挿入するためのトータルエネルギーを決定する方法 Download PDFInfo
- Publication number
- JP6341922B2 JP6341922B2 JP2015540192A JP2015540192A JP6341922B2 JP 6341922 B2 JP6341922 B2 JP 6341922B2 JP 2015540192 A JP2015540192 A JP 2015540192A JP 2015540192 A JP2015540192 A JP 2015540192A JP 6341922 B2 JP6341922 B2 JP 6341922B2
- Authority
- JP
- Japan
- Prior art keywords
- energy
- anion
- cation
- determining
- pair
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 54
- 239000003990 capacitor Substances 0.000 title claims description 39
- 150000002500 ions Chemical group 0.000 claims description 60
- 238000001179 sorption measurement Methods 0.000 claims description 43
- 238000004088 simulation Methods 0.000 claims description 36
- 150000001768 cations Chemical class 0.000 claims description 34
- 150000001450 anions Chemical class 0.000 claims description 33
- 239000002071 nanotube Substances 0.000 claims description 31
- 239000011148 porous material Substances 0.000 claims description 30
- 239000003463 adsorbent Substances 0.000 claims description 18
- 230000009881 electrostatic interaction Effects 0.000 claims description 16
- 238000004807 desolvation Methods 0.000 claims description 12
- 239000002156 adsorbate Substances 0.000 claims description 11
- 238000000926 separation method Methods 0.000 claims description 10
- 238000012216 screening Methods 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 5
- 238000012900 molecular simulation Methods 0.000 claims description 5
- 239000002904 solvent Substances 0.000 claims description 5
- 238000000324 molecular mechanic Methods 0.000 claims description 4
- 102000004190 Enzymes Human genes 0.000 claims description 3
- 108090000790 Enzymes Proteins 0.000 claims description 3
- 229910021536 Zeolite Inorganic materials 0.000 claims description 3
- HNPSIPDUKPIQMN-UHFFFAOYSA-N dioxosilane;oxo(oxoalumanyloxy)alumane Chemical group O=[Si]=O.O=[Al]O[Al]=O HNPSIPDUKPIQMN-UHFFFAOYSA-N 0.000 claims description 3
- 102000004169 proteins and genes Human genes 0.000 claims description 3
- 108090000623 proteins and genes Proteins 0.000 claims description 3
- 239000010457 zeolite Substances 0.000 claims description 3
- 125000002057 carboxymethyl group Chemical group [H]OC(=O)C([H])([H])[*] 0.000 claims 1
- 230000003993 interaction Effects 0.000 description 23
- 238000005411 Van der Waals force Methods 0.000 description 5
- 239000003792 electrolyte Substances 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 3
- 238000013459 approach Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- -1 electrode Chemical class 0.000 description 3
- 238000005381 potential energy Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 238000000329 molecular dynamics simulation Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- UINDRJHZBAGQFD-UHFFFAOYSA-O 2-ethyl-3-methyl-1h-imidazol-3-ium Chemical compound CCC1=[NH+]C=CN1C UINDRJHZBAGQFD-UHFFFAOYSA-O 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 238000000418 atomic force spectrum Methods 0.000 description 1
- 229910001423 beryllium ion Inorganic materials 0.000 description 1
- 239000002041 carbon nanotube Substances 0.000 description 1
- 229910021393 carbon nanotube Inorganic materials 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 150000003949 imides Chemical class 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000000302 molecular modelling Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000027756 respiratory electron transport chain Effects 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- CBXCPBUEXACCNR-UHFFFAOYSA-N tetraethylammonium Chemical compound CC[N+](CC)(CC)CC CBXCPBUEXACCNR-UHFFFAOYSA-N 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 125000001889 triflyl group Chemical group FC(F)(F)S(*)(=O)=O 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01G—CAPACITORS; CAPACITORS, RECTIFIERS, DETECTORS, SWITCHING DEVICES, LIGHT-SENSITIVE OR TEMPERATURE-SENSITIVE DEVICES OF THE ELECTROLYTIC TYPE
- H01G11/00—Hybrid capacitors, i.e. capacitors having different positive and negative electrodes; Electric double-layer [EDL] capacitors; Processes for the manufacture thereof or of parts thereof
- H01G11/22—Electrodes
- H01G11/26—Electrodes characterised by their structure, e.g. multi-layered, porosity or surface features
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01G—CAPACITORS; CAPACITORS, RECTIFIERS, DETECTORS, SWITCHING DEVICES, LIGHT-SENSITIVE OR TEMPERATURE-SENSITIVE DEVICES OF THE ELECTROLYTIC TYPE
- H01G11/00—Hybrid capacitors, i.e. capacitors having different positive and negative electrodes; Electric double-layer [EDL] capacitors; Processes for the manufacture thereof or of parts thereof
- H01G11/22—Electrodes
- H01G11/30—Electrodes characterised by their material
- H01G11/32—Carbon-based
- H01G11/36—Nanostructures, e.g. nanofibres, nanotubes or fullerenes
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01G—CAPACITORS; CAPACITORS, RECTIFIERS, DETECTORS, SWITCHING DEVICES, LIGHT-SENSITIVE OR TEMPERATURE-SENSITIVE DEVICES OF THE ELECTROLYTIC TYPE
- H01G11/00—Hybrid capacitors, i.e. capacitors having different positive and negative electrodes; Electric double-layer [EDL] capacitors; Processes for the manufacture thereof or of parts thereof
- H01G11/54—Electrolytes
- H01G11/58—Liquid electrolytes
- H01G11/62—Liquid electrolytes characterised by the solute, e.g. salts, anions or cations therein
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/13—Energy storage using capacitors
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Electric Double-Layer Capacitors Or The Like (AREA)
- Chemical & Material Sciences (AREA)
- Power Engineering (AREA)
- Fixed Capacitors And Capacitor Manufacturing Machines (AREA)
- Crystallography & Structural Chemistry (AREA)
- Nanotechnology (AREA)
- Materials Engineering (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
本発明は、電気二重層キャパシタの分野に関し、より一般的には、電気によって帯電された吸着体と電気によって帯電された吸着物との間の静電相互作用をシミュレーションする方法に関する。
特に、本発明は、イオン性溶液のアニオン−カチオンペアを脱溶媒和しこのペアを電気二重層キャパシタの電極の中に挿入するためのトータルエネルギーを決定する方法である。本発明は、電気二重層キャパシタのイオン性溶液の材料／カチオン−アニオンのペアのスクリーニングも可能にする。実際、このトータルエネルギーは、電気二重層キャパシタの静電容量に関連し、よいスクリーニング基準である。
電気エネルギー蓄積は、ハイブリッドおよび／または電気自動車の発達の状況において、特に重要な面である。電気二重層キャパシタ（または、電気化学キャパシタ）は、その（非常に）高いパワーにおいて、従来のキャパシタと異なる。それらは、１分のオーダの非常に短い時間にわたる複数のエネルギーパルスを必要とするアプリケーションに特によく適合する。
電気二重層キャパシタは、２つの多孔性の電極で構成され、この電極は、一般的には活性炭で作られていて電解液で充満されており、この２つの電極は、浸透性の隔離膜（イオン導電性を有する）によって分けられている。電気二重層は、各電極−電解質インタフェースで発達し、電気二重層キャパシタは、２つのキャパシタの直列接続として模式的に表すことができ、その２つのうちの一方が陽極を有し、他方が陰極を有する。キャパシタの接続の法則によれば、直列接続の静電容量は、２つのキャパシタのうちの低い方よりも常に低い。
バッテリと異なり、電子伝達を招く化学反応は電気二重層キャパシタでは生じず、電気二重層キャパシタは物理現象によって支配される。
充電期間の間、電位が２つの電極に印加され、電解質内のイオンは分離されて電気化学二重層を形成し、陰極が陽イオンを引きつけ、陽極が陰イオンを引きつける。カレントコレクタが使用されるとすぐに、キャパシタは放電し、二重層は壊れる。
キャパシタの静電容量（C）は、
に従って、蓄積された電荷の量（Q）と印加された電位差（ΔV）によって決定される。
電圧差は、
の電荷の関数として表現される。
ここで、Eは電界の値、εは媒体の誘電率、Aは電極／電解液インタフェースの表面積、そして、dは二重層の厚さである。
ここで、Eは電界の値、εは媒体の誘電率、Aは電極／電解液インタフェースの表面積、そして、dは二重層の厚さである。
我々は、このようにして、
を有する。
電極においては、その静電容量は、インタフェースの表面積が増加したときおよび二重層の厚さが減少したときに増加する。その厚さは、イオンが溶媒分子（水、アセトニトリン）または反対イオンに奪われたときに減少するであろう。イオン性溶液族（１００℃より低い温度の塩類溶液）は、それらの反対イオンに奪われる可能性／性質を有する。
イオン性溶液を使う電気二重層キャパシタの静電容量を最適化するために、イオン性溶液の多くの材料／カチオン−アニオンペアがテストされる必要がある。
そのようなスクリーニングを実行する方法は、ペアごとに、イオン性溶液（アニオン＋カチオン）ペアを脱溶媒和しそしてそれを電極の孔に吸着するためのトータルエネルギー（ΔEtot）を決定することを有する。このトータルエネルギーは、静電容量に反比例するようになる。
ΔEtotの計算を容易にするために、ΔEtotを、その合計がΔEtotとなる３つのエネルギー項に分割できる熱力学サイクルを使うことがよく知られている（図１）。
１．脱溶媒和エネルギー（ΔEdesolv）
それは、イオン性溶液（カチオン＋アニオン）ペアを凝縮相（純粋な、または、アセトニトリンのような溶媒内に溶けたイオン性溶液）から気相に移すために必要なエネルギーに対応する。この技術は、Koddermann, D. Paschek, R. Ludwig, ChemPhysChem 2008, 9, 549.という文献に明白に記載されている。
２．分離エネルギー（ΔEdiss）
それは、カチオンとアニオンを分離するために必要なエネルギーを表す。この技術は、Shimizu, K.; Tariq, M.; Costa Gomes, M.F.; Rebelo, L.P.N.; Canongia Lopes, J.N. The Journal of Physical Chemistry, B. 2010, 114, 5831-4という文献に明白に記載されている。
３．吸着エネルギー（ΔEads）
それは、イオンが電極の孔の中に入り込むときに放出される。電気二重層キャパシタの構成内において、アニオンの吸着エネルギー（ΔEads-an）とカチオンの吸着エネルギー（ΔEads-cat）が計算される。
１．脱溶媒和エネルギー（ΔEdesolv）
それは、イオン性溶液（カチオン＋アニオン）ペアを凝縮相（純粋な、または、アセトニトリンのような溶媒内に溶けたイオン性溶液）から気相に移すために必要なエネルギーに対応する。この技術は、Koddermann, D. Paschek, R. Ludwig, ChemPhysChem 2008, 9, 549.という文献に明白に記載されている。
２．分離エネルギー（ΔEdiss）
それは、カチオンとアニオンを分離するために必要なエネルギーを表す。この技術は、Shimizu, K.; Tariq, M.; Costa Gomes, M.F.; Rebelo, L.P.N.; Canongia Lopes, J.N. The Journal of Physical Chemistry, B. 2010, 114, 5831-4という文献に明白に記載されている。
３．吸着エネルギー（ΔEads）
それは、イオンが電極の孔の中に入り込むときに放出される。電気二重層キャパシタの構成内において、アニオンの吸着エネルギー（ΔEads-an）とカチオンの吸着エネルギー（ΔEads-cat）が計算される。
図１は、カチオンとアニオンを凝縮相から取り出しそれらを電極の孔の中に吸着するために必要なトータルエネルギーを計算することを可能にする方針を表す。
しかしながら、この熱力学サイクルを使った吸着エネルギーの決定は、イオン（ゲストとして称する）と電極（ホスト）との間の短い作用距離および長い作用距離（short- and long-range）の静電相互作用を考慮する必要がある。ところで、それらは必須であるけれども、計算形式によって課せられた制約は、これらの相互作用が考慮されることを許さない。
したがって、本発明の目的は、イオンと電極との間の静電相互作用を考慮しつつ、電気二重層キャパシタ内のアニオン−カチオンペアの吸着エネルギーを決定する方法に関する。
本発明は、アニオン−カチオンペアをスクリーニングする方法、および、電気二重層キャパシタの電極の孔の大きさを決める方法にも関する。
概括的に言えば、本発明は、電気によって帯電された吸着体と電気によって帯電された吸着物との間の吸着エネルギーを、吸着物と吸着体との間の静電相互作用を考慮して決定する方法に関する。本方法は、以下のステージを有する：
前記吸着体（A、B）と前記吸着物（i+、i-）を、同じタイプで反対の電荷の他の吸着体（A’、B’）および同じタイプで反対の電荷の他の吸着物（i’+、i’-）と同様に含むシミュレーションボックス（BSA）を、前記シミュレーションボックスがゼロの電荷を有するように構成することと、
前記第１のシミュレーションボックス（BSA）内の前記吸着物（i+、i-およびi’）の吸着エネルギーを、分子シミュレーションによって、エワルド法を用いて決定し、そこから前記吸着体（A、B）上の前記吸着物（i+、i-）の吸着エネルギーを推測すること。
前記吸着体（A、B）と前記吸着物（i+、i-）を、同じタイプで反対の電荷の他の吸着体（A’、B’）および同じタイプで反対の電荷の他の吸着物（i’+、i’-）と同様に含むシミュレーションボックス（BSA）を、前記シミュレーションボックスがゼロの電荷を有するように構成することと、
前記第１のシミュレーションボックス（BSA）内の前記吸着物（i+、i-およびi’）の吸着エネルギーを、分子シミュレーションによって、エワルド法を用いて決定し、そこから前記吸着体（A、B）上の前記吸着物（i+、i-）の吸着エネルギーを推測すること。
本発明によれば、吸着体は、ゼオライトまたはナノチューブまたは酵素または電極（A、B）とすることができ、吸着物は、イオン（i+、i-）またはタンパク質とすることができる。
１つの実施形態によれば、イオン性溶液のアニオン−カチオンペアの吸着エネルギーが、イオンと電極との間の静電相互作用を考慮することによって、電気二重層キャパシタの２つの電極に基づき、以下のステージを実行することによって決定される。
前記陽極Aと少なくとも１つのアニオンを、反対の電荷の電極A’および反対の電荷の少なくとも１つのイオンと同様に含む第１のシミュレーションボックス（BSA）を、前記シミュレーションボックスがゼロの電荷を持つように構成することと、
前記陰極Bと少なくとも１つのカチオンを、反対の電荷の電極B’と反対の電荷の少なくとも１つのイオンと同様に含む第２のシミュレーションボックス（BSB）を、前記シミュレーションボックスがゼロの電荷を持つように構成することと、
エワルド法を使って前記第１のシミュレーションボックス（BSA）内の前記イオンの吸着エネルギーの静電的寄与を決定し、そこからアニオンの吸着エネルギーを推測することと、
エワルド法を使って前記第２のシミュレーションボックス（BSB）内の前記イオンの吸着エネルギーの静電的寄与を決定し、そこからカチオンの吸着エネルギーを推測する。
前記陽極Aと少なくとも１つのアニオンを、反対の電荷の電極A’および反対の電荷の少なくとも１つのイオンと同様に含む第１のシミュレーションボックス（BSA）を、前記シミュレーションボックスがゼロの電荷を持つように構成することと、
前記陰極Bと少なくとも１つのカチオンを、反対の電荷の電極B’と反対の電荷の少なくとも１つのイオンと同様に含む第２のシミュレーションボックス（BSB）を、前記シミュレーションボックスがゼロの電荷を持つように構成することと、
エワルド法を使って前記第１のシミュレーションボックス（BSA）内の前記イオンの吸着エネルギーの静電的寄与を決定し、そこからアニオンの吸着エネルギーを推測することと、
エワルド法を使って前記第２のシミュレーションボックス（BSB）内の前記イオンの吸着エネルギーの静電的寄与を決定し、そこからカチオンの吸着エネルギーを推測する。
我々は、溶媒から前記アニオン−カチオンペアを脱溶媒和し前記ペアを電気二重層キャパシタの２つの電極内に挿入するためのトータルエネルギーΔEtotを、以下のステージを実行することによって決定できる。
前記アニオン−カチオンペアの脱溶媒和エネルギー（ΔEdesolv）を決定することと、
前記アニオン−カチオンペアの分離エネルギー（ΔEdiss）を決定することと、
本発明に従った方法で前記アニオン−カチオンペアの吸着エネルギー（ΔEads）を決定することと、
該脱溶媒和エネルギーと該分離エネルギーと該吸着エネルギーを合計することによってトータルエネルギー（ΔEtot）における変化を決定すること。
前記アニオン−カチオンペアの脱溶媒和エネルギー（ΔEdesolv）を決定することと、
前記アニオン−カチオンペアの分離エネルギー（ΔEdiss）を決定することと、
本発明に従った方法で前記アニオン−カチオンペアの吸着エネルギー（ΔEads）を決定することと、
該脱溶媒和エネルギーと該分離エネルギーと該吸着エネルギーを合計することによってトータルエネルギー（ΔEtot）における変化を決定すること。
本発明によれば、脱溶媒和エネルギー（ΔEdesolv）は、第１の分子力学シミュレーションを実行して所定温度での凝縮相の平均トータルエネルギーを計算し、第２の分子力学シミュレーションを実行して１つのイオンペアの平均トータルエネルギーを計算することによって決定可能であり、そして、分離エネルギー（ΔEdiss）は、アニオン−カチオンペアのエネルギー、カチオンのエネルギーおよびアニオンのエネルギーを決定することで決定可能である。
１つの実施形態によれば、電気二重層キャパシタの電極を形成する材料のスクリーニングは、以下のステージ：
イオン性溶液のためにアニオン−カチオンペアを選択し、
本発明に従った方法で電極の異なる孔サイズごとに、これらのイオンのトータルエネルギーΔEtotを決定し、
最小のトータルエネルギーΔEtotに対応する孔サイズを選択することによって最大の静電容量を得ることを可能にする孔サイズを決定すること、
を実行することによって実行される。
イオン性溶液のためにアニオン−カチオンペアを選択し、
本発明に従った方法で電極の異なる孔サイズごとに、これらのイオンのトータルエネルギーΔEtotを決定し、
最小のトータルエネルギーΔEtotに対応する孔サイズを選択することによって最大の静電容量を得ることを可能にする孔サイズを決定すること、
を実行することによって実行される。
電気二重層キャパシタのイオン性溶液のカチオン−アニオンペアのスクリーニングは、以下のステージ：
電極のために孔サイズを選択し、
本発明に従った方法で前記孔サイズについて異なるアニオン−カチオンペアのトータルエネルギーΔEtotを決定し、
前記孔サイズで最小のトータルエネルギー（ΔEtot）を持つペアを選択することによって最大の静電容量を得ることを可能にするアニオン−カチオンペアを選択すること、
を実行することによって実行可能である。
本発明に従った方法の他の特徴および効果は、添付図面に関して、非限定的な例の目的で与えられた実施形態の以下の記載を読むことから明らかであろう。
電極のために孔サイズを選択し、
本発明に従った方法で前記孔サイズについて異なるアニオン−カチオンペアのトータルエネルギーΔEtotを決定し、
前記孔サイズで最小のトータルエネルギー（ΔEtot）を持つペアを選択することによって最大の静電容量を得ることを可能にするアニオン−カチオンペアを選択すること、
を実行することによって実行可能である。
本発明に従った方法の他の特徴および効果は、添付図面に関して、非限定的な例の目的で与えられた実施形態の以下の記載を読むことから明らかであろう。
本発明は、電気によって帯電された吸着体と電気によって帯電された吸着物との間の静電（クーロン力）相互作用をシミュレーションする方法である。
その例は、
ゼオライト内のイオン、
MOF（金属有機構造体）内のイオン、
ナノチューブ内のイオン、
酵素内のタンパク質である。
その例は、
ゼオライト内のイオン、
MOF（金属有機構造体）内のイオン、
ナノチューブ内のイオン、
酵素内のタンパク質である。
特に、本発明は、イオン性溶液のアニオン−カチオンペアを脱溶媒和し、このペアを電気二重層キャパシタの電極の中に挿入するためのトータルエネルギーを決定する方法である。
図１に示したように、トータルエネルギーΔEtotを、合計がΔEtotとなる３つのエネルギー項に分割することを可能にする熱力学サイクルが使用される。
それゆえ、本方法は、以下のステージを有する：
１．脱溶媒和エネルギー（ΔEdesolv）を決定すること
２．分離エネルギー（ΔEdiss）を決定すること
３．吸着エネルギー（ΔEads）を決定すること
１．脱溶媒和エネルギー（ΔEdesolv）を決定すること
２．分離エネルギー（ΔEdiss）を決定すること
３．吸着エネルギー（ΔEads）を決定すること
１．脱溶媒和エネルギー（ΔE desolv ）を決定すること
脱溶媒和エネルギーは、イオン性溶液（カチオン＋アニオン）ペアを、凝縮相（純粋なまたはアセトニトリルのような溶媒の中に溶けているイオン性溶液）から気相に移すために必要なエネルギーに対応する。
脱溶媒和エネルギーは、イオン性溶液（カチオン＋アニオン）ペアを、凝縮相（純粋なまたはアセトニトリルのような溶媒の中に溶けているイオン性溶液）から気相に移すために必要なエネルギーに対応する。
それを計算するために、２つのシミュレーションが必要となる。
１）分子動力学シミュレーションが実行されて、所定の温度での凝縮相（ｎ個のイオンペアを有する）の平均トータルエネルギーE1を計算する。このような技術は、例えば以下の文献に記載されている。
Allen, M. P.; Tildesley, D. J. Computer simulation of liquids; Oxford: Clarendon press, Ed.; 1987th ed.; Oxford University Press: Oxford, 1987.
Frenkel, D.; Smit, B. Understanding Molecular Simulation; 2nd ed.; Academic Press: London, 2002; p. 638.
２）分子動力学シミュレーションが実行されて、１つのイオンペア（カチオン＋アニオン）の平均トータルエネルギーE2を計算する。
１）分子動力学シミュレーションが実行されて、所定の温度での凝縮相（ｎ個のイオンペアを有する）の平均トータルエネルギーE1を計算する。このような技術は、例えば以下の文献に記載されている。
Allen, M. P.; Tildesley, D. J. Computer simulation of liquids; Oxford: Clarendon press, Ed.; 1987th ed.; Oxford University Press: Oxford, 1987.
Frenkel, D.; Smit, B. Understanding Molecular Simulation; 2nd ed.; Academic Press: London, 2002; p. 638.
２）分子動力学シミュレーションが実行されて、１つのイオンペア（カチオン＋アニオン）の平均トータルエネルギーE2を計算する。
このように、脱溶媒和エネルギーは、ペアの数（n）によって割られたE1とE2との間の差、すなわち、ΔE desolv ＝E1/n - E2である。
このような技術は、例えば、以下の文献に記載されている。
Shimizu, K.; Tariq, M.; Costa Gomes, M. F.; Rebelo, L. P. N.; Canongia Lopes, J. N. The journal of physical chemistry. B 2010, 114, 5831-4.
Shimizu, K.; Tariq, M.; Costa Gomes, M. F.; Rebelo, L. P. N.; Canongia Lopes, J. N. The journal of physical chemistry. B 2010, 114, 5831-4.
２．分離エネルギー（ΔE diss ）を決定すること
分離エネルギーは、カチオンとアニオンを分離するために必要なエネルギーを表す。それは、３つの量子計算から計算される。
１）カチオン＋アニオンペアのエネルギー：E3
２）カチオンのエネルギー：E4
３）アニオンのエネルギー：E5
分離エネルギーは、カチオンとアニオンを分離するために必要なエネルギーを表す。それは、３つの量子計算から計算される。
１）カチオン＋アニオンペアのエネルギー：E3
２）カチオンのエネルギー：E4
３）アニオンのエネルギー：E5
我々は、以降、ΔE diss ＝E3−E4−E5と記載する。
このような技術は、例えば、以下の文献に記載されている。
Fernandes, A. M.; Rocha, M. A. A.; Freire, M. G.; Marrucho, I. M.; Coutinho, J. A. P.; Santos, L. M. N. B. F. The journal of physical chemistry. B 2011, 115, 4033-41.
Fernandes, A. M.; Rocha, M. A. A.; Freire, M. G.; Marrucho, I. M.; Coutinho, J. A. P.; Santos, L. M. N. B. F. The journal of physical chemistry. B 2011, 115, 4033-41.
３．吸着エネルギー（ΔE
ads
）を決定すること
吸着エネルギーは、イオンが電極（本発明によればカーボンナノチューブによって表される）の孔の中に挿入されたときリリースされるエネルギーである。
吸着エネルギーは、分子シミュレーションにて計算され、そこでは、ショートレンジおよびロングレンジの静電相互作用がエワルド法を使って以下のように計算される
ここで、Eelectrode+ion、EelectrodeおよびEionは、電極とイオンのエネルギー、電極のみのエネルギーおよびイオンのみのエネルギーにそれぞれ対応する。このエネルギーは、力場（force field）で計算されたエネルギーであり、力場のパラメータはイオン性溶液を十分に表すために最適化されている。この力場は、以下の文献に基づくことができる。
Canongia Lopes, J. N.; Deschamps, J.; Padua, A. A. H. The Journal of Physical Chemistry B 2004, 108, 11250.
Canongia Lopes, J. N.; Padua, A. A. H. The Journal of Physical Chemistry B 2006, 110, 19586-19592.
Canongia Lopes, J. N.; Padua, A. A. H. The Journal of Physical Chemistry B 2004, 108, 16893-16898.
De Andrade, J.; Bo, E. S.; Stassen, H. Journal of Physical Chemistry B 2002, 3546-3548.
Kaminski, G. A.; Jorgensen, W. L. Journal of the Chemical Society, Perkin Transactions 2 1999, 2365-2375.
Canongia Lopes, J. N.; Deschamps, J.; Padua, A. A. H. The Journal of Physical Chemistry B 2004, 108, 11250.
Canongia Lopes, J. N.; Padua, A. A. H. The Journal of Physical Chemistry B 2006, 110, 19586-19592.
Canongia Lopes, J. N.; Padua, A. A. H. The Journal of Physical Chemistry B 2004, 108, 16893-16898.
De Andrade, J.; Bo, E. S.; Stassen, H. Journal of Physical Chemistry B 2002, 3546-3548.
Kaminski, G. A.; Jorgensen, W. L. Journal of the Chemical Society, Perkin Transactions 2 1999, 2365-2375.
各ポテンシャルエネルギー（Epot）は、２つの重要な項：分子内エネルギー（Eintra）と分子間エネルギー（Einter）に分けることができる。例えば、
第１の項（Eintra）は、原子の手（bond）、結合角または二面角にてリンクされた原子間の相互作用を考慮している。第２の項（Einter）は、いわゆる非結合相互作用である、ファンデルワークス力と静電気力を含む。
一般論として、種（イオン、電極、イオンペア）のポテンシャルエネルギーは、
によって定義される。
ここで、
li＝原子の手（bond）iの長さ
li,0＝この原子の手のための基準距離
θi＝結合角iの大きさ
θi,0＝この結合角のための基準
ω＝二面角
Vn＝その二面角のための定数
γ＝その二面角のための定数
εij＝原子iおよびjのペアのための基準エネルギー
σij＝原子iおよびjのペアのための基準距離
ここで、
li＝原子の手（bond）iの長さ
li,0＝この原子の手のための基準距離
θi＝結合角iの大きさ
θi,0＝この結合角のための基準
ω＝二面角
Vn＝その二面角のための定数
γ＝その二面角のための定数
εij＝原子iおよびjのペアのための基準エネルギー
σij＝原子iおよびjのペアのための基準距離
この定義は、Leach, A. R. Molecular Modelling: Principles and Applications; 2nd ed.; Prentice Hall, 2001.に、明白に記載されている。
静電相互作用は、（非常に）長い作用距離（long range）を有する。このため、長い距離で分離された複数の帯電原子（または複数の帯電体）でさえ、静電相互作用を受ける。この相互作用を考慮するために、エワルド法が非常に一般的に適用される。エワルド法は、分子シミュレーションのために広く使われている方法である。それは、実空間における短い距離と逆空間における長い距離での静電相互作用を評価することを可能にし、その合計が、ある値に収束する。
エワルド法は、周期系の相互作用エネルギー、より具体的には静電エネルギーを計算するための方法である。静電エネルギーは、短い作用距離および長い作用距離の相互作用項の両方を有しており、相互作用ポテンシャルを、短い作用距離項（その合計が実空間で実行される）と長い作用距離項（その合計がフーリエ空間（逆空間）で実行される）に分離することは非常に興味深い。このアプローチの長所は、長い作用距離の相互作用の場合において、フーリエ空間での合計の収束が、実空間でのそれと比較して速いということである。
カチオンとアニオンの全ての原子は、部分的な実効電荷を持っている。このため、各原子は、「点電荷」として考えられる。そして、その系での各「点電荷」は、以下の式で計算されるトータル静電エネルギーに寄与する。
ここで、
第１の項は、実空間での寄与を表し、
第２の項は、逆空間での寄与を表し、
第３の項は、それ自身の各ガウス関数の間の相互作用の補正を表し、
ここで（主要変数の定義）、
Eelec：静電エネルギー（２つの電気によって帯電されたパーティクル間の相互作用）、
qiまたはqj：原子iおよびjの実効原子電荷
rij：電荷qiおよびqjをそれぞれ運ぶ原子iおよびjの間の距離
N：電荷を運ぶ原子のトータル数（電極の原子＋イオン性溶液の原子）
n = (nxL, nyL, nzL) nx, ny, nzは整数
L：ボックスの長さ
erfc：相補誤差関数：
a：ガウス関数
の「幅」を定義
k：2πn/L（逆格子ベクトル）
k：ベクトルkのノルム
第１の項は、実空間での寄与を表し、
第２の項は、逆空間での寄与を表し、
第３の項は、それ自身の各ガウス関数の間の相互作用の補正を表し、
ここで（主要変数の定義）、
Eelec：静電エネルギー（２つの電気によって帯電されたパーティクル間の相互作用）、
qiまたはqj：原子iおよびjの実効原子電荷
rij：電荷qiおよびqjをそれぞれ運ぶ原子iおよびjの間の距離
N：電荷を運ぶ原子のトータル数（電極の原子＋イオン性溶液の原子）
n = (nxL, nyL, nzL) nx, ny, nzは整数
L：ボックスの長さ
erfc：相補誤差関数：
k：2πn/L（逆格子ベクトル）
k：ベクトルkのノルム
エワルド法は、静電部分を計算する。しかしながら、この方法は、シミュレーションボックス内のトータルの電荷がゼロ（中性）にする必要があり、他の状態では、式（３）内の第２項が収束しない。シミュレーションボックスは、事象が調べられる物理環境の再現（representation）、モデルである。
ここで、イオンが電極に挿入されたときの物理系のダイレクトシミュレーションは、電極のみが帯電されているか、または、電極とイオンの系が帯電されているかを含意する。その結果、エワルド法は、もはや適用できない。
この問題を克服するために、本発明による方法は、図２に示したように、各電極が、正確に反対の電荷を受け取ることによって二重化されている。
このため、本発明によれば、２つのシミュレーションボックスBSAおよびBSBが物理系（SP）で構成されている。
物理系内の陽極Aとアニオン（i-）は、電極Aとアニオン（i-）を、反対電荷の電極A’と反対電荷のイオン（i’+）と共に有するシミュレーションボックスBSA内で調べられる。
物理系内の陰極Bとカチオン（i+）は、電極Bとカチオン（i+）を、反対電荷の電極B’と反対電荷のイオン（i’-）と共に有するシミュレーションボックスBSB内で調べられる。
各シミュレーションボックスにおいて電気的中性を確実にするために、カチオンの数と反対電荷のカチオン（カチオン’）の数は同じであり、トータルの電荷がゼロになる。同様に、電極AとA’の電荷が互いにキャンセルされる。
このため、エワルド法を適用できる。
電極A内のアニオンの吸着エネルギーを得るために、我々は、シミュレーションボックスBSAで計算された吸着エネルギー
を２つに分ける。
電極B内のカチオンの吸着エネルギーを得るために、我々は、シミュレーションボックスBSBで計算された相互作用エネルギーを２つに分ける。
電極B内のカチオンの吸着エネルギーを得るために、我々は、シミュレーションボックスBSBで計算された相互作用エネルギーを２つに分ける。
このように、より詳細には、各電極は、トータルの電荷がゼロになるように、距離ｄだけ離れ反対であるが同じ絶対値の電荷を正確に有する無限長の複数のナノチューブにてモデル化される。距離ｄは、静電力およびファンデルワークス力相互作用が無視できるように、２つのナノチューブが互いから十分に離れるように選択されている。各ナノチューブは、同じ性質であるが反対の原子電荷を有する同じ数（1、2、3…）のイオンで満たされている。１以上のアニオンは、陽電荷を有するナノチューブの中に挿入可能であり、第２のナノチューブには、カチオンになるように、同じアニオンであるが反対の原子電荷を有する負電荷が挿入可能である。
この説明は、図３に示されており、上のナノチューブは2-の電荷を保持し、下のナノチューブは、同じ絶対値で反対の電荷（2+）を保持する。下のナノチューブは、３つの（ヘキサフルオロホスフェイト）[PF6]-アニオンでチャージされ、このため、このサブシステムのトータルの電荷は−１となり、一方、上のナノチューブは、反対電荷を持つ３つのアニオン[PF6]+でチャージされ、このサブシステムのトータルの電荷は＋１となり、シミュレーションボックスのトータル電荷はゼロである。
我々は、原子電荷のみを対比し、力場の他の全てのパラメータ（パーティクルの系のポテンシャルエネルギーを表す数式とパラメータのセット）は依然として同一であり、トータルの内部エネルギーが２つの種（アニオンとカチオン）で正確に同一となる、ということを、特に言及する。このように、複数のナノチューブは、２つのイオンのように反対の電荷を持つので、電気的中性はトータルの系に常に提供される。このため、エワルド法は、静電相互作用を評価するために使用可能となる。
このアプローチを使用することは、いくつかの利点を含む。
１．電気的中性がトータルの系に常に提供される。その理由は、複数のナノチューブは２つのイオンのように反対の電荷をもつからである。エワルド法が、静電相互作用の評価に使用可能になる。
２．ある数のイオンで満たされたナノチューブのトータルの電荷は、ゼロである必要がない。その理由は、この実効電荷は、同じ数のイオン（しかし反対の電荷）を持つ他方のナノチューブの電荷で補償されるからである。このアプローチは、以下のことを可能にする。
ａ）ナノチューブの電荷に関わらず、各ナノチューブ内のイオンの数を変更すること。このため、ナノチューブの一定の電荷で、ナノチューブに入れることができるイオン（アニオンとカチオン）の最大数が調査可能になる。
ｂ）挿入されたイオンの数に関わらず、ナノチューブの電荷を変更すること。もし、ナノチューブの直径がイオンのサイズに対して「非常に小さく」てナノチューブの電荷を増加することによってフィットさせることができないと、静電相互作用は、イオンの変形エネルギーを補償できる。
１．電気的中性がトータルの系に常に提供される。その理由は、複数のナノチューブは２つのイオンのように反対の電荷をもつからである。エワルド法が、静電相互作用の評価に使用可能になる。
２．ある数のイオンで満たされたナノチューブのトータルの電荷は、ゼロである必要がない。その理由は、この実効電荷は、同じ数のイオン（しかし反対の電荷）を持つ他方のナノチューブの電荷で補償されるからである。このアプローチは、以下のことを可能にする。
ａ）ナノチューブの電荷に関わらず、各ナノチューブ内のイオンの数を変更すること。このため、ナノチューブの一定の電荷で、ナノチューブに入れることができるイオン（アニオンとカチオン）の最大数が調査可能になる。
ｂ）挿入されたイオンの数に関わらず、ナノチューブの電荷を変更すること。もし、ナノチューブの直径がイオンのサイズに対して「非常に小さく」てナノチューブの電荷を増加することによってフィットさせることができないと、静電相互作用は、イオンの変形エネルギーを補償できる。
３．システムのシミュレーションは、イオンとナノチューブの間の平均相互作用エネルギー（ナノチューブごとの平均エネルギー）の２倍に対応する計算をもたらすので、平均吸着エネルギーの計算の統計も改良される。
例
吸着エネルギー
図４において、吸着エネルギー（ΔEads）は、ファンデルワールス力相互作用（Evdw）と静電相互作用（Eelectro）とからの寄与に分解されている。ファンデルワールス力相互作用（Evdw）は円で、静電相互作用（Eelectro）は三角で、吸着エネルギー（(ΔEads）は四角で表されている。
吸着エネルギー
図４において、吸着エネルギー（ΔEads）は、ファンデルワールス力相互作用（Evdw）と静電相互作用（Eelectro）とからの寄与に分解されている。ファンデルワールス力相互作用（Evdw）は円で、静電相互作用（Eelectro）は三角で、吸着エネルギー（(ΔEads）は四角で表されている。
これらのエネルギーは、ナノチューブの内径（DIN）の関数としてプロットされ、この直径は、テトラエチルアンモニウム（TEA+）とエチルメチルイミダゾリウム（EMIM+）のカチオンと、テトラフルオロホウ塩酸（BF4 -）とビス（トリフルオロメタンスルホニル）イミド（TFSI-）のアニオンとのそれぞれについて、電極の孔サイズをモデル化している。
結果は、一定の電荷（ナノチューブごとに±2e）では、静電の寄与は、電極の孔サイズおよびイオンの種類と関わりないことを示し、実際、イオンとナノチューブ（常に±2e）のトータル電荷（常に±1e）によってもっぱら決定される。
一方、吸着エネルギーは、最小値を有する曲線を示し、その形状はファンデルワールスエネルギーの寄与を表す曲線の形状によって影響を受けている。
図４に示した吸着エネルギー曲線の最小値は、イオンのサイズに依存している。幾何学サイズが小さいイオン（BF4およびEMIM）ほど、小さい径（5〜6オングストローム）を有するナノチューブに中に入ることができる。すなわち、吸着エネルギーは、DINが4.5オングストロームよりも大きい限りはマイナスであり、イオンが大きくなるほど中に入れなくなる。すなわち、DINが5.5オングストローム（TFSI）または6オングストロームあたりより小さくなるとすぐ、吸着エネルギーはプラスになる（図４）。
イオンとホストの間の相互作用エネルギー（井戸の深さ）は、大きい数の原子からなるイオン（TFSI：15個の原子、EMIM：19個の原子、NEt4：29個の原子）では、小さい数の原子を有するイオン（BF4：5個の原子）よりも大きくなる。この相互作用は、各原子がファンデルワールス相互作用の総和に寄与するため高い。イオンの形状も重要な役割を果たすことを言及する。TFSIのような伸長した形状を有するイオンは、NEt4のようなより球形状を有するイオンよりも、円柱（ナノチューブ）に対してより高い相互作用を有する。
方法の適用の例
静電エネルギーを考慮してイオン性溶液のアニオン＋カチオンペアを脱溶媒和しそれを電極の孔に吸着するためのトータルエネルギー（ΔEtot）を決定するために本発明に応じた方法を使うことは、電極の孔サイズ、または、イオン性溶液内のアニオン−カチオンペアの関数として、電気二重層キャパシタの静電容量の予測を可能にする。
そして、所定のイオン性溶液について前記電気二重層キャパシタの最大静電容量を得ることを許容する電極孔サイズを決定することを可能にする。
所定の電極孔サイズについて前記電気二重層キャパシタの最大静電容量を得ることを許容するアニオン−カチオンペアを決定することも可能にする。
孔サイズの関数として電気二重層キャパシタの静電容量を予測すること
（測定された）静電容量曲線と、（孔サイズ分布を考慮して）本発明に従って決定されたトータルエネルギー（ΔEtot）を表す曲線は、実質的に同じ孔サイズで、最大と最小をそれぞれ持つ。
このように、トータルエネルギーを表す曲線は、電極の最適の孔サイズを評価するための非常に良い指標である
そして、電気二重層キャパシタの電極を形成する材料をスクリーニングする方法は、以下のステージを有する：
イオン性溶液のためにアニオン−カチオンペア（電気二重層キャパシタの電解液）を選択すること、
本発明に従った方法で、異なる電極孔サイズについてこれらのイオンのトータルエネルギー（ΔEtot）を決定すること、
最小エネルギー（ΔEtot）に応じた孔サイズを選択することによって、最大の静電容量を得ることを可能にする孔サイズを決定すること。
イオン性溶液のためにアニオン−カチオンペア（電気二重層キャパシタの電解液）を選択すること、
本発明に従った方法で、異なる電極孔サイズについてこれらのイオンのトータルエネルギー（ΔEtot）を決定すること、
最小エネルギー（ΔEtot）に応じた孔サイズを選択することによって、最大の静電容量を得ることを可能にする孔サイズを決定すること。
例
以下のアニオン−カチオンペアが選択される：EMIM/TFSI。
本発明に従った方法を適用することによって、これらのイオンのトータルエネルギーが、異なる電極孔サイズについて決定される。
トータルエネルギー曲線は、P. Simon et al. (Science (New York, N.Y.) 2006, 313, 1760-3)にて実験的に測定された正規化静電容量（normalized capacitance）と類似している。これは、一方では、実験的に測定された電気二重層キャパシタ（EMIM/TFSIペアを有する）の正規化された静電容量（NC−ひし形を有する曲線）の展開を、孔サイズ（PS）の関数として表し、他方では、トータルエネルギー（ΔEtot−円を有する曲線）の展開を孔サイズの関数として示す図５において見ることができる。
「静電容量」曲線の最大とΔEtot曲線の最小は実質的に0.7nmの孔サイズにところにあることが観測される。
この孔サイズは、EMIM/TFSIペアで動作する電気二重層キャパシタにとって最適の孔サイズである。
電気二重層キャパシタの静電容量をアニオン−カチオンペアの関数として予測すること
反対に、所定の電極の孔サイズについて、前記電気二重層キャパシタの最大静電容量を得ることを可能にするアニオン−カチオンペアを決定する方法は、以下のステージを有する。
電極のための孔サイズを選択すること、
本発明に従った方法で前記孔サイズについて異なるアニオン−カチオンペアのトータルエネルギー（ΔEtot）を決定し、
前記孔サイズについて最小のトータルエネルギー（ΔEtot）を有するペアを選択することで、最大の静電容量を得ることが可能なアニオン−カチオンペアを選択すること。
電極のための孔サイズを選択すること、
本発明に従った方法で前記孔サイズについて異なるアニオン−カチオンペアのトータルエネルギー（ΔEtot）を決定し、
前記孔サイズについて最小のトータルエネルギー（ΔEtot）を有するペアを選択することで、最大の静電容量を得ることが可能なアニオン−カチオンペアを選択すること。
Claims (5)
- 電気によって帯電された吸着体と電気によって帯電された吸着物との間の吸着エネルギーを、前記吸着体と前記吸着物との間の静電相互作用を考慮するとともに、
前記吸着体（A、B）と前記吸着物（i+、i-）とを含み、同じタイプで反対の電荷の他の吸着体（A’、B’）と同じタイプで反対の電荷の他の吸着物（i’+、i’-）とをさらに含むシミュレーションボックス（BSA）を、前記シミュレーションボックスがゼロの電荷を有するように構成する工程と、
前記シミュレーションボックス（BSA）内の前記吸着物（i+、i-およびi’）の吸着エネルギーを、分子シミュレーションによって、エワルド法を使って決定し、そこから前記吸着体（A、B）上の前記吸着物（i+、i-）の吸着エネルギーを推測する工程と、を実行することで決定する方法であって、
前記吸着体は、ゼオライトまたはナノチューブまたは酵素または電極（A、B）であり、前記吸着物は、イオン（i+、i-）またはタンパク質であり、
電気二重層キャパシタの２つの電極について、イオン性溶液のアニオン−カチオンペアの吸着エネルギーが、イオンと電極の間の静電相互作用を考慮するとともに、
前記陽極Aと少なくとも１つのアニオンとを含み、反対の電荷の電極A’と反対の電荷の少なくとも１つのイオンとをさらに含む第１のシミュレーションボックス（BSA）を、前記シミュレーションボックスがゼロの電位を有するように構成する工程と、
前記陰極Bと少なくとも１つのカチオンとを含み、反対の電荷の電極B’と反対の電荷の少なくとも１つのイオンとをさらに含む第２のシミュレーションボックス（BSB）を、前記シミュレーションボックスがゼロの電位を有するように構成する工程と、
エワルド法を使って前記第１のシミュレーションボックス（BSA）内の前記イオンの吸着エネルギーの静電寄与を決定し、そして、そこからアニオンの吸着エネルギーを推定する工程と、
エワルド法を使って前記第２のシミュレーションボックス（BSB）内の前記イオンの吸着エネルギーの静電寄与を決定し、そして、そこからカチオンの吸着エネルギーを推定する工程と、を実行することで決定される、方法。 - 溶媒から前記アニオン−カチオンペアを脱溶媒和し前記ペアを電気二重層キャパシタの２つの電極の中に挿入するためのトータルエネルギーΔEtotが、
前記アニオン−カチオンペアの脱溶媒和エネルギー（ΔEdesolv）を決定する工程と、
前記アニオン−カチオンペアの分離エネルギー（ΔEdiss）を決定する工程と、
請求項１に記載の方法で前記アニオン−カチオンペアの吸着エネルギー（ΔEads）を決定する工程と、
前記脱溶媒和エネルギー、前記分離エネルギーおよび前記吸着エネルギーを合計することによってトータルエネルギー（ΔEtot）の変化を決定するする工程と、を実行することで決定される、請求項１に記載の方法。 - 前記脱溶媒和エネルギー（ΔEdesolv）が、所定温度での凝縮相の平均トータルエネルギーを計算する第１の分子力学シミュレーションと、１つのイオンペアでの平均トータルエネルギーを計算する第２の分子力学シミュレーションと、を実行することによって決定され、
前記分離エネルギー（ΔEdiss）が、アニオン−カチオンペアのエネルギー、前記カチオンのエネルギーおよび前記アニオンのエネルギーを決定することによって決定される、請求項２に記載の方法。 - 電気二重層キャパシタの電極を構成する材料のスクリーニングが、
前記イオン性溶液のためのアニオン−カチオンペアを選択する工程と、
請求項２または３に記載の方法で、前記電極の異なる孔サイズについて、これらのイオンについての前記トータルエネルギーΔEtotを決定する工程と、
最小のトータルエネルギーΔEtotに対応する孔サイズを選択することによって、最大の静電容量を得ることを可能にする孔サイズを決定するする工程と、を実行することで実行される、請求項２または３に記載の方法。 - 電気二重層キャパシタのイオン性溶液のカチオン−アニオンペアのスクリーニングが、
前記電極の孔サイズを選択する工程と、
請求項２または３に記載の方法で、前記孔サイズについて、異なるアニオン−カチオンペアのトータルエネルギーΔEtotを決定する工程と、
前記孔サイズでの最小のトータルエネルギー（ΔEtot）を有するペアを選択することによって、最大の静電容量を得ることを可能にするアニオン−カチオンペアを選択するする工程と、を実行することで実行される、請求項２または３に記載の方法。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR1202895 | 2012-10-29 | ||
| FR1202895A FR2997316B1 (fr) | 2012-10-29 | 2012-10-29 | Procede pour determiner l'energie totale pour desolvater un couple anion-cation et l'inserer dans l'electrode d'un supercondensateur |
| PCT/FR2013/052137 WO2014068204A1 (fr) | 2012-10-29 | 2013-09-17 | Procédé pour déterminer l'énergie totale pour désolvater un couple anion-cation et l'insérer dans l'électrode d'un supercondensateur |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2015537380A JP2015537380A (ja) | 2015-12-24 |
| JP6341922B2 true JP6341922B2 (ja) | 2018-06-13 |
Family
ID=47750750
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2015540192A Expired - Fee Related JP6341922B2 (ja) | 2012-10-29 | 2013-09-17 | アニオン−カチオンペアを脱溶媒和しそれを電気二重層キャパシタの電極に挿入するためのトータルエネルギーを決定する方法 |
Country Status (8)
| Country | Link |
|---|---|
| US (1) | US20150347718A1 (ja) |
| EP (1) | EP2912674B1 (ja) |
| JP (1) | JP6341922B2 (ja) |
| KR (1) | KR20150079831A (ja) |
| CN (1) | CN104756119B (ja) |
| CA (1) | CA2888435A1 (ja) |
| FR (1) | FR2997316B1 (ja) |
| WO (1) | WO2014068204A1 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10606969B2 (en) * | 2017-10-12 | 2020-03-31 | The Boeing Company | Predicting electrostatic charges in a liquid container |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005027159A1 (ja) * | 2003-09-11 | 2005-03-24 | Matsushita Electric Industrial Co., Ltd. | 電気二重層コンデンサの製造方法 |
| US20050118091A1 (en) * | 2003-12-01 | 2005-06-02 | Cooper Alan C. | Hydrogen storage utilizing carbon nanotube materials |
| EP1746612B1 (en) * | 2004-05-13 | 2012-04-04 | Panasonic Corporation | Electric double-layer capacitor and production method therefor |
| JP5207338B2 (ja) * | 2006-06-05 | 2013-06-12 | トヨタ自動車株式会社 | 多孔質炭素材料、多孔質炭素材料の製造方法、電気二重層キャパシタ |
| CN101108194B (zh) * | 2007-06-30 | 2010-12-01 | 广西壮族自治区化工研究院 | 一种除去右旋糖酐铁络合物水溶液中氯化钠的方法及装置 |
| US20110111284A1 (en) * | 2008-06-24 | 2011-05-12 | Panasonic Corporation | Activated carbon for electrochemical element and electrochemical element using the same |
-
2012
- 2012-10-29 FR FR1202895A patent/FR2997316B1/fr not_active Expired - Fee Related
-
2013
- 2013-09-17 WO PCT/FR2013/052137 patent/WO2014068204A1/fr active Application Filing
- 2013-09-17 KR KR1020157013855A patent/KR20150079831A/ko not_active Application Discontinuation
- 2013-09-17 EP EP13779253.7A patent/EP2912674B1/fr not_active Not-in-force
- 2013-09-17 CN CN201380056488.3A patent/CN104756119B/zh not_active Expired - Fee Related
- 2013-09-17 US US14/654,635 patent/US20150347718A1/en not_active Abandoned
- 2013-09-17 CA CA 2888435 patent/CA2888435A1/fr not_active Abandoned
- 2013-09-17 JP JP2015540192A patent/JP6341922B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2015537380A (ja) | 2015-12-24 |
| KR20150079831A (ko) | 2015-07-08 |
| US20150347718A1 (en) | 2015-12-03 |
| CN104756119B (zh) | 2017-11-17 |
| CN104756119A (zh) | 2015-07-01 |
| EP2912674A1 (fr) | 2015-09-02 |
| FR2997316B1 (fr) | 2014-11-21 |
| EP2912674B1 (fr) | 2016-08-24 |
| CA2888435A1 (fr) | 2014-05-08 |
| FR2997316A1 (fr) | 2014-05-02 |
| WO2014068204A1 (fr) | 2014-05-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | MXene‐reduced graphene oxide aerogel for aqueous zinc‐ion hybrid supercapacitor with ultralong cycle life | |
| Lozovoi et al. | Ab initio simulation of charged slabs at constant chemical potential | |
| Hutzenlaub et al. | Three-dimensional electrochemical Li-ion battery modelling featuring a focused ion-beam/scanning electron microscopy based three-phase reconstruction of a LiCoO2 cathode | |
| Zhang et al. | Modelling electrochemical systems with finite field molecular dynamics | |
| Breitsprecher et al. | Electrode models for ionic liquid-based capacitors | |
| Lenze et al. | Simulation-supported analysis of calendering impacts on the performance of lithium-ion-batteries | |
| Wang et al. | Electric potential calculation in molecular simulation of electric double layer capacitors | |
| Neto et al. | Differential capacitance and energetics of the electrical double layer of graphene oxide supercapacitors: Impact of the oxidation degree | |
| Yang et al. | Reliability of Constant Charge Method for Molecular Dynamics Simulations on EDLCs in Nanometer and Sub‐Nanometer Spaces | |
| Fileti et al. | Investigating the asymmetry in the EDL response of C 60/graphene supercapacitors | |
| Shin et al. | Metal–organic framework supercapacitors: challenges and opportunities | |
| Kong et al. | Molecular dynamics for the charging behavior of nanostructured electric double layer capacitors containing room temperature ionic liquids | |
| JP6341922B2 (ja) | アニオン−カチオンペアを脱溶媒和しそれを電気二重層キャパシタの電極に挿入するためのトータルエネルギーを決定する方法 | |
| Ers et al. | Double layer in ionic liquids: Temperature effect and bilayer model | |
| Kumbhar et al. | Mesoscale analysis of Li‐ion battery microstructure using sequential coupling of discrete element and finite element method | |
| Balagopal et al. | Effect of anode conductivity degradation on the Thevenin Circuit Model of lithium ion batteries | |
| Houck et al. | Parameterization and modeling protocols for ultra-fast charging Wadsley-Roth lithium-ion batteries from coin to pouch cells | |
| Romero-Becerril et al. | Reduced order dynamical model for supercapacitors | |
| Nasimi et al. | Investigation of different effects on the capacity of supercapacitor comprising an ionic liquid between graphene oxide nanosheets | |
| Marrocchelli et al. | Molecular dynamics simulations of electrochemical energy storage devices | |
| Singh et al. | Double layer capacitance for different thicknesses of carbon overlay and types of electrolytes: a circuit analogue approach | |
| German et al. | Electrochemical double layer capacitors (supercapacitors) ageing impacts and comparison on different impedance models | |
| Pavaskar et al. | Modeling electric double-layer capacitors using charge variation methodology in gibbs ensemble | |
| Patra et al. | Electric double layer at a metal/electrolyte interface: A density functional approach | |
| He et al. | Equivalent circuit model for electrosorption with redox active materials |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160916 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20170706 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20170822 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20171121 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20180424 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20180515 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6341922 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |