Given a impulse responce can we derive a state space model?
you might need some procedure like the prony method to convert your impulse response (as a bandlimited and uniformly-sampled signal) into a sum of decaying exponentials and from that obtain a sum of partial fractions of first and second-order transfer functions, and from that derive poles and zeros.  from that can be derived a state-space model.