/usr/share/psi/samples/props1/input.dat is in psi4-data 1:0.3-5.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 | #! RHF STO-3G dipole moment computation, performed by applying a finite electric field and numerical differentiation.
memory 250 mb
molecule h2o {
0 1
O
H 1 r
H 1 r 2 a
r=0.98944
a=100.047
}
pert = 0.001
lambdas = [pert, -pert, 2.0*pert, -2.0*pert]
set {
basis sto-3g
e_convergence 10
scf_type pk
}
energies = []
for l in lambdas:
set scf perturb_h true
set scf perturb_with dipole_z
set scf perturb_magnitude $l
energies.append(energy('scf'))
# The nuclear dipole is returned in a zero-based vector, containing [x, y, z]
nuc = nuclear_dipole(h2o)
nuc_z = nuc[2]
# Now use 3- and 5-point finite difference formulae to compute the dipole
dm_z_3point = (energies[0] - energies[1]) / (2.0*pert)
dm_z_5point = (8.0*energies[0] - 8.0*energies[1] - energies[2] + energies[3]) / (12.0*pert)
# The a.u. to Debye conversion factor is automatically available in Psithon as psi_dipmom_au2debye
set scf perturb_h false
# So we can get the analytic result to compare to
energy('scf')
# Tabulate the results of the energy computation
for val in range(len(lambdas)):
print_out("Perturbation strength = %7.4f, computed energy = %16.10f\n" % (lambdas[val], energies[val]))
compare_values(ref_energies[val], energies[val], 8, "Energy for displacement %d" % val)
print_out("nuclear z component = %10.6f ea0, %10.6f Debye\n" % (nuc_z, nuc_z*psi_dipmom_au2debye))
print_out("Electronic contributions:\n")
print_out("3 Point formula: mu(z) = %10.6f ea0, %10.6f Debye\n" % (dm_z_3point, dm_z_3point*psi_dipmom_au2debye))
print_out("5 Point formula: mu(z) = %10.6f ea0, %10.6f Debye\n" % (dm_z_5point, dm_z_5point*psi_dipmom_au2debye))
dm_z_3point = dm_z_3point + nuc_z
dm_z_5point = dm_z_5point + nuc_z
print_out("Total dipoles\n")
print_out("3 Point formula: mu(z) = %10.6f ea0, %10.6f Debye\n" % (dm_z_3point, dm_z_3point*psi_dipmom_au2debye))
print_out("5 Point formula: mu(z) = %10.6f ea0, %10.6f Debye\n" % (dm_z_5point, dm_z_5point*psi_dipmom_au2debye))
compare_values(ref_3pt, dm_z_3point, 8, "Z dipole, using 3 point formula")
compare_values(ref_5pt, dm_z_5point, 8, "Z dipole, using 5 point formula")
|