Publications Search

Molecular dynamics simulation was implemented using LAMMPS simulation package (1) to study the diffusivity of He³ and CO₂ in NaCl aqueous solution. To simulate at infinite dilute gas concentration, we placed one He³ or CO₂ molecule in an initial simulation box of 24x24x33Å³ containing 512 water molecules and a certain number of NaCl molecules depending on the concentration. Initial configuration was set up by placing water, NaCl, and gas molecules into different regions in the simulation box. Calculating diffusion coefficient for one He or CO₂ molecule consistently yields poor results. To overcome this, for each simulation at specific conditions (i.e., temperature, pressure, and NaCl concentration), we conducted 50 simulations initiated from 50 different configurations. These configurations are obtained by performing the simulation starting from the initial configuration mentioned above in the NVE ensemble (i.e., constant number of particles, volume, and energy). for 100,000 time steps and collecting one configuration every 2,000 times step. The output temperature of this simulation is about 500K. The collected configurations were then equilibrated for 2ns in the NPT ensemble (i.e., constant number of particles, pressure, and temperature) followed by 9ns simulations in the NVT ensemble (i.e., constant number of particles, volume, and temperature). The time step is 1fs for all simulations.

More Details

TYPE Other Report YEAR 2016

OSTI DOI