A nonlocal constitutive model generated by matrix functions for polyatomic periodic linear chains