<div dir="ltr">Hi,<div>I'm trying to calculate AlCl dipole moment at equilibrium using finite field by below input:</div><div><br></div><div><br><div><div>***,AlCl finite field calculations</div><div>ANGSTROM</div><div>r=[2.130143]                !set geometry parameters</div><div> geometry={Al;Cl,Al,r}</div><div>basis=aug-cc-pv6z                        !define default basis</div><div>field=[0,0.0003,-0.0003]             !define finite field strengths</div><div>$method=[hf,casscf,ccsd,ccsd(t)]</div><div><br></div><div>k=0</div><div> do i=1,#field                      !loop over fields</div><div> dip,,,field(i)                   !add finite field to H</div><div>do m=1,#method                   !loop over methods</div><div>  k=k+1</div><div>  $method(m)                     !calculate energy</div><div>!expec,dm   !calculate dipole moments for ACPF</div><div>   e(k)=energy                    !save energy </div><div>enddo</div><div>enddo</div><div>k=0</div><div> n=#method</div><div>do m=1,#method</div><div> k=k+1</div><div> energ(m)=e(k)</div><div> dipmz(m)=(e(k+n)-e(k+2*n))/(field(2)-field(3))   !dipole moment as</div><div>first energy derivative</div><div>  dpolz(m)=(e(k+n)+e(k+2*n)-2*e(k))/((field(2)-field(1))*(field(3)-field(1)))</div><div> !polarizability as second der.</div><div>enddo</div><div><br></div><div>table,method,energ,dipmz,dpolz</div><div> title,results for AlCl, r=$R</div></div></div><div><br></div><div><br></div><div>however, the output gives me below error:</div><div><div> SETTING E(6)           =      -701.56161745  HARTREE                         </div><div> SETTING K              =         0.00000000                                  </div><div> SETTING N              =         2.00000000                                  </div><div><br></div><div> DO M                   =         1.00000000   </div><div> SETTING K              =         1.00000000                                  </div><div> SETTING ENERG(1)       =      -701.51812000  HARTREE                         </div><div> SETTING DIPMZ(1)       =        -0.57529364  HARTREE                         </div><div> ? Error</div><div> ? Command FIRST not found</div><div> ? The problem occurs in get_commandset</div></div><div><br></div><div><br></div><div>can anyone please help me to solve this problem</div><div><br></div><div>thank you</div><div><br></div></div>