Solving a Graphite Structure Factor Calculation Problem

In summary, the conversation discusses a problem with calculating the structure factor of graphite and finding the correct sum in the structure factor. The questioner has tried multiple calculations and coordinates, but their result does not match the result given in a paper. After further investigation, they realize their result is correct and the paper may have typos. They ask for confirmation from others who have calculated the structure factor of graphite.
  • #1
spf
17
0
Hi,

this is my first post on this forum and I hope that I am in the right category. My question concerns an independent study and not a homework assignment. However, it's a standard homework question in introductory solid state physics courses which makes me both the more frustrated that I cannot figure it out and the more surprised that I cannot find the answer on the web in spite of extensive googling.

Homework Statement



The structure factor of graphite needs to be calculated.

The following article might be helpful to answer this question but is not necessary as I give all the values and formulas below: Chung, JMS 2002 (Journal of Materials Science), p.2, bottom right
http://www.springerlink.com/content/3u4j6xc32h49de1a/

More precisely, the problem lies in the calculation of the following sum in the structure factor

[itex]\sum_{j} e^{- i \vec{G}_m \vec{\rho}_j}[/itex]

with [itex]\vec{G}_m[/itex] being the reciprocal lattice vector and [itex]\vec{\rho}_j[/itex] being the coordinates of the four atoms of the graphite unit cell (assuming the standard ABAB stacking).

The correct result (according to Chung, JMS 2002) should be:

[itex](1 + e^{- i \frac{2\pi}{3}(2m_1-m_2)})(1 + e^{- i \pi m_3})[/itex]

which is the same as

[itex](1 + e^{i \frac{2\pi}{3}(m_1+m_2)})(1 + e^{- i \pi m_3})[/itex]

because

[itex]1 = e^{i 2\pi} = e^{i \frac{2\pi}{3} \cdot 3 m_1}[/itex]

My problem is that I get

[itex](1 + e^{i \frac{2\pi}{3}(m_1+m_2)})(e^{- i \frac{2\pi}{3}(m_1+m_2)} + e^{- i \pi m_3})[/itex]

The result from Chung JMS 2002 seems more reasonable as it explains the experimentally confirmed 001 Bragg spot extinctions which my result does not. But I cannot find any error in my calculations nor in the start values. Can anybody help me out?

Homework Equations



The direct space basis vectors of the graphite lattice are

[itex]\vec{a}_1 = a (\sqrt{3}/2, - 1/2, 0)[/itex]
[itex]\vec{a}_2 = a (\sqrt{3}/2, 1/2, 0)[/itex]
[itex]\vec{a}_3 = c (0, 0, 1)[/itex]

a = 0.246 nm being the graphite in-plane lattice parameter, c = 0.671nm the graphite out-of-plane lattice constant (distance between two A planes in a ABAB stacking) and [itex]a/\sqrt{3}[/itex] being the next-neighbour carbon distance of 0.142nm.

A possible choice of the coordinates of the four atoms of the graphite unit cell (assuming the standard ABAB stacking) are

[itex]\vec{\rho}_A = (0, 0, 0)[/itex]
[itex]\vec{\rho}_B = (\frac{a}{\sqrt{3}}, 0, 0)[/itex]
[itex]\vec{\rho}_A' = (0, 0, \frac{c}{2})[/itex]
[itex]\vec{\rho}_B' = (- \frac{a}{\sqrt{3}}, 0, \frac{c}{2})[/itex]


One gets the reciprocal vectors with the following formula:

[itex]\vec{b}_1 = 2\pi \frac{\vec{a}_2 \times \vec{a}_3}{\vec{a}_1 \cdot (\vec{a}_2 \times \vec{a}_3)}[/itex]
[itex]\vec{b}_2 = 2\pi \frac{\vec{a}_3 \times \vec{a}_1}{\vec{a}_2 \cdot (\vec{a}_3 \times \vec{a}_1)}[/itex]
[itex]\vec{b}_3 = 2\pi \frac{\vec{a}_1 \times \vec{a}_2}{\vec{a}_3 \cdot (\vec{a}_1 \times \vec{a}_2)}[/itex]

The reciprocal lattice vector is then

[itex]\vec{G}_m = m_1 \vec{b}_1 + m_2 \vec{b}_2 + m_3 \vec{b}_3 [/itex]

The sum in the structure factor is then calculated by

[itex]\sum_{j} e^{- i \vec{G}_m \vec{\rho}_j}[/itex]


The Attempt at a Solution



By using the formulas and values given above, one gets

[itex]\vec{b}_1 = \frac{4 \pi}{\sqrt{3} a} (+ 1/2, - \sqrt{3}/2, 0)[/itex]
[itex]\vec{b}_2 = \frac{4 \pi}{\sqrt{3} a} (+ 1/2, + \sqrt{3}/2, 0)[/itex]
[itex]\vec{b}_3 = \frac{2 \pi}{c} (0, 0, 1)[/itex]

and thus

[itex]\vec{G}_m = ( \frac{2 \pi}{\sqrt{3} a} (m_1 + m_2), \frac{2 \pi}{a} (m_2 - m_1), \frac{2 \pi}{c} m_3)[/itex]

and thus the following sum in the structure factor

[itex](1 + e^{i \frac{2\pi}{3}(m_1+m_2)})(e^{- i \frac{2\pi}{3}(m_1+m_2)} + e^{- i \pi m_3})[/itex]

instead of

[itex](1 + e^{i \frac{2\pi}{3}(m_1+m_2)})(1 + e^{- i \pi m_3})[/itex]

Where is the error? Any help would be appreciated.

Thanks,

spf
 
Physics news on Phys.org
  • #2


OK, up to now I tried many things to figure it out, have redone this calculation many times, rechecked the values and have tried out three different sets of coordinates of the A, B, A' and B' atoms in the unit cell and I always get the same result. I went the other way round by asking myself which coordinates the A, B, A' and B' atoms should have within the unit cell in order to produce the result from the paper and I got AA graphite and not the standard ABAB graphite.

In conclusion, I start to think that my result is correct and the one in the paper has a typo. That's not impossible as (despite peer-review) I spotted three other typos in the paper (sign error in the direct space basis vectors (a), sign error in the reciprocal lattice vector (G), missing indices in the text. I was wrong by writing that my result does not explain the 001 extinctions (or more generally the 00l extinctions for l=2n+1). By inserting m1=0 and m2=0 into my result, it predicts these extinctions as required. (By the way, you can predict also these extinctions based on the fact that graphite (the standard ABAB one) has [itex]6_3[/itex] screw symmetry (which means that if you shift it by c/2 in [itex]\vec{a}_3[/itex] direction and rotate it by 60°), it looks exactly the same again.)

So, I guess that's problem solved then. However, there must be lots of people out there who have calculated the structure factor of ABAB graphite in their introductory solid state physics or crystallography classes, either as students or as teachers, so if one of them could just confirm me that my result is correct, it would mean a lot to me. If not, that's OK as well, I am fairly sure of my result by now.
 
Last edited:
  • #3


Dear spf,
I got the same results as yours when I teach the exercise course in solid state physics. your results is correct!
this results can also predict that the (101) diffraction really exist, while in that paper, (101) is extinction, which is totally wrong from experimental results!
 

FAQ: Solving a Graphite Structure Factor Calculation Problem

1. What is a graphite structure factor calculation problem?

A graphite structure factor calculation problem involves determining the intensity of X-rays scattered by a crystal lattice of graphite. This information can be used to understand the arrangement of atoms in the crystal and the distances between them.

2. Why is solving a graphite structure factor calculation problem important?

Solving a graphite structure factor calculation problem is important because it provides valuable information about the structure and properties of graphite. This can be used in materials science, engineering, and other fields to develop new products or improve existing ones.

3. What factors affect the accuracy of a graphite structure factor calculation?

Several factors can affect the accuracy of a graphite structure factor calculation, including the quality of the X-ray source, the crystal structure of the graphite, and the mathematical models used to analyze the data. It is important to carefully control these factors to obtain reliable results.

4. How is a graphite structure factor calculation problem solved?

A graphite structure factor calculation problem is typically solved using X-ray diffraction techniques. This involves directing a beam of X-rays onto the crystal and measuring the intensity of the scattered rays at different angles. This data is then used to calculate the structure factor using mathematical equations.

5. What are the applications of solving a graphite structure factor calculation problem?

The applications of solving a graphite structure factor calculation problem are numerous. In addition to providing information about the structure and properties of graphite, it can also be used in the development of new materials, such as carbon fibers, coatings, and batteries. It is also important in studying the behavior of graphite under different conditions, such as temperature and pressure.

Back
Top