QuantumATK Forum
QuantumATK => General Questions and Answers => Topic started by: ramkrishna on April 20, 2011, 09:21
-
Dear Sir ,
I wish to let you know that we are trying to find out the band gap of Hydrogen passivated rectangular Silicon nanowire. In this context there are few suggestions and I am extremely grateful for their suggestions (and the scripts attached).
However, the following scripts does not tend to produce the actual band gap as is generally observed in the literature.
I would be extremely grateful if you correct/modify/suggest me how to make changes in the scripts attached.
Thanking you,
Ramkrishna Ghosh
-
Maybe, most of the band gap calculations of silicon nanowires in literature were done with the first-principles, rather than the Huckel model. As the size of silicon nanowire in your case seems large, you may try the LCAO calculator with the smallest basis set to achieve a higher accuracy than that of the Huckel model.
-
Dear Sir
I have checked the all possible ways (for large as well as for smaller system) to get the accuracy but it does not give any proper result and I am working on that for last few months but whenever I ask in forum, I never get any appropriate answer . Can you please provide me a proper script so that I can work on this?
Regards
Ramkrishna
-
The attached script files are modified for the LCAO calculator with SZP basis set. Please try them.
-
Dear ramkrishna,
First of all, both structures you configured are inadequate.
In the case of SiNW_1.py, you should define the structure as not a bulk but a molecule.
In the case of silicon_wire.py, the length of C vector of the unit cell is wrong.
The correct one is "5.4306" angstrom.
Could you try "silicon_wire-1_modified.py" and verify that the appropriate calculation can be done with it?
-
I've figured it out that the passivation of dangling bonds in the Si nanowire is incomplete.
That's why you can't obtain proper result with "silicon_wire-1_modified2.py".
Sorry for it.
I modified the script again and confirmed that reasonable band structure can be obtained.
Could you try "silicon_wire-1_modified3.py" ?
-
Dear Sir ,
Thank you for your valuable reply . Please let me know the script to do the Hydrogen passivation of dangling bonds so that I can do it properly in lower size like 1X1nm^2 or 2X2nm^2 also . Please also let me know what should I do for changing the bond length along X and Y direction.
Regards
Ramkrishna
-
Dear Sir
I am waiting for your replies. Please help me.
Regards
Ramkrishna
-
Please let me know the script to do the Hydrogen passivation of dangling bonds
I've made the script to do it with Custom Builder but it's at my office now.
My office is closed until 8th May that's why I will able to give it to you next week.
what should I do for changing the bond length along X and Y direction.
Which do you mention, Si-H bond or Si-Si bond?
-
Dear Sir,
I will wait for that and for changing the bond length I want to mention Si-Si bond.
Regards
Ramkrishna
-
for changing the bond length I want to mention Si-Si bond.
Perhaps you can make it by the following steps:
1. drop bare silicon nanowire onto Builder,
2. switch coordinate type to fractional,
3. change the unit cell size against x or y direction.
-
Dear Sir,
I have got the script for the Hydrogen passivation from your post in "Scripts, Tutorials and Application" and applying that I have got a reasonable band structure and band gap value for higher cross section (5X5 nm^2) but for lower cross section (like 3X3 nm^2) it gives almost zero band gap again but it should be around 2eV according to the literature. I can't understand whether there is any limitation of the software for lower cross section or not. Please find the attachment and suggest me how to overcome this problem.
Regards
Ramkrishna
-
In this particular script, the Si atom in the corner with smallest XYZ coordinates is not hydrogenated properly. Maybe there are others too, I just noticed that one since it's easy to spot. This creates a dangling bond which may change the electronic structure, at least a little bit.
There is no limitation in the software that makes certain sizes of system behave differently. However, depending on the computational parameters or model used, the results may of course vary to other studies.
-
I can't check your structure now but it's expected that the bare Si nanowire you prepared would be something wrong.
You should prepare a bare Si nanowire in which all Si-Si bond length are about 2.35angstrom in order to achieve proper sp3 passivation with the script.
-
Dear Sir,
I have made the bare Si nanowire by the following mathod :
custom->nanowire->wire direction(100)->radius 10A-> (100) surface energy 0.1
Now it is difficult to check form the structure that every Si-Si bonds are minimum 2.35 angstrom so that I can get proper sp3 passivation. Please suggest me, how can I achieve this ?
Regards
Ramkrishna
-
Nanowire Wulff constructor has a bug now, that's why you have to modify a unit cell length along z direction before sp3 passivation.
For instance of Si nanowire(100), you should change the length of C vector from "4.70304 angstrom" to "5.4306 angstrom".
-
Dear Sir,
If you have check my script which is attached in the previous post, there you can see that I have done the change of the length of C vector from "4.70304 angstrom" to "5.4306 angstrom" , although I have not got proper result and I also want to verify one thing that, in the nanowire Wulff constructor window, how the radius slot is affecting the length and breadth ? how can I specify the dimension of the rectangular nanowire?
Regards
Ramkrishna
-
If you really do only the step as "custom->nanowire->wire direction(100)->radius 10A-> (100) surface energy 0.1", you should obtain the Si nanowire in which Si-Si bond length is 2.35152 angstrom and only the unit cell size along z direction is wrong.
But in your structure Si-Si bond length shrink to 2.19488 angstrom, that means you did some additional step to get the structure.
I'm not sure what you did, but this is because you obtained the wrong result.
-
Dear Sir,
thank you for your reply and please give the information about my second query in the previous post.
Regards
Ramkrishna
-
I'm afraid that I don't have an answer for the second question because I don't know what algorithm is used for Nanowire Wulff constructor.
-
Dear Sir,
I am getting proper result now. Thank you so much for your kind help.
Regards
Ramkrishna