Example: QM/FDE/FQ: excitations of acrolein in water¶
Three examples: a QM/FDE/FQ calculation with the frozen-layer response modeled with FQ, one without the response, and one with no explicit FQ layer.
#!/bin/bash
# This example shows how to perform an excited-state calculation
# of an organic molecule in aqueous solution using the QM/FDE/FQ model
# In this example the response of the FDE layer is also included
# by means of the FQ method
# The calculation starts with the solute embedded within an FQ environment
AMS_JOBNAME="acrolein" "$AMSBIN/ams" <<EOF
System
Atoms
C 26.163000 27.087000 27.854000
O 25.922000 26.598000 28.918000
H 25.425000 27.661000 27.251000
C 27.541000 26.980000 27.309000
C 27.781000 27.369000 26.027000
H 28.266000 26.451000 27.925000
H 28.823000 27.287000 25.727000
H 27.066000 27.797000 25.330000
End
End
Task SinglePoint
Engine ADF
Title Acrolein PBE/DZ single point
Symmetrytolerance 1e-2
Relativity
Level None
End
Basis
Type DZ
Core None
End
XC
GGA PBE
End
ExactDensity
EPrint
SFO NOEIG NOOVL NOORBPOP
SCF NOPOP
End
NoPrint BAS FUNCTIONS
NoSymFit
QMFQ
FORCEFIELD FQ
MolCharge 0.00
AtomType
Symbol O
Chi 0.189194
Eta 0.523700
SubEnd
AtomType
Symbol H
Chi 0.012767
Eta 0.537512
SubEnd
Coords
O 28.826000 25.183000 30.762000 1
H 28.995000 26.089000 31.149000 1
H 28.144000 24.706000 31.315000 1
O 27.763000 31.250000 26.749000 2
H 27.469000 30.950000 25.842000 2
H 26.974000 31.279000 27.363000 2
O 24.563000 24.510000 25.661000 3
H 25.293000 24.179000 25.064000 3
H 23.710000 24.040000 25.436000 3
O 24.219000 23.640000 28.272000 4
H 23.531000 22.942000 28.076000 4
H 24.474000 24.103000 27.423000 4
O 22.979000 24.696000 30.549000 5
H 23.726000 24.648000 29.886000 5
H 23.222000 24.170000 31.364000 5
O 22.709000 27.494000 30.324000 6
H 23.077000 26.584000 30.512000 6
H 22.782000 27.689000 29.346000 6
O 29.009000 27.562000 32.091000 7
H 28.267000 27.944000 31.541000 7
H 29.871000 28.003000 31.840000 7
O 26.815000 24.369000 24.266000 8
H 27.032000 25.081000 23.598000 8
H 27.472000 24.407000 25.018000 8
O 26.727000 23.089000 28.975000 9
H 26.980000 23.432000 29.879000 9
H 25.761000 23.284000 28.803000 9
O 25.883000 31.425000 28.603000 10
H 25.156000 31.260000 29.269000 10
H 26.624000 31.936000 29.038000 10
O 29.317000 29.573000 27.978000 11
H 30.130000 29.270000 27.481000 11
H 28.782000 30.188000 27.398000 11
O 24.110000 29.703000 24.837000 12
H 23.685000 28.874000 25.201000 12
H 23.551000 30.064000 24.091000 12
O 26.708000 30.327000 24.312000 13
H 25.716000 30.396000 24.419000 13
H 26.920000 29.651000 23.606000 13
O 30.565000 24.816000 28.707000 14
H 30.237000 24.502000 27.816000 14
H 29.793000 24.898000 29.337000 14
O 28.442000 23.780000 26.550000 15
H 27.844000 23.801000 27.351000 15
H 28.454000 22.857000 26.166000 15
SubEnd
End
EndEngine
EOF
mv acrolein.results/adf.rkf t21.acrolein
# The QM water layer is then simulated
#H2O B3LYP/DZP
AMS_JOBNAME=H2O "$AMSBIN/ams" <<EOF
System
Atoms
O 23.219000 27.243000 25.752000
H 22.533000 27.083000 25.042000
H 23.738000 26.403000 25.909000
O 27.812000 28.937000 30.049000
H 28.442000 28.831000 29.280000
H 26.915000 28.573000 29.797000
O 23.304000 28.960000 27.797000
H 23.864000 29.646000 27.332000
H 23.228000 28.146000 27.220000
End
End
Task SinglePoint
Engine ADF
Title H2O PBE/DZ single point
Symmetrytolerance 1e-2
Relativity
Level None
End
Basis
Type DZ
Core None
End
XC
GGA PBE
End
ExactDensity
EPrint
SFO NOEIG NOOVL NOORBPOP
SCF NOPOP
End
NoPrint BAS FUNCTIONS
NoSymFit
QMFQ
FORCEFIELD FQ
MolCharge 0.00
AtomType
Symbol O
Chi 0.189194
Eta 0.523700
SubEnd
AtomType
Symbol H
Chi 0.012767
Eta 0.537512
SubEnd
Coords
O 28.826000 25.183000 30.762000 1
H 28.995000 26.089000 31.149000 1
H 28.144000 24.706000 31.315000 1
O 27.763000 31.250000 26.749000 2
H 27.469000 30.950000 25.842000 2
H 26.974000 31.279000 27.363000 2
O 24.563000 24.510000 25.661000 3
H 25.293000 24.179000 25.064000 3
H 23.710000 24.040000 25.436000 3
O 24.219000 23.640000 28.272000 4
H 23.531000 22.942000 28.076000 4
H 24.474000 24.103000 27.423000 4
O 22.979000 24.696000 30.549000 5
H 23.726000 24.648000 29.886000 5
H 23.222000 24.170000 31.364000 5
O 22.709000 27.494000 30.324000 6
H 23.077000 26.584000 30.512000 6
H 22.782000 27.689000 29.346000 6
O 29.009000 27.562000 32.091000 7
H 28.267000 27.944000 31.541000 7
H 29.871000 28.003000 31.840000 7
O 26.815000 24.369000 24.266000 8
H 27.032000 25.081000 23.598000 8
H 27.472000 24.407000 25.018000 8
O 26.727000 23.089000 28.975000 9
H 26.980000 23.432000 29.879000 9
H 25.761000 23.284000 28.803000 9
O 25.883000 31.425000 28.603000 10
H 25.156000 31.260000 29.269000 10
H 26.624000 31.936000 29.038000 10
O 29.317000 29.573000 27.978000 11
H 30.130000 29.270000 27.481000 11
H 28.782000 30.188000 27.398000 11
O 24.110000 29.703000 24.837000 12
H 23.685000 28.874000 25.201000 12
H 23.551000 30.064000 24.091000 12
O 26.708000 30.327000 24.312000 13
H 25.716000 30.396000 24.419000 13
H 26.920000 29.651000 23.606000 13
O 30.565000 24.816000 28.707000 14
H 30.237000 24.502000 27.816000 14
H 29.793000 24.898000 29.337000 14
O 28.442000 23.780000 26.550000 15
H 27.844000 23.801000 27.351000 15
H 28.454000 22.857000 26.166000 15
SubEnd
End
EndEngine
EOF
mv H2O.results/adf.rkf t21.H2O
# Finally, the QM/FDE/FQ calculation is performed
# The response of the FDE layer is activated with the FDERESP keyword
AMS_JOBNAME=FDE "$AMSBIN/ams" <<EOF
System
Atoms
C 26.163000 27.087000 27.854000 adf.f=frag1
O 25.922000 26.598000 28.918000 adf.f=frag1
H 25.425000 27.661000 27.251000 adf.f=frag1
C 27.541000 26.980000 27.309000 adf.f=frag1
C 27.781000 27.369000 26.027000 adf.f=frag1
H 28.266000 26.451000 27.925000 adf.f=frag1
H 28.823000 27.287000 25.727000 adf.f=frag1
H 27.066000 27.797000 25.330000 adf.f=frag1
O 23.219000 27.243000 25.752000 adf.f=frag2 imol=1
H 22.533000 27.083000 25.042000 adf.f=frag2 imol=1
H 23.738000 26.403000 25.909000 adf.f=frag2 imol=1
O 27.812000 28.937000 30.049000 adf.f=frag2 imol=2
H 28.442000 28.831000 29.280000 adf.f=frag2 imol=2
H 26.915000 28.573000 29.797000 adf.f=frag2 imol=2
O 23.304000 28.960000 27.797000 adf.f=frag2 imol=3
H 23.864000 29.646000 27.332000 adf.f=frag2 imol=3
H 23.228000 28.146000 27.220000 adf.f=frag2 imol=3
End
End
Task SinglePoint
Engine ADF
Title Acrolein-H2O PBE/PW91K/DZ FDE single point with interaction energy
SymmetryTolerance 1e-2
STOFit
Relativity
Level None
End
XC
GGA PBE
End
ExactDensity
FDE
PW91K
GGAPOTXFD PBEx
GGAPOTCFD PBEc
FULLGRID
ENERGY
End
Fragments
frag1 t21.acrolein
frag2 t21.H2O type=FDE
End
EPrint
SFO NOEIG NOOVL NOORBPOP
SCF NOPOP
End
NoPrint BAS FUNCTIONS
NoSymFit
Excitations
LOWEST 4
ONLYSING True
End
QMFQ
FORCEFIELD FQ
MolCharge 0.00
FDERESP
AtomType
Symbol O
Chi 0.189194
Eta 0.523700
SubEnd
AtomType
Symbol H
Chi 0.012767
Eta 0.537512
SubEnd
Coords
O 28.826000 25.183000 30.762000 1
H 28.995000 26.089000 31.149000 1
H 28.144000 24.706000 31.315000 1
O 27.763000 31.250000 26.749000 2
H 27.469000 30.950000 25.842000 2
H 26.974000 31.279000 27.363000 2
O 24.563000 24.510000 25.661000 3
H 25.293000 24.179000 25.064000 3
H 23.710000 24.040000 25.436000 3
O 24.219000 23.640000 28.272000 4
H 23.531000 22.942000 28.076000 4
H 24.474000 24.103000 27.423000 4
O 22.979000 24.696000 30.549000 5
H 23.726000 24.648000 29.886000 5
H 23.222000 24.170000 31.364000 5
O 22.709000 27.494000 30.324000 6
H 23.077000 26.584000 30.512000 6
H 22.782000 27.689000 29.346000 6
O 29.009000 27.562000 32.091000 7
H 28.267000 27.944000 31.541000 7
H 29.871000 28.003000 31.840000 7
O 26.815000 24.369000 24.266000 8
H 27.032000 25.081000 23.598000 8
H 27.472000 24.407000 25.018000 8
O 26.727000 23.089000 28.975000 9
H 26.980000 23.432000 29.879000 9
H 25.761000 23.284000 28.803000 9
O 25.883000 31.425000 28.603000 10
H 25.156000 31.260000 29.269000 10
H 26.624000 31.936000 29.038000 10
O 29.317000 29.573000 27.978000 11
H 30.130000 29.270000 27.481000 11
H 28.782000 30.188000 27.398000 11
O 24.110000 29.703000 24.837000 12
H 23.685000 28.874000 25.201000 12
H 23.551000 30.064000 24.091000 12
O 26.708000 30.327000 24.312000 13
H 25.716000 30.396000 24.419000 13
H 26.920000 29.651000 23.606000 13
O 30.565000 24.816000 28.707000 14
H 30.237000 24.502000 27.816000 14
H 29.793000 24.898000 29.337000 14
O 28.442000 23.780000 26.550000 15
H 27.844000 23.801000 27.351000 15
H 28.454000 22.857000 26.166000 15
SubEnd
End
EndEngine
EOF
Download FQQM_FDE_noFDERESP.run
#!/bin/bash
# This example shows how to perform an excited-state calculation
# of an organic molecule in aqueous solution using the QM/FDE/FQ model
# In this example the response of the FDE layer is not included
# The calculation starts with the solute embedded within an FQ environment
AMS_JOBNAME="acrolein" "$AMSBIN/ams" <<EOF
System
Atoms
C 26.163000 27.087000 27.854000
O 25.922000 26.598000 28.918000
H 25.425000 27.661000 27.251000
C 27.541000 26.980000 27.309000
C 27.781000 27.369000 26.027000
H 28.266000 26.451000 27.925000
H 28.823000 27.287000 25.727000
H 27.066000 27.797000 25.330000
End
End
Task SinglePoint
Engine ADF
Title Acrolein PBE/DZ single point
Symmetrytolerance 1e-2
Relativity
Level None
End
Basis
Type DZ
Core None
End
XC
GGA PBE
End
ExactDensity
EPrint
SFO NOEIG NOOVL NOORBPOP
SCF NOPOP
End
NoPrint BAS FUNCTIONS
NoSymFit
QMFQ
FORCEFIELD FQ
MolCharge 0.00
AtomType
Symbol O
Chi 0.189194
Eta 0.523700
SubEnd
AtomType
Symbol H
Chi 0.012767
Eta 0.537512
SubEnd
Coords
O 28.826000 25.183000 30.762000 1
H 28.995000 26.089000 31.149000 1
H 28.144000 24.706000 31.315000 1
O 27.763000 31.250000 26.749000 2
H 27.469000 30.950000 25.842000 2
H 26.974000 31.279000 27.363000 2
O 24.563000 24.510000 25.661000 3
H 25.293000 24.179000 25.064000 3
H 23.710000 24.040000 25.436000 3
O 24.219000 23.640000 28.272000 4
H 23.531000 22.942000 28.076000 4
H 24.474000 24.103000 27.423000 4
O 22.979000 24.696000 30.549000 5
H 23.726000 24.648000 29.886000 5
H 23.222000 24.170000 31.364000 5
O 22.709000 27.494000 30.324000 6
H 23.077000 26.584000 30.512000 6
H 22.782000 27.689000 29.346000 6
O 29.009000 27.562000 32.091000 7
H 28.267000 27.944000 31.541000 7
H 29.871000 28.003000 31.840000 7
O 26.815000 24.369000 24.266000 8
H 27.032000 25.081000 23.598000 8
H 27.472000 24.407000 25.018000 8
O 26.727000 23.089000 28.975000 9
H 26.980000 23.432000 29.879000 9
H 25.761000 23.284000 28.803000 9
O 25.883000 31.425000 28.603000 10
H 25.156000 31.260000 29.269000 10
H 26.624000 31.936000 29.038000 10
O 29.317000 29.573000 27.978000 11
H 30.130000 29.270000 27.481000 11
H 28.782000 30.188000 27.398000 11
O 24.110000 29.703000 24.837000 12
H 23.685000 28.874000 25.201000 12
H 23.551000 30.064000 24.091000 12
O 26.708000 30.327000 24.312000 13
H 25.716000 30.396000 24.419000 13
H 26.920000 29.651000 23.606000 13
O 30.565000 24.816000 28.707000 14
H 30.237000 24.502000 27.816000 14
H 29.793000 24.898000 29.337000 14
O 28.442000 23.780000 26.550000 15
H 27.844000 23.801000 27.351000 15
H 28.454000 22.857000 26.166000 15
SubEnd
End
EndEngine
EOF
mv acrolein.results/adf.rkf t21.acrolein
# The QM water layer is then simulated
#H2O B3LYP/DZP
AMS_JOBNAME=H2O "$AMSBIN/ams" <<EOF
System
Atoms
O 23.219000 27.243000 25.752000
H 22.533000 27.083000 25.042000
H 23.738000 26.403000 25.909000
O 27.812000 28.937000 30.049000
H 28.442000 28.831000 29.280000
H 26.915000 28.573000 29.797000
O 23.304000 28.960000 27.797000
H 23.864000 29.646000 27.332000
H 23.228000 28.146000 27.220000
End
End
Task SinglePoint
Engine ADF
Title H2O PBE/DZ single point
Symmetrytolerance 1e-2
Relativity
Level None
End
Basis
Type DZ
Core None
End
XC
GGA PBE
End
ExactDensity
EPrint
SFO NOEIG NOOVL NOORBPOP
SCF NOPOP
End
NoPrint BAS FUNCTIONS
NoSymFit
QMFQ
FORCEFIELD FQ
MolCharge 0.00
AtomType
Symbol O
Chi 0.189194
Eta 0.523700
SubEnd
AtomType
Symbol H
Chi 0.012767
Eta 0.537512
SubEnd
Coords
O 28.826000 25.183000 30.762000 1
H 28.995000 26.089000 31.149000 1
H 28.144000 24.706000 31.315000 1
O 27.763000 31.250000 26.749000 2
H 27.469000 30.950000 25.842000 2
H 26.974000 31.279000 27.363000 2
O 24.563000 24.510000 25.661000 3
H 25.293000 24.179000 25.064000 3
H 23.710000 24.040000 25.436000 3
O 24.219000 23.640000 28.272000 4
H 23.531000 22.942000 28.076000 4
H 24.474000 24.103000 27.423000 4
O 22.979000 24.696000 30.549000 5
H 23.726000 24.648000 29.886000 5
H 23.222000 24.170000 31.364000 5
O 22.709000 27.494000 30.324000 6
H 23.077000 26.584000 30.512000 6
H 22.782000 27.689000 29.346000 6
O 29.009000 27.562000 32.091000 7
H 28.267000 27.944000 31.541000 7
H 29.871000 28.003000 31.840000 7
O 26.815000 24.369000 24.266000 8
H 27.032000 25.081000 23.598000 8
H 27.472000 24.407000 25.018000 8
O 26.727000 23.089000 28.975000 9
H 26.980000 23.432000 29.879000 9
H 25.761000 23.284000 28.803000 9
O 25.883000 31.425000 28.603000 10
H 25.156000 31.260000 29.269000 10
H 26.624000 31.936000 29.038000 10
O 29.317000 29.573000 27.978000 11
H 30.130000 29.270000 27.481000 11
H 28.782000 30.188000 27.398000 11
O 24.110000 29.703000 24.837000 12
H 23.685000 28.874000 25.201000 12
H 23.551000 30.064000 24.091000 12
O 26.708000 30.327000 24.312000 13
H 25.716000 30.396000 24.419000 13
H 26.920000 29.651000 23.606000 13
O 30.565000 24.816000 28.707000 14
H 30.237000 24.502000 27.816000 14
H 29.793000 24.898000 29.337000 14
O 28.442000 23.780000 26.550000 15
H 27.844000 23.801000 27.351000 15
H 28.454000 22.857000 26.166000 15
SubEnd
End
EndEngine
EOF
mv H2O.results/adf.rkf t21.H2O
# Finally, the QM/FDE/FQ calculation is performed
AMS_JOBNAME=FDE "$AMSBIN/ams" <<EOF
System
Atoms
C 26.163000 27.087000 27.854000 adf.f=frag1
O 25.922000 26.598000 28.918000 adf.f=frag1
H 25.425000 27.661000 27.251000 adf.f=frag1
C 27.541000 26.980000 27.309000 adf.f=frag1
C 27.781000 27.369000 26.027000 adf.f=frag1
H 28.266000 26.451000 27.925000 adf.f=frag1
H 28.823000 27.287000 25.727000 adf.f=frag1
H 27.066000 27.797000 25.330000 adf.f=frag1
O 23.219000 27.243000 25.752000 adf.f=frag2 imol=1
H 22.533000 27.083000 25.042000 adf.f=frag2 imol=1
H 23.738000 26.403000 25.909000 adf.f=frag2 imol=1
O 27.812000 28.937000 30.049000 adf.f=frag2 imol=2
H 28.442000 28.831000 29.280000 adf.f=frag2 imol=2
H 26.915000 28.573000 29.797000 adf.f=frag2 imol=2
O 23.304000 28.960000 27.797000 adf.f=frag2 imol=3
H 23.864000 29.646000 27.332000 adf.f=frag2 imol=3
H 23.228000 28.146000 27.220000 adf.f=frag2 imol=3
End
End
Task SinglePoint
Engine ADF
Title Acrolein-H2O PBE/PW91K/DZ FDE single point with interaction energy
SymmetryTolerance 1e-2
STOFit
Relativity
Level None
End
XC
GGA PBE
End
ExactDensity
FDE
PW91K
GGAPOTXFD PBEx
GGAPOTCFD PBEc
FULLGRID
ENERGY
End
Fragments
frag1 t21.acrolein
frag2 t21.H2O type=FDE
End
EPrint
SFO NOEIG NOOVL NOORBPOP
SCF NOPOP
End
NoPrint BAS FUNCTIONS
NoSymFit
Excitations
LOWEST 4
ONLYSING True
End
QMFQ
FORCEFIELD FQ
MolCharge 0.00
AtomType
Symbol O
Chi 0.189194
Eta 0.523700
SubEnd
AtomType
Symbol H
Chi 0.012767
Eta 0.537512
SubEnd
Coords
O 28.826000 25.183000 30.762000 1
H 28.995000 26.089000 31.149000 1
H 28.144000 24.706000 31.315000 1
O 27.763000 31.250000 26.749000 2
H 27.469000 30.950000 25.842000 2
H 26.974000 31.279000 27.363000 2
O 24.563000 24.510000 25.661000 3
H 25.293000 24.179000 25.064000 3
H 23.710000 24.040000 25.436000 3
O 24.219000 23.640000 28.272000 4
H 23.531000 22.942000 28.076000 4
H 24.474000 24.103000 27.423000 4
O 22.979000 24.696000 30.549000 5
H 23.726000 24.648000 29.886000 5
H 23.222000 24.170000 31.364000 5
O 22.709000 27.494000 30.324000 6
H 23.077000 26.584000 30.512000 6
H 22.782000 27.689000 29.346000 6
O 29.009000 27.562000 32.091000 7
H 28.267000 27.944000 31.541000 7
H 29.871000 28.003000 31.840000 7
O 26.815000 24.369000 24.266000 8
H 27.032000 25.081000 23.598000 8
H 27.472000 24.407000 25.018000 8
O 26.727000 23.089000 28.975000 9
H 26.980000 23.432000 29.879000 9
H 25.761000 23.284000 28.803000 9
O 25.883000 31.425000 28.603000 10
H 25.156000 31.260000 29.269000 10
H 26.624000 31.936000 29.038000 10
O 29.317000 29.573000 27.978000 11
H 30.130000 29.270000 27.481000 11
H 28.782000 30.188000 27.398000 11
O 24.110000 29.703000 24.837000 12
H 23.685000 28.874000 25.201000 12
H 23.551000 30.064000 24.091000 12
O 26.708000 30.327000 24.312000 13
H 25.716000 30.396000 24.419000 13
H 26.920000 29.651000 23.606000 13
O 30.565000 24.816000 28.707000 14
H 30.237000 24.502000 27.816000 14
H 29.793000 24.898000 29.337000 14
O 28.442000 23.780000 26.550000 15
H 27.844000 23.801000 27.351000 15
H 28.454000 22.857000 26.166000 15
SubEnd
End
EndEngine
EOF
#!/bin/bash
# In this example, the FDE method is used to simulate the absorption
# spectrum of acrolein in water.
# Normally, the FDE method neglects to include the response of the frozen
# layer, however this assumption can lead to large systematic errors
# One way to include the response of the environment in a very cost-efficient
# way is to resort to the FQ method, as shown below
AMS_JOBNAME="acrolein" "$AMSBIN/ams" <<EOF
System
Atoms
C 26.163000 27.087000 27.854000
O 25.922000 26.598000 28.918000
H 25.425000 27.661000 27.251000
C 27.541000 26.980000 27.309000
C 27.781000 27.369000 26.027000
H 28.266000 26.451000 27.925000
H 28.823000 27.287000 25.727000
H 27.066000 27.797000 25.330000
End
End
Task SinglePoint
Engine ADF
Title Acrolein PBE/DZ single point
Symmetrytolerance 1e-2
Relativity
Level None
End
Basis
Type DZ
Core None
End
XC
GGA PBE
End
ExactDensity
EPrint
SFO NOEIG NOOVL NOORBPOP
SCF NOPOP
End
NoPrint BAS FUNCTIONS
NoSymFit
EndEngine
EOF
mv acrolein.results/adf.rkf t21.acrolein
#H2O B3LYP/DZP
AMS_JOBNAME=H2O "$AMSBIN/ams" <<EOF
System
Atoms
O 23.219000 27.243000 25.752000
H 22.533000 27.083000 25.042000
H 23.738000 26.403000 25.909000
O 27.812000 28.937000 30.049000
H 28.442000 28.831000 29.280000
H 26.915000 28.573000 29.797000
O 23.304000 28.960000 27.797000
H 23.864000 29.646000 27.332000
H 23.228000 28.146000 27.220000
End
End
Task SinglePoint
Engine ADF
Title H2O PBE/DZ single point
Symmetrytolerance 1e-2
Relativity
Level None
End
Basis
Type DZ
Core None
End
XC
GGA PBE
End
ExactDensity
EPrint
SFO NOEIG NOOVL NOORBPOP
SCF NOPOP
End
NoPrint BAS FUNCTIONS
NoSymFit
EndEngine
EOF
mv H2O.results/adf.rkf t21.H2O
# Having converged the two QM layers the FDE calculation can be carried out
# The response of the frozen layer is modeled using the FQ method, note the use
# of the FDERESP keyword in the QMFQ block
AMS_JOBNAME=FDE "$AMSBIN/ams" <<EOF
System
Atoms
C 26.163000 27.087000 27.854000 adf.f=frag1
O 25.922000 26.598000 28.918000 adf.f=frag1
H 25.425000 27.661000 27.251000 adf.f=frag1
C 27.541000 26.980000 27.309000 adf.f=frag1
C 27.781000 27.369000 26.027000 adf.f=frag1
H 28.266000 26.451000 27.925000 adf.f=frag1
H 28.823000 27.287000 25.727000 adf.f=frag1
H 27.066000 27.797000 25.330000 adf.f=frag1
O 23.219000 27.243000 25.752000 adf.f=frag2 imol=1
H 22.533000 27.083000 25.042000 adf.f=frag2 imol=1
H 23.738000 26.403000 25.909000 adf.f=frag2 imol=1
O 27.812000 28.937000 30.049000 adf.f=frag2 imol=2
H 28.442000 28.831000 29.280000 adf.f=frag2 imol=2
H 26.915000 28.573000 29.797000 adf.f=frag2 imol=2
O 23.304000 28.960000 27.797000 adf.f=frag2 imol=3
H 23.864000 29.646000 27.332000 adf.f=frag2 imol=3
H 23.228000 28.146000 27.220000 adf.f=frag2 imol=3
End
End
Task SinglePoint
Engine ADF
Title Acrolein-H2O PBE/PW91K/DZ FDE single point with interaction energy
SymmetryTolerance 1e-2
STOFit
Relativity
Level None
End
XC
GGA PBE
End
ExactDensity
FDE
PW91K
GGAPOTXFD PBEx
GGAPOTCFD PBEc
FULLGRID
ENERGY
End
Fragments
frag1 t21.acrolein
frag2 t21.H2O type=FDE
End
EPrint
SFO NOEIG NOOVL NOORBPOP
SCF NOPOP
End
NoPrint BAS FUNCTIONS
NoSymFit
Excitations
LOWEST 4
ONLYSING True
End
QMFQ
FORCEFIELD FQ
MolCharge 0.00
FDERESP
AtomType
Symbol O
Chi 0.189194
Eta 0.523700
SubEnd
AtomType
Symbol H
Chi 0.012767
Eta 0.537512
SubEnd
End
EndEngine
EOF