add a member function Op *=cd
switch to Hz in function omega
add a notice for CSA constant
nothing new except some code optimization
modify the HZ ham,since the omega function in Hz
nothing