Using the spectral theory on the S-spectrum it is possible to define the fractional powers of a large class of vector operators. This possibility leads to new fractional diffusion and evolution problems that are of particular interest for nonhomogeneous materials where the Fourier law is not simply the negative gradient operator but it is a nonconstant coefficients differential operator of the form T=∑ℓ=13eℓaℓ(x)∂xℓ,x=(x1,x2,x3)∈Ω¯,where, Ω can be either a bounded or an unbounded domain in R3 whose boundary ∂Ω is considered suitably regular, Ω ¯ is the closure of Ω and eℓ, for ℓ= 1 , 2 , 3 are the imaginary units of the quaternions H. The operators Tℓ:=aℓ(x)∂xℓ, for ℓ= 1 , 2 , 3 , are called the components of T and a1, a2, a3: Ω ¯ ⊂ R3→ R a...