\section{Introduction}

Since its first outbreak in Wuhan (China) in December 2019, the coronavirus disease (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has spread rapidly globally.
On March 2, 2020, the World Health Organization (WHO) declared COVID-19 a pandemic~\footnote{https://www.who.int/director-general/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19---11-march-2020}.
Even in 2022, the COVID-19 pandemic shows no signs of cessation.
In response to the COVID-19 pandemic, various public health and social measures have been implemented globally~\cite{estrada2020covid}.
The large-scale preventative measures implemented, such as the imposition of travel restrictions, closure of schools and businesses, prohibition of group gatherings, and implementation of lockdowns, have severely disrupted social and economic activities.
To limit the spread of the virus, washing of hands, social distancing, and mask usage have become the new prerequisites for partaking in social activities.
Thus, it is important to investigate the effectiveness of mask use in suppressing the virus transmission in various situations.

To understand the spread of diseases beyond the population level, it is necessary to consider a contact network of individuals~\cite{barabasi2016network,newman2018networks}.
Numerous studies have shown that the structure of an underlying contact network significantly impacts the propagation of epidemics (see review~\cite{pastor2015epidemic} and references therein).
For instance, an epidemic can propagate even with a small infection rate if the underlying network is highly heterogeneous~\cite{pastor2001epidemic,moreno2002epidemic}.
It should further be noted that population heterogeneity in terms of socio-demographic and health factors also affects disease dynamics~\cite{markovivc2021socio}.
Network scientists have explored network-based strategies to efficiently suppress epidemics, including optimal vaccine allocation~\cite{gallos2007improving, wang2016statistical}, rapid isolation of infected individuals~\cite{takaguchi2014suppressing,hasegawa2017efficiency,strona2018rapid}, social distancing~\cite{gutin2021effect,maharaj2012controlling,valdez2012intermittent,maheshwari2020network,nande2021dynamics}, and community lockdown~\cite{gosak2021community}.
Recent studies (since COVID-19) have discussed the effectiveness of mask use in limiting the spread of the virus in networks~\cite{li2020mask,eikenberry2020mask,leung2020contact,yagan2021modeling,tian2021analysis,ying2021modelling,sridhar2021leveraging,lee2021epidemic,tian2021role}.
Ya{\u{g}}an et al.~\cite{yagan2021modeling} proposed the implementation of a susceptible-infected-recovered (SIR) model in a network consisting of mask wearers and non-mask wearers, based on a multiple-strain model, where a virus may mutate into a different strain~\cite{alexander2010risk,eletreby2020effects}.
They referred to it as the {\it mask model}.
Tian et al. derived the expected size and occurrence probability of large outbreaks in the mask model~\cite{tian2021analysis}.
Additionally, Lee and Zhu considered the same model to estimate the infection and recovery rates from John Hopkins University's COVID-19 database.
They further assessed the mask coverage using the maximum likelihood estimation method~\cite{lee2021epidemic}.

Homophily, the tendency for individuals to connect with similar others, has been found in a wide variety of attributes such as race, age, and gender~\cite{mcpherson2001birds}.
An experimental study by Centola~\cite{centola2011experimental} discerned homophily in the adoption of health behaviors.
Since COVID-19, some studies have reported homophily in mask use~\cite{haischer2020wearing,woodcock2021role}.
Despite the guidelines from the United States Centers for Disease Control and Prevention~\footnote{https://www.cdc.gov/coronavirus/2019-ncov/prevent-getting-sick/diy-cloth-face-coverings.html} and the WHO~\footnote{https://www.who.int/emergencies/diseases/novel-coronavirus-2019/advice-for-public/when-and-how-to-use-masks}, some members of the population do not wear masks, or even if they do, they do not wear highly protective masks, such as surgical masks and N95 masks.
Homophily in mask use will be factored in as {\it assortative mixing}~\cite{newman2003mixing} in networks.
Assortative mixing here means a tendency for (non-)mask wearers to prefer to connect with other (non-)mask wearers and disassortative mixing means a tendency for mask-wearers to prefer to connect with non-mask wearers.
Whereas previous studies~\cite{yagan2021modeling,tian2021analysis,sridhar2021leveraging} assumed that mask wearers and non-mask wearers are randomly distributed, the present study considers mixing pattern by mask-wearing to discuss the impact of assortative mixing on the effectiveness of mask use in suppressing the propagation of epidemics.

In this paper, we analyze the mask model in assortatively mixed random networks that follow an arbitrary degree distribution.
We derive the occurrence probability of large outbreaks, mean size of large outbreaks, and average epidemic size of the mask model in assortatively mixed random networks using probability generating functions.
Applying our analysis to the mask model in Poisson random networks, we establish that the assortative (disassortative) mixing by mask-wearing decreases (increases) the epidemic threshold.
It means that disassortative mixing is effective in containing epidemics in that a high transmissibility is required for large outbreaks.
On the other hand, the result for the case of high transmissibility shows that the mask use is most effective in reducing the probability and mean size of large outbreaks when the mixing pattern is strongly assortative.
The separation of mask wearers and non-mask wearers, which occurs in strongly assortative networks, reduces the probability and degree of a large outbreak in high-transmissibility cases, although it allows a large outbreak to occur even in low-transmissibility cases.
For scale-free networks, we establish that mask use is most effective when the mixing pattern is strongly assortative and maximally disassortative (provided that the mask coverage is not extremely low).
We also performed extensive simulations to confirm that our analytical results are in full agreement with the simulation results.

The paper is organized as follows.
In Sec.~\ref{sec:model}, we introduce the mask model, which is a generalization of a discrete-time SIR model, and assortatively mixed random networks proposed in~\cite{newman2003mixing}.
In Sec.~\ref{sec:analysis}, we derive the occurrence probability and mean size of large outbreaks, as well as the average epidemic size, for the mask model in assortatively mixed random networks with an arbitrary degree distribution and mixing matrix.
In Sec.~\ref{sec:results}, we discuss our results for the mask model in two representative networks: a Poisson network and a scale-free network.
A comparison between the analytical treatment and simulation results is shown.
Section~\ref{sec:summary} is devoted to the summary.

\section{Model \label{sec:model}}

In this study, we treat a discrete-time SIR-type model as follows~\cite{newman2002spread}.
Let us consider a network of $n$ nodes and edges, where the nodes represent individuals and the edges represent contact between them.
We assume that each node can be in one of three states: susceptible, infected, or recovered.
The transmission of infections to a susceptible node occurs, with a constant probability, when the latter is adjacent to an infected node.
The system evolves in discrete times as follows.
The dynamics start with an initial state where a single seed node is randomly selected and is set to infected, while the other nodes are set to susceptible.
The seed node infects each of its adjacent nodes independently with probability $T$.
Here $T$ is called the {\it transmissibility} of the disease~\cite{newman2002spread}.
Thereafter, the seed node spontaneously recovers (recovered node), i.e., it can neither infect nor be infected.
In a next step, a newly infected node infects each of its adjacent susceptible nodes with probability $T$ and subsequently recovers.
This process continues until no newly infected node appears (or equivalently, only susceptible and recovered nodes exist in the network).

The scale of outbreaks in a network depends entirely on whether the transmissibility, $T$, exceeds the epidemic threshold, $T_c$, or not.
For $T<T_c$, epidemics starting from a single seed node always terminate at an early stage and the fraction of nodes eventually infected is negligibly small.
For $T>T_c$, a single infected seed node can cause a large outbreak (an outbreak of finite relative size in a large network) with nonzero probability.
The average epidemic size $s_{\rm av}$, defined as the average fraction of recovered nodes in the final state, is employed to characterize this phase transition: in the limit $n \to \infty$, $s_{\rm av}=0$ when $T \le T_c$ and $s_{\rm av}>0$ when $T>T_c$.
It is known~\cite{newman2002spread} that the epidemic threshold of the SIR model in random networks with degree distribution $p(k)$ is
\begin{equation}
T_c = \frac{\langle k \rangle}{\langle k^2 \rangle -\langle k \rangle}, \label{eq:knownThreshold}
\end{equation}
where $\langle A(k) \rangle = \sum_k A(k) p(k)$.

The mask model was proposed as a variant of the above epidemic model, where some individuals wear masks, while others do not~\cite{yagan2021modeling,tian2021analysis}.
Each node in the mask model is classified into one of the following two types: type-1 nodes representing mask wearers and type-2 nodes representing non-mask wearers.
We denote the mask coverage by $p_1$ and assume that the $p_1$ fraction of randomly chosen nodes is of type 1, while the other $p_2(=1-p_1)$ fraction is of type 2.
In the mask model, there are four types of contacts between infected and susceptible nodes, depending on whether each of them wears a mask or not.
We designate $T_{xy}$ as the transmissibility from an infected node of type $x$ ($x=1,2$) to a susceptible node of type $y$ ($y=1,2$).
As the previous study~\cite{yagan2021modeling} pointed out, mask use is effective in reducing the droplets expelled from an infected individual, but it is less effective in preventing an individual from inhaling droplets from an infected non-mask wearer~\cite{brienen2010effect, lai2012effectiveness, davies2013testing,ueki2020effectiveness}; thus, we suppose that $T_{11}<T_{12}<T_{21}<T_{22}$.

The objective of this study is to assess the impact of assortative mixing by mask-wearing on the effectiveness of mask use in suppressing the propagation of epidemics.
The mixing pattern is characterized by mixing matrix ${\bm e}$, whose element $e_{xy}$ denotes the fraction of edges connecting a type-$x$ node and a type-$y$ node~\cite{newman2003mixing}.
For a network in which nodes can be one of two types (wearing or non-wearing group), ${\bm e}=(e_{xy})$ satisfies:
\begin{equation}
e_{xy}=e_{yx},\; p_x = \sum_{y=1,2}e_{xy}=\sum_{y=1,2}e_{yx}, \; \sum_{x,y=1,2} e_{xy} = 1.
\end{equation}
The tendency of (non-)mask wearers to be connected to other (non-)mask wearers is quantified by the {\it assortativity coefficient} $r$~\cite{newman2003mixing}.
For the present case, the assortativity coefficient of mask-wearing is determined only by $p_1$ ($0 < p_1 < 1$) and $e_{12}$ ($0 \le e_{12} \le {\rm min}\{p_1,1-p_1\}$):
\begin{equation}
r
=\frac{\sum_{x} e_{xx} - \sum_x p_x^2}{1-\sum_x p_x^2}
=\frac{p_1(1-p_1) -e_{12}}{p_1(1-p_1)},
\end{equation}
where the range of $r$ is
\begin{equation}
r_{\rm min} \le r \le 1,
\end{equation}
$r_{\rm min}$ being
\begin{equation}
r_{\rm min} = \begin{cases}
-p_1/(1-p_1) & (0 < p_1 \le 1/2) \\
-(1-p_1)/p_1 & (1/2 \le p_1 < 1)
\end{cases}.
\label{eq:rmin}
\end{equation}
The network is {\it assortative}, meaning that mask-wearers are likely to be connected to other mask-wearers, {\it neutral}, and {\it disassortative}, meaning that mask-wearers are unlikely to be connected to other mask-wearers, when $r>0$, $r \approx 0$, and $r<0$, respectively.

We introduce a procedure~\cite{newman2003mixing} for generating assortatively mixed and random networks with an arbitrary degree distribution $p(k)$ and mixing matrix ${\bm e}$.
Let the number of edges $m$, degree distribution $p(k)$, and mixing matrix ${\bm e}$ be predetermined.
The procedure is as follows:
(i) we prepare $m$ edges so that the type of ends of each edge follow ${\bm e}$.
Later, $m e_{11}$ edges connect two type-$1$ nodes, $m(e_{12}+e_{21})$ edges connect a type-$1$ node and a type-$2$ node, and $m e_{22}$ edges connect two type-$2$ nodes.
The total degrees of type-1 nodes, $m_{1}^{(\rm e)}$, and type-2 nodes, $m_{2}^{(\rm e)}$, are expressed as $m_{1}^{(\rm e)}=m (2 e_{11} + e_{12})$ and $m_{2}^{(\rm e)}=m (e_{21} + 2 e_{22})$, respectively.
Subsequently, we determine the number, $n_x$, of type-$x$ nodes by $n_x=m_x^{(\rm e)}/\langle k \rangle$ ($x=1,2$).
(ii)
For each type $x$, we assign $k_i$ stubs to each node $i$ according to the distribution $p(k)$.
We select a type-$x$ node at random and reassign stubs to it from $p(k)$, repeating the process until the total degree of type-$x$ nodes is equal to $m_x^{(\rm e)}$.
(iii)
We repeatedly join a randomly selected pair of stubs to form an edge, thereby realizing a network.
Here, multiple edges and self-loops are prohibited.
By implementing this procedure while tuning $p_1$ and $e_{12}$, we obtain assortatively mixed random networks with a given distribution $p(k)$ and level of assortative mixing $r$.

Figure~\ref{fig:Snapshot} illustrates assortatively mixed random networks with $r = -0.95$, $0.0$, and $0.95$.
The more $r$ is negative, the number of edges connecting a mask wearer and a non-mask wearer increases and the network becomes closer to bipartite.
The more positive $r$ is, the number of edges connecting a mask wearer and a non-mask wearer decreases and the separation of groups of mask-wearers and non-mask-wearers becomes more pronounced.
In particular, the perfectly disassortative case ($r=-1$) results in a bipartite network of mask wearers and non-mask wearers, and the perfectly assortative case ($r=1$) results in two disconnected random networks, one consisting of mask wearers and the other consisting of non-mask wearers.
The case of $r=0$ corresponds to the mask model on random networks analyzed in~\cite{yagan2021modeling,tian2021analysis}.

 \begin{figure}[tb]
 \centering
  \includegraphics[width=12cm]{./figures/fig1.eps}
 \caption{
Illustration of a disassortative network ($r = -0.95$), a neutral network ($r = 0.00$), and an assortative network ($r = 0.95$), from left to right.
The mask coverage $p_1$ is set to $p_1=0.5$.
All networks have $n \approx 200$ nodes and $m = 400$ edges, and follow a Poisson degree distribution with $\langle k \rangle \approx 4$.
The red nodes and green nodes represent mask wearers and non-mask wearers, respectively.
The more $r$ is negative, the number of edges connecting a mask wearer and a non-mask wearer increases and the network becomes closer to bipartite.
The more positive $r$ is, the number of edges connecting a mask wearer and a non-mask wearer decreases and the separation of groups of mask-wearers and non-mask-wearers becomes more pronounced.
 }
 \label{fig:Snapshot}
 \end{figure}

\section{Analysis \label{sec:analysis}}

The mask model can be mapped onto the heterogeneous bond percolation of a multitype network~\cite{allard2009heterogeneous}.
In Ref.~\cite{allard2009heterogeneous}, Allard et al. derived a number of percolation properties for heterogeneous bond percolation on multitype networks with arbitrary degree distribution, using a generating function formalism.
In the context of epidemiology, their generating functions yield the expected size of small outbreaks, the mean size of large outbreaks, the occurrence probability of large outbreaks, and the critical transmissibility set, over which a seed node leads to a large outbreak with nonzero probability.
Similar to Ref.~\cite{allard2009heterogeneous}, we obtain the relative size of the giant component, $s_{\rm L}$, of heterogeneous bond percolation on assortatively mixed random networks.
This giant component size $s_{\rm L}$ corresponds to the mean size of large outbreaks.
Additionally, we derive the occurrence probability of large outbreaks $P_{\rm L}$ using the probability generating functions.
Thereafter, the average epidemic size $s_{\rm av}$ is given by $P_{\rm L}$ and $s_{\rm L}$.
It should be noted that in heterogeneous bond percolation, the mean size of large outbreaks $s_{\rm L}$, given as the probability that a randomly chosen node is reached from an infinite connected component does not always coincide with the occurrence probability $P_{\rm L}$, given as the probability that a randomly chosen node does not link only to finite components~\cite{allard2009heterogeneous}.
The difference between $s_{\rm L}$ and $P_{\rm L}$ has been known for the SIR model on semidirected networks~\cite{meyers2006predicting} and the SIR model with heterogeneity of infectious periods, infectivity, and susceptibility~\cite{kenah2007second, miller2007epidemic} prior to~\cite{allard2009heterogeneous}.

Let us consider an infinitely large, assortatively mixed, and random network with degree distribution $p(k)$ and mixing matrix ${\bm e}$.
A node reached by following a randomly chosen edge (or, a randomly chosen neighbor of a node in a random network) has $k$ other edges with probability $q(k)=(k+1)p(k+1)/\langle k \rangle$, where $q(k)$ represents the excess degree distribution.
The generating functions for $p(k)$ and $q(k)$ are defined as
\begin{equation}
G_p(x)=\sum_k p(k) x^k
\quad {\rm and} \quad
G_q(x)=\sum_k q(k) x^k,
\end{equation}
respectively.

Now, we consider the heterogeneous bond percolation~\cite{allard2009heterogeneous} with probability $\{T_{xy}\}$ ($T_{xy}$ is the occupation probability of an edge from a type-$x$ node to a type-$y$ node).
We denote by $u_x$ the probability that a type-$x$ node ($x=1,2$) does not receive the infection (occupied edge) from a neighbor in the giant component.
The probability, $p_{xy}$, that a randomly chosen neighbor of a type-$x$ node is of type $y$ is given by $p_{xy} = e_{xy}/p_x$.
Noting that a type-$y$ neighbor belongs to the giant component is $1-\sum_k q(k) u_y^k = 1 - G_q(u_y)$, the probability that a neighbor of a type-$x$ node is of type $y$, belongs to the giant component, and infects the type-$x$ node is given by $p_{xy} T_{yx} (1-G_q(u_y))$.
Thus, in the infinitely large, assortatively mixed, random network with ${\bm e}$, $u_1$ and $u_2$ obey the following equations:
\begin{subequations}
\begin{align}
u_1 &= 1 - p_{11} T_{11}(1-G_q (u_1)) - p_{12} T_{21}(1-G_q (u_2)), \label{eq:u1}\\
u_2 &= 1 - p_{21} T_{12}(1-G_q (u_1)) - p_{22} T_{22}(1-G_q (u_2)). \label{eq:u2}
\end{align} \label{eq:u}
\end{subequations}
A randomly chosen node is of type $x$ with probability $p_x$, and it does not belong to the giant component with probability $\sum_k p(k) u_x^k = G_p(u_x)$.
The relative size of the giant component $s_{\rm L}$ is one minus the probability that a randomly chosen node is not connected to the giant component:
\begin{equation}
s_{\rm L} = 1- p_1 G_p(u_1) - p_2 G_p(u_2), \label{eq:S}
\end{equation}
where $u_1$ and $u_2$ are given as the solutions of Eq.~(\ref{eq:u}).
Since this giant component can be regarded as the set of nodes that will eventually be infected in an epidemic if a large outbreak occurs, $s_{\rm L}$ corresponds to the mean size of large outbreaks (the expected size of an outbreak conditioned on it being a large one).

Furthermore, we derive the critical transmissibility set from Eq.~(\ref{eq:u}).
For this, we consider the stability of the trivial solution of Eq.~(\ref{eq:u}), i.e., $(u_1, u_2)=(1,1)$, assuming that $u_1 \approx 1-\epsilon_1$ and $u_2 \approx 1-\epsilon_2$.
By expanding Eq.~(\ref{eq:u}) to the leading order in $\epsilon_1$ and $\epsilon_2$, we have
\begin{equation}
{\bm \epsilon} = A {\bm \epsilon},
\end{equation}
where
\begin{equation}
{\bm \epsilon}=
\begin{pmatrix}
\epsilon_1 \\
\epsilon_2
\end{pmatrix}
\quad {\rm and} \quad
A=
G_q' (1)
\begin{pmatrix}
p_{11} T_{11} & p_{12} T_{21} \\
p_{21} T_{12} & p_{22} T_{22}
\end{pmatrix}.
\end{equation}
Thus, the critical transmissibility set is obtained from ${\rm det} (A-I)=0$; more specifically,
\begin{equation}
(p_{11}T_{11} G_q'(1) -1)(p_{22} T_{22} G_q'(1) -1)-p_{12} p_{21} T_{21}T_{12} G_q'(1)^2 =0. \label{eq:thresholdCondition}
\end{equation}
It should be noted that condition (\ref{eq:thresholdCondition}) reduces to the epidemic threshold (\ref{eq:knownThreshold}) for the SIR model with transmissibility $T$ when mask usage is completely ineffective, i.e., $T_{11}=T_{12}=T_{21}=T_{22}=T$.

The probability $p_{xy}$ is written in terms of $p_1$ and $r$ as $p_{11} = 1-(1-p_1)(1-r)$, $p_{12}=(1-p_1)(1-r)$, $p_{21}=p_1(1-r)$, and $p_{22}=1-p_1(1-r)$.
Equation~(\ref{eq:thresholdCondition}) is then rewritten as
\begin{equation}
((r+p_1(1-r)) T_{11} G_q'(1) -1)((1-p_1(1-r))T_{22} G_q'(1) -1) - p_1(1-p_1)(1-r)^2 G_q'(1)^2 T_{12} T_{21} =0. \label{eq:thresholdCondition2}
\end{equation}
We should mention the following cases.
(i) For the neutral case ($r=0$), Eq.~(\ref{eq:thresholdCondition2}) reduces to
\begin{equation}
(p_1 T_{11} G_q'(1)-1)((1-p_1) T_{22} G_q'(1)-1) - p_1(1-p_1) T_{12} T_{21} G_q'(1)^2 =0,
\end{equation}
which is identical to that of the mask model in random networks already reported in~\cite{tian2021analysis}.
(ii) For the perfectly assortative case ($r=1$), Eq.~(\ref{eq:thresholdCondition2}) reduces to
\begin{equation}
(T_{11} G_q'(1)-1)(T_{22} G_q'(1)-1) = 0. \label{eq:thresholdAssortative}
\end{equation}
In this case, the critical transmissibility set is not affected by $T_{12}$ and $T_{21}$ because no edges exist between mask wearers and non-mask wearers.
Supposing that $T_{11}<T_{22}$ (i.e., mask has an effect in reducing the transmissibility), the epidemic threshold is given by $T_{22} G_q'(1)=0$.
(iii) For the maximally disassortative case ($r=r_{\rm min}$), Eq.~(\ref{eq:thresholdCondition2}) reduces to
\begin{equation}
\begin{cases}
(1-2p_1) T_{22} G_q'(1) - p_1 T_{12}T_{21} G_q'(1)^2 = 1-p_1 & (0<p_1<1/2) \\
(2p_1-1) T_{11} G_q'(1) +(1-p_1) T_{12}T_{21} G_q'(1)^2 = p_1 & (1/2 \le p_1 <1)
\end{cases}
.
\end{equation}
The lack of $T_{11}$ ($T_{22}$) for $p_1<1/2$ ($p_1 \ge 1/2$) reflects that in the maximally disassortative case, all neighbors of mask wearers (non-mask wearers) are non-mask wearers (mask wearers).
The perfectly disassortative case ($r = -1$) is realized only when $p_1=1/2$.
In this case, the above condition becomes $T_{12}T_{21} G_q'(1)^2 =1$, so neither $T_{11}$ nor $T_{22}$ exists.

Next, we derive the probability $P_{\rm L}$ that a single seed node causes a large outbreak.
We designate $P_n^{(x)}(n_1,n_2)$ as the probability that a seed node of type $x$ infects $n_1$ neighbors of type 1 and $n_2$ neighbors of type 2.
This probability, $P_n^{(x)}(n_1,n_2)$, is given by
\begin{equation}
P_n^{(x)}(n_1,n_2) = \sum_{k_1=n_1}^\infty \sum_{k_2=n_2}^\infty p(k_1+k_2) \binom{k_1+k_2}{k_1} \prod_{y = 1,2} p_{xy}^{k_y}
\binom{k_y}{n_y} T_{xy}^{n_y} (1-T_{xy})^{k_y-n_y}.
\end{equation}
By introducing the generating function for $P_n^{(x)}(n_1,n_2)$ as $F_p^{(x)}(y,z) = \sum_{n_1} \sum_{n_2} P_n^{(x)} (n_1,n_2) y^{n_1} z^{n_2}$, we have
\begin{subequations}
\begin{align}
F_p^{(1)}(y,z) & = G_p(1+p_{11}T_{11}(y-1)+p_{12}T_{12}(z-1)), \\
F_p^{(2)}(y,z) & = G_p(1+p_{21}T_{21}(y-1)+p_{22}T_{22}(z-1)).
\end{align}
\end{subequations}
Likewise, we denote by $Q_n^{(x)}(n_1,n_2)$ the probability that a type-$x$ infected node that is reached by following a randomly chosen edge and infected through it newly infects $n_1$ type-1 neighbors and $n_2$ type-2 neighbors.
The generating function for $Q_n^{(x)}(n_1,n_2)$ is $F_q^{(x)}(y,z)= \sum_{n_1} \sum_{n_2} Q_n^{(x)} (n_1,n_2) y^{n_1} z^{n_2}$, which is written as
\begin{subequations}
\begin{align}
F_q^{(1)}(y,z) &= G_q(1+p_{11}T_{11}(y-1)+p_{12}T_{12}(z-1)), \\
F_q^{(2)}(y,z) &= G_q(1+p_{21}T_{21}(y-1)+p_{22}T_{22}(z-1)).
\end{align}
\end{subequations}
We further introduce the probability, $P_{s}^{(x)}(s_1,s_2)$, that a seed node of type $x$ causes an outbreak, where $s_1$ type-1 nodes and $s_2$ type-2 nodes are eventually infected, and the probability, $Q_{s}^{(x)}(s_1,s_2)$, that the infection of a type-$x$ node reached by following a randomly chosen edge results in the secondary infection of $s_1$ type-1 nodes and $s_2$ type-2 nodes.
We define the generating function for $P_{s}^{(x)}(s_1,s_2)$ and $Q_{s}^{(x)}(s_1,s_2)$ as $H_{p}^{(x)}(y,z) = \sum_{s_1} \sum_{s_2} P_s^{(x)} (s_1,s_2) y^{s_1} z^{s_2}$ and $H_{q}^{(x)}(y,z) = \sum_{s_1} \sum_{s_2} Q_s^{(x)} (s_1,s_2) y^{s_1} z^{s_2}$, respectively.
These generating functions obey the following equations:
\begin{subequations}
\begin{align}
H_{p}^{(1)}(y,z) &= y F_p^{(1)} (H_q^{(1)}(y,z), H_q^{(2)}(y,z)),
\quad
H_{p}^{(2)}(y,z) = z F_p^{(2)} (H_q^{(1)}(y,z), H_q^{(2)}(y,z)),
\\
H_{q}^{(1)}(y,z) &= y F_q^{(1)} (H_q^{(1)}(y,z), H_q^{(2)}(y,z)),
\quad
H_{q}^{(2)}(y,z) = z F_q^{(2)} (H_q^{(1)}(y,z), H_q^{(2)}(y,z)).
\end{align}
\end{subequations}
A randomly chosen node is of type $x$ with probability $p_x$, and it causes an epidemic with finite infections with probability $\sum_{s_1} \sum_{s_2} P_s^{(x)}(s_1,s_2)=H_p^{(x)}(1,1)$.
Thus, the occurrence probability of a large outbreak $P_{\rm L}$, which is the probability that an epidemic starting from a single seed node does not terminate with finite infections, is
\begin{equation}
P_{\rm L}
= 1- p_1 H_p^{(1)}(1,1) - p_2 H_p^{(2)}(1,1)
= 1- p_1 F_p^{(1)}(v_1,v_2) - p_2 F_p^{(2)}(v_1,v_2), \label{eq:pG}
\end{equation}
where $v_1= H_q^{(1)}(1,1)$ and $v_2 = H_q^{(2)}(1,1)$ are the solutions of the following equations:
\begin{equation}
v_1 = F_q^{(1)}(v_1,v_2),
\quad
v_2 = F_q^{(2)}(v_1,v_2).
\end{equation}

It should be stated that the critical transmissibility set is related to the reproductive number $R_0$, which is the average number of infections caused by a newly infected individual.
The critical transmissibility set is given as $R_0=1$ in that an epidemic persists for $R_0>1$ and dies out for $R_0<1$.
Following \cite{tian2021analysis,tian2021role}, $R_0$ is defined as as the spectral radius of matrix $B = (B_{xy})$,
\begin{equation}
R_0 = \rho(B),
\end{equation}
where $B_{xy}$ is the expected number of type-$y$ nodes infected by a newly infected type-$x$ node.
Since the expected number of neighbors of a newly infected node is $G_q'(1)$ and each neighbor is of type $y$ with probability $p_{xy}$ and will be infected with probability $T_{xy}$,
we have
\begin{equation}
B_{xy} = G_q'(1) T_{xy} p_{xy}. \label{eq:mat}
\end{equation}
Substituting Eq.~(\ref{eq:mat}) into $R_0 = \rho(B) = 1$, we obtain the critical transmissibility set (\ref{eq:thresholdCondition}).

Finally, we determine the average epidemic size $s_{\rm av}$ using $s_{\rm L}$ and $P_{\rm L}$.
An infection starting from a single seed node evolves into a large epidemic outbreak of expected size $s_{\rm L}$ with probability $P_{\rm L}$; otherwise, the epidemic propagation becomes small so that the final fraction of recovered nodes is negligible.
Thus, the average epidemic size $s_{\rm av}$ is
\begin{equation}
s_{\rm av} = P_{\rm L} s_{\rm L}
=(1- p_1 F_p^{(1)}(v_1,v_2) - p_2 F_p^{(2)}(v_1,v_2))(1- p_1 G_p(u_1) - p_2 G_p(u_2)). \label{eq:sAverage}
\end{equation}

\section{Results \label{sec:results}}

To discuss the effect of assortative mixing by mask-wearing, we consider the mask model in two typical random networks with tunable assortative mixing: a Poisson network and a scale-free network.
With reference to an experiment~\cite{ueki2020effectiveness}, which measured the effectiveness of mask use in filtering expelled and inhaled droplets, the present study employs the following parameters for $T_{xy}$,
\begin{equation}
\begin{pmatrix}
T_{11} & T_{12} \\
T_{21} & T_{22} \\
\end{pmatrix}
=
T
\begin{pmatrix}
0.24 & 0.42 \\
0.50 & 1.00 \\
\end{pmatrix}. \label{eq:maskEffect}
\end{equation}
We present analytical results as well as Monte Carlo (MC) simulation results.
In the simulations, we generated $10^2$ network realizations of $n = 10^5$ nodes, and implemented the mask model with random placements of a single seed and non-mask wearers $10^3$ times on each realization.
In each run, the size of an epidemic was given by the fraction of recovered nodes in the final state.
Afterward, the average epidemic size $s_{\rm av}$ was given as the epidemic size averaged over all samples.
In the present simulation, we state that a large outbreak occurs if the epidemic size exceeds 0.01.
We regard the frequency of samples in which a large outbreak occurred as the occurrence probability of large outbreaks $P_{\rm L}$, and the epidemic size averaged over those samples in which a large outbreak occurred as the mean size of large outbreaks $s_{\rm L}$.

 \begin{figure}[tb]
 \centering
 (a)
  \includegraphics[width=4.75cm]{./figures/fig2a.eps}
 (b)
  \includegraphics[width=4.75cm]{./figures/fig2b.eps}
 (c)
  \includegraphics[width=4.75cm]{./figures/fig2c.eps}
 \caption{
Results of the mask model on the Poisson network with no assortative mixing ($r=0$).
(a) Occurrence probability of large outbreaks $P_{\rm L}$ (black), mean size of large outbreaks $s_{\rm L}$ (red), and average epidemic size $s_{\rm av}$ (blue).
In panel (a), the mask coverage is fixed at $p_1=0.5$.
The symbols represent the results of MC simulations of $n=10^5$ nodes, and the lines represent the analytical results obtained from Eq.~(\ref{eq:pG}),~(\ref{eq:S}), and~(\ref{eq:sAverage}) for $P_{\rm L}$, $s_{\rm L}$, and $s_{\rm av}$, respectively.
The analytical lines are in excellent agreement with the simulation results.
(b) The cumulative epidemic size distribution, $\Pi(s)=\sum_{s'=s}^\infty \pi(s')$, obtained from MC simulations with $p_1=0.5$ and with $T=0.3114<T_c$ (red-dashed line), $T=0.4114 \approx T_c$ (black-solid line), and $T=0.5114>T_c$ (blue-dotted line).
Here $\pi(s)$ denotes the ratio of the samples whose epidemic size is $s$ in $10^2 \times 10^3$ samples.
Panel (b) indicates that $\Pi(s)$ quickly decreases with $s$ for $T<T_c$, $\Pi(s) \propto s^{2-\tau}$ at $T=T_c$, and $\Pi(s)$ shows the bimodality of $\pi(s)$ for $T>T_c$.
(c) Impact of the mask coverage ($p_1=0.0,\, 0.25,\, 0.5,\, 0.75$, and $1.0$) on the average epidemic size $s_{\rm av}$.
The symbols represent the simulation results, and the lines represent the analytical results based on Eq.~(\ref{eq:sAverage}).
The increase in mask coverage decreases the average epidemic size and increases the epidemic threshold.
 }
 \label{fig:Examples}
 \end{figure}

First, we treat the mask model on the Poisson network whose degree distribution is $p(k)=\langle k \rangle^k e^{-\langle k \rangle}/k!$ and average degree is $\langle k \rangle=4$.
Substituting Eq.~(\ref{eq:maskEffect}) and $G_q'(1)=\langle k^2 - k \rangle/\langle k \rangle=\langle k \rangle$ into Eq.~(\ref{eq:thresholdCondition}), it follows that the epidemic threshold $T_c$ is given as the solution of $(0.96(r+p_1(1-r))T_c-1)(4(1-p_1(1-r))T_c-1)-3.36 p_1(1-p_1)(1-r)^2 T_c^2=0$.
We begin by considering the case of no mixing ($r = 0$).
Figure~\ref{fig:Examples} (a) shows the average epidemic size $s_{\rm av}$, occurrence probability of large outbreaks $P_{\rm L}$, and mean size of large outbreaks $s_{\rm L}$ when the mask coverage is fixed at $p_1=0.5$.
The epidemic threshold is $T_c \approx 0.4114$ when $p_1=0.5$ and $r=0$.
Evidently, all the theoretical lines drawn from equations in the previous section are in excellent agreement with the simulation data (symbols).

As shown in Fig.~\ref{fig:Examples} (a), all the $s_{\rm av}$, $s_{\rm L}$, and $P_{\rm L}$ parameters exhibit non-zero values when $T$  is greater than the epidemic threshold $T_c $.
Further, we observe from the distribution of the epidemic sizes that a phase transition occurs at $T=T_c$.
Figure~\ref{fig:Examples} (b) shows the cumulative distribution of epidemic sizes, $\Pi(s)=\sum_{s'=s}^\infty \pi(s')$, where $\pi(s)$ is the frequency of samples in which the epidemic size is $s$, when $T =0.3114 \;(< T_c)$, $T=0.4114 \;(\approx T_c)$, and $T=0.5114 \;(>T_c)$.
The figure indicates that $\Pi(s)$ for $T<T_c$ quickly decreases with $s$;
$\Pi(s)$ for $T>T_c$ stands for the bimodality of $\pi(s)$ with peaks at small $s$ ($s < 10^{-3}$) and at large $s$ ($s > 10^{-1}$) values;
and $\Pi(s)$ for $T=T_c$ follows the power law, $\Pi(s) \propto s^{2-\tau}$.
In the present case, exponent $\tau$ is given as $\tau=5/2$, which coincides with that in the mean-field theory of percolation \cite{stauffer2018introduction} (and the SIR model is known to be in percolation universality class \cite{grassberger1983critical,tome2010critical,de2011new}).

Figure~\ref{fig:Examples} (c) indicates the impact of the mask coverage in suppressing the average epidemic size $s_{\rm av}$.
As expected, the increase in the number of mask wearers decreases the average epidemic size and also increases the epidemic threshold.
Furthermore, a high mask coverage results in lowered $P_{\rm L}$ and $s_{\rm L}$ (not shown), i.e., the higher the number of mask wearers, the lower the occurrence probability of a large outbreak, and even if it does occur, the mask wearers will reduce its propagation.

 \begin{figure}[tb]
 \centering
(a)
  \includegraphics[width=6.5cm]{./figures/fig3a.eps}
(b)
  \includegraphics[width=6.5cm]{./figures/fig3b.eps}
 \caption{
(a) Occurrence probability $P_{\rm L}$ of the mask model in the Poisson network with different assortative mixing patterns: assortative ($r=0.95$; red), neutral ($r=0$; black), and disassortative cases ($r=-0.95$; blue).
The symbols represent the simulation results of $n = 10^5$ nodes, and the lines represent the analytical results based on Eq.~(\ref{eq:pG}).
The epidemic threshold above which $P_{\rm L}$ becomes nonzero is the lowest (highest) in the assortative (disassortative) case.
On the other hand, $P_{\rm L}$ for high-transmissibility cases ($T \approx 1$) is smaller in the assortative case than in the disassortative case (the situation is the same for $s_{\rm L}$ and $s_{\rm av}$).
Assortative mixing reduces the probability and degree of a large outbreak in high-transmissibility cases, although it allows a large outbreak to occur even in low-transmissibility cases.
(b) Epidemic threshold $T_c$ as a function of $r$ when $p_1=0.1,\, 0.3,\, 0.5,\, 0.7$, and $0.9$.
Here $T_c$ is evaluated from Eq.~(\ref{eq:thresholdCondition}).
It is obvious that $T_c$ continuously decreases as $r$ increases.
 }
 \label{fig:severalMaskPatterns}
 \end{figure}

Next, we discuss the impact of the assortative mixing by mask-wearing on epidemic spreadings, fixing the mask coverage at $p_1=0.5$ (where the assortativity coefficient $r$ is in the range of $[-1,1]$).
Figure~\ref{fig:severalMaskPatterns} (a) shows a comparison of the occurrence probabilities for the assortative ($r=0.95$), neutral ($r=0$), and disassortative cases ($r=-0.95$).
As indicated by both the theoretical line and simulation data, the epidemic threshold is the lowest (highest) when the mixing based on mask-wearing is assortative (disassortative).
Disassortative mixing is effective in suppressing epidemics in that a high transmissibility is required for large outbreaks.
On the other hand, $P_{\rm L}$ for high-transmissibility cases ($T \approx 1$) is smaller in the assortative case than in the disassortative case (the situation is the same for $s_{\rm L}$ and $s_{\rm av}$).
In the assortative case where non-mask wearers tend to be connected, large outbreaks are more likely to occur in the groups of non-mask wearers.
Simultaneously, in the assortative case, mask wearers are often surrounded by mask wearers, and there is a good probability of containing the epidemic in the early stages, if the seed node is a mask wearer, even with high transmissibility.
The separation of mask wearers and non-mask wearers reduces the probability and degree of a large outbreak in high-transmissibility cases, although it allows a large outbreak to occur even in low-transmissibility cases.
Note that the value of $P_{\rm L}$ (also $s_{\rm L}$ and $s_{\rm av}$) in a high $T$ case is affected by the mask effectiveness, in addition to assortative mixing by mask-wearing, and mask coverage.

 \begin{figure}[tb]
 \centering
  (a)
  \includegraphics[height=6.5cm]{./figures/fig4a.eps}
  (b)
  \includegraphics[height=6.5cm]{./figures/fig4b.eps}
 \caption{
Impact of assortative mixing on (a) epidemic threshold $T_c$ and (b) average epidemic size $s_{\rm av}$ with $T=1$ in the $(p_1, r)$ plane.
Here, the color-coded values of $T_c$ and $s_{\rm av}$ are evaluated from Eq.~(\ref{eq:thresholdCondition}) and Eq.~(\ref{eq:sAverage}), respectively.
The black curved lines represent $r_{\rm min}(p_1)$ given as Eq.~(\ref{eq:rmin}).
In panel (a), $T_c$ has a minimum (maximum) value at $r=1$ ($r=r_{\rm min}$) for a given mask coverage $p_1$ ($0<p_1<1$).
Panel (a) illustrates that $T_c$ is generally increased by decreasing $r$ or by increasing $p_1$.
The red circles (the green squares) in panel (b) indicate the position of $r$ at which $s_{\rm av}$ reaches the minimum (maximum) for a given $p_1$.
On the right side of the white dashed line, a large outbreak does not occur even when $T=1$ (i.e., $T_c=1$ and $s_{\rm av}=0$).
Irrespective of $p_1$, $s_{\rm av}$ reaches the minimum in the strongly assortative case ($r \approx 1$).
On the other hand, $s_{\rm av}$ reaches the maximum in a moderately assortative case when the mask coverage is high ($p_1 \gtrsim 0.55$).
}
 \label{fig:colorMap}
 \end{figure}

Assortative mixing by mask-wearing alters the epidemic threshold.
Figure~\ref{fig:severalMaskPatterns} (b) plots the epidemic threshold $T_c$ as a function of the assortativity coefficient $r$ with different mask coverages: $p_1 = 0.1,\, 0.3,\, 0.5,\, 0.7$, and $0.9$.
It is obvious that the epidemic threshold continuously decreases with increasing $r$.
As $r$ increases, the separation of mask wearers and non-mask wearers becomes more pronounced, resulting in large outbreaks among non-mask wearers more easily.
At $r = 1$, regardless of the value of $p_1(<1)$, non-mask wearers are completely detached from mask wearers and form their own giant component~\footnote{In assortatively mixed random networks, even if the fraction of non-mask wearers is small, as long as $p_1<1$, the number of non-mask wearers is infinite and the degree distribution of non-mask wearers' component follows $p(k)$, thus a giant component can be formed.}, so the epidemic threshold becomes determined by whether the infection spreads among non-mask wearers, which is consistent with the known one (\ref{eq:knownThreshold}) (see also Eq.~(\ref{eq:thresholdAssortative})).
In addition, we find that the impact of assortative mixing on the epidemic threshold increases as the mask coverage increases.
Figure~\ref{fig:colorMap} (a) illustrates the color-coded values of $T_c$ in the $(p_1,r)$ plane.
The figure shows that the epidemic threshold is generally increased by disassortative mixing by mask-wearing or by increasing the mask coverage.
For each $p_1$ value ($0<p_1<1$), $T_c$ has a minimum (maximum) value at $r=1$ ($r=r_{\rm min}$).

Although the assortative mixing lowers the epidemic threshold, it suppresses $P_{\rm L}$, $s_{\rm L}$, and $s_{\rm av}$ when the transmissibility is high, as already mentioned in Fig.~\ref{fig:severalMaskPatterns} (a).
Figure~\ref{fig:colorMap} (b) shows the average epidemic size $s_{\rm av}$ at $T=1$.
Irrespective of $p_1$, $s_{\rm av}$ reaches the minimum in the strongly assortative case ($r \approx 1$).
That is, the mask use is most effective in terms of reducing $P_{\rm L}$ when $r \approx 1$.
 This, however, does not imply that the average epidemic size for high $T$ always decreases more when the mixing is assortative than when it is neutral or disassortative.
When the mask coverage is high ($p_1 \gtrsim 0.55$), $s_{\rm av}$ reaches the maximum in a moderately assortative case (see the green line in Fig.~\ref{fig:colorMap} (b)).
The same holds for $P_{\rm L}$ and $s_{\rm L}$ (not shown).

Finally, we consider the example of heterogeneous networks.
We employ a scale-free random network having $p(k) = k^{-\gamma} / \sum_{k'=k_{\rm min}}^{k_{\rm max}} k'^{-\gamma}$ ($k=k_{\rm min}, k_{\rm min}+1, \ldots, k_{\rm max}$).
Here, we set $k_{\rm min}=2$, $k_{\rm max}=500$, and $\gamma=2.625$ so that $\langle k \rangle \approx 4$.
Although the epidemic threshold for a scale-free network with finite $k_{\rm max}$ is comparatively low regardless of the assortative mixing pattern, it is the lowest in the assortative case (Figs.~\ref{fig:heterogeneousNetwork} (a) and (b)).
(Note that this small but nonzero epidemic threshold is due to the finiteness of $k_{\rm max}$. If $k_{\rm max}$ is unbounded, $T_c$ becomes zero, regardless of the assortative mixing pattern, since $G_q'(1)$ diverges.)
Figures~\ref{fig:heterogeneousNetwork} (a) and (b) reveal that the occurrence probability and the average epidemic size in high $T$ cases are almost identical in the assortative and disassortative cases.
We see from Fig.~\ref{fig:heterogeneousNetwork} (c) that the average epidemic size at $T=1$ is minimized in the strongly assortative case ($r \approx 1$), but it also exhibits an equally small value in the maximally disassortative case ($r \approx r_{\rm min}$), provided that the mask coverage is not low.
Conversely, the mask use is less effective when the network is moderately assortative (as evidenced by the large $s_{\rm av}$), as in the case of Poisson networks.
We conclude that the effectiveness of mask use is enhanced when the mixing pattern is strongly assortative so that the network is split into two groups of mask wearers and non-mask wearers, or when the mixing pattern is strongly disassortative, which reduces the number of high-risk edges between non-mask wearers significantly.

 \begin{figure}[tb]
 \centering
 \vspace{3cm}
(a)
 \includegraphics[width=4.75cm]{./figures/fig5a.eps}
(b)
 \includegraphics[width=4.75cm]{./figures/fig5b.eps}
(c)
\includegraphics[width=4.0cm]{./figures/fig5c.eps}
 \caption{
Results of the mask model in a scale-free network:
(a) occurrence probability of large outbreaks $P_{\rm L}$ and (b) average epidemic size $s_{\rm av}$ as a function of $T$ in the assortative ($r=0.95$; red), neutral ($r=0$; black), and disassortative cases ($r=-0.95$; blue).
The symbols represent the simulation results of $n = 10^5$ nodes, and the lines represent the analytical results based on Eq.~(\ref{eq:pG}) for panel (a) and Eq.~(\ref{eq:sAverage}) for panel (b).
The degree distribution is $p(k) = k^{-2.625} / \sum_{k'=k_{\rm min}}^{k_{\rm max}} k'^{-2.625}$ for $k_{\rm min} \le k \le k_{\rm max}$, where $k_{\rm min}=2$ and $k_{\rm max}=500$.
The epidemic threshold is the lowest (highest) in the assortative (disassortative) case, although it is comparatively low regardless of the mixing pattern.
The occurrence probability and the average epidemic size in high $T$ cases are almost identical in the assortative and disassortative cases.
(c) Color-coded values of $s_{\rm av}$ at $T=1$ in the $(p_1,\;r)$ plane.
The black curved lines in panel (c) represent $r_{\rm min}(p_1)$, which is given as Eq.~(\ref{eq:rmin}).
The red circles (the green squares) in panel (c) indicate the position of $r$ at which $s_{\rm av}$ reaches the minimum (maximum) for a given $p_1$.
The average epidemic size at $T=1$ is minimized in the strongly assortative case ($r \approx 1$), but it also exhibits an equally small value in the maximally disassortative case ($r \approx r_{\rm min}$), provided that the mask coverage is not low.
 }
 \label{fig:heterogeneousNetwork}
 \end{figure}

\section{Summary \label{sec:summary}}

In this work, we studied the mask model, which is a SIR-type epidemic model involving mask wearers and non-mask wearers, in assortatively mixed random networks to discuss the impact of assortative mixing by mask-wearing on the effectiveness of mask use.
We derived the occurrence probability and mean size of large outbreaks, epidemic threshold, and average epidemic size using the probability generating functions.
In both the Poisson and scale-free random networks, we established that the assortative mixing by mask-wearing decreases the epidemic threshold.
It means that large outbreaks can occur more easily as a network is assortative.
In high-transmissibility cases, on the other hand, the average epidemic size is minimized when the mixing pattern is strongly assortative.
Strongly assortative mixing, resulting in the separation of mask wearers and non-mask wearers, reduces the probability and degree of a large outbreak in high-transmissibility cases, although it allows a large outbreak to occur even in low-transmissibility cases.
In scale-free networks, the mask use in the maximally disassortative case is equally effective in terms of reducing the probability and mean size of large outbreaks as in the strongly assortative case (when the mask coverage is not low).
We also performed MC simulations of the mask model in the Poisson and scale-free networks, which validated our analytical treatment results.

This study has the following limitations:
(i) the epidemic model employed is a discrete-time SIR model with constant transmissibility, which is very simple.
The effectiveness of the mask use should be discussed further in more general epidemic models.
(ii) We assumed random networks.
Real networks have structural properties, such as degree correlation, high clustering, and community structure.
A recent study~\cite{duh2019assortativity} reported that degree correlation in a network affects human behaviors, which may be related to epidemic control.
A further study of how the combination of network structure and mixing pattern by mask-wearing affects the spread of the disease should be conducted.
Furthermore, we may consider multiplex or multilayer networks~\cite{boccaletti2014structure, bianconi2018multilayer} or simplicial complexes~\cite{torres2020simplicial, bianconi2021higher}.
Using recent bond percolation literatures~\cite{kryven2017finite,kryven2019bond,bianconi2019percolation,balogh2020networks}, we could investigate mask models on such networks.
(iii) The present study discussed the impact of assortative mixing by mask-wearing.
One can consider the assortative mixing based on other attributes, such as age, gender, obesity, and vaccination~\footnote{After the submission of this manuscript, the impacts of vaccination homophily (assortative mixing by vaccination adoptions) on disease dynamics have been reported~\cite{hiraoka2021herd,burgio2022homophily}.}.
Regarding the mask use, one may further consider the mask quality (cotton mask, surgical mask, or N-95 mask) as an attribute rather than usage~\footnote{In preparation of this manuscript, Tian et al.~\cite{tian2021role} proposed a multi-type mask model, where each node wears a different type of mask.}.
The analysis provided in this study will be easily generalized to discuss some of the limitations mentioned above.
It is hoped that the present study will contribute to a better understanding of epidemic propagation patterns in our society.

\section*{Acknowledgement}
This work was supported by JSPS KAKENHI Grant Numbers JP19K03648 and JP21H03425.


\section{Introduction}

The use of standards is a fundamental step for reusability of digital resources. In Terminology, the reference ISO standard ISO 30042:2019,\footnote{ISO 30042:2019 Management of terminology resources - TermBase eXchange (TBX), \url{https://www.iso.org/standard/62510.html}} describes an XML-based terminology exchange format, designed to make terminology databases easier and safer to maintain, distribute, and use. Nevertheless, this standard is not immediately useful when a researcher want to expose the data structured as Linked Data.
Attempts to fill this gap has been proposed in the literature \cite{dibuono_et_al_2020a,dibuono_et_al_2020b}; however, there are still some manual corrections or interventions that are needed in order to make the conversion perfect~\cite{reineke_2014}, as well as problems in the depth of the information that can be represented with other models more concept-oriented~\cite{reineke_romary_2019}.

In this work, we want to discuss the problem from a different perspective. We make a step back and look at these issues as a data modeling problem. Our hypothesis is that if the terminological data are modeled correctly at an abstract level, then the data themselves can be exported or exposed in any format and maintain both the linguistic dimension and its relationship to the conceptual dimension. For this purpose, we will use the Entity-Relationship (ER) diagram as our tool for the conceptual model of the problem and the explanation of the associations between the object of the study.

The paper is organized as follows: in Section~\ref{sec:conceptual} we introduce the basic elements of the conceptual data modeling and the pictorial representation that we use, the Entity-Relationship (ER) diagram. In Section~\ref{sec:model}, we analyze step by step four approaches in Terminology, onomasiological, semasiological, ontoterminological, frame-based, and we design the corresponding conceptual data model. In Section~\ref{sec:final}, we give our final comments and a view on the current and future work.

\section{Conceptual Data Modeling: ER Diagrams}\label{sec:conceptual}

A conceptual data model represents the objects of the domain of interest and the relationships among them. It is a (not-so) complicated visual language which is effective for communicating and describing how the knowledge is organized at an abstract level without any specific reference to how the data will be physically implemented. The Entity-Relationship (ER) Model~\cite{chen_1976} is a model of this kind, and its pictorial representation is called ER diagram. This diagram uses three basic elements, shown in Figure~\ref{fig:er_schema} that describes all the important semantic information about the domain of interest: entities, associations, attributes. An entity is the class of elements that are described with a set of attributes, some of the attributes are `identifiers' (distinguished with a black dots from the non-dientifiers) and distinguish uniquely each element of the class.  An association describes a class of elements that represent the relationships between entities.

\begin{figure}[t]
\begin{center}

\includegraphics[width=\textwidth]{figure/ER_constructs.pdf}

\caption{ER diagram basic constructs: entity, association, attributes}
\label{fig:er_schema}
\end{center}
\end{figure}

To make an example, let us suppose to model the following problem described by Roche~\cite{roche_2012}:
\begin{quote}
[...] unlike classical terminology ontoterminology preserves the diversity of language between different communities of practice since they share the same domain and standardised conceptualisation. In point of fact, two different terms can denote the same concept whose identifier should be built so that we understand the right place of the concept in the ontology.
\end{quote}

The previous paragraph describes a particular requirement of the domain of interest where a concept can be designated by more than one term. The corresponding ER diagram is shown in Figure~\ref{fig:er_concept}. There are two entities: the entity \textsc{concept} which represents the class of all the elements (the concepts) that are uniquely specified by an identifier, and the entity \textsc{term} which represents the class of elements (the terms) that have a unique identifier and a designation, i.e. the sequence of characters that compose the (multi-)word that designates the term.
The association Denoted describes the relationship between pairs of elements (concept, term). We have an additional information in this diagram, a pair of numbers that indicates a constraint on the participation of an entity to an association. This information is called cardinality and it tells the minimum and the maximum number of times that an element of an entity can participate to the relationship.
In this case, on the segment that connects \textsc{concept} to Denoted we have a cardinality $(1, n)$ which means that an element of \textsc{concept} must be denoted by at least $1$ term and at most by $n$ terms (we use $n$ when there are no specified numbers). On the other side of the association, each element of \textsc{term} must be associate to 1 and only 1 element of \textsc{concept}.

\begin{figure}[t]
\begin{center}

\includegraphics[width=\textwidth]{figure/ER_concept_term.pdf}

\caption{ER diagram for concepts denoted by terms}
\label{fig:er_concept}
\end{center}
\end{figure}

This simple example already contains some interesting food for thoughts: while the cardinality of the \textsc{concept} derives directly from the textual description of the problem, the cardinality $(1, 1)$ of \textsc{term} was our implicit choice to say that \textit{homographs} must be uniquely identified even though the designation is the same. For example, the sequence of characters \textit{bank} is the value of the attribute designation of two different (and uniquely identified) elements of \textsc{term} that refer to the concept that describes the side of a river and the other concept that refers to the financial institution.

\begin{figure}[t]
\begin{center}

\includegraphics[width=\textwidth]{figure/ER_concept_term_extensional.pdf}

\caption{Extensional representation of concepts denoted by terms}
\label{fig:er_concept_extensional}
\end{center}
\end{figure}

To conclude this section, we show in Figure~\ref{fig:er_concept_extensional} the extensional representation (sets of elements are explicitly depicted) of the ER diagram. In the background and greyed out, we have the same ER diagram of Figure~\ref{fig:er_concept}, in the foreground, we have the sets (classes) of elements and the explicit links to the association. Concept \textsc{c1} is denoted by two terms, \textsc{t1} and \textsc{t2}, while Concept \textsc{c2} is denoted by term \textsc{t3} only.

\section{From The General Theory of Terminology to Terminological Data Modeling Approach}\label{sec:model}

Despite the title of this section, it is not our intent to enter into the debates and the criticisms to the General Theory of Terminology of Eugene W\"uster.
We will rather propose an analysis of the requirements of the different theories and approaches to terminology summarised in~\cite{cabre_2003}, that starts from the theory proposed by W\"uster and the design of its corresponding ER diagram.

\subsection{Onomasiological approach}

The main aspect of the General Theory of Terminology proposed by W\"uster  revolves around the notion of concepts, while terms are a sort of ``consequence'' of the organization of the concepts. In the onomasiological approach, concepts are delineated by delimiting characteristics~\cite{suonuuti_2012} while terms are assigned permanently to concepts; in particular, term denotes only one concept. Concepts can be placed in a concept system that describes how concepts are related to each other.

If we analyze these requirements and compare them with the previous example shown in Figure~\ref{fig:er_concept}, there are at least two new classes of objects involved in the problem: the characteristics of a concept, and the relations among concepts to describe the conceptual system. If we carry out a short analysis of requirements, we find that, according to the ISO 1087:2019~\cite{iso_1087_2019}, there are `type of characteristics' that are categories of characteristics that are grouped (for the purpose of the terminological analysis), essential, non-essential, and delimiting characteristics. Moreover, there are `hierarchical' or `generic' relations among concepts.

\begin{figure}[t]
\begin{center}

\includegraphics[width=\textwidth]{figure/ER_onomasiological.pdf}

\caption{Simplified ER diagram for the onomasiological, semasiological, ontoterminological, frame-based approaches.}
\label{fig:er_onomsaiological}
\end{center}
\end{figure}

In Figure~\ref{fig:er_onomsaiological}, we show the resulting ER diagram that includes both the characteristics of a concept and the hierarchical relations among them. For space reasons, the generic relations are not displayed (it would be a similar association that connects \textsc{concept} with itself). The entity \textsc{characteristic} represents the class of elements identified by the `name' (we assume that the name of the characteristic is unique) and has an attribute `variety' to describe whether the element belongs to the type, essential, non-essential, or delimiting characteristic. The association Group describes the fact that: 1) on the one hand, an element of \textsc{characteristic} can be part (see here that it is not mandatory for a characteristic to be part of a group) of one or more types of characteristics (the segment \textit{is part of}), 2) on the other hand, the segment \textit{Has elements} indicate an element `type of characteristic' that can group one or more characteristics. The association Hierarchical models the hierarchy of sub-ordinate and super-ordinate concepts. The cardinality show that an element of \textsc{concept} can (it is optional) have sub-ordinate concepts, while every concept must have at most one superordinate concept (it is not the goal of this paper to discuss whether a concept can belong to different hierarchies or the situation where the concept-of-concepts has as a super-ordinate concept itself).

\subsection{Semasiological approach}

As Santos and Costa argue in~\cite{santos_costa_2015}:
\begin{quote}
That the concept is a central element is not, for us, the question. The question is, knowing if in terminological work it is always the starting point as W\"uster argues or if, depending on the circumstances, it can be the point of arrival.
\end{quote}
The semasiological approach stems from a linguistic analysis of the texts to study the language `in action' and identify terms before concepts. Despite being clearly a sort of `opposite direction' flow of analysis (it is not our intent to oppose the two approaches), the relative conceptual data model is pretty much identical to the one shown in Figure~\ref{fig:er_onomsaiological}. We may want to add additional entities and relationships to document, for example, the \textsc{text} and the \textsc{collection} from which the term was extracted, as well as new types of relationships such as `connected to', `consists of', 'is a', to build the concept map of a domain of knowledge (we would strongly suggest to do it in a real case scenario,). Nevertheless, we witness the fact that the core structure of the conceptual data model is the same to that of the onomasiological one.

\subsection{Ontoterminological approach}

Continuing with the idea of this dichotomy of flow of analysis about how to build knowledge, we analyze the paradigm called ontoterminology~\cite{roche_2012}: ``a terminology whose conceptual system is a formal ontology relying on epistemological principles''. The aim of this paradigm derives from the need to make the operationalization on terminology easier and formally correct, especially from a description logic point of view.

One of the most important features of this paradigm is the double semantic triangle that maps the linguistic knowledge extracted from text to the ontologies derived from domain experts. As a consequence, the definition for a term and for a concept are different in the sense that the former derives from the natural language explanation, while the latter is a formal definition in terms of the characteristics.

Once again, whereas this approach is different from the other two,
the conceptual data model is still the same. There is one important aspect that we must take into consideration, though: the definition of a term should be included (probably as an attribute of \textsc{term}) to the ER diagram in Figure~\ref{fig:er_onomsaiological}.

\subsection{Frame-based approach}

A similar perspective on the problem of knowledge driven versus corpus driven terminology and on how to enhance the linguistic properties of an ontology is provided by frame semantics~\cite{fillmore_1985}. In the frame based terminology approach~\cite{faber_2015}, the notion of semantic frames which is a structural background knowledge able to construct the meaning of terms is applied. A semantic frame ``models a given situation; situations comprise participants, props, and other conceptual elements, which constitute its frame elements''\cite{lhomme_2018}.

Even in this approach, the core structure of the conceptual data model remains the same. Frames and frame elements could be represented as additional entities related to terms, and their relationships, that describe framework for the definition and representation of specialized knowledge units.

\section{Final Considerations and Future Work}\label{sec:final}

In the previous section, we have presented an investigation of how to model the different theories in terminology. We do not intend to create a new integrated theory, we rather want to study whether a unified conceptual data model that accepts all the different facets that experts in terminology have proposed exists. In this regard, our proposal is  similar to the Cabr\'e's``theory of doors''\cite{cabre_2003} where the terminological unit is the object under observation that can be accessed from different perspectives (the ``doors'') in such a way that the central object, the terminological unit, is directly addressed whether starting from the concept or the term (or even `the situations' as described in her work). In our design of the conceptual data model, the core data structure is the object that can be accessed and the ``approaches'' are only a different view of the same data. Of course, this does not mean that different approaches will produce different data for the same object under study, but the idea is that different views can enhance the information about the term (ontological, linguistic, practical information) or simply access portions of data.

Moreover, the ER model can be used as a basis for unification of different views of data: the network model, the relational model, and the entity set model. This fact will be very useful for what comes next. The positive consequence of this homogeneous conceptual data model is that we can decide what `direction' we want (or need) for the implementation of the structure that will represent the data and how we want to export them.
For example, we can:
\begin{itemize}
\item transform the ER diagram into the corresponding relational database and export the standard TermBase Exchange (TBX) file~\cite{vezzani_dinunzio_2020};
\item provide an experimental infrastructure for the study of the properties of terms and their relations with distributional semantics\cite{bonato_et_al2021}, or the study of the quality of machine translations~\cite{cambedda_et_al_2021};
\item make a comparative qualitative quantitative analysis of the use of specialized languages in narrative~\cite{vezzani_dinunzio_2019};
\item study the performance of query variation for technology assisted medical systematic review systems~\cite{dinunzio_vezzani_2021}.
\end{itemize}

In addition, we may want to increase the level of interoperability by exposing the data with a Linked Open Data approach. Since the ER diagram can be mapped to a graph model, we can directly map the overall conceptual data model to a graph and expose it, for example in RDF \cite{dibuccio_et_al_2013,dibuccio_et_al_2014}.

\bibliographystyle{splncs04}