[SOLVED] Are there using coordinate transf in caculating RF?

Phonons, DFPT, electron-phonon, electric-field response, mechanical response…

Moderators: mverstra, joaocarloscabreu

Locked
lzl8181
Posts: 5
Joined: Thu Jan 27, 2011 12:31 am

[SOLVED] Are there using coordinate transf in caculating RF?

Post by lzl8181 » Thu Jan 27, 2011 1:27 am

we know there exist three sets of coordinates in solid physics or crystal physics: the conventional unit (basal vector of this coordinate: A, B,C), the primitive unit (basal vector :a,b,c), and the Cartesian coordinate (basal vector : x,y,z), but in usual caculation, we often adopt the primitive unit (basal vector :a,b,c) or self-defined unit cell (in supercell case) to reduce the weight of caculation, so the confusion comes:
when we set the direction of perturbation (rfdir), which coordinate of the above three it adopts? is that of primitive unit ?
if rfdir=100, which direction it refers to, the direction of vector a, or of vector A, or of vector X? As i know, often the macroscopic symmetry of a crystal is anylyzed in the conventional unit, and many physical properties, such as dielectric function, elastic constants,etc, must correlate with the macroscopic symmetry.
if it's true, does there exist coordinate transformation inside Abinit to get those symmtery-imposed results using the coordinate of primitive unit ? how can this coordinate just give out symmtery-imposed results as we can see from the results of course of elasic tensor caculation?
sorry for my poor english, best wishs!
Last edited by lzl8181 on Thu Jan 27, 2011 9:45 pm, edited 1 time in total.

User avatar
jzwanzig
Posts: 504
Joined: Mon Aug 17, 2009 9:25 am

Re: Are there using coordinate transformation in caculating

Post by jzwanzig » Thu Jan 27, 2011 4:01 pm

rfdir refers to cartesian directions. rprim is the map from cartesian to crystal (abc) directions. The output of the perturbation calculations rfelfd (electric field) and rfstrs (strains, thus elastic constants) are in cartesian directions. How this works will depend a bit on how you define your input structure; for example, an FCC crystal could be orientated such that its conventional axes lie along x/y/z cartesian directions, then rprim will be (0,1/2,1/2),(1/2,0,1/2),(1/2,1/2,0). It could also be imagined to be arbitrarily rotated, and then rprim will reflect this; however, this would be a very hard way to solve the problem. Easiest to imagine your starting crystal as close to an obvious x/y/z orientation as possible.
Josef W. Zwanziger
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com

lzl8181
Posts: 5
Joined: Thu Jan 27, 2011 12:31 am

Re: Are there using coordinate transformation in caculating

Post by lzl8181 » Thu Jan 27, 2011 9:42 pm

Dear prof. Zwanziger,
I'm very glad to receive you timely reply! Thank you for your help! That question has confused me for a long time, many articles give different blurred expressions, it is sloved today!

Locked