The manifold of unitary and symmetric matrices: characterization, Riemannian optimization and application to BD-RIS design
Pith reviewed 2026-05-08 10:24 UTC · model grok-4.3
The pith
Riemannian optimization on unitary symmetric matrices improves BD-RIS MIMO performance.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The geometry of the set of unitary symmetric matrices admits a Riemannian structure with explicit tangent space, retraction, and geodesics, which supports the design of efficient optimization algorithms for BD-RIS scattering matrices in MIMO systems that achieve superior sum-gain and rates compared to existing techniques while exploiting low-rank structure when the RIS is large.
What carries the argument
The manifold of unitary symmetric matrices U_s, with its derived tangent space, simple retraction, and closed-form geodesics that enable line-search and geodesic phase-optimization algorithms.
Load-bearing premise
The proposed retraction and closed-form geodesics remain computationally tractable and numerically stable for practical BD-RIS MIMO dimensions, and local optima found by the algorithms are close enough to the global optimum for the non-convex objectives.
What would settle it
Running the Riemannian line-search and phase-optimization algorithms on a BD-RIS MIMO setup with more elements than antennas and verifying if the optimized matrix has rank equal to the antenna count, or comparing the achieved performance metrics against baseline methods to see if they are higher.
Figures
read the original abstract
This paper proposes and analyzes Riemannian optimization algorithms on the manifold of unitary and symmetric matrices, denoted ${\cal {U}}_s$, which naturally models the scattering matrices of passive and reciprocal devices such as beyond-diagonal reconfigurable intelligent surfaces (BD-RISs). Despite its relevance, the geometry of ${\cal {U}}_s$ has remained largely unexplored, and existing BD-RIS optimization methods either ignore the symmetry constraint or rely on costly Takagi-based parameterizations. We first provide a rigorous geometric characterization of ${\cal {U}}_s$, deriving its tangent space, a simple retraction, and closed-form expressions for geodesics. Building on these results, we develop two Riemannian manifold optimization (MO) algorithms tailored to ${\cal {U}}_s$: a line-search (LS) based scheme and a phase-optimization (PO) update along geodesics. We then apply the proposed framework to BD-RIS-assisted multiple-input multiple-output (MIMO) links, addressing sum-gain maximization, rate maximization, and minimum mean-square error problems, where they outperform existing approaches. Furthermore, we show that when the number of BD-RIS elements exceeds the total number of antennas, the optimal scattering matrix is low-rank, which motivates and enables efficient low-rank variants of the proposed algorithms.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims to rigorously characterize the manifold of unitary and symmetric matrices (denoted U_s) by deriving its tangent space, a simple retraction, and closed-form geodesic expressions. Building on this, it develops two Riemannian manifold optimization algorithms (a line-search scheme and a phase-optimization update along geodesics) and applies them to BD-RIS-assisted MIMO systems for sum-gain maximization, rate maximization, and MMSE problems, reporting outperformance over existing methods. It further shows that when the number of BD-RIS elements exceeds the total number of antennas, the optimal scattering matrix is low-rank, motivating efficient low-rank algorithm variants.
Significance. If the geometric derivations hold and the algorithms prove tractable, this work offers a specialized Riemannian framework for symmetric unitary constraints that arise in reciprocal passive devices. The closed-form primitives and low-rank insight could improve scalability and performance in BD-RIS MIMO optimization compared with Takagi factorizations or unconstrained relaxations, with potential impact on practical system design.
major comments (3)
- Abstract and low-rank section: the claim that 'the optimal scattering matrix is low-rank' when the number of BD-RIS elements exceeds the total number of antennas requires precise mathematical clarification. Unitary matrices are full-rank by definition, so the statement must specify whether it refers to a reduced-dimensional submanifold, a low-rank factorization preserving symmetry and unitarity, or the rank of an effective operator; without this, the motivation for the efficient variants is unclear.
- Algorithm development and complexity analysis: the retraction and closed-form geodesics rely on matrix functions (exponentials, square roots) whose O(N^3) cost and numerical stability for practical BD-RIS sizes (N=64–256) are not quantified. The manuscript should provide explicit complexity counts, runtime benchmarks, and floating-point error analysis to confirm tractability.
- Optimization results sections: for the non-convex sum-gain and rate objectives, the paper should supply evidence that the Riemannian local optima are sufficiently close to global solutions. This could be addressed by reporting statistics over multiple random initializations or comparisons against convex relaxations or global solvers.
minor comments (2)
- Notation consistency: the manifold symbol U_s and related objects (tangent space, retraction) should be introduced once with clear definitions and used uniformly.
- Experimental tables: performance comparisons should report standard deviations or confidence intervals alongside mean gains to allow assessment of variability.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable suggestions. We have carefully considered each major comment and provide point-by-point responses below, along with our plans for revisions to the manuscript.
read point-by-point responses
-
Referee: Abstract and low-rank section: the claim that 'the optimal scattering matrix is low-rank' when the number of BD-RIS elements exceeds the total number of antennas requires precise mathematical clarification. Unitary matrices are full-rank by definition, so the statement must specify whether it refers to a reduced-dimensional submanifold, a low-rank factorization preserving symmetry and unitarity, or the rank of an effective operator; without this, the motivation for the efficient variants is unclear.
Authors: We appreciate the referee's observation regarding the potential ambiguity in our low-rank claim. Upon reflection, the statement in the abstract and low-rank section refers to the fact that, when the number of BD-RIS elements exceeds the total number of antennas, there exists an optimal scattering matrix that admits a low-rank factorization while preserving the required symmetry and unitarity constraints. Specifically, the optimization can be restricted to a reduced-dimensional submanifold parameterized by lower-rank factors. This provides the motivation for developing efficient low-rank variants of our algorithms. We will revise the abstract and the corresponding section to include this precise mathematical clarification. revision: yes
-
Referee: Algorithm development and complexity analysis: the retraction and closed-form geodesics rely on matrix functions (exponentials, square roots) whose O(N^3) cost and numerical stability for practical BD-RIS sizes (N=64–256) are not quantified. The manuscript should provide explicit complexity counts, runtime benchmarks, and floating-point error analysis to confirm tractability.
Authors: We agree that quantifying the computational aspects is important for assessing practicality. The proposed retraction and geodesic computations involve matrix square roots and exponentials, which have a complexity of O(N^3) per operation. In the revised manuscript, we will include explicit complexity analyses for the individual primitives and the complete algorithms. Furthermore, we will add runtime benchmarks and floating-point error analyses based on our existing simulations for BD-RIS sizes up to N=256, demonstrating that the methods remain tractable and numerically stable for the considered problem dimensions. revision: yes
-
Referee: Optimization results sections: for the non-convex sum-gain and rate objectives, the paper should supply evidence that the Riemannian local optima are sufficiently close to global solutions. This could be addressed by reporting statistics over multiple random initializations or comparisons against convex relaxations or global solvers.
Authors: We recognize the value of providing evidence on the quality of the local optima for these non-convex problems. To address this, we will augment the optimization results sections with statistics (such as average performance and standard deviation) obtained from multiple random initializations of our Riemannian algorithms. This will illustrate the consistency of the solutions. Where computationally feasible, we will also include comparisons with convex relaxations on smaller-scale instances to support that our local solutions are close to global optima. revision: yes
Circularity Check
No significant circularity: geometric primitives derived from first principles
full rationale
The paper's core contribution is a direct characterization of the manifold U_s (tangent space, retraction, closed-form geodesics) followed by standard Riemannian algorithms (LS and PO) applied to BD-RIS objectives. These steps do not reduce to fitted parameters renamed as predictions, self-citations as load-bearing premises, or ansatzes smuggled via prior work. The low-rank property when N exceeds antenna count follows from the unitarity/symmetry constraints and dimension counting rather than being presupposed. Performance comparisons are external to the derivation chain. The work is self-contained against external benchmarks with no evidence of the enumerated circular patterns.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The set of unitary symmetric matrices forms a smooth embedded submanifold of the unitary group
- standard math Standard Riemannian optimization theory (line search, geodesic updates) applies directly once the manifold geometry is known
Reference graph
Works this paper leans on
-
[1]
P. A. Absil, R. Mahony, and R. Sepulchre,Optimization Algorithms on Matrix Manifolds. Princeton University Press, 2008
2008
-
[2]
Boumal,An Introduction to Optimization on Smooth Manifolds
N. Boumal,An Introduction to Optimization on Smooth Manifolds. Cambridge University Press, 2023. [Online]. Available: https://www. nicolasboumal.net/book
2023
-
[3]
The geometry of algorithms with orthogonality constraints,
A. Edelman, T. A. Arias, and S. T. Smith, “The geometry of algorithms with orthogonality constraints,”SIAM Journal on Matrix Analysis and Applications, vol. 20, no. 2, pp. 303–353, 1998
1998
-
[4]
Ramírez, I
D. Ramírez, I. Santamaría, and L. Scharf,Coherence: In Signal Pro- cessing and Machine Learning. Springer Nature, 2023
2023
-
[5]
Precoder design for massive MIMO downlink with matrix manifold optimization,
R. Sun, C. Wang, A. Lu, X. Gao, and X.-G. Xia, “Precoder design for massive MIMO downlink with matrix manifold optimization,”IEEE Trans. Signal Process., vol. 72, pp. 1065–1080, 2024
2024
-
[6]
Steepest descent algo- rithms for optimization under unitary matrix constraint,
T. E. Abrudan, J. Eriksson, and V . Koivunen, “Steepest descent algo- rithms for optimization under unitary matrix constraint,”IEEE Trans. Signal Process., vol. 56, no. 3, pp. 1134–1147, 2008
2008
-
[7]
Conjugate gradient algorithm for optimization under unitary matrix constraint,
T. Abrudan, J. Eriksson, and V . Koivunen, “Conjugate gradient algorithm for optimization under unitary matrix constraint,”Signal Processing, vol. 89, no. 9, pp. 1704–1714, 2009. [Online]. Available: https://www.sciencedirect.com/science/article/pii/S0165168409000814
2009
-
[8]
projUNN: efficient method for training deep networks with unitary matrices,
B. T. Kiani, R. Balestriero, Y . LeCun, and S. Lloyd, “projUNN: efficient method for training deep networks with unitary matrices,” in Proceedings of the 36th International Conference on Neural Information Processing Systems, ser. NIPS ’22, Red Hook, NY , USA, 2022
2022
-
[9]
Short notes: Completion of a symmetric unitary matrix,
R. F. Mathis, “Short notes: Completion of a symmetric unitary matrix,” SIAM Review, vol. 11, no. 2, pp. 261–263, 1969. [Online]. Available: http://www.jstor.org/stable/2028117
-
[10]
The threefold way. algebraic structure of symmetry groups and ensembles in quantum mechanics,
F. J. Dyson, “The threefold way. algebraic structure of symmetry groups and ensembles in quantum mechanics,”J. Math. Phys., vol. 3, no. 6, pp. 1199–1215, 1962
1962
-
[11]
D. M. Pozar,Microwave Engineering, 4th Ed.John Wiley and Sons, 2011
2011
-
[12]
A tutorial on beyond- diagonal reconfigurable intelligent surfaces: Modeling, architectures, system design and optimization, and applications,
H. Li, M. Nerini, S. Shen, and B. Clerckx, “A tutorial on beyond- diagonal reconfigurable intelligent surfaces: Modeling, architectures, system design and optimization, and applications,”IEEE Communica- tions Surveys & Tutorials, vol. 28, pp. 4086–4126, 2026
2026
-
[13]
Beyond diagonal reconfigurable intelli- gent surfaces: From transmitting and reflecting modes to single-, group- , and fully-connected architectures,
H. Li, S. Shen, and B. Clerckx, “Beyond diagonal reconfigurable intelli- gent surfaces: From transmitting and reflecting modes to single-, group- , and fully-connected architectures,”IEEE Trans. Wireless Commun., vol. 22, no. 4, pp. 2311–2324, 2022
2022
-
[14]
Modeling and architecture design of reconfigurable intelligent surfaces using scattering parameter network analysis,
S. Shen, B. Clerckx, and R. Murch, “Modeling and architecture design of reconfigurable intelligent surfaces using scattering parameter network analysis,”IEEE Trans. Wireless Commun., vol. 21, pp. 1229–1243, 2022
2022
-
[15]
Capacity maximization for MIMO channels assisted by beyond-diagonal RIS,
E. Björnson and T. Demir, “Capacity maximization for MIMO channels assisted by beyond-diagonal RIS,” in19th European Conference on Antennas and Propagation, Stockholm, Sweden, 2025
2025
-
[16]
MIMO channel shaping and rate maximization using beyond-diagonal RIS,
Y . Zhao, H. Li, B. Clerckx, and M. Franceschetti, “MIMO channel shaping and rate maximization using beyond-diagonal RIS,”IEEE Trans. on Signal Processing, vol. 73, pp. 4397–4414, 2025
2025
-
[17]
E. Illi and M. Qaraqe, “Beyond-diagonal RIS for enhanced secrecy and sensing gains in secure ISAC networks: An optimization framework,” arXiv: 10.48550/arXiv.2604.04480, 2026
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2604.04480 2026
-
[18]
A low-complexity beamforming design for beyond-diagonal RIS aided multi-user networks,
T. Fang and Y . Mao, “A low-complexity beamforming design for beyond-diagonal RIS aided multi-user networks,”IEEE Comm. Letters, vol. 28, no. 1, pp. 203–207, 2024
2024
-
[19]
MIMO ca- pacity maximization with beyond-diagonal RIS,
I. Santamaria, M. Soleymani, E. Jorswieck, and J. Gutiérrez, “MIMO ca- pacity maximization with beyond-diagonal RIS,” inIEEE International Workshop on Signal Processing Advances in Wireless Communications (SPAWC), Lucca, Italy, 2024, pp. 936–940
2024
-
[20]
Statistical CSI- enabled optimization for beyond diagonal SIM,
Q. Xia, J. Zhang, K. Xu, S. Ma, S. Jin, and C. Yuen, “Statistical CSI- enabled optimization for beyond diagonal SIM,”IEEE Wireless Comm. Letters, vol. 14, no. 9, pp. 2972–2976, 2025
2025
-
[21]
Interfer- ence minimization in beyond-diagonal RIS-assisted MIMO interference channels,
I. Santamaria, M. Soleymani, E. Jorswieck, and J. Gutierrez, “Interfer- ence minimization in beyond-diagonal RIS-assisted MIMO interference channels,”IEEE Open Journal of Vehicular Technology, vol. 6, pp. 1005–1017, 2025. 12
2025
-
[22]
On an algebraic problem related to an analytic theorem of Carathéodory and Fejér and on an allied theorem of Landau,
T. Takagi, “On an algebraic problem related to an analytic theorem of Carathéodory and Fejér and on an allied theorem of Landau,”Jpn. J. Math. Trans. Abstr., vol. 1, pp. 83–93, 1924
1924
-
[23]
Sur les matrices hypohermitiennes et sur les matrices unitaires,
L. Autonne, “Sur les matrices hypohermitiennes et sur les matrices unitaires,”Ann. Univ. Lyon, vol. 38, pp. 1–77, 1915
1915
-
[24]
Reciprocal beyond-diagonal reconfigurable intelligent surface (BD-RIS): Scattering matrix design via manifold optimization,
M. Fidanovski, I. A. M. Sandoval, H. S. Rou, G. T. F. de Abreu, and E. Björnson, “Reciprocal beyond-diagonal reconfigurable intelligent surface (BD-RIS): Scattering matrix design via manifold optimization,”
-
[25]
Available: https://arxiv.org/abs/2509.20246
[Online]. Available: https://arxiv.org/abs/2509.20246
-
[26]
RIS optimization on the complex circle manifold for interference mitigation in interference channels,
M. A. ElMossallamy, K. G. Seddik, W. Chen, L. Wang, G. Y . Li, and Z. Han, “RIS optimization on the complex circle manifold for interference mitigation in interference channels,”IEEE Trans. on Veh. Tech., vol. 70, no. 6, pp. 6184–6189, 2021
2021
-
[27]
Sum-rate maximization for RIS-assisted integrated sensing and communication systems with manifold optimization,
E. Shtaiwi, H. Zhang, A. Abdelhadi, A. Swindlehurst, Z. Han, and H. V . Poor, “Sum-rate maximization for RIS-assisted integrated sensing and communication systems with manifold optimization,”IEEE Trans. on Communications, vol. 71, pp. 4909–4923, 2023
2023
-
[28]
Manifold-based optimizations for RIS-aided massive MIMO systems,
W. D. S. Junior, D. W. M. Guerra, J. C. M. Filho, T. Abrão, and E. Hossain, “Manifold-based optimizations for RIS-aided massive MIMO systems,”IEEE Open Journal of the Communications Society, vol. 5, pp. 7913–7940, 2024
2024
-
[29]
Riemannian optimization on the manifold of unitary and symmetric matrices with application to BD-RIS-assisted systems,
I. Santamaria, M. Soleymani, E. Jorswieck, J. Gutierrez, and C. Beltran, “Riemannian optimization on the manifold of unitary and symmetric matrices with application to BD-RIS-assisted systems,” inIEEE In- ternational Conference on Acoustics, Speech and Signal Processing (ICASSP), Barcelona, Spain, 2026
2026
-
[30]
Postnikov,Lie groups and Lie algebras
M. Postnikov,Lie groups and Lie algebras. MIR Moscow: Lectures in Geometry 5, 1994
1994
-
[31]
Cayley differential unitary space-time codes,
B. Hassibi and B. M. Hochwald, “Cayley differential unitary space-time codes,”IEEE Trans. on Information Theory, vol. 48, no. 6, pp. 1485– 1503, Jun. 2002
2002
-
[32]
B. C. Hall,Lie Groups, Lie Algebras, and Representations: An Elemen- tary Introduction. Springer, 2015
2015
-
[33]
SNR maximization in beyond diagonal RIS-assisted single and multiple antenna links,
I. Santamaria, M. Soleymani, E. Jorswieck, and J. Gutiérrez, “SNR maximization in beyond diagonal RIS-assisted single and multiple antenna links,”IEEE Signal Proc. Letters, vol. 30, pp. 923–926, 2023
2023
-
[34]
Some metric inequalities in the space of matrices,
K. Fan and A. J. Hoffman, “Some metric inequalities in the space of matrices,”Proceedings of the American Mathematical Society, vol. 6, no. 1, pp. 111–116, 1955. [Online]. Available: http://www.jstor.org/stable/2032662
-
[35]
Sur quelques propriétés des déterminants gauches,
A. Cayley, “Sur quelques propriétés des déterminants gauches,”J. Reine Angew. Math., vol. 32, pp. 119–123, 1846
-
[36]
Special paraunitary matrices, cayley transform, and multidimensional orthogonal filter banks,
J. Zhou, M. N. Do, and J. Kovacevic, “Special paraunitary matrices, cayley transform, and multidimensional orthogonal filter banks,”IEEE Trans. on Image Processing, vol. 15, no. 2, pp. 511–519, Feb. 2006
2006
-
[37]
Beyond diagonal reconfig- urable intelligent surfaces utilizing graph theory: Modeling, architecture design, and optimization,
M. Nerini, S. Shen, H. Li, and B. Clerckx, “Beyond diagonal reconfig- urable intelligent surfaces utilizing graph theory: Modeling, architecture design, and optimization,”IEEE Trans. on Wireless Commun., vol. 23, no. 8, pp. 9972–9985, Aug. 2024
2024
-
[38]
Closed-form global optimization of beyond diagonal reconfigurable intelligent surfaces,
M. Nerini, S. Shen, and B. Clerckx, “Closed-form global optimization of beyond diagonal reconfigurable intelligent surfaces,”IEEE Trans. Wireless Commun., vol. 23, no. 2, pp. 1037–1051, 2024
2024
-
[39]
D. P. Bertsekas,Nonlinear Programming, 2nd ed. Athena Scientific, 1999
1999
-
[40]
Convergence of a block coordinate descent method for nondifferentiable minimization,
P. Tseng, “Convergence of a block coordinate descent method for nondifferentiable minimization,”Journal of Optimization Theory and Applications, vol. 109, pp. 475–494, 2001
2001
-
[41]
A unified convergence analysis of block successive minimization methods for nonsmooth optimization,
M. Razaviyayn, M. Hong, and Z. Q. Luo, “A unified convergence analysis of block successive minimization methods for nonsmooth optimization,”SIAM J. Opt., vol. 23, no. 2, pp. 1126–1153, 2013
2013
-
[42]
Optimizing power consumption, energy efficiency, and sum-rate using beyond diago- nal RIS—a unified approach,
Y . Zhou, Y . Liu, H. Li, Q. Wu, S. Shen, and B. Clerckx, “Optimizing power consumption, energy efficiency, and sum-rate using beyond diago- nal RIS—a unified approach,”IEEE Trans. on Wireless Communications, vol. 23, no. 7, pp. 7423–7438, 2024
2024
-
[43]
Physically consistent modeling of wireless links with reconfigurable intelligent sur- faces using multiport network analysis,
J. A. Nossek, D. Semmler, M. Joham, and W. Utschick, “Physically consistent modeling of wireless links with reconfigurable intelligent sur- faces using multiport network analysis,”IEEE Wireless Communications Letters, vol. 13, no. 8, pp. 2240–2244, Aug. 2024
2024
-
[44]
A novel Q-stem connected architecture for beyond-diagonal reconfigurable intelligent surfaces,
X. Zhou, T. Fang, and Y . Mao, “A novel Q-stem connected architecture for beyond-diagonal reconfigurable intelligent surfaces,” inIEEE Inter- national Conference on Communications (ICC), 2025, pp. 6880–6885
2025
-
[45]
X. Zhou, T. Fang, Y . Mao, and B. Clerckx, “Generalized beyond-diagonal RIS architectures: Theory and design via structured-oriented symmetric unitary projection,”arXiv: https://arxiv.org/pdf/2509.17804v2, 2025
-
[46]
Beyond-diagonal RIS in multiuser MIMO: Graph theoretic modeling and optimal architectures with low complex- ity,
Z. Wu and B. Clerckx, “Beyond-diagonal RIS in multiuser MIMO: Graph theoretic modeling and optimal architectures with low complex- ity,”IEEE Trans. on Inf. Theory, vol. 71, no. 11, pp. 8506–8523, 2025
2025
-
[47]
Optimal symmetric low-rank BD-RIS configuration maximizing the determinant of a MIMO link
I. Santamaria, M. Soleymani, J. Gutiérrez, and E. Jorswieck, “Optimal symmetric low-rank BD-RIS configuration maximizing the determinant of a MIMO link,”arXiv: https://doi.org/10.48550/arXiv.2604.09335, 2026
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2604.09335 2026
-
[48]
Gradient of mutual information in linear vector Gaussian channels,
D. P. Palomar and S. Verdú, “Gradient of mutual information in linear vector Gaussian channels,”IEEE Trans. Inf. Theory, vol. 52, no. 1, pp. 141–154, 2006
2006
-
[49]
Capacity characterization for intelligent reflecting surface aided MIMO communication,
S. Zhang and R. Zhang, “Capacity characterization for intelligent reflecting surface aided MIMO communication,”IEEE J. Sel. Areas Commun., vol. 38, no. 8, pp. 1823–1838, 2020
2020
-
[50]
MMSE design of RIS-aided communications with spatially-correlated channels and electromagnetic interference,
W. X. Long, M. Moretti, A. Abrardo, L. Sanguinetti, and R. Chen, “MMSE design of RIS-aided communications with spatially-correlated channels and electromagnetic interference,”IEEE Trans. Wireless Com- mun., vol. 23, no. 11, pp. 16 992–17 006, 2024
2024
-
[51]
MIMO transceiver design via majorization theory,
D. P. Palomar and Y . Jiang, “MIMO transceiver design via majorization theory,”Found. Trends Commun. Inf. Theory, vol. 3, no. 4, 2006
2006
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.