BLW

return to Tutorial n°4

BLW - must know

INPUT (for new file roblw.src)

The BLW input is that of Gamess supplemented by two groups to define the localization of the electrons.

$BLW group

BLW refers to the Block-Localized Wave function (BLW) method, which is a variant of the ab initio Valence Bond (VB) theory. This code can perform RBLW calculations for close-shell systems and ROBLW calculations for open-shell systems. It can also be combined with the grid-DFT calculations. The algorithm to derive self-consistent BLW is based on:

• 1) E. Gianinetti; M. Raimondi; E. Tornaghi, Int. J. Quant. Chem. 60, 157-166(1996).
• 2) A. Famulari, E. Gianinetti, M. Raimondi, M. Sironi, I. Vandoni, Theor. Chem. Acc. 99, 358-365(1998).
• 3) Y. Mo, S. D. Peyerimhoff, J. Chem. Phys. 109, 1687-1697 (1998).
• 4) Y. Mo, J. Gao, S. D. Peyerimhoff, J. Chem. Phys. 112, 5530-5538 (2000).
• 5) L. Song, Y. Lin, Y. Mo, J. Phys. Chem. A, 111, 8291-8301 (2007).

NBLOCK = number of blocks, the definition of each block shall be listed in $BLWDAT group. If nblock=1, it’s just a regular RHF/ROHF calculation.

ITER = maximum number of BLW cycles (default=50)

SCFCOV = NONE, do not take any actions during the SCF iterations DAMP, use damping of the Fock matrix if energy rises DIIS , selects Pulay's DIIS interpolation (based on FOCK matrix) DDIIS, selects Pulay's DIIS interpolation (based on density matrix) (default=DIIS)

IFZB(1) = an array that lists the blocks to be frozen (orbitals kept unchanged) during the optimization. This is useful to examine individual polarization effects, for example.

To switch DFT calculation to full accuracy since the very first step (for the guess orbitals), de-activate coarse grid and SCF pre-optimization. GAMESS always starts an SCF (Hartree-Fock) and a DFT with coarse grid to accelerate the calculation. But if one needs to obtain the polarized energy, the energy of the first iteration shall be exact. To do so, add the following $DFT group. to the .inp file.

 $DFT NRAD0=96 NLEB0=302 NTHE0=12 NPHI0=24  SWOFF=0.0 $END

$BLWDAT group (required by $BLW)

For each block: NE NBF IFLAG

List of basis functions

NE is the number of electrons of that block. NBF is the number of basis functions of that block. IFLAG is an input control. IFLAG = 0, all the basis functions shall be specified; =1, the basis functions are sequential, and only the first one and the last one need be specified; =-1, this block has all the rest basis functions, thus no further specification.

The initial guess can be setup after the definitions of all blocks. It can be useful if the program fails to obtain convergence or one needs to derive the polarized energy term for the BLW-ED (energy decomposition) analysis. Note that the coefficients of MOs correspond to the list of basis functions in the definition lines of blocks. An blank line separates each blocks.

For example:
--------------------------------
$BLW  NBLOCK=2 $END
$BLWDAT
4 4 0       ← block 1 has 4 electrons and 4 basis functions
1 2 3 4    ← bfs of block 1
4 6 1    ← block 2 has 4 electrons and 6 basis functions
5 10     ← bfs of block 2 is from 5 to 10
1.0 0.0 0.0 0.0
0.0 1.0 0.0 0.0                ← initial guess for the 2 occupied orbitals(4/2) of block 1 
                             ← a blanck line at the end of the block
1.0 0.0 0.0 0.0 0.0 0.0   
0.0 1.0 0.0 0.0 0.0 0.0    ← initial guess for the 2 occupied orbitals(4/2) of block 2
$END
---------------------------------