element(s): ['C', 'H'] AFLOW prototype label: AB_cI16_199_a_a Parameter names: ['a', 'x1', 'x2'] model type (only 'standard' supported at this time): standard number of parameter sets: 1 Parameter values for parameter set 0: ['4.3206', '0.088952473', '0.44338402'] model name: Sim_LAMMPS_ReaxFF_ChenowethVanDuinGoddard_2008_CHO__SM_584143153761_001 ==== Building ASE atoms object with: ==== representative atom symbols = ['C', 'H'] representative atom coordinates = [[0.08895247 0.08895247 0.08895247] [0.44338402 0.44338402 0.44338402]] spacegroup = 199 cell = [[4.3206, 0, 0], [0, 4.3206, 0], [0, 0, 4.3206]] ========================================= Step Time Energy fmax BFGS: 0 12:16:18 -87.873036 1.709038 BFGS: 1 12:16:18 -88.027668 1.448738 BFGS: 2 12:16:19 -88.232716 1.101834 BFGS: 3 12:16:20 -88.365121 1.055347 BFGS: 4 12:16:20 -88.412694 0.856407 BFGS: 5 12:16:21 -88.506627 0.722248 BFGS: 6 12:16:21 -88.614015 0.759988 BFGS: 7 12:16:22 -88.720045 0.691494 BFGS: 8 12:16:22 -88.814236 0.588367 BFGS: 9 12:16:23 -88.893978 0.468593 BFGS: 10 12:16:23 -88.954834 0.356187 BFGS: 11 12:16:24 -88.998286 0.239917 BFGS: 12 12:16:25 -89.023671 0.123768 BFGS: 13 12:16:25 -89.030696 0.039220 BFGS: 14 12:16:26 -89.030747 0.006301 BFGS: 15 12:16:26 -89.030748 0.002024 BFGS: 16 12:16:26 -89.030749 0.000018 BFGS: 17 12:16:27 -89.030749 0.000001 BFGS: 18 12:16:27 -89.030749 0.000000 Minimization converged after 18 steps. Maximum force component: 1.2030085029459912e-09 eV/Angstrom Maximum stress component: 3.732933129127767e-10 eV/Angstrom^3 ==== Minimized structure obtained from ASE ==== symbols = ['C', 'C', 'C', 'C', 'C', 'C', 'C', 'C', 'H', 'H', 'H', 'H', 'H', 'H', 'H', 'H'] basis = [[0.07899897 0.07899897 0.07899897] [0.42100103 0.92100103 0.57899897] [0.92100103 0.57899897 0.42100103] [0.57899897 0.42100103 0.92100103] [0.57899897 0.57899897 0.57899897] [0.92100103 0.42100103 0.07899897] [0.42100103 0.07899897 0.92100103] [0.07899897 0.92100103 0.42100103] [0.43059741 0.43059741 0.43059741] [0.06940259 0.56940259 0.93059741] [0.56940259 0.93059741 0.06940259] [0.93059741 0.06940259 0.56940259] [0.93059741 0.93059741 0.93059741] [0.56940259 0.06940259 0.43059741] [0.06940259 0.43059741 0.56940259] [0.43059741 0.56940259 0.06940259]] cellpar = Cell([4.164100454529655, 4.164100454529653, 4.164100454529655]) forces = [[ 1.20300850e-09 1.20300850e-09 1.20300850e-09] [-1.20300850e-09 -1.20300850e-09 1.20300850e-09] [-1.20300850e-09 1.20300850e-09 -1.20300850e-09] [ 1.20300850e-09 -1.20300850e-09 -1.20300850e-09] [ 1.20300850e-09 1.20300850e-09 1.20300850e-09] [-1.20300850e-09 -1.20300850e-09 1.20300850e-09] [-1.20300850e-09 1.20300850e-09 -1.20300850e-09] [ 1.20300850e-09 -1.20300850e-09 -1.20300850e-09] [ 2.40430966e-10 2.40430966e-10 2.40430966e-10] [-2.40430966e-10 -2.40430966e-10 2.40430966e-10] [-2.40430966e-10 2.40430966e-10 -2.40430966e-10] [ 2.40430966e-10 -2.40430966e-10 -2.40430966e-10] [ 2.40430966e-10 2.40430966e-10 2.40430966e-10] [-2.40430966e-10 -2.40430966e-10 2.40430966e-10] [-2.40430966e-10 2.40430966e-10 -2.40430966e-10] [ 2.40430966e-10 -2.40430966e-10 -2.40430966e-10]] stress = [3.73293313e-10 3.73293313e-10 3.73293313e-10 0.00000000e+00 0.00000000e+00 0.00000000e+00] energy per atom = -5.510303768071336 =============================================== Parameter sets [0] relaxed to duplicate structures, attempting to write only one of them. Successfully added property instance for parameter set 0