VirialVLE:Virial Basic Layout
From Wiketomica
| Introduction |
| Background |
Basic Layout
|
| Examples |
| Problems |
| Run Simulation |
| Assessment |
| References |
| Credits |
| Etomica Modules |
When the Virial component of the application is started, it looks like this:
The left side of the display is used for numerical data input (B at different temperatures). At the bottom are
three sliders and input fields which are used to adjust the basic parameters of the LJ potential (
(“epsilon”) and σ (“sigma”) according to Eq. (4)) and the quadrupole moment Q. At the center of the interface is the
Function Display where the analytical B(T ) calculated from Eq. (7) and additional quadrupole terms is
displayed together with the data points.
Data points that are entered in the data input field are automatically added as points in the Function
Display. If any of the sliders is moved the function B(T) will change accordingly. By adjusting the potential
parameters one can bring the calculated B(T) curve to agree with the experimental data, as shown in the
screenshot below:
This way, the LJ parameters and the quadrupole moment can be obtained from experimentally determined
virial coefficients B. These potential parameters can be used for subsequent Molecular Dynamics (MD) or
Monte Carlo (MC) simulations.
The unit for σ is
m. The parameter
corresponds to the depth of the LJ potential well
and therefore measures the interaction energy at the molecular level. The usual energy unit in SI is Joule
(J) but the value for the interaction between two molecules is tiny compared to a Joule. The Boltzmann
constant k is of “atomic dimension”; therefore
is often given in units of kT and one simply specifies a
temperature T instead. Therefore, the module reads the
energy value in Kelvin. The dimension of the
quadrupole moment is
. Again, SI units are adapted to macroscopic objects and lead to very
small numerical values at the molecular scale. Here, the quadrupole moment is given in units of
with 1 debye =
esu·cm. The unit esu is a unit of electrical charge and comes from the old CGS unit system (Centimeter–Gram–Second). The corresponding SI unit is
.


