Setup a Hamiltonian for the Triangular lattice.
(4) Compute the upper limit of the number of Transfer & Interaction and malloc them.
35 int isite, jsite, kCell, ntransMax, nintrMax;
38 double complex Cphase;
44 fp = fopen(
"lattice.gp",
"w");
47 fprintf(stdout,
" @ Lattice Size & Shape\n\n");
61 StdI->
tau[0][0] = 0.0; StdI->
tau[0][1] = 0.0; StdI->
tau[0][2] = 0.0;
65 fprintf(stdout,
"\n @ Hamiltonian \n\n");
70 if (strcmp(StdI->
model,
"spin") == 0 ) {
143 if (strcmp(StdI->
model,
"hubbard") == 0 ) {
153 fprintf(stdout,
"\n @ Numerical conditions\n\n");
159 if (strcmp(StdI->
model,
"kondo") == 0 ) StdI->
nsite *= 2;
162 if (strcmp(StdI->
model,
"spin") == 0 )
164 else if (strcmp(StdI->
model,
"hubbard") == 0 )
167 for (iL = 0; iL < StdI->
nsite / 2; iL++) {
174 if (strcmp(StdI->
model,
"spin") == 0 ) {
175 ntransMax = StdI->
nsite * (StdI->
S2 + 1 + 2 * StdI->
S2);
177 * (3 * StdI->
S2 + 1) * (3 * StdI->
S2 + 1);
180 ntransMax = StdI->
NCell * 2 * (2 * StdI->
NsiteUC + 6 + 6 + 6);
181 nintrMax = StdI->
NCell * (StdI->
NsiteUC + 4 * (3 + 3 + 3));
183 if (strcmp(StdI->
model,
"kondo") == 0) {
184 ntransMax += StdI->
nsite / 2 * (StdI->
S2 + 1 + 2 * StdI->
S2);
185 nintrMax += StdI->
nsite / 2 * (3 * StdI->
S2 + 1) * (3 * StdI->
S2 + 1);
193 for (kCell = 0; kCell < StdI->
NCell; kCell++) {
195 iW = StdI->
Cell[kCell][0];
196 iL = StdI->
Cell[kCell][1];
201 if (strcmp(StdI->
model,
"kondo") == 0 ) isite += StdI->
NCell;
203 if (strcmp(StdI->
model,
"spin") == 0 ) {
209 if (strcmp(StdI->
model,
"kondo") == 0 ) {
218 StdFace_SetLabel(StdI, fp, iW, iL, 1, 0, 0, 0, &isite, &jsite, 1, &Cphase, dR);
220 if (strcmp(StdI->
model,
"spin") == 0 ) {
230 StdFace_SetLabel(StdI, fp, iW, iL, 0, 1, 0, 0, &isite, &jsite, 1, &Cphase, dR);
232 if (strcmp(StdI->
model,
"spin") == 0 ) {
242 StdFace_SetLabel(StdI, fp, iW, iL, 1, - 1, 0, 0, &isite, &jsite, 1, &Cphase, dR);
244 if (strcmp(StdI->
model,
"spin") == 0 ) {
254 StdFace_SetLabel(StdI, fp, iW, iL, 2, - 1, 0, 0, &isite, &jsite, 2, &Cphase, dR);
256 if (strcmp(StdI->
model,
"spin") == 0 ) {
266 StdFace_SetLabel(StdI, fp, iW, iL, 1, 1, 0, 0, &isite, &jsite, 2, &Cphase, dR);
268 if (strcmp(StdI->
model,
"spin") == 0 ) {
278 StdFace_SetLabel(StdI, fp, iW, iL, - 1, 2, 0, 0, &isite, &jsite, 2, &Cphase, dR);
280 if (strcmp(StdI->
model,
"spin") == 0 ) {
290 StdFace_SetLabel(StdI, fp, iW, iL, 2, 0, 0, 0, &isite, &jsite, 3, &Cphase, dR);
292 if (strcmp(StdI->
model,
"spin") == 0) {
302 StdFace_SetLabel(StdI, fp, iW, iL, 0, 2, 0, 0, &isite, &jsite, 3, &Cphase, dR);
304 if (strcmp(StdI->
model,
"spin") == 0) {
314 StdFace_SetLabel(StdI, fp, iW, iL, 2, -2, 0, 0, &isite, &jsite, 3, &Cphase, dR);
316 if (strcmp(StdI->
model,
"spin") == 0) {
325 fprintf(fp,
"plot \'-\' w d lc 7\n0.0 0.0\nend\npause -1\n");
void StdFace_PrintVal_i(char *valname, int *val, int val0)
Print a valiable (integer) read from the input file if it is not specified in the input file (=214748...
double V2
Anisotropic Coulomb potential (1st), input parameter.
double J0pp[3][3]
Isotropic, diagonal/off-diagonal spin coupling (3rd Near.), input parameter J0''x, J0''y, J0''z, J0''xy, etc. or set in StdFace_InputSpin().
double Jp[3][3]
Isotropic, diagonal/off-diagonal spin coupling (2nd Near.), input parameter J'x, J'y, J'z, J'xy, etc.
double complex t2p
Anisotropic hopping (2nd), input parameter.
double J[3][3]
Isotropic, diagonal/off-diagonal spin coupling (1st Near.), input parameter Jx, Jy, Jz, Jxy, etc.
void StdFace_GeneralJ(struct StdIntList *StdI, double J[3][3], int Si2, int Sj2, int isite, int jsite)
Treat J as a 3*3 matrix [(6S + 1)*(6S' + 1) interactions].
void StdFace_HubbardLocal(struct StdIntList *StdI, double mu0, double h0, double Gamma0, double U0, int isite)
Add intra-Coulomb, magnetic field, chemical potential for the itenerant electron. ...
double D[3][3]
Coefficient for input parameter D. Only D[2][2] is used.
double J1p[3][3]
Isotropic, diagonal/off-diagonal spin coupling (2nd Near.), input parameter J1'x, J1'y...
void StdFace_PrintGeometry(struct StdIntList *StdI)
Print geometry of sites for the pos-process of correlation function.
double J1ppAll
Anisotropic, diagonal spin coupling (3rd Near), input parameter J1''.
double J2pp[3][3]
Isotropic, diagonal/off-diagonal spin coupling (3rd Near.), input parameter J2''x, J2''y, J2''z, J2''xy, etc. or set in StdFace_InputSpin().
double JpAll
Isotropic, diagonal spin coupling (2nd Near), input parameter Jp.
void StdFace_InputHopp(double complex t, double complex *t0, char *t0name)
Input hopping integral from the input file, if it is not specified, use the default value(0 or the is...
double J1[3][3]
Isotropic, diagonal/off-diagonal spin coupling (1st Near.), input parameter J1x, J1y, J1z, J1xy, etc. or set in StdFace_InputSpinNN().
double J2p[3][3]
Isotropic, diagonal/off-diagonal spin coupling (2nd Near.), input parameter J2'x, J2'y...
void StdFace_Hopping(struct StdIntList *StdI, double complex trans0, int isite, int jsite, double *dR)
Add Hopping for the both spin.
double complex t
Nearest-neighbor hopping, input parameter.
void StdFace_MallocInteractions(struct StdIntList *StdI, int ntransMax, int nintrMax)
Malloc Arrays for interactions.
double JAll
Isotropic, diagonal spin coupling (1st Near.), input parameter J.
int S2
Total spin |S| of a local spin, input from file.
int NsiteUC
Number of sites in the unit cell. Defined in the beginning of each lattice function.
void StdFace_InputSpin(double Jp[3][3], double JpAll, char *Jpname)
Input spin-spin interaction other than nearest-neighbor.
double J1All
Anisotropic, diagonal spin coupling (1st Near), input parameter J1.
double V0p
Anisotropic Coulomb potential (2nd), input parameter.
double J1pp[3][3]
Isotropic, diagonal/off-diagonal spin coupling (3rd Near.), input parameter J1''x, J1''y, J1''z, J1''xy, etc. or set in StdFace_InputSpin().
void StdFace_InitSite(struct StdIntList *StdI, FILE *fp, int dim)
Initialize the super-cell where simulation is performed.
double V2p
Anisotropic Coulomb potential (2nd), input parameter.
double V2pp
Anisotropic Coulomb potential (3rd), input parameter.
double J2ppAll
Anisotropic, diagonal spin coupling (3rd Near), input parameter J2''.
double complex tpp
3rd-nearest hopping, input parameter
char model[256]
Name of model, input parameter.
void StdFace_NotUsed_J(char *valname, double JAll, double J[3][3])
Stop HPhi if variables (real) not used is specified in the input file (!=NaN).
double Gamma
Transvars magnetic field, input parameter.
double JppAll
Isotropic, diagonal spin coupling (3rd Near), input parameter J''.
double complex t1pp
Anisotropic hopping (3rd), input parameter.
double V1
Anisotropic Coulomb potential (1st), input parameter.
double J0All
Anisotropic, diagonal spin coupling (1st Near), input parameter J0.
double J0[3][3]
Isotropic, diagonal/off-diagonal spin coupling (1st Near.), input parameter J0x, J0y, J0z, J0xy, etc. or set in StdFace_InputSpinNN().
double J0ppAll
Anisotropic, diagonal spin coupling (3rd Near), input parameter J0''.
double U
On-site Coulomb potential, input parameter.
double complex t0p
Anisotropic hopping (2nd), input parameter.
int ** Cell
[StdIntList][3] The cell position in the fractional coordinate. Malloc and Set in StdFace_InitSite()...
double complex t2pp
Anisotropic hopping (3rd), input parameter.
double phase[3]
Boundary phase, input parameter phase0, etc.
double V1pp
Anisotropic Coulomb potential (3rd), input parameter.
double J0pAll
Anisotropic, diagonal spin coupling (2nd Near), input parameter J0'.
double length[3]
Anisotropic lattice constant, input parameter wlength, llength, hlength.
int * locspinflag
[StdIntList::nsite] LocSpin in Expert mode, malloc and set in each lattice file.
double complex tp
2nd-nearest hopping, input parameter
double Jpp[3][3]
Isotropic, diagonal/off-diagonal spin coupling (3rd Near.), input parameter J''x, J''y...
double complex t1
Anisotropic hopping (1st), input parameter.
double V
Off-site Coulomb potential (1st), input parameter.
double complex t0
Anisotropic hopping (1st), input parameter.
double direct[3][3]
The unit direct lattice vector. Set in StdFace_InitSite().
void StdFace_MagField(struct StdIntList *StdI, int S2, double h, double Gamma, int isite)
Add longitudinal and transvars magnetic field to the list.
void StdFace_NotUsed_d(char *valname, double val)
Stop HPhi if a variable (real) not used is specified in the input file (!=NaN).
double complex t0pp
Anisotropic hopping (3rd), input parameter.
double V0
Anisotropic Coulomb potential (1st), input parameter.
void StdFace_Coulomb(struct StdIntList *StdI, double V, int isite, int jsite)
Add onsite/offsite Coulomb term to the list StdIntList::Cinter and StdIntList::CinterIndx, and increase the number of them (StdIntList::NCinter).
void StdFace_InputCoulombV(double V, double *V0, char *V0name)
Input off-site Coulomb interaction from the input file, if it is not specified, use the default value...
void StdFace_SetLabel(struct StdIntList *StdI, FILE *fp, int iW, int iL, int diW, int diL, int isiteUC, int jsiteUC, int *isite, int *jsite, int connect, double complex *Cphase, double *dR)
Set Label in the gnuplot display (Only used in 2D system)
double Vpp
Off-site Coulomb potential (3rd), input parameter.
double complex t1p
Anisotropic hopping (2nd), input parameter.
double complex t2
Anisotropic hopping (1st), input parameter.
double mu
Chemical potential, input parameter.
void StdFace_PrintVal_d(char *valname, double *val, double val0)
Print a valiable (real) read from the input file if it is not specified in the input file (=NaN)...
int NCell
The number of the unit cell in the super-cell (determinant of StdIntList::box). Set in StdFace_InitSi...
double Vp
Off-site Coulomb potential (2nd), input parameter.
double J1pAll
Anisotropic, diagonal spin coupling (2nd Near), input parameter J1'.
void StdFace_NotUsed_c(char *valname, double complex val)
Stop HPhi if a variable (complex) not used is specified in the input file (!=NaN).
int nsite
Number of sites, set in the each lattice file.
double J2All
Anisotropic, diagonal spin coupling (1st Near), input parameter J2.
double V0pp
Anisotropic Coulomb potential (3rd), input parameter.
void StdFace_NotUsed_i(char *valname, int val)
Stop HPhi if a variable (integer) not used is specified in the input file (!=2147483647, the upper limt of Int).
double h
Longitudinal magnetic field, input parameter.
double a
The lattice constant. Input parameter.
double J2[3][3]
Isotropic, diagonal/off-diagonal spin coupling (1st Near.), input parameter J2x, J2y, J2z, J2xy, etc. or set in StdFace_InputSpinNN().
double J2pAll
Anisotropic, diagonal spin coupling (2nd Near), input parameter J2'.
double ** tau
Cell-internal site position in the fractional coordinate. Defined in the beginning of each lattice fu...
double J0p[3][3]
Isotropic, diagonal/off-diagonal spin coupling (2nd Near.), input parameter J0'x, J0'y...
void StdFace_InputSpinNN(double J[3][3], double JAll, double J0[3][3], double J0All, char *J0name)
Input nearest-neighbor spin-spin interaction.
double V1p
Anisotropic Coulomb potential (2nd), input parameter.
double K
4-spin term. Not used.