July 08, 2026
Decentralized federated learning (DFL) dispenses with the central server of classical FL by utilizing peer-to-peer model exchanges among edge devices. This server-free architecture enables ad-hoc, flexible distributed learning in large device-to-device (D2D) networks. However, wireless DFL converges slowly because peer-to-peer model aggregation incurs high delays and errors. Each DFL training round involves many-to-many gradient sharing over wireless channels, resulting in uncoordinated channel access, large communication errors from stragglers, and slow model consensus, especially in large-scale D2D networks with pronounced clustering structures. We address these aggregation bottlenecks by provisioning a few reliable backhaul links at straggling nodes to enhance network connectivity. Building on this idea, our budget-aware, cluster-centric DFL framework first partitions the network into densely connected clusters, and then allocates the limited backhaul budget to selected cluster heads. The resulting two-tier protocol executes fast, parallel model aggregation within clusters and infrequent inter-cluster exchanges among the heads, yielding an \(\mathcal{O}(1/t)\) convergence rate in \(t\) iterations. Numerical experiments on image-classification tasks confirm that our approach accelerates convergence compared to state-of-the-art DFL baselines with only a few strategically placed backhaul links.
Decentralized federated learning, device-to-device network, gossip algorithm, over-the-air computation, device clustering.
Federated learning (FL) has emerged as a compelling framework for distributed training of artificial-intelligence (AI) models at the network edge, where both data and computation are distributed over edge devices [1]. Classical FL relies on a central server (e.g., a base station or edge server) to coordinate model aggregation: edge devices train local models and periodically upload their models or gradients to the server for global consensus. Although effective, this server-centric architecture becomes impractical whenever device-to-server connections are unavailable, unreliable, or undesirable, particularly in ad hoc networks [2]. Decentralized FL (DFL) overcomes this limitation by replacing the server with peer-to-peer gossip exchanges [3], thereby aligning with the connectivity patterns and privacy requirements of large-scale device-to-device (D2D) networks and the Internet-of-Things (IoT) [4].
A growing body of work confirms that model communication and aggregation, rather than local model computation, dominate the runtime of FL systems [5]. Exchanging high-dimensional model parameters across rate-limited links is extremely costly for large networks. The challenge is amplified in DFL as every device must alternately transmit to and receive from its neighbors, turning the “many-to-one" uplink of classical FL into “many-to-many" communication [6]. The absence of a coordinating server further complicates decentralized multiple-access control, leading to collisions, long delays, and ultimately slower convergence.
To mitigate the model aggregation bottlenecks, recent research advocates over-the-air (OTA) computation as a scalable solution. By exploiting the signal-superposition property of the wireless multiple-access channel, OTA computation sums gradients “in the air" by transmit scaling, allowing all devices to transmit simultaneously to a single receiver over the same channel [7]. When extended to DFL [8]–[10], OTA aggregation reduces the peer-to-peer gossip delay among \(N\) devices from \(\mathcal{O}(N^{2})\) transmissions to \(\mathcal{O}(N)\).
Despite its promise, OTA-based DFL still faces three key hurdles: (a) without a server, devices struggle to schedule concurrent transmissions and receptions [8]; (b) OTA aggregation relies on transmitter-side channel inversion to align the received gradients, forcing all transmitters to scale to the stragglers with weakest links, which inflates noise and slows convergence [11]; and (c) large D2D networks typically exhibit “small-world’’ structure, where dense clusters bridged by only a few inter-cluster edges. As a consequence, the sparsely connected clusters throttle global consensus [6].
These observations suggest harnessing the natural clustering of “small-world" networks to localize OTA aggregation within dense clusters and reinforce the few critical long-range edges. This insight motivates supplementing DFL with a few extra reliable backhaul links to assist the sparsely connected stragglers, bridge otherwise isolated clusters, and, in turn, accelerate model aggregation. Specifically, we propose a budget-aware, cluster-centric DFL framework that (a) partitions the network into densely connected clusters, (b) identifies straggling devices whose poor connections dominate aggregation error, and (c) deploys a small set of extra reliable backhaul links between selected cluster heads. The resulting two-tier protocol performs fast, head-centric OTA aggregation inside each cluster and low-frequency inter-cluster exchanges among the heads, achieving an \(\mathcal{O}(1/t)\) convergence rate in the optimality-gap bound after \(t\) iterations. Finally, we develop a low-complexity clustering algorithm that minimizes the convergence rate bound under link-deployment constraints using only long-term channel statistics. Notably, our approach generalizes the spirit of hierarchical FL to a fully decentralized setting, removing the need for dedicated cloud servers when partitioning users. Moreover, it naturally adapts to heterogeneous ad-hoc networks: clusters can leverage existing infrastructure, such as routers or factory controllers, as pre-established inter-cluster links, further trimming costs.
In this work, we consider a DFL system operating over a D2D network, where \(N\) single-antenna IoT devices collaboratively train a shared AI model using their local datasets. We assume that these devices form a connected network through wireless connections where there is a routing path connecting any two devices. The global objective is to minimize \[\min_{\xv \in \mathbb{R}^d} f(\xv)=\min_{\xv \in \mathbb{R}^d} \sum_{n=1}^N \frac{|\mathcal{D}_n|}{D} f_n(\xv), \label{eq:32grand32objective}\tag{1}\] where \(\xv \in \mathbb{R}^d\) represents the model parameters of dimension \(d\), \(\mathcal{D}_n\) is the dataset on device \(n\), \(D=\sum_n |\mathcal{D}_n|\) is the total number of data, and \(f_n(\cdot)\) is the local loss of device \(n\). Specifically, the local loss is evaluated over all samples in the local dataset \(\mathcal{D}_n\) as \[f_n(\xv)=\frac{1}{|\mathcal{D}_n|} \sum_{i=1}^{|\mathcal{D}_n|} f_n(\xv;\xi_{n,i}), \label{eq:32local32objective}\tag{2}\] where \(\xi_{n,i}\) denotes the \(i\)-th sample in \(\mathcal{D}_n\).
DFL proceeds iteratively with on-device local training and inter-device model aggregation. Due to the network connectivity and data privacy concerns, we consider model communication occurs only among neighboring devices over wireless fading channels, with no central server involved. This fully decentralized architecture aligns naturally with dense D2D IoT systems, where devices primarily share local connections rather than relying on a base station.
Conventional DFL algorithms, e.g., [6], [10], typically employ gossip algorithms to disseminate and aggregate local information to enforce global model consensus. For each training round \(t=1,2,\ldots,T\), a standard gossip-based DFL procedure comprises the following steps:
Local gradient computation: Given its current model \(\xv_n^{(t)}\), device \(n\) draws a mini-batch \(\mathcal{B}_{n}^{(t)}\subset\mathcal{D}_n\) and runs mini-batch stochastic gradient descent (SGD) to obtain the stochastic gradient vector \({\nabla} f_{n}(\xv_n^{(t)};\mathcal{B}_{n}^{(t)})\).
Gradient exchange: Let \(\mathcal{N}_n\) denote the neighbors of device \(n\) that share direct wireless links with it. Device \(n\) transmits its model update with respect to \(\hat{\xv}_n^{(t)}\) to all \(m\in\mathcal{N}_n\) and simultaneously receives their updates. Due to communication impairments (e.g., wireless fading, compression loss), the gradient received from \(m\) arrives perturbed. We denote the received model update from \(m\) to \(n\) by \(\gv_{m\to n}^{(t)}\).
Model aggregation: After collecting \(\{\gv_{m\to n}^{(t)}\}_{m\in\mathcal{N}_n}\) from its neighbors, each device \(n\) updates its model via (cf. [6]) \[\begin{align} &\xv_n^{(t+1)}= \xv_n^{(t+\frac{1}{2})} + \gamma \sum_{m\in\mathcal{N}_n \cup \{n\}}w_{nm} (\widehat\xv_m^{(t+1)}- \widehat\xv_n^{(t+1)}),\label{Eq3} \end{align}\tag{3}\] where \(w_{nn}\geq 0\) and \(w_{nm}\geq 0\) are the mixing weights controlling the balance between local and neighboring model information and satisfying \(w_{nn}+\sum_{m\in \mathcal{N}_n} w_{nm}=1\), \(\gamma\) denotes the consensus stepsize, \(\xv_n^{(t+\frac{1}{2})}=\xv_n^{(t)}-\eta_t {\nabla} f_{n}(\xv_n^{(t)};\mathcal{B}_{n}^{(t)})\) is the locally updated model with stepsize \(\eta_t\), and \(\widehat \xv_m^{(t)}=\xv_m^{(0)}-\sum_{t^\prime=1}^{t-1} \gv_{m\to n}^{(t^\prime)}\) is the running estimate of \(m\)’s model reconstructed at device \(n\) from the received model updates.
It has been extensively documented that exchanging gradient information constitutes the critical bottleneck in DFL systems. Iterative transmission and reception of local gradients over every pair of neighboring devices impose a prohibitively large latency and introduce communication errors that slow down model consensus and hinder convergence. Our work aims to address this bottleneck through an efficient device scheduling and gradient aggregation algorithm tailored to fully decentralized networks.
OTA computation is a scalable “many-to-one" gradient uploading scheme for dense networks. Since our proposed method employs OTA schemes to facilitate inter-device aggregation, we summarize the underlying channel model and OTA aggregation methods as follows.
For every existing wireless link between devices \(m\) and \(n\), the instantaneous channel coefficient at iteration \(t\) is assumed to follow a Rayleigh fading model, as \[\begin{align} h_{mn}^{(t)}=\sqrt{\alpha_{mn}}\,\beta_{mn}^{(t)}, \qquad n\in[N],\; m\in\mathcal{N}_n,\label{eq95channel} \end{align}\tag{4}\] where \(\alpha_{mn}\) is the large-scale fading (assumed constant during training) and \(\beta_{mn}^{(t)}\sim \mathcal{CN}(0,1)\) is the small-scale fading. Because large-scale coefficients vary slowly, \(\{\alpha_{mn}\}\) are treated as known a priori to the DFL system design. In contrast, the fast-varying coefficients \(\{\beta_{mn}^{(t)}\}\) are estimated at the start of each training iteration. Therefore, instantaneous channel state information (CSI) is available for the OTA transceiver design but not for DFL device scheduling.
Consider a fixed iteration \(t\), where the index \(t\) is omitted when clear. Device \(n\) aims to receive and aggregate the stochastic gradients \(\{{\nabla} f_{m}(\xv_m;\mathcal{B}_{m})\}_{m\in\mathcal{N}_n}\) from its neighbors. To achieve this, each sender normalizes its gradient by the empirical mean and standard deviation, scales the normalized vector with a transmit
scalar, and transmits the vectors entry-by-entry over the shared channel. The transmitted signal from device \(m\in\mathcal{N}_n\) to \(n\) is given by \[\begin{align} \zv_{m\to n}=q_{m\to n}\, \frac{{\nabla} f_{m}(\xv_m;\mathcal{B}_{m})-\texttt{mean}\cdot\boldsymbol{1}}{\texttt{std}},\label{eq5}
\end{align}\tag{5}\] where \(q_{m\to n}\in\mathbb{C}\) is a complex scaling factor and mean, std are the common empirical statistics of all the gradients in this iteration.1 The average per-entry power is constrained by \(\Ebb[\lVert \zv_{m\to n}\rVert^2/d]=|q_{m\to n}|^{2}\le P_{0}\).
Device \(n\) receives a noisy superposition \(\yv_{n}=\sum_{m\in\mathcal{N}_n}h_{mn}\zv_{m\to n}+\nv_{n}\), where \(\nv_{n}\) is additive white Gaussian noise. By applying a linear receive scalar \(\omega_{n}\), the device estimates \[\begin{align} \sum_{m\in\mathcal{N}_n}w_{nm} \gv^{(t)}_{m\to n} \approx \frac{\texttt{std}}{\omega_{n}}\yv_{n}+\texttt{mean}\cdot\boldsymbol{1},\label{eq95otareceiver} \end{align}\tag{6}\] which is then used in the local update step in 3 .
Following [7], [11], the transceiver scalars that minimize the mean-squared error (MSE) of the aggregated gradients under the maximum power \(P_0\) are given by \[\omega_{n}=\sqrt{P_{0}}\min_{m\in\mathcal{N}_n}\frac{|h_{mn}|}{w_{nm}}, q_{m\to n}=\omega_{n}w_{nm}\frac{h_{mn}^{*}}{|h_{mn}|^{2}},\label{eq6}\tag{7}\] where \(h_{mn}^{*}\) denotes the complex conjugate of \(h_{mn}\).
Applying OTA aggregation to fully decentralized settings introduces several obstacles. In DFL, each device must both transmit to and receive from its neighbors, significantly complicating the gradient exchange procedure. Concretely, OTA-based DFL faces three key challenges, which hamper consensus and slow convergence:
Multiple-access control: Classic OTA aggregation assumes a single server that orchestrates channel access. In DFL, no such coordinator exists, making it difficult to schedule the channel access, so that every device can fully transmit to, and receive from, all of its neighbors without collisions.
Gradient aggregation error: According to 7 , OTA transmitters scale their gradients with transmit power inversely proportional to the channel gain to ensure that the superimposed signals arrive with the desired weights. If a link suffers from deep fading, all the other devices sharing that channel must reduce transmit power, which amplifies the overall aggregation error and decelerates convergence [10].
Network connectivity bottleneck: The consensus speed of the approach in 3 critically depends on the edge connectivity reflected in mixing weights \(\Wv=[w_{nm}]\). Real-world D2D topologies often follow a “small-world’’ pattern: dense clusters connected by only a handful of inter-cluster edges. Consequently, some poorly connected nodes at the cluster edge become bottlenecks that dominate the global consensus rate.
We develop a cluster-centric DFL model aggregation scheme to tackle the above limitations. The core idea is to provision a small number of extra reliable communication links between carefully selected devices (e.g., by backhaul connections), thereby strengthening the weakest parts of the network and accelerating global consensus. Because the placement of these links depends only on slow-varying large-scale channel statistics (e.g., large-scale fading coefficients), it can be optimized a priori and kept fixed throughout training.
As shown in Fig. 1, the original D2D graph is viewed as a collection of densely connected intra-cluster subgraphs bridged by sparsely connected inter-cluster edges. We exploit this “small-world" structure in two steps:
Device clustering: Using large-scale fading information, we partition the devices into clusters and pinpoint the few nodes whose poor connectivity in the intersections of clusters creates network bottlenecks.
Link enhancement and hybrid aggregation: A limited budget of high-quality links is allocated to these stragglers, yielding robust shortcuts for inter-cluster gossip. Model aggregation is then decomposed into (a) noisy OTA exchanges within each cluster and (b) reliable gossip-based exchanges across clusters.
This hybrid design simultaneously (a) alleviates the connectivity bottleneck by fortifying inter-cluster paths and (b) simplifies multiple-access scheduling, as gradient transmission now occurs alternately in intra-cluster and inter-cluster phases.
We represent the \(N\)-device D2D network by an undirected weighted graph \(\mathcal{G}\) whose edge weights correspond to the large-scale fading coefficients. The goal is to partition \(\mathcal{G}\) into \(K\) disjoint clusters, where both \(K\) and the cluster assignments are design variables. Denote by \(\mathcal{C}_k\) the node set of cluster \(k\) and by \(v_k\in\mathcal{C}_k\) its cluster head (centroid). We jointly compute \(K\), and \(\{\mathcal{C}_k,v_k\}_{k=1}^{K}\) subject to the following requirements:
Star-type intra-cluster topology: For every \(k\), the subgraph induced by \(\mathcal{C}_k\) must contain a star with \(v_k\) at its center; that is, each node in \(\mathcal{C}_k\) already shares a direct wireless link with \(v_k\). This ensures that intra-cluster aggregation can adopt a server-based FL protocol with \(v_k\) acting as the local “server."
Link-enhancement budget: We augment the connectivity among the heads \(\{v_k\}\) by provisioning a limited number of reliable links (e.g., wired backhaul). Let \(\mathcal{H}\) be the \(K\)-node weighted graph whose edges correspond to these links. Suppose that deploying a reliable transceiver at a node costs \(L_{N}\), and creating a link of length \(e(v_i,v_j)\) between \(v_i\) and \(v_j\) costs \(L_{D}\cdot e(v_i,v_j)\). The total expense must satisfy \[\begin{align} \label{eq8} K L_{N}+L_{D}\sum_{1\le i<j\le K}e(v_i,v_j)\le B, \end{align}\tag{8}\] where \(B\) is the overall budget and \(e(v_i,v_j)\) is the edge weight in \(\mathcal{H}\). Moreover, each reliable link must obey the physical reach limit \[\begin{align} \label{eq9} \max_{1\le i<j\le K}e(v_i,v_j)\le \tau, \end{align}\tag{9}\] where \(\tau\) is the maximum feasible distance.
Learning-rate maximization: Subject to the topology and budget constraints above, we select \(K\), \(\{\mathcal{C}_k,v_k\}\), and \(\mathcal{H}\) so that, under the hybrid aggregation protocol described next, the resulting DFL procedure achieves the largest possible learning rate.
Before formally formulating and solving the clustering problem, we first outline the OTA-based hybrid aggregation scheme over \(\mathcal{G}\) and \(\mathcal{H}\) in the sequel.
Given the clustering solution \(\{\mathcal{C}_k,v_k\}\) and its head–link graph \(\mathcal{H}\), the objective 1 is equivalent to \[\min_{\xv\in\mathbb{R}^d}\; \sum_{k=1}^{K} p_k\,f_{\mathcal{C}_k}(\xv), \label{eq:32grand32objective2}\tag{10}\] where \(p_k=D_{\mathcal{C}_k}/D\) with \(D_{\mathcal{C}}=\sum_{i\in\mathcal{C}}|\mathcal{D}_i|\) and \[f_{\mathcal{C}_k}(\xv)=\sum_{i\in\mathcal{C}_k} \frac{|\mathcal{D}_i|}{D_{\mathcal{C}_k}}\,f_i(\xv). \label{eq:32grand32objective3}\tag{11}\] During any training iteration \(t\) (index omitted for brevity), the model aggregation proceeds in two steps.
Intra-cluster aggregation. Every device \(i\in\mathcal{C}_k\) computes the stochastic gradient \({\nabla} f_{i}(\xv_k;\mathcal{B}_{i})\) and uploads it to the head \(v_k\) via the OTA scheme in 5 . Using 6 , the head recovers a noisy estimate of the weighted sum \(\sum_{i\in\mathcal{C}_k}\frac{|\mathcal{D}_i|}{D_{\mathcal{C}_k}}\, {\nabla} f_{i}(\xv_k;\mathcal{B}_{i})\), abbreviated as \(\widehat{\gv}_k\), with the mixing weight \(w_{iv_k}=\frac{|\mathcal{D}_i|}{D_{\mathcal{C}_k}}\). Then, it updates its intra-cluster model \(\xv_k\) by using the estimated aggregated gradient, and broadcasts it to all members of \(\mathcal{C}_k\).
Inter-cluster aggregation. For every \(H>0\) training iterations, the heads perform a noiseless gossip step over the reliable links in \(\mathcal{H}\) once: \[\begin{align} \xv_k \;\gets\; \sum_{l\in\mathcal{N}_{\mathcal{H}}(v_k)} p_l\,\xv_l + (1-\sum_{l\in \mathcal{N}_\mathcal{H}(v_k)} p_l) \xv_k, \label{eq12} \end{align}\tag{12}\] where \(\mathcal{N}_{\mathcal{H}}(v_k)\) denotes the neighbors of \(v_k\) in \(\mathcal{H}\). Restricting gossip to the \(K\) heads, rather than to all the \(N\) devices as in conventional DFL [6], [10], markedly accelerates global consensus. The inter-cluster communication interval \(H\) is introduced to reduce communication delay and boost the overall efficiency. The impact of \(H\) is investigated in Sections 4.1 and 5.
The complete procedure is summarized in Algorithm 2. The proposed approach speeds up convergence by combining server-based OTA aggregation within clusters and infrequent, noiseless gossip among a small set of well-connected heads, at the cost of a modest additional deployment budget for the reliable links.
This section first establishes the convergence rate of the proposed clustered DFL algorithm and then uses the results to optimize the device clustering variables of \(K,\mathcal{H}\), and \(\{\mathcal{C}_k,v_k\}\).
We begin with two standard assumptions, which are widely adopted in the distributed learning literature.
Assumption 1. Each local objective \(f_n\) in 2 is \(L\)-smooth and \(\mu\)-strongly convex, guaranteeing the existence and uniqueness of the global optimum \(\xv^*\) of 1 .
Assumption 2. For every device \(n\) and iteration \(t\), the mini-batch gradient \({\nabla} f_{n}(\xv_n^{(t)};\mathcal{B}_{n}^{(t)})\) is unbiased and has bounded variance and second moment: \[\begin{align} &\Ebb \| {\nabla} f_{n}(\xv_n^{(t)};\mathcal{B}_{n}^{(t)}) - \nabla f_n (\xv_n^{(t)};\mathcal{D}_n) \|^2 \leq \sigma_n^2,\nonumber\\ &\Ebb \| {\nabla} f_{n}(\xv_n^{(t)};\mathcal{B}_{n}^{(t)}) \|^2 \leq G^2. \end{align}\]
Theorem 1 quantifies the MSE of the network-wide average model relative to \(\xv^*\); its proof is omitted here.
Theorem 1. Suppose the instantaneous CSI \(\{h_{mn}^{(t)}\}_{\forall m,n,t}\) is given at each iteration. Let the SGD stepsize be \(\eta_t=\frac{2}{\mu(t+a)}\) with a given constant \(a=\Omega(H)>0\), and define \(\overline{\xv}^{(t)}=\frac{1}{N}\sum_{n=1}^{N}\xv_n^{(t)}\) as the network-wide average model computed by Algorithm 2. It follows that2 \[\begin{align} \Ebb \| \overline{\xv}^{(t)} & - \xv^* \|^2 \leq \frac{4 ( \sum_{k}p_k^2 \zeta_k^2 + \frac{1}{D^2}\sum_{n} |\mathcal{D}_n|^2 \sigma_n^2 ) t}{\mu^2 (t+a-2)^{2 } }\nonumber \\ & + \frac{a^2 \| \xv^{(0)} - \xv^* \|^2}{(t+a-2)^{2 }} + \Ocal \left(H^2\frac{\log t}{t^2} \right), \label{theorem32main} \end{align}\tag{13}\] where expectation is over mini-batch sampling and communication noise; and \(\zeta_k^2 \triangleq \Ebb \left\| \widehat \gv_k^{(t)} - \sum_{i\in\mathcal{C}_k}\frac{|\mathcal{D}_i|}{D_{\mathcal{C}_k}}\,{\nabla} f_{i}(\xv_k;\mathcal{B}_{i}) \right \|^2\) is MSE caused by the intra-cluster OTA aggregation (cf. Line 8 of Algorithm 2).
Theorem 1 shows that the dominant error terms decay as \(\Ocal(1/t)\) and are shaped by the OTA aggregation term \(\sum_{k}p_k^2 \zeta_k^{2}\). Under the transceiver design in 7 , this error term depends on the instantaneous CSI \(\{h_{mn}\}\) and, by [11], can be expressed as \[\begin{align} \sum_{k=1}^K p_k^2 \zeta_k^2&\propto \sum_{k=1}^K \left(\min_{i\in\mathcal{C}_k} \frac{|h_{iv_k}^{(t)}|^2}{|\mathcal{D}_i|^2}\right)^{-1}. \end{align}\] Since \(h_{iv_k}^{(t)}\!\sim\!\mathcal{CN}(0,\alpha_{iv_k})\), \(\min_{i}|h_{iv_k}^{(t)}|^{2}/|\mathcal{D}_i|^{2}\) is exponentially distributed with rate \(\sum_{i\in\mathcal{C}_k}|\mathcal{D}_i|^{2}/\alpha_{iv_k}\). Motivated by this, without access to instantaneous CSI, we design the clusters by minimizing the following aggregation MSE statistic subject to the budget constraints in 8 –9 : \[\begin{align} \sum_{k=1}^K \left(\Ebb\left[\min_{i\in\mathcal{C}_k} \frac{|h_{iv_k}^{(t)}|^2}{|\mathcal{D}_i|^2}\right]\right)^{-1}=\sum_{k=1}^K\sum_{i\in \mathcal{C}_k}\frac{|\mathcal{D}_i|^2}{\alpha_{iv_k}}.\label{cut32left32tails} \end{align}\tag{14}\]
We present a low-complexity, sub-optimal solution to the minimization of 14 as follows:
Construct \(\mathcal{H}\) for a given \(K\) and clustering. With fixed \(K\) and \(\{\mathcal{C}_k,v_k\}\), the objective 14 is independent of \(\mathcal{H}\). We therefore build \(\mathcal{H}\) at minimum cost by computing the minimum spanning tree (MST) of the \(K\) heads, which takes \(\Ocal(K^{2}\log K)\) time [12].
Cluster formation for a given \(K\). For any feasible \(K\!\ge\!1\), the constrained clustering of \(\{\mathcal{C}_k,v_k\}\) is NP-hard. We adopt a heuristic method to solve for a sub-optimal solution. When the head vector \(\boldsymbol{v}=[v_1,\dots,v_K]\) is fixed, each node \(i\) can be greedily assigned to the cluster that minimizes \(|\mathcal{D}_i|^{2}/\alpha_{iv_k}\). We then refine \(\boldsymbol{v}\) iteratively via Gibbs sampling. Denote the head selection state in iteration \(r\) as \(\boldsymbol{v}_r\). Consider all feasible neighboring states that differ from \(\boldsymbol{v}_r\) in exactly one head index; we draw a new state \(\boldsymbol{v}_{r+1}\) from the neighboring state with probability proportional to \(e^{-{\rm Obj}(\boldsymbol{v}_r')/\Lambda_r}\), where \({\rm Obj}(\cdot)\) is the value of 14 evaluated at the neighboring state \(\boldsymbol{v}_r'\), and \(\Lambda_r=\rho\Lambda_{r-1}\) (\(0<\rho<1\)) is a decreasing “temperature" parameter. Here, we adopt a cooling schedule for \(\Lambda_r\), which balances the trade-off between exploration and exploitation. Finally, the best state encountered in \(R\) iterations is returned.
Line search over \(K\). Starting from \(K=1\), we iteratively increment \(K\), apply the above two steps, and stop once the budget constraints are violated. Finally, the best feasible solution is selected.
The complete procedure is summarized in Algorithm 3.

Figure 4: Device clustering results and the objective (in dB) in Algorithm 3 versus the value of \(K \in \{2,3, 4, 5\}\). Different colors represent clustering groups \(\{\mathcal{C}_k\}\). The reliable links in \(\mathcal{H}\) is visualized via solid black lines with the link distances marked above..
None
Figure 5: Average training loss (top) and test accuracy (bottom) versus number of iterations for the Fashion-MNIST image classification with the D2D network in Fig. 4..
In this section we examine the proposed approach through simulations. A D2D network with \(N=50\) devices is generated over a square region \([-100~\text{m},100~\text{m}]^{2}\) using a stochastic-block model with four blocks, where the intra- and inter-block connection probabilities are set to \(0.99\) and \(0.01\). The large-scale fading is set to \(\alpha_{mn}=\alpha_0 d^{-\gamma}\) with the reference path loss \(\alpha_0=0\) dB and path-loss exponent \(\gamma=3.76\). The small-scale fading in \(h_{mn}^{(t)}\) is drawn i.i.d. in each training iteration according to 4 . The transmit-power limit is \(P_0=2\) W, and the communication noise power is \(0\) dBW.
Impact of the number of clusters. We first investigate the device clustering performance of Algorithm 3 with budget parameters \(\tau=120\), \(B=530\), \(L_N=50\), and \(L_D=1\) in 8 –9 . Fig. 4 reports the clustering resulting and the objective value in 14 representing the convergence-rate error term versus \(K\). Among all feasible solutions, \(K=4\) minimizes the objective. The plot also shows that the objective is not monotone in \(K\): if \(K\) grows too large, the budget constraints force cluster heads to lie closer together, which subsequently degrades the clustering performance. Note that the line search in Algorithm 3 correctly identifies \(K=4\) as optimal.
Learning performance. We next test the proposed clustered DFL method in Algorithm 2 with \(K=4\) on the Fashion-MNIST classification task [13]. The \(60{,}000\) training images are i.i.d. distributed among the \(N=50\) devices, while the \(10{,}000\) test images are used for evaluating the accuracy of the network-wide average model. We train a \(3\)-layer multi-layer perceptron (MLP) model of \(d=235,146\) parameters with stepsize \(\eta_t=\tfrac{30}{1000+t}\), and local batch size \(100\). We compare the proposed method with the baseline method in [10], which relies on the OTA gossip scheme in 3 for model aggregation. Fig. 4 depicts the training loss and test accuracy. Thanks to the reliable links between cluster heads, the proposed scheme mitigates inter-cluster stragglers and achieves markedly faster convergence than the baseline. Moreover, the proposed method is insensitive to the value of the inter-cluster communication interval \(H\). This observation corroborates the analysis in Theorem 1, where \(H\) does not appear in the leading \(\mathcal{O}(1/t)\) error term of the optimality gap. From a practical standpoint, a larger \(H\) is therefore desirable to reduce communication overhead without sacrificing accuracy.
In this work, we addressed the straggler problem that limits DFL model aggregation efficiency in large-scale D2D networks. We proposed a budget-aware, cluster-centric framework, which identifies sensitive straggling devices and provisions a small set of reliable backhaul links between carefully chosen devices. This design converts decentralized model aggregation into (i) fast device-to-server model aggregation within each cluster, and (ii) low-frequency gossip exchange among cluster heads. We established an \(\Ocal(1/t)\) convergence bound on the optimality gap, which explicitly captures the impact of model aggregation error. We then developed a budget-constrained device-clustering algorithm that minimizes the leading error term without requiring instantaneous CSI. Experiments on image classification verified that the proposed method markedly outperforms the state-of-the-art DFL approach in learning convergence. These results demonstrate that modest, strategically placed link enhancements can unlock substantial efficiency gains for large-scale DFL.