動物細胞上有無數的離子通道，每種通道會控制細胞膜電壓差異並允許對應的離子通過，藉此維持生理功能正常。由於離子通道實在太小了，要觀察並記錄它的各種數據需要大量的時間和昂貴的設備，因此我們使用
泊松-能斯特-普朗克-費米 (PNPF) 模型模擬陽離子通過此通道時的情形。
該方法能分析通道蛋白質通道內的分子表面、奇異電荷並展現出模擬離子流動情形上的重要特徵，如最佳集中、有效率的非線性迭代和物理特性。PNPF 模型模擬 CaVab 通道在不同電壓下所產生的電流與實驗的電流-電壓取得數據上的一致。另一個模擬方法 MOL 則是將通道內重要的點位加上離子，經過計算後得出離子是否在該位置會穩定存在，用以展現出CaVab 通道的推擠與離子進出細胞的方式，在結果上也與實現數據相當的吻合。

Animal cells have numerous ion channels, each channel controls cellular membrane voltage difference and allows corresponding ions to pass through the channel by this way to keep the normal physiological function. Because the size of ion channels is too tiny, observing and recording its various data needs a lot of time and expensive equipment. Because of that, we simulate the cation transport through CaVab by using the Poisson-Nernst-Planck-Fermi (PNPF) model. The method can analyze molecular surfaces and singular charges of channel proteins and exhibit important features in flow simula-tions such as optimal convergence, efficient nonlinear iterations, and physical conservation. The PNPF currents are in accord with the experimental I-V (current-voltage) data of the CaVab channel with various voltage. Another simulation method using atomic PF theory adds binding sites in the channel with ions. After calculation, it is concluded that the ion will be stable or not at the binding site. It shows that the theory can describe a "knock-off" mech-anism of calcium selectivity and the transport of ions through the channel.

[1] Bazant, M. Z., Storey, B. D., and Kornyshev, A. A. Double Layer in Ionic Liquids: Overscreening versus Crowding. Physical Review Letters 106, 4 (2011), 6–9.
[2] Berman, H. M., Battistuz, T., Bhat, T. N., Bluhm, W. F., Bourne, P. E., Buckhardt, K., Feng, Z., Gilliland, G. L., Iype, L., Jain, S., Fagan, P., Marvin, J., Padilla, D., Ravichan-dran, V., Schneider, B., Thanki, N., Weissig, H., Westbrook,J. D., and Zardecki, C. The Protein Data Bank. Acta crystallograph-ica D 58 (2002), 899–907.
[3] Chen, D., and Eisenberg, R. Charges, currents, and potentials in ionic channels of one conformation. Biophysical journal 65, 2 (1993), 727–46.
[4] Liu, J. L. Numerical methods for the Poisson –Fermi equation in electrolytes. Journal of Computational Physics 247 (2013), 88–99.
[5] Liu, J. L., and Eisenberg, B. Correlated ions in a calcium channel model: A Poisson-Fermi Theory. Journal of Physical Chemistry B 117, 40 (2013), 12051–12058.
[6] Liu, J. L., and Eisenberg, B. Poisson-Nernst-Planck-Fermi theory for modeling biological ion channels. The Journal of Chemical Physics 141, 22 (2014), 22D532.
[7] Liu, J. L., and Eisenberg, B. Numerical methods for a Poisson-Nernst-Planck-Fermi model of biological ion channels. Physical Review E 92, 1 (2015), 12711.
[8] Tang, L., Tamer M. Gamal El-Din, J. P., Gilbert Q. Martinez, T. M. H., Todd Scheuer, N. Z., and Catterall, W. A. [10]2014 (Nature) Structural basis for Ca selectivity. Nature, 1 (2014), 1–5.
[9] Todd J., D., Paul, C., Li, H., Nielsen, J. E., Jensen, J. H.,
Klebe, G., and Baker, N. A. PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simula-tions. Nucleic Acids Research 35, SUPPL.2 (2007), 522–525.