Molecular Simula-on studies of Organic salts inside carbon nanopores Harsha Dissanayake Dr. Francisco R. Hung Dr. Joshua D. Monk Ramesh Singh Cain Department of Chemical Engineering, Louisiana State University, Baton Rouge, LA, USA
What are Ionic liquids? An ionic liquid is a salt in liquid phase at room temperature. Made up of ions and have ionic bonds. Posses low vapor pressure and are moderate or low electric conductors. Highly tunable properies; 10 9 10 18 ILs could be formed by varying caions, anions. 1,3 dimethylimidazolium chloride [DMIM+][Cl ] Vast number of applicaions due to disinct properies in many fields such as Organic chemistry, Engineering, Electro chemistry, Catalysis and Physical chemistry. (Solvents, electrolytes, etc)
Importance of Ionic Liquids(IL s) IL s can be synthesized to possess fluorescent, magneic and many other qualiies according to necessity. Ionic liquids can be synthesized from safe, FDA approved compounds; therefore could be used in Bio medical applicaions. Nanomaterials with fluorescent proper-es can be used to label Issues and cells for biomedical imaging Nanomaterials with magne-c proper-es can be used in cancer treatment (magneic hyperthermia). Aber acaching them to cancer cells, the paient can be exposed to a safe magne-c field and the heat generated by the magneic nanoparicles will destroy the respecive cancer cells that they are acached to. Image courtesy of the Warner group [PF 6ˉ ] 1 butyl 2,3 methylimidazolium hexafluorophosphate, [bm 2 im + ][PF 6ˉ ]
Why molecular simulaion? The properies of ionic liquids confined inside nanopores can be very different from what we see in day to day usage (bulk systems). It is extremely difficult to extract experimental data regarding such systems Molecular dynamics (MD) simulaion is a very useful technique in such cases to give insights to experimental scienists and to complement their observaions regarding properies and characterisics of molecules at nanolevel.
Overview of simula-on methods TIME /s 10 0 (ms)10-3 (µs)10-6 (ns)10-9 (ps)10-12 (fs)10-15 Atomistic simulation methods Semi-empirical methods Ab initio methods Mesoscale methods tight-binding MNDO, INDO/S Lattice Monte Carlo Brownian dynamics Dissipative particle dynamics Monte Carlo Molecular dynamics Continuum Methods 5 10-10 10-9 10-8 10-7 10-6 10-5 10-4 (nm) (µm) LENGTH/meters
Molecular dynamics: a quick introduc-on In molecular dynamics (MD), we follow the posiions and velociies of each molecule in Ime by solving Newton s equaions of moion: i=1,2,3, N U = inter/intra molecular potenial (i.e., interacions between atoms) r i = posiion vector of atom i F i = force acing over atom i m i = mass of atom i In MD these equations are integrated numerically to obtain the time evolution of the system 6
Pair Poten-als: Large, Flexible Mols. Total pair energy breaks into a sum of terms: DNA U str stretch U bend bend U tors torsion U disp dispersion (van der Waals) U elec electrostaic U pol polarizaion 7
System to be studied 1,3 dimethylimidazolium chloride [DMIM+][Cl ] Molar mass 132.5 mols/gram MulI walled carbon nanotube hcp://www.photon.t.u tokyo.ac.jp/~maruyama/agallery/ nanotubes/mwnt.gif In my summer project I will be invesigaion various properies of the ionic liquid [DMIM+][Cl ] when it is confined in a muli walled carbon nanotube. I will be using Gromacs MD sobware to perform MD simulaions of [DMIM+][Cl ] Maolin Sha, Guozhong Wu, Haiping Fang, Guanglai Zhu, and Yusheng Liu J. Phys. Chem. C, 2008, 112 (47), 18584 18587 DOI: 10.1021/jp8079183 Publica?on Date (Web): 30 October 2008
[DMIM+][Cl ] molecules inserted into carbon nanotube System without [Cl ] atoms
IniIal configuraion Minimized Energy state Melted state The [Cl ] ion has not been depicted for visualizaion purposes.
0 200 400 600 800 Density Density 10 per. Mov. Avg.(Density) 0 0.2 0.4 0.6 0.8 1 1.2 0 700 1400 2100 2800 3500 4200 4900 5600 6300 7000 7700 rmsd rmsd 10 per. Mov. Avg. (rmsd) 0 0.5 1 1.5 2 2.5 3 3.5 4 0 0.078 0.156 0.234 0.312 0.39 0.468 0.546 0.624 0.702 0.78 0.858 0.936 1.014 1.092 1.17 1.248 1.326 1.404 1.482 rdf rdf 10 per. Mov. Avg.(rdf) Aber running the system under bath condiions for approximately 10ns we can extract the data and analyze it to understand the molecular level properies and characterisics of the system. This informaion can be useful to experimentalists.
End
Molecular dynamics: Introduc-on MD simula5ons are similar to real experiments In a basic MD simulaion program: Prepare the sample : IniIal energy (or temperature) Number of paricles N Box size (N/L 3 = density) Force field equaions, parameters Time step dt for integraing equaions of moion IniIalize posiions and velociies of all paricles Do the experimental measurement : Compute forces on all paricles Integrate Newton s equaions (F = ma) Update paricle posiions and velociies Calculate instantaneous properies Stop aber iteraing t max steps Repeat unil we reach an appropriate simulated Ime 13