Let xy lie in rad(H), while x does not lie in rad(H). This means (xy)n lies in H, but no power of x lies in H. Seen another way, xnyn lies in H, but xn does not, hence some power of yn lies in H, a power of y lies in H, and y lies in rad(H). Therefore rad(H) is a prime ideal.
Though H has its own prime radical, a prime ideal may be the radical of several different primary ideals. For example the primary ideals generated by Pn in Z all have the same prime radical P.
Assume P = rad(H) and H is primary. Now if xy is in H and x is not in H then yn is in H, and y is in P. We normally say xy in P means at least one factor is in P, but now, if xy is in H at least one factor is in P. Think of H as a concentrator for P.
Conversely, assume P = rad(H), and let H be a concentrator for P. If x is not in H then y is in P is in rad(H), yn lies in H, and H is primary.
When H is primary and P = rad(H) we say H is P primary. Here is a simple theorem about M primary ideals, where M is maximal.
Since M = rad(H), H is M primary.
As a corollary, all powers of M are M primary. We proved this for a pid, but it holds in any ring.
Let J contain xy, but not x. Thus there is some primary ideal Hi that contains xy, but not x. Since Hi contains a power of y, P contains y. Since rad(J) = P, J contains a power of y, and J is P primary.
Let R = Z[x], and let 2 and x generate P, while 4 and x generate Q. Since R/P = Z2, P is maximal. Since R/Q = Z4, Q is primary. The image of P in R/Q gives the nil radical, so P = rad(Q), and Q is P primary. Yet P2, generated by 4, 2x, and x2, lies properly inside Q, hence Q is not a power of P.
Let R = K[x,y,z] mod z2-xy, and let x and z generate P. The quotient R/P is K[y], an integral domain, so P is indeed prime. The product xy lies in P2, being equal to z2. However, x is not in P2, and no power of y lies in P2, hence P2 is not primary.