Skip to content

Commit 64b8bbe

Browse files
committed
Merge branch 'ml-dft-sandia-master'
2 parents 544c1bf + c9ef077 commit 64b8bbe

File tree

219 files changed

+57862
-1
lines changed

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

219 files changed

+57862
-1
lines changed

ml-dft-sandia/README.md

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,3 +1,3 @@
11
# ML-DFT-SANDIA
22

3-
The Sandia code goes here.
3+
This is the original code used in [1].
Lines changed: 38 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,38 @@
1+
This example shows how to generate the SNAP bispectrum on a grid.
2+
3+
Usage:
4+
5+
python3.7 ./bgrid.py > bgrid.out
6+
7+
Requirements:
8+
9+
Python 3.6 or 3.7
10+
LAMMPS, with compute sna/grid, from here:
11+
https://github.com/athomps/lammps/tree/compute-grid
12+
13+
Before running, you will need to do the following:
14+
15+
1. Compile the LAMMPS shared library with the SNAP package installed
16+
2. Make the shared library visible e.g.
17+
18+
ln -s $LAMMPS/src/liblammps.so .
19+
(for other ways see https://lammps.sandia.gov/doc/Python_shlib.html)
20+
21+
3. Make the LAMMPS lammps.py visible e.g.
22+
23+
ln -s $LAMMPS/python/lammps.py .
24+
(for better ways see https://lammps.sandia.gov/doc/Python_install.html)
25+
26+
The files that are produced are as follows:
27+
28+
bgrid.out: Screen output from LAMMPS and bgrid.py
29+
log.lammps: LAMMPS log file
30+
b.dump: per-atom bispectrum components
31+
bgrid.dat: grid bispectrum components as LAMMPS array
32+
bgrid_np.dat: grid bispectrum components as Numpy array
33+
pgrid_np.dat: grid power spectrum components as Numpy array
34+
bgridnorm_np.dat: normalized bispectrum
35+
pgridnorm_np.dat: normalized power spectrum
36+
37+
You can compare these to the corresponding files
38+
with names containing a timestamp e.g. b.25Oct2019.dump
Lines changed: 11 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,11 @@
1+
ITEM: TIMESTEP
2+
0
3+
ITEM: NUMBER OF ATOMS
4+
2
5+
ITEM: BOX BOUNDS pp pp pp
6+
0.0000000000000000e+00 3.3159999999999998e+00
7+
0.0000000000000000e+00 3.3159999999999998e+00
8+
0.0000000000000000e+00 3.3159999999999998e+00
9+
ITEM: ATOMS id x y z c_b[1] c_b[2] c_b[3] c_b[4] c_b[5]
10+
1 0 0 0 108.173 3.21778 0.712238 7.06634 1.04273
11+
2 1.658 1.658 1.658 108.173 3.21778 0.712238 7.06634 1.04273
Lines changed: 11 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,11 @@
1+
ITEM: TIMESTEP
2+
0
3+
ITEM: NUMBER OF ATOMS
4+
2
5+
ITEM: BOX BOUNDS pp pp pp
6+
0.0000000000000000e+00 3.3159999999999998e+00
7+
0.0000000000000000e+00 3.3159999999999998e+00
8+
0.0000000000000000e+00 3.3159999999999998e+00
9+
ITEM: ATOMS id c_b[1] c_b[2] c_b[3] c_b[4] c_b[5]
10+
1 3375 846.026 36.811 8.51924 0.247118
11+
2 3375 846.026 36.811 8.51924 0.247118
Lines changed: 128 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,128 @@
1+
x, y, z = 0 0 0
2+
108.173 3.21778 0.712238 7.06634 1.04273
3+
x, y, z = 1.658 0 0
4+
174.4 5.59324 1.35171 13.8088 2.13628
5+
x, y, z = 0 0.829 0
6+
183.116 3.64679 1.14854 30.7642 7.26691
7+
x, y, z = 1.658 0.829 0
8+
168.973 6.62458 1.46561 11.0631 1.61839
9+
x, y, z = 0 1.658 0
10+
174.4 5.59324 1.35171 13.8088 2.13628
11+
x, y, z = 1.658 1.658 0
12+
174.4 5.59324 1.35171 13.8088 2.13628
13+
x, y, z = 0 2.487 0
14+
183.116 3.64679 1.14854 30.7642 7.26691
15+
x, y, z = 1.658 2.487 0
16+
168.973 6.62458 1.46561 11.0631 1.61839
17+
x, y, z = 0 0 0.4145
18+
189.401 2.43052 0.962196 44.4583 12.4353
19+
x, y, z = 1.658 0 0.4145
20+
172.181 6.01397 1.39915 12.5923 1.90608
21+
x, y, z = 0 0.829 0.4145
22+
180.705 4.08676 1.22596 27.1833 6.05503
23+
x, y, z = 1.658 0.829 0.4145
24+
170.986 6.15471 1.45038 13.3359 2.10603
25+
x, y, z = 0 1.658 0.4145
26+
172.181 6.01397 1.39915 12.5923 1.90608
27+
x, y, z = 1.658 1.658 0.4145
28+
175.99 5.14283 1.33312 18.0104 3.26169
29+
x, y, z = 0 2.487 0.4145
30+
180.705 4.08676 1.22596 27.1833 6.05503
31+
x, y, z = 1.658 2.487 0.4145
32+
170.986 6.15471 1.45038 13.3359 2.10603
33+
x, y, z = 0 0 0.829
34+
183.116 3.64679 1.14854 30.7642 7.26691
35+
x, y, z = 1.658 0 0.829
36+
168.973 6.62458 1.46561 11.0631 1.61839
37+
x, y, z = 0 0.829 0.829
38+
175.638 5.11301 1.3681 19.6169 3.72319
39+
x, y, z = 1.658 0.829 0.829
40+
175.638 5.11301 1.3681 19.6169 3.72319
41+
x, y, z = 0 1.658 0.829
42+
168.973 6.62458 1.46561 11.0631 1.61839
43+
x, y, z = 1.658 1.658 0.829
44+
183.116 3.64679 1.14854 30.7642 7.26691
45+
x, y, z = 0 2.487 0.829
46+
175.638 5.11301 1.3681 19.6169 3.72319
47+
x, y, z = 1.658 2.487 0.829
48+
175.638 5.11301 1.3681 19.6169 3.72319
49+
x, y, z = 0 0 1.2435
50+
175.99 5.14283 1.33312 18.0104 3.26169
51+
x, y, z = 1.658 0 1.2435
52+
172.181 6.01397 1.39915 12.5923 1.90608
53+
x, y, z = 0 0.829 1.2435
54+
170.986 6.15471 1.45038 13.3359 2.10603
55+
x, y, z = 1.658 0.829 1.2435
56+
180.705 4.08676 1.22596 27.1833 6.05503
57+
x, y, z = 0 1.658 1.2435
58+
172.181 6.01397 1.39915 12.5923 1.90608
59+
x, y, z = 1.658 1.658 1.2435
60+
189.401 2.43052 0.962196 44.4583 12.4353
61+
x, y, z = 0 2.487 1.2435
62+
170.986 6.15471 1.45038 13.3359 2.10603
63+
x, y, z = 1.658 2.487 1.2435
64+
180.705 4.08676 1.22596 27.1833 6.05503
65+
x, y, z = 0 0 1.658
66+
174.4 5.59324 1.35171 13.8088 2.13628
67+
x, y, z = 1.658 0 1.658
68+
174.4 5.59324 1.35171 13.8088 2.13628
69+
x, y, z = 0 0.829 1.658
70+
168.973 6.62458 1.46561 11.0631 1.61839
71+
x, y, z = 1.658 0.829 1.658
72+
183.116 3.64679 1.14854 30.7642 7.26691
73+
x, y, z = 0 1.658 1.658
74+
174.4 5.59324 1.35171 13.8088 2.13628
75+
x, y, z = 1.658 1.658 1.658
76+
108.173 3.21778 0.712238 7.06634 1.04273
77+
x, y, z = 0 2.487 1.658
78+
168.973 6.62458 1.46561 11.0631 1.61839
79+
x, y, z = 1.658 2.487 1.658
80+
183.116 3.64679 1.14854 30.7642 7.26691
81+
x, y, z = 0 0 2.0725
82+
175.99 5.14283 1.33312 18.0104 3.26169
83+
x, y, z = 1.658 0 2.0725
84+
172.181 6.01397 1.39915 12.5923 1.90608
85+
x, y, z = 0 0.829 2.0725
86+
170.986 6.15471 1.45038 13.3359 2.10603
87+
x, y, z = 1.658 0.829 2.0725
88+
180.705 4.08676 1.22596 27.1833 6.05503
89+
x, y, z = 0 1.658 2.0725
90+
172.181 6.01397 1.39915 12.5923 1.90608
91+
x, y, z = 1.658 1.658 2.0725
92+
189.401 2.43052 0.962196 44.4583 12.4353
93+
x, y, z = 0 2.487 2.0725
94+
170.986 6.15471 1.45038 13.3359 2.10603
95+
x, y, z = 1.658 2.487 2.0725
96+
180.705 4.08676 1.22596 27.1833 6.05503
97+
x, y, z = 0 0 2.487
98+
183.116 3.64679 1.14854 30.7642 7.26691
99+
x, y, z = 1.658 0 2.487
100+
168.973 6.62458 1.46561 11.0631 1.61839
101+
x, y, z = 0 0.829 2.487
102+
175.638 5.11301 1.3681 19.6169 3.72319
103+
x, y, z = 1.658 0.829 2.487
104+
175.638 5.11301 1.3681 19.6169 3.72319
105+
x, y, z = 0 1.658 2.487
106+
168.973 6.62458 1.46561 11.0631 1.61839
107+
x, y, z = 1.658 1.658 2.487
108+
183.116 3.64679 1.14854 30.7642 7.26691
109+
x, y, z = 0 2.487 2.487
110+
175.638 5.11301 1.3681 19.6169 3.72319
111+
x, y, z = 1.658 2.487 2.487
112+
175.638 5.11301 1.3681 19.6169 3.72319
113+
x, y, z = 0 0 2.9015
114+
189.401 2.43052 0.962196 44.4583 12.4353
115+
x, y, z = 1.658 0 2.9015
116+
172.181 6.01397 1.39915 12.5923 1.90608
117+
x, y, z = 0 0.829 2.9015
118+
180.705 4.08676 1.22596 27.1833 6.05503
119+
x, y, z = 1.658 0.829 2.9015
120+
170.986 6.15471 1.45038 13.3359 2.10603
121+
x, y, z = 0 1.658 2.9015
122+
172.181 6.01397 1.39915 12.5923 1.90608
123+
x, y, z = 1.658 1.658 2.9015
124+
175.99 5.14283 1.33312 18.0104 3.26169
125+
x, y, z = 0 2.487 2.9015
126+
180.705 4.08676 1.22596 27.1833 6.05503
127+
x, y, z = 1.658 2.487 2.9015
128+
170.986 6.15471 1.45038 13.3359 2.10603

0 commit comments

Comments
 (0)