LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] ReaxFF form in Ashraf and van Duin (JPCA 2016, 10.1021/acs.jpca.6b12429)
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] ReaxFF form in Ashraf and van Duin (JPCA 2016, 10.1021/acs.jpca.6b12429)


From: Luca Gelisio <luca.gelisio@...24...>
Date: Fri, 11 Aug 2017 23:51:26 +0200

Hey Ray,
Thanks a lot!
Luca

On Fri, Aug 11, 2017 at 10:51 PM, Ray Shan <rshan@...1795...> wrote:

Hi Luca, Michal,

 

Indeed, removing Delta_boc from lines 152/153 in reaxc_valence_angles.cpp would reproduce the Ashraf_2017 paper for modified valence angle term, but simply removing this variable would break backward compatibility from all previous ReaxFF descriptions.

 

A better way to do this is to add a keyword, e.g. modval, in the pair_style reax/c argument – as shown in the attached pair_reaxc.html.  With the keyword set to yes, the valence angle takes the modified form. 

 

The attached tar ball contains all necessary changes, including the doc page, reax/c source code for the original, kokkos, and omp versions.  Untar it in the top level LAMMPS directory to include the changes.

 

Thanks,

Ray

 

--

T. Ray Shan, Ph.D.   |   Support Scientist

MATERIALS DESIGN - San Diego, California, USA

Direct +1.760.495.4924, ext. 219www.materialsdesign.com

 

 

From: Luca Gelisio <luca.gelisio@...33....24...>
Date: Friday, August 11, 2017 at 8:50 AM
To: Giacomo Fiorin <giacomo.fiorin@...24...>
Cc: LAMMPS Users Mailing List <lammps-users@...396...sourceforge.net>
Subject: Re: [lammps-users] ReaxFF form in Ashraf and van Duin (JPCA 2016, 10.1021/acs.jpca.6b12429)

 

Hi Giacomo and Michal,

Yes, they've changed the valence angle term and I think Michal is right, that patch should be enough.

Thanks a lot,

luca

 

 

On Fri, Aug 11, 2017 at 4:08 PM, Giacomo Fiorin <giacomo.fiorin@...24...> wrote:

I see, thanks for the correction.

 

On Fri, Aug 11, 2017 at 9:57 AM, Michał Kański <michal.kanski@...2460...> wrote:

Hi Giacomo,

 

They’ve changed the valence angle term by removing the undercoordination term. See the equations 3 and 4 in the paper.

The modification isn’t big and I think one only needs to delete “workspace->Delta_boc[j]” from lines 152 and 153 in reaxc_valence_angles.cpp in order to introduce it. Someone with a larger experience with ReaxFF should confirm this though.

 

Michal

 

 

 

Od: Giacomo Fiorin
Wysłano: piątek, 11 sierpnia 2017 15:07

Do: Luca Gelisio
DW: LAMMPS Users Mailing List
Temat: Re: [lammps-users] ReaxFF form in Ashraf and van Duin (JPCA 2016, 10.1021/acs.jpca.6b12429)

 

Hi Luca, it looks like only the parameters set has been optimized, but the energy function is still the same.  Can you try running models with the force field file in the SI of this paper?

Giacomo

 

On Fri, Aug 11, 2017 at 8:14 AM, Luca Gelisio <luca.gelisio@...24...> wrote:

Dear LAMMPS users and developers,

 

According to LAMMPS reax/c webpage (http://lammps.sandia.gov/doc/pair_reaxc.html), the Reax version currently implemented in LAMMPS is the one reported in Chenoweth et al (JPCA 2008, 10.1021/acs.jpca.6b12429).

By chance, has someone already experimented with the version reported in Ashraf and van Duin (JPCA 2016, 10.1021/acs.jpca.6b12429)?

 

Thanks for any comment,

luca

 


------------------------------------------------------------------------------
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@...396...sourceforge.net
https://lists.sourceforge.net/lists/listinfo/lammps-users




--

Giacomo Fiorin

Associate Professor of Research, Temple University, Philadelphia, PA

Contractor, National Institutes of Health, Bethesda, MD

http://goo.gl/Q3TBQU
https://github.com/giacomofiorin

 




--

Giacomo Fiorin

Associate Professor of Research, Temple University, Philadelphia, PA

Contractor, National Institutes of Health, Bethesda, MD