The [paper](https://www.sciencedirect.com/science/article/pii/S0012365X07007960) "On the arithmetic product of combinatorial species" introduces a combinatorial interpretation for the arithmetic product of two combinatorial species $M,N$. 

In short, let $F:Bij\to Set$ be a (finite) combinatorial species and define its (modified) Dirichlet series to be
$$ D_F(s) = \sum_{n=1}^\infty\frac{|F[n]|}{n!n^s}$$ where $|F[n]|$ is the cardinality of the set $F[n]$. Then, the arithmetic product $\boxdot$ is a monoidal operation on the category of species such that $D_{M\boxdot N}(s)=D_M(s)D_N(s)$.

Explicitly, one takes an (exponential) power series $\sum \frac{a_n}{n!}X^n$ for the species $M$, and one $\sum \frac{b_n}{n!}X^n$ for $N$, and their arithmetic product is 
$$ \sum \frac{c_n}{n!}X^n \qquad c_n := \sum_{d\mid n} \begin{Bmatrix}n\\d\end{Bmatrix}a_db_{n/d}$$ (there is an unfortunate clash of notation between this notation and Stirling numbers of the 2nd kind here, but $\left\{\begin{smallmatrix}n\\ d\end{smallmatrix}\right\}$ is just the "modified binomial" $\frac{n!}{d!(n/d)!}$). 

With sparse remarks here and there, the paper proves what boils down to the fact that the category of species has a symmetric monoidal structure given by $\boxdot$. The only thing is, usually Dirichlet series and their convolution are not defined in this way, but instead without the $\frac 1{n!}$. At any rate, it seems to me that with the appropriate substitution/transport of structure the monoidal category $(Species,+,\boxdot)$ can be regarded as a categorified version of the semiring of Dirichlet series with nonnegative integer coefficients.

My question is: is there a description of this monoidal structure that looks concise to a category theorist, or at least as concise as all other monoidal structures (Cartesian, Day convolution, substitution...) do? What is the universal property, if any, of the $\boxdot$ structure? (Compare this question with: the substitution product of species is the unique monoidal structure making $Species$ monoidally equivalent to the category of analytic functors + composition, and it's called "substitution product" because it coincides with the composition of formal power series $(\sum a_n t^n)\circ g = \sum a_n g(t)^n$).

My understanding is that this structure is *related* to the substitution product (obtained as iterated convolution) which is what the authors (and other people studying species with more combinatorial jargon) call "assemblies": indeed, they say: 
> The most interesting combinatorial construction associated to the arithmetic product is the assembly of cloned
structures. Informally, an assembly of cloned G-structures is an assembly of G-structures in the above sense, where
all structures in the assembly are isomorphic replicas of the same structure. Moreover, information about ‘homologous
vertices’ or ‘genetic similarity’ between each pair in the assembly is also provided. The structures of F ⊡ G have
some resemblance with the structures of the substitution F(G). An element of F ⊡ G can be represented as a cloned
assembly of G-structures together with an external F-structure (an F-assembly of cloned G-structures). Because
of the symmetry F ⊡ G = G ⊡ F this structures can also be represented as G-assemblies of cloned F-structures.
For example, if L+ denotes the species of non-empty lists, the structures of F ⊡ L+ could be thought of either as
F-assemblies of cloned lists, or as lists of cloned F-structures.

But I'm having a really hard time even understanding what is the notion of a "rectangle" that the authors employ to give to $\boxdot$ a combinatorial interpretation. I reason way better with coends and Kan extensions...