[Pw_forum] input for a polymer crystal
Axel Kohlmeyer
akohlmey at cmm.chem.upenn.edu
Sun Jan 11 05:12:49 CET 2009
On Sat, Jan 10, 2009 at 10:12 PM, Jonas Baltrusaitis <jasius_1 at yahoo.com> wrote:
> I am trying to optimize a hydrogen bonded polymer crystal structure. I have experimental a,b,c and all three angles of the cell, as well as the point group. Also fractional coordinates
>
> what I don't have is an idea how to prepare an input. I read Doc/Input_pw closely and couldn't find how to utilize any of those (except atomic positions).
>
> CELL_PARAMETERS should accommodate those somwehow but it's only 3x3 matrix input in the manual instead. also, specifying ibrav, celldm, nat, ntyp and nbnd are not clear as it's nothing like a simple Al or Si structure. also, what do I put for K_POINTS? Gamma only?
if you have all you this information, ibrav, celldm, nat, ntyp are all defined.
CELL_PARAMETERS is only needed if you set ibrav=0. nbnd depends on
your pseudopotentials
and whether you use smearing or not. probably the default set by pw.x
is sufficient and
you can probably get away with fixed occupations and gamma, but that
just needs to be
tested. there is a wealth of material on how to set up proper jobs
with pw.x and other codes,
also discussing convergence issues (which cannot be predefined) at:
http://www.quantum-espresso.org/wiki/index.php/Tutorials_and_Examples
perhaps you shoud study this more closely. if all these rather
elementary parameters
tell you nothing, than you should work on your understanding of the
methods you use
_before_ you use them for a production calculation...
cheers,
axel.
>
> thanks for the help
>
> Jonas Baltrusaitis
> University of Iowa
>
>
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
>
--
=======================================================================
Axel Kohlmeyer akohlmey at cmm.chem.upenn.edu http://www.cmm.upenn.edu
Center for Molecular Modeling -- University of Pennsylvania
Department of Chemistry, 231 S.34th Street, Philadelphia, PA 19104-6323
tel: 1-215-898-1582, fax: 1-215-573-6233, office-tel: 1-215-898-5425
=======================================================================
If you make something idiot-proof, the universe creates a better idiot.
More information about the users
mailing list