We introduce a prime density triangle, a multiplicative array constructed by a modified rule of indices. Two display formats are given: a right-angle triangle with converted and unconverted entries, and an equilateral triangle obtained by mirroring the right-angle triangle about the central term xn. We derive the row-sum generating function and provide worked examples using both direct computation and the generating function. Using a four-rule division scheme applied to the unconverted rows we obtain the sequence whose integer parts sum to the magic numbers for two spin orientations. Replacing the repeated end integers by 1,1 yields whose integer parts give the standard harmonic oscillator,HO and in conjunction with the 2 spin orientation nucleon magic numbers gives precisely the standard nuclear magic numbers 2, 8, 20, 28, 50, 82, 126, 184, . . . Subshell filling up to x9, deductions, and evidence from nuclear physics are included. A recurrence M = m ± [(n1n2) + 2] is presented for generating further magic numbers. The sequences appear as OEIS A005897 and OEIS A018226 and OEIS A007290