nextuppreviouscontentsindex
MOLCAS manual:

Next: 8.16 ffpt Up: 8. Programs Previous: 8.14 expbas

Subsections



8.15 falcon


8.15.1 Description

FALCON calculates total energy of the large system based on the fragment approach. Total energy of the whole system is calculated from total energies of fragments as follows,

\begin{displaymath}
E^{whole}=\sum C_i^{fragment} E_i^{fragment},
\end{displaymath}

where Eifragment is the total energy of fragment i, and Cifragment is its coefficient.

In addition to the total energy, FALCON can calculate orbitals of the whole system. Fock matrix and overlap matrix of the whole system are calculated from ones of fragments using following equations,

\begin{displaymath}
F^{whole}=\sum C_i^{fragment} \textbf{F}_i^{fragment},
\end{displaymath}

and

\begin{displaymath}
S^{whole}=\sum C_i^{fragment} \textbf{S}_i^{fragment},
\end{displaymath}

where Fi and Si are the Fock matrix and overlap matrix, respectively, of fragment i.

Then

\begin{displaymath}
\textbf{FC}=\textbf{SC$\varepsilon$}
\end{displaymath}

is solved to obtain the orbitals, C, and orbitals energies, $\varepsilon$.


8.15.2 Input

Below follows a description of the input to FALCON.

The input for each module is preceded by its name like:

  &FALCON

Argument(s) to a keyword, either individual or composed by several entries, can be placed in a separated line or in the same line separated by a semicolon. If in the same line, the first argument requires an equal sign after the name of the keyword.

8.15.2.1 Keywords

KeywordMeaning
TITLeOne-line title.
FRAGmentTakes one, two or three argument(s). The first value (integer) defines the fragment number, the second value (real) determines coefficient, and the third value (integer) is the fragment number that is equivalent to this fragment when translational symmetry is used. A default for the second value is 1.0 where the first and third values have no default. Other keyword(s) specific to this fragment must follow this keyword.
OPERatorA real value following this keyword represents a coefficient, Cifragment, of fragment i (current fragment), where i is a value specified by FRAGMENT keyword. This keyword is equivalent with the second value of keyword, FRAGMENT.
EQUIvalenceAn integer, j, following this keyword declares that current fragment is translationally equivalent with fragment j, and information provided for fragment j are tranfered to current fragment. This keyword is equivalent with the third value of keyword, FRAGMENT.
TRANslateThree real numbers following this keyword specifies the translational vector by which the current fragment is translated to give new coordinate. A unit of either bohr or angstrom can follow. The default unit is angstrom. This keyword takes effect only when the equivalent fragment is specified.
RUNFileFollowing this keyword specifies the name of RunFile file for the corresponding fragment.
ONEIntFollowing this keyword specifies the name of OneInt file for the corresponding fragment.
NFRAgmentAn integer following this keyword specifies the number of fragments. If this keyword is not given, the largest fragment number given by FRAGMENT keyword is set to be the number of fragment.
NIRRepAn integer following this keyword specifies the number of irreducible representation of point group symmetry.
OCCUpationA list of integer(s) following this keyword specifies the number of occupied orbitals in each symmetry representation in the unfragmented system.
DISTanceA real number following this keyword specifies the distance of two atoms that are equivalent to each other, followed by a unit that is eather angstrom or bohr. Default is angstrom.
NEARA real number following this keyword specifies the distance of two atoms within which atoms are considered to be too close each other. An unit that is eather angstrom or bohr can follow. Default is angstrom.
PRINtAn integer following this keyword specifies the format of orbital print out.
ORBEneA real number follwing this keyword stands for the threshold for orbital print out. The orbitals with orbital energy below this value are print out.
ORBOccA real number follwing this keyword stands for the threshold for orbital print out. The orbitals with occupation number above this value are print out.

8.15.2.2 Input examples

Below shows an example of input file for the three fragment system of which energy, Ewhole, is written as

Ewhole= E1fragment + E2fragment - E3fragment,

by fragment energies, E1fragment, E2fragment, and E3fragment.



  &FALCON
  Fragment=1,  1.0
  Fragment=2,  1.0
  Fragment=3,  -1.0

which can be simplified as,



  &FALCON
  Fragment=3,  -1.0

The next example is a two fragment system in which fragment 1 and fragment 2 are equivalent except for their positons. When their difference in position is described by a vector, (1.0, 1.0, -1.0), a translational symmetry can be used and the input becomes as follows,



  &FALCON
  Fragment=2,  1.0,  1
  Translate=1.0,  1.0,  -1.0

If the total energy of the whole system is given by the sum of total energies of three fragment,

Ewhole= E1fragment + E2fragment + E3fragment,

input is simplly as follows,



  &FALCON
  nFragment=3


next up previous contents index
Next: 8.16 ffpt Up: 8. Programs Previous: 8.14 expbas