18.11.2012 Views

ASE Manual Release 3.6.1.2825 CAMd - CampOS Wiki

ASE Manual Release 3.6.1.2825 CAMd - CampOS Wiki

ASE Manual Release 3.6.1.2825 CAMd - CampOS Wiki

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

<strong>ASE</strong> <strong>Manual</strong>, <strong>Release</strong> 3.6.1.2828<br />

nwannier: The number of Wannier functions you wish to construct. This must be at least<br />

half the number of electrons in the system and at most equal to the number of bands in the<br />

calculation.<br />

calc: A converged DFT calculator class. If file arg. is not provided, the calculator must<br />

provide the method get_wannier_localization_matrix, and contain the wavefunctions<br />

(save files with only the density is not enough). If the localization matrix is read<br />

from file, this is not needed, unless get_function or write_cube is called.<br />

Optional arguments:<br />

nbands: Bands to include in localization. The number of bands considered by Wannier can<br />

be smaller than the number of bands in the calculator. This is useful if the highest bands of<br />

the DFT calculation are not well converged.<br />

spin: The spin channel to be considered. The Wannier code treats each spin channel independently.<br />

fixedenergy / fixedstates: Fixed part of Heilbert space. Determine the fixed part of<br />

Hilbert space by either a maximal energy or a number of bands (possibly a list for multiple<br />

k-points). Default is None meaning that the number of fixed states is equated to nwannier.<br />

file: Read localization and rotation matrices from this file.<br />

initialwannier: Initial guess for Wannier rotation matrix. Can be ‘bloch’ to start from<br />

the Bloch states, ‘random’ to be randomized, or a list passed to calc.get_initial_wannier.<br />

seed: Seed for random initialwannier.<br />

verbose: True / False level of verbosity.<br />

get_centers(scaled=False)<br />

Calculate the Wannier centers<br />

pos = L / 2pi * phase(diag(Z))<br />

get_function(index, repeat=None)<br />

Get Wannier function on grid.<br />

Returns an array with the funcion values of the indicated Wannier function on a grid with the size of<br />

the repeated unit cell.<br />

For a calculation using k-points the relevant unit cell for eg. visualization of the Wannier orbitals is<br />

not the original unit cell, but rather a larger unit cell defined by repeating the original unit cell by the<br />

number of k-points in each direction. Note that for a Γ-point calculation the large unit cell coinsides<br />

with the original unit cell. The large unitcell also defines the periodicity of the Wannier orbitals.<br />

index can be either a single WF or a coordinate vector in terms of the WFs.<br />

get_functional_value()<br />

Calculate the value of the spread functional.<br />

Tr[|ZI|^2]=sum(I)sum(n) w_i|Z_(i)_nn|^2,<br />

where w_i are weights.<br />

get_hamiltonian(k=0)<br />

Get Hamiltonian at existing k-vector of index k<br />

dag<br />

H(k) = V diag(eps ) V<br />

k k k<br />

get_hamiltonian_kpoint(kpt_c)<br />

Get Hamiltonian at some new arbitrary k-vector<br />

7.22. Density Functional Theory 159

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!