Electronic transport properties of molecules, nanowires, and bulks such as superlattice structures can be calculated based on a nonequilibrium Green function (NEGF) method within the collinear and noncollinear DFT methods. The features and capabilities are listed below:
The details of the implementation can be found in Ref. [73]. First the usage of the functionalities for the collinear case is explained in the following subsections. After that, the noncollinear case will be discussed.
System we consider
In the current implementation of OpenMX Ver. 3.9, a system shown in Fig. 39(a) is treated by the NEGF method. The system consists of a central region connected with infinite left and right leads, and the two dimensional periodicity spreads over the bcplane. Considering the two dimensional periodicity, the system can be cast into an onedimensional problem depending on the Bloch wave vector shown in Fig. 39(b). Also, the Green function of the region is selfconsistently determined in order to take account of relaxation of electronic structure around the interface between the central region and the region . It should be noted that the electronic transport is assumed to be along the aaxis in the current implementation. Thus, users have to keep in mind the specification when the geometrical structure is constructed. See also the subsection 'Step 1: The calculations for leads'.

Computational flow
The NEGF calculation is performed by the following three steps:
Step 1 Step 2 Step 3
Each step consists of
The band structure calculations are performed for the left and right leads using a program code 'openmx'. The calculated results will be used to represent the Hamiltonian of the leads in the NEGF calculation of the step 2.
The NEGF calculation is performed for the structure shown in Fig. 39 under zero or a finite bias voltage using a program code 'openmx', where the result in the step 1 is used for the construction of the leads.
By making use of the result of the step 2, the transmission, charge/spin current density, and the eigenchannel are calculated by a program code 'openmx'.
An example: carbon chain
As a first trial, let us illustrate the three steps by employing a carbon chain.
Step 1
%./openmx LeadChain.dat  tee leadchain.stdA file 'negfchain.hks' is generated by the step 1.
Step 2
%./openmx NEGFChain.dat  tee negfchain.stdA file 'negfchain.tranb' is generated by the step 2.
Step 3
openmx
starts the step 3 immediately after it finishes the step 2.
If we perform separately the step 2 and the step 3, we run openmx
as follows:
%./openmx NEGFChain.dat  tee negfchain.std
In the step 3, openmx
reads a file 'negfchain.tranb' and,
calculates the transmission, current, and eigen channels.
After the calculation, the following files:
negfchain.conductance negfchain.tranec0_0_0_2_r.cube negfchain.current negfchain.tranec0_0_0_3_i.cube negfchain.tran0_0 negfchain.tranec0_0_0_3_r.cube negfchain.tranec0_0_0_0_i.cube negfchain.tranec0_0_0_4_i.cube negfchain.tranec0_0_0_0_r.cube negfchain.tranec0_0_0_4_r.cube negfchain.tranec0_0_0_1_i.cube negfchain.traneval0_0_0 negfchain.tranec0_0_0_1_r.cube negfchain.tranevec0_0_0 negfchain.tranec0_0_0_2_i.cubeare generated by the step 3.
The calculations can be traced by using the input files stored in a directory of 'work/negf_example'. By plotting the sixth column in 'negfchain.tran0_0' as a function of the fourth column, you can see a transmission curve as shown in Fig. 40.
