1. Not finding help here? Sign up for a free 30min tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Vibrational frequency of C2

  1. Dec 7, 2009 #1
    1. The problem statement, all variables and given/known data
    The proposed problem is to find the vibrational frequency of a carbon dimer C2. Then I have to write and run an MD simulation to find the period of oscillations and compare the two.

    This is part of a final project for a computational/numerical methods course I am taking. Up until now, we have done nothing but numerical computations, but the professor decided it would be good for us to know what it is like to solve a real physics problem. The only prerequisite for this course is experience in programming, so I am not very far along in my physics (through Modern I: intro to quantum and classical dynamics). I am a bit overwhelmed, so any help is greatly appreciated.

    2. Relevant equations
    [tex]E_{tot}=E_{bs}+E_{rep}[/tex]

    where [tex]E_{bs}[/tex] is the sum of electronic eigenvalues over all occupied states, and
    [tex]E_{rep}[/tex] is a short-ranged repulsive energy.

    [tex]E_{rep}=\sum_{i}f\left(\sum{j}\phi\left(r_{ij}\right)\right)[/tex]

    where [tex]\phi\left(r_{ij}\right)[/tex] is a pairwise potential between atoms i and j, and f is a functional expressed as a 4th-order polynomial with argument [tex]\sum{j}\phi\left(r_{ij}\right)[/tex].

    [tex]s\left(r\right)=\left(r_{0}/r\right)^{n}exp\left(n\left[-\left(r/r_{c}\right)^{n_{c}}+\left(r_{0}/r_{c}\right)^{n_{c}}\right] \right)[/tex]

    [tex]\phi\left(r\right)=\phi_{0}\left(d_{0}/r\right)^{m}exp\left(m\left[-\left(r/d_{c}\right)^{m_{c}}+\left(d_{0}/d_{c}\right)^{m_{c}}\right] \right)[/tex]

    where [tex]r_{0}[/tex] denotes the nearest-neighbor atomic separations, and [tex]n, n_{c}, r_{c}, \phi_{0}, m, d_{c}, and m_{c}[/tex] are parameters that need to be determined.

    These equations are from A transferable tight-binding potential for carbon by C H Xu et al. They describe the process of finding the total energy for diamond.


    3. The attempt at a solution
    In order to calculate the vibrational frequency, I need to calculate the total energy of the molecule. Once I have the total energy, I can plot the energy as a function of the inter-atomic distance and take the second derivative to find the spring constant and calculate the vibrational frequency.

    For Ebs, the electronic eigenvalues can be obtained by solving an empirical tight-binding Hamiltonian [tex]H_{TB}[/tex]. The off-diagonal are described by a set of orthogonal sp3 two-center hopping parameters, [tex]V_{ss\sigma}[/tex],
    [tex]V_{sp\sigma}[/tex], [tex]V_{pp\sigma}[/tex], and [tex]V_{pp\pi}[/tex], scaled with interatomic separation r as a function of s(r); and the on-site elements are the atomic orbital energies of the corresponding atom.

    This is all I know. I know what is necessary to get Ebs, but I don't know how to construct the matrix. I don't know what dimensions it should be, but I suspect 8x8 because carbon is tetravalent, and I don't really understand what the elements should be. Are the diagonal elements found from the Hamiltonian? What are the hopping parameters for the off-diagonal elements?

    If anyone has a good article or webpage that might help me, it would also be appreciated.
     
  2. jcsd
  3. Dec 8, 2009 #2
    sorry, I am still having issues with this, so i am bumping the thread to make it more visible
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook