1
General Questions and Answers / Re: broken bonds in carbon nanotube when optimization
« on: September 16, 2013, 19:21 »
Thanks Dr. Nordland.
Thanks Dr. Blom.
Thanks Dr. Blom.
This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.
Method used? Input file would help.
So in essence you are describing the effect of hot electrons, which usually is very small. In the tutorial I referred to, the main effect is observed when the electrodes is at different temperatures, but also the thermionic effect is larger for larger tunnel distances.
ATK 13.8 can compute phonon frequencies and phonon transport, but electron-phonon interaction is not yet included.
In theory you can do DFT-MD but it's extremely time-consuming, so the sensible choice for mechnical properties is surely classical MD. ATK 13.8 introduces many new functionalities in this area, such as the Tersoff potential, NPT, and a few other things.
If you are assuming Cu-C bonds as covalent ones the optimized bond length 2.0 Ang is reasonable.
My idea to relax the initial structure :
The Cu tube radius of your structure is too big, I mean the strain between CNT and Cu are too tight.
Making a smaller radius of Cu tube, you can construct more relaxed structure as the starting point of the optimization.
And also to enlarge the CNT radius is an alternative idea to release the strain in the system.Hi Mr. Usui, thanks for your answer.
I am not sure about the initial distance. I did use 1.5 Angstrom as the initial one, but after relaxation, for device with less copper atoms, it always arrive at the value of around 2 Angstrom.
I am not sure if that's OK for these device results.
I mean, since it will always go through optimization, finally it would arrive at the equilibrium bond length.
Am I understanding it right?
Secondly, what do you mean by "You should make more relaxed structure as the initial configuration for the optimization"?Please check the covalent radius of Cu and C.
Approximate values of them are
Cu = 1.32 Ang.
C = 0.73 Ang (sp2)
by Wiki.
So the initial bond length of the Cu-C for the optimization calculation should be around 1.32 + 0.73 =~ 2.0 Ang.
Looking your initial configuration for the optimization some Cu-C length are too short ( around 1.46 Ang)
and too much forces could be applied to the system and you found the C-C bond breaking.
You should make more relaxed structure as the initial configuration for the optimization.
That's a bit too big a question to answer for me, so I will not.
But as a general comment for any study, I would say that unless it's performed with the atoms in the equilibrium position, then one should be able to explain why not. And if the equilibrium positions cannot even be determined since the optimization does not work, well, then it's a bit of a problem for sure. I mean, you can dream up any system, but if it's realistic or relevant is a very different story.
I don't see any principal problem with the calculation. Most likely it's simply not a physical situation to have the Cu atoms that close to each other inside the nanotube. The trajectory clearly shows that in the very first steps the Cu atoms experience very large forces, as I interpret it due to the repulsion between Cu-Cu, and so they move outwards which in turn creates large forces on the C atoms, and then it all breaks. So Shinji is correct, the initial structure is the problem, it's not realistic.
Not enough information available. Method used, DFT?
It's probably not really a suitable way to study the temperature effect as the temperature-dependent of the conductance is more due to the interaction of phonons and electrons than the difference in conductance due to some random fluctuation of the atom position with temperature. One could do an MD simulation and take snapshots at different times, do statistics, and the compare the average at different temperatures, but it seems that's not how you do it. In fact, I'm not sure I understand how you take the temperature into account at all in the optimization (are you doing MD?), hence the question about method etc.
Please check the covalent radius of Cu and C.
Approximate values of them are
Cu = 1.32 Ang.
C = 0.73 Ang (sp2)
by Wiki.
So the initial bond length of the Cu-C for the optimization calculation should be around 1.32 + 0.73 =~ 2.0 Ang.
Looking your initial configuration for the optimization some Cu-C length are too short ( around 1.46 Ang)
and too much forces could be applied to the system and you found the C-C bond breaking.
You should make more relaxed structure as the initial configuration for the optimization.
Method used? Input file would help.