abstract: While 5G mobile communication systems are currently in deployment, researchers around the world have already started to discuss 6G technology and funding agencies started their first programs with a 6G label. Although it may seem like a good idea from a historical point of view with returning generations every decade, this contribution will show that there is a great risk of introducing 6G labels at this time. While the reasons to not talk about 6G yet are manifold, some of the more dominant ones are $i.)$ there exists a lack of real technology advancements introduced by a potential 6G system; $ii.)$ the flexibility of the 5G communication system introduced by softwarization concepts, such as in the Internet community, allows for daily updates; and $iii.)$ introducing widespread 6G discissions can have a negative impact on the deployment and evolution of 5G with completely new business cases and customer ecosystems compared to its predecessors. Finally, as we do not believe that 5G is the end of our journey, we will provide an outlook on the future of mobile communication systems, independent of the current mainstream discussion.
author: Frank H. P. Fitzek and Patrick Seeling [^1] [^2] [^3]
title: Delay Impacts on EEG-Based Determination of the Human Visual Interface QoE for Virtual and Augmented Realities; Why We Should NOT Talk about 6G

# Introduction – The Evolution of Mobile Communication Systems

Whether there is a clear need for a sixth generation (6G) of mobile communication systems can only be answered through an initial recognition of the history of mobile communications and the motivations underpinning the continuous development efforts. Mobile communication systems are commonly categorized in generations. (We refer to mobile communication systems throughout this contribution with a focus on *cellular* mobile communication systems.) So far, a new generation has been introduced in each decade. The first generation (1G) of mobile communication systems allowed a small group of privileged users to experience mobile voice services employing purely analog technology . The second generation (2G) democratized these services for the masses exploiting digital technology. Both generations have been straightforward extensions of the existing public switched telephone networks (PSTNs) of their times to mobile services. The main services of 2G, following the global system for mobile communications (GSM) definition, were speech services, mobility, and security. Message and data service implementations were only added to the network afterwards, featuring limited data rates and capacities . With the introduction of the third generation (3G), an extension of mobile communication systems to the Internet was targeted. Even though the first realizations, such as i-mode in Japan or the wireless application protocol (WAP) in Europe, had been great failures, the fourth generation (4G) finally enabled the mobile Internet for users with great success . The introduction of a new generation was always motivated by and aligned with changes in the underpinning technical details. While 1G was based purely on analog technology, the digital era was introduced with 2G. Furthermore, 2G in Europe was based on time division multiple access (TDMA), while code division multiple access (CDMA) was deployed in the USA and other countries alongside TDMA. With 3G, CDMA was introduced as a global standard, significantly unifying the formerly divided technology landscape. Finally, 4G uses orthogonal division multiple access (OFDM). With that, each generation could have been identified by the predominant access scheme. All generations so far have in common that they targeted the consumer market's main demands over time: to enable voice services and mobility in the beginning (starting in the 1980s) up to video streaming and social networking towards the end (the 2010s). Nevertheless, each generation has its entitlement in a technical development stage. The fifth generation (5G) is different , as $i.)$ it continues to use OFDM as it provides technical benefits over other access technologies and there are simply no novel access technologies left, $ii.)$ its standardization efforts from the beginning onward did not just follow a philosophy of pure capacity increases but already addressed new technical parameters, such as latency, resilience, heterogeneity, or energy, $iii.)$ it already has a strong impact on the overall network architecture rather than focusing on a new air interface with some radio access network (RAN) elements, and $iv.)$ the target use cases are mainly in the domain of machine control rather than the consumer market, a major shift with respect to prior generations.

As illustrated in Figure , the initial four generations have been addressing the wireless domain with some additional RAN elements. 5G will follow this approach on the wireless side, but is strongly engaged with the wired domain. One evidence for that is the large number of bilateral meetings from the 3rd Generation Partnership Project (3GPP) and the Internet Engineering Task Force (IETF) to define 5G technologies. These meetings are embossed by a cultural clash where one side is dominated by hardware while the other is dominated by software. While the generations on the mobile communication system side are a result of changing hardware in the time frame of one decade, the Internet domain is realized by software and allows for daily updates. Once the two worlds are truly merged in an holistic manner in the future, software will play a dominant role, while specialized hardware is needed only for special tasks, such as energy efficient computing or energy efficient and secure radio front-ends with antennas.

5G is a highly disruptive communication system, as can already be derived from only some of the main differences to its predecessors, as highlighted in the following.

Low Latency as a Service:

:
    Besides pure data rate, 5G also targets low latency communications. This enables communication for time-critical scenarios. In other words, 5G combines control and communication theory. While the predecessors target the consumer market for communication needs of humans, 5G is mainly targeting machine-type communications. Even though 5G is currently introduced as a bit pipe, with upcoming releases 16 and onward, the new service area will be the target. Therefore, 5G is able to become the communication network for a broad variety of use cases, such as transportation, Industry 4.0, Agriculture 4.0, construction, health care, or the upcoming Tactile Internet. The Tactile Internet has the most severe latency requirements with $1$ ms end-to-end delay including the embedded system, the wireless link, and the computing within the wired network. As given in , the most challenging entity for low latency are network functions realized in software.

Softwarization:

:
    5G is not just an agnostic bit pipe as the generations before. This is achieved by the concept of softwarization. While its predecessors came in the form of dedicated hardware, 5G is mainly based on generic hardware with software solutions. Different concepts, such as software defined radio (SDR), software defined networks (SDN), and network function virtualization (NFV), enable network operators to keep operational expenses (OPEX) as well as capital expenses (CAPEX) low. Additionally, softwarization enables network operators to dynamically reorganise their networks for different current and future use cases. New concepts such as mobile edge cloud and network slicing are only possible due to the introduction of the softwarization concept.

5G as Campus Solutions:

:
    Due to the last point, the realization of so-called 5G campus solutions is now possible. In some countries, industry or government entities can apply for 5G frequencies for a given local area. A 5G campus solution operates in an isolated geographic area and provides localized services. Mobility support is not needed in this application scenario, but ultra-low latency communication is.

These three main differentiators of 5G with respect to prior generations already highlight how it completely differs from its predecessors. It is the first generation to target machine-type communications, even though the support of multimedia services is still an option. 5G will support different technical parameters such as low latency, security, heterogeneity, and more. More importantly, the replacement of specialized hardware with generic hardware in combination with dedicated software solutions is such a fundamental change, that it will have long-term implications for the future of communication networks.

Based on the paradigm shift towards network as a software solution, the architecture of 5G will be extended by campus solutions. Those campus solutions do not necessarily support mobility in the cellular WAN context, but rather offer mobility in the LAN context, which is more nomadic. While earlier generations defined communication as a transport of bits leading to the end-to-end paradigm of communication, with 5G we enter a new definition where communication is defined as transport, storage, and computing of bits within the network. The disruptive nature of 5G is that it elevates the network from data to information.

Despite these revolutionary changes, it is highly likely that 5G will face immense problems in the beginning of its deployment phase. Classic services, such as video streaming, voice calls, etc., are part of the general update of capabilities from 4G to 5G – but focusing on their contribution to business success might be more than disappointing. This should not surprise, as the superiority of 5G lies in low latency communications, which will become part of implementation rollouts with the standard release 16 and onward. Once the implmentations find their way into associated use cases, the business side of the success story of 5G will start as well.

# The Atom – A Holistic View on 5G

In this section, we briefly discuss the 5G landscape from a holistic viewpoint by employing the *5G Atom*. Figure  illustrates this view, which deviates from the *massiveness* in the International Mobile Telecommunications (IMT) vision  of massive IoT, massive multimedia, and massive low latency, by following 5G use cases defining higher layers, up to novel applications.

The *massive* use cases identified in IMT-2020 are part of the initial wave of 5G use cases, but they do not account for all of them. One can readily separate the use cases by their focus of communication, whereby the machine-to-machine communication will be expressed in roll-outs in connected autonomous cars, Industry 4.0, Agriculture 4.0, energy grids, as well as construction, health care, or education – in what is heavily marketed today as the Internet of Things (IoT). In contrast to these, human-machine-type communications that are summarized as the Tactile Internet  will feature ultra-reliable low latency communications (URLLC)  enabling human-machine collaborations across all sectors of industry and in mixed realities. The Tactile Internet has the potential to usher in a new paradigm shift from a network of information to a network of skills  with a focus on layer 8 of the network.

A readily notable technical requirement resulting from these different use cases is the handling of increased heterogeneity of communication types and challenges when considering different communication actors and device types. For example, connecting massive numbers of IoT devices in the billions requires new approaches to data handling and storage. Additionally, 5G networks need to feature security as well as resilience, not just with respect to adverse actors threatening the network, privacy concerns, or interference. Maintaining quality of service (QoS) levels to higher degrees than before for at least a slice of the network is of utmost importance as use cases feature control loops that are comprised of machines as well as humans, with potentially dire consequences as results of service degradation. For example, URLLC requires significant reductions of the latency present in today's networks down to a millisecond, while simultaneously bounding delay variations (jitter). As the case with all prior generations, increases of throughput are common requests in use case evolutions. Some of these move from today's dominant multimedia streaming services to fully immersive experience streaming services requiring significantly more data. With such a broad landscape of devices and use cases, energy consumption is an additional factor, e.g., to reduce battery depletion rates and required maintenance for devices in the field.

The following concept layer in Figure  highlights approaches to resolve these often diametrically opposed technical requirements. The 5G new radio air interface will significantly reduce latency, offer higher throughput, and enable massive numbers of securely connected devices. In addition, new concepts that offer extensions of the 5G network through meshed networks as well as multi-path communications can aid significantly to realize the technical requirements as part of a softwarized solution. In addition to networking-centric approaches, the softwarization inherent in 5G also is tied to new approaches to managing the network, such as focusing on information rather than just data, and to dynamically allocate the network resources (including computation and storage) through network slicing. This approach enables an ability to dynamically move these resources at the network's edge, as mobile edge cloud (MEC).

These solutions make heavy use of softwarization technologies across the 5G network by employing software-defined radio (SDR), software-defined networking (SDN), network function virtualization (NFV), and service function chaining (SFC). Jointly, these technologies allow a dynamic allocation of all resources within a 5G network and their logical combination. This makes the common computing resources (communication, computation, and storage) flexible and versatile to perform flexible service composition.

The developed services can readily be altered to reflect changing use cases as well as new use cases providing network operators with long-term flexibility as well. The outer layer application examples in Figure  that will likely become part of these future services include block chains, machine learning, network coding, and compressed sensing. While a detailed discussion of these are out of scope here, we note that all of these have been extensively researched in conjunction with networked applications. It is therefore foreseeable that they will become part of future service offerings – or not, as the flexibility offered by the softwarized network can readily support others not even considered yet.

# Mythbusting 6G

After the introduction of 5G, in this section, we present several theses which question the necessity of a 6G technology.

## Nothing New in 6G

Currently, there is nothing to little novelty proposed for 6G that can not be adopted or realized in 5G networks and its evolution plan. Often new frequencies, such as terahertz frequencies, are labeled as 6G. But the 5G evolution foresees already higher frequencies up to 60 GHz and to include even higher frequencies would not require a new generation. Even a change of hardware or software components is not necessarily leading to a new generation – 2G already supported channel bundling (software) and novel modulation schemes (hardware). Some research labs are even proposing artificial intelligence or mobile edge cloud as a 6G feature, but this has been already introduced in 5G. Oftentimes, lower latencies than promised in 5G are also referred to when discussing 6G. Some researchers claim that higher data rates are needed to achieve this goal, others rely on a novel medium access protocol in their claims. In general, the data rates of 5G are already so high that the propagation delay increases in importance when compared to the relative delay, as  have shown for simple M/M/1 queues. In conclusion, there is little in individual technologies that would require a new generation. General evolution steps as well as disruptive ideas can be realized with software changes within 5G networks, as introduced beforehand.

## The Internet Has No Generation, Why Mobile Networks?

Given the evolutionary changes in the roll-outs of cellular network generations to date, a common thread is the combination of typically specialized hardware combined with specialized and provider- or operator-based services targeting cellular end-systems used by humans. This development increased capacities on the network side to support mobile Internet based over-the-top (OTT) services at the end-units (e.g., smartphones for media consumption). The fixation on continuously evolving services for individual human customers meant that architectural changes were required to fully enable Internet-based services following the end-to-end paradigm.

At the same time as mobile communication network operators switched out generations of different custom hardware and software solutions moving from 3G to 4G, the smartphones they were connecting to the Internet were already softwarized. The software-oriented approach inherent to common smartphones enables the continuous addition of new features based on a common set of hardware components via software updates.

Throughout the history of the Internet, the closest notion of generational change can be seen in the versions of protocols, predominantly the migration from IPv4 to IPv6. The reason that the sunset of IPv4 takes so long is also based on hardware legacy. Nevertheless, researchers in standardization groups such as IETF would never use the term generation as software enables agile testing and implementation. The overall nature of "Everything over IP and IP over Everything" abstracts or even virtualizes the underlying complexity for upper layer protocols. In turn, changes can be made on all layers above IP as needed in a fluid fashion.The softwarized nature of even a significant portion of the core network components means that within reason, new features can now be added via system software upgrades on general-purpose hardware and no longer require full system hardware swaps, i.e., a new generation.

As the network role is elevated from data to information, it becomes clear that only individual components need changes over time, as 5G already contains the built-in flexibility of softwarization of the network. Thus, it is difficult to justify individual component tweaks to modify network capabilities to be labeled as a generational change. Indeed, intuitive examples such as replacing an older general-purpose server in a rack with a more powerful could increase the in-network computation speed (e.g., to support more services offered through a MEC or in a nomadic 5G campus container), but hardly qualify as a new generation. Similarly, significantly changing TCP's congestion algorithm did not result in new protocol versions despite significant performance changes. The network softwarization of 5G enables continuous, fluid change that reduces the need for generational changes.

## 5G Engaged Vertical Players – 6G Scares Them Away

5G was very successful to engage vertical players even before the start of the technology itself. While 3G failed in identifying possible markets and 4G took over the promise of 3G, 5G breaks new grounds for a wide range of heterogeneous markets. Some markets such as agriculture or construction happily engaged with the 5G technology. Other market segments such as connected cars or Industry 4.0 had alternatives such as IEEE 802.11p or field busses. It took, and still takes, some effort to explain the advantages of 5G over existing solutions. This was successful for most use cases and is ongoing for others. The upcoming discussions of the 6G technology creates some nervousness among the management in industry. When considered by merit and potential, there should be a Fear Of Missing Out (FOMO), but these discussions more likely begin to create a Fear Of Joining In (FOJI). Some might think that 5G is already an old technology and 6G is the holy grail now. This has the potential to destroy the momentum 5G has. Distracting the industry from 5G could lead to a decreased market adoption. Subsequently, this would put both network operators and network equipment suppliers in great financial distress. Without the revenues earned from 5G, a new generation, possibly 6G, would not even see the light of day.

## Only Even Numbered Generations Are a Success

In several keynotes about 5G and 6G, the speakers claim that only even numbered generations are a success. Such statements are not only ridiculous but also technically wrong. These statements are based on the difficulties 3G had in the beginning of its introduction, while 2G and 4G were a success story from the beginning. However, this argumentation has several flaws. The first generation 1G was undoubtedly a success even if not available for the masses. Furthermore, the problem with such a statement is the statistical confidence one can achieve with only five samples. Finally, if the statement were true, we should just use even numbered generations from the beginning.

## The Use Case, Stupid

Making a numbers-only argument additionally obfuscates a significant shortcoming in the overall interplay of 5G and 4G with respect to use cases. Through complete softwarization of the network, components can be changed dynamically and without a need to call for a new number. It is noteworthy to consider the scenarios for which these last two generations were designed: mobile Internet user-centered in 4G (over-the-top media consumption by humans) and expanding with (mobile) mission-critical with ultra-low latency in 5G (massive machine-type communications, including the Tactile Internet). Some use cases for 5G, such as the outlined 5G campus, offer the seamless integration of solutions into a locally managed 5G solution, but no mobility. Jointly, these cases cover the currently imaginable scenarios for implementation fairly well, as listed in Table .

```latex
\begin{table*}
\centering
\caption{Minimum requirements (in terms of generations and use cases) for human- and machine-type communications. \label{tab:one}}
\begin{tabularx}{\textwidth}{|l|l|X|}\hline
Mobility target    & Human                     & Machine  \\\hline
Mobile              &3G/4G for browsing   & 5G for cars, drones (for both mobility is needed), energy grids (coverage is needed), also for the Tactile Internet \\
                    &4G/5G for over-the-top services      & \\
                    &5G for the Tactile Internet      & \\\hline
Nomadic/stationary  & WiFi                      & 5G campus, OpenRAN, also for the Tactile Internet\\
                    &5G, also for the Tactile Internet & \\\hline
\end{tabularx}
\end{table*}
```

Differently put, there are no current new scenarios that are not already covered through 5G or with a potential incorporation of non-cellular approaches. The remaining generational updates could be attributed to radio access network upgrades. While additional arguments could be made that increased radio access network capabilities would yield better front-end latency characteristics, the question is how much is enough for use cases? As usual when moving into business considerations, one also needs to carefully consider the diminishing rates of return for increased resource expenditures. Realizing higher throughput to reduce transmission delays (assuming that propagation delays cannot be influenced in the near future, e.g., with quantum networking), higher frequency bands are needed. As a direct downside, this requires significant attenuation handling and will reduce the current small cell sizes for going beyond mmWave bands even further, which might impose the need for a significant and cost-prohibitive increase in the number of cells needed for sufficient coverage.

# So, What *IS* the Future of Mobile Communication Networks?

While we advocate against the 6G label for the moment, we still believe in fundamental changes in future communication systems. As illustrated in Figure , the first (cellular) mobile communication systems leading to 4G target humans and their mobility. Beginning with 5G, the focus shifts to being able to support communications with and between machines in real-time. As the 5G standard continues to evolve jointly with applied use cases, its evolution will continue with future releases 16, 17, …, simultaneously incorporating additional features for deployments.

In contrast to the common prior focus on use cases incorporating mobility, 5G represents a turning point for the support of stationary, localized communications. Similar to prior generations, 5G will continue to support the mobility of humans as well as machines, such as cars or drones. The URLLC components of 5G can even be seen as one enabler for future mobility scenarios with connected autonomous cars. Mobility could also refer to the scope of information needs, which can readily go beyond single-cell sizes when considering the organization of large, distributed deployments, such as the energy grid. While information mobility support is highly needed for these types of use cases, the emergence of human-machine co-habitation for the future of work and life represents a new scenario. Here, URLLC will support the Tactile Internet-based and localized cooperation of humans and robots within their similarly networked but transparently operating local environment.

In these localized scenarios, mobility is just needed in a nomadic context common to WiFi networks, e.g., as provided by EduRoam for academic users that move between campuses, and the focus will be on URLLC. 5G will always be superior over WiFi (i.e., the IEEE 802.11 WLAN standard range) products, even WiFi6, when it comes to dedicated latency requirements. WiFi's core philosophy of "listen before talk" clashes with deterministic latencies. Similarly, potential enhancements to WiFi protocols still require CSMA to operate in the general ISM bands, see, e.g., , which will continue to result in a no-guarantee upper timing bound unsuitable for mission-critical application scenarios requiring URLLC. Nevertheless, the cost efficiency of WiFi networks that maintain their ISM band reliance is still remarkable and readily serves as a target for future communication systems. This future might well be in the combination of the strengths of the two aforementioned technologies, namely 5G and WiFi6.

Subsequently, it is highly likely that we will witness two main trends, namely the evolution of 5G technologies and the rise of a new communication system that is not categorized in a generation. This new communication system will target specific use cases that do not need mobility support and national coverage to a large extent. Here we see Industry 4.0 or human–machine interaction as possible use cases. This can be seen as a low cost variant of the 5G campus idea, where the wireless technology is realized by OpenRAN solutions.

The flexibility of OpenRAN in combination with network slicing ideas of the core network will enable ultra-low latency communications that is also resilient against hostile agents. For example, adverse activities could include denial of service through extensive jamming of individual cells or modifying the data exchanged. The latter case is particularly interesting in the context of URLLC, as sophisticated mitigation techniques might readily impact the ability to remain at the ultra-low latency level, e.g., as required for the 1 ms round-trip (including all sensing, computing, and actuating next to network transport) deadlines imposed by the Tactile Internet. In the context of collaborative robotics in Industry 4.0 use cases, one result could be that adversaries become enabled to tamper with the actual payload without detection, e.g., through knowledge of the actual information content, for an individual robot's control and cause undetectable extensive damage. One opportunity that OpenRAN offers with software-defined radios is to change the underlying MAC protocols at a high frequency. (We note that this is not the same approach as security through obfuscation, which is commonly considered harmful, but a new approach to security enabled by softwarization.) A similar trend to enable the lower latency bounds imposed by use cases is the need to move computation from the cloud to the edge of the network. Here, significant computational resources will need to be dynamically allocated and re-partitioned based on use cases. Again, this can only be realized through softwarization. Further latency reductions for actual use cases that take the shift to information-based mobile communications into account will, in turn, more likely be based on commodity hardware upgrades at the network edge and algorithmic improvements in software, as ultra-low latency dependent services will not be able to leave an operator's network. Even more stringent latency requirements could be realized in 5G campus offerings with the possibility of significant embedded computational resources at its center (if an individual single cell) or edge (in general).

Considering the combination of flexible, commodity hardware paired with fully softwarized access and core networks as well as OpenRAN for the last hop, it becomes clear that 5G generates the ability for network operators to offer everything as a service (\*aaS). For example, general latency as a service can be realized through flexible resource slicing on the network and computation side. Upgraded, guaranteed latencies can be provided through software-based OpenRAN changes that implement a bus-based polling protocol of subscribed stations, resulting in upper latency bounds. Storage, in-network computation, resource pooling, and, ultimately, services as a service (SaaS) all become enabled through elevation of the network to an information level, where an interplay of wireless access technologies could be orchestrated through 5G.

# Conclusion

In this article, we made a case for future mobile communication network research needs. Research and development must focus on the evolutionary advancement of the 5G roadmap. This roadmap provides many improvements in the areas of spectrum usage, integration of IETF software solutions for the mobile edge cloud, network slicing, and artificial intelligence. In parallel, there will be a disruptive evolution in the architectures of future mobile communications systems, targeting very specific applications, such as human–machine cooperation and Industry 4.0. The envisioned mobile radio systems will no longer be cellular and will only provide local coverage. The latter will allow low latency applications, novel security concepts, and extremely low cost. This approach could be seen as a symbiosis of the existing 5G campus solutions and WiFi6. Nevertheless, we have also argued in this article that the use of the 6G label is not only counterproductive, but it also is not even technically motivated. Using 6G is counterproductive where new markets have just been opened up by 5G. The players in these new markets could be very nervous about investing in 5G, when new technology may soon be available. For example, the mere use of the term 6G may prevent the 5G market from fully developing. From a technical point of view, there is little reason to talk about a further generation in mobile communications because with 5G we leave the area of hardware behind us and now mostly realize everything in software. Thus, mobile communication systems are reaching the status that the Internet community has held for decades: continuous evolution and improvement through software. Nowadays, the 6G term is more a marketing tool that might end up in a 6G label without technical needs.

[^1]: F. H. P. Fitzek is with the Centre for Tactile Internet with Human-in-the-Loop and the Deutsche Telekom Chair of Communication Networks, Technische Universität Dresden, Dresden, Saxony, Germany, email: email@example.com

[^2]: Patrick Seeling is with the Department of Computer Science at Central Michigan University, Mount Pleasant, MI, USA, email: email@example.com

[^3]: Funded by the German Research Foundation (DFG, Deutsche Forschungsgemeinschaft) as part of Germany's Excellence Strategy – EXC 2050/1 – Project ID 390696704 – Cluster of Excellence "Centre for Tactile Internet with Human-in-the-Loop" (CeTI) of Technische Universität Dresden.


author: S.E. van Terwisga; E.F. van Dishoeck; M. Ansdell; N. van der Marel; L. Testi; J.P. Williams; S. Facchini; M. Tazzari; M.R. Hogerheijde; L. Trapman; C.F. Manara; A. Miotello; L.T. Maud; D. Harsono
bibliography: V1094_bib.bib
title: V1094 Sco: a rare giant multi-ringed disk around a T Tauri star

To characterize the structure of an extended, multi-ringed disk discovered serendipitously in the ALMA Lupus disk survey and put it in the context of the Lupus disk population. ALMA observations in Band 6 at 234 GHz and Band 7 at 328 GHz at 0.3$''$ resolution toward the K6 star V1094 Sco in Lupus III are presented, and its disk structure analyzed. The spectral index $\alpha_{\rm{mm}}$ is determined in the inner $150$ AU of the disk. The ALMA continuum data show a very extended disk with two gap/ring pairs. The gaps are located at 100 AU and 170 AU, the bright rings at 130 AU and 220 AU. Continuum emission is detected out to a 300 AU distance, similar to IM Lup but a factor of 5 larger than typically found for Lupus disks at this sensitivity and resolution. The bright central region of the disk (within 35 AU) is possibly optically thick at 1 mm wavelengths, and has a brightness temperature of only 13 K. The spectral index increases between the inner disk and the first ring, at the location of the first gap. Due to the low temperature of the disk midplane, snow lines can be excluded as the drivers behind the ring and gap formation in this disk. Disks the size of V1094 Sco are rare, and only $2.1 \pm 1.5 \%$ of disks in Lupus show continuum emission beyond 200 AU. Possible connections between the large primordial disk population, transition disks, and exoplanets are discussed.

# Introduction

One of the key discoveries of the Atacama Large Millimeter/submillimeter Array (ALMA) has been the unexpectedly wide variety of structures seen in protoplanetary disks, ranging from (asymmetric) dust cavities to multiple rings at a range of radii. These multiple rings have been found in disks around both Herbig Ae/Be and T Tauri stars (). Of these T Tauri stars, HL Tau, Elias 2-24, and AS 209 are very young objects ($< 1$ Myr), while the massive TW Hya disk is unexpectedly old (up to 10 Myr). The bulk of the disks imaged in continuum surveys of nearby star-forming regions such as Lupus does not seem to show substructure at 0.2 – 0.3$''$ resolution, apart from transition disks with large inner dust cavities ($> 20$ AU) .

The existence of multi-ringed disks raises important questions on the evolution of disks, and therefore on planet formation, since dust particles trapped in a ring may grow efficiently without being accreted onto the star . Thus, rings can provide the material necessary for possible future planet formation. So far, several mechanisms have been proposed in the literature to generate ring-like continuum structures: (secular) gravitational instabilities, the effect of snow lines of various molecules, dead zones and, most tantalizingly, the presence of already-formed planets in the intervening gaps . Our understanding of the physics behind these multi-ringed disks is limited by the small number of known examples ($N = 5$ around T Tauri stars), as well as the difficulty of placing them in the context of a well-studied, similar population. Moreover, it is still not clear how common it is for disks to show multiple continuum rings, whether they correlate with stellar or disk parameters, or when in their evolution features like these are most likely to be present.

Another, less obvious structure found in many protoplanetary disks is an intensity profile with a bright central region surrounded by a fainter, extended outer region, resembling a "fried egg". This structure can be seen both in multi-ringed disks, such as AS 209 and HD 163296 (), and in full disks where no rings have been identified, such as V883 Ori . In the well-studied case of V883 Ori, the origin of the morphology is due to a steep (spatially unresolved) increase in optical depth at the location of the knee in the intensity profile. This increase may be linked to the presence of a water snowline, inside of which the originally large icy grains fragment into many smaller bare silicate grains. These smaller particles have smaller drift velocities causing a 'traffic jam' and thus a dust pile-up in the inner disk increasing the optical depth . More generally, this type of intensity profile must trace the underlying grain properties (number density, composition and size, or both).

In this paper, we present ALMA Band 6 and Band 7 observations of V1094 Sco, a newly-identified giant multi-ringed disk in Lupus. V1094 is a K6 star with a luminosity of 1.7 (1.95) L$_{\odot}$ and an accretion rate of $\log( \dot{M}_{\rm{acc}} / (M_{\odot} \rm{yr}^{-1}) ) = -8.12\,(-7.7)$ . Throughout this article, we take a 150 pc distance to the source, consistent with the latest *Gaia* distance estimate of $153 \pm 1$ pc ; the values between brackets correspond to the published values assuming a pre-*Gaia* 200 pc source distance. The V1094 disk had previously been recognized as an enigmatic object, in particular by , who inferred from spectral energy distribution (SED) modeling that it was a cold and massive disk. Photometry from the *Herschel Space Observatory* shows a large far-infrared excess, with an SED classification inconsistent with that of a transition disk . Previous direct imaging searches for companions, either stellar or of planetary mass, have not found any down to $\sim 10\,M_J$ beyond 75 AU .

Several key aspects of this disk make it an especially valuable addition to the so far small sample of multi-ringed disks around T Tauri stars: first, the population of Class II Young Stellar Objects in the Lupus clouds is exceptionally well-studied, with a $100\%$ coverage in ALMA Band 6 and 7 continuum at $0.3''$ resolution and $0.3$ mJy beam$^{-1}$ sensitivity . $80\%$ of these sources, including all the objects discussed here, also have VLT/X-Shooter spectroscopy data available resulting in accurate stellar types and accretion properties . This enables us to compare the disk and its host star to a large, homogenous sample, and calculate the occurrence rate of this type of object at the age of the Lupus cloud. Moreover, at $\sim 3$ Myr, the age of this disk places it in the middle between the young ($<1$ Myr) objects, such as Elias 2-24 and HL Tau, and $10$ Myr-old TW Hya disk .

In §2, the ALMA observational details are discussed, whereas continuum images showing a bright central millimeter emission core with two outer gaps and rings are presented and analyzed in §3. The availability of data at two frequencies allow the spectral index as function of radius to be determined. The V1094 Sco disk is compared with other similar disk structures in §4.1, and put in the context of the full Lupus disk sample in §4.2, demonstrating that such large disks are rare.

# Observations and data reduction

The data presented here are part of the Lupus completion survey (ALMA proposal 2016.1.01239.S, P.I.: S.E. van Terwisga) which aimed to provide data on the sources not previously included in the ALMA survey of Lupus protoplanetary disks and for which no data at similar spectral and resolution settings were available at the time. The sources included are EX Lup, GQ Lup, RX J1556.1-3655, Sz 102, Sz 76, Sz 77, and V1094 Sco.

Our observations cover Band 6 and Band 7 continuum, and the $^{12}$CO $2-1$, $^{13}$CO $3-2$ and $2-1$, and C$^{18}$O $3-2$ and $2-1$ lines, as well as CN $3-2$. The Band 6 observations were performed on 7 July 2017, using 44 antennae and $2600 - 16700$ m baselines, and a PWV of around 0.64 mm. The sample was observed in Band 7 on 21 May 2017, using 41 antennae and baselines between 1100 m and 15100 m at PWV columns of $\sim$ 0.7 mm. On-source integration times were 3 and 4.5 min in Band 6 and 7, respectively.

The Band 6 data were flux-calibrated using J1427-4206, and had J1610-3958 as phase- and J1517-2422 as bandpass calibrators. The Band 7 data used J1517-2422 as flux- and bandpass calibrator and J1607-3331 and J1610-3958 as phase calibrator. While no additional flux calibration on Solar-system objects was done, the Band 6 data were observed on the same day as their flux calibrator, and therefore have a small absolute flux error ($\sim 10 \%$) typical of ALMA. However, the Band 7 flux calibrator was both highly variable (by a factor of 1.5 in the month following the observations) and not observed close to the date on which our data were obtained. As such, its flux — and therefore the flux calibration of all Band 7 data — is very uncertain, complicating the calculation of spectral indices and disk masses.

Fortunately, the Band 7 data in this proposal also cover GQ Lup, which has previously been observed in that band at high $S/N$ and with excellent flux calibration (on Titan and Ceres) by . In addition to their accurate fluxes (a less than $10\%$ systematic uncertainty is quoted), these authors do not find any evidence for variability of this source at mm-wavelengths, retrieving fluxes consistent with previous studies in the $1997 - 2015$ period (). Because our observations covered GQ Lup and V1094 Sco on the same day, it is safe to compare these data with the  results to re-calibrate the Band 7 fluxes. This gives a flux ratio of $F_{\rm{MacGregor\,et\,al.}} / F_{\rm{this\,work}} = 1.3 \pm 0.009$, which we will implicitly apply to the Band 7 data throughout this paper.

Phase-only self-calibration was used to maximize the $S/N$ of the observations in both bands. The Band 6 data allowed self-calibration down to the integration interval ($6$ s), whereas for Band 7 we found useable solutions down to a $15$ s interval.

Imaging the continuum data using the `CASA` `tclean` task with Briggs weighting with a robust parameter of 0.5, the beam shapes become $0.25'' \times 0.22''$ and $0.19'' \times 0.18''$ in Band 6 and Band 7 with rms noise values of 0.11 and 0.3 mJy beam$^{-1}$ respectively for V1094 Sco.

The discussion in this article will focus on the high-$S/N$ continuum observations. Fluxes for the CO (isotopologue) lines can be found in , and CN is discussed in Van Terwisga et al. (in prep.). The gas emission covers a (very) extended region, so that it is partly resolved-out. Also, the $S/N$ of the line data is low, with just a few minutes of integration time, so they are not useful for a spatially resolved analysis.

# Results and analysis

Figure  presents the Band 6 and Band 7 continuum images resulting from the data reduction procedure described above. The continuum rings and bright central core are immediately obvious, especially in the Band 6 image. Also evident is the large radial extent of the source, with continuum emission detected up to 300 AU when placed at 150 pc. This is comparable to the radius of the well-known massive IM Lup disk , one of the largest and most massive disks in Lupus. The full integrated flux (inside of a $2.16''$ radius aperture, beyond which the image is entirely dominated by noise) in Band 6 is $180 \pm 14$ mJy, while in Band 7 it is $390 \pm 19$ mJy. The fluxes and errors here were derived using aperture photometry with a curve-of-growth method to set the aperture radius, identical to the method used in .

The key features are clearly the same in both images: the disk has a large radial extent, with two bright rings superimposed over a fainter, extended continuum, and a bright and at this resolution apparently featureless central region. Like other multi-ringed disks, V1094 Sco does not show significant azimuthal asymmetries at this sensitivity and resolution. Assuming the disk is indeed purely radially symmetric, fitting a two-dimensional Gaussian in the image plane with the `imfit` task in CASA led to a beam-deconvolved PA$= 109 \pm 1.9 ^{\circ}$ and $i = 53 \pm 1.3 ^{\circ}$, and a J2000 phase center of 16h08m36.17s, -39$^{\circ}$`<!-- -->`{=html}23m02.87s.

Comparing the Band 6 and Band 7 data in the image plane, however is not entirely straightforward. Due to the lower $S/N$ and worse weather conditions affecting the Band 7 data, the phase-calibration is only partially successful. As a result, the Band 7 phase noise is higher than that of the Band 6 observations. Also, the $u,v$-plane coverage of the Band 7 data is different. These conditions particularly affect the emission in the North-East and South-West of the disk at $\sim 1''$ from the disk center, and makes the presence of the rings less obvious. The Band 6 data, on the other hand suffer from a low-amplitude North-South ripple, since the longest baselines are only sampled in a single direction. However, its impact is relatively minimal. Therefore we choose the Band 6 data as the basis for the subsequent analysis.

## Disk structure: core, gaps and rings

In order to compare the continuum structure of V1094 Sco with that of other, similar disks around both earlier and similar spectral-type stars, we model the intensity profile using an approach similar to that of : the intensity at radius $R$ is considered as the sum of Gaussian contributions, multiplied by a cosine. This model is then Fourier-transformed and fit to the derotated, deprojected data in the $u,v$-plane, after binning the data to $10$ k$\lambda$ bins. The advantage of this method is that it uses a minimal number of parameters and can easily formally confirm the presence of structures, and precise locations of the brightness maxima and minima. Moreover, it does not rely on a model-dependent temperature profile or assumptions about the dust properties. The $u,v$-plane data show four local maxima between 0 and 700 k$\lambda$; therefore, we use up to $n=4$ terms, for a total of 13 free parameters.

$$\label{eq:zhangmodel}
\begin{split}
    I(\theta) = & \frac{1}{\sqrt{2 \pi} \sigma_0} \exp \left( - \frac{\theta^2}{2\sigma_0^2} \right) \\
    & + \sum_1^3 \cos (2\pi\theta\rho_j) \times \frac{a_j}{\sqrt{2\pi}\sigma_j} \exp \left( -\frac{\theta^2}{2\sigma_j^2} \right).
\end{split}$$

In this model, $\theta$ is the angular separation from the source center, $\sigma_i$ the width of the $i$th component, $a_i$ its intensity and $\rho_i$ the spatial frequency of the cosinusoidal perturbation. No baselines longer than $700$ k$\lambda$ are included in the fit, since the data are noise-dominated beyond that point. A Monte Carlo Markov Chain (MCMC) fitting procedure was used to infer the best fit parameters, using the `pymc` package. Since the model is degenerate if the ranges for $\rho_i$ overlap, their allowed values were restricted based on the location of the peaks in the u,v-plane, but always at least $100$ k$\lambda$ wide in order to sample the full width of each peak. We confirmed that the behavior of the Markov chain was proper, by ensuring the walkers were burnt-in and inspecting their sampling of the parameter space, and checking that the posterior distributions for the parameters showed a single clear maximum. Since we start our model with values for the disk inclination, position angle, and coordinates, derived with the `imfit` task, the resulting best-fit parameters are conditional on these priors. We confirmed their impact on the resulting fit was minimal by re-running the fitting code with different initial choices of inclination, position angle, and source center within their error bars.

Figure  presents the resulting model of the intensity distribution for the Band 6 data, with the best-fit values in Table  and inferred positions for the gaps and rings in Table . The best-fit ring and gap positions are also shown on top of the Band 6 data in Figure . The results show that our 4-component model is a good description of the data, fitting all local maxima within this baseline range. The parameter values also indicate our description of the disk, as a bright core with an extended outer region featuring two rings, is accurate: the 'core' of the disk is contained within the $i=1$ component, having a characteristic radius of $35 \pm 2$ AU.

In the best-fit model, the two bright rings and two gaps are also present in the image plane, arising from the combination of the various components in Equation . They are identified by taking local maxima and minima of $I(\theta)$. Taking into account the errors on the fit parameters, the locations of the bright rings become $130 \pm 9$ AU and $220 \pm 16$ AU, while the gaps are found at $100 \pm 6$ AU and $170 \pm 13$ AU. This closely matches the intensity profile of the disk after deprojecting and derotating in the image plane (see also § ). The best-fit model gives an intensity ratio of $2.3$ for the first ring/gap combination and $1.4$ for the second ring and gap. It is difficult to constrain the width of the gaps and rings from this fit, as there may be structure at smaller scales that is not probed, due to insufficient $S/N$ at large baselines and the resolution limit.

```latex
\begin{table*}[!ht]\caption{Parameter values and inferred errors from the MCMC fit to the Band 6 $u,v$-data.}
\label{tab:uvfit}

    \centering
    \begin{tabular}{l l l l l l}
    \hline \hline
                          & $i=0$       & $i=1$        & $i=2$          & $i=3$         & $i=4$          \\
    \hline
    $a_i$ [Jy]            &  -          &$1.8 \pm 0.3$ &$0.64 \pm 0.096$&$0.46 \pm 0.05$&$0.52 \pm 0.18$ \\
    $\sigma_i$ [AU]       &$161 \pm 8.3$&$35 \pm 1.9$  &$57 \pm 5.0$    &$130 \pm 37$   &$120 \pm 29$    \\
    $\rho_i$ [k$\lambda$] &  -          &$20.7 \pm 0.1$&$446 \pm 12$    &$146 \pm 15$   &$255 \pm 9$     \\
    $ \lambda / \rho_i $ [AU] & -       &$1495 \pm 7$  &$69\pm1.9$      &$212\pm22$     &$121\pm4$       \\
    \hline
    \end{tabular}
\end{table*}
```
```latex
\begin{table*}[!ht]\caption{Ring and gap positions from the MCMC fit to the Band 6 $u,v$-data.}
\label{tab:imres}

    \centering
    \begin{tabular}{l l l l l l}
    \hline \hline
             &  Gap 1  &  Ring 1 &   Gap 2  & Ring 2   \\
    \hline
    R [AU] &$100\pm6$&$130\pm9$&$170\pm13$&$220\pm16$\\
    \hline
    \end{tabular}
\end{table*}
```

## Spectral index analysis

The continuum observations of V1094 Sco show that in both Band 6 and Band 7, the structure is dominated by two azimuthally symmetric regions: the bright, inner 'core' and a fainter, more extended outer disk, where the rings are located. Directly linking intensity to dust properties is difficult, but necessary for tying these observations to models of ring formation. Thanks to the range in frequencies covered by our observations (between 234 and 328 GHz), the spectral index $\alpha_{\rm mm}$ of the continuum emission can be determined. In particular, if the intensity at a given radius is given by $I_{\nu} = B_{\nu}(T)(1 - e^{-\tau_{\nu}})$ and $\tau \sim \nu^\beta$, then if $\tau \gg 1$ and $h\nu / k_B T \ll 1$, $\alpha_{\rm mm} = 2$, while if $\tau \ll 1$, $\alpha_{\rm mm} = 2 + \beta$. However, these simple relations only hold in the Rayleigh-Jeans regime. If this does not apply, an extra (negative) temperature-dependent term appears, and values of $\alpha_{\rm{mm}}$ lower than 2 can occur.

The resolution and sensitivity of our data allow us to determine the spectral index as a function of radius. However, we cannot immediately use the results presented in Figure , since their imaging beam sizes and the $u,v$-plane coverage of the underlying data differ. To ensure that $\alpha_{\rm{mm}}$ contains only variations due to actual spectral index variations and image noise, the data were re-imaged with a circular, $0.25''$-radius beam, using only baselines with lengths between $17 - 1260$ k$\lambda$. This is the largest baseline range with which it is possible to cover all the continuum spectral windows in both Band 6 and Band 7. The profile of $\alpha_{\rm{mm}}$ was subsequently derived by deprojecting and derotating the images in both bands, calculating the value of $\alpha_{\rm{mm}}$ and its standard error in each pixel, and radially averaging these values. Full maps of $\alpha_{\rm{mm}}$ and $\sigma(\alpha_{\rm{mm}})$ are included in Appendix , but here only the radially averaged profiles will be used.

The uncertainty in $\alpha_{\rm{mm}}$ at any point is determined primarily by two components: first, the noise properties of the images in both bands, and second, the absolute flux calibration errors in the images.

The image noise consists of the random noise term to $\alpha_{\rm{mm}}$ which affects each beam separately, but is also affected by residual phase noise. This is particularly relevant for the Band 7 data. Figure  shows that the effect of this type of error does propagate into the $\alpha_{\rm{mm}}$-map. However, the result is that the radial spectral index profile becomes much noisier in those low $S/N$ regions where the azimuthal asymmetries are located, which happens at distances beyond $\sim 1''(150$ AU). Therefore, we restrict the range of radii over which we discuss the spectral index variations to within this region.

In contrast to the random noise, the absolute flux calibration error adds a constant at all positions. Since our data have been self-calibrated, and the rms phase noise on the shortest baselines is less than $15\%$, assuming a $10\%$ flux calibration error is appropriate in both bands (; see also § ). This uncertainty leads to an 0.42 error on $\alpha_{\rm{mm}}$. A third source of error, specific to the radial profile of the spectral index, is the impact of our self-calibration models on the resulting image properties. It is difficult to quantify the result of slightly different source models. However, in general, their main influence should be constrained to changes in the steepness of $\alpha_{\rm{mm}}(R)$ in the bright inner disk region (the inner $38$ AU), where the self-calibration's impact is the largest, and — partly — in extra noise in the radial profile.

The profile of the spectral index $\alpha_{\rm{mm}}(R)$, shown in Figure  shows several key features: a low spectral index in the inner disk, and a peak and subsequent drop at the location of the first gap identified in § . Here we explore these features in more detail.

### A cold inner disk

In the inner beam of our observations (dominated by the bright central component with characteristic radius $38\,$AU in the model fit in § , see Figure ), the spectral index appears to be consistent with $\alpha_{\rm{mm}} = 2.09 \pm 0.03$, increasing outward. The dominant source of error in this region is the flux calibration uncertainty in both bands, resulting in an absolute uncertainty in $\alpha_{\rm{mm}}$ over the entire image of 0.42.

If we assume the emission is in the Rayleigh-Jeans regime, then $\beta \lesssim 0.5$ - a value that can only be reached if the particles have grown significantly, possibly down to gray ($\beta = 0$) opacity indicating grains with sizes much larger than the wavelengths of these observations. On the other hand, values for $\alpha_{\rm{mm}} \leq 2$ can be explained if the emission in this part of the disk is optically thick, or if the Rayleigh-Jeans approximation does not apply. In both of these cases, the intensity of the central part of the disk becomes sensitive to the temperature of the dust. The observed Band 6 peak intensity of 0.022 Jy beam$^{-1}$ in the central 0.25$''$ beam corresponds to a brightness temperature of only 12.9 K if we assume the emission is optically thick. This suggests that the disk midplane could be cold enough for CO to be frozen out within a $\sim 20$ AU radius.

To test the hypothesis of an optically thick inner disk, an independent dust temperature estimate is necessary. Previously, a temperature profile for V1094 Sco has been given by , who fit a power-law temperature structure to the disk's spectral energy distribution (SED). In their best-fit model, $T(R) = 99.3\,\rm{K} \left( \frac{R}{1\,\rm{AU}} \right)^{-0.68}$, giving a temperature of 13 K at 20 AU. An additional check of a plausible range of midplane temperatures was done with the physical-chemical disk modeling code `DALI` , which was used to perform dust radiative transfer calculations for a set of 45 disk models to obtain their dust temperatures, raytrace their emission, and image the data with the same $0.25''$ beam used for the observations. These models were based on a parametric surface density structure and dust properties with values in a range that produces a near- to mid-infrared SED similar to that of V1094 Sco. The parameter ranges and sampling for these models are given in Appendix , Table .

In both approaches, the underlying disk density structure is smooth and does not contain rings or gaps, and both therefore have radially decreasing midplane temperature profiles. In the  model, the midplane temperature at $20$ AU is $13$ K. In the grid of `DALI` models, we find midplane temperatures of $20 \pm 4$ K at this radius. The resulting raytraced images of the disk show peak brightness temperatures between 17 and 10 K and a median temperature of $12.7$ K, very close to the observed peak temperature. This shows that our observations are consistent with a cold inner disk for V1094 Sco.

An important consequence of the low midplane temperatures predicted by these models is that the midplane CO snowline, which is expected to be around 20 K for pure CO ice , is not or in the most optimistic cases barely resolved. The gaps and rings are located far outside the inner $20$ AU, meaning that another explanation for the presence of these features than the influence of a snowline is needed.

If the mm emission in this region is indeed optically thick, this implies that additional substructures could be present in the inner disk, not revealed by these data. High-angular resolution $\geq 3$ mm data are needed to probe this region, and will also allow for a large reduction in systematic uncertainty in $\alpha_{\rm{mm}}$ across the disk.

### Radial spectral index variations in the first gap/ring-pair

The high $S/N$ achieved in the image plane in both bands allows us to describe the radial variations of $\alpha_{\rm{mm}}$ out to a radius of $\sim 150$ AU. Crucially, this region contains the first gap/ring pair in the disk, allowing us to place some constraints on the grain properties in these features.

As can be seen from Figure , $\alpha_{\rm{mm}}$ increases from $2.1 \pm 0.03$ in the inner disk, by $\sim 0.9$, towards $3.0 \pm 0.2$ in the first gap at $100$ AU. Subsequently, the spectral index drops by $\sim 0.7$ to $2.3 \pm 0.3$ in the first ring at $130$ AU. Some tentative hints of a subsequent increase in $\alpha_{\rm{mm}}$ can be seen towards the location of the second gap, but are very small compared to the (large) increase in image noise at those radii, and will therefore not be further discussed.

Due to the previously mentioned absolute flux calibration error, two possible interpretations for the behavior of $\alpha_{\rm{mm}}$ between $20 - 130$ AU can be found. On the one hand, the inner region may be optically thick. In this scenario, the material in the gap is optically thin with a $\beta \sim 1$, consistent with grain growth up to mm sizes in this region , and the material in the first ring may be optically thick as well, or, interestingly, have a very low $\beta (\sim 0.3)$ due to a population dominated by grains with sizes larger than a few cm, even if the emission is not in the Rayleigh-Jeans regime.

Alternatively, the entire disk could be optically thin, so that $\beta$ becomes even higher (up to $\sim 1.4$) in the gap, requiring grains in the this part of the disk to be small ($\mu$m-sized), while larger grain sizes are found in the rings and central disk. Indeed, in simulations, local pressure maxima in the disk can reproduce both rings and gaps in disks as well as this type of spectral index profile, as shown in, amongst others, .

Intriguingly, similar radial spectral index variations have so far also been identified observationally in both TW Hya and HL Tau , suggesting that this behavior occurs more often in rings and gaps in protoplanetary disks. A decreasing $\beta$ with radius is also present in HD 163296, although these observations lack the $S/N$ to resolve the rings .

# Discussion

## Comparing multi-ringed disks

Besides studying the detailed properties of a single object, the comparison of V1094 Sco with other multi-ringed disks, as well as with other disks in the Lupus star forming region, can provide insight into the mechanisms that drive this peculiar emission morphology.

Despite the small size of the sample of multi-ringed disks found in the literature, there is quite a large spread in the temperatures and luminosities of the stars around which they are found. For example, HD 163296 has an effective temperature of $9333$ K and a luminosity of $36$ $L_{\odot}$  while V1094 Sco is both much cooler at $T_{\rm{eff}} = 4205$ K and an order of magnitude fainter with $L_{\star} = 1.95$ $L_{\odot}$ .

To facilitate the comparison between other disks in the literature, Figure  presents the normalized radial intensity profile of V1094 Sco compared to several other protoplanetary disks. The left panel compares the V1094 disk with two of the other largest and brightest Lupus disks, both around relatively cool stars: IM Lup ($T_{\rm{eff}} = 4350$ K) and Sz 98 ($T_{\rm{eff}} = 4060$ K) . The right panel compares disk around V1094 with that of TW Hya, an old, massive, multi-ringed disk around a cool star with $T_{\rm{eff}} = 3800$ K , and the previously mentioned Herbig Ae star HD 163296. All radial distances of the Lupus disks have been calculated assuming a uniform 150 pc distance (as opposed to the 200 pc distance, used in ).

Several conclusions can be drawn from Figure . First, it is clear that the continuum emission from V1094 Sco is similar in extent to that of IM Lup , with the caveat that the IM Lup observations shown here are with a larger beam width of $0.4''$, and the outer radius might thus be smaller. This disk also seems to show a continuum break, at similar radii of $\sim 100$ AU and relative intensity to V1094 Sco. Such a break is also seen in Sz 98, albeit at a higher intensity relative to the peak flux and a smaller radius ($\sim 80$ AU). Intriguingly,  found ring-like residuals for this disk when fitting it with a smooth decreasing surface density profile, implying there may be additional structure present in this object. The older TW Hya disk, in contrast, is much more compact than the Lupus disks compared here. This source also shows substructure down to much smaller scales, which would not be resolved with the beam used for the other sources, demonstrating that rings in protoplanetary disks are not necessarily limited to the most radially extended objects.

It is illustrative to compare the observed intensity profile of V1094 with that HD 163296 as derived by . While HD 163296 is more compact than V1094 Sco, the two disks share remarkably similar intensity profiles, with a central core and more extended outer region, which hosts two rings. The more compact central core of the HD 163296 disk is at least in part due to the smaller beam of the  observations ($0.2''$ at a distance of $122$ pc).

The difference in ring locations around both stars is significant, since it provides an additional, independent indication that snow lines are an unlikely explanation for the rings in V1094 Sco: in the hotter HD 163296 system, the rings are located further inward. This contrasts with the result in , who found an overlap in the (inner) ring and gap of HD163296 and the CO snowline, as well as coinciding rings and CO snowlines in the T Tauri stars in their sample (HL Tau and TW Hya). However, both of these T Tauri disks are much more compact than the V1094 Sco disk.

An alternative explanation for the formation of rings, the presence of a magnetohydrodynamic dead zone, also does not appear to be sufficient. Although the models presented by , for instance, show a ring and gap at the outer edge of the dead zone, their model does not lead to a double ring structure, such as in these disks. However, the total (gas$+$dust) column density where the  model predicts a gap and ring is similar (at $\sim 3$ g cm$^{-2}$) to that of the best-fit V1094 Sco model of  at 100 AU, indicating that one of the rings could be linked to the dead zone edge. Alternatively, the change in turbulence inside the dead zone may be an explanation for the bright inner region of the disk and the 'fried-egg' morphology of the emission. The volatile snow lines, located inside 20 AU, may contribute to the bright inner core region, as per .

A (secular) gravitational instability (GI) or the effect of planets are both plausible mechanisms for ring formation at large radii. Secular GI effects can form rings at 100 AU scales around Solar-like stars . Likewise, it is well established that one (or more) previously formed planets of sufficient mass can form gaps and rings of dust in protoplanetary disks ; however, the nature of the planets (their mass, and the number of planets per gap) is difficult to constrain based on these data alone. Observations of V1094 Sco by the *SEEDS* survey have not led to the detection of a companion ; however, the median sensitivity of those observations is between $4 - 10$ $M_J$ in the $100-300$ AU region. This means that even a quite massive planet could easily escape detection using current instruments. Alternatively, perhaps several lower-mass planets are responsible for one or both of the observed gaps.

## The occurrence rate of large dust disks

Since our survey of disks in Lupus is — with the addition of this sample from the Lupus completion survey (see § ) — complete, it is now for the first time possible to place concrete upper limits on the occurrence rate of bright, radially extended disks. Here, we look at the radii of disks based on their mm-continuum emission. Although $^{12}$CO emission is often more extended than the continuum,  show that the average ratio $R{\rm{gas}}/R_{\rm{dust}}$ is within a factor $1.5 - 3$ in Lupus. In total, there are $95$ Class-II sources in the Lupus clouds. Of these disks, IM Lup and V1094 Sco are by far the most extended with radii up to $\sim 300$ AU, assuming a 150 pc distance, with the next largest disk at a radius of $\sim160$ AU. This implies an occurrence rate of $2.1 \pm 1.5 \%$ for these largest disks.

Given the high surface densities at large radii of disks like V1094 Sco and IM Lup, and the finding of rings in V1094 Sco that could be sculpted by planet-disk interactions, it is interesting to compare the statistics of large disks to the number of massive planets at large orbital radii.  find an occurrence rate of $<2.1 \%$ for planets in orbits between $100 - 1000$ AU with masses of $5 - 13\,M_J$, around $5-300$ Myr-old stars across all spectral types.

In contrast, the number of full (non-transitional) disks in Lupus where continuum emission was detected beyond $100$ AU is significantly larger: $8 \pm 3.0 \%$ of the full population . Since we might not detect disks with fainter continuum emission at large radii, this is a lower limit to the actual number of these disks.

If we also consider the Lupus transition disk population, the number of disks with continuum emission beyond 100 AU further increases, to $18 \pm 4.3 \%$. However, these transition disks are consistent with clearing by planets interior to the gap radius, which is $<100$ AU for the entire sample . Thus, the transition disk population is not necessarily comparable to the statistics of massive exoplanets in the $100-1000$ AU interval.

The exoplanet occurrence rates discussed here are, in either case, much smaller than the number of disks with significant amounts of continuum emission at similar radii. However, the statistics of the extrasolar planet population at these orbital radii are very incomplete, especially for lower-mass planets ($<5\,M_J$) around fainter stars. Thus, it is possible that the discrepancy in the occurrence rate of planets at large orbital radii is due to the disks forming mostly lower-mass planets. Indeed, no particularly massive planets are needed to open detectable gaps at large orbital distances, since the Hill sphere radius increases linearly in orbital distance but only as the third root of mass. Both  and , for instance, find planet masses below $5\,M_J$ are sufficient to form gaps in HD 163296 and Elias 2-24 respectively. Alternatively, the large difference in planet detection statistics could be due to significant planet migration moving most of the planets that formed at large distances from their stars to smaller orbits.

Apart from comparing the population of large disks directly to the exoplanet population, relating large primordial disks to the entire Lupus transition disk sample is also informative. Interestingly, although transition disks with large cavities make up $\gtrsim 11 \pm 3.7 \%$ of the disks in Lupus , they are strongly overrepresented in the sample of largest disks: $50\%$ of disks with continuum radii in the upper octile show evidence of central cavities in their continuum emission. This led  to propose a two-mode evolutionary sequence for protoplanetary disks. In this view, large primordial disks like V1094 Sco and IM Lup are the progenitors of transition disks, while the majority of lower-mass, more compact disks may evolve through a different pathway, either without a transitional stage or as a scaled-down version of the pathway of the large disks. The massive primordial disks may trap dust by some means (possibly through planet-disk interactions, or by some other mechanism), which prevents its inward drift. This leads to the formation of rings and, later, the emergence of a large, possibly multi-ringed, transition disk with a large dust cavity.

In the context of this hypothesis, it is possible to interpret the ratio of large primordial disks to transition disks as an indication of the speed at which this evolution occurs. Since half the disks larger than $120$ AU are transition disks at the $1 - 3$ Myr age of Lupus, this process should take place on $\sim 1$ Myr timescales. Importantly, this means we have a testable prediction for the occurrence rate of large primordial disks in other star forming regions: they should be rarer in more evolved regions, and more common in younger systems, if the absolute rate of evolution across regions is the same. In terms of disk mass, too, it is then expected that the disk mass distribution in more evolved regions is dominated by a small number of transition disks, while younger regions should have an upper mass distribution with mostly large primordial disks. Although statistics across other regions at similar sensitivity and resolution are still lacking, this framework is at least consistent with the observation of several large, multi-ringed systems in very young systems: HL Tau, Elias 2-24 and AS 209 are all at most 1 Myr old . However, such young systems tend to be brighter, leading to an observational bias for high-resolution imaging campaigns towards these disks, which in turn results in their overrepresentation in the already small sample of multi-ringed disks.

# Conclusions

In this work the nature of V1094 Sco as a giant, multi-ringed protoplanetary disk is identified for the first time, based on data from the Lupus survey of protoplanetary disks. Analyzing its structure and comparing this disk to the rest of the Lupus disk population allows us to exclude some pathways for ring formation, and provides a possible link to a general disk formation picture.

-   Continuum emission from V1094 Sco extends out to 300 AU from the central star, making it larger than the other Lupus disks by a factor $\sim 5$, but similar in size to the well-studied IM Lup disk.

-   Two bright rings are identified at 130 and 220 AU, separated by gaps at 100 and 170 AU, based on a u,v-plane fit of the intensity profile.

-   Disks with continuum radii of this size ($> 200$ AU) are rare: only $\sim 2\%$ of the Class II objects in Lupus are giant disks similar to V1094 Sco.

-   Radial variations in $\alpha_{\rm{mm}}$ are consistent with optically thick emission in the bright disk core. The spectral index increases at the location of the first gap, and subsequently decreases again towards the first ring, possibly as a result of efficient grain growth in this ring.

-   Based on the temperature of the optically thick emission and on independent disk models, the midplane dust temperature at the ring locations is too low for snow lines to be a valid explanation, while a dead zone would only be expected to cause one ring. However, gravitational instabilities or the presence of planets up to $\sim 10$ $M_J$ in the gaps cannot be excluded.

-   The number of primordial disks with continuum emission at radii $> 100$ AU in Lupus seems to be much larger than the occurrence rate of giant planets with $5 - 13\,M_J$ between $100-1000$ AU. The large disk/transition disk ratio may be a probe of disk evolution, if large multi-ringed disks evolve into (multi-ringed) transition disks with large cavities on $\sim$Myr timescales.

```latex
\begin{appendix}

%\include{V1094_SED}
\renewcommand\thefigure{\thesection.\arabic{figure}}
\section{Spectral index map of V1094 Sco}
\setcounter{figure}{0}
\label{app:alphaims}

Figure~\ref{fig:alphamap} shows the full maps of $\alpha_{\rm{mm}}$ and $\sigma(\alpha_{\rm{mm}})$ that were used to construct the radial $\alpha_{\rm{mm}}$ profile discussed in Section~\ref{sec:alpha}. Pixels in both Band 6 and Band 7 with $S/N < 2.5$ were masked.

\begin{figure*}[!h]
\begin{center}
        \includegraphics[width=1\textwidth]{V1094_Sco_alphamap_noradcut2.png}
        \caption{\textit{Left}: $\alpha_{\rm{mm}}$, deprojected and derotated. \textit{Right}: Standard error on $\alpha_{\rm{mm}}$, not including the $\pm 0.42$ offset due to flux calibration errors. The Band 6 data are overplotted in black, with contours at 3, 5, 7.5, 10, 15, 20 and 30$\sigma$. The locations of the rings (R1 and R2), gaps (G1 and G2) and the bright core are indicated with white and gray arcs.}
        \label{fig:alphamap}
\end{center}
\end{figure*}

\section{Midplane temperatures of V1094 Sco disk models}
    \setcounter{figure}{0}
    \label{app:sed}
    In Table~\ref{tab:sedpar}, the parameters chosen to for our \texttt{DALI} mini-grid of V1094 Sco-like disks are listed. \texttt{DALI}~\citep{dali1,dali2} is a physical-chemical code that was used to perform dust radiative transfer calculations for the range of dust distributions, and to obtain the dust temperatures and raytrace their emission.

    The disk models themselves are based on the parametric disk density structure introduced in~\citet{andrews11}, here truncated at a radius $R_c$, so that $\Sigma(R) = \Sigma_c \left(\frac{R}{R_c}\right)^{-\gamma} \exp \left( (R/R_c)^{2-\gamma}\right)$. $R_c$ was set to 250\,AU based on the Band 6 continuum image, and we fix $\gamma$ to $0.6$. For a full disk, the key parameters are then the total mass (gas + dust) $M_d$, the scale height $h_c$ (the standard deviation of the vertical Gaussian distribution of the dust), flaring angle $\psi$ (defined so that $h(R) = h_c (R/R_c)^{\psi}$), the large grain fraction $f_{ls}$, and the relative settling factor $\chi$, which sets the large grain scale height by $h_{c\,\rm{large}} = \chi h_c$. Large grains have sizes up to $1000\,\mu$m, while small grains do not exceed $1\,\mu$m in size. The total mass was fixed at $0.1\,M_{\odot}$, reproducing the Band 6 continuum flux assuming a gas to dust ratio of 100.

\begin{table}[hb]
\caption{Parameter ranges and number of samples for the disk models}
\label{tab:sedpar}
    \centering
    \begin{tabular}{l l l l l}
    \hline \hline
    Parameter & minimum & maximum & $N_{\rm{samp}}$ \\
    \hline
    Scale height $h_c$ & 0.02 & 0.1 & 4 \\
    Flaring angle $\psi$ & 0.1 & 0.3 & 3 \\
    Large grain fraction $f_{ls}$ & 0.85 & 0.99 & 2 \\
    Settling parameter $\chi$ & 0.2 & 1.0 & 2 \\
    \hline
    \end{tabular}
\end{table}

\begin{table}[]
\caption{Stellar properties used for the disk models}
\label{tab:starpar}
    \centering
    \begin{tabular}{l l l l l l}
    \hline \hline
    %Parameter & minimum & maximum & $N_{\rm{samp}}$ \\
    $M_{\star}$\tablefootmark{a} & $R_{\star}$   & $L_{star}$    & $T_{\rm{eff}}$ & SpT & $d$ \\
    $[M_{\odot}]$                & $[R_{\odot}]$ & $[L_{\odot}]$ & $[\rm{K}]$     &     & pc  \\
    \hline
    0.92          & 1.9           & 1.7           & 4205           & K6  & 150 \\
    \hline
    \end{tabular}
    \tablefoot{
    Values derived from~\citet{alcala17} assuming a 150\,pc distance to the disk.\\
    \tablefoottext{a}{Mass derived with the \citet{siess00} model grid.}
    }
\end{table}

\begin{figure}[!htb]
\begin{center}
        \includegraphics[width=0.5\textwidth]{DALI_midplane_temps_convolved_hist.png}
        \caption{\texttt{DALI} model midplane temperatures at 20\,AU (orange) and Band 6 peak brightness temperatures (blue) after convolution with an $0.25''$ beam for a grid of disks. The black line indicates the observed peak brightness temperature of V1094 Sco. For pure CO ice in typical disk midplane conditions, the freeze-out temperature is 20 K ({\"O}berg et al. 2005).}
        \label{fig:sedpar}
\end{center}
\end{figure}

\end{appendix}
```