[ase-users] another nearest neighbor list question

Jens Jørgen Mortensen jensj at fysik.dtu.dk
Mon Feb 21 15:28:15 CET 2011


On Thu, 2011-02-17 at 10:19 -0500, John Kitchin wrote:
> That doesn't make sense to me. There are 12 atoms that are equidistant
> from the origin. The halflist shouldn't apply or change the results of
> this case because there is only one atom in the unit cell.
> 
> If I make a 3x3x3 bulk fcc cell,  and get the neighbor list of the
> central atom, then it returns only 8 atoms even with bothways=True. I
> would have expected 12 atoms.

Yes, you should get 12.  Fixed in trunk.

Jens Jørgen

> John
> 
> -----------------------------------
> John Kitchin
> Assistant Professor
> Doherty Hall A207F
> Department of Chemical Engineering
> Carnegie Mellon University
> Pittsburgh, PA 15213
> 412-268-7803
> http://kitchingroup.cheme.cmu.edu
> 
> 
> 
> 
> 
> On Thu, Feb 17, 2011 at 9:49 AM, Tao Jiang <tjiang at fysik.dtu.dk> wrote:
> > I think this is the half list as Jakob mentioned before. By doing so,
> > the calculator such as emt does not need to go through the whole
> > list twice to calculate the total energy.
> >
> > /TJ
> >
> >
> > On Feb 17, 2011, at 3:25 PM, John Kitchin wrote:
> >
> >> for the nearest neighbor list, would you expect to get 12 results for
> >> the atom at the origin of an FCC bulk lattice (not counting the atom
> >> itself)? I only get 6 neighbors. Am I misunderstanding the results, or
> >> how to setup the list?
> >>
> >> from ase.lattice.cubic import FaceCenteredCubic
> >> from ase.calculators.neighborlist import *
> >>
> >> atoms = FaceCenteredCubic(directions=[[0,1,1], [1,0,1], [1,1,0]],
> >>                          size=(1,1,1), symbol='Cu', pbc=(1,1,1))
> >>
> >>
> >> nl = NeighborList([2.6/2.],bothways=False,self_interaction=False)
> >> nl.update(atoms)
> >> indices, offsets = nl.get_neighbors(0)
> >>
> >> print indices
> >> print offsets
> >>
> >> [0 0 0 0 0 0]
> >> [[ 0  0  1]
> >> [ 0  1 -1]
> >> [ 0  1  0]
> >> [ 1 -1  0]
> >> [ 1  0 -1]
> >> [ 1  0  0]]
> >>
> >>
> >> John
> >>
> >> -----------------------------------
> >> John Kitchin
> >> Assistant Professor
> >> Doherty Hall A207F
> >> Department of Chemical Engineering
> >> Carnegie Mellon University
> >> Pittsburgh, PA 15213
> >> 412-268-7803
> >> http://kitchingroup.cheme.cmu.edu
> >> _______________________________________________
> >> ase-users mailing list
> >> ase-users at listserv.fysik.dtu.dk
> >> https://listserv.fysik.dtu.dk/mailman/listinfo/ase-users
> >
> >
> 
> _______________________________________________
> ase-users mailing list
> ase-users at listserv.fysik.dtu.dk
> https://listserv.fysik.dtu.dk/mailman/listinfo/ase-users




More information about the ase-users mailing list