\section{Introduction}
	\vspace{-0.1cm}

The increasingly demanding objectives for 6th Generation (6G) communications have spurred recent research activities on novel transceiver hardware architectures and relevant communication algorithms. Such hardware architectures comprise large numbers of ElectroMagnetic (EM) radiating elements, thus paving the way for \ac{mmimo} communications.
A \ac{mmimo} system with arbitrarily large number of antenna elements can provide substantial gains in spectral efficiency with relatively simple signal processing algorithms. This potential has motivated the incorporation of \ac{mmimo} technology in the 5G New Radio (NR) interface
, and \ac{mmimo} transceivers with extremely large number of antennas are considered to continue being one of the key technologies for 6G communications \cite{saad2019vision}.

With the widespread deployment of Internet of Things (IoT) devices,
the number of  nodes connected over wireless media is expected to reach the order of tens of billions in the next few years. To address these massive connectivity, high peak device rates and increased  throughput requirements, future wireless networks are expected to transit into dense deployments of coordinating extreme \ac{mmimo} transceivers (namely, \acp{bs} and access points), particularly in urban and indoor environments. While the theoretical gains of densely deployed \ac{mmimo} systems are still being unveiled \cite{bjornson2019massive},
implementing such systems in practice is a challenging task. In particular, realizing \acp{bs} with hundreds, or even thousands, of antenna elements being able to simultaneously serve multiple users, gives rise to a multitude of practical difficulties for conventional  sub-6 GHz bands as well  millimeter wave and THz bands. Among those challenges are  high fabrication cost, increased power consumption, constrained physical size and shape, and deployment limitations.

Over the last few years, metamaterials have emerged as a powerful technology with a broad range of applications, including wireless communications. Metamaterials comprise a class of artificial materials whose physical properties, and particularly their permittivity and permeability, can be engineered to exhibit various desired characteristics  \cite{smith2004metamaterials}. When deployed in planar structures (a.k.a. metasurfaces), their effective parameters can be tailored to realize a desired transformation on the transmitted, received, or impinging EM waves \cite{DSmith-2017PRA}. Such structures have been lately  envisioned  as a revolutionary means to transform any naturally passive wireless communication environment (the set of objects between a transmitter and a receiver constitute the wireless environment) to an active one \cite{huang2019reconfigurable,di2019smart}.
Their extremely small hardware footprint enables their cost-effective embedding in various 3D components of the environment (e.g., building facades and room walls/ceilings).

\acp{dma} have been recently proposed as an efficient realization of massive antenna arrays for wireless communications \cite{Yoo2018TCOM, shlezinger2019dynamic}. They provide beam tailoring capabilities \textcolor{NewColor}{and facilitate processing of the transmitted and received signals in the analog domain in a flexible and dynamically configurable manner using simplified transceiver hardware. In addition, \ac{dma}-based architectures require much less power and cost compared with conventional antenna arrays (i.e., those based on patch arrays and phase shifters) eliminating the need for complicated corporate feed and / or active phase shifters.} \acp{dma} may comprise a large number of tunable metamaterial antenna elements that can be packed in small physical areas \cite{Akyildiz-2016NCN} for a wide range of operating frequencies. This feature makes them an appealing technology for the extreme \ac{mmimo} transceivers of 6G wireless networks. In contrast to passive metasurfaces that have received extensive attention recently \cite{huang2019reconfigurable,di2019smart}, the potential and capabilities of metasurfaces as active \ac{mmimo} antenna arrays, as well as their associated challenges and future research directions, have not yet been properly treated in the wireless communication literature.

In this article, we discuss the promising application of \acp{dma} as \ac{mmimo} transceivers for future
wireless communications. We commence with a brief introduction on metasurfaces considered  for wireless communications, highlighting the differences between their nearly passive and dynamic counterparts. We then present the unique properties of \acp{dma} as transceiving structures with Analog and Digital (A/D) processing capabilities. We  elaborate on the relationship between \acp{dma} and conventional hybrid A/D \ac{mmimo} architectures, which are based on phase shifters for analog processing \cite{ioushua2019family} and discuss
the main advantages and drawbacks of using \ac{dma}-based \acp{bs} for \ac{mmimo} communications. A simulation study evaluating the ability of \acp{dma} to reliably communicate \textcolor{NewColor}{with a reduced number of Radio Frequency (RF) chains by exploiting their inherent analog signal processing flexibility} is detailed. Next, we present experimental results demonstrating the strong potential benefits of \acp{dma} for wireless communication. We conclude the article with a description of open problems in this area of research and a discussion of future directions towards unveiling the potential of \acp{dma} for 6G wireless communications.

\section{Metasurfaces for Wireless Communications}
\label{sec:Metasurfaces}
\textcolor{NewColor}{
Metamaterials comprise a class of artificial materials in which macroscopic, structured elements mimic the atoms or molecules. Each metamaterial element behaves as an electric or magnetic polarizable dipole, the collection of which can often be characterized by an effective permittivity and permeability. These dipole moments can be engineered so as to achieve desired EM properties} \cite{smith2004metamaterials}. The underlying idea behind metamaterials is to introduce tailored inclusions in a host medium to emulate diverse targeted responses.
\textcolor{NewColor}{
When metamaterial elements are distributed over a planar surface, they are often referred to as metasurfaces} \cite{HMIMOS}. Such 2-Dimensional (2D) structures support the individual tuning of each metamaterial element, allowing the metasurface to carry out different functionalities, such as radiation, reflection, beamforming, and reception of propagating waves \cite{DSmith-2017PRA}. \textcolor{NewColor}{The ability to stack a large number of elements in a limited surface area allows metasurfaces to achieve highly directed signalling enabling holographic beamforming \cite{bjornson2019massive}.}
The properties of each element can often be externally controlled, yielding a dynamic metasurface, also referred to as a reconfigurable intelligent surface \cite{huang2019reconfigurable}.

Two main applications of dynamically tuned metasurfaces have been considered recently in the context of wireless communications, as illustrated in Fig.~\ref{fig:PassiveVsActive1}. These are passive reflective surfaces \cite{ huang2019reconfigurable} and active transceiver antenna arrays \cite{shlezinger2019dynamic}. The former application builds upon the capability of metasurfaces to generate reconfigurable reflection patterns. In this case, metasurfaces deployed in urban settings can facilitate and improve  communication between the \ac{bs} and multiple users by effectively modifying EM signal propagation. \textcolor{NewColor}{The addition of a metasurface enables the communications system as a whole to overcome} harsh non-line-of-sight conditions and improve coverage, when the metasurface is placed within small distances from the BS or the users, without increasing transmission power.
	\begin{figure*}
	\centering
	\includegraphics[width = \linewidth]{PassiveVsActive3.png}
	\vspace{-0.6cm}
	\caption{Two applications of reconfigurable metasurfaces in the downlink (gray arrows) and uplink (blue arrows) of multi-user mMIMO wireless communications: a) as nearly passive reflective surfaces; and b) as active transceiver antenna arrays.
	}
	\label{fig:PassiveVsActive1}
\end{figure*}
Nearly passive reconfigurable intelligent surfaces require some level of control to alter the impinging EM wave in light of the dynamic wireless environment. This is achieved by embedding a digital control unit, which is capable of tuning the metamaterial elements to obtain desired reflection patterns. Such reflective surfaces do not implement conventional relaying techniques (i.e., neither power amplification nor baseband processing) \cite{huang2019reconfigurable}, but only reflect the impinging signal in a controllable manner.

\acp{dma} comprise an additional recent usage of metasurfaces for communication as mMIMO antenna structures. This application exploits their ability to realize planar, compact, and low cost dynamically tunable massive antenna arrays \cite{shlezinger2019dynamic}, which can be deployed in current and future BSs and access points. As such, a DMA-based BS will consist of a multitude of radiating metamaterial elements that can transmit and receive communication signals over the wireless channel. By dynamically tuning  the EM properties of the DMA elements, one can control the analog beampattern for transmission and reception. Unlike passive reflective metasurfaces, which are  the focus of considerable research attention lately, the usage of metasurfaces as active mMIMO antennas is a relatively new area of research. To understand the potential of \acp{dma} for mMIMO systems, we discuss in the following section the architecture, properties, advantages, and drawbacks of using \acp{dma} for transmitting and receiving communication signals.

\section{\acp{dma} for Massive MIMO Communications}
\label{sec:DMAs}
In this section, we discuss the application of metamaterial-based planar antenna arrays in wireless communication systems. We first detail the architecture of DMAs and their main characteristics in the context of communications, and then, elaborate on its operation when deployed for mMIMO \acp{bs}. Finally, some representative numerical evaluations are presented together with experimental results.

 \subsection{\ac{dma} Hardware Architecture}
 \label{subsec:DMAsModel}
 \acp{dma} consist of a multitude of reconfigurable metamaterial radiating elements that can be used both as transmit and receive antennas. Those elements are placed on a waveguide through which the signals to be transmitted, and the received waveforms intended for information decoding, are transferred. The transceiver digital processor, which generates the outgoing signals and processes the received signals, is connected to the waveguide through dedicated input and output ports, respectively. In general, \acp{dma} can use 2D waveguides connected to several input/output ports \cite{Yoo2018TCOM}. However, the common \ac{dma} architecture, on which we focus in this article, consists of multiple separate waveguide-fed element arrays, referred to as {\em 1-Dimensional (1D) waveguides}, each connected to a single input/output port, as illustrated in Fig. \ref{fig:DMA_Diagram1}.
 %

Such waveguides can accommodate a large number of radiating elements, which are commonly sub-wavelength spaced, allowing each input/output port to feed a multitude of possibly coupled radiators. Since this waveguide is typically designed to be single mode and the wave can only propagate along one line, its analysis is much easier than 2D waveguides (such as parallel plate waveguides), where a scattered wave from each element propagates in all directions. Furthermore, ensuring isolation between different ports is easier in 1D waveguides than in multiple ports of a 2D waveguide. A common implementation of 1D waveguides is based on \textit{microstrips}, as illustrated in Fig.~\ref{fig:DMA_Diagram1}.

	\begin{figure*}
	\centering
	\includegraphics[width = \linewidth]{DMA_Diagram3.png}
	\vspace{-0.6cm}
	\caption{An $N$-element \ac{dma} consisting of $M$ micropstrips, where each microstrip is implemented as an 1D waveguide. The upper right part and the lower right part of the figure illustrate their operation during transmission  and reception, respectively, along with their equivalent signal path models.  }
	\label{fig:DMA_Diagram1}
\end{figure*}

When \acp{dma} are deployed as receive antennas, the signals captured at each metamaterial element propagate through the corresponding waveguide to the output port, where they are acquired and forwarded to the digital unit for baseband processing.
In a \ac{dma}-based transmit antenna array, the signals to be radiated from its metamaterial elements are fed to each waveguide via its input port.	The relationships among the radiating signals and those captured/fed at the input/output port of each waveguide are determined by the following two properties arising from the \ac{dma} architecture:
	\begin{enumerate}
		\item \label{itm:P1} Each metamaterial element acts as a resonant electrical circuit whose frequency response is described by a Lorentzian equation \cite{DSmith-2017PRA}. The parameters of the resonant circuit, i.e., its oscillator strength, damping factor, and resonance frequency, are externally controllable and can be configured in run-time for each element individually. An illustration of the normalized magnitude of the element response achieved for several different resonant frequencies is depicted in Fig.~\ref{fig:Lorentz1Mag}. This figure demonstrates that the elements can be tuned to exhibit diverse responses, varying from bandpass to frequency flat filters.

		\begin{figure}
	    \centering
		\scalebox{0.58}{\includegraphics{LorentzMag3.pdf}}

		\caption{The normalized element response magnitudes for different settings of the Lorentzian resonant frequency (Res. freq) as functions of the operating frequency.}
		\label{fig:Lorentz1Mag}
	\end{figure}

		\item \label{itm:P2} In each waveguide, the signal has to travel between the feed port and each metamaterial element. Consequently, the signals propagating along the waveguide accumulate different phases for each metamaterial element.
	\end{enumerate}
	An illustration of the resulting equivalent model for each DMA waveguide, which arises from the aforementioned properties, is depicted in the  right part of Fig.~\ref{fig:DMA_Diagram1}. As shown in this figure, the input/output port is connected to a set of radiating elements via a reconfigurable filter encapsulating the joint effect of the propagation inside the waveguide and the externally controllable frequency response of the metamaterial element.

 \subsection{\ac{dma}-Based mMIMO Transceivers}
 \label{subsec:DMAsBSs}
Similar to conventional massive antenna arrays, \acp{dma} interface the electronic circuitry of a \ac{mmimo} transceiver with the EM waveforms propagating over the wireless channel. However, their application in \ac{mmimo} BSs induces several unique characteristics, which follow from the \ac{dma} structure discussed in the previous section. First, the number of independent data streams that can be processed by a \ac{dma}-based BS in the digital domain is typically much smaller than the number of its metamaterial antenna elements. In particular, a \ac{dma} consisting of $M$ waveguides (e.g., the 1D waveguides in Fig.~\ref{fig:DMA_Diagram1} termed as microstrips), each embedded with $L$  metamaterial antenna elements, allows the transceiver to process $M$ data streams, while utilizing in total $N=ML$ metamaterial elements. Consequently, \ac{dma}-based transceivers implement a form of hybrid A/D beamforming, since part of the processing of the transmitted and received signals is carried out in the analog domain, as an inherent byproduct of the waveguide-fed metamaterial array architecture. Such hybrid processing, which requires additional dedicated hardware in conventional \ac{mmimo} architectures, allows the \ac{bs} to utilize an amount of metamaterial elements that is much larger than the number of the digitally processed data streams. This inherent expansion upon transmission and compression in reception is typically desirable in \ac{mmimo} \acp{bs} as a method to reduce the number of costly RF chains \cite{ioushua2019family}, while achieving increased beamforming gain, as well as to facilitate efficient operation under low quantization requirements \cite{shlezinger2018hardware}.

In DMA-based transceivers, each metamaterial antenna element can exhibit a broad range of frequency responses with various combinations of amplitude and phase values, ranging from frequency selective to frequency flat profiles. These profiles are externally controllable for each element separately. This implies that \ac{dma}-based mMIMO \acp{bs} implement dynamically reconfigurable and frequency selective hybrid A/D beamforming by tuning their  elements to impart different levels of attenuation and phase shift on the transmitted and received signals. Consequently, \ac{dma}-based BSs can be treated as hybrid A/D beamforming systems which do not require additional dedicated analog combining circuitry, while offering increased flexibility, compared to conventional hybrid architectures consisting of interconnections of frequency static phase shifters and switches \cite{ioushua2019family}.

Furthermore, DMAs provide similar beamforming capabilities to those achievable with typical phased array antennas, but with much lower power consumption and cost. To better illustrate this point, we note that the radiation pattern formed by a DMA is the superposition of the radiated field from many metamaterial radiators, whose complex amplitude are determined by two factors: the tunable resonance response of the elements and the phase accumulated by the guided wave. Using simple holographic techniques, the tuning states of the metamaterial elements can be designed to form beams in any direction of interest. Since the tuning of metamaterial elements is usually accomplished with simple components, such as varactors, it requires minimal additional power for beam steering. This is in contrast to many common antenna arrays which use active phase shifters that consume significant power. The beamforming capability of DMAs is heavily dependent on the subwavelength spacing of metamaterial elements. This is due to the fact that the metamaterial elements amplitude and phase cannot be tuned separately. To augment this limited design space, one can utilize the phase accumulation of the feed wave as it travels between different elements. A densely sampled waveguide thus provides enough degrees of freedom to form any desired beamforming pattern.

 \subsection{Numerical and Experimental Results}
 \label{subsec:DMA_Prototype}
 The \ac{dma} architecture facilitates the implementation of massive amounts of metamaterial elements in simple, energy-efficient, compact, and low-profile design configurations, irrespective of the operating frequency. The planar physical shape of \acp{dma} makes them suitable for installations in urban and indoor environments, and the inherent expansion/compression induced by their waveguide-fed architecture reduces the number of required expensive RF chain circuits.
However, this expansion/compression, which stems from the fact that the transceiver can access the signals only at each waveguide's input/output ports, reduces the  achievable performance compared to a fully-digital transceiver which is capable of simultaneously accessing the signals transmitted and received from each antenna element separately. In particular, this performance reduction follows from the DMA requirement for fewer digital streams than antenna elements, and is a common characteristic of hybrid A/D beamforming architectures.

To evaluate the effect of the aforementioned properties on the performance of \ac{dma}-based \ac{mmimo} \acp{bs}, we depict in Fig.~\ref{fig:SEvsSNR_IFreq} the achievable uplink sum-rate performance of narrowband communications at $3.5$GHz carrier frequency as a function of the operating \ac{snr}. The achievable sum-rate is computed by treating the uplink wireless channel as a multiple access channel, and thus the number of RF chains is not restricted to be larger than the number of users.
In this figure, the same performance with a fully-connected hybrid A/D architecture is included, as well as the sum-capacity of the uplink channel when the number of \ac{bs} antennas equals the number $N$ of receiving elements and when it equals the number $M$ of receive RF chains. The corresponding simulation scenario, based on the setup detailed in \cite[Sec. IV]{shlezinger2019dynamic}, consists of $64$ user terminals uniformly distributed in a single cell of $400$m radius that simultaneously communicate in the uplink direction with a mMIMO \ac{bs} having $N=160$ antenna elements. For the DMA architecture, $L=10$ unit metamaterial elements are coated upon each of the $M=16$ waveguides (microstrips in this case). With the fully connected hybrid A/D beamforming architecture, all antenna elements are attached via phase-shifter networks to each of the $M=16$ RF chains.
As observed in Fig.~\ref{fig:SEvsSNR_IFreq}, the DMA architecture yields sum-rate performance closer to the sum-capacity with $M=N$ than conventional hybrid A/D architectures based on fully-connected networks of phase shifters. This advantage follows from the previously discussed capability of \acp{dma} to implement different forms of analog signal processing. The  rate gains of \acp{dma} with $M$ RF chains and $N=ML$ elements over the sum-capacity when the \ac{bs} has $M$ antennas indicates that, connecting each RF chain to a microstrip with multiple elements can notably improve the communication rate compared to using each RF chain to feed a single antenna. Similar results have been reported for the downlink case
as well as
for communications with low-resolution analog-to-digital converters \cite{wang2020dynamic}.

 	\begin{figure}
		\centering
		\scalebox{0.48}{\includegraphics{SEvsSNR_IFreq2_16_64_2.pdf}}

		\caption{Achievable uplink sum-rate performance in bps/Hz versus the SNR in dB for a $160$-antenna mMIMO BS with $16$ RF chains serving simultaneously $64$ users within a cell with $400$m radius \textcolor{NewColor}{based on \cite{shlezinger2019dynamic}}. For the DMA architecture, each of the $M=16$ microstrips includes $L=10$ unit metamaterial elements. In the fully connected hybrid A/D beamforming architecture, the phase-shifter network interconnects each antenna element to all RF chains. }
		\label{fig:SEvsSNR_IFreq}
	\end{figure}

Despite the growing interest in using reconfigurable metasurfaces (either as an intelligent reflector or a transceiver) to augment wireless communication systems, very few experimental studies have been conducted. Most experimental studies focus on verifying the possibility to model each tunable metamaterial radiator as a polarizable dipole and, optimize their tuning states to form desired beams. For instance, the \ac{dma} configuration illustrated in Fig.~\ref{fig:DMA_Diagram1} consists of an array of microstrip, each comprised of a 1D waveguide-fed metasurface. An experimental implementation of such a dynamic 1D waveguide-fed metasurface, based on \cite{sleasman2017experimental,sleasman2016design}, is shown in Fig.~\ref{fig:DMA_sample}. This waveguide is fashioned with metamaterial elements with two different (interleaved) resonant frequencies. Each element is loaded with two PIN diodes, which render the metamaterial element radiating and non-radiating, as shown in Fig.~\ref{fig:DMA_sample}. A close-up view of a fabricated metamaterial is also shown. In addition, the figure includes experimental results for the ability of the metasurfaces to form beams in different directions, showcasing the potential of DMA-based mMIMO transceivers in forming multiple directed beams.

 	\begin{figure}
		\centering
		\scalebox{0.5}{\includegraphics{Comm_magazine.png}}

		\caption{\textcolor{NewColor}{An experimental implementation of a dynamic 1D waveguide-fed metasurface \cite{sleasman2017experimental,sleasman2016design}:} a) Detailed circuitry of a reconfigurable metamaterial element. b) The simulated response of the metamaterial element, where the impact of the PIN diodes to render the element in radiating ($e_{1,\rm on}$ and $e_{2,\rm on}$) and non-radiating ($e_{1,\rm off}$ and $e_{2,\rm off}$) states are shown. c) A close-up view of a sample 1D DMA with metamaterial elements having two different resonance frequencies. d) Beamforming capability of the fabricated 1D DMA.}
		\label{fig:DMA_sample}
	\end{figure}

\section{Open Research Challenges}
\label{sec:Challenges}
As previously discussed, the hardware architecture of \acp{dma} enables  efficient implementation of extreme mMIMO transceivers, which renders this concept a strong candidate technology for the physical layer of 6G communications. However, \acp{dma} come with certain design challenges that need to be carefully addressed, while giving rise to new opportunities. In the following, we list some of the most important related research challenges to date.

\textbf{Frequency-Selective Analog Beamforming:} Current algorithmic designs for \acp{dma} focus either on narrowband communications or ignore their capability to dynamically configure the frequency-selective profile of each unit metamaterial element. This unique property, which does not exist in any conventional hybrid A/D architecture, provides increased flexibility for the design of wideband operation by matching the spectral behavior of each element to optimize the equivalent wideband channel. Consequently, the true potential of \ac{dma}-based mMIMO systems in achieving reliable and ultra high rate communications is not yet fully explored.

\textbf{Wireless Channel Estimation and Tracking:} To date, studies on \ac{dma}-based communications assume that the transceiver has full channel knowledge. In practice, however, the channel coefficients need to be efficiently estimated, which is a challenging task with hybrid A/D receivers. When channel estimation is carried out in a time-division duplexing manner, \acp{dma} offer the possibility of tuning their elements to facilitate channel estimation via pilot signals, and to adapt in a manner which optimizes data reception in light of the estimated channel. The design and analysis of efficient algorithms for \ac{dma}-based \acp{bs}, which have to estimate features of the wireless channel and reliably communicate have not yet been properly treated in the literature.

\textbf{Hardware Design and Experimentation:} A large body of fabricated designs and experimental work is still required in order to transit the DMA concept into an established technology for future extreme mMIMO \acp{bs} with hundreds or thousands of metamaterial elements. The experimental studies detailed in the previous section  are an important first step which demonstrates the feasibility of this concept. Yet, the performance and other implications of using \ac{dma}-based \acp{bs} need be tested in a broad range of wireless setups. To date, most experimental studies of \acp{dma} are geared toward imaging and radar systems.

In particular, several key properties of \acp{dma} should be thoroughly tested to fully understand their potential for mMIMO communications. First, one  needs to quantify the level of correlation among the signals received by different metamaterial elements on the same waveguide as well as different waveguides. The impact of this correlation on the overall response needs to be investigated. In the example discussed in the previous section, the tuning of the metamaterial elements is optimized based on certain communication criteria. However, the resulting pattern may introduce more correlation or result in low \ac{snr} (for example, in a scenario where the resulting pattern is directed toward directions other than the intended users). The co-design of the metamaterial elements tuning states with the resulting pattern is another important factor that needs to be examined. When used as a transmitter, the possibility of the high amplitude RF signal (the carrier) to push the tunable component of each metamaterial element into a nonlinear regime needs to be investigated. This nonlinearity can cause intersymbol or intercarrier interference.
Finally, although conceptually feasible, \acp{dma} for communications in the high millimeter wave (around 60 GHz) and THz (above 100 GHz) bands are still in relatively early stages of research.

\textbf{Hybrid Passive and Active Metasurfaces:}
As discussed in the introduction, the emerging concept of reconfigurable metasurfaces, with both its passive and active counterparts, is lately gaining increased interest for both controllable reflection and transmission/reception. It is reasonable to envision hybrid passive and active reconfigurable metasurfaces, notably strengthening the design flexibility for such surfaces and their performance gains in programmable wireless environments. For instance, having such a hybrid metasurface acting as a receiving device can significantly facilitate channel estimation \cite{GA2020_RIS}, which is currently a major challenge and a source of substantial overhead in purely passive metasurfaces. In addition, hybrid metasurfaces will enable more advanced relaying strategies, overcoming the dominating impact of pathloss in their passive versions.

\textbf{Use Cases and Applications:} The use cases and applications where \acp{dma} can provide substantial improvement compared to current architectures have not yet been thoroughly identified. For example, their planar shape and reduced size for a small number of RF chains facilitate their deployment in indoor environments, like buildings, malls, train stations, and airports. In such setups,  \acp{dma} are expected to communicate with multiple users in close line-of-sight conditions, possibly operating in the near-field regime. Such near-field scenarios bring forth the possibility of beam focusing,
namely, the ability to focus the signal towards a specific location, instead of a specific direction as in far-field conditions via conventional beamforming, allowing multiple orthogonal links to exist in the same indoor or outdoor area.

\section{Conclusion}
\label{sec:Conclusion}
\acp{dma} are an attractive radiating technology for next generation wireless systems, allowing to realize dynamically controllable \ac{mmimo} antennas of reduced cost and power consumption compared to conventional arrays. In this article we surveyed some of the key properties of \acp{dma} in the context of \ac{mmimo} systems, including their operation model during their transmission and reception, tunable frequency selective profile, beam steering capabilities, as well as advantages and drawbacks compared to conventional arrays. We concluded with a set of key open research directions, which are expected to pave the way in unveiling the full potential of using active metasurfaces in 6G wireless communications.

	\bibliographystyle{IEEEtran}


\section{Main Results}}

The irrational rotation C*-algebra $A_\theta$ is generated by unitary operators $U,V$ satisfying the Heisenberg commutation
\[
VU = e^{2\pi i\theta} UV
\]
where $\theta$ is an irrational number. The Flip $\Phi,$ Fourier transform $\sigma,$ and Cubic transform $\kappa,$ are the canonical automorphisms of $A_\theta$ defined by
\begin{align*}
 \Phi(U) &= U^{-1},  &  \sigma(U) &= V^{-1}, & \kappa(U) &= e(-\tfrac\theta2)U^{-1}V,
\\
\Phi(V) &= V^{-1}, &  \sigma(V) &= U,	 & \kappa(V) &= U^{-1},
\end{align*}
where we shall sometimes use the conventional notation $e(\theta) := e^{2\pi i\theta}$. The Flip orbifold is the fixed point C*-subalgebra of $A_\theta$ given by $A_\theta^\Phi = \{x\in A_\theta: \, \Phi(x) = x\}$. The Flip was studied extensively, e.g. in \ccite{BEEKa} \ccite{BEEKb} \ccite{SW-JLMS}, and it is known that the Flip orbifold is approximately finite dimensional for irrational $\theta$ (\ccite{BK} and \ccite{SW-CMP}).
\medskip

The well-known modular action of the group SL$(2,\mathbb Z)$ on $A_\theta$ (due to Watatani \ccite{YW} and Brenken \ccite{BB}) is given by associating to the integral matrix
$\left[\smallmatrix a & c \\ b & d \endsmallmatrix\right]$ (with $ad-bc=1$) the automorphism $\alpha$ defined by
\begin{equation}\label{modaut}
\alpha(U) = e(\tfrac\theta2 ab)\, U^a V^b, \qquad
\alpha(V) = e(\tfrac\theta2 cd)\, U^c V^d.
\end{equation}
We call these modular automorphisms of the rotation algebra. It is straightforward to check that the unitaries on the right sides of these equalities satisfy the commutation relation $VU = e(\theta) UV,$ guaranteeing that the automorphism $\alpha$ exists by the universal property of this relation.
\medskip

It's easy to check that all these modular automorphisms commute with the Flip, being given by minus the identity matrix $\left[\smallmatrix -1 & 0 \\ 0 & -1 \endsmallmatrix\right]$. Therefore, the modular group SL$(2,\mathbb Z)$ also acts on the Flip orbifold $A_\theta^\Phi,$ leaving it invariant.
\medskip

\begin{dfn}
A projection $e$ in a C*-algebra $A$ is called {\bf approximately central} (or AC, for short) if $e=e_n$ is a sequence depending on an integer parameter $n$ such that
\[
\lim_{n\to \infty} \|x e_n - e_n x\| \ = \ 0
\]
for each $x\in A$.  Two AC projections $e=e_n$ and $f=f_n$ are said to be {\bf centrally equivalent} or {\bf AC-equivalent} in $A$ if $\forall \epsilon > 0$ and each finite subset $F$ of $A,$ there exists a partial isometry $u \in A$ such that $uu^* = e_n, \ u^*u = f_n$ (for large enough $n$) and $\|x u - u x\| < \epsilon$ for each $x\in F$. We use the notation
\[
e  {\text{\ \Large$\sim_{{}_{\!\!\!\!\!\!\!\!c}}$\ \,}}   f
\]
to indicate central equivalence of $e$ and $f$ in $A$.
\end{dfn}

(Of course, this last condition is just usual Murray-von Neumann equivalence with the added condition that the partial isometry, or unitary, is approximately central in the algebra $A$.)
\medskip

\begin{dfn}
A projection $e$ in the rotation algebra is called {\bf smooth} when it belongs to the canonical smooth dense *-subalgebra $A_\theta^\infty$ consisting of rapidly decreasing Laurent series $\rmsumop a_{mn}U^mV^n$. An automorphism is called {\bf smooth} if it maps $A_\theta^\infty$ onto itself.
\end{dfn}
\medskip

Of course, all modular automorphisms are smooth. The smooth Flip orbifold $A_\theta^{\infty,\Phi} = A_\theta^\infty \cap A_\theta^\Phi$ is also a dense *-subalgebra of $A_\theta^\Phi$ and is closed under the holomorphic functional calculus (as is the case for $A_\theta^\infty$ and $A_\theta$).
\medskip

We can now state the main results.

\begin{thm}\label{mainthm}
Let $\alpha$ be a modular automorphism of the irrational rotation C*-algebra $A_\theta$ arising from the canonical action of SL$(2,\mathbb Z)$. For any approximately central projection $e$ in the Flip orbifold $A_\theta^\Phi,$ the projection $\alpha(e)$ is centrally equivalen in $A_\theta^\Phi$ to one of the projections in its $S_3$-orbit
\begin{equation}\label{orbite}
e, \quad \sigma(e), \quad \kappa(e), \quad \kappa^2(e), \quad \sigma\kappa(e), \quad \sigma\kappa^2(e).
\end{equation}
\end{thm}
\bigskip

The symmetric group $S_3$ is here envisaged as generated by $\sigma$ and $\kappa$ subject to the relation $\kappa\sigma = \sigma\kappa^2$ (and $\sigma^2 = \kappa^3 = 1$). Some of the projections in the $S_3$-orbit of $e$ could be centrally equivalent to one another. Indeed, they can all be according to the following result.
\medskip

\begin{thm}\label{prop}
Let $\alpha$ be a smooth automorphism of the irrational rotation C*-algebra $A_\theta$ commuting with the Flip $\Phi$. Then for any AC smooth projection $g$ in $A_\theta$ such that $g \Phi(g) = 0,$ the projections $e = g + \Phi(g)$ and $\alpha(e)$ are centrally equivalent in $A_\theta^\Phi$. In particular, the $S_3$-orbit projections of $e$ \eqref{orbite} are all centrally equivalent in $A_\theta^\Phi$.
\end{thm}
\medskip

(Note that the automorphism in this theorem is not required to be modular.) In a forthcoming paper \ccite{SWb} we construct an AC Powers-Rieffel projection such that its six $S_3$-orbit projections \eqref{orbite} are pairwise not centrally equivalent (Theorem 1.4 of \ccite{SWb}).
\medskip

The orbifold $A_\theta^{\infty,\Phi}$ is known to have four basic unbounded traces denoted
$\phi_{jk}$ (see equation \eqref{traces}). The following result shows that smooth automorphisms act on them in a very specific manner.

\begin{thm}\label{mainB}
Let $\beta$ be an automorphism of the smooth Flip orbifold $A_\theta^{\infty,\Phi}$. One has the linear combinations
\[
\phi_{00}\beta = \phi_{00} + b\phi_{01} +  c\phi_{10} + d \phi_{11}
\]
\[
\phi_{jk}\beta = b_{jk}\phi_{01} +  c_{jk} \phi_{10} + d_{jk} \phi_{11}
\]
for $jk\in\{01,10,11\},$ where $b, c, d, b_{jk}, c_{jk}, d_{jk} \in \tfrac12\mathbb Z$ are half-integers. In particular, the latter gives a group representation $\text{Aut}(A_\theta^{\infty,\Phi}) \to \text{SL}_{\pm}(3,\mathbb R)$.
\end{thm}
\medskip

(The last representation is discussed in Remark 3.6 below.) From this we are able to say when an AC projection is centrally equivalent to a smooth automorphic image.

\begin{thm}\label{mainC}
Let $\beta$ be a smooth automorphism of the Flip orbifold that induces the identity on $K_0(A_\theta^{\infty,\Phi})$. Then
\[
\beta(e)  {\text{\ \Large$\sim_{{}_{\!\!\!\!\!\!\!\!c}}$\ \,}}   e
\]
are centrally equivalent in $A_\theta^{\Phi}$ for any approximately central projection $e$ in $A_\theta^{\infty,\Phi}$.
\end{thm}

\begin{cor}
Let $\alpha, \beta$ be smooth automorphisms of the Flip orbifold $A_\theta^{\infty,\Phi}$ such that $\alpha_* = \beta_*$ on $K_0(A_\theta^{\Phi})$ (e.g., if they are homotopic). Then
\[
\alpha(e)  {\text{\ \Large$\sim_{{}_{\!\!\!\!\!\!\!\!c}}$\ \,}}   \beta(e)
\]
in $A_\theta^{\Phi}$ for any approximately central projection $e$ in $A_\theta^{\infty,\Phi}$.
\end{cor}

\begin{rmk}
The canonical toral automorphism $\gamma_1(U) = -U,\  \gamma_1(V)=V$ does not induce the identity on $K_0(A_\theta^{\Phi})$ of the orbifold (see \eqref{gamma1Q}), but to show that $\gamma_1(e)$ and $e$ are not centrally equivalent for some AC projection $e$ in the Flip orbifold is not trivial. This, however, is shown to follow from \ccite{SWb} for an AC Powers-Rieffel projection $e$.  (The same is also shown to occur for the Fourier and Cubic transforms on the Flip orbifold.)
\end{rmk}

\medskip

We recall the basic unbounded $\Phi$-traces $\phi_{jk}$ on the rotation algebra $A_\theta$ are defined by
\begin{equation}\label{traces}
\phi_{jk}(U^mV^n)\ =\ e(-\tfrac{\theta}2 mn)\,\updelta_2^{m-j} \updelta_2^{n-k}
\end{equation}
for $jk = 00, 01,10,11,$ where $\updelta_2^n$ is the divisor delta function defined to be 1 if $n$ is even, and 0 if $n$ is odd.
\medskip

By definition, a $\Phi$-trace on $A_\theta$ is a linear functional $\phi$ (usually discontinuous) defined on the canonical smooth subalgebra $A_\theta^\infty$ such that
\[
\phi(xy) = \phi(\Phi(y) x)
\]
for all smooth elements $x,y \in A_\theta^\infty$. It is well-known (e.g., from \ccite{SW-JLMS}) that the complex vector space of all $\Phi$-traces on $A_\theta^\infty$ is 4-dimensional and has the four functionals \eqref{traces} as basis.  It is easy to see that $\phi\alpha$ is a $\Phi$-trace for each $\Phi$-trace $\phi$ and each smooth automorphism $\alpha$ commuting with the Flip. So, in particular, $\phi\alpha$ is a linear combination of the $\phi_{jk}$. As $\Phi$-traces on $A_\theta$, the maps $\phi_{jk}$ restrict to actual (unbounded) trace maps on the smooth Flip orbifold, and hence induce real-valued group morphisms on $K_0(A_\theta^{\infty,\Phi}) = K_0(A_\theta^{\Phi})$.
\medskip

The canonical trace $\uptau$ and the unbounded traces $\phi_{jk}$ together form the Connes-Chern character map for the Flip orbifold given by
\begin{equation}
\bold T : K_0(A_\theta^\Phi) \to \mathbb R^5, \qquad
\bold T(x) = (\uptau(x); \phi_{00}(x), \phi_{01}(x), \phi_{10}(x), \phi_{11}(x))
\end{equation}
since each positive $K_0$ class $x$ has a smooth representative. In \ccite{SW-JLMS} (Proposition 3.2) this map was shown to be injective for irrational $\theta$.\footnote{In \cite{SW-JLMS} we worked with the crossed product algebra $A_\theta \times_\Phi \mathbb Z_2$, but since this is strongly Morita equivalent to the fixed point algebra, the injectivity follows for the fixed point subalgebra.}  The ranges of the traces $\phi_{jk}$ on projections are known to be half-integers: $\tfrac12\mathbb Z$, while the canonical trace has range $(\mathbb Z + \mathbb Z\theta) \cap [0,1]$ on projections in $A_\theta^\Phi$.
\medskip

If $m,n$ is any pair of integers, we shall write $\phi_{mn} = \phi_{\conj m\,\conj n}$ where $\conj m = 0, 1$ according to whether $m$ is even or odd (respectively).
\medskip

The zeroth cyclic cohomology (see \ccite{AC}) $HC^0(A_\theta^{\infty,\Phi})$ of the Flip orbifold is a 5-dimensional complex vector space spanned by the (bounded) canonical trace $\tau$ and the four unbounded traces $\phi_{jk}$. Indeed, this can be seen in much the same way as how the zeroth cyclic cohomology was calculated for the Fourier transform in \ccite{SWChern} (Proposition 2.2 and Theorem 2.3). Thus $HC^0(A_\theta^{\infty,\Phi}) \cong \mathbb C^5.$

\begin{cor}\label{coro}
A smooth automorphism of the Flip orbifold that induces the identity map on $K_0(A_\theta^{\Phi}) = \mathbb Z^6$ necessarily induces the identity map on the zeroth cyclic cohomology group $HC^0(A_\theta^{\infty,\Phi}) = \mathbb C^5$.
\end{cor}
\medskip

\noindent{\bf Acknowledgement.}
This paper and \ccite{SWb} were written at about the time the author retires, he is therefore most grateful to his home institution of 26 years, the University of Northern British Columbia, for many years of research support. The author also expresses his nontrivial gratitude to the many referees who made helpful review reports over the years (including critical ones).

{\Large \section{Proof Of Theorem \ref{mainthm}}}

We shall let $\alpha$ be a fixed modular automorphism, as in equation \eqref{modaut}. We require the following lemma.

\medskip

\begin{lem}\label{lemma}
Let $\alpha$ be a modular automorphism of the rotation algebra $A_\theta$ induced by the integral matrix $\left[\smallmatrix a & c \\ b & d \endsmallmatrix\right]$ (where $ad-bc=1$). Then
\[
\phi_{00}\alpha = \phi_{00}
\]
and its action on the remaining unbounded $\Phi$-traces is given as follows
\[
\text{For even } a: \quad \phi_{01}\alpha = \phi_{10}, \quad \phi_{10}\alpha = \phi_{d,1}, \quad \phi_{11}\alpha = \phi_{1-d, 1}
\]
\[
\text{For even } b: \quad \phi_{01}\alpha = \phi_{c,1}, \quad \phi_{10}\alpha = \phi_{10}, \quad \phi_{11}\alpha = \phi_{1-c, 1}
\]
\[
\text{For even } c: \quad \phi_{01}\alpha = \phi_{01}, \quad \phi_{10}\alpha = \phi_{1,b}, \quad \phi_{11}\alpha = \phi_{1,1-b}
\]
\[
\text{For even } d: \quad \phi_{01}\alpha = \phi_{1,a}, \quad \phi_{10}\alpha = \phi_{01}, \quad \phi_{11}\alpha = \phi_{1,1-a}
\]
\end{lem}
\begin{proof} A straightforward computation verifies that
\[
(U^k V^\ell)^m = e(\tfrac{\theta}2 k\ell m(m-1))\ U^{km} V^{\ell m}.
\]
In view of \eqref{modaut}, this gives
\[
\alpha(U^m V^n) = e(\tfrac{\theta}2[abm^2 + cdn^2])\ e(\theta bcmn) \ U^{am+cn} V^{bm+dn}.
\]
From this it follows that $\alpha$ transforms the $\Phi$-trace functionals according to
\[
\phi_{jk}\alpha (U^mV^n)\
=\ e(-\tfrac{\theta}2 mn)\, \updelta_2^{am+cn-j} \updelta_2^{bm+dn-k}.
\]
Since $ad-bc=1$, one of the integers $a,b,c,d$ has to be even, so it will be convenient to express $\phi_{jk}\alpha$ in terms of the basic $\Phi$-traces $\phi_{jk}$ by considering the four cases where one of these integers is even.
\medskip

{\bf Case 1:} $a$ is even. Here, $b,c$ are odd so they can be replaced by 1's in the delta functions, and of course $a$ disappears, so we have
\[
\phi_{jk}\alpha (U^mV^n)\ =\ e(-\tfrac{\theta}2 mn)\, \updelta_2^{m+dn-k} \updelta_2^{n-j}
=\ e(-\tfrac{\theta}2 mn)\, \updelta_2^{m+dj-k} \updelta_2^{n-j}
\]
(since from the second delta, $n$ can be replaced by $j$ in the first delta), so
\[
\phi_{jk}\alpha = \phi_{k-dj, j}		\tag{$a$ even}
\]
(where, as we noted, the subscript index $k-dj$ is reduced mod 2). Explicitly,
\[
\phi_{00}\alpha = \phi_{00}, \quad \phi_{01}\alpha = \phi_{10}, \quad
\phi_{10}\alpha = \phi_{d,1}, \quad \phi_{11}\alpha = \phi_{1-d, 1}.
\]
In similar vain we proceed to deal with the remaining cases.

\medskip

{\bf Case 2:} $b$ is even. Here $a,d$ are odd so they can be replaced by 1's so we have
\[
\phi_{jk}\alpha (U^mV^n)\ =\ e(-\tfrac{\theta}2 mn)\, \updelta_2^{m+ck-j} \updelta_2^{n-k}
\]
from which
\[
\phi_{jk}\alpha = \phi_{j-ck, k}		\tag{$b$ even}
\]
\medskip

{\bf Case 3:} $c$ is even. Here $a,d$ are odd so they can be replaced by 1's so we have
\[
\phi_{jk}\alpha (U^mV^n)\
= e(-\tfrac{\theta}2 mn)\, \updelta_2^{m-j} \updelta_2^{n+ bj-k}
\]
from which
\[
\phi_{jk}\alpha = \phi_{j, k-bj}.		\tag{$c$ even}
\]
\medskip

{\bf Case 4:} $d$ is even. Here $b,c$ are odd, so
\[
\phi_{jk}\alpha (U^mV^n)\
=\ e(-\tfrac{\theta}2 mn)\, \updelta_2^{m-k} \updelta_2^{n + ak-j}
\]
from which
\[
\phi_{jk}\alpha = \phi_{k, j-ak}.		\tag{$d$ even}
\]
From each of these four cases we obtain the results stated in the lemma.
\end{proof}

\medskip

Note that the actions on the three unbounded traces $\phi_{01}, \phi_{10}, \phi_{11}$ give rise to all their cyclic permutations, exhausting the symmetric group $S_3$ on these three objects. For instance, in the first case where $a$ is even, if $d$ is also even we have
\[
\phi_{01}\alpha = \phi_{10}, \quad \phi_{10}\alpha = \phi_{01}, \quad \phi_{11}\alpha = \phi_{11}
\]
so $\alpha$ acts like a transposition on the first two (like the Fourier transform), and when $d$ is odd we have
\[
\phi_{01}\alpha = \phi_{10}, \quad
\phi_{10}\alpha = \phi_{11}, \quad \phi_{11}\alpha = \phi_{01}
\]
so $\alpha$ induces a cyclic permutation (like the inverse Cubic transform $\kappa^{-1}$). One gets the other transpositions and cyclic permutations by similarly examining the other cases for $b,c,d$.
\medskip

This shows that modular automorphisms act on $\phi_{01}, \phi_{10}, \phi_{11}$ simply by  permutations, and that all permutations arise from the symmetric group $S_3$ on three letters generated by the order 2 element $\sigma$ (Fourier transform on the Flip orbifold) and the order 3 element $\kappa$ (Cubic transform) satisfying the commutation $\kappa\sigma = \sigma\kappa^2$.

\medskip

Having fixed the automorphism $\alpha$, by Lemma \ref{lemma} choose $\beta$ in the symmetric group
\[
S_3 = \{ I, \sigma, \kappa, \kappa^2, \sigma\kappa, \sigma\kappa^2\}
\]
such that the $\gamma := \beta^{-1}\alpha$ fixes all the unbounded traces: $\phi_{jk}\gamma = \phi_{jk}$. It follows that $\bold T\gamma(x) = \bold T(x)$ on $K_0$ and hence $\gamma_* = id_*$ is the identity on $K_0(A_\theta^\Phi)$.
\medskip

Now suppose $[P]$ is a class from any given (but fixed) basis for $K_0(A_\theta^\Phi)$ where $P$ is a projection in $A_\theta^\Phi$. The cutdown of $P$ by the projection $\gamma(e)$ is the projection $\chi(\gamma(e) P \gamma(e))$ (defined since $e$ is approximate central), where $\chi$ is the characteristic function of the interval $[\tfrac12,\infty)$.  Writing $\phi = \tau, \phi_{jk}$ for short, the topological invariant of this cutdown is
\[
\phi \chi(\gamma(e)P\gamma(e)) = \phi\gamma \chi(e \gamma^{-1}(P) e)
= \phi \chi(e \gamma^{-1}(P) e)
\]
since $\phi\gamma = \phi$. As $[\gamma^{-1}(P)] = [P]$ in $K_0(A_\theta^\Phi),$ they are unitarily equivalent: $\gamma^{-1}(P) = uPu^*$ for some Flip invariant unitary $u$ (by \ccite{SW-JLMS}, Theorem 5.3 or Corollary 5.6). Hence,
\begin{equation*}
\phi \chi(\gamma(e) P \gamma(e))
= \phi \chi(e uPu^* e)
= \phi( u\chi(eP e) u^* )
= \phi( \chi(eP e) )
\end{equation*}
\vskip5pt

\noindent using the trace property of $\phi$ (and $u$ being Flip invariant). This shows that the cutdown projections $\chi(\gamma(e) P \gamma(e))$ and $\chi(eP e)$ have the same canonical trace and same unbounded traces so that their Connes-Chern characters $\bold T$ are equal. Therefore, since $\bold T$ is injective on $K_0$, their classes in $K_0(A_\theta^\Phi)$ are equal.  By Kishimoto's\footnote{See Remark 2.9 in \ccite{AK} where it is noted that his Theorem 2.1 applies equally to simple AT-algebras of real rank zero, which include the irrational rotation C*-algebra, and its canonical orbifolds under the canonical automorphisms of order 2, 3, 4, and 6. Further, the $K_1$ part of his theorem trivially holds since the $K_1$ group of the Flip orbifold is 0, being an AF-algebra.}Theorem 2.1 \ccite{AK}, the projections $\gamma(e)$ and $e$ are centrally equivalent $\gamma(e)  {\text{\ \Large$\sim_{{}_{\!\!\!\!\!\!\!\!c}}$\ \,}}   e$, so that $\gamma(e) = \beta^{-1}\alpha(e) = vev^*$ for some approximately central unitary $v$ in the Flip orbifold. Applying $\beta,$ we see that $\alpha(e)  {\text{\ \Large$\sim_{{}_{\!\!\!\!\!\!\!\!c}}$\ \,}}   \beta(e)$ are centrally equivalent in $A_\theta^\Phi$. Since $\beta$ is one of the six basic automorphisms mentioned earlier, it follows that $\alpha(e)$ is centrally equivalent to one of the projections in the $S_3$-orbit of $e$:
\begin{equation}
e, \quad \sigma(e), \quad \kappa(e), \quad \kappa^2(e), \quad \sigma\kappa(e), \quad \sigma\kappa^2(e).
\end{equation}
In particular, for any modular automorphism $\alpha$, the AC projection $\alpha(e)$ is centrally equivalent to one of these projections in the Flip orbifold. This proves Theorem \ref{mainthm}.
\medskip

We now prove Theorem \ref{prop}.

\begin{proof} (Of Theorem \ref{prop}.) One can construct AC projections $g$ in the irrational rotation C*-algebra $A_\theta$ orthogonal to their Flip: $g\Phi(g) = 0$ (see \ccite{SWmathscand}). The resulting Flip-invariant projection $e = g + \Phi(g)$ is also approximately central in $A_\theta$. Now fix a class $[P]$ from any given (but fixed) basis for $K_0(A_\theta^\Phi)$ where $P$ is a (smooth) projection in $A_\theta^\Phi$. Then since $\alpha^{-1}(P)$ and $P$ have the same trace, they are unitarily equivalent in the rotation algebra $A_\theta$ (see \ccite{MR}, Corollary 2.5, based on the highly nontrivial results of \ccite{PV} and \ccite{MR1981}). So let $w \in A_\theta$ be a unitary such that $\alpha^{-1}(P) = wPw^*$.
Since $\tau\alpha = \tau,$ we have
\[
\tau \chi(\alpha(e) P \alpha(e) ) = \tau\alpha \chi( e \alpha^{-1}(P) e) =
\tau \chi( e wPw^* e) = \tau (w\chi( e P e)w^*) = \tau (\chi( e P e))
\]
since the projections $\chi( e wPw^* e)$ and $w\chi( e P e)w^*$ are close (as $e$ is AC) so therefore have the same trace.

Next, let $Q = \alpha^{-1}(P)$ and note that from the norm approximation $gQ \Phi(g) \approx 0$ (where $a \approx b$ is short for $\|a-b\|\to 0$), we have
\[
e Q e = (g + \Phi(g)) Q (g + \Phi(g)) \approx gQg + \Phi(g)Q \Phi(g) = gQg + \Phi(gQg)
\]
is a sum of two orthogonal positive elements, each norm-close to their squares (since $g$ is AC), thus giving
\[
\chi(e Q e) \approx \chi(gQg) + \Phi\chi(gQg)
\]
where both sides are $\Phi$-invariant projections. On the right side one has $\phi \chi(g Q g) =  0$ for each $\Phi$-trace $\phi$ (and $\phi \Phi\chi(gQg) = 0$). To see this, we observe that if $f = h + \Phi(h),$ where $h \Phi(h) = 0,$ then $\phi(f) = 0$ for each $\phi$ since $\phi(h) = \phi(hh) = \phi(\Phi(h)h) = 0$. Therefore, $\phi\chi(e Q e) = 0$ and we have
\[
\phi_{jk}\chi(\alpha(e) P \alpha(e)) = \phi_{jk} \alpha \chi(e Q e) = 0
\]
since $\phi_{jk} \alpha$ is also a $\Phi$-trace ($\alpha$ being $\Phi$-invariant). This shows that the projections $\chi(\alpha(e) P \alpha(e))$ and $\chi(e P e)$ have the same canonical trace and both have vanishing unbounded traces, so they yield the same class in $K_0(A_\theta^\Phi),$ and thus, by Kishimoto's Theorem, $\alpha(e)$ and $e$ are centrally equivalent in $A_\theta^\Phi$ (bearing in mind the comment in footnote 2 regarding $K_1=0$).
\end{proof}

{\Large \section{Smooth Orbifold Automorphisms}}

\medskip

In this section we show how a general smooth automorphism of the Flip orbifold transforms unbounded traces (Theorem \ref{mainB}). Thereby we prove that when it induces the identity on $K_0$ of the Flip orbifold it transforms AC projections to ones that are centrally equivalent to them (Theorem \ref{mainC}).
\medskip

In the proof below we shall make use of the toral automorphisms $\gamma_1, \gamma_2, \gamma_3$ of the rotation algebra, which are also automorphisms of the Flip orbifold (since they commute with the Flip), given by
\begin{align}\label{gammas}
\gamma_1(U) &= -U, \quad \gamma_1(V)=V	\notag
\\
\gamma_2(U) &= U, \quad \ \ \ \gamma_2(V)=-V
\\
\gamma_3(U) &= -U, \quad \gamma_3(V)=-V		\notag
\end{align}
where $\gamma_3=\gamma_1\gamma_2$. (It's easy to check that these, aside from the identity, are the only toral automorphisms, arising from the usual action of the torus $\mathbb T^2,$ that commute with the Flip.) They transform the unbounded traces according to
\begin{equation}\label{gammaphis}
\phi_{jk} \gamma_1 = (-1)^j \phi_{jk}, \qquad
\phi_{jk} \gamma_2 = (-1)^k \phi_{jk}, \qquad
\phi_{jk} \gamma_3 = (-1)^{j+k} \phi_{jk}.
\end{equation}

\bigskip

\begin{thm}\label{betaspan}
Let $\beta$ be an automorphism of the smooth Flip orbifold $A_\theta^{\infty,\Phi}$. One has the linear combination
\[
\phi_{jk}\beta = a\phi_{00} + b\phi_{01} +  c\phi_{10} + d \phi_{11}
\]
where $b,c,d$ are half-integers. Further, $a = 1$ when $jk=00,$ and $a=0$ otherwise.
\end{thm}
\begin{proof} Fix $j,k$. From the paragraph just before Corollary \ref{coro}, since $\phi_{jk}\beta$ is a trace map on $A_\theta^{\infty,\Phi},$ it is a linear combination
\begin{equation}\label{phibetacomb}
\phi_{jk}\beta = K \tau + a \phi_{00} + b\phi_{01} +  c\phi_{10} + d \phi_{11}
\end{equation}
where $K,a,b,c,d$ are complex constants (depending on $j,k$).

By Lemma 3.2 in \ccite{SamHouston2018}, there are Flip-invariant Powers-Rieffel projections $Q$ with Connes-Chern character $\bold T(Q) = (\theta'; \tfrac12, \tfrac12, \tfrac12, \tfrac12)$ for a dense set of irrationals $\theta' \in \mathbb (Z + \mathbb Z\theta) \cap (0,\tfrac12)$. Pick such a $Q$ with $\theta'$ in the range $\tfrac13 < \theta' < \tfrac49 < \tfrac12,$ so that $0 < 3(3\theta'-1) < 1$.

Applying the toral automorphisms $\gamma_j$ to $Q$ gives the following Connes-Chern characters (in view of \eqref{gammaphis})
\begin{align}\label{gamma1Q}
\bold T(Q) & = (\theta'; \tfrac12, \tfrac12, \tfrac12, \tfrac12)		\notag
\\
\bold T(\gamma_1Q) &= (\theta'; \tfrac12,  \tfrac12, -\tfrac12, -\tfrac12)
\\
\bold T(\gamma_2 Q) &= (\theta'; \tfrac12, - \tfrac12, \tfrac12, - \tfrac12)	\notag
\\
\bold T(\gamma_3 Q) &= (\theta'; \tfrac12, - \tfrac12, - \tfrac12, \tfrac12) 	\notag
\end{align}
By Theorem 1.2 of \ccite{SWmathscand}, which requires the condition $0 < 3(3\theta'-1) < 1,$ there exists a Flip invariant projection $f_0$ such that $\bold T(f_0) = (2(3\theta'-1); 0,0,0,0)$. Now we evaluate equation \eqref{phibetacomb} at the projections $1, f_0, Q, \gamma_1Q, \gamma_2Q, \gamma_3Q,$ and recall that $\phi_{jk}$ is always a half-integer on projections, we obtain the respective equations
\[
K + a = \phi_{jk}(1) \in \{0,1\}
\]
\[
12K\theta' - 4K = N_1
\]
\[
K\theta' + \tfrac12 a  + \tfrac12 b +  \tfrac12 c + \tfrac12 d = \tfrac12 N_2
\]
\[
K\theta' + \tfrac12 a + \tfrac12 b - \tfrac12 c - \tfrac12 d = \tfrac12 N_3
\]
\[
K\theta' + \tfrac12 a  - \tfrac12 b + \tfrac12 c - \tfrac12 d = \tfrac12 N_4
\]
\[
K\theta' + \tfrac12 a  - \tfrac12 b - \tfrac12 c + \tfrac12 d = \tfrac12 N_5
\]
for some integers $N_1,\dots,N_5$. Adding the last four of these equations gives
\[
4K\theta' + 2a = \tfrac12 M
\]
where $M = N_2 + N_3 + N_4 + N_5$ is an integer. From $a = \phi_{jk}(1) - K$ we get
\[
4K\theta' - 2K  = \tfrac12 M - 2\phi_{jk}(1).
\]
But since also $12K\theta' - 4K = N_1$ is an integer, after eliminating $K\theta'$ using these two equations, we see that $K$ must be rational. So from $4K(3\theta' - 1) = N_1$ it follows that $K=0$ (or else $\theta'$ would be rational). Therefore $K = N_1 = 0$ and hence also $a=\phi_{jk}(1)$ and $M = 4\phi_{jk}(1)$. Further, the linear equations above involving $a,b,c,d$ become
\begin{align*}
b +  c +  d &=  N_2 - a
\\
b -  c -  d &=  N_3 - a
\\
-  b +  c -  d &=  N_4 - a
\\
-  b -  c +  d &=  N_5 - a
\end{align*}
which imply that $b,c,d$ are half integers since $a=\phi_{jk}(1)$ is 1 for $jk=00$ and is 0 otherwise. Equation \eqref{phibetacomb} has now become
\begin{equation}
\phi_{jk}\beta = a \phi_{00} + b\phi_{01} +  c\phi_{10} + d \phi_{11}
\end{equation}
where $a=\phi_{jk}(1)$ and $b,c,d\in \tfrac12\mathbb Z,$ hence the result.
\end{proof}

\medskip

\begin{qst} As noted above, the invariance relation $\phi_{00}\alpha = \phi_{00}$ holds for any modular $\alpha$ as well as for $\gamma_j$'s. Is it true for all smooth automorphisms $\alpha$ of the Flip orbifold? Or for any smooth automorphism of $A_\theta$  commuting with the Flip?
\end{qst}

\medskip

\begin{cor}
A smooth automorphism of the Flip orbifold that induces the identity map on $K_0(A_\theta^{\infty,\Phi}) = \mathbb Z^6$ necessarily induces the identity map on the zeroth cyclic cohomology group $HC^0(A_\theta^{\infty,\Phi}) = \mathbb C^5$.
\end{cor}

\noindent(Automorphisms homotopic to the identity induce the identity on $K_0$.)

\begin{proof}
Let $\beta$ denote a smooth automorphism for which, by hypothesis, $[\beta(P)] = [P]$ for each projection $P$. In particular, $\phi_{jk}\beta(P) = \phi_{jk}(P)$ for each projection $P,$ and we want to show that $\phi_{jk}\beta = \phi_{jk}$ as maps on the orbifold $A_\theta^{\infty,\Phi}$.\footnote{Although the linear span of projections is norm dense we cannot conclude $\phi_{jk}\beta = \phi_{jk}$ from this since $\phi_{jk}$ are not norm continuous.} From Theorem \ref{betaspan} we have
\[
\begin{bmatrix} \phi_{00}\beta \\ \phi_{01}\beta \\ \phi_{10}\beta \\ \phi_{11}\beta
\end{bmatrix}
= M \begin{bmatrix} \phi_{00} \\ \phi_{01} \\ \phi_{10} \\ \phi_{11} \end{bmatrix}
\]
for some $4\times4$ matrix $M$ (with half-integer entries). Evaluating this at the projections $P$ mentioned in \eqref{gamma1Q}, and using $\phi_{jk}\beta(P) = \phi_{jk}(P),$ we get
\[
\begin{bmatrix} \phi_{00}(P) \\ \phi_{01}(P) \\ \phi_{10}(P) \\ \phi_{11}(P)
\end{bmatrix}
= M \begin{bmatrix} \phi_{00}(P) \\ \phi_{01}(P) \\ \phi_{10}(P) \\ \phi_{11}(P) \end{bmatrix}.
\]
Stacking the left and right sides for the four projections $P = Q, \gamma_1Q, \gamma_2Q, \gamma_3Q$ (and multiplying by a factor of 2) we get
\[
\begin{bmatrix} 1 & 1 & 1 & 1 \\ 1 & 1 & -1 & -1 \\ 1 & -1 & 1 & -1 \\ 1 & -1 & -1 & 1
\end{bmatrix}
= M \begin{bmatrix} 1 & 1 & 1 & 1 \\ 1 & 1 & -1 & -1 \\ 1 & -1 & 1 & -1 \\ 1 & -1 & -1 & 1
\end{bmatrix}.
\]
The matrix appearing on the left and right side has determinant $-16$, so is invertible, hence $M$ is the identity matrix, showing that $\phi_{jk}\beta = \phi_{jk}$ as maps on the smooth orbifold.
\end{proof}

\medskip

This, together with the same proof argument used in Section 2, gives us the following result.

\begin{thm}
Let $\beta$ be a smooth automorphism of the Flip orbifold that induces the identity on $K_0(A_\theta^{\infty,\Phi})$. Then $\beta(e)  {\text{\ \Large$\sim_{{}_{\!\!\!\!\!\!\!\!c}}$\ \,}}   e$ for any approximately central projection $e$ in $A_\theta^{\infty,\Phi}$.
\end{thm}

\smallskip

\begin{rmk}
The smooth automorphism $\gamma_1$ (see \eqref{gammas}) is homotopic to the identity on the rotation algebra $A_\theta,$ but it is not homotopic to the identity on the Flip orbifold since $\gamma_1$ does not induce the identity map on $K_0(A_\theta^\Phi)$ (which is clear from \eqref{gamma1Q}). To show that $\gamma_1(e)$ and $e$ are not centrally equivalent for some AC projection $e$ in the Flip orbifold is not trivial, but can be shown from \ccite{SWb} for an AC Powers-Rieffel projection $e$  (where the same is also shown to occur for the Fourier and Cubic transforms on the Flip orbifold).
\end{rmk}

\begin{rmk}
From Theorem \ref{betaspan}, one obtains a group representation
\[
M: \text{Aut}(A_\theta^{\infty,\Phi}) \to \text{SL}_{\pm}(3,\mathbb R)
\]
into the modular subgroup $\text{SL}_{\pm}(3,\mathbb R)$ of $\text{GL}(3,\mathbb R)$ consisting of matrices with determinant $\pm1,$ such that the $3\times3$ invertible matrix $M(\beta)$ has half-integer entries for each $\beta,$ defined by noting that any given smooth automorphism $\beta$ (of the Flip orbifold) leaves the subspace spanned by $\phi_{01}, \phi_{10}, \phi_{11}$ invariant.  (In fact, all its powers $M(\beta)^n = M(\beta^n)$ also have half-integer entries for each $n\in \mathbb Z,$ from which it is easy to see that $\det\, M(\beta) = \pm1$.) Do the entries all have to be integers? We do not have an example where $M(\beta)$ has non-integer entries (for smooth $\beta$). For instance, for automorphisms $\alpha_1, \alpha_2$ mentioned in Lemma \ref{lemma} (for even $a,d,$ and for even $a$ and odd $d,$ respectively), one has cyclic permutation matrices
\[
M(\alpha_1) = \begin{bmatrix}
0 & 1 & 0 \\
1 & 0 & 0 \\
0 & 0 & 1
\end{bmatrix}, \qquad
M(\alpha_2) = \begin{bmatrix}
0 & 1 & 0 \\
0 & 0 & 1 \\
1 & 0& 0
\end{bmatrix}.
\]
(The first one has determinant $-1,$ so $M(\beta)$ can have negative determinant.)
For $\gamma_1$ it is easy to see that its matrix is $M(\gamma_1)= $ diag$(1,-1,-1)$. It would be interesting to know the homomorphic image of $M$ (which contains the subgroup generated by the permutation matrices and the indicated diagonals).
\end{rmk}

\medskip