===== Exploing the chiral field =====

<WRAP box>
<WRAP column 500px>
{{:ncpchem:hlgap.jpg?490}}
</WRAP>
<WRAP column 500px>
The CHFClBr equilibrium geometry, optimised at the scalar relativistic CCSD(T) level (taken from [[http://link.aps.org/doi/10.1103/PhysRevA.65.042508|here]]), was scaled in the range 0.5 - 2.0 and the total and EPV energy
was calculated at the B3LYP level.
As pointed out by Peter the bond stretching will inevitably lead to a breakdown of the SCF model. Indeed one observes that the HOMO-LUMO gap closes
and the SCF did not converge beyond scaling factor 1.4
</WRAP>
</WRAP>

<WRAP box>
<WRAP column 500px>
{{:ncpchem:nucrep.jpg?490}}
</WRAP>
<WRAP column 500px>
The nuclear repulsion energy is evidently not affected by the breakdown of the wavefunction model and drops off as 1/s, where s is the scaling factor.
</WRAP>
</WRAP>

<WRAP box>
<WRAP column 500px>
{{:ncpchem:toterg.jpg?490}}
</WRAP>
<WRAP column 500px>
The total energy can be fitted to a 5th-order polynomial within the range of validity of the model and has, as expected, a minimum near the scaling factor 1.0.
</WRAP>
</WRAP>



<WRAP box>
<WRAP column 500px>
{{:ncpchem:totalepv.jpg?490}}
</WRAP>
<WRAP column 500px>
The total E<sub>PV</sub> energy is not the most interesting quantity since the change in the chiral field for all atoms except the central carbon atom is rather complicated.
</WRAP>
</WRAP>

<WRAP box>
<WRAP column 500px>
{{:ncpchem:epv_c.jpg?490}}
</WRAP>
<WRAP column 500px>
More interesting is the contribution to the E<sub>PV</sub> energy from the central carbon atom. Surprisingly it increases within the range of validity of the model and goes through a minimum around s=0.82.
</WRAP>
</WRAP>