[ase-users] From CIF files to XYZ

[Eric Hermes]

On Fri, 2017-12-08 at 18:48 +0000, Landi, Alessandro via ase-users
wrote:
> Hello.
> 
> Obviously, I typed it wrong.
> 
> Here the right command lines I use:
> 
> atoms_ase=ase.io.read('filename.cif')
> ase.io.write('filename.xyz', atoms_ase)
> 
> Anyway, I have the problem that the xyz file is wrong as I explained
> in my first e-mail: the molecules need to be reconnected (see files
> attached).

I don't understand what you mean by "reconnected". Neither xyz file
that you uploaded contains any information regarding connectivity, and
ASE doesn't have any conception of connectivity anyway (aside from a
distance-based metric used in the ASE GUI to draw sticks between nearby
atoms). The atom ordering produced by ASE is more consistent than the
"realstructure.xyz" file as compared to the original cif ordering. It
might be that whatever program you are using with the xyz files expects
connected atoms to be contiguous; to the best of my knowledge this is
not typical. If you want to ensure this ordering within ASE, you will
need to manually rearrange the Atoms object.

Eric

> 
> "pentacene.cif" is the cif file that I give as input to ASE;
> "pentacene.xyz" is the xyz file written by ASE;
> "realstructure.xyz" is the structure of the reconnected molecules
> (obtained using CCDC Mercury)
> 
> Any suggestions?
> 
> Best regards,
> Alessandro
> 
> Da: Oscar Xavier Guerrero <oscarxavier.ox at gmail.com>
> Inviato: venerdì 8 dicembre 2017 18:36:10
> A: Landi, Alessandro
> Cc: ase-users at listserv.fysik.dtu.dk
> Oggetto: Re: [ase-users] From CIF files to XYZ
>  
> Hello,
> 
> When you read the structure an Atoms object is created. You need to
> store that object on a variable and then you can write that one. As
> so:
> 
> First read the data:
> from ase.io import read, write
> atoms = read('filename.cif')
> Then write:
> write('filename.xyz', atoms)
> or you can use:
> atoms.write('filename.xyz')
> 
> 
> 
> 
> 2017-12-08 12:19 GMT+00:00 Landi, Alessandro via ase-users <ase-users
> @listserv.fysik.dtu.dk>:
> > Hello everyone.
> > 
> > I am using ASE to read some CIF files and (among other things) to
> > convert the structure in xyz format.
> > 
> > However using the following commands:
> > 
> > ase.io.read('filename.cif', atoms_ase)
> > ase.io.write('filename.xyz', atoms_ase)
> > 
> > I obtain a wrong structure, since the molecules need to be
> > reconnected. This happens for all the molecules I have tested
> > (about 20).
> > 
> > Could you please explain me what I am doing wrong?
> > 
> > If could be of any help, I attached an example:
> > "pentacene.cif" is the cif file that I give as input to ASE;
> > "pentacene.xyz" is the xyz file written by ASE;
> > "realstructure.xyz" is the structure of the reconnected molecules
> > (obtained using CCDC Mercury)
> > 
> > Best regards,
> > Alessandro
> > 
> > 
> > _______________________________________________
> > ase-users mailing list
> > ase-users at listserv.fysik.dtu.dk
> > https://listserv.fysik.dtu.dk/mailman/listinfo/ase-users
> > 
> 
> _______________________________________________
> ase-users mailing list
> ase-users at listserv.fysik.dtu.dk
> https://listserv.fysik.dtu.dk/mailman/listinfo/ase-users