Hello everyone,

The stress increases unphysically for Tobermorite-14A simulation using CLAYFF potential. I am using Morse potential for bond_style. But, it should fracture pretty fast somewhere along 0.1-0.2 strain. I have water inside tobermorite. I have the water to see the difference. But it doesn’t make any difference. Here’s a part of my code:

units real

atom_style full

boundary p p p

# Bonded Interaction

bond_style morse

angle_style harmonic

dihedral_style harmonic

improper_style cvff

# Calculation of stresses and strains # per-atom stress is sum of diagonal components in stress tensor

variable multiplication_factor equal 101.325*1e-6 # stress converstion from atm to GPa
variable f equal {multiplication_factor}
variable sigmax equal (v_f)*(c_p[1]/{Tvol})
variable sigmay equal (v_f)*(c_p[2]/{Tvol})
variable sigmaz equal (v_f)*(c_p[3]/{Tvol})

# Applying axial and shear strain by fix deform

#fix 8 all nvt temp {temp} {temp} 5.0 drag 1

fix 8 solid npt temp {temp} {temp} 5.0 x 0 0 100 y 0 0 100 drag 1

fix 9 all deform 1 z erate ${erate1} units box remap x

Regards,

Baig