PairEx.c File Reference

Calculating the pair excited state generated by the pair operator

\[ \sum_{i,j, \sigma_1, \sigma_2} h_{i,\sigma_1, j, \sigma_2} c_{i\sigma_1} c_{j\sigma_2} (a_{i\sigma_1} a_{j\sigma_2})\]

, where \( c_{i\sigma_1} (a_{i\sigma_1})\) indicates a creation (anti-creation) operator at \(i\)-th site with \( \sigma_1 \) spin. More...

#include "bitcalc.h"
#include "mltplyCommon.h"
#include "PairEx.h"
#include "PairExHubbard.h"
#include "PairExSpin.h"

Go to the source code of this file.

Functions
int	GetPairExcitedState (struct BindStruct *X, double complex *tmp_v0, double complex *tmp_v1)
	Calculating the pair excited state by the pair operator; \[ \sum_{i,j, \sigma_1, \sigma_2} h_{i,\sigma_1, j, \sigma_2} c_{i\sigma_1} c_{j\sigma_2} (a_{i\sigma_1} a_{j\sigma_2})\] , where \( c_{i\sigma_1} (a_{i\sigma_1})\) indicates a creation (anti-creation) operator at \(i\)-th site with \( \sigma_1 \) spin. Target System: Hubbard, Kondo, Spin. More...

Detailed Description

Calculating the pair excited state generated by the pair operator

\[ \sum_{i,j, \sigma_1, \sigma_2} h_{i,\sigma_1, j, \sigma_2} c_{i\sigma_1} c_{j\sigma_2} (a_{i\sigma_1} a_{j\sigma_2})\]

, where \( c_{i\sigma_1} (a_{i\sigma_1})\) indicates a creation (anti-creation) operator at \(i\)-th site with \( \sigma_1 \) spin.

Version

1.2

Author

Kazuyoshi Yoshimi (The University of Tokyo)

Definition in file PairEx.c.

Function Documentation

GetPairExcitedState()

int GetPairExcitedState	(	struct BindStruct *	X,
		double complex *	tmp_v0,
		double complex *	tmp_v1
	)

Calculating the pair excited state by the pair operator;

\[ \sum_{i,j, \sigma_1, \sigma_2} h_{i,\sigma_1, j, \sigma_2} c_{i\sigma_1} c_{j\sigma_2} (a_{i\sigma_1} a_{j\sigma_2})\]

, where \( c_{i\sigma_1} (a_{i\sigma_1})\) indicates a creation (anti-creation) operator at \(i\)-th site with \( \sigma_1 \) spin. Target System: Hubbard, Kondo, Spin.

Parameters

X	[in] define list to get and put information of calculation
tmp_v0	[in, out] Result \( v_0 += H v_1 \).
tmp_v1	[in] The wave vector \( v_1 \) to be mltiplied by the Hamiltonian \( H v_1 \).

Returns

TRUE: Normally finished

FALSE: Unnormally finished

Author

Kazuyoshi Yoshimi (The University of Tokyo)

Version

1.2

Definition at line 47 of file PairEx.c.

References FALSE, GetPairExcitedStateHubbard(), GetPairExcitedStateHubbardGC(), GetPairExcitedStateSpin(), GetPairExcitedStateSpinGC(), GetSplitBitByModel(), GetSplitBitForGeneralSpin(), and X.

Referenced by GetExcitedState().

{
    int iret;
    long unsigned int irght,ilft,ihfbit;

  //  i_max = X->Check.idim_max;
    if(X->Def.iFlgGeneralSpin == FALSE) {
        if (GetSplitBitByModel(X->Def.Nsite, X->Def.iCalcModel, &irght, &ilft, &ihfbit) != 0) {
            return -1;
        }
    }
    else {
        if (GetSplitBitForGeneralSpin(X->Def.Nsite, &ihfbit, X->Def.SiteToBit) != 0) {
            return -1;
        }
    }

  X->Large.i_max    =  X->Check.idim_maxOrg;
  X->Large.irght    = irght;
  X->Large.ilft     = ilft;
  X->Large.ihfbit   = ihfbit;
  X->Large.mode=M_CALCSPEC;

    switch(X->Def.iCalcModel){
  case HubbardGC:
      iret=GetPairExcitedStateHubbardGC(X, tmp_v0, tmp_v1);
    break;

  case KondoGC:
  case Hubbard:
  case Kondo:
      iret=GetPairExcitedStateHubbard(X, tmp_v0, tmp_v1);
    break;

    case Spin: // for the Sz-conserved spin system
      iret =GetPairExcitedStateSpin(X, tmp_v0, tmp_v1);
      break;

    case SpinGC:
      iret=GetPairExcitedStateSpinGC(X,tmp_v0, tmp_v1);
      break;

    default:
        iret =FALSE;
      break;
    }

    return iret;
}