|
Size: 2824
Comment:
|
Size: 4211
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 7: | Line 7: |
| = resp = | = RESP module for response properties based on HF and DFT = |
| Line 10: | Line 10: |
| {{{ Response properties based on DFT/HF theory. }}} == Keywords == |
== Keywords for general information == |
| Line 22: | Line 18: |
| === CHCK === Check the interface with several external packages. === CTHRD === == Keyworks for processing excited-state information == |
|
| Line 26: | Line 30: |
| === CHCK === | === NFILES === |
| Line 28: | Line 32: |
| Check the interface with several external packages. | Linked with '''istore''' value in TD-DFT input for loading output. |
| Line 30: | Line 34: |
| === NFILES === | == Keyword for geometric derivatives == |
| Line 35: | Line 39: |
| == Keywords for linear response calculations == | |
| Line 37: | Line 42: |
| Enable linear response == REDUCED == Solve the response equation in its reduced form [(A-B)(A+B)-w2](X+Y)=Rvo+Rov (not preferred). |
|
| Line 39: | Line 48: |
| Polarizabiity: '''<<A;B>>(wB)''', where the operators A and B can be dipole (DIP), quadruple (QUA), SOC (HSO), EFG. == Keywords for quadratic response calculations == |
|
| Line 41: | Line 53: |
| Enable quadratic response function (QRF) calculations | |
| Line 43: | Line 56: |
| Hyperpolarizability: <<A;B,C>>(wB,wC) === SINGLE:STATES === Single residue of QRF, STATES can be used to specify the number of states followed by a detailed specification via the triple (ifile,isym,istate). === DOUBLE: PAIRS === Double residue of QRF, PAIRS can be used to specify the number of pairs followed by a detailed specification via two triples (ifile,isym,istate,jfile,isym,jstate). |
|
| Line 45: | Line 65: |
=== SINGLE:STATES === === DOUBLE: PAIRS === |
First-order nonadiabatic couplings |
| Line 51: | Line 68: |
| Neglect the response part of transition density matrix in DOUBLE and FNAC calculations (recommended) == Keywords for finite difference calculations === |
|
| Line 53: | Line 73: |
| Enable finite difference calculations |
|
| Line 54: | Line 76: |
| followed by a real number for the step size, default 0.001 [unit]. |
|
| Line 55: | Line 79: |
| The default unit is angstrom, to use bohr. This keyword must be specified. | |
| Line 56: | Line 81: |
| == REDUCED == | === IGNORE === Ignore the recomputation of excitation energies for check consistency. |
| Line 58: | Line 84: |
| == IGNORE == | |
| Line 60: | Line 85: |
| === CTHRD === == Quick guides by examples == |
= Quick guides by examples = |
| Line 65: | Line 88: |
| === Example: first-order NAC === | == Example: first-order NAC == |
RESP module for response properties based on HF and DFT
Contents
Keywords for general information
IPRT
Print level, >1 gives more information, >2 give more information about integral evaluations.
NPRT
CHCK
Check the interface with several external packages.
CTHRD
Keyworks for processing excited-state information
METHOD
=1, ground state gradients; =2, excited-state calculations which will load TD-DFT output.
NFILES
Linked with istore value in TD-DFT input for loading output.
Keyword for geometric derivatives
GEOM: NORDER
GEOM enables geometric derivatives, NORDER=1, gradient and fo-NACMEs; =2, hessian (not implemented yet.)
Keywords for linear response calculations
LINE
Enable linear response
REDUCED
Solve the response equation in its reduced form [(A-B)(A+B)-w2](X+Y)=Rvo+Rov (not preferred).
POLA: AOPER, BOPER, BFREQ
Polarizabiity: <<A;B>>(wB), where the operators A and B can be dipole (DIP), quadruple (QUA), SOC (HSO), EFG.
Keywords for quadratic response calculations
QUAD
Enable quadratic response function (QRF) calculations
HYPE: AOPER, BOPER, BFREQ, COPER, CFREQ
Hyperpolarizability: <<A;B,C>>(wB,wC)
SINGLE:STATES
Single residue of QRF, STATES can be used to specify the number of states followed by a detailed specification via the triple (ifile,isym,istate).
DOUBLE: PAIRS
Double residue of QRF, PAIRS can be used to specify the number of pairs followed by a detailed specification via two triples (ifile,isym,istate,jfile,isym,jstate).
FNAC
First-order nonadiabatic couplings
NORESP
Neglect the response part of transition density matrix in DOUBLE and FNAC calculations (recommended)
== Keywords for finite difference calculations ===
FDIF
Enable finite difference calculations
STEP
followed by a real number for the step size, default 0.001 [unit].
BOHR
The default unit is angstrom, to use bohr. This keyword must be specified.
IGNORE
Ignore the recomputation of excitation energies for check consistency.
Quick guides by examples
The following examples give the minimal inputs for starting response calculations:
Example: first-order NAC
$COMPASS Title nh3 Basis sto-3g Geometry C 0.00000000 -1.20809142 -1.14173975 C 0.00000000 -1.20797607 0.25342015 C 0.00000000 0.00000000 0.95085852 C -0.00000000 1.20797607 0.25342015 C -0.00000000 1.20809142 -1.14173975 C 0.00000000 0.00000000 -1.83922155 H 0.00000000 -2.16045397 -1.69142002 H 0.00000000 -2.16044427 0.80300713 H -0.00000000 2.16044427 0.80300713 H -0.00000000 2.16045397 -1.69142002 H 0.00000000 0.00000000 -2.93882555 F 0.00000000 0.00000000 2.30085848 End geometry skeleton group c(1) nosym $END $xuanyuan direct schwarz $end $scf RHF charge 0 spin 1 THRESHCONV 1.d-10 1.d-8 OPTSCR 1 iaufbau 0 $end $tddft imethod 1 isf 0 iexit 2 itda 1 idiag 1 istore 1 crit_e 1.d-10 crit_vec 1.d-8 lefteig AOKXC DirectGrid $end $resp iprt 1 QUAD FNAC single states 1 1 1 2 double pairs 1 1 1 1 1 1 2 norder 1 method 2 nfiles 1 FDIF step 0.001 ignore 1 noresp $end
To use finite-difference, a script fdiff.py should be used as
./fbdiff.py run.sh input.inp > log
After the calculation is done, an output file input.out will present in the current directory. The log file saves the information during the calculations.
Note: If FDIF is omitted, the analytic calculation will be carried out by simply using the run.sh script.