Liouvillian EnsembleCCE::create_spin_liouvillian(const Hamiltonian& hami0, const Hamiltonian hami1)
{
Liouvillian lv0(hami0, SHARP);
Liouvillian lv1(hami1, FLAT);
Liouvillian lv = lv0 - lv1;
return lv;
}
Liouvillian SingleSampleCCE::create_spin_liouvillian(const Hamiltonian& hami0, const Hamiltonian hami1)
{/*{{{*/
Liouvillian lv0(hami0, SHARP);
Liouvillian lv1(hami1, FLAT);
Liouvillian lv = lv0 - lv1;
return lv;
}/*}}}*/
void prepare_data(string filename)
{/*{{{*/
// create defect center
double x = 0.0, y = 0.0, z = 0.89175;
vec coord; coord << x << y << z;
NVCenter nv(NVCenter::N14, coord);
double magBx = 0.0, magBy = 0.0, magBz = 1e-5;
nv.set_magB(magBx, magBy, magBz);
nv.make_espin_hamiltonian();
cout << nv.get_eigen_state(0) << endl;
cout << nv.get_eigen_state(1) << endl;
cSPIN espin=nv.get_espin();
PureState st0(nv.get_eigen_state(0));
PureState st1(nv.get_eigen_state(1));
espin.set_coordinate(coord);
cout << "espin coordinate = " << espin.get_coordinate() << endl;
// create bath spins
cSpinSourceFromFile spin_file(filename);
cSpinCollection spins(&spin_file);
spins.make();
vector<cSPIN> sl = spins.getSpinList();
cout << sl[0].get_coordinate() << sl[0].get_gamma() << endl;
cout << sl[1].get_coordinate() << sl[1].get_gamma() << endl;
vec magB; magB << magBx << magBy << magBz;
SpinZeemanInteraction zee(sl, magB);
SpinDipolarInteraction dip(sl);
DipolarField hf_field0(sl, espin, st0);
DipolarField hf_field1(sl, espin, st1);
Hamiltonian hami0(sl);
hami0.addInteraction(zee);
hami0.addInteraction(dip);
hami0.addInteraction(hf_field0);
hami0.make();
Hamiltonian hami1(sl);
hami1.addInteraction(zee);
hami1.addInteraction(dip);
hami1.addInteraction(hf_field1);
hami1.make();
cout << hami0.getMatrix() << endl;
cout << hami1.getMatrix() << endl;
Liouvillian lv0(hami0, SHARP);
Liouvillian lv1(hami1, FLAT);
Liouvillian lvH = lv0 - lv1;
double rate = 2.0*datum::pi*1e4;
vec axis; axis << 1.0 << 1.0 << 1.0;
SpinDephasing dephasing(sl, rate, normalise(axis));
LiouvilleSpaceOperator dephaseOperator(sl);
dephaseOperator.addInteraction(dephasing);
dephaseOperator.make();
QuantumOperator lv = lvH + dephaseOperator;
lv.saveMatrix("lv");
vec _bath_polarization = zeros<vec>(3);
SpinPolarization p(sl, _bath_polarization);
DensityOperator ds(sl);
ds.addStateComponent(p);
ds.make();
ds.makeVector();
cout << ds.getVector() << endl;
cout << ds.getMatrix() << endl;
}/*}}}*/
