How can I calculate the lattice parameter for building Argon atoms using LAMMPS?

  • Context: Undergrad 
  • Thread starter Thread starter positroncascade
  • Start date Start date
  • Tags Tags
    Flow
Click For Summary
SUMMARY

This discussion focuses on calculating the lattice parameter for modeling Argon atoms using LAMMPS in a 3D Poiseuille flow simulation. The user is attempting to determine the correct number of Argon atoms needed for their model, referencing the Knudsen number and bulk density equations. The key variables include a diameter (D) of 103 Angstroms and a bulk density of Argon at 1.6 kg/m³. The user seeks assistance in calculating the lattice parameter to ensure the expected number of atoms aligns with the simulation results.

PREREQUISITES
  • Understanding of molecular dynamics (MD) simulations
  • Familiarity with LAMMPS commands, specifically the lattice fcc command
  • Knowledge of the Knudsen number and its application in gas dynamics
  • Basic concepts of bulk density and mean free path in gases
NEXT STEPS
  • Research how to calculate the lattice parameter for Argon using the fcc structure in LAMMPS
  • Explore the implications of the Knudsen number on molecular dynamics simulations
  • Learn about the relationship between bulk density and number density in gas simulations
  • Investigate the use of the mean free path formula in determining atom counts in simulations
USEFUL FOR

This discussion is beneficial for molecular dynamics researchers, LAMMPS users, and anyone involved in simulating gas flows and atomic structures in computational materials science.

positroncascade
Messages
2
Reaction score
0
TL;DR
I try to model Poiseuille flow of Argon gas confined between two parallel Platinum plates. However, I am not able to obtain the density of Ar at the end of my simulation. So, probably, I am doing something wrong and I need help to calculate lattice parameter.
Hello everyone. I am newbie at MD. I started with Poiseuille flow.

I try to model 3D Poiseuille flow of Argon gas confined between two parallel Platinum plates by using LAMMPS. I have one reference paper, as given in [1] below. I built my Platinum structure correctly and I placed Ar atoms and run my simulation. At the end of simulation, I checked the density and it was totally wrong. I suspect that, I am not able to build sufficient number of atoms in my model. I try to use equations given below to obtain number of atoms but I couldn't decide the number of Ar atoms in my model, still. So, I need your help.

In the paper, they obtain the number of atoms by using Knudsen number. So, from the definition of Knudsen,

$$
Kn = {\lambda \over D}
$$

and

$$
\lambda = {1 \over \sqrt{2} \Pi d^2 {n_0}}
$$

where ## {n_0} ## is bulk number density. The two things that I know are D, which is 103 Angstrom and the bulk density of Argon, which is 1.6 ## kg \over {m^3} ## .

I want to build Ar atoms using lattice fcc command of Lammps and for that, I need lattice parameter. How can I calculate it and match the number of atoms that I expect and that I get from Lammps?

I would be very happy if you can help me which means a lot for me as I am stuck at this point.

Thanks for your help, indeed.The reference that I use:

[1] https://doi.org/10.1080/15567265.2016.1215364
 
Last edited:
Physics news on Phys.org
Usual formula for mean free path is
{\displaystyle \ell ={\frac {k_{\text{B}}T}{{\sqrt {2}}\pi d^{2}p}},}
where ##k_B## is the Boltzmann constant. Comparing it with your fomula it seems
n_0=\frac{p}{k_B T}=\frac{N}{V}
Does it meet your "bulk density" though its dimension is ##L^{-3}## number density not ##ML^{-3}## mass density?
 
Last edited:

Similar threads

  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 6 ·
Replies
6
Views
5K
Replies
1
Views
5K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 15 ·
Replies
15
Views
6K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 7 ·
Replies
7
Views
5K
  • · Replies 75 ·
3
Replies
75
Views
7K
  • · Replies 3 ·
Replies
3
Views
1K
Replies
9
Views
3K