Multiple ion layer formation associated with molecular orientation alteration and hydrogen bonding in a supercapacitor model with 1-ethyl-3-methylimidazolium bis-(trifluoromethylsulfonyl)imide ([C 2 mim][NtF 2 ]) as a roomtemperature ionic liquid (RTIL) electrolyte has been investigated using potential-bias classical molecular dynamics (CMD) simulations. A 50-ion pair model was used to observe the molecular ion layer formation of [C 2 mim][NtF 2 ] of the simulated electrified graphite electrode. After applying the potential bias, multiple layers of [C 2 mim] + and [NtF 2 ] − formed and became more pronounced.An orientation analysis indicated that the higher potential differences play a significant role in [C 2 mim] + , especially near the negative surface. The π−π stacking interaction between [C 2 mim] + rings and graphene electrodes was clearly perceived. Principal component analysis (PCA) was introduced to address the characteristics of the layer formation, respectively, to the potential bias applied to the simulated model. PCA scores from the applied potential (ΔE = 0.0, 1.0, 2.0, 3.0, and 4.0 V) were used to group the characteristics of the atom pairing associated with the number of H-bonds between H atoms in the [C 2 mim] + cation and all atoms in the [NtF 2 ] − anion.