Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
documentation:language_reference:objects:responsefunction:functions:calculateselfenergy [2024/09/27 11:26] – Aleksandrs Zacinskisdocumentation:language_reference:objects:responsefunction:functions:calculateselfenergy [2024/10/30 18:43] (current) – Aleksandrs Zacinskis
Line 3: Line 3:
  
 ### ###
-alligned paragraph text+ResponseFunction.CalculateSelfEnergy($G_0$,$G$) calculates  
 +\begin{equation}  
 +\Sigma = G_0^{-1} - G^{-1} 
 +\end{equation}
 ### ###
 +===== Input ===== 
 +  * $G_0$ : ResponseFunction object 
 +  * $G$ : ResponseFunction object 
 +  * Possible options are: 
 +      * "epsilon" :  minimal distance between two poles to be considered different in energy (default value 2.3E-13) 
 +===== Output ===== 
 +  * $\Sigma$ : ResponseFunction object
 ===== Example ===== ===== Example =====
- 
 ### ###
-G = GPESDn[1] + ResponseFunction.InvertEnergy(GPESDn[1])### +###
 ==== Input ==== ==== Input ====
-<code Quanty Example.Quanty> +<code Quanty CalculateSelfEnergy.Quanty> 
--- some example code+ 
 +-- H0, H1, Npsi, TPes, TIPes are defined beforehand 
 +psiList = Eigensystem(H0,StartRestrictions, Npsi) 
 + 
 +G0PESIPES_Spectra, G0PESIPES_ResponseFunction = CreateSpectra(H0, {TPes, TIPes}, psiList[1], {{"Emin",-1}, {"Emax",9}, {"NE",1000}, {"Gamma",0.25}}) 
 +G0 = G0PESIPES_ResponseFunction[2]  + ResponseFunction.InvertEnergy(G0PESIPES_ResponseFunction[1]) 
 + 
 +Eigensystem(H0+H1, psiList[1]) -- Use H0 groundstate as Ansatz for Full Hamiltonian groundstate calculation 
 +--For memory efficiency, this way of calling Eigensystem overwrites Ansatz wavefunction with a new one.   
 + 
 +GPESIPES_Spectra, GPESIPES_ResponseFunction = CreateSpectra(H0+H1, {TPes, TIPes}, psiList[1], {{"Emin",-1}, {"Emax",9}, {"NE",1000}, {"Gamma",0.25}}) 
 +G = GPESIPES_ResponseFunction[2]  + ResponseFunction.InvertEnergy(GPESIPES_ResponseFunction[1]) 
 + 
 +Sigma = ResponseFunction.CalculateSelfEnergy(G0,G) 
 + 
 +print(ResponseFunction.ToTable(Sigma))
 </code> </code>
  
 ==== Result ==== ==== Result ====
 <file Quanty_Output> <file Quanty_Output>
-text produced as output +{ { 0.5 , 6.9374060711003e-16 , 3.0413812651491 } ,  
-</file> +  { 0.085601012694643 , 0.16439898730536 } , 
 +  name = Self energy , 
 +  mu = 0 , 
 +  type = ListOfPoles } 
 + </file>
 ===== Table of contents ===== ===== Table of contents =====
 {{indexmenu>.#1|msort}} {{indexmenu>.#1|msort}}
  

You are here: Quanty - a quantum many body script language » Documentation » Language Reference » Objects » Response function » Functions » CalculateSelfEnergy

Print/export