Group Abstract Group Abstract

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
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