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

  • Thread starter Thread starter positroncascade
  • Start date Start date
  • Tags Tags
    Flow
Click For Summary
The discussion revolves around calculating the lattice parameter for modeling Argon atoms in LAMMPS while simulating 3D Poiseuille flow between Platinum plates. The user is struggling with incorrect density results, suspecting an insufficient number of Argon atoms in their model. They reference the Knudsen number and provide equations for calculating the mean free path and bulk number density, seeking guidance on determining the correct number of Argon atoms. The user specifically wants to use the lattice fcc command in LAMMPS and is looking for help in calculating the lattice parameter to align their expected atom count with simulation results. Assistance is requested to resolve these issues and advance their simulation work.
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
5K
  • · Replies 15 ·
Replies
15
Views
5K
Replies
6
Views
2K
  • · Replies 75 ·
3
Replies
75
Views
7K
Replies
9
Views
3K
Replies
7
Views
9K
  • · Replies 11 ·
Replies
11
Views
4K