CHMAS: A Coupled Hierarchical Framework for Multi-Agent Reinforcement Learning


Abstract

Multi-agent reinforcement learning (MARL) systems face fundamental challenges in balancing global coordination with local execution across different temporal scales. This paper introduces the Coupled Hierarchical Multi-Agent System (CHMAS), a novel framework that decomposes multi-agent decision-making into centralized strategic planning and distributed tactical execution with bidirectional information flow. The strategic layer integrates all agents’ states with an exclusive global environmental state to generate guidance actions every \(T\) timesteps, while tactical agents execute distributed policies augmented by strategic guidance and local neighborhood observations. Unlike existing hierarchical approaches with unidirectional control, CHMAS establishes a feedback mechanism where accumulated tactical rewards influence strategic objectives through a coupling coefficient \(\lambda\), ensuring strategic plans remain grounded in tactical feasibility. To address the non-stationarity inherent in hierarchical learning, we propose an asynchronous update protocol where strategic parameters update every \(N_f\) tactical episodes, allowing tactical policies to converge to quasi-stationary points between strategic changes. We present both a general bi-level formulation capturing full system dynamics and a tractable additive approximation enabling rigorous analysis. Theoretical analysis proves that this asynchronous scheme achieves \(\mathcal{O}(\log K/\sqrt{K})\) convergence for the strategic layer after \(K\) strategic updates under standard assumptions. Experimental validation in a multi-agent foraging domain demonstrates successful learning of spatially partitioned exploration strategies, with both layers converging stably despite hierarchical coupling.

1 INTRODUCTION↩︎

Multi-agent reinforcement learning (MARL) has emerged as a fundamental paradigm for autonomous coordination in complex environments, with applications spanning robotics [1], traffic management [2], and resource allocation [3]. Despite significant progress, real-world deployment faces a critical challenge: coordinating distributed agents across multiple temporal and spatial scales while maintaining computational tractability.

The curse of dimensionality renders purely centralized approaches intractable. For \(N\) agents with state spaces \(\mathcal{S}_i\) and action spaces \(\mathcal{A}_i\), the joint space grows as \(O(\prod_{i=1}^N |\mathcal{S}_i| \cdot |\mathcal{A}_i|)\) [4]. While distributed methods address scalability through local decision-making [5], they sacrifice coordination quality for long-horizon objectives. Traditional approaches like QMIX [6] and MADDPG [7] leverage centralized training with decentralized execution, but assume all agents access relevant global information during training, which often fails in practice due to partial observability and communication constraints [8].

Furthermore, multi-agent systems inherently operate across multiple temporal scales. Strategic decisions (formation control, task allocation) require global information including environmental states inaccessible to individual agents, operating on extended timescales. Tactical responses (collision avoidance, local adaptation) demand rapid reactions based on immediate observations. While hierarchical reinforcement learning offers a natural framework for multi-scale decision-making [9], [10], existing hierarchical MARL approaches [11], [12] employ unidirectional control where high-level policies dictate low-level actions, missing opportunities for tactical feedback to inform strategic planning.

The interaction between hierarchical levels creates complex learning dynamics. Strategic changes fundamentally alter the environment perceived by tactical agents, violating stationarity assumptions [8]. Conversely, evolving tactical policies change execution characteristics that strategic planners rely upon. Previous hierarchical MARL work [13], [14] addresses this through reward engineering or option learning but requires domain expertise and often fails to generalize.

This paper introduces CHMAS, which addresses these limitations through bidirectional coupling between strategic and tactical layers at different temporal scales. The strategic level maintains centralized access to both agents’ states and an exclusive global environmental state \(\mathcal{S}^{\text{env}}\) capturing system-wide properties essential for coordination. Strategic guidance, generated every \(T\) timesteps, broadcasts to tactical agents who execute distributed policies based on augmented local observations. Tactical rewards accumulate and influence strategic objectives through coupling coefficient \(\lambda\), creating feedback that grounds strategic planning in tactical realities. An asynchronous protocol with strategic updates every \(N_f\) tactical episodes manages non-stationary dynamics, allowing tactical convergence between strategic changes.

Our contributions are threefold. First, we formalize the hierarchical MARL problem with bidirectional coupling and information asymmetry, providing both a general bi-level formulation and tractable additive approximation. Second, we develop the Asynchronous Hierarchical Policy Gradient (AHPG) algorithm managing non-stationarity through temporal separation while maintaining \(\mathcal{O}(\log K/\sqrt{K})\) convergence guarantees for the strategic layer. Third, we validate the framework empirically, demonstrating stable bidirectional learning in a cooperative multi-agent foraging task.

2 PRELIMINARIES↩︎

We consider a cooperative multi-agent system with \(N\) agents operating in a shared environment. The system employs a hierarchical architecture with a strategic layer that accesses global information at low temporal frequency and a tactical layer that operates distributedly at high frequency using only local information.

2.1 Multi-Agent Communication Graph↩︎

The tactical agents interact through a time-invariant communication graph \(\mathcal{G} = (\mathcal{V}, \mathcal{E})\) where \(\mathcal{V} = \{1, 2, \ldots, N\}\) represents the agent set and \(\mathcal{E} \subseteq \mathcal{V} \times \mathcal{V}\) denotes the edge set [15], [16]. For each agent \(i\), the neighborhood \(\mathcal{N}_i = \{j \in \mathcal{V} : (i,j) \in \mathcal{E}\}\) is the set of agents that can directly communicate with agent \(i\). This graph constrains information flow at the tactical level, ensuring distributed execution remains scalable.

2.2 Strategic MDP↩︎

Definition 1 (Strategic MDP). The strategic layer evolves as \(\mathcal{M}^{\text{str}} = (\mathcal{S}^{\text{str}}, \mathcal{A}^{\text{str}}, P^{\text{str}}, R^{\text{str}}, \gamma_{\text{str}})\), where: (i) the state \(\mathcal{S}^{\text{str}} = \mathcal{S}^{\text{env}} \times \prod_{i=1}^N \mathcal{S}_i\) combines the exclusive global environmental state \(\mathcal{S}^{\text{env}}\) (resource distributions, global objectives) with all agent states, so that at epoch \(k\) the strategic state is \(s^{\text{str}}_k = (s^{\text{env}}_k, s_{1,k}, \ldots, s_{N,k})\); (ii) actions \(\mathcal{A}^{\text{str}} = \prod_{i=1}^N \mathcal{A}^{\text{str}}_i\) provide personalized guidance \(a^{\text{str}}_{i,k}\) to each agent, persisting for the interval \([kT,(k+1)T)\); and (iii) the reward integrates global objectives with tactical feedback: \[\label{eq:strategic95reward} R^{\text{str}}(s^{\text{str}}_k, \mathbf{a}^{\text{str}}_k) = R^{\text{global}}(s^{\text{env}}_k, \mathbf{a}^{\text{str}}_k) + \lambda \sum_{t=kT}^{(k+1)T-1} \sum_{i=1}^N r^{\text{tac}}_{i,t},\qquad{(1)}\] where \(\lambda \in [0,1]\) controls the influence of tactical execution quality on strategic evaluation.

The strategic policy \(\pi^{\text{str}}: \mathcal{S}^{\text{str}} \rightarrow \mathcal{P}(\mathcal{A}^{\text{str}})\) maps complete state observations to guidance distributions, encoding high-level coordination logic that individual agents cannot compute from local information alone.

2.3 Tactical Multi-Agent MDP↩︎

Definition 2 (Tactical MAMDP). Given fixed strategic guidance \(\mathbf{a}^{\text{str}}_k\), the tactical layer evolves as \(\mathcal{M}^{\text{tac}} = (\{\mathcal{S}_i\}, \{\mathcal{O}^{\text{tac}}_i\}, \{\mathcal{A}^{\text{tac}}_i\}, P^{\text{tac}}, \{R^{\text{tac}}_i\}, \gamma_{\text{tac}})\), where: (i) each agent \(i\) observes \[\label{eq:tactical95obs} o^{\text{tac}}_{i,t} = \phi_i(s_{i,t}, \{s_{j,t}\}_{j \in \mathcal{N}_i}, a^{\text{str}}_{i,k}) \in \mathcal{O}^{\text{tac}}_i,\qquad{(2)}\] combining local state, neighbor states, and strategic guidance (\(k = \lfloor t/T \rfloor\)); (ii) the joint transition factorizes over the communication graph: \(P^{\text{tac}}(\mathbf{s}_{t+1} | \mathbf{s}_t, \mathbf{a}^{\text{tac}}_t) = \prod_{i=1}^N P_i(s_{i,t+1} | s_{i,t}, \{s_{j,t}, a^{\text{tac}}_{j,t}\}_{j \in \mathcal{N}_i \cup \{i\}})\), maintaining tractability as \(N\) grows [4]; and (iii) the per-agent reward \[\label{eq:tactical95reward} R^{\text{tac}}_i(s_{i,t}, a^{\text{tac}}_{i,t}, \{s_{j,t}, a^{\text{tac}}_{j,t}\}_{j \in \mathcal{N}_i}, a^{\text{str}}_{i,k})\qquad{(3)}\] enables local optimization and neighbor coordination while incorporating strategic objectives.

2.4 Bidirectional Coupling↩︎

Unlike traditional top-down hierarchies [11], CHMAS establishes feedback in both directions. Downward coupling (strategic \(\to\) tactical): guidance \(a^{\text{str}}_{i,k}\) shapes tactical behavior through observation augmentation ?? , reward modulation ?? , and temporal persistence over \(T\) timesteps. Upward coupling (tactical \(\to\) strategic): accumulated rewards \(R^{\text{tac,cum}}_k = \sum_{t=kT}^{(k+1)T-1} \sum_{i=1}^N r^{\text{tac}}_{i,t}\) influence strategic rewards via \(\lambda\) in ?? , so that strategic decisions account for tactical feasibility. Setting \(\lambda = 0\) recovers a purely top-down hierarchy. The framework maintains asymmetric information access: the strategic layer observes complete state \(s^{\text{str}}_k\) while tactical agents access only \(\mathcal{N}_i\), with strategic guidance serving as an information bottleneck [17].

3 PROBLEM FORMULATION↩︎

3.1 General Coupled Formulation↩︎

The bidirectional coupling creates a bi-level optimization structure [18], [19]. The strategic objective: \[\label{eq:strategic95objective} J^{\text{str}}(\theta^{\text{str}} | \{\theta^{\text{tac}}_i\}) = \mathbb{E}\!\left[\sum_{k=0}^{\infty} \gamma^k_{\text{str}} R^{\text{str}}(s^{\text{str}}_k, \mathbf{a}^{\text{str}}_k)\right]\tag{1}\] depends on tactical policies through ?? . Each tactical objective: \[\label{eq:tactical95objective} J^{\text{tac}}_i(\theta^{\text{tac}}_i | \theta^{\text{str}}) = \mathbb{E}\!\left[\sum_{t=0}^{\infty} \gamma^t_{\text{tac}} r^{\text{tac}}_{i,t}\right]\tag{2}\] depends on strategic parameters through ?? and ?? .

Definition 3 (Coupled Hierarchical Problem). The coupled optimization problem is the bi-level program: \[\begin{align} \max_{\theta^{\text{str}}} \quad & J^{\text{str}}(\theta^{\text{str}} | \{\theta^{\text{tac}*}_i(\theta^{\text{str}})\}) \\ \text{s.t.} \quad & \theta^{\text{tac}*}_i(\theta^{\text{str}}) = \arg\max_{\theta^{\text{tac}}_i} J^{\text{tac}}_i(\theta^{\text{tac}}_i | \theta^{\text{str}}), \quad \forall i \in \mathcal{V}, \end{align}\] where the strategic optimization accounts for tactical best responses.

Remark 1. We operate in the cooperative MARL setting where agents share a common team objective. Accordingly, the per-agent best responses \(\theta^{\text{tac}*}_i\) are aligned with the joint team optimum, and no equilibrium-seeking across agents is required at the tactical level [4].

This Stackelberg structure is computationally intractable in general [20], motivating the approximation below.

3.2 Additive Approximation for Tractable Optimization↩︎

Definition 4 (Additive Hierarchical Problem). Under the additive approximation, the global objective is: \[\label{eq:additive} J^{\text{global}}(\theta^{\text{str}}, \{\theta^{\text{tac}}_i\}) = J^{\text{str}}(\theta^{\text{str}}) + \sum_{i=1}^N J^{\text{tac}}_i(\theta^{\text{tac}}_i),\qquad{(4)}\] where each term is evaluated with fixed parameters of the other layer.

Remark 2. Equation ?? is an approximation of the true coupled objective: the cross-layer dependencies in 12 are not eliminated but handled implicitly through the asynchronous update schedule. The approximation error is controlled by the tactical optimality gap \(\delta_k\) which will be defined in Assumption 2 at each epoch, which vanishes as learning progresses.

The separable gradient structure, \[\label{eq:gradient} \nabla_{\theta} J^{\text{global}} = \bigl[\nabla_{\theta^{\text{str}}} J^{\text{str}},\; \nabla_{\theta^{\text{tac}}} J^{\text{tac}}\bigr]^T,\tag{3}\] allows updates to be treated as coordinate steps in the joint parameter space, enabling rigorous convergence analysis via two-timescale stochastic approximation [21], [22].

4 ASYNCHRONOUS LEARNING ALGORITHM↩︎

4.1 Algorithm Design Principles↩︎

AHPG leverages three principles for stable hierarchical learning.

Temporal separation. Strategic decisions occur every \(T\) timesteps and persist throughout the interval. Strategic parameters follow a decaying step size schedule \(\eta^{\text{str}}_k = \alpha/\sqrt{k}\) and update only every \(N_f\) tactical episodes, operating at a slower effective timescale than tactical learning. This allows tactical policies to converge to near-optimal responses between strategic changes, transforming the non-stationary hierarchical problem into alternating stationary subproblems.

Structured information flow. Downward flow provides coordination signals through strategic guidance without requiring centralized control at each timestep. Upward flow aggregates tactical performance \(R^{\text{tac}}_k = \sum_{t,i} r^{\text{tac}}_{i,t}\) into strategic rewards, keeping high-level plans grounded in execution realities. The strategic layer maintains exclusive access to \(s^{\text{env}}_k\), preserving the information asymmetry that enables scalable distributed execution.

Distributed tactical learning. Each tactical agent maintains its own policy \(\pi^{\text{tac}}_i\) and experience buffer \(\mathcal{D}^{\text{tac}}_i\), updating parameters using only local gradients without inter-agent communication during learning, so that tactical computation scales linearly with \(N\).

4.2 Algorithm Specification↩︎

Algorithm 1 presents the complete AHPG procedure. Strategic guidance is generated at the start of each new strategic epoch (every \(N_f\) episodes), and strategic parameters are updated at the end of each epoch using the accumulated tactical performance. Tactical parameters update every episode.

Figure 1: Asynchronous Hierarchical Policy Gradient (AHPG)

During the \(N_f\) episodes between strategic updates, tactical agents operate in a stationary environment with fixed guidance, converging toward near-optimal policies. This ensures strategic updates observe stable tactical performance that accurately estimates the tactical contribution to strategic rewards.

For learning components, the strategic layer employs centralized policy gradient algorithms (e.g., PPO [23]) leveraging complete state access, while tactical agents use distributed policy gradient or actor-critic variants operating on local observations. The framework naturally extends to continuous control through continuous guidance vectors and to partially observable settings through recurrent policies.

5 CONVERGENCE ANALYSIS↩︎

We establish convergence for AHPG under the additive approximation \(J^{\text{global}}(\theta) = J^{\text{str}}(\theta^{\text{str}}) + \sum_{i=1}^N J^{\text{tac}}_i(\theta^{\text{tac}}_i)\).

Assumption 1.

  1. Smoothness: \(J^{\text{str}}\) and each \(J^{\text{tac}}_i\) are \(L\)-smooth with \(L\)-Lipschitz gradients.

  2. Boundedness: \(J^{\text{str}} \le J^{\text{str}*}\) and \(J^{\text{tac}}_i \le J^{\text{tac}*}_i\) for all \(i\).

  3. Bounded variance: Stochastic gradient estimates satisfy \(\mathbb{E}[\|g^{\text{str}}_k - \nabla J^{\text{str}}_k\|^2] \le \sigma^2_{\text{str}}\) and \(\mathbb{E}[\|g^{\text{tac}}_{i,e} - \nabla J^{\text{tac}}_i\|^2] \le \sigma^2_{\text{tac}}\).

  4. Biased strategic gradient: \(\mathbb{E}[g^{\text{str}}_k \mid \theta^{\text{str}}_k] = \nabla J^{\text{str}}(\theta^{\text{str}}_k) + b_k\), where \(b_k\) captures the error from evaluating strategic gradients with non-converged tactical parameters.

Assumption 2.

  1. PL condition: Each \(J^{\text{tac}}_i\) satisfies the \(\mu\)-PL inequality: \(\|\nabla J^{\text{tac}}_i(\theta)\|^2 \ge 2\mu(J^{\text{tac}*}_i - J^{\text{tac}}_i(\theta))\).

  2. Warm-start: At the beginning of each strategic epoch \(k\), tactical agents are initialized from the converged state of the previous epoch. The initial optimality gap satisfies \(\mathbb{E}[J^{\text{tac}*}_i - J^{\text{tac}}_i(\theta^{\text{tac}}_{i,kN_f})] \le \delta_k\), where \(\delta_k = \mathcal{O}(\eta^{\text{str}}_k) = \mathcal{O}(1/\sqrt{k})\), reflecting the magnitude of the strategic parameter change between successive epochs.

  3. Round-robin updates: Within each tactical phase, agents update sequentially with others’ parameters held fixed, so that the PL condition applies to each per-agent landscape independently.

Remark 3. The PL condition in Assumption 2 holds for log-linear (softmax) policies in tabular settings [24] and locally for over-parameterized neural networks near stationary points. The round-robin assumption is a standard device to decouple per-agent landscapes; it can be relaxed to simultaneous updates at the cost of a more involved coupling analysis.

Assumption 3 (Coupling structure). The strategic gradient is \(L_b\)-Lipschitz in the tactical parameters: \[\|\nabla J^{\text{str}}(\theta^{\text{str}}, \boldsymbol{\theta}^{\text{tac}}) - \nabla J^{\text{str}}(\theta^{\text{str}}, \boldsymbol{\theta}^{\text{tac}*})\| \le L_b \|\boldsymbol{\theta}^{\text{tac}} - \boldsymbol{\theta}^{\text{tac}*}\|.\]

Theorem 1 (AHPG Convergence). Under Assumptions 13, with \(\eta^{\text{str}}_k = \alpha/\sqrt{k}\) where \(\alpha \le 1/(8L)\), \(\eta^{\text{tac}} = \beta/\sqrt{K}\) where \(\beta < 1/(2L)\), and \(N_f = c\sqrt{K}\) for a constant \(c > 0\), after \(K\) strategic updates AHPG achieves: \[\min_{k \in \{1,\ldots,K\}} \mathbb{E}[\|\nabla J^{\text{str}}(\theta^{\text{str}}_k)\|^2] = \mathcal{O}\!\left(\frac{\log K}{\sqrt{K}}\right),\] \[\frac{1}{KN_f}\sum_{k=1}^{K}\sum_{e=kN_f}^{(k+1)N_f-1} \sum_{i=1}^N \mathbb{E}[\|\nabla J^{\text{tac}}_i(\theta^{\text{tac}}_{i,e})\|^2] = \mathcal{O}\!\left(\frac{1}{\sqrt{K}}\right).\]

Proof. Strategic layer. By \(L\)-smoothness of \(J^{\text{str}}\) and the update \(\theta^{\text{str}}_{k+1} = \theta^{\text{str}}_k + \eta^{\text{str}}_k g^{\text{str}}_k\): \[J^{\text{str}}_{k+1} \ge J^{\text{str}}_k + \eta^{\text{str}}_k \langle\nabla J^{\text{str}}_k,\, g^{\text{str}}_k\rangle - \tfrac{L}{2}(\eta^{\text{str}}_k)^2 \|g^{\text{str}}_k\|^2.\] Taking expectation with \(\mathbb{E}[g^{\text{str}}_k] = \nabla J^{\text{str}}_k + b_k\), we bound the two terms. For the inner product, applying Young’s inequality \(\langle a, b\rangle \ge -\tfrac{1}{4}\|a\|^2 - \|b\|^2\): \[\mathbb{E}[\langle\nabla J^{\text{str}}_k, g^{\text{str}}_k\rangle] \ge \tfrac{3}{4}\|\nabla J^{\text{str}}_k\|^2 - \|b_k\|^2.\] For the squared norm, using \(\mathbb{E}[\|g^{\text{str}}_k\|^2] \le \|\nabla J^{\text{str}}_k + b_k\|^2 + \sigma^2_{\text{str}} \le 2\|\nabla J^{\text{str}}_k\|^2 + 2\|b_k\|^2 + \sigma^2_{\text{str}}\): \[\begin{align} \mathbb{E}[J^{\text{str}}_{k+1}] &\ge J^{\text{str}}_k + \!\left(\tfrac{3}{4}\eta^{\text{str}}_k - L(\eta^{\text{str}}_k)^2\right) \|\nabla J^{\text{str}}_k\|^2 \\ &\quad - \!\left(\eta^{\text{str}}_k + L(\eta^{\text{str}}_k)^2\right)\|b_k\|^2 - \tfrac{L}{2}(\eta^{\text{str}}_k)^2 \sigma^2_{\text{str}}. \end{align}\] Since \(\eta^{\text{str}}_k \le 1/(8L)\), we have \(L(\eta^{\text{str}}_k)^2 \le \eta^{\text{str}}_k/8\), so \(\tfrac{3}{4}\eta^{\text{str}}_k - L(\eta^{\text{str}}_k)^2 \ge \tfrac{5}{8}\eta^{\text{str}}_k \ge \tfrac{1}{8}\eta^{\text{str}}_k\) and \(\eta^{\text{str}}_k + L(\eta^{\text{str}}_k)^2 \le \tfrac{9}{8}\eta^{\text{str}}_k \le 2\eta^{\text{str}}_k\). Rearranging and telescoping over \(k=1,\ldots,K\): \[\begin{align} \label{eq:str95telescope} \frac{1}{8}\sum_{k=1}^K \eta^{\text{str}}_k \mathbb{E}[\|\nabla J^{\text{str}}_k\|^2] &\le \Delta J^{\text{str}} + 2\sum_{k=1}^K \eta^{\text{str}}_k \|b_k\|^2 \nonumber\\ &\quad + \frac{L\sigma^2_{\text{str}}}{2} \sum_{k=1}^K (\eta^{\text{str}}_k)^2, \end{align}\tag{4}\] where \(\Delta J^{\text{str}} = J^{\text{str}*} - J^{\text{str}}_1\).

Bias bound. By Assumption 3, \(\|b_k\|^2 \le L_b^2 \|\boldsymbol{\theta}^{\text{tac}}_{(k+1)N_f} - \boldsymbol{\theta}^{\text{tac}*}\|^2\). Running \(N_f = c\sqrt{K}\) SGD steps under the PL condition (Assumption 2) with step size \(\beta/\sqrt{K}\) and warm-start gap \(\delta_k\) yields contraction factor \(\rho = (1 - 2\mu\beta/\sqrt{K})^{c\sqrt{K}} \to e^{-2\mu\beta c} < 1\), so the residual tactical deviation satisfies \(\|b_k\|^2 = \mathcal{O}(\delta_k + 1/\sqrt{K}) = \mathcal{O}(1/\sqrt{k})\).

Completing the strategic bound. With \(\eta^{\text{str}}_k = \alpha/\sqrt{k}\): \[\begin{align} \sum_{k=1}^K \eta^{\text{str}}_k\|b_k\|^2 &= \mathcal{O}\!\left(\alpha\sum_{k=1}^K \tfrac{1}{\sqrt{k}} \cdot \tfrac{1}{\sqrt{k}}\right) = \mathcal{O}(\log K),\\ \sum_{k=1}^K(\eta^{\text{str}}_k)^2 &= \alpha^2\sum_{k=1}^K \tfrac{1}{k} = \mathcal{O}(\log K). \end{align}\] Dividing 4 by \(\sum_{k=1}^K \eta^{\text{str}}_k = \mathcal{O}(\alpha\sqrt{K})\) gives \(\min_k \mathbb{E}[\|\nabla J^{\text{str}}_k\|^2] = \mathcal{O}(\log K/\sqrt{K})\).

Tactical layer. Within epoch \(k\), \(\theta^{\text{str}}_k\) is fixed, so \(J^{\text{tac}}_i\) is stationary. Standard SGD under the \(\mu\)-PL condition over \(N_f\) steps gives per-epoch average gradient norm \(\mathcal{O}(\delta_k + 1/\sqrt{K})\). Averaging over \(K\) epochs and using \(K^{-1}\sum_{k=1}^K\delta_k = K^{-1}\sum_{k=1}^K\mathcal{O}(k^{-1/2}) = \mathcal{O}(K^{-1/2})\), the overall tactical average is \(\mathcal{O}(1/\sqrt{K})\), completing the proof. ◻

Remark 4. Theorem 1 is established for policy gradient implementations. The DQN-based experiments in Section 6 serve as an empirical proof of concept; extending the convergence analysis to Q-learning variants is left for future work.

6 EXPERIMENTAL VALIDATION↩︎

6.1 Experiment Setting↩︎

We evaluate CHMAS on a 25\(\times\)​25 GridWorld environment with \(N=4\) agents performing cooperative foraging, instantiating the key features of hierarchical multi-agent coordination: information asymmetry, temporal abstraction, and bidirectional coupling.

Figure 2: The 25\times​25 GridWorld environment. Four agents (A0-blue, A1-orange, A2-purple, A3-green) operate within their assigned 9\times​9 regions (dashed boxes) determined by strategic guidance. Colored areas indicate visited cells tracked in the global state s^{\text{env}}, while agents pursue local resources (red apple).

The global environmental state \(s^{\text{env}} \in \{0,1\}^{625}\) tracks cell visitation, representing system-wide coverage that individual agents cannot observe. The strategic layer accesses the complete state \(s^{\text{str}} = (s^{\text{env}}, s_1, s_2, s_3, s_4)\), while each tactical agent \(i\) observes only local information \(o^{\text{tac}}_i \in \mathbb{R}^{20}\) comprising its position, nearby resources, and strategic guidance.

Strategic actions assign each agent a 9\(\times\)​9 region center \(a^{\text{str}}_i = (g^x_i, g^y_i) \in [4,20]^2\) at the start of each epoch (\(T=5\) timesteps). Tactical agents select discrete actions \(a^{\text{tac}}_i \in \{0,1,2,3,4,5\}\) (movement, collection, wait).

The reward structure implements bidirectional coupling: tactical agents receive \(r^{\text{tac}}_{i,t} = r^{\text{collect}}_{i,t} - 0.01 - 0.01 \cdot \mathbf{1}[\text{outside region}]\), while the strategic layer exclusively receives global coverage reward \(R^{\text{global}}_k = 0.01 \times |\text{newly visited cells}|\) computed from \(s^{\text{env}}\). Strategic rewards integrate both: \(R^{\text{str}}_k = R^{\text{global}}_k + \lambda \sum_{t,i} r^{\text{tac}}_{i,t}\) with \(\lambda = 0.5\).

The implementation is publicly available at https://github.com/EricDmWang/Hierarchy_RL.git.

6.2 Results and Analysis↩︎

We implement the tactical layer using Multi-Agent DQN (MA-DQN) with parameter sharing, dueling architecture, double Q-learning, and prioritized experience replay. The strategic layer employs a centralized DQN on the global state. Training runs for 400 episodes with 5 independent seeds.

a
b
c
d

Figure 3: Tactical layer learning dynamics.. a — Tactical Q-values, b — Tactical returns, c — Tactical critic loss, d — Tactical policy loss

Figure 3 shows the tactical layer’s learning progression. Q-values (a) grow rapidly before stabilizing near episode 300. Returns (b) improve from approximately \(-80\) to \(-10\), showing agents learn efficient resource collection while respecting region constraints. Critic loss (c) and policy loss (d) peak early during exploration and then decrease as policies stabilize.

a
b
c
d

Figure 4: Strategic layer learning dynamics with \(N_f=5\).. a — Strategic Q-values, b — Strategic returns, c — Strategic critic loss, d — Strategic actor loss

Figure 4 shows the strategic layer’s asynchronous learning with updates every \(N_f=5\) episodes. Strategic Q-values (a) improve from \(-8\) to near \(0\), and strategic returns (b) increase from \(-10\) to approximately \(15\), demonstrating successful optimization of global coverage objectives. The delayed onset of learning (first 10–20 episodes) reflects the asynchronous protocol allowing tactical policies to stabilize before strategic updates commence. Low variance across runs confirms algorithmic stability.

6.3 Qualitative Analysis↩︎

Figure 5 shows the execution trajectory of a trained AHPG model over a complete episode. The strategic layer learns to assign non-overlapping regions that collectively tile the grid, initially placing agents in corners and then adapting assignments to unexplored areas. Agents balance region compliance with opportunistic resource collection. By frame (e), full resource collection with maintained spatial coordination validates the effectiveness of bidirectional hierarchical control.

a
b
c
d
e

Figure 5: Evolution of hierarchical coordination during execution. Agents operate within their assigned 9\(\times\)​9 regions (dashed boxes) while collectively achieving grid coverage. Red apples represent resources collected opportunistically.. a — \(t=0\), b — \(t=6\), c — \(t=9\), d — \(t=16\), e — \(t=19\)

7 CONCLUSION↩︎

This paper introduced CHMAS, a hierarchical MARL framework based on bidirectional coupling between strategic and tactical layers. By enabling strategic planning with exclusive global state access while maintaining distributed tactical execution, CHMAS provides a new architectural paradigm for coordinating agents across multiple temporal scales. The asynchronous update protocol, with a decaying strategic step size schedule, manages the non-stationarity inherent in hierarchical learning and admits \(\mathcal{O}(\log K/\sqrt{K})\) convergence guarantees for the strategic layer.

The bidirectional coupling mechanism offers a principled approach to constrained reinforcement learning: strategic guidance encodes global constraints while tactical feedback ensures feasibility, with \(\lambda\) providing tunable constraint satisfaction.

References↩︎

[1]
T. T. Nguyen, N. D. Nguyen, and S. Nahavandi, “Deep reinforcement learning for multiagent systems: A review of challenges, solutions, and applications,” IEEE transactions on cybernetics, vol. 50, no. 9, pp. 3826–3839, 2020.
[2]
M. A. Wiering et al., “Multi-agent reinforcement learning for traffic light control,” in Machine learning: Proceedings of the seventeenth international conference (ICML’2000), 2000, pp. 1151–1158.
[3]
G. Qu, Y. Lin, A. Wierman, and N. Li, “Scalable multi-agent reinforcement learning for networked systems with average reward,” Advances in Neural Information Processing Systems, vol. 33, pp. 2074–2086, 2020.
[4]
K. Zhang, Z. Yang, and T. Başar, “Multi-agent reinforcement learning: A selective overview of theories and algorithms,” Handbook of Reinforcement Learning and Control, pp. 321–384, 2021.
[5]
K. Zhang, Z. Yang, H. Liu, T. Zhang, and T. Basar, “Fully decentralized multi-agent reinforcement learning with networked agents,” in International conference on machine learning, 2018, pp. 5872–5881.
[6]
T. Rashid, M. Samvelyan, C. S. De Witt, G. Farquhar, J. Foerster, and S. Whiteson, “Monotonic value function factorisation for deep multi-agent reinforcement learning,” Journal of Machine Learning Research, vol. 21, no. 178, pp. 1–51, 2020.
[7]
R. Lowe, Y. I. Wu, A. Tamar, J. Harb, P. Abbeel, and I. Mordatch, “Multi-agent actor-critic for mixed cooperative-competitive environments,” in Advances in neural information processing systems, 2017, vol. 30.
[8]
P. Hernandez-Leal, B. Kartal, and M. E. Taylor, “A survey and critique of multiagent deep reinforcement learning,” Autonomous Agents and Multi-Agent Systems, vol. 33, no. 6, pp. 750–797, 2019.
[9]
R. S. Sutton, D. Precup, and S. Singh, “Between MDPs and semi-MDPs: A framework for temporal abstraction in reinforcement learning,” Artificial Intelligence, vol. 112, no. 1–2, pp. 181–211, 1999.
[10]
P.-L. Bacon, J. Harb, and D. Precup, “The option-critic architecture,” in Proceedings of the AAAI conference on artificial intelligence, 2017, vol. 31.
[11]
S. Pateria, B. Subagdja, A. Tan, and C. Quek, “Hierarchical reinforcement learning: A comprehensive survey,” ACM Computing Surveys, vol. 54, no. 5, pp. 1–35, 2021.
[12]
A. S. Vezhnevets et al., “Feudal networks for hierarchical reinforcement learning,” in International conference on machine learning, 2017, pp. 3540–3549.
[13]
T. D. Kulkarni, K. Narasimhan, A. Saeedi, and J. Tenenbaum, “Hierarchical deep reinforcement learning: Integrating temporal abstraction and intrinsic motivation,” in Advances in neural information processing systems, 2016, vol. 29.
[14]
O. Nachum, S. S. Gu, H. Lee, and S. Levine, “Data-efficient hierarchical reinforcement learning,” Advances in Neural Information Processing Systems, vol. 31, 2018.
[15]
R. Olfati-Saber, J. A. Fax, and R. M. Murray, “Consensus and cooperation in networked multi-agent systems,” Proceedings of the IEEE, vol. 95, no. 1, pp. 215–233, 2007.
[16]
A. Nedic, A. Ozdaglar, and P. A. Parrilo, “Constrained consensus and optimization in multi-agent networks,” IEEE Transactions on Automatic Control, vol. 55, no. 4, pp. 922–938, 2010.
[17]
N. Tishby, F. C. Pereira, and W. Bialek, “The information bottleneck method,” arXiv preprint physics/0004057, 2000.
[18]
B. Colson, P. Marcotte, and G. Savard, “An overview of bilevel optimization,” Annals of Operations Research, vol. 153, no. 1, pp. 235–256, 2007.
[19]
A. Sinha, P. Malo, and K. Deb, “A review on bilevel optimization: From classical to evolutionary approaches and applications,” IEEE Transactions on Evolutionary Computation, vol. 22, no. 2, pp. 276–295, 2018.
[20]
M. Hong, H.-T. Wai, Z. Wang, and Z. Yang, “A two-timescale stochastic algorithm framework for bilevel optimization: Complexity analysis and application to actor-critic,” SIAM Journal on Optimization, vol. 33, no. 1, pp. 147–180, 2023.
[21]
V. S. Borkar, “Stochastic approximation with two time scales,” Systems & Control Letters, vol. 29, no. 5, pp. 291–294, 1997.
[22]
Y. Nesterov, “Efficiency of coordinate descent methods on huge-scale optimization problems,” SIAM Journal on Optimization, vol. 22, no. 2, pp. 341–362, 2012.
[23]
J. Schulman, F. Wolski, P. Dhariwal, A. Radford, and O. Klimov, “Proximal policy optimization algorithms,” arXiv preprint arXiv:1707.06347, 2017.
[24]
A. Agarwal, S. M. Kakade, J. D. Lee, and G. Mahajan, “Optimality and approximation with policy gradient methods in markov decision processes,” in Conference on learning theory, 2020, pp. 64–66.

  1. This work was supported by the National Science Foundation under Grant ECCS-2129949.↩︎

  2. The authors are with the Department of Electrical and Computer Engineering, University of California, Riverside, CA 92521, USA. {wdong025, jxu150, yzhan831, wei.ren}@ucr.edu.↩︎