From b82fa2fb83c1f84e9e446e1fb498614b62f40a6f Mon Sep 17 00:00:00 2001 From: angelariva Date: Sun, 28 Apr 2024 16:14:31 +0200 Subject: [PATCH] documented fermionic functions --- src/models.jl | 23 ++++++++++++++++++++++- 1 file changed, 22 insertions(+), 1 deletion(-) diff --git a/src/models.jl b/src/models.jl index 6c5ea07..ba42fea 100644 --- a/src/models.jl +++ b/src/models.jl @@ -704,8 +704,20 @@ end """ tightbinding_mpo(N, ϵd, chainparams1, chainparams2) + Generate MPO for a tight-binding chain of N fermionic sites with a single impurity site (fermionic as well) + of energy ϵd at the center. The impurity is coupled to two leads, each described by a set of chain parameters. + The interactions are nearest-neighbour, with the first N/2-1 sites corresponding to the first lead, + the Nth site corresponding to the impurity, and the rest of the sites corresponding to the second + lead. + # Arguments + * `N::Int`: number of sites in the chain + * `ϵd::Real`: energy of the impurity site at the center, as Ed - μ, where μ is the chemical potential + * chainparams1::Array{Real,1}: chain parameters for the first lead + * chainparams2::Array{Real,1}: chain parameters for the second lead + + The chain parameters are given in the standard form: `chainparams` ``=[[ϵ_0,ϵ_1,...],[t_0,t_1,...],c_0]``. """ @@ -816,10 +828,19 @@ end """ interleaved_tightbinding_mpo(N, ϵd, chainparams1, chainparams2) + Generate MPO for a tight-binding chain of N fermionic sites with a single impurity site (fermionic as well) + of energy ϵd. The impurity is coupled to two leads, each described by a set of chain parameters. + The interactions are next-nearest-neighbour, with the first site corresponding to the impurity, and the + two chains organised in an interleaved fashion. + # Arguments + * `N::Int`: number of sites in the chain + * `ϵd::Real`: energy of the impurity site at the first site, as Ed - μ, where μ is the chemical potential + * chainparams1::Array{Real,1}: chain parameters for the first lead + * chainparams2::Array{Real,1}: chain parameters for the second lead - + The chain parameters are given in the standard form: `chainparams` ``=[[ϵ_0,ϵ_1,...],[t_0,t_1,...],c_0]``. """ function interleaved_tightbinding_mpo(N, ϵd, chainparams1, chainparams2)