Message Boards Message Boards

[WSS20]Design & Implementation of Data Structure of Crystalline Atomic Data

frontimage

POSTED BY: Claudio Chaib
6 Replies

Could this be used to make the plesiohedra?

POSTED BY: Ed Pegg

No, the fractional lattice coordinates of the atoms and the symmetry elements of the space group are part of the input data in the CIF file. It could be adapted to tile 3D space with plesiohedra given the coordinates of their vertices, lists of vertices comprising the faces, and the elements of the corresponding space group.

For those not familiar with plesiohedra, see also the Wikipedia page.

POSTED BY: Robert Nachbar

enter image description here -- you have earned Featured Contributor Badge enter image description here

Your exceptional post has been selected for our editorial column Staff Picks http://wolfr.am/StaffPicks and Your Profile is now distinguished by a Featured Contributor Badge and is displayed on the Featured Contributor Board. Thank you!

POSTED BY: EDITORIAL BOARD

Update 2 - EARTH’S MOST COMPLEX MINERAL:

  • Ewingite: Earth’s most complex mineral

Ewingite is the most complex mineral on Earth, so I thought it would be interesting to use the CrystalLattice function to reproduce it’s visualization.

But there is still no Ewingite .CIF file available. So to make it possible, I manually edited and developed a .CIF file with the Ewingite information (this .CIF is attached for those interested). The information on this mineral was extracted from few articles (the links are at the end of this Update 2 and in the attach).

Image 1 (basic information - see link 1):

1

Image 2 (place of discovery - see link 1):

2

Image 3 (real image of the mineral - golden yellow):

3

  • Codes (.CIF + using CrystalLattice function):

First, what the .CIF of the mineral Ewingite looks like without taking into account the mineral's natural water molecules (NW – No Water):

4

ewNW = CrystalLattice[
  "C:\\Users\\quimi\\OneDrive\\Área de \
Trabalho\\Ewingite\\EwengiteNW.cif"]

5

After creating Ewingite's CrystalLattice object, we apply CrystalLatticeModify to perform spatial symmetry in a unit cell:

ewNWsym = CrystalLatticeModify[ewNW, {"UnitCellCopy", {0, 0, 0}}]

6

After that, we check the density of the object with the density of the X-ray crystallography (note that the value is very close to that of the article):

valueD = CrystalLatticeValue[ewNWsym, 
  "CellDensity"]; UnitConvert[valueD, {"Grams"/"Centimeters"^3}]

7

So, we generated the 3D graph using CrystalLatticePlot3D, note that even without the water molecules this seems to be a very complex mineral:

CrystalLatticePlot3D[ewNWsym]

8

Refining the visualization using the CrystalLatticeModify function in conjunction with the “UnitCellGather” argument. We can better see its structure inside the unit cell:

ewNWsymGather = 
 CrystalLatticeModify[ewNWsym, 
  "UnitCellGather"]; CrystalLatticePlot3D[ewNWsymGather]

9

Now, comparing with the symmetry given in the article (we can see that the information proceeds):

Image 4 (Space Group symmetry - in attachment 1):

10

Finally, to have a glimpse of the real complexity, the coordinates of the oxygen atoms from the water molecules found in the mineral Ewingite were inserted in the .CIF file. Even so, the visualization still lacks the +3500 hydrogen atoms of this water that was not added due to difficulties in the data information:

ew = CrystalLattice[
  "C:\\Users\\quimi\\OneDrive\\Área de \
Trabalho\\Ewingite\\Ewingite.cif"]; ewSym = 
 CrystalLatticeModify[ew, {"UnitCellCopy", {0, 0, 0}}]; plot1 = 
 CrystalLatticePlot3D[ewSym]; ewSymGather = 
 CrystalLatticeModify[ewSym, "UnitCellGather"]; plot2 = 
 CrystalLatticePlot3D[ewSymGather]; {ewSym, plot1, plot2} // Column

11

  • Links:

Ewingite: Earth’s most complex mineral Travis A. Olds1 , Jakub Plasil2 , Anthony R. Kampf3 , Antonio Simonetti1 , Luke R. Sadergaski1 , Yu-Sheng Chen4 and Peter C. Burns1,5: file in attach*.

https://zh.mindat.org/min-47791.html (images 1 and 2)

Thank you.

POSTED BY: Claudio Chaib

It should be emphasized that the data for the structure of pentadiamond are from an electronic structure calculation (reference above) for this theoretical carbon allotrope. Its properties are therefore predicted and yet to be confirmed experimentally.

Moreover, the data were not in the form of a CIF file and were compiled from information in the article and manually encoded as a CrystalLattice object in Wolfram Language. Additional coding was required to generate the 192 symmetry equivalent positions of the F m -3 m space group, which were critical for the construction of the theoretical crystalline network.

The central cavity of the structure shown in the third figure is composed of 24 5-membered rings and 6 8-membered rings, and has a total of 56 carbon atoms. Its point group symmetry is Oh. While the overall shape is spherical, it's best visualized as a cube with an 8-membered ring on each face (shown below), a triquinacene moiety at each vertex (shown below), and a pentalene moiety on each edge (shown below) of the cube. The diameter of the cavity is 6.50 Angstroms from pentalene to pentalene, and are large enough to contain a helium atom. These cavities are interconnected by cubic chambers composed solely of six 8-membered rings. The diameter of the chamber is 5.15 Angstroms.

Face view: enter image description here

Vertex view: enter image description here

Edge view: enter image description here

POSTED BY: Robert Nachbar

Update NEW CARBON ALLOTROPE

  • A new form of carbon allotrope was developed (30 June 2020 - link at the end of this post), harder and lighter than diamond; the pentadiamond! And I decided to update my work to demonstrate the use of the functions already developed, CrystalLattice, CrystalLatticeModify, CrystalLatticePlot3D.., for the visualization of this new crystal.

i1

  • First, the CrystalLattice object is created, and the symmetry group is checked:

An attached notebook contains the pentadiamond object!

1a

Then:

pentadiamond["SpaceGroupName"]
  • The symmetry group is the same as in the article (as expected!):

1b

There are only 3 unequal carbons in position:

2

  • Using CrystalLatticeModify with "UnitCellCopy" {0,0,0} we obtain the symmetries together:

    pentadiamondSymmetry = 
     CrystalLatticeModify[pentadiamond, {"UnitCellCopy", {0, 0, 0}}]
    

i3

  • Then CrystalLatticeModify is used again with the argument "UnitCellGather" to have the same display as the article of "pentadiamond":

    pentadiamondGather = 
      CrystalLatticeModify[pentadiamondSymmetry, "UnitCellGather"];
    CrystalLatticePlot3D[pentadiamondGather]
    

i4

The article is: https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.125.016001 from Yasumaru Fujii, Mina Maruyama, Nguyen Thanh Cuong, and Susumu Okada Phys. Rev. Lett. 125, 016001 – Published 30 June 2020

I would like to point out that Robert Nachbar (Wolfram), my mentor, who provided the fundamental pieces to obtain the coordinates of this allotrope. Thank you!

Thank you all!

Attachments:
POSTED BY: Claudio Chaib
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract