void CCE::create_bath_spins()
{
string method = _cfg.getStringParameter("SpinBath", "method");
if( method.compare("TwoDimLattice") == 0 )
{
double lattice_const = _cfg.getDoubleParameter("Lattice", "lattice_const");
string isotope = _cfg.getStringParameter("Lattice", "isotope");
int range_i = _cfg.getIntParameter("Lattice", "full_range");
_lattice = TwoDimFaceCenterLattice(lattice_const, isotope);
_lattice.setRange(range_i);
cSpinSourceFromLattice spin_on_lattice(_lattice);
_bath_spins = cSpinCollection(&spin_on_lattice);
_bath_spins.make();
}
else
{
cSpinSourceFromFile spin_file(_bath_spin_filename);
_bath_spins = cSpinCollection(&spin_file);
_bath_spins.make();
if(_my_rank == 0)
cout << _bath_spins.getSpinNum() << " spins are read from file: " << _bath_spin_filename << endl << endl;
}
}
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;
}/*}}}*/
cSpinSourceFromFile create_spin_bath_from_file(const po::variables_map& para)
{/*{{{*/
string input_filename = INPUT_PATH + para["input"].as<string>();
cSpinSourceFromFile spin_file(input_filename);
return spin_file;
}/*}}}*/
