0.8.2

Released 31 of March 2017.

Contributors

A total of 1 people contributed to this release. People with a “+” by their names contributed a patch for the first time.

  • Nick Papior

Pull requests merged

A total of 0 pull requests were merged for this release.

  • Fixed reading _hr.dat from Wannier90, now the band-structure of SrTiO3 (Junquera’s test example) is correct.

  • Speeded up tbtrans.py analyzing methods enourmously by introducing faster sparse iterators. Now one can easily perform data-analysis on systems in excess of 10.000 atoms very fast.

  • Added the TBT.AV.nc file which is meant to be created by sisl from the TBT.nc files (i.e. create the k-averaged output). This enables users to run tbtrans, create the k-averaged output, and then delete the old file to heavily reduce disk-usage.

    An example:

    tbtrans RUN.fdf > TBT.out sdata siesta.TBT.nc –tbt-av rm siesta.TBT.nc

    after this siesta.TBT.AV.nc exists will all k-averaged quantites. If one is not interested in k-resolved quantities this may be very interesting.

  • Updated the TBT.nc sile for improved readability.

  • Easier script data-extraction from TBT.nc files due to easier conversion between atomic indices and pivoting orbitals.

    For this:

    • a2p returns the pivoting indices for the given atoms (complete set)

    • o2p returns the pivoting indices for the given orbitals

    • Added atom keyword for retrieving DOS for a given set of atoms

    • sdata and TBT.nc files now enable the creation of the TBT.AV.nc file which is the k-averaged file of TBT.nc

  • Faster bond-current algorithms (faster iterator)

  • Initial template for TBT.Proj files for sdata processing

  • Geometry:

    • Enabled multiplying geometries with integers to emulate repeat or tile functions:

      >>> geometry * 2 == geometry.tile(2, 0).tile(2, 1).tile(2, 2)
      >>> geometry * [2, 1, 2] == geometry.tile(2, 0).tile(2, 2)
      >>> geometry * [2, 2] == geometry.tile(2, 2)
      >>> geometry * ([2, 1, 2], 'repeat') == geometry.repeat(2, 0).repeat(2, 2)
      >>> geometry * ([2, 1, 2], 'r') == geometry.repeat(2, 0).repeat(2, 2)
      >>> geometry * ([2, 0], 'r') == geometry.repeat(2, 0)
      >>> geometry * ([2, 2], 'r') == geometry.repeat(2, 2)
      

      This may be considered an advanced feature but useful nonetheless.

    • Enabled “adding” geometries in a similar way as multiplication I.e. the following applies:

      >>> A + B == A.add(B)
      >>> A + (B, 1) == A.append(B, 1)
      >>> A + (B, 2) == A.append(B, 2)
      >>> (A, 1) + B == A.prepend(B, 1)
      
    • Added origo and atom argument to rotation functions. Previously this could be accomblished by:

      rotated = geometry.move(-origo).rotate(…).move(origo)

      while now it is:

      rotated = geometry.rotate(…, origo=origo)

      The origo argument may also be a single integer in which case the rotation is around atom origo.

      Lastly the atom argument enables only rotating a sub-set of atoms.

    • Geometry[..] is now calling axyz if is pure indices, if it is a slice it does not work with super-cell indices

    • Added rij functions to the Geometry for retrieving distances between two atoms (orij for orbitals)

    • Renamed iter_linear to iter

    • Added argument to iter_species for only looping certain atomic indices

    • Added iter_orbitals which returns an iterator with atomic _and_ associated orbitals. The orbitals are with respect to the local orbital indices on the given atom

      ` >>> for ia, io in Geometry.iter_orbitals(): >>>     Geometry.atom[ia].R[io] `

      works, while

      ` >>> for ia, io in Geometry.iter_orbitals(local=False): >>>     Geometry.atom[ia].R[io] `

      does not work because io is globally defined.

    • Changed argument name for coords, atom instead of the old idx.

    • Renamed function axyzsc to axyz

  • SparseCSR:

    • Added iter_nnz(i=None) which loops on sparse elements connecting to row i (or default to loop on all rows and columns).

    • ispmatrix to iterate through a scipy.sparse.*_matrix (and the SparseCSR matrix).

  • Hamiltonian:

    • Added iter_nnz which is the Hamiltonian equivalent of SparseCSR.iter_nnz. It enables explicit looping on atomic couplings, or orbital couplings. I.e. one may specify a subset of atoms or orbitals to loop over.

    • Preliminary implementation of the non-collinear spin-case. Needs testing.