<div dir="ltr">Hi, I am wondering how to calculate octupole moments at the CCSD(T) level of theory via the finite field method? I found an example which explicitly does this for quadrupole moments (snippet below), but could not find anything for octupole moments. Any help would be appreciated, thank you!<div>
<br></div><div><span style="background-color:rgb(255,255,255)"><span style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">k=0</span><br style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">
<span style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">n=#method</span><br style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">
<span style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">do m=1,#method</span><br style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">
<span style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">  k=k+1</span><br style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">
<span style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">  energ(m)=e(k)</span><br style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">
<span style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">  quadzz(m)=(e(k+n)-e(k+2*n))/(field(2)-field(3))  ! quadrupole moment as first energy derivative</span><br style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">
<span style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">  dpolz(m)=(e(k+n)+e(k+2*n)-2*e(k))/((field(2)-field(1))*(field(3)-field(1)))  ! polarizability as second derivative</span><br style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">
<span style="color:rgb(0,0,0);font-family:arial,sans,sans-serif;font-size:13px;font-style:italic;white-space:pre-wrap">enddo</span><br clear="all"></span><div><br></div><div><br></div><div>Best,</div><div>Chris</div><div dir="ltr">
<div><div></div></div></div>
</div></div>