---------------------------------------------------------------------------- c a u c h y - b o r n c o n s t i t u t i v e m o d e l This program computes the finite deformation constitutive response of an arbitrary multilattice single crystal. The code is based on the local quasicontinuum (QC) method (http://qcmethod.org). The program uses interatomic potentials compatible with the Knowledgebase of Interatomic Models (KIM) project (https://openkim.org). Copyright (C) 2015 E. B. Tadmor ---------------------------------------------------------------------------- Test calculation: Stillinger-Weber silicon ground state and compression **setting factor** * Default globaltol set to 0.10000E-15 CG global tolerance = 0.10000E-15 **reading material data** Material name 'Al-Ni(L12)' material information Bravais lattice vectors: a1 = 4.00000 0.00000 0.00000 a2 = 0.00000 4.00000 0.00000 a3 = 0.00000 0.00000 4.00000 Basis atoms: Number of basis atoms = 4 spec mass reference shift vector 13 1.00 0.00000 0.00000 0.00000 28 1.00 0.00000 2.00000 2.00000 28 1.00 2.00000 0.00000 2.00000 28 1.00 2.00000 2.00000 0.00000 **reading constitutive information** * Error: at line 169 in mod_poten_kim.F90 Message: kim_api_string_init KIM_STATUS_MSG: unsuccessful completion KIM Model name: IMD_EAM_Schopf_AlNiCo_A__MO_122703700223_002 Match failed; removing process_d2Edr2 and trying again. Moduli calculation disabled. Preconditioning disabled in lattice parameter optimization. Maximum cut-off distance: 7.2000000000 **global parameter settings** Number of spatial dimensions = 3 Number of dofs per node = 3 Number of nodes per element = 4 Shuffles optimization = ACTIVATED Prefer multilattice stress correction = OFF Optimize F and shuffles concurrently = ACTIVATED Use preconditioned CG optimization = OFF **check equilibrium of multilattice** Input lattice vectors for CheckEquilibrium 4.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 4.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 4.0000000000000000 Input basis atom positions for CheckEquilibrium 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 2.0000000000000000 2.0000000000000000 2.0000000000000000 0.0000000000000000 2.0000000000000000 2.0000000000000000 2.0000000000000000 0.0000000000000000 scaling parameter: 0.9316432643371826 Updated lattice vector after uniform optimization of lattice parameter: 3.7265730573487303 0.0000000000000000 0.0000000000000000 0.0000000000000000 3.7265730573487303 0.0000000000000000 0.0000000000000000 0.0000000000000000 3.7265730573487303 Updated basis atom positions after uniform optimization of lattice parameter: 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000051714 1.8632865287242122 1.8632865286875797 1.8632865286631897 0.0000000000051711 1.8632865286224738 1.8632865287242175 1.8632865286631950 -0.0000000000051714 dwde for CheckEquilibrium 0.0000000000000000 -0.0000000000688929 -0.0000000007094164 0.0000000000000000 0.0000000009742882 -0.0000000000688937 0.0000000000000000 0.0000000007656732 -0.0000000010306955 -0.0000000000688929 -0.0000000007094164 0.0000000009744395 **optimize lattice parameter** Input lattice vectors 4.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 4.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 4.0000000000000000 Input basis atom positions 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 2.0000000000000000 2.0000000000000000 2.0000000000000000 0.0000000000000000 2.0000000000000000 2.0000000000000000 2.0000000000000000 0.0000000000000000 Relaxed lattice vectors 3.72657305724183229856 0.00000000000000000461 0.00000000000000000170 0.00000000000000000461 3.72657305724183007811 0.00000000000000000475 0.00000000000000000170 0.00000000000000000475 3.72657305724183185447 Relaxed basis atom positions 0.00000000000000000000 0.00000000000000000000 0.00000000000000000000 -0.00000000000000002152 1.86328652862091503906 1.86328652862091592723 1.86328652862091614928 0.00000000000000003495 1.86328652862091592723 1.86328652862091614928 1.86328652862091503906 0.00000000000000000097 **deform crystal** Method: relax shifts Mode: single step Deformation gradient 1.00000000 0.00000000 0.00000000 0.00000000 1.00000000 0.00000000 0.00000000 0.00000000 1.00000000 Shuffle vectors and relaxed displacements 1 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 2 -0.00000000 1.86328653 1.86328653 0.00000000 0.00000000 0.00000000 3 1.86328653 0.00000000 1.86328653 0.00000000 0.00000000 0.00000000 4 1.86328653 1.86328653 0.00000000 0.00000000 0.00000000 0.00000000 energy density (per reference volume) = -0.3576888966 Kirchhoff stress -0.90101E-15 0.21453E-17 0.21453E-17 0.00000E+00 -0.16980E-14 -0.42905E-17 0.21453E-17 -0.53632E-17 0.33681E-15 * CB-KIM calculation successfully completed