#!/usr/bin/perl -w

use System;

my $a=1.0;
my $m=2.3266804e-25;

$cell=[
  $a,$a,0,
  $a,0,$a,
  0,-$a,-$a
];

my $atoms=[
  [ "Ce", $m ],
];

my $sites=[
  [ 0,0,0,             0 ],
];

my $bonds=[
  [ 0,0, 1,1,0,                   4.3726,  4.5798,  0.0000,
                                  4.5798,  4.3726,  0.0000,
                                  0.0000,  0.0000, -0.2264 ],

  [ 0,0, 2,0,0,                  -2.3562,  0.0000,  0.0000,
                                  0.0000,  0.0773,  0.0000,
                                  0.0000,  0.0000,  0.0773 ],

  [ 0,0, 2,1,1,                   0.2058, -0.0496, -0.0496,
                                 -0.0496,  0.3169, -0.0547,
                                 -0.0496, -0.0547,  0.3169 ],

  [ 0,0, 2,2,0,                   0.1231,  0.1505,  0.0000,
                                  0.1505,  0.1231,  0.0000,
                                  0.0000,  0.0000,  0.0114 ],

  [ 0,0, 3,1,0,                  -0.0525,  0.0193,  0.0000,
                                  0.0193, -0.0992,  0.0000,
                                  0.0000,  0.0000, -0.1044 ],

  [ 0,0, 2,2,2,                  -0.3316, -0.2194, -0.2194,
                                 -0.2194, -0.3316, -0.2194,
                                 -0.2194, -0.2194, -0.3316 ],

  [ 0,0, 3,2,1,                   0.1057, -0.0068,  0.0763,
                                 -0.0068, -0.1138,  0.0050,
                                  0.0763,  0.0050,  0.0263 ],

  [ 0,0, 4,0,0,                  -0.0009,  0.0000,  0.0000,
                                  0.0000,  0.2219,  0.0000,
                                  0.0000,  0.0000,  0.2219 ],
];

System::write($cell,$atoms,$sites,$bonds,"cubic");
