Register to reply 
PhononPhonon Coupling (Different Wave Vectors) 
Share this thread: 
#1
Jul2814, 04:20 AM

P: 5

In a perfect periodic crystal, do lattice vibrational modes corresponding to different wave vectors (kvectors) interact with each other? I'm asking with reference to anharmonic latticedynamics calculations of a timeindependent nature.
For example: consider a linear, monatomic chain with a primitive basis of one atom. Build a supercell containing 3 primitive cells. Calculate lattice vibrations using periodic boundary conditions at the two edges of the supercell (lattice constant is now triple the primitive lattice constant). The primitive cell had 1 degree of freedom, and the supercell will have 3 (31 translation = 2 vibrational dof). One of the supercell's vibrational modes will correspond to wave vector [itex]k=\frac{1}{3}[/itex], and the other to [itex]k=\frac{2}{3}b[/itex], where [itex]b[/itex] is the reciprocal lattice vector for the chain. [itex]k[/itex] and [itex]b[/itex] are defined relative to the primitive cell. Note that in speaking of modes in the supercell which "correspond to" particular wave vectors, I could have more formally said that there is a set of wave vectors which are commensurate with the supercell. Further clarification: if we had a linear, diatomic chain then the 3xsupercell would have 2 vibrational modes each for the wave vectors: [itex]0b, \frac{1}{3}b [/itex] & [itex]\frac{2}{3}b[/itex], and one could imagine modecoupling between the two vibrational modes for [itex]k=\frac{1}{3}b[/itex], say. Does one also include coupling between modes corresponding to different wave vectors? The Question: In a calculation including modecoupling (evaluating the potential energy surface of a nucleus as a function of 2 vibrational mode coordinates [itex]V=V(Q_1,Q_2)[/itex]  [itex]Q_1,Q_2[/itex] are vibrational normal coordinates  rather than just [itex]V=V(Q_1)+V(Q_2)[/itex]), is it correct to couple the modes from different wave vectors? If not, is there a physical/mathematical reason why one refrains from coupling different wave vectors? It seems reasonable to me, as I just imagine coupling modes of different wave vectors as superposing two plane waves. To be clear: I'm trying to avoid a discussion of multiphonon scattering, because I'm only interested in an "equilibrium" snapshot in time / the "mean" vibrational state of the perfect crystal. 


#2
Jul2814, 04:46 AM

Sci Advisor
P: 3,628

As the quadratic part of the Hamiltonian is diagonalized yielding the phonons, the lowest order of phononphonon coupling is trilinear, i.e. it involves 3 phonons. That's an effect of anharmonicity, i.e. also an equilibrium effect.
With a diatomic basis, I don't see why you don't also consider k=0 for the optical phonons. So you could have k1=1/3, k2=2/3 and k3=0 (in units of b). I don't quite understand your mode coupling example. 


#3
Jul2814, 05:36 AM

P: 5

[itex]V(Q_1,Q_2,...,Q_N) \approx \sum_{i=1}^{N}V(Q_i)[/itex] [itex]V(Q_i) := V(Q_i,\left\{Q_l\right\}) , Q_l = 0, \forall l \ne i, l \in \left\{1..N\right\} [/itex] If we consider pairwise coupling, then: [itex]V(Q_1,Q_2,...,Q_N) \approx \sum_{i<j}V(Q_i,Q_j)[/itex] [itex]V(Q_i,Q_j) := V(Q_i,Q_j,\left\{Q_l\right\}) , Q_l = 0, \forall l \ne i,j, l \in \left\{1..N\right\} [/itex] And so on in that fashion for higher order mode couplings. Of course, that's for a molecule. For a supercell in a periodic solid, I'm asking whether or not we should consider mode couplings when the modes correspond to different wave vectors. That is, should we write (for pairwise coupling): [itex]V(Q_1,Q_2,...,Q_N) \approx \sum_{k}\sum_{i^{k}<j^{k}}V(Q_i^{k},Q_j^{k})[/itex] [itex]V(Q_i^{k},Q_j^{k}) := V(Q_i^{k},Q_j^{k},\left\{Q_l\right\}) , Q_l = 0, \forall l \ne i^k,j^k, l \in \left\{1..N\right\} [/itex] ? (where the superscript [itex]k[/itex] indicates that supercell mode corresponds to wave vector [itex]k[/itex]  the sum over [itex]k[/itex] includes all wave vectors commensurate with the supercell) That would be not coupling modes from different wave vectors. If I were to couple different wave vectors, I would retain the molecular formula, only with periodic boundary conditions. I hope that clarifies my question :) If not, what the above equations mean (in words): For no coupling, one constructs a potential energy surface with displacements of nuclei along one vibrational normal coordinate at a time, with all other vibrational coordinates held at their "equilibrium" values. For pairwise coupling, one allows two vibrational coordinates to vary at the same time (take all pairs). In the [itex]N[/itex] mode coupling situation, one would construct the entire potential energy surface by allowing all [itex]N[/itex] coordinates to vary at the same time. Then the approximation would become an equality. 


#4
Jul2914, 01:13 AM

Sci Advisor
P: 3,628

PhononPhonon Coupling (Different Wave Vectors)
This certainly gives you couplings of q's with different k.



#5
Jul2914, 10:16 PM

P: 5

DrDu, that was my first thought as well. To complicate matters, I've found a respectable paper where this appears to be contradicted, and I've been having trouble with my calculations also (results tend to become unphysical when I do couple different kvectors). This is why I'm asking, because there seems to be more to it.
Reference: Physical Review B 78, 224108 (2008) "Effect of Quantization of Vibrations On The Structural Properties of Crystals" Therein, under the heading, "Choice of Vibrational coordinates in the VSCF"  VSCF meaning Vibrational SelfConsistent Field methodology for calculating vibrational frequencies  the following was stated, without justification: "... We now analyze the question of which vibrational coordinates are optimal for the VSCF approximation in the case of the fouratom cell. The three modes are zonecenter optical phonon, which exhibits a double well, and the zoneboundary acoustic and optical phonons, which are practically harmonic. In principle, any linear combination of these modes is possible. Nevertheless, since the modes at different kvectors don't mix, we fix the doublewell mode and optimize the choice of modes in the twodimensional boundary subspace." That was with reference to a 1 dimensional FH (hydrogen bonded) chain, on which VSCF calculations were performed. The "modes" mentioned therein, are represented as eigenvectors in the VSCF; the eigenvectors come from ordinary harmonic approximation (dynamical matrix diagonalisation) calculations. The 4 atom cell is the periodic unit, because the structure is chosen like this: ... F  H     F  H     F  H     F H ... Perhaps they mean something different by the mixing of modes than what I mean by the coupling of eigenvectors, though it sounds exactly the same. Perhaps there is some reason that the different wave vectors don't mix for a 1 dimensional wave that is obvious to everyone else? Or maybe there is actually an important general principle here. Sorry for the longwinded explanations. If you think I should start a new thread, I'd be happy to, but I'm struggling to make this question more concise. Thanks a lot for reading all this. I'm very grateful you've taken the time. 


#6
Jul3014, 01:50 AM

Sci Advisor
P: 3,628

Probably I was wrong or misled by your notation. I would split the index of the Q into two indices where the first one labels the unit cell (or molecule). So let's try to clear this up. We can introduce ##q_{kl}=\sum_j Q_{jl} \exp i j k ## or ##Q_{jl}=1/2 \pi \int dk q_{kl}\exp ikj##. Now your quadratic potential is ##V=\sum_{iljm} V_{iljm} Q_{il} Q_{jm} ## where ## V_{iljm} \propto \delta_{ij}##. Doing the summations you get something like ##V=\sum_{klm}\tilde{V}_{klm}q_{kl}q_{km}=\sum_{klm}\tilde{V}_{klm}q_{kl}q^*_{km}## which is diagonal in k, due to the ##\delta_{ij}##. This can be generalized to potentials which depend only on ji, i.e.translation invariant couplings between different molecules.



#7
Jul3114, 03:36 AM

P: 5

Thanks very much for your reply, which is an introduction to a completely new idea for me. I've done my best to follow it, but I'm not confident I've understood entirely. Here's what I have:
[itex] q_{kl}=\sum_{j}Q_{jl}exp(Ijk) [/itex] [itex] Q_{jl} = \frac{1}{2\pi}\int dkq_{kl}exp(Ikj) [/itex] where [itex]q_{kl}[/itex] is some discrete Fourier transform of [itex]Q_{jl}[/itex], which (in turn), is the eigenvector for vibrational mode [itex]l[/itex] of a primitive cell, broken up into the primitive cell indices [itex]j[/itex]. The index [itex]l[/itex] runs over all vibrational modes in a primitive cell, and the index [itex]j[/itex] runs over all primitive cells in the supercell. [itex]k[/itex] is the wave vector (1 D wavenumber in this case) and [itex]exp(Ijk)[/itex] is just the 1 D version of the phase factor, since [itex]k[/itex] is expressed in fractional coordinates of the reciprocal lattice vectors (reciprocal to the primitive cell real space lattice vectors). I've used capital I as [itex]\sqrt{1}[/itex]. All good so far? You then write the quadratic potential as the summation over [itex]i,j,l,m[/itex] of [itex]V_{ijlm}Q_{il}Q_{jm}[/itex], where [itex]V_{ijlm}[/itex] is ...? Is this a 2nd order partial derivative, or is it [itex]V(Q_{il},Q_{jm})[/itex] with all other [itex]Q_{ab} \ne Q_{il},Q_{jm}[/itex] held at equilibrium? The Kronecker delta: because you've split up the eigenvectors into their primitive cell components, simultaneous displacement along eigenvectors belonging to different primitive cells breaks primitive cell translational invariance (though not supercell translational invariance), and hence the contribution to the potential must be zero for [itex]i \ne j[/itex]. I'd really like clarification here, because I feel I'm starting to stretch a bit. Simplifying gives: [itex]V = \sum_{ilm}V_{ilm}Q_{il}Q_{im} [/itex] where [itex]V_{ilm} := V_{ilim}[/itex] Now jumping to your formula, I can see that (by analogy with earlier discrete Fourier transform): [itex]q_{kl}q_{km} = \sum_{j}Q_{jl}exp(Ijk)Q_{jm}exp(Imk) [/itex] and hence: [itex]Q_{il}Q_{im} = \int dkq_{kl}q_{km}exp(Ik(im)) [/itex] giving: [itex]V = \sum_{ilm}\int dk V_{ilm}q_{kl}q_{km}exp(Ik(im)) [/itex] Since it is only the difference between cells [itex]i[/itex] and [itex]m[/itex] which is important, one can use [itex]m' := im [/itex] and, since the supercell is a finite approximation to the infinite lattice, one must replace the integral over [itex]k[/itex] with a summation over wave vectors commensurate with the supercell. [itex]V = \frac{1}{N_k}\sum_{klm'} \tilde{V}_{klm'} q_{kl}q_{km'} [/itex] with [itex] \tilde{V}_{klm'} := Vexp(ikm') [/itex] and [itex]N_k[/itex] the number of wave vectors commensurate with the supercell. Returning from [itex]m'[/itex] back to [itex]m[/itex] (we remember our reference cell [itex]i[/itex] or consider it shifting by a phase factor), I get something proportional to your final expression: [itex]V = \frac{1}{N_k}\sum_{klm} \tilde{V}_{klm} q_{kl}q_{km}^{*} [/itex] which is indeed diagonal in [itex]k[/itex]. What this means: we have built a paircoupled potential due to the vibrations of different primitive cells in a supercell, without mixing wave vectors, whilst maintaining translational invariance? 


#8
Aug114, 03:33 PM

Sci Advisor
P: 3,628

This sounds all quite reasonable. Take in mind that my post was only a crude sketch, so it is up to you to work out all the 2π's and sums instead of integrals. V_ilmn was meant to be the Taylor coeficient for the second order term.



Register to reply 
Related Discussions  
Phonon Coupling  Mechanical Engineering  0  
Phononphonon interaction in solid lattice  Atomic, Solid State, Comp. Physics  0  
Pseudogap and phonon gap. Does phonon (pseudo)gap exist?  Atomic, Solid State, Comp. Physics  1  
Interaction of a phonon with evanescent wave  Atomic, Solid State, Comp. Physics  0  
Do you know the meaning of the Zero phonon or phononassisted?  Atomic, Solid State, Comp. Physics  2 