after opening the correct folder in the terminal and running the code, lmp_serial -in in.lmp, the file is not opening. is there a specific way the files should be organized so that lammps is able to read it or is there another issue here.

Describe the bug
As https://icme.hpc.msstate.edu/mediawiki/index.php.1.html was archived in favour of https://www.cavs.msstate.edu/icme/ the link to certain files no longer work.

To Reproduce
https://icme.hpc.msstate.edu/mediawiki/index.php/Fe_110_sig3.txt

Expected behavior
It could be replaced with a link to the Github repository:
https://raw.githubusercontent.com/mrkllntschpp/lammps-tutorials/master/Fe_110_sig3.txt

Additional context
Maybe other links are effected as well.

Hello,
Thanks for all these great tutorials! I got stuck with the delete_atoms command. I am looking at the LAMMPS documentation and comparing codes and I can't see where I am going wrong, yet I continue to get the error "ERROR: Could not find delete_atoms group ID (../delete_atoms.cpp:306)
Last command: delete_atoms overlap 0.35 region lower upper"
I get this error right after it starts the displacing atoms section. I copied my code directly from your code posted on this tutorial, so I am confused with what is happening.

Dear sir,

Thanks for your input file to calculate the energy of grain boundary.
I have an question in tutorial-5.
Why the minimumenergy is euqal to -3.36? This item should be the atomic energy in bulk. How did you get this value?

Kind regards,
Meng Zhang

Hello,

Thank you for providing such comprehensive list of tutorials for LAMMPS. I am just starting out with LAMMPS and I am finding your tutorial very useful.

In tutorial 1, you show how to perform energy minimization and obtain the optimized lattice parameter for Al lattice.

I was wondering, if there is way to perform energy minimization in LAMMPS where only the atoms are free to relax while the lattice shape and volume doesn't change. If you are familiar with VASP input, this would be equivalent to setting ISIF = 2.

I tried looking it up online and so far I haven't found any direct method of doing it.

I want to know if such energy minimizations are possible in LAMMPS and if so could you please give me some hints on how to do it.

Thank you!

Thank you very much for this tutorial. According to this tutorial, I have realized the tensile mechanical properties of BCC(body centered cubic) Fe block (x - < 100 >, Y - < 010 >, Z - < 001 >).

But I couldn't achieve it when I tried the tensile mechanical properties of another crystal oriented Fe block (x - < 110 >, Y - < 110 >, Z - < 001 >). I have tried to modify the system temperature, cutoff
and other parameters, but I still can't. I have to ask you for help. I hope you can help me.

Thank you!

Dear Mark,

I am currently learning to utilize Lammps software for atomistic simulations to extract material parameters necessary for continuum modeling of materials and alloys, such as employing crystal plasticity. I am writing to seek your assistance regarding some doubts I have encountered while interpreting the simulation results from the Lammps log file. Your Tutorial - 01 has been immensely helpful in understanding the intricacies of molecular dynamics simulations.

However, I have come across a few sections in the log file that I am having trouble comprehending. I was hoping you could kindly shed some light on them:

1. The following portion of the output appears in the log file, and I am unsure about its significance:

Nlocal:              2 ave           2 max           2 min
Histogram: 4 0 0 0 0 0 0 0 0 0
Nghost:            459 ave         459 max         459 min
Histogram: 4 0 0 0 0 0 0 0 0 0
Neighs:            140 ave         151 max         129 min
Histogram: 1 0 1 0 0 0 0 1 0 1

2. In Tutorial-01, a total of 8 atoms were modeled as indicated in the Lammps log file:

replicate 1 1 1
Replicating atoms ...
  orthogonal box = (0 0 0) to (5.6568542 5.6568542 4)
  2 by 2 by 1 MPI processor grid
  8 atoms
  replicate CPU = 0.004 seconds

However, later in the output file, it shows that the number of neighbors is equal to 560:

Total # of neighbors = 560
Ave neighs/atom = 70
Neighbor list builds = 0
Dangerous builds = 0

I am puzzled as to how the number of neighbors can be 560 for only 8 atoms in the simulation.

To facilitate your assistance, I have attached the input file and log file to this message.

Thank you for taking the time to review my inquiries. Your guidance will be immensely valuable to me in gaining a better understanding of these simulation results.

Best Regards,
Vikram Roy
Country: India
log.lammps.txt
calc_fcc.in.txt

How to write input script for single polymer chain of polycarbonate in for MD simulation in lammps?

Thank you for the tutorial. Have tried so far, but how to achieve polycrystalline stretching? It has always been a problem.

The polycrystalline structure uses random grains created by Atomsk software. The above is the structure file. I have tried calculations, but the performance at the grain boundary is not very good. So ask you for help.

Write to Fe_unitcell.xsf

Bcc Fe oriented X=[100] Y=[010] Z=[001].

CRYSTAL
PRIMVEC
2.85530000 0.00000000 0.00000000
0.00000000 2.85530000 0.00000000
0.00000000 0.00000000 2.85530000
CONVVEC
2.85530000 0.00000000 0.00000000
0.00000000 2.85530000 0.00000000
0.00000000 0.00000000 2.85530000
PRIMCOORD
2 1
26 0.00000000 0.00000000 0.00000000
26 1.42765000 1.42765000 1.42765000

Fe_voronoiRandom.txt
box 100 100 25
random 10

Fe_polycrystal.lmp is generated by the following command

atomsk --polycrystal Fe_unitcell.xsf Fe_voronoiRandom.txt Fe_polycrystal.cfg lmp -wrap

Happy New Year!
Sincerely

Before I begin, thanks a lot for making these tutorials available. I've a doubt here. In tutorial 3, after the initial equilibration, the simulation box is deformed

in the x-direction at a strain rate of 10^10 1/s

That is implemented in the LAMMPS script in the following lines

units 		metal
...
timestep 0.001  #in picosecond 
... 
variable srate equal 1.0e10
variable srate1 equal "v_srate / 1.0e12"
fix 2 all deform 1 x erate ${srate1} units box remap x

How do you calculate the factor of 1.0e12 ? Isn't the timestep dt involved here in some way?