element(s): ['Ca', 'Si'] AFLOW prototype label: AB2_hR6_166_c_2c Parameter names: ['a', 'c/a', 'x1', 'x2', 'x3'] model type (only 'standard' supported at this time): standard number of parameter sets: 1 Parameter values for parameter set 0: ['3.8848956', '7.8159275', '0.91996424', '0.81742379', '0.65322287'] model name: Sim_LAMMPS_ReaxFF_ManzanoMoeiniMarinelli_2012_CaSiOH__SM_714124634215_000 ==== Building ASE atoms object with: ==== representative atom symbols = ['Ca', 'Si', 'Si'] representative atom coordinates = [[0. 0. 0.41996424] [0. 0. 0.31742379] [0. 0. 0.15322287]] spacegroup = 166 cell = [[3.8849, 0, 0], [-1.94245, 3.3644220911622, 0], [0, 0, 30.3641]] ========================================= Step Time Energy fmax BFGS: 0 16:22:10 -69.334117 2.913354 BFGS: 1 16:22:11 -69.643444 2.847716 BFGS: 2 16:22:11 -69.985475 2.771288 BFGS: 3 16:22:11 -70.304705 2.693717 BFGS: 4 16:22:11 -70.606134 2.611946 BFGS: 5 16:22:11 -70.892373 2.522781 BFGS: 6 16:22:12 -71.164414 2.422613 BFGS: 7 16:22:12 -71.421994 2.307299 BFGS: 8 16:22:12 -71.663782 2.172286 BFGS: 9 16:22:12 -71.887547 2.013136 BFGS: 10 16:22:13 -72.090476 1.826536 BFGS: 11 16:22:13 -72.269757 1.611739 BFGS: 12 16:22:13 -72.423404 1.371937 BFGS: 13 16:22:13 -72.551118 1.114753 BFGS: 14 16:22:14 -72.654691 0.851121 BFGS: 15 16:22:14 -72.737646 0.592816 BFGS: 16 16:22:14 -72.804208 0.417424 BFGS: 17 16:22:14 -72.858195 0.545900 BFGS: 18 16:22:15 -72.902486 0.655555 BFGS: 19 16:22:15 -72.939653 0.736376 BFGS: 20 16:22:15 -72.966314 0.726267 BFGS: 21 16:22:15 -73.000884 0.630671 BFGS: 22 16:22:15 -73.028489 0.516323 BFGS: 23 16:22:16 -73.052880 0.441194 BFGS: 24 16:22:16 -73.077600 0.399324 BFGS: 25 16:22:16 -73.103270 0.376558 BFGS: 26 16:22:16 -73.129553 0.391938 BFGS: 27 16:22:17 -73.099763 0.443664 BFGS: 28 16:22:17 -73.137321 0.413713 BFGS: 29 16:22:17 -73.171144 0.413336 BFGS: 30 16:22:17 -73.201604 0.419362 BFGS: 31 16:22:18 -73.229190 0.422833 BFGS: 32 16:22:18 -73.254345 0.423561 BFGS: 33 16:22:18 -73.277449 0.421410 BFGS: 34 16:22:19 -73.298820 0.416289 BFGS: 35 16:22:19 -73.318716 0.408141 BFGS: 36 16:22:19 -73.337334 0.396936 BFGS: 37 16:22:19 -73.354826 0.382666 BFGS: 38 16:22:20 -73.371293 0.365334 BFGS: 39 16:22:20 -73.386801 0.344949 BFGS: 40 16:22:20 -73.401384 0.321522 BFGS: 41 16:22:20 -73.415048 0.295048 BFGS: 42 16:22:21 -73.427781 0.265505 BFGS: 43 16:22:21 -73.439553 0.232823 BFGS: 44 16:22:21 -73.450325 0.196848 BFGS: 45 16:22:22 -73.460052 0.157255 BFGS: 46 16:22:22 -73.468683 0.139600 BFGS: 47 16:22:22 -73.476156 0.142709 BFGS: 48 16:22:23 -73.481507 0.124724 BFGS: 49 16:22:23 -73.484288 0.091503 BFGS: 50 16:22:23 -73.486243 0.041512 BFGS: 51 16:22:24 -73.486722 0.024337 BFGS: 52 16:22:24 -73.486870 0.014112 BFGS: 53 16:22:24 -73.486923 0.006198 BFGS: 54 16:22:24 -73.486937 0.002077 BFGS: 55 16:22:25 -73.486932 0.000819 BFGS: 56 16:22:25 -73.486931 0.000390 BFGS: 57 16:22:25 -73.486931 0.000365 BFGS: 58 16:22:26 -73.486932 0.000222 BFGS: 59 16:22:26 -73.486932 0.000091 BFGS: 60 16:22:26 -73.486932 0.000011 BFGS: 61 16:22:26 -73.486932 0.000001 BFGS: 62 16:22:27 -73.486932 0.000000 BFGS: 63 16:22:27 -73.486932 0.000000 Minimization converged after 63 steps. Maximum force component: 2.969683046744984e-09 eV/Angstrom Maximum stress component: 1.0616876912751443e-10 eV/Angstrom^3 ==== Minimized structure obtained from ASE ==== symbols = ['Ca', 'Ca', 'Ca', 'Ca', 'Ca', 'Ca', 'Si', 'Si', 'Si', 'Si', 'Si', 'Si', 'Si', 'Si', 'Si', 'Si', 'Si', 'Si'] basis = [[3.38035849e-32 5.04626786e-31 4.17315368e-01] [1.96357103e-17 1.31995688e-16 5.82684632e-01] [6.66666667e-01 3.33333333e-01 7.50648701e-01] [6.66666667e-01 3.33333333e-01 9.16017966e-01] [3.33333333e-01 6.66666667e-01 8.39820343e-02] [3.33333333e-01 6.66666667e-01 2.49351299e-01] [1.37393736e-17 0.00000000e+00 3.18762574e-01] [3.90131959e-17 1.00925357e-30 6.81237426e-01] [6.66666667e-01 3.33333333e-01 6.52095907e-01] [6.66666667e-01 3.33333333e-01 1.45707592e-02] [3.33333333e-01 6.66666667e-01 9.85429241e-01] [3.33333333e-01 6.66666667e-01 3.47904093e-01] [4.35433706e-18 1.57695870e-31 1.52517037e-01] [4.83982325e-17 5.04626786e-31 8.47482963e-01] [6.66666667e-01 3.33333333e-01 4.85850371e-01] [6.66666667e-01 3.33333333e-01 1.80816296e-01] [3.33333333e-01 6.66666667e-01 8.19183704e-01] [3.33333333e-01 6.66666667e-01 5.14149629e-01]] cellpar = Cell([[3.6101853413585743, 2.2648954378131503e-17, -1.0767736236975243e-15], [-1.8050926706792867, 3.126512217986714, 7.53700282381232e-16], [-9.143259342222432e-15, 2.638256566412998e-15, 32.3885924380764]]) forces = [[-6.27698131e-26 1.81120173e-26 2.22352426e-10] [ 6.27696549e-26 -1.81118803e-26 -2.22352426e-10] [-6.27698280e-26 1.81120173e-26 2.22352426e-10] [ 6.27696846e-26 -1.81118974e-26 -2.22352426e-10] [-6.27698230e-26 1.81120173e-26 2.22352426e-10] [ 6.27696599e-26 -1.81118803e-26 -2.22352426e-10] [-3.22246087e-25 9.29824479e-26 1.14150533e-09] [ 3.22246008e-25 -9.29831330e-26 -1.14150533e-09] [-3.22245691e-25 9.29825849e-26 1.14150533e-09] [ 3.22246087e-25 -9.29835440e-26 -1.14150533e-09] [-3.22245770e-25 9.29821738e-26 1.14150533e-09] [ 3.22246482e-25 -9.29834070e-26 -1.14150533e-09] [-8.38337987e-25 2.41899546e-25 2.96968305e-09] [ 8.38338422e-25 -2.41900025e-25 -2.96968305e-09] [-8.38337947e-25 2.41899751e-25 2.96968305e-09] [ 8.38338303e-25 -2.41900094e-25 -2.96968305e-09] [-8.38337789e-25 2.41899614e-25 2.96968305e-09] [ 8.38338501e-25 -2.41900162e-25 -2.96968305e-09]] stress = [-9.23033281e-11 -9.23033281e-11 -1.06168769e-10 1.48302119e-26 -5.77230102e-26 -5.16938786e-25] energy per atom = -3.99862116389727 =============================================== Parameter sets [0] relaxed to duplicate structures, attempting to write only one of them. Successfully added property instance for parameter set 0