File: tut_ffpt.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 (45 lines) | stat: -rw-r--r-- 1,632 bytes parent folder | download | duplicates (3)
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
.. index::
   single: Program; FFPT
   single: FFPT
   single: Properties; Finite-field PT
   single: Properties; FFPT

.. _TUT\:sec\:ffpt:

:program:`FFPT` --- A Finite Field Perturbation Program
========================================================

Many molecular properties of wave functions can be computed using the
:program:`FFPT` program module in |molcas|. It adds the requested operator to
the integrals computed by the :program:`seward` module. This must be done
before the |molcas| module calculating the required wave function is requested
so the :program:`FFPT` module is best run directly after the :program:`seward`
module.

The :kword:`TITLe` keyword behaviors in a similar fashion to other
|molcas| modules.
The sample input below contains the :program:`FFPT` input
requesting that the dipole moment operator be added to the integrals
using the :kword:`DIPOle` keyword.
The size and direction is specified using the :kword:`COMP` keyword
which accepts free format input. We can compute the dipole of the
molecule by numerical determination of the gradient of the energy
curve determined for several values of the dipole operator. From the second
derivative we can obtain the polarizability component.

Sample input requesting the FFPT module to
include a dipole moment operator in the integral file: ::

  &FFPT
  Title= Finite Perturbation with a dipole in the x negative of strength 0.1 au
  FFPT
  Dipole
   Comp
   X -0.1

:program:`ffpt` Output
----------------------

The :program:`ffpt` section of the output is short and self
explanatory. The :file:`ONEINT` file is updated with the requested
operator.