cancel
Showing results for 
Search instead for 
Did you mean: 
cancel
Showing results for 
Search instead for 
Did you mean: 

Community Tip - Did you get an answer that solved your problem? Please mark it as an Accepted Solution so others with the same problem can find the answer easily. X

Fermi integral (01)

payman
14-Alexandrite

Fermi integral (01)

Hello All,

I am trying to calculate the Fermi integral of the order of 1/2 and its derivative. In the attached worksheet, the integral cannot be evaluated for x>0. I used an approximated model as is given in the following article but the result is not similar to what the article obtained.

Koroleva, Olga & V Mazhukin, A & Mazhukin, Vladimir & Breslavskiy, Pavel. (2016). APPROXIMATION OF FERMI-DIRAC INTEGRALS OF DIFFERENT ORDERS USED TO DETERMINE THE THERMAL PROPERTIES OF METALS AND SEMICONDUCTORS. MATHEMATICA MONTISNIGRI. 36. 37-53.  

 

Could you please help me with this problem. 

Thank you very much in advance.

Payman

ACCEPTED SOLUTION

Accepted Solutions
Werner_E
25-Diamond I
(To:payman)

To get single values you may consider using symbolic evaluation with "float"

Bild1.png

For numeric approximation (for plotting) you may also use the same integral with a rather high upper value (not sure in which range this would be an acceptable approximation)

Bild2.png

And to make your plot to math the one you try to duplicate you will have to

  • not use the secondary y-axis (otherwise you will have different y scales)
  • set the plot-limits manually
  • notice, that the plot in the article uses a logarithmic y-scale

Bild3.png

 

 

 

View solution in original post

6 REPLIES 6
Werner_E
25-Diamond I
(To:payman)

To get single values you may consider using symbolic evaluation with "float"

Bild1.png

For numeric approximation (for plotting) you may also use the same integral with a rather high upper value (not sure in which range this would be an acceptable approximation)

Bild2.png

And to make your plot to math the one you try to duplicate you will have to

  • not use the secondary y-axis (otherwise you will have different y scales)
  • set the plot-limits manually
  • notice, that the plot in the article uses a logarithmic y-scale

Bild3.png

 

 

 

Werner_E
25-Diamond I
(To:Werner_E)

This is the second time I noticed it - I am not offered the option to edit my post in the 3 dots menu. Anybody else having the same problem?

 

Just wanted to add, that you don't need to define eta as a vector and create the vectors for FF, etc.. A simple range is all thats neeed and if you are happy with the default range from -10 to 10 of the quickplot feature, you may also omit the range defintion altogether.

Bild4.png

 

payman
14-Alexandrite
(To:Werner_E)

Thank you very much Werner !

I don't know why, but my worksheet is still unable to perform the integration as you showed in your post. Could you please take a look?

Werner_E
25-Diamond I
(To:payman)

I have no idea whats going on. You file failed for me, too. But when I retyped the definition for F10 the very same way you did, it worked - no clue why. I attach the file - hope it works at your end, too.

Pic.PNG

Werner_E
25-Diamond I
(To:Werner_E)

I really wonder what the reason for the error is if we feed the function defined with upper limit infinity with values higher than the one shown in the picture below

Pic2.PNG

payman
14-Alexandrite
(To:Werner_E)

Thank you, Werner, I really appreciate your help. I can now proceed forward.
Announcements

Top Tags