Bilayer Intelligent Omni-Surface-Assisted Full-Duplex Systems: Simultaneous Self-Interference Cancellation and Sum Rate Maximization

Recently, reconfigurable intelligent surfaces (RISs) have attracted increasing attentions in the design of full-duplex (FD) systems due to their novel capability of propagation environment reconfiguration. However, in conventional RIS-assisted FD systems, the beamforming for self-interference cancellation (SIC) and sum rate maximization (SRM) are highly coupled during RIS optimization, which significantly degrades the system performance. To tackle this issue, we exploit a novel bilayer intelligent omni-surface (BIOS) structure in FD systems. Compared with the conventional RIS designs, the BIOS provides independent beams on both sides, thus enabling more flexible achievement of SRM and SIC. For the BIOS-assisted FD system, we first formulate an optimization problem to achieve SRM and efficient SIC simultaneously. Then, we exploit the relationship between the SRM and mean square error (MSE), and propose a weighted MSE minimization with SIC algorithm to solve the problem. Specifically, we jointly design the beamforming at the base station and the BIOS with manifold optimization while guaranteeing an SIC constraint. Furthermore, we theoretically derive a lower band for the BIOS size required for efficient SIC in FD systems. Simulation results indicate that the BIOS outperforms the conventional RIS designs in FD systems, and verify the accuracy of the derived lower bound for the BIOS size.


Introduction
With the explosive growth of wireless communication devices, the scarcity of available spectrum resources has become increasingly acute.Consequently, effectively utilizing existing spectrum resources has emerged as a pivotal factor in advancing wireless communication technology [1][2][3].Recently, in-band full-duplex (FD) systems have drawn notable interest in the design of next-generation wireless communication systems due to their capability of transmitting and receiving data on the same frequency band concurrently.Compared to traditional half-duplex (HD) systems, such as time division duplex or frequency division duplex, FD systems can significantly enhance spectrum utilization and reduce the time delay of wireless communication systems [4,5].Additionally, FD systems can significantly improve confidentiality while maintaining the same data rate [4].Therefore, FD technology is critical for communication in ultra-dense network environments in the future.
Despite the advantages mentioned above, a prominent problem in FD systems is the inevitable self-interference (SI) emanating from the transmitter to the receiver, as its high intensity significantly interferes with uplink (UL) signals.To tackle this issue, three typical categories of SI cancellation (SIC) methods are discussed in the literature: transmissiondomain processing, analog-domain processing, and digital-domain signal processing [4].The transmission-domain SIC scheme involves suppressing SI through hardware design, such as physical isolation, shielding materials, and antennas with different polarization directions, thereby minimizing the power of SI signals before they reach the receiver [6][7][8].The analog-domain SIC scheme focuses on suppressing SI in the radio frequency (RF) chain by adding isolators or designing circuits to generate signals that counteract the SI signals [9][10][11][12][13].In [9], the authors added an additional RF chain at the transmitter, connected to the receiver.They estimated the SI channel and used the results to digitally preprocess the transmitted signal, generating a counteracting signal with the same amplitude but opposite phase to the SI, thus canceling it at the receiver.Similarly, in [10], an active RF chain was utilized to generate a signal compensating for the inevitable leakage from passive analog SIC.The digital-domain SIC scheme employs digital signal processing techniques on baseband signals to cancel residual SI, exploiting the spatial degrees of freedom provided by multi-antenna systems.Typically, digital SIC methods involve beamforming techniques such as zero-forcing (ZF), minimum mean square error (MMSE), and least squares [14][15][16][17].In [14], the authors used ZF precoding at the transmitter and maximum ratio combining at the receiver to effectively cancel SI.In [15], the authors employed the MMSE algorithm to jointly design transmitting and receiving beamforming, aiming to minimize the mean square error (MSE) and reduce the SI.However, despite numerous attempts, these SIC methods come with high hardware expenses, substantial complexity, and constrained performance.Thus, the reconfigurable intelligent surface (RIS) is considered a promising technology and has attracted increasing attention in the design of FD systems.Specifically, RIS is a two-dimensional (2D) metasurface consisting of programmable, low-cost, and passive elements [18,19].Due to its flexibility and low-cost nature, RIS has been extensively applied in wireless communication systems [20][21][22][23][24].
For the FD systems, RIS has been proven a promising technology to tackle the SIC and thus is widely investigated in the literature [25][26][27][28][29].In [26], intelligent reflecting surface (IRS), as the most famous type of RIS, was exploited in FD systems to achieve efficient SIC.Specifically, the IRS smartly reconfigured the cascade signal emitted from the base station (BS) transmitter (BS Tx ) to the IRS and then reflected to the BS receiver (BS Rx ) to create destructive interference with the SI signal.Moreover, taking into account the quantization noise of analog-to-digital converters (ADCs), the authors of [27] considered a more practical scenario and proposed to jointly design the downlink (DL) precoding matrix and the IRS coefficients.However, due to hardware limitations, the IRS can only reflect the signal back to the same side [30].Therefore, it can only service user equipment (UE) on the same side of the BS.To tackle this constraint, the authors of [28] first introduced intelligent omnisurface (IOS) to FD systems, which is named simultaneously transmitting and reflecting reconfigurable intelligent surface (STAR-RIS) in some studies [31].Compared with the IRS, the IOS can simultaneously reflect and refract the incident signal to both sides [30,31].Thus, it can extend the wireless coverage to the full space while suppressing the SI signal in FD systems.
Although previous studies have focused on the performance of RIS-assisted FD systems, there remains a notable gap in the literature regarding the simultaneous SIC and sum rate maximization (SRM) of DL and UL.Ref. [28] only addressed SI constraints targeting DL rate optimization, while [25] assumed perfect SIC when maximizing the sum rate (SR).However, minimizing the residual SI is crucial for effective air interface SIC as excessive SI will lead to the saturation of ADCs, preventing it from functioning properly [7].Thus, it is necessary to investigate beamforming optimization for the SRM under practical SI constraints in RIS-assisted FD systems.Moreover, it is notable that although current IOS models offer both reflection and refraction capabilities, they cannot independently adjust reflecting and refracting phase shifts [30][31][32].This limitation results in coupled beamforming on both sides, which leads to a serious performance loss in the FD system as the IOS is expected to simultaneously suppress the SI signal and enhance the system data rate.
To deal with these issues, in this paper, we introduce a novel bilayer intelligent omnisurface (BIOS) structure to simultaneously suppress the SI and enhance the system SR in FD systems.This BIOS structure deploys two neighboring IOS layers, and thus can ensure independent beamforming design on both sides [33].Unlike [33], which utilizes both reflected and refracted beams for communication, our design employs the reflected signals for SIC and the refracted signal for communication, thus achieving both objectives simultaneously.Our approach is grounded in a practical scenario that considers the power constraint of the SI signal.The main contributions are summarized as follows:

•
We introduce a novel BIOS in a multi-user FD system.With the unique bilayer structure, the BIOS can provide flexible beamforming on both sides and thus can simultaneously improve the SR and cancel the SI more effectively than the conventional RIS designs.

•
For the BIOS-assisted multi-user FD system, we jointly optimize the beamforming at the BS and BIOS to maximize the SR of UL and DL.We ensure effective over-theair SIC by constraining the residual SI power below a threshold, marking a distinct approach from that of other studies.To solve this problem, we establish the proof of equivalence between the SRM and a weighted MSE minimization (WMMSE) and then propose a novel SIC-WMMSE algorithm.

•
To further guide the practical deployment of BIOS-assisted FD systems, we theoretically derive a lower bound of the BIOS size required for efficient SIC, which is proven to be quite close to the accurate size required in practical systems.

•
We provide simulation results to show that the BIOS-assisted FD system significantly outperforms that of the conventional IOS due to its capability of independent beamforming on both sides.
Notations: Matrices and vectors are denoted by boldface upper case and boldface lower case, respectively.(•) T , (•) H , (•) −1 , and (•) * , respectively, denote the transpose, Hermitian transpose, inverse, and conjugate.[•] ij and [•] i denote the (i, j)-th element of a matrix and the i-th element of a vector.C M×N denotes the space of M × N complex-valued matrices.The trace, l 1 norm, and Frobenius norm are denoted by Tr(•), | • | 1 , and • F , respectively.R{•} and E(•), respectively, denote the real-part and expectation operation.• denotes Hadamard products.Furthermore, I N denotes the N × N identity matrix.diag(a) is a diagonal matrix with the elements of a on its main diagonal and diag(A) represents the extraction of the diagonal of A. CN (0, K) denotes the circularly symmetric complex Gaussian distribution with zero mean and covariance matrix K.

BIOS Structure
As mentioned in Section 1, hardware limitations prevent IOS from independently controlling the phases of refracted and reflected signals, resulting in symmetric beams on both sides.However, in the proposed FD system, we aim to utilize the reflected signals to cancel SI and for the refracted signals to communicate.Thus, the symmetric beams will degrade the system performance.Specifically, the cascaded reflected signal, which is emitted from the BS Tx into the IOS and then reflected to the BS Rx , serves exclusively for SIC while the refracted signals aim to maximize SR.Moreover, refraction beamforming must align DL signals with the DL UEs and enhance UL reception, complicating the optimization problem due to increased coupling and non-convexity.Consequently, performance tradeoffs are necessary to find an optimal balance.
To mitigate the constraints, we propose a BIOS-assisted FD system as shown in Figure 1, which can achieve simultaneous SIC and SRM by flexibly controlling the reflection and refraction beams on both sides of the BIOS [33,34].The BIOS consists of two neighboring IOSs (referred to as IOS 1 and IOS 2 ).Situated proximately to the BS, IOS 1 functions like a traditional IOS, managing both reflection and refraction, while IOS 2 is set to the fulltransmitting mode.With this arrangement, the reflected cascaded signals can be smartly designed by controlling the IOS 1 to suppress SI signals.Furthermore, the SR can be maximized by independently controlling the reflection and refraction beams on both sides of the BIOS.Thus, the BIOS can achieve the dual-function target of SIC and SRM in FD systems better than the conventional RIS designs.Specifically, compared to traditional IRS, BIOS allows the BS to achieve greater service coverage by reflecting and refracting signals into the full space.Moreover, the capability of independent beamforming on both sides enables separate and effective management for SIC and SRM.Compared to single-layer IOS, the addition of IOS 2 enables independent beams on both sides, allowing for more flexible adjustment.Additionally, compared to multilayer refractive RIS investigated in the previous studies [35,36], the BIOS is distinguished by its specific dual function of simultaneous and decoupled reflection and refraction, which enhances its performance in FD systems.
Without loss of generality, we assume that both IOS 1 and IOS 2 are uniform square planar arrays consisting of M = M y × M z elements, where M y = M z .Thus, the effective reflecting and refracting coefficient matrices of BIOS can be respectively presented as follows: where Φ 1 ∈ C M×M and Φ 2 ∈ C M×M denote the diagonal coefficient matrices of IOS 1 and IOS 2 .ξ r and ξ t are the respective fractions of reflection and refraction signal power relative to the total power of the incoming signal for IOS 1 , satisfying the power constraint ξ r + ξ t = 1.(Following the assumption in [32], we presume negligible penetration loss for IOS and BIOS.Nevertheless, the penetration loss will be an interesting and important issue that deserves further investigation.)Similarly, r and t denote those for IOS 2 with r + t = 1.According to the setup, r = 0 and t = 1.G ∈ C M×M denotes the channel between IOS 1 and IOS 2 .

System Model
Consider a narrowband multi-input single-output (MISO) FD system comprising a BS, a BIOS, K D DL UEs, and K U UL UEs in Figure 1.We assume that all the UEs are single-antenna, BS Tx employs N T transmit antennas, and BS Rx has N R receive antennas.The BIOS is deployed significantly close to the BS to ensure the reflected signals are strong enough to efficiently suppress the SI signals, so we assume that there is no line-of-sight (LoS) path between the BS and UEs [28].
For this FD system, the DL signal from the BS is denoted as denotes the original symbol vector, and x U,l is the symbol transmitted by the l-th UL UE with the power constraint E{|x U,l | 2 } = P U .The channels between the BS Tx -IOS 1 , BS Tx -BS Rx , BS Rx -IOS 1 , IOS 2 -UE D,k , and IOS 2 -UE U,l are denoted, respectively, as , and h U,l ∈ C M×1 .Therefore, the DL received signal at the k-th DL UE and the UL received signal at the BS Rx can be respectively represented as follows: where n MU D,k represents the interference from other DL UEs for the k-th DL UE and n MU U,l represents the interference from other UL UEs for the l-th UL UE. n SI is the SI at the BS.n D,k denotes the Gaussian noise at the k-th DL UE with mean zero and variance σ 2 D , and n U denotes the noise at the BS satisfying the distribution of n U ∼ CN 0, σ 2 U I .

Channel Model
Before problem formulation, we first present detailed channel modeling.As the BIOS is deployed at the BS to achieve efficient SIC, H TB , H TR , and H BR are considered nearfield channels without non-line-of-sight (NLoS) paths.Thus, considering the properties of spherical waves, we model the LoS path gain from the n-th receive antenna of BS Tx to the m-th element of IOS 1 as follows [37]: where a is the size of scattering elements, λ is the carrier wavelength, and d m,n is the distance between the n-th antenna of BS Tx and the m-th element of IOS 1 .θ r m,n and φ r m,n , respectively, denote the elevation and azimuth angle of the received signal at the m-th element of IOS 1 .Moreover, G TB = 2 and F TB (θ, φ) = |cos θ|, respectively, denote parameters related to the equivalent antenna gain and the normalized power radiation pattern between BS Tx and IOS 1 .
Similarly, the channel gain of H BR , H TR , and G can be modeled as follows: where Considering the NLoS paths, h D,k and h U,l can be modeled as Rician channels.Specifically, the channel gain between the k-th DL UE and m-th element of IOS 2 can be modeled as follows: where κ is the Rician factor.h LoS m and h NLoS m , respectively, denote the LoS and NLoS path gain from the m-th element of IOS 2 to the k-th DL UE. d m denotes the distance between the m-th element of IOS 2 and the k-th DL UE. α 0 denotes the pass loss of the reference distance d 0 with the path loss exponent β. h s m 1 ,n ∼ CN (0, 1)denotes the cumulative effect of a large number of scattered paths.Moreover, G UE = 2, F UE (θ, φ) = | cos θ|.Similarly, h U,l can be modeled in the same way.

Problem Formulation
According to (2) and ( 3), the data rate of the k-th DL UE can be expressed as Similarly, the UL data rate of the l-th UL UE is given by where Ĥ = In this paper, we jointly optimize the beamforming design at the BS and BIOS to maximize the SR of UL and DL while effectively managing SIC.It is worth noting that excessive residual SI will cause the saturation of ADC, preventing it from functioning properly.However, most of the existing works are based on perfect SIC.To tackle this issue, we take into account the hardware limitations and then impose a constraint to maintain the residual SI power below a threshold, which is different from other papers.Specifically, the problem can be formulated as follows: max Tr FF H ≤ P s .
Constraint (13b) maintains the residual SI power below a threshold to achieve efficient SIC.Constraints (13c) and (13d) result from the constant modulus constraint of the BIOS scattering elements and the power limitation at the BS, respectively.

Beamforming Design for BIOS-FD Systems
Problem ( 13) is difficult to solve directly due to the severely coupled variables and non-convex constraints.Thus, we propose an SIC-WMMSE method to transform the primal problem into an equivalent weighted MSE minimization problem, and then employ alternating optimization (AO) to solve it.

SIC-WMMSE Approach
Due to the multiple coupled variables and non-convex constraints, problem ( 13) is difficult to solve directly.However, by leveraging the relationship between MSE and data rate, we convert the primal problem (13) into an equivalent WMMSE problem as follows: min where a D,k and a U,l are the weighting coefficients for the k-th DL UE and the l-th UL UE, respectively.u * D,k is the receive coefficient at the k-th DL UE and u H U,l is the receive combiner for the l-th UL UE at the BS Rx , respectively.e D,k and e U,l , respectively, denote the MSEs of the k-th DL UE and the l-th UL UE as follows: where ĥH D,k = represents the equivalent channel of the l-th UL UE for simplification and ĤU = ĥU,1 , • • • , ĥU,K U .The proof of the equivalence between problems ( 13) and ( 14) can be found in Appendix A. For a more concise representation, problem ( 14) can be rewritten as follows by transforming the scalars into matrices: min where respectively, denote diagonal matrices containing the MSEs of each DL and UL UE.Although we have transformed the SRM problem into an equivalent WMMSE problem (17), it is still difficult to solve it directly due to the highly non-convex constraints.Thus, we propose an efficient way to deal with the residual SI power constraint (13b).Specifically, move the SI suppression term into the objective function with a penalty factor ρ, resulting in a new optimization problem as follows: min where F is the power of the SI signal.In order to facilitate the optimization of each variable, we substitute all variables into objective function (18) to derive a more detailed expression of f , which can be expressed as follows: where

Optimize A
According to the WMMSE design approach [28] and the derivation in Appendix A, A D , A U , U D , and U U all have closed-form solutions when other variables are fixed, which can be expressed as follows:

Optimize F with Fixed
When A D , A U , U D , U U , Φ 1 , and Φ 2 are fixed, problem ( 18) is reduced into a subproblem for the optimization of F with the power constraint, which can be rewritten as follows by ignoring constant terms: It can be observed that sub-problem ( 27) is convex.Thus, with other parameters fixed, the optimal F can be obtained by utilizing the Karush-Kuhn-Tucker (KKT) conditions as where δ is the KKT coefficient whose optimal solution can be obtained through a one-dimensional search based on the power constraint.The follow two subsections focus on the optimization of BIOS, i.e., Φ 1 , Φ 2 .

Optimize
In this subsection, we will optimize Φ 1 by fixing other variables.By further defin- , the objective function in (18) related to Φ 1 can be simplified as follows: Due to the constant modulus constraint of BIOS, sub-problem (29) remains highly nonconvex.Thus, we propose to apply the manifold optimization (MO) algorithm to address this challenge.Similar to the gradient descent method, the MO algorithm aims to find the most rapid descent direction on the manifold which also minimally deviates from it, namely, the Riemann conjugate gradient (RCG).It can be derived that the RCG is the projection of the Euclidean gradient in the tangent space, which is the set of directions that least deviate from the manifold surface.Consequently, the key of MO lies in the derivation of the Euclidean gradient.To obtain the gradient of f 1 (Φ 1 ), we first extract the main diagonal elements of Φ 1 containing useful information to form vector φ 1 , i.e., φ 1 = diag(Φ 1 ), and transform the objective function f 1 (Φ 1 ) into f 1 (φ 1 ) as follows: where A and B are equivalent matrices as and equality (a) holds due to the following properties: Tr(C H PCQ) = c H (P Q T )c, Tr(PC) = p T c for any matrices P, Q, and diagonal matrix C, with p = diag(P) and c = diag(C).Thus, the Euclidean gradient of f 1 (φ 1 ) can be obtained as Thus, the RCG can be derived to facilitate gradient descent with an appropriate step size.Subsequently, the optimized φ 1 is projected back onto the manifold surface to ensure the constant modulus constraint.Consequently, the optimal φ 1 can be obtained by iteratively descending along the RCG until a predefined stopping criterion is satisfied.

It is important to notice that the term Tr
TB results from the SI power constraint, while the other terms are associated with the SR.Consequently, the optimization of Φ 1 takes into account SRM and SIC simultaneously.
Similarly, we can also employ MO to solve (33).Extract the main diagonal elements of Φ 2 to form vector φ 2 ; the objective f 2 (Φ 2 ) can be transformed into f 2 (φ 2 ) as follows: where φ 2 = diag(Φ 2 ), and C and D are equivalent matrices as Thus, the Euclidean gradient of f 2 (φ 2 ) can be obtained as It is worth noting that, despite the similar form, the optimization of Φ 1 and Φ 2 has different meanings for this system.As indicated in the UL rate expression (12), the SI is tied to Φ 1 while not directly related to the optimization of Φ 2 .Therefore, we can ignore the penalty item related to SIC in (19) during the Φ 2 optimization.Furthermore, as the BIOS can independently control the reflection and refraction beams on both sides, the additional freedom enables Φ 2 to focus on the beamforming for SRM without affecting the SIC quality.
The proposed SIC-WMMSE algorithm for solving the sub-problem is summarized in Algorithm 1.It is worth noting that the application of the KKT conditions and the MO algorithm enables a reduction in the objective function for each variable within each iteration.When the decrement in the weighted MSE between two iterations is below a predefined threshold, the iteration is terminated.Since the objective function decreases after each step of updating, Algorithm 1 can converge to a locally optimal solution.To analyze its computational complexity, we count the number of complex multiplications of the key steps in one iteration turn.Assuming that N T = N T = N, K D = K U = K, and M > N > K, the complexity can be summarized as follows: updating 2 and F (0) and set i = 0.

Analysis of BIOS Size
To further guide the practical deployment of BIOS-assisted FD systems, a lower bound of BIOS size required for efficient SIC is derived here for the single-input single-output (SISO) scenario.Specifically, the channels between the BS Tx -BS Rx , BS Tx -IOS 1 , and BS Rx -IOS 1 are denoted as h TR , h TB , and h BR , respectively.Theorem 1.To achieve efficient SIC, the number of elements in one layer of the BIOS must exceed a threshold M that satisfies the following inequality: where d BR m , d TB m , respectively, denote the distances from the m-th element of the IOS 1 to the BS Rx and BS Tx , and d BB and d TR denote the distance between BS-BIOS and BS Rx -BS Tx .
Proof of Theorem 1.To ensure effective SIC in the SISO scenario, the residual SI should satisfy the following power constraint: Notably, for effective SIC, the threshold γ should be minimized, and the power constraint (38) should be universally valid for any s D .This requires that the residual interference approaches zero, which can be expressed as To clarify, the necessary condition for efficient SIC is that the modulus of the cascaded channel √ ξ r h H BR Φ 1 h TB should exceed that of the SI channel h TR , which can be formulated as follows: We define g(M) as the maximum gain of the cascade channel generated by BIOS with M scattering elements in one layer, which is given by It can be observed that maximizing the cascade channel gain occurs when the phase shift has an evident closed-form solution as By substituting the channel settings, we have According to the derivations mentioned above, the lower bound of the element number can be derived through the inequality constraint The proof is thus completed.
Theorem 1 provides a lower bound of BIOS size required for SIC, which is of practical importance in the design of FD systems.As will be shown in Section 5, the derived lower bound is rather close to the required BIOS size in the simulation.

Simulation Results
This section presents numerical results to illustrate the performance evaluation of the proposed BIOS-assisted FD system.For comparison, an IOS-assisted FD system is employed as a benchmark, which is characterized by identical phase matrices for reflecting and refracting with ξ r = ξ t = 1/2.Throughout simulations, the centers of the BS and BIOS are set at the origin of a three-dimensional coordinate system.The distance between two neighboring IOSs of the BIOS is set as d II = 0.05 m.The DL and UL UEs are randomly distributed within two respective circles of 3m radius centered at (20, −10, 0) and (10,10,0).Considering the NLoS paths, h D,k and h U,l are modeled as Rician channels with a Rician factor κ = 4.Moreover, α 0 = 10 −4 is the pass loss of the reference distance d 0 = 1 m with the path loss exponent β = 2.2.Furthermore, λ = 0.05 m and σ 2 U = σ 2 D = −80 dBm. Figure 2 shows the SR performance as a function of the BIOS/IOS size for different designs in a multi-user MISO system with N T = N R = 3 and K D = K U = 2. Specifically, the solid lines illustrate the SR of FD systems, and the dashed lines indicate the rate of HD systems, with blue for BIOS assistance and red for IOS assistance, respectively.Meanwhile, as the BIOS may be too small to enable effective SIC, we switch the system from the FD mode to the HD one with only the DL transmission when the SIC constraint cannot be satisfied, which is indicated by the shaded area.It is worth noting that as the channels are randomly generated, the shaded area means that in most of the simulation realizations, the system has to switch to the HD mode as the size of the BIOS/IOS is too small to achieve the SIC threshold.It can be observed that for all systems, the data rate monotonically increases with M y for both BIOS and IOS due to the enhanced beamforming gains.For FD systems, as M y increases, the SR of BIOS-assisted systems significantly outperforms that of HD systems.However, the IOS-assisted FD system demonstrates only a slight enhancement over HD systems.This is because the IOS beams for SIC and SRM are tightly coupled, limiting the system's efficiency.Specifically, if the IOS is designed to cancel the SI signal on the reflecting side, it cannot be efficiently optimized for SRM on the refracting side.Conversely, the BIOS-assisted FD system can provide flexible beamforming on both sides, thereby enhancing communication performance and achieving SIC simultaneously.To further investigate the beamforming performance of the BIOS and IOS in FD systems, Figure 3 shows the detailed UL and DL rates of the BIOS/IOS-assisted system with different M. As can be observed, the system switches from the FD mode to the HD mode with M y < 9.In this range, the data rate, i.e., the DL rate, consistently increases with M. This is because the expansion of the array scale provides higher degrees of freedom, thus enhancing the beamforming gain.However, when M y changes from 8 to 9, the DL rate unexpectedly decreases while the UL rate substantially increases.This is because the power of the cascaded signal generated by the BIOS/IOS is strong enough to facilitate effective SIC, thus leading to a trade-off between the SRM and SIC for beamforming.It is worth noting that, compared with the IOS-assisted system, the DL rate in the BIOS-assisted system is less affected due to the independent beamforming on both sides.Consequently, the SR performance of the BIOS-assisted FD system significantly surpasses that of the IOS-assisted system.Furthermore, Figure 4 shows the performance comparison between the proposed BIOS and an IOS with the same scattering elements.As the element spacing is kept unchanged, the IOS with 2M elements has a larger aperture size, it can be observed that the IOS with 2M elements slightly outperforms the IOS with M elements due to the additional beamforming gain from the increase in the element number and the aperture size.Similarly, it outperforms the BIOS when M y is small for the same reasons.However, with the increase in M y , the performance of the IOS with 2M elements remains significantly inferior to that of the BIOS due to the coupling between the reflected and refracted beams.Finally, to visualize the SIC capability of the BIOS and verify the theoretical lower bound discussed in Section 4, Figure 5 shows the residual SI power as a function of M y in SISO FD systems with different d BB .Specifically, the solid line represents the residual SI power, and the dashed line of the same color represents the theoretical lower bound of M for efficient SIC in the system with the same d BB .It can be observed that the residual SI power with the proposed scheme monotonically decreases as M increases due to the beamforming gain.Meanwhile, despite the approximation, the theoretical lower bound calculated according to (37) is consistent with the simulation result.

Conclusions
We proposed a novel BIOS-assisted FD system for simultaneous SIC and SRM, and formulated the beamforming optimization problem aiming at maximizing the SR with the SIC requirement constraint.The BS precoder and the BIOS phase shifts were jointly optimized by the proposed SIC-WMMSE approach.Furthermore, a theoretical lower bound of BIOS size required for efficient SIC was derived and proved to be very tight through simulation.The simulation results also showed that the BIOS-assisted FD system significantly outperforms the conventional IOS-assisted system, as the BIOS can flexibly control the beamforming for SIC and SRM.
Proof of Theorem A1.In the DL scenario, considering the independence between UEs, we take the k-th DL UE as an example to prove the equivalence between SRM and WMMSE.The proof is to demonstrate that maximizing the rate of the k-th DL UE is equivalent to solving the following optimization problem: min According to (A4), the optimal a D,k can be rewritten as By taking (A8) back to g D,k , we have Comparing (11) with (A9), we can find that the logarithmic term in (A9) matches the expression of the achievable rate.Consequently, the optimal solution Φ 1 , Φ 2 , F that maximizes the rate in (11) is also the one that minimizes the WMMSE as described in (A9).Similarly, we take the l-th UL UE as an example to prove the equivalence between SRM and WMMSE in the UL scenario.This equivalence is shown by transforming the rate maximization problem of the l-th UL UE into the following minimization problem min Φ 1 ,Φ 2 ,F g U,l = a U,l e U,l − log 2 a U,l . (A10) With Φ 1 , Φ 2 , F fixed, the optimal a U,l , u U,l can be derived as ∇ a U,1 (g U,l ) = 0 ⇒ a U,l = e −1 U,l , (A11)

Figure 1 .
Figure 1.System model of the BIOS-assisted FD system.
respectively, denote diagonal matrices containing the weighting coefficients of each DL and UL UE.U D = diag u D,1 , • • • , u D,K D denotesthe diagonal matrix containing the receive coefficients of each DL UE, and U U = u U,1 , . . ., u U,K U denotes the matrix containing the receive combiners of each UL UE.

3. 5 .
Optimize Φ 2 with Fixed F, Φ 1 , A D , A U , U D , U U Similarly, we optimize Φ 2 with other variables fixed in this subsection.By defining equivalent matrices

Figure 3 .
Figure 3. DL/UL rate performance of BIOS/IOS-assisted systems with different M y .

Figure 5 .
Figure 5. Residual SI power of BIOS-assisted FD systems with different M y .

Φ 1 ,
Φ 2 ,F g D,k = a D,k e D,k − log 2 a D,k .(A3)Itcan be observed that when variables Φ 1 , Φ 2 , and F are fixed, the optimization problem (A3) becomes convex with respect to A D,k , u D,k .Thus, the optimal A D,k can be directly obtained through the first-order conditions∇ a D,k (g D,k ) = 0 ⇒ a D,k = e −1 D,k .(A4)Similarly, the optimal u D,k can be derived as∇ u * D,k (g D,k ) = a D,k ∂ u * D,k u D,k ĥH D,k FF H ĥD,k − 2 Re u * D,k ĥH D,k f k + σ 2 D,k u * D,k u D,k + 1 ∂u * D,k = a D,k u D,k ĥH D,k FF H ĥD,k − ĥH D,k f k + σ 2 D,k u D,k , (A5) ∇ u D,k * (g D,k ) = 0 ⇒ u D,k = ĥH D,k FF H ĥD,k + σ 2 D,k −1 ĥH D,k f k .(A6)By substituting u D,k into e D,k , the minimized MSE can be derived as follows: