Momentum density of silicon valence band using bandstructure

  • Context: Graduate 
  • Thread starter Thread starter Samisdabes
  • Start date Start date
  • Tags Tags
    Silicon Valence band
Click For Summary
SUMMARY

The discussion focuses on analyzing the momentum density of the silicon valence band using CASTEP software. The user encountered unexpected results when plotting data extracted from the .bands file into Python. It was clarified that to obtain meaningful band-gradients in momentum space, one should utilize the optical matrix elements by setting the task to 'spectral' and specifying 'spectral_task' as 'optics' in CASTEP. Additionally, ensuring the correct k-points are included in the spectral_kpoints_list is crucial for accurate results.

PREREQUISITES
  • Familiarity with CASTEP software for electronic structure calculations
  • Understanding of band structure and optical properties in solid-state physics
  • Proficiency in Python for data manipulation and visualization
  • Knowledge of k-point sampling in Brillouin zone integration
NEXT STEPS
  • Learn how to configure CASTEP for spectral calculations and optical matrix elements
  • Explore the significance of k-point sampling in band structure analysis
  • Investigate Python libraries for advanced data visualization of band structures
  • Study the relationship between band gradients and optical properties in semiconductors
USEFUL FOR

Researchers and students in solid-state physics, computational chemists, and materials scientists focusing on electronic properties and optical characteristics of materials.

Samisdabes
Messages
1
Reaction score
0
TL;DR
I am attempting to use python to calculate the gradient of the bandstructure for the 3 bands in the valence band for a bandstructure, at a discrete set of points and then plot them in a histogram to get the electron momentum density curve. However I have problems.
I have included the resulting histogram that emerges when I do this procedure and it does not look as expected. I am using CASTEP and have read the .bands file as a txt file into python where i have extracted the relevant data.
Screenshot 2023-12-13 at 20.16.28.png
 
Physics news on Phys.org
Why do you need the gradients? Your plot doesn't really make any sense to me.
 
You don't give us enough information to interpret your plot. However, if you want the band-gradients in momentum space then they are basically the optical matrix elements, so just ask CASTEP for them:

task : spectral
spectral_task : optics

and then make sure your spectral_kpoints_list has the k-points you want.
 
  • Like
Likes   Reactions: Spinnor

Similar threads

  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 4 ·
Replies
4
Views
7K
  • · Replies 8 ·
Replies
8
Views
3K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 3 ·
Replies
3
Views
11K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 8 ·
Replies
8
Views
9K