Skip to Content.
Sympa Menu

forum - Re: [abinit-forum] About the rprim parameter

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: [abinit-forum] About the rprim parameter


Chronological Thread 
  • From: "Anglade Pierre-Matthieu" <anglade@gmail.com>
  • To: forum@abinit.org
  • Subject: Re: [abinit-forum] About the rprim parameter
  • Date: Mon, 30 Apr 2007 19:27:04 +0200
  • Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=beta; h=received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=UT+CoO34Algn3aXu64BIdmZq1+KZm6ni4jkz1dIY9V6B5yvJnrmFVHuTyE1Wb6KmrhY1hH7+w2VCZ8wkJPfjPRk8QJCp6GDasEb67Y2hrKOqt91nMI738T7bsWGHQHSAYk2FXb6v6eGqLyGuUuPOqu9ZSuPu6p6jd7g0UIXILz0=

Dear Chol-Jun Yu
You should probably take the transpose of the matrix you have as input
for rprim.
indeed Abinit is doing the following:
10*(1,0,0)
10*(0,1,0)
14*(1/2,1/2,1/2)
If you input the transpose you'll probably get:
10*(1,0,1/2)
10*(0,1,1/2)
14*(0,0,1/2)
regards

PMA

On 4/30/07, Chol-Jun Yu <yucj@ghi.rwth-aachen.de> wrote:
Dear ABINIT users,

I am confusing about the rprim parameter. If I want to creat the
body-center-tetragonal (bct) primitive unit cell, I will have to make the
primitive lattice vectors R1=(a, 0, 0), R2=(0,a,0), R3=(0.5a, 0.5a, 0.5c). In
order to do that, I can use the rprim as (1, 0, 0), (0, 1, 0), (0.5, 0.5,
0.5), where acells are, for example, 10, 10, 14. At this time it is
expected that R1=(10, 0, 0) R2=(0, 10, 0), R3=(5, 5, 7). However, the results
are different from this expectation in R3 vector: R3=(7, 7, 7). What do I
misunderstand?

Best regards.

Yu Chol Jun
-----------------------------------------------------------------------------
Computational Materials Engineering(CME)
Center for Computational Engineering Science(CCES)
Institute of Mineral Engineering(GHI), RWTH Aachen
Mauerstrasse 5, D-52064 Aachen, Germany
Tel: ++49 241 80 94969
Fax: ++49 241 80 92271
e-mail: yucj@ghi.rwth-aachen.de
-----------------------------------------------------------------------------




--
Pierre-Matthieu Anglade



Archive powered by MHonArc 2.6.16.

Top of Page