Property Changes in a Packed Bed

I am newer to LIGGGHTS so I am still learning my way around assigning properties to particles. In my simulations the particles are spheres, and they are in a packed bed. Over time the particle's properties need to change based upon their position, in particular radius, Youngs Modulus and Poissons ratio. I am using the granular hertzian pair style without history.

So far as I have seen I can create a material type with specific properties, and then assign each sphere that material type. But I was wondering if there was a way to output data that includes each sphere's youngs modulus and poissons ratio, and then to alter them which would allow each sphere to have slightly different material properties (and radius) and then re-read them back in or use them for the liggghts calculations.

Thanks!