diff --git a/potentials/README b/potentials/README index d4dba824c..599432ce8 100644 --- a/potentials/README +++ b/potentials/README @@ -1,103 +1,104 @@ This directory contains potential files for different elements and alloys, as used by LAMMPS for various pair styles. See the description of the "pair_style" and "pair_coeff" commands for details of the file formats and the various styles in LAMMPS that read these files. IMPORTANT NOTE: These files are provided primarily to demonstrate the different types of interatomic potentials that LAMMPS supports. Each file has a header line with a date for when it was added to the LAMMPS distribution. Also a citation and contact info for the person who contributed it to LAMMPS (if we remember who that is). This info is not meant to "guarantee" that the potential is correct. I.e. that the contributor transcribed the info from the paper correctly or that the paper itself had no errors. In many cases (but not all), we or other LAMMPS users have confirmed that when the potential file is used with the current version of LAMMPS, it reproduces results in the cited publication. In some cases, this accuracy check may require other parameters not contained in the potential file to be specified as part of a LAMMPS simulation, e.g. a distance cutoff. Also, for particular materials and applications modeled with a pair style coded in LAMMPS, a different potential file may be more suitable than the one provided here. For best results when choosing a potential, you should do a thorough search of published literature and on-line databases such as the Interatomic Potentials Repository Project (NIST) or the Knowledgebase of Interatomic Models (KIM). Whatever potential you choose for your application, you should verify that you have defined it and are using it correctly in LAMMPS, by comparing with published results for that potential. 2nd IMPORTANT NOTE: The DATE field in the first line of each of these files is printed to the screen and log file when it is read by a LAMMPS input script. If an updated or corrected version of the same potential file is later added to the LAMMPS distribution, then a new DATE will be added to the file. This means you can "diff" an old and new log file and see that the potential file changed, which could affect your simulation results. A small amount of metadata is included in the first line of each file in order to track the provenance of each file. The metadata is indicated by a keyword followed by white space, followed by the metadata, followed by whitespace. The metadata is intended to be straightforward and human-readable, while still conforming to a standard format. DATE: Format is "yyyy-mm-dd". This indicates the date of a significant change to the file. Multiple entries can appear in reverse chronological order. As described above, the first of these will be printed to the screen and log file when it is read by a LAMMPS input script. CONTRIBUTOR: Format is "name[, email address]". This indicates the person who contributed the file and/or who is best able to provide more details about its provenance. CITATION: Format is "surname[[, surname] and surname], Publication abbreviation with spaces and no periods, volume, page[-page], (year)" COMMENT: This one is optional and is used to hold any other text that can not go elsewhere. If the first line of the file is always skipped by the file reader, then the first line should begin with the DATE keyword. If the file format supports comment lines, then the first line should be a comment line with the metadata e.g. "# DATE: 2010-01-01..." If the first line of the file is required to begin with data, then the metadata will be appended to the first line e.g. "7 DATE: 2010-01-01..." The prefix of each file indicates the element(s) it is parameterized for. An additional lower-case identification tag may be appended. Si = Silicon SiC = Silicon and Carbon Au_u3 = Gold universal 3 The suffix of each file indicates the pair style it is used with: adp ADP angular dependent potential airebo AI-REBO and REBO potentials bop.table BOP potential, tabulated form cdeam concentration-dependent EAM comb COMB potential comb3 COMB3 potential eam embedded atom method (EAM) single element, DYNAMO funcfl format eam.alloy EAM multi-element alloy, DYNAMO setfl format eam.fs Finnis-Sinclair EAM format (single element or alloy) edip EDIP potential for silicon-based materials eim embedded-ion method (EIM) potential lcbop LCBOP long-range bond-order potential meam modified EAM (MEAM) library and individual elements/alloys meam.spline modified EAM (MEAM) spline potential meam.sw.spline modified EAM (MEAM) Stillinger-Weber spline potential nb3b.harmonic nonbonded 3-body harmonic potential poly polymorphic 3-body potential reax ReaxFF potential (see README.reax for more info) +smtbq Second Moment Tight Binding - QEq potential snap SNAP potential snapcoeff SNAP potential snapparam SNAP potential streitz Coulombic portion of Streitz-Mintmire potential sw Stillinger-Weber potential tersoff Tersoff potential tersoff.mod modified Tersoff potential tersoff.zbl Tersoff with ZBL core vashishta Vashishta 2-body and 3-body potential diff --git a/potentials/ffield.smtbq.Al b/potentials/ffield.smtbq.Al new file mode 100755 index 000000000..b8fd478a0 --- /dev/null +++ b/potentials/ffield.smtbq.Al @@ -0,0 +1,34 @@ +# +# SMTBQ parameter for Al-Al interaction. +# Edited by N. Salles Univ. Bourgogne and E. Maras from Aalto Univ. +# year: 2014 +# ========================================================================= +' Nombre.de.type.d.atome..........:' 1 +' ====== atomic parameters ======= ' +' Cation.de.l.oxyde..Stoechio.....:' 'Al' 1 +' Qform.....masse.................:' 3.0 26.98 +' Param.QEq.(ne,.Chi,.J,.R_eff)...:' 3 1.19258 11.05345 0.57701 +' Nbre.d.etats.partage.par.cation.:' 3 +' ===== potential Parameter ======' +' Atom1..atom2..potential..mode....' 'Al' 'Al' 'second_moment' 'metal' +' Potentiel.Cat-Ox.(A,.p,.Ksi,.q).:' 0.1221 8.612 1.316 2.516 +' Pot..Cat-Ox.(rc1,.rc2,.r0).Iota.:' 5.6 7.0 2.863 2.0 +' ======== Parametre tab ========= ' +' Rcoul...........................:' 11.1714 +' rmin...dr.......................:' 1.18845 0.001 +' ======== IFQM Parameter ======== ' +' Frenquency.Q.resolution..........' 0 +' loopmax.-.precision..............' 5000 0.0002 +' ==== Coordination parameters ====' +' .r1n................r2n..........' 2.5 3.2 +' ========== QInitMode ========= ' +' QInitMode....QInit(if.needed)....' 'false' 0.0 +' ======== Mode for QEq ======== ' +' mode(see.end.of.this.file.)......' 'QEqAllParallel' +' parameters.for.mode..............' +' ========== Verbose ============ ' +' Verbose(true.or.false)...........' 'false' +' Print.Energy.components..........' 'false' 300.0 +' Print.electroneg...components....' 'false' 300.0 +# =========================== END's parameters ========================= + diff --git a/potentials/ffield.smtbq.Al2O3 b/potentials/ffield.smtbq.Al2O3 new file mode 100755 index 000000000..9c26ade9e --- /dev/null +++ b/potentials/ffield.smtbq.Al2O3 @@ -0,0 +1,56 @@ +# DATE: 2015-10-22 CONTRIBUTOR: Nicolas Salles, nsalles@laas.fr CITATION: N. Salles, O. Politano, E. Amzallag and R. Tetot, Comput. Mater. Sci. 111 (2016) 181-189 +# SMTBQ parameter for AlO interaction with a limit length 'rc2sm=dc2**2'. +# Edited by N. Salles from Univ Bourgogne and E. Maras from Aalto Univ. +# +# Presentation atom : nature +# q, qmin, qmax, masse +# parameter QEq : Chi, J, R_eff +# Parameter SM : A, p, Ksi, q +# CutOff SM : dc1, dc2, r0 +# ========================================================================= +# -------------------------- Begin's parameters --------------------------- +' Number.of.atoms.type............:' 2 +' ====== atomic parameters ======= ' +' 1st.element.(Oxygen).Stoechio...:' 'O' 3 +' Qform.....mass..................:' -2.0 16.00 +' Param.QEq.(ne,.Chi0,.JiO).......:' 2 6.57 10.22 +' coordBB.coordB.coordS.rBB.rB.rS.:' 6. 4. 3.00 0.529 0.529 0.529 +' Number.of.shared.state.by.ions..:' 3 +' -------------------------------- ' +' 2nd.element.(metal).Stoechio....:' 'Al' 2 +' Qform.....mass..................:' 3.0 26.98 +' Param.QEq.(ne,.Chi0,.Ji0,.R_eff):' 3 1.19009 11.1903 0.56619 +' Number.of.shared.state.by.ions..:' 4 +' ===== potential Parameter ======' +' Atom1..atom2..potential..mode....' 'Al' 'O' 'second_moment' 'oxide' +' Pot.CatOx.(A,.p,.Ksi(ref=O),.q).:' 0.18176 8.80041 0.26044 1.58851 +' Pot..Cat-Ox.(rc1,.rc2,.r0)......:' 4. 5.6 1.91 +' -------------------------------- ' +' atom1..atom2..potential..........' 'O' 'O' 'buck' +' Potentiel.O-O...(C,.Rho)........:' 580.440 0.3540 +' ======== Parametre tab ========= ' +' Rcoul=a*rc(SMASH)...............:' 11.1714 +' rmin...dr.......................:' 1.18845 0.001 +' ======== IFQM Parameter ======== ' +' Nevery.charge.calculation........' 1 +' loopmax....precision.............' 7000 0.000001 +' ==== Coordination parameters ====' +' .r1n................r2n..........' 2.5 3.2 +' ========== QInitMode ========= ' +' QInitMode....QInit(if.needed)....' 'false' -1.8 +' ======== Mode for QEq ======== ' +' mode(see.end.of.this.file.)......' 'QEqAll' +' parameters.for.mode..............' +' ========== Verbose ============ ' +' Verbose(true.or.false)...........' 'false' +' Print.Energy.components..........' 'false' 300.0 +' Print.electroneg...components....' 'false' 300.0 +# =========================== FIN des parametres ========================= + +#Possible QInit modes +# true (then initialize all the oxygen charges to QOxInit and set the cation charge in order to keep the charge balance (neutrality of the box) +#any other name would lead to either 0 charges or charge read from the lammps atomic position file +#Possible QEq modes | parameters +# QEqAll | no parameters +# QEqAllParallel | no parameters +# Surface | zlim (QEq only for z>zlim) diff --git a/potentials/ffield.smtbq.TiO2 b/potentials/ffield.smtbq.TiO2 new file mode 100755 index 000000000..4c8d9df53 --- /dev/null +++ b/potentials/ffield.smtbq.TiO2 @@ -0,0 +1,53 @@ +# DATE: 2015-10-22 CONTRIBUTOR: Nicolas Salles, nsalles@laas.fr CITATION: N. Salles, O. Politano, E. Amzallag and R. Tetot, Comput. Mater. Sci. 111 (2016) 181-189 +# ======================================================================== +# SMTBQ parameter for Ti-O interaction with a limit length 'rc2sm=dc2**2'. +# Edited by N. Salles from Univ Bourgogne and E. Maras from Aalto Univ +# september 2014 +# ======================================================================== +' Number.of.atoms.type............:' 2 +' ====== atomic parameters ======= ' +' 1st.element.(Oxygen).Stoechio...:' 'O' 2 +' Qform.....mass..................:' -2.0 16.00 +' Param.QEq.(ne,.Chi0,.JiO).......:' 2 6.57 10.22 +' coordBB.coordB.coordS.rBB.rB.rS.:' 6. 3. 2.00 0.52 0.54348 0.58 +' Number.of.shared.state.by.ions..:' 3 +' -------------------------------- ' +' 2nd.element.(metal).Stoechio....:' 'Ti' 1 +' Qform.....mass..................:' 4.0 26.98 +' Param.QEq.(nq,.Chi0,.Ji0,.R_eff):' 3 0.00 10.572 0.734 +' Number.of.shared.state.by.ions..:' 5 +' ===== potential Parameter ======' +' Atom1..atom2..potential..mode....' 'Ti' 'O' 'second_moment' 'oxide' +' Pot.CatOx.(A,.p,.Ksi(ref=O),.q).:' 0.134 12.609 0.5434 2.0965 +' Pot..Cat-Ox.(rc1,.rc2,.r0)......:' 3.6 6.0 1.95 +' -------------------------------- ' +' atom1..atom2..potential..........' 'O' 'O' 'buckPlusAttr' +' Potential.O-O...(C,.Rho)........:' 580.440 0.3540 +' Potential.O-O...(D.B.r1OO.r2OO).:' -20.86 -0.916 1.4 1.8 +' ======== Tab Parameter ========= ' +' Rcoul=a*rc(SMASH)...............:' 12.1744 +' rmin...dr.......................:' 1.0675 0.001 +' ======== IFQM Parameter ======== ' +' Nevery.charge.calculation........' 1 +' loopmax....precision.............' 7000 0.000001 +' ==== Coordination parameters ====' +' .r1n................r2n..........' 2.0 3.5 +' ========== QInitMode ========= ' +' QInitMode....QInit(if.needed)....' 'false' -1.0 +' ======== Mode for QEq ======== ' +' mode(see.end.of.this.file.)......' 'QEqAll' +' parameters.for.mode..............' +' ========== Verbose ============ ' +' Verbose(true.or.false)...........' 'false' +' Print.Energy.components..........' 'false' 300.0 +' Print.electroneg...components....' 'false' 300.0 +# =========================== END's parameters ========================= + +#Possible QInit modes +# true (then initialize all the oxygen charges to QOxInit and set the cation charge in order to keep the charge balance (neutrality of the box) +#any other name would lead to either 0 charges or charge read from the lammps atomic position file +#Possible QEq modes | parameters +# QEqAll | no parameters +# QEqAllParallel | no parameters +# Surface | zlim (QEq only for z>zlim) +# BulkFromSlab | zlim1 zlim2 (QEq only for zlim1