File: expbas.rst

package info (click to toggle)
openmolcas 25.02-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 170,204 kB
  • sloc: f90: 498,088; fortran: 139,779; python: 13,587; ansic: 5,745; sh: 745; javascript: 660; pascal: 460; perl: 325; makefile: 17
file content (179 lines) | stat: -rw-r--r-- 6,533 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
.. index::
   single: Program; ExpBas
   single: ExpBas

.. _UG\:sec\:expbas:

:program:`expbas`
=================

.. only:: html

  .. contents::
     :local:
     :backlinks: none

.. xmldoc:: <MODULE NAME="EXPBAS">
            %%Description:
            <HELP>
            This module is used to manipulate orbital files. By default it is
            used to expand the orbital file to a larger basis set.
            This feature is available only for generally contracted basis sets (ANO or CC).
            The original orbital file INPORB is needed.
            By default, the Runfile, RUNFIL1,
            from the small calculation is also needed. One runs also the
            GATEWAY module with
            the larger basis set to obtain the corresponding runfile, RUNFIL2. The
            program will produce the file EXPORB, which can be used as the input
            orbital file for calculations using the larger basis set.
            The DESYmmetrize keyword is used to transform orbitals with symmetry to
            C1 symmetry.
            If one wants to desymmetrize the orbitals without expanding the basis set, the keyword
            NOEXpand must be used.
            When the DESY keyword is used, in combination with NOEXpand, only the INPORB file is needed. One runs
            the GATEWAY module
            with the same basis set of the original orbital file to obtain the corresponding
            runfile. The program will produce in output the DESORB file.
            In any subsequent calculation using the DESORB file, modified Cartesian
            coordinates are required, the ones internally generated by Molcas and printed
            out in the GATEWAY section of the output.
            </HELP>

This module is used to manipulate orbital files. By default it is
used to expand the orbital file to a larger basis set.
This feature is available only for generally contracted basis sets (ANO or CC).
The original orbital file :file:`INPORB` is needed.
By default, the Runfile, :file:`RUNFIL1`,
from the small calculation is also needed. One runs also the
:program:`GATEWAY` module with
the larger basis set to obtain the corresponding runfile, :file:`RUNFIL2`. The
program will produce the file :file:`EXPORB`, which can be used as the input
orbital file for calculations using the larger basis set.
The :kword:`DESYmmetrize` keyword is used to transform orbitals with symmetry to
C1 symmetry.
If one wants to desymmetrize the orbitals without expanding the basis set, the keyword
:kword:`NOEXpand` must be used.
When the :kword:`DESYmmetrize`
keyword is used, in combination with :kword:`NOEXpand`, only the :file:`INPORB`
file is needed. One runs the :program:`GATEWAY`
with the same basis set of the original orbital file to obtain the corresponding
runfile. The program will produce the :file:`DESORB` file.
In any subsequent calculation using the `DESORB` file, modified Cartesian
coordinates are required, the ones internally generated by |molcas| and printed
out in the :program:`GATEWAY` section of the output.

.. index::
   pair: Dependencies; ExpBas

.. _UG\:sec\:expbas_dependencies:

Dependencies
------------

The :program:`EXPBAS` requires one InpOrb file and two Runfiles produced by other programs.
If the keyword :kword:`NOEXpand` is used, only the :file:`INPORB` file is needed.

.. index::
   pair: Files; ExpBas

.. _UG\:sec\:expbas_files:

Files
-----

Input files
...........

:file:`INPORB`, an orbital file for the original basis set.

Two more files are needed only if the basis set expansion is to be performed:
:file:`RUNFIL1`, the runfile for the small basis set.
:file:`RUNFIL2`, the runfile for the larger basis set.

Output files
............

:program:`EXPBAS` generates the files:
:file:`EXPORB` which contains orbitals, orbital energies, occupation numbers, and
type indices for the larger basis set.
:file:`DESORB` which contains orbitals, orbital energies, occupation numbers, and
type indices for the desymmetrized basis set.

.. index::
   pair: Input; ExpBas

.. _UG\:sec\:expbas_input:

Input
-----

::

  &EXPBAS

Input example
.............

::

  >>COPY name.InpOrb INPORB
  >>COPY name.RunFil1 RUNFIL1
  >>COPY name.Runfil2 RUNFIL2
   &EXPBAS

Optional general keywords
.........................

.. class:: keywordlist

:kword:`NOEXpand`
  Deactivates the basis set expansion.

  .. xmldoc:: <KEYWORD MODULE="EXPBAS" NAME="NOEX" APPEAR="Do not expand basis set" LEVEL="BASIC" KIND="SINGLE">
              %%Keyword: NOEX <basic>
              <HELP>
              Deactivates the basis set expansion.
              </HELP>
              </KEYWORD>

:kword:`DESYmmetrize`
  Activates the basis set desymmetrization.
  The output is sorted first by orbital kind ascendingly (frozen, inactive, RAS1, ...),
  second by occupation number descendingly (2.0, 2.0, 1.x, 0.0, ...),
  third by energy ascendingly (-3.0, -2.0, -2.0, 0.0, 1.0, ...),
  and fourth by irrep ascendingly (1, 2, ..., 8)

  .. xmldoc:: <KEYWORD MODULE="EXPBAS" NAME="DESY" APPEAR="Desymmetrize basis set" LEVEL="BASIC" KIND="SINGLE">
              %%Keyword: DESY <basic>
              <HELP>
              Activates the basis set desymmetrization.
              </HELP>
              </KEYWORD>

:kword:`FILEorb`
  The next line specifies the filename containing the original orbitals that will
  be expanded to the larger basis set. By default a file named :file:`INPORB` will be used.

  .. xmldoc:: <KEYWORD MODULE="EXPBAS" NAME="FILE" APPEAR="Orbitals file" KIND="STRING" LEVEL="BASIC">
              %%Keyword: FileOrb <basic>
              <HELP>
              The next line specifies the filename containing the original orbitals that will
              be expanded to the larger basis set. By default a file named INPORB will be used.
              </HELP>
              </KEYWORD>

**Comments**. For the basis set expansion
procedure, this is the strategy:

#. Do a calculation (SCF, DFT, CASSCF, etc.) with a small basis set. Save the
   produced orbital file and the runfile as: :file:`name.Inporb` and :file:`name.Runfil1`.

#. Run a :program:`GATEWAY` to create the second runfile, which is then saved as
   :file:`name.Runfil2`.

#. Run :program:`EXPBAS` as indicated above. The file :file:`name.ExpOrb` is produced.
   This file can be used as input for large basis set calculations. Note that the
   inactive and active orbitals will have the same place in the list as with the
   smaller basis set, so they can be easily identified.

.. xmldoc:: </MODULE>