[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

From: |
John Smith <johnson21smith@...24...> |

Date: |
Fri, 23 Jun 2017 17:26:48 -0500 |

Yes I ran both of them and there was a difference in the output for both methods (~0.3 using vacf and ~0.2 using msd). However, in my calculations, they don't even have the same number of significant digits. I followed the exact same procedure as given in the examples/DIFFUSE directory and took into consideration the 2D factors as well.

What exactly do you mean by long-time snapshots?

Sincerely,

John

On 23 June 2017 at 09:38, Steve Plimpton <sjplimp@...24...> wrote:

Stevewith VACF and proper averaging and statistics, etc. MSD isway to compute a diffusion coeff. You have to worry about timescalesabout those than your scripts. I will say that VACF is a much trickierHave you tried running the MSD and VACF input scriptin examples/DIFFUSION? If so, try to see what is different

simple. There is little that can go wrong. And you can validatewith just a couple long-time snapshots, that you have a validMSD number as produced by LAMMPS.On Thu, Jun 22, 2017 at 4:55 PM, John Smith <johnson21smith@...24...> wrote:I also tried the same for the mixture simulation and calculated the VACF for the non-SRD particles. In this case, there are interactions between the particles and hence velocity correlations should exist, right? However, this is not what I am observing as the diff coefficient calculated using VACF and MSD dont match.I am once again attaching the input script and two graphs.Sincerely,JohnOn 22 June 2017 at 08:39, Steve Plimpton <sjplimp@...24...> wrote:Steveand how to interpret it, is a different Q.is the same. What that means for an SRD fluidat some time delta. I'm guessing you will find itone or more dump files and compute the MSD yourselfLAMMPS is producing a valid number. Just outputMSD is incorrect, you can check yourself ifSince SRD particles do not interact with each other,I don't think a VACF is meaningful. Effectively

there are no velocity correlations.

When you sayOn Wed, Jun 21, 2017 at 5:40 PM, John Smith <johnson21smith@...24...> wrote:------------------------------Hello,I have tried to modify the example script given under SRD to calculate vacf and mean squared displacement and I am attempting to calculate the diffusion coefficient from this. However, the values are not matching and the range is different. I also included the factors necessary for a 2D simulation (details given in examples/DIFFUSE) but this is still giving an incorrect result. I am attaching my input script as well as the output graphs. Please advice!Sincerely,John------------------------------ ------------------

Check out the vibrant tech community on one of the world's most

engaging tech sites, Slashdot.org! http://sdm.link/slashdot

_______________________________________________

lammps-users mailing list

lammps-users@...655....net

https://lists.sourceforge.net/lists/listinfo/lammps-users

**Follow-Ups**:**Re: [lammps-users] Diffusion Coefficient in 2D***From:*Steve Plimpton <sjplimp@...24...>

**References**:**[lammps-users] Diffusion Coefficient in 2D***From:*John Smith <johnson21smith@...24...>

**Re: [lammps-users] Diffusion Coefficient in 2D***From:*Steve Plimpton <sjplimp@...24...>

**Re: [lammps-users] Diffusion Coefficient in 2D***From:*John Smith <johnson21smith@...24...>

**Re: [lammps-users] Diffusion Coefficient in 2D***From:*Steve Plimpton <sjplimp@...24...>

- Prev by Date:
**Re: [lammps-users] Question about Pair Coeff for Beads of same atoms types** - Next by Date:
**[lammps-users] Different outcomes with same *.data & *.in while changing the parallel core numbers** - Previous by thread:
**Re: [lammps-users] Diffusion Coefficient in 2D** - Next by thread:
**Re: [lammps-users] Diffusion Coefficient in 2D** - Index(es):