Example: damped second hyperpolarizability: LiH

Download LiH_DampedGamma.run

#! /bin/sh

# If the subkey lifetime and GAMMA or CUBIC is included in the key AORESPONSE,
# the damped (frequency dependent) second hyperpolarizability is calculated.
# This test example consists of three calculations calculations: two with GAMMA
# and one with CUBIC.

# The the first two examples subkey EFIOR is used, which means the electric
# field induced optical rectification gamma(0;omega,-omega,0). In the third
# example two-photon absorption (TPA) cross sections are calculated, the reduced
# form of gamma(-omega;omega,omega,-omega). Examples can easily be modified to
# calculate the static case gamma(0;0,0,0), the optical Kerr effect
# gamma(-omega;omega,0,0), the intensity dependent refractive index
# gamma(-omega;omega,omega,-omega), the electric field induced second harmonic
# generation gamma(-2omega;omega,omega,0), the third harmonic generation
# gamma(-3omega;omega,omega,omega) use the subkey THG, or the general case
# gamma(-(omega1+omega2+omega3);omega1,omega2,omega3).

# Note: results will be physically meaningless due to small basis set. Purpose
# of this job is to provide a test case for the second hyperpolarizability
# implementation

# In the first example the second hyperpolarizability is calculated with the
# subkey GAMMA, for EFIOR, the electric field induced optical rectification
# gamma(0;omega,-omega,0).


####################
# Damped EFIOR, 2n+1
####################

$ADFBIN/adf <<eor
title Damped EFIOR of LiH, 2n+1

basis
 Type SZ
 core None
end

atoms 
Li    0.00000000    0.00000000    3.49467000
H     0.00000000    0.00000000    1.89402300
end

symmetry nosym
allpoints

numericalquality good

aoresponse
 scf iterations 50
 frequency 3 0.1000 -0.1000 0.0000 Hartree
 gamma 
 efior
 ALDA
 lifetime 0.0034
end

eor

rm TAPE21 logfile

# In the second example the second hyperpolarizability is calculated with the
# subkey CUBIC, again for EFIOR, the electric field induced optical
# rectification gamma(0;omega,-omega,0).

#####################
# Damped EFIOR, Cubic
#####################

$ADFBIN/adf <<eor
title Damped EFIOR of LiH, Cubic

basis
 Type SZ
 core None
end

atoms 
Li    0.00000000    0.00000000    3.49467000
H     0.00000000    0.00000000    1.89402300
end

symmetry nosym
allpoints

numericalquality good

aoresponse
 scf iterations 50
 frequency 3 0.1000 -0.1000 0.0000 Hartree
 cubic 
 efior
 ALDA
 lifetime 0.0034
end

eor

rm TAPE21 logfile

# In the third example the subkeyword TPA is be used to calculate the gammagamma
# corresponding to the two photon absorption process (i.e., the reduced form of
# gamma(-omega;omega,omega,-omega)), however, TPA can ONLY be used with keyword
# GAMMA.

##################
# Damped TPA, 2n+1
##################

$ADFBIN/adf <<eor
title Damped TPA of LiH, 2n+1

basis
 Type SZ
 core None
end

atoms 
Li    0.00000000    0.00000000    3.49467000
H     0.00000000    0.00000000    1.89402300
end

symmetry nosym
allpoints

numericalquality good

aoresponse
 scf iterations 50
 frequency 3 0.1000 0.1000 -0.1000 Hartree
 gamma 
 tpa
 ALDA
 lifetime 0.0034
end

eor

rm TAPE21 logfile

# For the static case gamma(0;0,0,0) use the subkey STATIC
#   aoresponse
#    ...
#    frequency 3 0.0000 0.0000 0.0000 Hartree
#    static
#   end

# For the optical Kerr effect gamma(-omega;omega,0,0) use the subkey OKE
#   aoresponse
#    ...
#    frequency 3 0.1000 0.0000 0.0000 Hartree
#    oke
#   end

# For the intensity dependent refractive index gamma(-omega;omega,omega,-omega)
# use the subkey IDRI
#   aoresponse
#    ...
#    frequency 3 0.1000 0.1000 -0.1000 Hartree
#    idri
#   end

# For the electric field induced second harmonic generation
# gamma(-2omega;omega,omega,0) use the subkey EFISHG
#   aoresponse
#    ...
#    frequency 3 0.1000 0.1000 0.0000 Hartree
#    efishg
#   end

# For the third harmonic generation gamma(-3omega;omega,omega,omega) use the
# subkey THG
#   aoresponse
#    ...
#    frequency 3 0.1000 0.1000 0.1000 Hartree
#    thg
#   end

# Or in the general case for gamma(-(omega1+omega2+omega3);omega1,omega2,omega3)
# choose three input frequencies omega1, omega2, and omega3
#   aoresponse
#    ...
#    frequency 3 omega1 omega2 omega3 Hartree
#   end