Another way to reach the mutual inductance is through the vector potential \( \mathbf{A} \) of the solenoid, using \( \Phi = \oint \mathbf{A} \cdot d\mathbf{l} \) instead of directly multiplying field by area.
Step 1: Vector potential of a long solenoid.
For a point inside the solenoid, at radial distance \( s \le r_1 \) from the axis, the vector potential circles the axis with magnitude \[A(s) = \frac{\mu_0 N_1 I_1}{2L} \, s\] obtained from \( \oint \mathbf{A} \cdot d\mathbf{l} = \Phi_{\text{enclosed}} \) applied to a circle of radius \( s \), enclosing flux \( B \pi s^2 \). For a point outside the winding (\( s > r_1 \)), all of the solenoid's flux \( B \pi r_1^2 \) is enclosed no matter how far out the circle is drawn, so \[A(s) = \frac{\mu_0 N_1 I_1}{2L} \, \frac{r_1^2}{s}\]
Step 2: Apply this at the outer coil.
The outer coil sits at radius \( r_2 > r_1 \), outside the solenoid winding, so \[A(r_2) = \frac{\mu_0 N_1 I_1}{2L} \, \frac{r_1^2}{r_2}\]
Step 3: Flux linked with one turn of the coil.\[\phi = \oint \mathbf{A} \cdot d\mathbf{l} = A(r_2) \cdot 2\pi r_2 = \frac{\mu_0 N_1 I_1}{2L} \cdot \frac{r_1^2}{r_2} \cdot 2\pi r_2\]The factor of \( r_2 \) cancels exactly:\[\phi = \mu_0 \frac{N_1 I_1}{L} \pi r_1^2\]the same result as before, but now it is clear why the coil's own radius \( r_2 \) drops out: the vector potential outside the solenoid falls off as \( 1/r_2 \) exactly fast enough to cancel the \( 2\pi r_2 \) circumference.
Step 4: Total flux and mutual inductance.
With \( N_2 \) turns on the outer coil,\[\Phi = N_2 \phi = \mu_0 \frac{N_1 N_2}{L} \pi r_1^2 I_1, \qquad M = \frac{\Phi}{I_1} = \mu_0 \frac{N_1 N_2 \pi r_1^2}{L}\]
On whether \( M_{12} = M_{21} \) holds: this identity does not depend on which circuit carries the current, it is a general consequence of the reciprocal nature of magnetic flux linkage between any two fixed loops in space. So yes, \( M_{12} = M_{21} \) is valid for this solenoid-coil pair, exactly as for any two coupled circuits.