Interfacial tension:Example 1

From Wiketomica
Basic Layout
Run Simulation
Etomica Modules

Simulation Example WITHOUT Surfactants

Re-initialize the simulation. If a simulation has not been run since the application has been started, the simulation will already be initialized and the first button in the Control panel will be labeled “Start”. If the simulation has been previously started, even if it is currently paused, push the Reinitialize button. The “Continue” button will change to “Start” again upon re–initialization.

Set the simulation parameters to :

  • Temperature () = 1.0
  • Box size = 12.0
  • Number of atoms = 643
  • Number of surfactants = 0.0
  • (Surfactant potential parameters not used)

Start the simulation by pushing the Start button. The button will change its label to “Pause” which can be used to pause and continue the simulation at any given time.

As the simulation proceeds, observe the atoms diffusing within the cubic simulation box at a liquid-like density. The density profile can be monitored (which should have an average value of around 0.54), as well as the potential energy and the other simulation metrics. Once the simulation has been allowed to run for approximately 100 ps, the liquid phase should be fairly well stabilized. On a modern PC, this should take approximately 1 minute. At this point, the fluid will not necessarily be fully equilibrated, but it should be stable enough to proceed to the next step, which begins the two-phase simulation.

At this point, the "Expand" button can be pressed, which increases the x-dimension of the simulation box by a scale-factor of 3. You can observe under the Configuration tab that the atoms from the liquid-like region begin evaporating into the newly-created void spaces in the expanded simulation box. It takes approximately 100 ps of simulation time in order to stabilize the two-phase system. Once this point has been reached (a total simulation time of 200 ps or more), the "Reset Averages" button can be pushed, which will reset the plots and the data collection for the averages will be reinitialized. The progress of the simulation will be unaffected by resetting the averages.

The “Error” printed for each property in the Metrics tab is the statistical error of the running average; it does not necessarily indicate if the results are converged. For example, there could be a steady drift in the surface tension despite small magnitudes of the Interfacial Tension Error listed in the Metrics tab. Below are some screen shots of what a simulation typically looks like after roughly 2000 ps (following the equilibration period), as well as representative data profile plots:

IT pic4.gif

  • Density Tab

VL interface.gif

  • Virial Profile


  • Tension Profile

ST profile.gif

It takes approximately 1500 ps (roughly 6 to 8 min on a modern PC) for the most drastic variations in the properties to disappear. Use more steps if necessary to diminish statistical errors in the computed properties and to make sure there are no systematic drifts in the calculated properties (this is easily revealed by the various plots). You will find that the vapor and the liquid density converge rather quickly (found under the Density tab) at temperatures well below the critical point, but that the virial, interfacial tension, and chemical potential profiles may require more steps to converge. One of the key features to examine is the symmetry of the various profile plots. Although there will still be fluctuations, they all should essentially become symmetrical as the simulation converges. Once converged, the statistics from the simulation can be documented, such as the surface tension, the liquid and gas phase densities, and the potential energy. These results are unique to the parameters specified for a given simulation, including the thermodynamic conditions and the intermolecular parameters.