News:

SMF - Just Installed!

Main Menu

Only one Framework atom taken from .cif file

Started by petersbk, April 01, 2019, 05:25:26 PM

Previous topic - Next topic

petersbk

Summary: Large number of FW atoms in .cif file, but simulation only uses 1 of them and at position 0,0,0. What is wrong with .cif or input file?
Warning: Stage 0 Newbie with respect to raspa.

I'm using a bunch of Argon atoms in a P1 40x40x40 box as a "Framework": just a bunch of fixed LJ atoms via the Garcia Perez FF. The simulation attempts gcmc of methane in this framework.  The .cif was made in Avogadro / openbabel after creating the atom positions elsewhere. Here is the top of the .cif:
--------------------------------------------------------------------------------------------
# CIF file generated by openbabel 2.3.2, see http://openbabel.sf.net
data_I
_chemical_name_common 'box6_ar.pdb'
_cell_length_a 40
_cell_length_b 40
_cell_length_c 40
_cell_angle_alpha 90
_cell_angle_beta 90
_cell_angle_gamma 90
_space_group_name_H-M_alt 'P 1'
_space_group_name_Hall 'P 1'
loop_
    _symmetry_equiv_pos_as_xyz
    'x,y,z'
loop_
    _atom_site_type_symbol
    _atom_site_label
    _atom_site_Cartn_x
    _atom_site_Cartn_y
    _atom_site_Cartn_z
     Ar   Ar1    19.87500    0.95800   36.89600
     Ar   Ar2    29.92900    6.71600   36.71700
     Ar   Ar3     6.11900    9.36600    1.67600
...
--------------------------------------------------------------------------------------------

And I'm using this simple input file (which started life as one of the Basic examples):

--------------------------------------------------------------------------------------------
SimulationType                MonteCarlo
NumberOfCycles                5000
NumberOfInitializationCycles  1000
PrintEvery                    1000

Forcefield                    GarciaPerez2006

RemoveAtomNumberCodeFromLabel yes
ReadCIFAsCartesian yes
UseChargesFromCIFFile no

Framework 0
FrameworkName box6_ar
ExternalTemperature 300.0
ExternalPressure 1e8

Component 0 MoleculeName             methane
            MoleculeDefinition       TraPPE
            TranslationProbability   0.5
            ReinsertionProbability   0.5
            SwapProbability          1.0
            CreateNumberOfMolecules  0
--------------------------------------------------------------------------------------------

When the simulation is run, this is the start of what is printed to stderr:

--------------------------------------------------------------------------------------------
_cell_length_a: 40.000000
_cell_length_b: 40.000000
_cell_length_c: 40.000000
_cell_length_alpha: 90.000000
_cell_length_beta: 90.000000
_cell_length_gamma: 90.000000
space group found from symmetry elements: 1 (nr elements: 1)
unknown: 19.87500
unknown: 0.95800
unknown: 36.89600
unknown: 29.92900
unknown: 6.71600
...
--------------------------------------------------------------------------------------------

There is a separate one of these "unknown" lines for each x, y, and z coordinate in the .cif file.

The simulation runs "fine", but from the _final.pdb etc. files in the Movies directory or from the output, one can see that only one framework atom (@ 0,0,0) was used in the simulation.

What is my problem (at least related to this simulation...)?
Thanks.

Pezhman Zarabadi-Poor

Dear Peter,

I think your issue is coming from the fact that you are providing Cartesian coordinates instead of fractional ones. You may convert it to fractional and then you will have

.
.
.
_atom_site_frac_x
_atom_site_frac_y
_atom_site_frac_z
Ar Ar1 frac_x frac_y frac_z
.
.
.

My best regards,
Pezhman

petersbk

Maybe. The manual says you can use Cartesian coordinates for P1 and I've used the "ReadCIFAsCartesian yes" command in the input file. It is off-topic, but I can't get openbabel to convert this file to fractional coordinates.  Thanks.

Pezhman Zarabadi-Poor

I have not used this option of RASPA but I guess it still needs the CIF in fractional coordinates and then RASPA internally converts it to Cartesian. You may use other software like VESTA to generate CIF in fractional coordinates.

My best,
Pezhman

petersbk

I converted to fractional coordinates by hand, removed the "ReadCIFAsCartesian yes" command, and the same behavior occurs: lots of "unknown" lines to stderr and only one FW atom appears @ 0,0,0. Any other advice? Thanks.

Pezhman Zarabadi-Poor

May you please send your CIF, force field definition, and input file?

petersbk

Thanks for looking into this. The input and .cif files can be taken from my original post (the three atoms shown give similar behavior as a larger number of atoms). This force-field.def file is the one supplied in the latest version of the download. I've done nothing to modify any .def file.

petersbk

In the interest of closure, I'll report that I was only able to solve my .cif file problems by using the header from a working P1 .cif file (e.g. from the structures directory in the distribution), cutting out the coordinates, pasting in fractional coordinates for my system of interest, and modifying the cell parameters. Error messages of the type: "unknown: x" where x is one of the x,y, or z coordinates in the .cif file can occur when something is wrong (not compatible with raspa) with the other lines in the .cif.