author: Javier Garcia-Bernardo; Jan Fichtner; Frank W. Takes; Eelke M. Heemskerk
bibliography: Zotero.bib
title: Uncovering Offshore Financial Centers: Conduits and Sinks in the Global Corporate Ownership Network

# Introduction

Multinational corporations use highly complex corporate structures of parents and subsidiaries to organize their global operations and ownership structure. For example, the Britain-based banking and financial services company *HSBC* is composed of at least 828 legal corporate entities in 71 countries. The largest brewing company in the world, *Anheuser-Busch InBev*, consists of at least 680 corporate entities involving 60 countries. These complex corporate structures purposefully span across countries and jurisdictions in order to increase competitive advantage by minimizing costs and accountability . Some jurisdictions are particularly popular, such as Bermuda, the British Virgin Islands, and the Cayman Islands . Often referred to as *Offshore Financial Centers* (OFCs), these jurisdictions attract financial activities from abroad through low taxation and lenient regulation. With an estimated 50% of the world's cross-border assets and liabilities (US\$21-US\$32 trillion) passing through OFCs  they have become dominant nodes in the transnational financial-economic network in which capital is stored and redistributed. The largest OFCs generally have well-developed regulatory institutions and comply with international laws on trade and money laundering . At the same time, the services that OFCs offer are increasingly under scrutiny for, amongst others, facilitating corporate tax avoidance. We consider that OFCs are not only places where capital is 'stored', but act as nodes in a complex network of international capital flows. We suggest a novel approach for identifying and classifying offshore financial centers based on the underlying large-scale granular firm level ownership data.

Corporations create complex corporate ownership structures for at least three reasons. First, corporations seek to increase legal protection. By organizing parts of their corporate structure in certain trusted territories with favorable legal conditions they can increase legal certainty for their operations or for joint-ventures. And by setting up subsidiaries in specific jurisdictions and using such subsidiaries to invest in other countries, multinationals can hedge their investment against decisions of governments. Second, favorable regulatory regimes in OFCs can be used by companies to avoid corporate accountability and public scrutiny of their operations, i.e. regulatory arbitrage. For instance, many of the opaque structured financial products that aggravated the global financial crisis since 2008 were created in OFCs . Third, complex corporate ownership structures help to minimize tax payments – especially for corporations that have many intangible assets, such as intellectual property rights. For instance, between 2007 and 2009 Google moved the majority of its profits generated outside the United States (US\$12.5 billion) to Bermuda through corporate entities in the Netherlands. As a result Google paid an effective tax rate of 2.4% on all operations . Similarly, Apple used Ireland to avoid US\$14.5 billion in taxes since 2003 , and Starbucks UK voluntarily payed £20 million after it came to light that it had paid virtually no taxes since establishing in the UK . In total, every year multinationals avoid paying US\$50-200 billion in taxes in the European Union using OFCs . In the United States, tax evasion by multinational corporations via offshore jurisdictions is estimated to be at least US\$130 billion per year. While reducing costs is in general a valid concern for corporations, these practices significantly diminish tax revenues and as such may inhibit the capacity of states. Moreover, tax avoidance provides a competitive advantage to multinationals over small and middle sized companies, which are taxed at national rates – ranging between 25 and 35% for most developed countries.

Given the policy-concerns caused by the use of OFCs in recent years, the identification of OFC jurisdictions has become a politicized and contested issue. The International Monetary Fund (IMF) and the Organization for Economic Co-operation and Development (OECD) have published lists of OFCs based on a qualitative assessment of the jurisdictional regulations and taxation frameworks. However, this approach says little about the real-world relevance of specific jurisdictions as OFCs and is also vulnerable to political influence. In contrast to these qualitative approaches, Zoromé has defined OFCs as jurisdictions that 'provide financial services to nonresidents on a scale that is incommensurate with the size and the financing of their domestic economies' . Zoromé  as well as Cobham et al.  used flow data on the export of financial services to calculate ratios that indicate how significantly a jurisdiction acts as an OFC. Fichtner further expanded this approach by using stock data on international banking assets, portfolio investment, and foreign direct investment (FDI) in relation to the gross domestic product (GDP) of a jurisdiction to calculate an 'offshore-intensity ratio' . While the offshore-intensity ratio approach identifies which countries have a disproportionate value of inward foreign investment, it is limited in at least two ways. First, aggregated national statistics such as FDI are subject to political preferences and influence . For instance, inward FDI is systematically underreported . Second, the existing methods cannot shed light on the position of a jurisdiction in the broader network of capital flows since they are not able to differentiate if the inward foreign investment reported by Bermuda originates in the Netherlands, or if in contrast it originates in Germany and is routed through the Netherlands.

We propose a novel network analytic approach to identify OFCs based on a country's position in the *network of global corporate ownership*. In contrast to prior work, our approach characterizes countries based on a distinction between sink-OFCs and conduit-OFCs. *Sink-OFCs* are countries that attract and retain foreign capital – territories in this category are usually characterized as tax havens, such as the British Virgin Islands, the Cayman Islands and Bermuda. Most sink-OFCs have small domestic economies and large values of foreign assets, which are attracted through low or zero corporate taxes. Because of this disparity between the local economy and external assets, the aforementioned offshore-intensity ratio approach is well suited for identifying these sink-OFCs . *Conduit-OFCs* on the other hand are 'countries that are widely perceived as attractive intermediate destinations in the routing of investments' . Conduit-OFCs typically have low or zero taxes imposed on the transfer of capital to other countries, either via interest payments, royalties, dividends or profit repatriation. In addition, such jurisdictions have highly developed legal systems that are able to cater to the needs of multinational corporations. Conduits play a key role in the global corporate ownership network by allowing the transfer of capital without taxation. In this way, profit from one country can be re-invested in another part of the world paying no or little taxes. Countries such as the Netherlands and Ireland have been criticized for these types of activities .

The building blocks of our method for identifying OFCs are what we call *global ownership chains* (GOCs), in which a series of companies are connected in a chain if for each two directly subsequent entities $A$ and $B$, it holds that firm $A$ is owned by firm $B$, i.e., there is a link between them in the ownership network. Under European Union regulations (Council Directive 2003/123/EC), transfers of capital without taxation are typically only allowed through ownership links from subsidiaries to parents, meaning that value can flow from $A$ to $B$. These EU regulations are expanded to other countries via tax treaties, allowing companies to transfer capital outside the EU through corporate structures. Based on the value going through these international ownership chains, we propose two new centrality measures specifically aimed at measuring the extent to which a jurisdiction is a sink-OFC or conduit-OFC. We furthermore introduce an entropy-based metric that can characterize the specialization of an OFC in terms of which countries it services.

The proposed network analytic approach to identifying OFCs has a number of advantages. First, it makes no *a priori* assumptions about the global economy and the countries involved; the possible identification of a country as an OFC is purely *data-driven*. Second, it does not rely solely on aggregated macroeconomic indicators that may introduce significant noise and deviations, but on *fine-grained data* of firm-level corporate ownership. Third, this firm-level data allows us for the first time to quantitatively identify and distinguish between both sink-OFCs and conduit-OFCs. Fourth, our approach is also suitable to classify and characterize the *specialization* of OFCs across geographic regions and industrial sectors. OFC specialization is a key issue from a policy and regulatory perspective, but cannot be answered by existing quantitative methods. If OFCs are not specialized and occupy structural and functional equivalent positions in the network of global corporate ownership, we can assume that firms can easily re-organize their corporate structures to other OFCs in the wake of regulation . On the other hand, if OFCs are specialized then there is likely more room for tailored policy interventions and regulation. Our research demonstrates that – contrary to the still prevalent conjecture that offshore finance resembles an 'atomized' marketplace in which a multitude of approximately equal jurisdictions compete and where regulation is therefore unfeasible – the corporate use of OFCs is in fact concentrated in a small number of key jurisdictions, most of which are highly-developed OECD countries.

# Methods

## Data extraction and quality assessment

We sourced company ownership data from the Orbis database (<http://orbis.bvdinfo.com>) in November 2015. Orbis is a unique and frequently used information provider that covers about 200 million public and private firms worldwide  compiled from official country registrars and other country collection agencies. For each available company, we extracted its operating revenue, country, city, sector, global ultimate owner (the parent firm who owns at least 50% of the company directly or indirectly and is not itself owned by any other firm) and all ownership relationships, with the direct and total ownership percentage. Moreover, since companies located in Isle of Man (IM), Jersey (JE) and Guernsey (GG) use the country code of the United Kingdom, every company located in IM, JE or GG cities was given the country code of the territories (see Supplementary Methods). The resulting dataset contains 71,201,304 distinct ownership relationships between 98,255,206 companies. Note that our data selection method does not include private wealth, only corporate structures are considered.

Data quality of Orbis differs across regions. Firm coverage is better for high-income countries than for low-income countries . In the United States, a significant number of companies registered in the state of Delaware are not covered because they are not required to file information – a problem shared between all corporate information providers. In general, poor data quality is associated with fiscal secrecy since these jurisdictions do not consistently report the companies registered in their territories. So, if we are missing data, then this data is more likely to be missing in an OFC than in a non-OFC. Therefore, the findings that we report likely represent a *lower bound* on the position of offshore financial centers in the corporate ownership network.

Two extra steps were done to ensure a correct analysis of the data. First, we deconsolidated financial accounts. Generally, only consolidated financial information is available for large companies – i.e., the revenue of all the subsidiaries is reported in both the subsidiaries and the parent company. We corrected this by recursively subtracting the operating revenue of all subsidiaries of companies with consolidated accounts along all paths of ownership relations (see Supplementary Methods). Second, since the information is collected by different country-level agencies and merged by Orbis, the sum of the ownership stakes was normalized as in Vitali et al.  to account for missing shareholders (see Supplementary Methods).

## From company data to global corporate ownership chains

In this section we outline our approach that from the ownership data constructs country chains which can ultimately be used to detect OFCs. For theoretical definitions of the the different concepts in each of these construction steps, the reader is referred to the Supplementary Methods.

### Ownership network

We considered the network of ownership relations as a directed graph with the firms as nodes and the links as ownership relations, where value flows from corporate entities to their owners. The structure of this network itself has been extensively studied and exhibits common properties of complex networks such as a power-law degree distribution and weight (strength) distribution. We can furthermore observe the emergence of a giant weakly connected component capturing the majority of corporations in the network, as well as a bow-tie structure featuring a smaller strongly connected component in the center. The remainder of this paper considers the derivation of chains from the ownership links, and is not directly concerned with further exploration of the network's macro level structure. As such, for more details the reader is referred to the excellent study of the global ownership network by Vitali et al. .

### Company level chains

From the set of over 71 million ownership relations we identified the global corporate ownership chains (GCOCs) as follows. For each node, we applied a depth-first search algorithm, exploring the resulting network as far as possible along each branch before backtracking, forming chains from the starting node. We continued adding nodes to a chain until the multiplicative ownership fell below 0.001 (for instance, four companies in a chain owning 10% of the next). Chains reaching the origin node (i.e., loops) or a node previously visited in the considered chain were ignored to avoid infinite loops. The results are robust to variations of the multiplicative ownership threshold (Supplementary Information). Our approach still reflects country round-tripping (where value flows from country A to B and to A again) since this strategy requires two different companies in country A. We repeated this process for all nodes in the network, which resulted in a set of 11,404,819 ownership chains.

### Multiplicative ownership

For each ownership chain we determined its value using the weighting method by Vitali et al. . It weighs the value $V_p$ of a chain $p$ in terms of the revenue of the initial company in the chain: $V_{C_1|C_2|C_3} = R_{C_1} \cdot \mathit{MO}_{C_1|C_2|C_3}$. Here $C_1|C_2|C_3$ corresponds to a chain of three companies in which $C_2$ owns $C_1$ and $C_3$ owns $C_2$, $R_{C_1}$ is the operating revenue of company $C_1$ and $\mathit{MO}$ is the multiplicative ownership, i.e., the product of the weights of the links between the subsequent firms in the chain.

### Aggregating at the country level

Next we determined the country of domicile for each corporate node in a GCOC. Given our goal of studying transnational links, we merged together adjacent nodes in a chain that are located in the same country. Finally, we divided each chain into *chunks* of length ${2,3, \ldots, |\textrm{chain}|}$, which resulted into 108,159,506 chunks. In order to avoid double counting revenue, we kept the maximum revenue for each group of chunks matching to the same country level chain, the same owner and the same upstream companies. For instance, we can have a chain $A_1|A_2|B_1|B_2$, where companies $A_1$ and $A_2$ are located in country $A$, and companies $B_1$ and $B_2$ are in country $B$. From this chain, we can create all chunks in Figure A, matching to the same country chain ($A|B$). Since the value originates in the same company ($A_1$), it has the same owner in country $B$ ($B_1$), and the flow happens through the same companies ($A_1|A_2$), we only keep the chunk with the largest value. In this way we obtained 16,448,469 chains.

After grouping all chains going through the same countries we obtained 377,098 different country-level ownership chains (e.g. $ES|NL|LU$). In these chains value flows from a company in the source (ES) to an owner in the sink (LU) (the sink owns the source either directly or through conduit countries). Next, we identify OFCs based on particular motifs in the ownership chains.

## Finding Sink and Conduit Offshore Financial Centers

### Avoiding double counting of sinks and conduits

In order to identify sinks and conduits we focus on two subsets of country chains. The first subset contains 52,655 chains of size three and is used for the conduit-OFC analysis. This selection strategy ensures that conduits (countries in the middle of an ownership chain) are not double counted (Fig. B–C). For instance, using chains of length four would result in double counting country C in Fig. C. The second subset contains 7,172 chains of size two and is used for the sink-OFC analysis. This selection strategy ensures that sources and sinks are not counted twice.

### Sink Offshore Financial Centers (sink-OFCs)

Sink-OFCs are jurisdictions that attract and retain foreign capital, i.e., jurisdictions in which GCOCs end. Thus, we define the sink centrality ($S_c$) of a country as the difference of value entering and leaving the country, divided by the sum of all value in the network. Moreover, since this difference is proportional to the size of the country, we normalized the centrality by the gross domestic product (GDP) of the country.

$$\label{eq:ss}
S_c =
\frac{\smashoperator[r]{\sum_{g \in G^2:\ g[1] = c}}V_{g} -\smashoperator[r]{\sum_{g \in G^2: g[0] = c}}V_{g}}
{\smashoperator[r]{\sum_{g \in G^2}}V_{g}}
\hspace{2mm}
\cdot
\hspace{2mm}
\frac{\sum_i{GDP_i}}
{GDP_c},$$

Here, $G^2$ is the set of country chains of size two and $g[i] = c$ means that the $i$-th country in country chain $g$ is country $c$. Furthermore, $GDP_c$ is the GDP of country $c$. Using the measure in Equation , we define as sink offshore financial centers (sink-OFCs) those countries that have a disproportional amount of value staying in the country, where a disproportional amount is set to 10 – i.e., the value staying in the country is 10 times higher than the value that would correspond to the country in terms of its GDP. The results are robust to variations of this threshold (Supplementary Information).

### Conduit Offshore Financial Centers (conduit-OFCs)

We define conduits as *jurisdictions that act as intermediate destinations to sink-OFCs*. The conduit centrality $C_c$ of a country $c$ is defined in two axes. The first axis, inward conduit centrality ($C_{c_{in}}$), measures the value of chains flowing from a sink-OFC, into the conduit, out to another country. The second axis, outward conduit centrality $C_{c_{out}}$ measures the value of chains flowing from any country, into the conduit, out to a sink-OFC. Since this flow is proportional to the size of the country, we normalized the centrality by the gross domestic product (GDP) of the country.

$$C_{c_{in}} =
\frac{\smashoperator[r]{\sum_{g \in G^3_{s1}:\ g[2] = c}}V_{g}} {\smashoperator[r]{\sum_{g \in G^3}}{V_{g}}}
\hspace{2mm}
\cdot
\hspace{2mm}\frac{\sum_i{GDP_i}}
{GDP_c},$$ $$C_{c_{out}} =
\frac{\smashoperator[r]{\sum_{g \in G^3_{s3}:\ g[2] = c}} V_{g}}{\smashoperator[r]{\sum_{g \in G^3}}{V_{g}}}
\hspace{2mm}
\cdot
\hspace{2mm}\frac{\sum_i{GDP_i}}
{GDP_c},$$

Here, $G^3$ are chains of length three, $G^3_{s1}$ is the subset of $G^3$ in which the first country in the chain is a sink-OFC. Analogously, $G^3_{s3}$ is the subset of $G^3$ in which the third and last country in the chain is a sink-OFC. Countries with a $C_c$ larger than 1 in both axes ($C_{c_{in}}$ and $C_{c_{out}}$) were considered conduit-OFCs. We empirically found that the division between conduit-OFCs and other countries occur naturally around $C_{c(in/out)} = 1$ (Supplementary Information).

In order to understand if a country $c$ serves as a conduit for many other countries or only a few, we calculate the Shannon entropy of the distribution of values leaving ($EC_{c_{out}}$) or entering ($EC_{c_{in}}$) the country in chains of size three. Shannon entropy is lower for conduits catering to few countries and more skewed value distributions – i.e., most of the value coming from the same country; and higher for conduits catering to many countries and with more even distributions. $$EC_{c_{in}} =-\sum_{i}
{(V_{i|c|s}\log{(V_{i|c|s})})},$$ $$EC_{c_{out}} =-\sum_{i}
{(V_{s|c|i}\log{(V_{s|c|i})})},$$ Where $V_{i|c|w}$ denotes the sum of value in all $G^3$ where the source is country $i$, the conduit is country $c$ and $s$ is any sink-OFC, normalized such that $\sum_i{(V_{i|c|s})} = 1$.

# Results

## Identifying Sink Offshore Financial Centers

We first identified sink and conduit Offshore Financial Centers. Figure A shows the sink-OFC centrality for every country, with the absolute value of the metric on the horizontal axis and the proper GDP-normalized value on the vertical axis (see Methods). Node size reflects the sum of all value entering or leaving the country. Orange colored countries receive more value than they send and are thus net sinks. Green countries cannot be sink-OFCs because they are net senders.

We denoted the 24 orange colored countries above the factor 10 threshold in Fig. A as sink OFCs (Tables  and S2). The sink-OFC centrality is proportional to FDI and as such mirrors the offshore intensity ratio approach  for identifying sink-OFCs (Supplementary Fig. S1). Our list of sink-OFCs is indeed associated with territories with low or zero corporate taxes, where capital accumulates. While Panama is well known as a tax haven, our approach does not identify it as a sink-OFC. Panama is mainly a tax haven for individuals and with relatively high corporate taxes (25%) is less attractive for corporate groups. In the notorious Panama Papers (<https://panamapapers.icij.org>) the large majority of the involved shell companies were actually domiciled in the British Virgin Islands (VG).

| ISO2 | Country name | $S_c$ | ISO2 | Country name | $S_c$ | ISO2 | Country name | $S_c$ |
|:---|:---|:---|:---|:---|:---|:---|:---|:---|
| VG | British Virgin Islands | 5235 | MH | Marshall Islands | 100 | BZ | Belize | 38 |
| TW | Taiwan | 2278 | MT | Malta | 100 | GI | Gibraltar | 34 |
| JE | Jersey | 397 | MU | Mauritius | 75 | AI | Anguilla | 27 |
| BM | Bermuda | 374 | LU | Luxembourg | 71 | LR | Liberia | 17 |
| KY | Cayman Islands | 331 | NR | Nauru | 67 | VC | St. Vincent & Granadines | 14 |
| WS | Samoa | 277 | CY | Cyprus | 62 | GY | Guyana | 14 |
| LI | Lichtenstein | 225 | SC | Seychelles | 60 | HK | Hong Kong | 14 |
| CW | Curaçao | 115 | BS | Bahamas | 40 | MC | Monaco | 11 |

List of sink-OFCs, ordered by sink centrality value.

In contrast to the results from the offshore-intensity ratio approach, our method identifies Taiwan (TW) as a very prominent sink-OFC. In-depth studies by tax specialists have suggested that Taiwan is an 'unnoticed tax haven' , since it has not signed the OECD Common Reporting Standard for the automatic exchange of financial information and maintains bank confidentiality. The prominence of Taiwan is driven by Taiwanese technological companies, which often own Chinese firms through Hong Kong (33%) and Caribbean Islands (20%), or own Hong Kong firms through Caribbean Islands (12%). Due to pressure by China, Taiwan does not participate in FDI statistics collected by the IMF and therefore was not detected by studies relying on aggregated international FDI data.

Our approach replicates the outcomes of the offshore-intensity ratio method, but is also able to identify hitherto 'unnoticed' sink-OFCs. In addition, our approach enables us to quantify the weight of each territory in corporate offshore finance. In conclusion, the five largest sink-OFCs in terms of non-normalized sink-OFC centrality are: Luxembourg, Hong Kong, the British Virgin Islands, Bermuda and Jersey (Supplementary Table S2).

## Identifying Conduit Offshore Financial Centers

While sink-OFCs store capital, conduit-OFCs facilitate the movement of capital between sink-OFCs and other countries. Figure B shows the inward versus outward conduit-OFC centrality (see Methods). We marked in red countries identified as sink-OFCs in Fig. A. Blue countries are not OFCs since they are not sink-OFCs and only a moderate sum of value in chains ending (or starting) in sink-OFCs goes through them. We denote all the jurisdictions in the upper right quadrant as conduit OFCs. As expected, the sink-OFCs (in red) are also important conduits.

Green countries are distinct conduit-OFCs (Table ). Our approach identifies five large conduit-OFCs that do not act as sink-OFCs: The Netherlands, the United Kingdom, Switzerland, Singapore and Ireland. These countries facilitate the transfer of value from and to sink-OFCs, and are used by a wide range of countries (Figure C). Importantly, these countries are also used extensively as conduits to non-OFCs (Figure D), indicating that conduit-OFCs are not used exclusively for the transfer of value to sink-OFCs. This contrasts with countries such as Russia, China or most sink-OFCs, which are used by companies as conduits to sink-OFCs more frequently than as conduits to non-OFCs (Figure D). Finally, our approach also identifies three small conduit-OFCs that do not act as sink-OFCs: Belgium, Panama and Guernsey. For a comparison of our approach with other rankings of OFCs, see Supplementary Information and Supplementary Table S5.

| ISO2 | Country name | Non-normalized $C_{c_{out}}$ | Non-normalized $C_{c_{in}}$ | $C_{c_{out}}$ | $C_{c_{in}}$ |
|:---|:---|:---|:---|:---|:---|
| NL | The Netherlands | $7.4\cdot10^{11}$ | $3.8\cdot10^{11}$ | 18.6 | 22.5 |
| GB | United Kingdom | $3.8\cdot10^{11}$ | $1.3\cdot10^{11}$ | 3.1 | 2.4 |
| CH | Switzerland | $2.2\cdot10^{11}$ | $2.7\cdot10^{10}$ | 6.9 | 2.0 |
| SG | Singapore | $7.2\cdot10^{10}$ | $2.2\cdot10^{10}$ | 5.1 | 3.8 |
| IE | Ireland | $6.4\cdot10^{10}$ | $3.3\cdot10^{10}$ | 5.9 | 7.2 |

List of conduit-OFCs, ordered by value flowing through the conduit toward sink-OFCs.

## Geographical specialization

Next, we investigated if there exists geographical specialization in the ownership network. Figure  shows the network of value flow between countries, where the color indicates the relative importance of the link in relation to a simple null model where the weights are set to the product of the GDP of each pair of countries (see Supplementary Methods). The size of the nodes is proportional to the fraction of chains where that country appears, and the color shows the sink-OFC centrality. The position of a country in the network is set through a force-directed layout, where well connected countries are close in space. This is reflected by European countries placed close to the Netherlands (NL) and Luxembourg (LU), while Asian countries are placed close to Hong Kong (HK) and other sink-OFCs, and United Kingdom (GB) acts as an integrator between Europe and Asia.

The geographical specialization observed at the country-by-country ownership network (Fig. ) can be further dissected at the global corporate ownership chain level. Figure  shows the countries that appear in the source and sink position of the chain. Since we are interested in conduit-OFCs, we restricted the analysis to those chains ending in a sink-OFC. For each country in the horizontal axis we visualized two columns. The leftmost column represents the value associated to the source countries, while the rightmost column shows the distribution of the value associated to sink-OFCs. The Netherlands and the United Kingdom are the largest conduits, with values going through them twice the value of the next largest country (Luxembourg). Each country is specialized in a geographical area: the United Kingdom serves as a conduit between European countries and LU, BM, JE, VG and KY. The Netherlands is the principal conduit between European companies and LU, CW, CY and BM. Importantly, the majority of investments from LU or HK do not require a conduit, and thus LU and HK companies invest directly in European countries and China (Fig.  inset). In contrast, investments from countries typically identified as tax havens (e.g., BM, VG or KY) do, and thus companies located in these jurisdictions invest in other OFCs (Fig.  inset).

Other countries are also specialized. Switzerland is used as a conduit to Jersey. Ireland is the route for Japanese and American companies to Luxembourg (note that since the US is not considered a sink-OFC, companies with the structure $EUcountry \rightarrow IE \rightarrow US$ (such as Apple) do not contribute to the conduit-OFC centrality of Ireland). Cyprus is primarily used by Russian companies owned from the British Virgin Islands. Belgium is used as a conduit for the company Euroclear (see Supplementary Information).

British Virgin Islands (VG), Cayman Islands (KY) and Bermuda (BM) are strongly interlinked, with many chains starting and ending in the three jurisdictions under British sovereignty (as evidenced in Fig. ). However they are used by different countries. VG is used for chains starting in China (CN) and Hong Kong (HK) and ending in Taiwan (TW) or HK. BM is used for European companies owned from the VG or LU. Cayman Islands (KY) chains start in TW, HK and China and end in TW and HK. Finally, HK and Singapore (SG) are small territories connecting Asia, Europe and sink-OFCs: HK acts as the conduit between China and TW, VG, KY; SG has an important role joining many different countries (mainly South-East Asia, CN, GB and NL) to HK, VG, KY and BM.

## Sector specialization

Next, we looked at whether specialization was also present at the sectoral level. We first converted each global corporate ownership chain in terms of the sectors involved by using the statistical classification of economic activities in the European community (NACE Rev. 2). For instance the chain $291\rightarrow 642\rightarrow 829$ indicates that a company in sector 291 (manufacture of motor vehicles) is owned by a company in sector 642 (holding company), which in turn is owned by a company in sector 829 (administrative and support activities sector). In order to explore the sectoral specialization, we classified the sectors by their dominant position in chains of size three (in terms of value), finding six categories: only source, only conduit, only sink, source+conduit, source+sink, conduit+sink and source+conduit+sink (see Supplementary Methods). For instance, sectors prevalent in the source position were assigned to the source group. Figure A summarizes our findings. We found that while the source category is relatively diversified in terms of sectors, this is not the case for conduit and sink sectors. The conduit positions in the chains are dominated by 'Holding companies' (prominent in the Netherlands). Sinks are specialized by sector as well, with 'Administrative' (829, prominent in Luxembourg) and 'Unknown' (prominent in sink-OFCs) sectors attracting the majority of the value. Two sectors are prominent in both source and conduit position, 'Head offices' (701) and 'Wholesale trade' (46) sectors.

The concentration is more pronounced at the end of the chains, where administrative and office support (829, driven by Luxembourg), oil industry (08/09, driven by Russia) and manufacture of computers (26, driven by Hong Kong and Taiwan) are the main sectors. As expected since sink-OFCs are usually associated with secrecy, companies with missing sectors (00, driven by sink-OFCs) also appear at the end of the chains. Finally, sectors that are evenly distributed among sectors correspond to other types of financial companies (sector 64). Different sectors appear in different parts of the chains, with only a small number of sectors appearing in the conduit and sink positions.

## Sectoral-geographical specialization

Finally, we investigated if the geographical specialization in conduits correlates with a sectoral specialization. Since conduit-OFCs have favorable legislation for the establishment of only a few concrete types of entities (e.g., head offices, or holding companies), we expected that the majority of the value going through a conduit-OFC would be concentrated in one or a few sectors. Figure B shows the sectoral specialization for different conduit-OFCs as well as Luxembourg and Hong Kong, given their relevance as sink-OFCs and their position in the European Union and in proximity of China. For each country in the horizontal axis we show three columns. The leftmost column represent the distribution of sectors of the companies located in the source country, the middle columns represent the distribution of sectors in the conduit country, while the right columns represent the distribution of sectors in the companies within the sink-OFCs. In general, we observe that while the source sector is generally known, this is not the case for the sink sector. Consequently with their classification as sink-OFCs, 40% of LU and 70% of HK companies involved in chains ending in sink-OFCs have missing sector information. We found that the Netherlands (NL), the United Kingdom (GB), Switzerland (CH) and Singapore (SG) are specialized in holding companies (code 642), while Ireland (IE) (together with GB) specializes in activities of head offices (code 701).

For chains using Luxembourg as a conduit, 90% of the value ends up in an unknown sector in a sink-OFC. Similar trends are observed for HK, SG and IE. For Switzerland, 70% of the value ends in sector 089 in Jersey. We observe the $641|642|641$ sandwich in GB, corresponding mainly to HSBC. Manufacturers of electronics and computers (261 and 262) use predominantly HK as conduit, while holding companies use the Netherlands and Luxembourg as conduits, and other financial sectors use Ireland. Finally, we observe that Dutch companies are owned by companies in sector 829 (concentrated in Luxembourg). Hence, we observe a clear sectoral specialization of conduit-OFCs.

# Discussion

Prior work on OFC identification used either qualitative assessments of policies and regulations or a quantitative approach based on ratios of foreign investment to GDP. Here we develop a novel method for OFC identification by analyzing the transnational ownership network based on global corporate ownership chains. By investigating which countries are used predominantly as owners in the chain we are able to find a set of 24 sink-OFCs. We show that the majority of investment from and to sink-OFCs occurs through only five jurisdictions (conduit-OFCs): the Netherlands, the United Kingdom, Switzerland, Ireland and Singapore. The Netherlands and the United Kingdom rank among the largest countries in global cross-border direct investment, according to IMF data. Thus, it could have been expected that these two countries would also appear in a firm-level analysis of the global corporate ownership network, such as ours. However, our method for the first time uncovers the role of both countries as dominant conduits within transnational corporate ownership chains. This granular analysis sheds new light on the outsized role of the Netherlands and the United Kingdom in global finance.

We find a clear geographical specialization in the offshore financial network : the Netherlands is the conduit between European companies and Luxembourg. The United Kingdom is the conduit between European countries and former members of the British Empire, such as Hong Kong, Jersey, Guernsey or Bermuda. Hong Kong and Luxembourg, being themselves sink-OFCs, also serve as the main countries in the route to typical tax havens (British Virgin Islands, Cayman Islands and Bermuda). The specialization is not only geographical, but also present at the sector level. For instance, the Netherlands is specialized in holding companies, the United Kingdom in head offices and fund management, Ireland is prominent in financial leasing and head offices, Luxembourg in support activities, Hong Kong and Switzerland are characterized under the "Other financial activities category", which encompasses commodity dealing, financial intermediation and derivatives dealing. Companies choose to centralize their investment apparatus in specific jurisdictions according to the tax regulations of the jurisdiction, its bilateral tax treaties, and its sectoral advantages. Our network approach thus sheds light on the geographic structure of the global ownership network, finding that only a small set of territories act as sinks of ownership chains (most of them under British sovereignty), and even a smaller subset act as conduits for ownership structures to sinks.

Our approach identifies, characterizes and ranks OFCs and as such helps to increase transparency of and insight in highly complex international corporate financial structures. The developed method to identify OFCs improves previous attempts such as the list of tax havens published by the European Union in 2015, where countries such as Luxembourg or the Netherlands – the most prominent sink-OFC and conduit-OFC (see Supplementary Table S5) – were not included. The European Union list also does not rank jurisdictions, giving the same status to the British Virgin Islands and to Anguilla, while in fact 170 times more value ends in the former than in the latter.

Since the financial crisis, the G20 and the OECD have increased pressure on tax evasion . However the effects of these efforts have been modest . Our contribution can help regulators target the policy to the sectors and territories where the offshore activity concentrates. While efforts usually focus on small exotic islands, we showed that the main sinks of corporate ownership chains are highly developed countries which have signed numerous tax treaty agreements (Luxembourg and Hong Kong). Moreover, we showed that only a small number of conduits canalize the majority of investments to typical sink-OFCs such as Bermuda, British Virgin Islands, Cayman Islands, or Jersey (which, in fact, are all under the sovereignty of the United Kingdom). Targeting conduit-OFCs could prove more effective than targeting sink-OFCs, since – while new territories with low or no corporate taxes are continuously emerging – the conditions for conduit-OFCs (numerous tax treaties, strong legal systems, good reputation) can only be found in a few countries.

Future work could investigate the resilience of the network, in order to find points of action for legislators. For instance, an estimation of the changes in corporate structures at the global scale following the implementation of a policy could be calculated by co-analyzing tax treaties, trade networks, and corporate ownership chains (in the same vein as Weyzig ). Future work can also compare the structure and resilience of the ownership network with other networks such as the airline  or trade networks , finding how trade (representing 'real' economy), ownership (representing flux of value) and airlines (representing connections between countries) are interlinked.

# Acknowledgements

This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement number 638946). We thank Arjan Reurink for useful suggestions and comments.

# Author contributions statement

All authors were involved in the research design. JGB collected the data, implemented the methods and performed the experiments. All authors wrote and reviewed the manuscript.

# Additional information

**Competing financial interests** The authors declare no competing financial interests.


abstract: We investigate experimentally magnetic frustration effects in thermally active artificial kagome spin ice. Starting from a paramagnetic state, the system is cooled down below the Curie temperature of the constituent material. The resulting magnetic configurations show that our arrays are locally brought into the so-called spin ice 2 phase, predicted by at-equilibrium Monte Carlo simulations and characterized by a magnetic charge crystal embedded in a disordered kagome spin lattice. However, by studying our arrays on a larger scale, we find unambiguous signature of an out-of-equilibrium physics. Comparing our findings with numerical simulations, we interpret the efficiency of our thermalization procedure in terms of kinetic pathways that the system follows upon cooling and which drive the arrays into degenerate low-energy manifolds that are hardly accessible otherwise.
author: I. A. Chioar$^{1,2}$, B. Canals$^{1,2}$, D. Lacour$^3$, M. Hehn$^3$, B. Santos Burgos$^4$, T. O. Menteş$^4$, A. Locatelli$^4$, F. Montaigne$^3$, and N. Rougemaille$^{1,2}$
date: 2024-09-30
title: Kinetic pathways to the magnetic charge crystal in artificial dipolar spin ice

Artificial spin ice (ASI) are systems designed to explore the intriguing physics observed in magnetically frustrated materials. Generally fabricated by using lithography techniques, they offer almost infinite possibilities to construct a wide variety of spin models which can be accessed experimentally in a controlled manner . ASI systems have been the subject of intense research in the last few years and have allowed the investigation of a rich physics and fascinating phenomena, such as the exploration of the extensively degenerate ground state manifolds of spin ice systems , the evidence of new magnetic phases in purely two-dimensional lattices and the observation of pseudo-excitations involving classical analogues of magnetic monopoles . Notably, artificial spin ices comprise very different types of systems, including macroscopic arrays of compass needles , Josephson junctions , superconducting loops , optical traps and colloidal systems .

Up until recently, most of the experimental realizations based on magnetic nanostructures were considered as purely athermal systems. Therefore, demagnetization protocols based on the slow decay of an oscillating field are often used to drag such systems into disordered magnetic phases . However, these protocols show severe limitations in bringing ASI into their predicted low-energy magnetic configurations, where exotic effects emerge, and several other routes have been suggested to make square or kagome ASI thermally active . Up to now, two main directions have been proposed. The first one consists in working close to the blocking temperature of the system, i.e. close to the ferromagnetic / superparamagnetic transition . The second approach consists in cooling the system down to the ferromagnetic state starting from a paramagnetic regime obtained above the Curie temperature ($T_{C}$) .

Following the second procedure, we show that thermally active artificial kagome arrays of nanomagnets can be locally brought into the magnetic charge crystal phase . Furthermore, we show that, within these magnetic charge crystallites, pairwise spin correlators are very similar to those expected for the so-called spin ice 2 phase, predicted by Monte Carlo simulations at low temperatures and characterized by a magnetic charge crystal embedded in a partially-ordered kagome spin lattice with preferential spin-loop configurations. In other words, besides the local observation of alternating $\pm$`<!-- -->`{=html}1 magnetic charges, the spin configurations are also consistent with those specific to the exotic spin ice 2 phase. However, by computing the pairwise spin and charge correlators on the array-scale, we find that the overall magnetic configurations are clearly out-of-equilibrium. Using a kinetic algorithm that models how our artificial arrays of nanomagnets behave when crossing the Curie temperature from the paramagnetic state, we manage to reproduce very well our experimental data. We thus interpret our experimental findings in terms of kinetic pathways that the system follows upon cooling and we further show how the kinetic process drives the arrays, locally, into low-energy degenerate manifolds that are hardly accessible using a field protocol.

Kagome arrays of 342 nanomagnets were fabricated by e-beam lithography and ion beam etching. Typical dimensions are 500$\times$`<!-- -->`{=html}100$\times$`<!-- -->`{=html}10 nm$^3$ and the nanomagnets are connected at the vertices (see Fig. 1). The constituent material is a ferrimagnetic CoGd alloy which has a Curie temperature adjustable over a wide temperature range by tuning the alloy composition. For the chosen composition (Co$_{0.7}$Gd$_{0.3}$), $T_{C}$ is close to 475 K. Therefore, the networks were annealed at about 500 K, before being cooled down below $T_{C}$ in the absence of the applied magnetic field. Since the remagnetization of the nanomagnets when crossing the Curie point is orders of magnitude faster than the cooling rate (ns and sec time-scales, respectively), the process can be considered as quasistatic. At room temperature, the nanomagnets are uniformly magnetized, with their magnetization aligned with the long axis of the magnetic elements, and can therefore be considered as Ising pseudo-spins. The resulting magnetic configurations were imaged using X-ray Magnetic Circular Dichroism PhotoEmission Electron Microscope (XMCD-PEEM) at the Nanospectroscopy beamline of the Elettra synchrotron radiation facility . A typical XMCD-PEEM image of a kagome array is shown in Fig. 1a. The magnetic configurations imaged after cooling the system through the Curie temperature show that the arrays are efficiently demagnetized and that the ice rule is strictly obeyed: among more than 3800 observed vertices (corresponding to 18 different arrays), no 3-in or 3-out configuration is observed. The magnetic configurations thus fall in the pseudo-spin ice manifold and all vertices are characterized by a $\pm$`<!-- -->`{=html}1 magnetic charge $Q$ .

From these images, the pairwise spin and charge correlators can be determined . First, we have measured the nearest neighbor charge correlator $\langle Q_{i}Q_{i+1} \rangle$ for the 18 arrays we studied. The values range between $-0.30$ and $-0.10$, with an average value of $-0.22$. Although far from the $-1$ value expected for the magnetic charge crystal phase, our measurements are consistent with Ref. and indicate that the magnetic charge has crystallized. In fact, we see the formation of charge domains, with a perfect alternation of positive and negative charges on adjacent vertices (see blue/red circles in Fig. 1b). In several cases, these charge domains can extend over a significant fraction of the array and can include more than 1/3 of the total number of vertices. To visualize the distribution of domain sizes, charge domains are colored in white and green in Fig. 1b.

To gain further insight into the underlying physics, we have also considered the pairwise spin correlators, which we computed by averaging within the crystallites of magnetic charges (locally) and over the entire array (globally). As an emergent charge order can be found in both the ground state configuration and the spin ice 2 phase, the local averages of the spin correlations can help discriminate between the two regimes. Both cases are reported in Fig. 2 up to the seventh neighbor (green and black diamonds, respectively) and each correlation value corresponds to an average performed over the 18 arrays studied. For comparison, the values for the ground state manifolds predicted by the nearest neighbor spin ice (SRSI) model and the dipolar spin ice (DSI) model are also reported (red and blue spheres, respectively) . This experimental data shows several remarkable features.

First, the measured pairwise spin correlators can have large values, even for higher order neighbors (see Fig. 2). By no means can these values be obtained with only a nearest neighbor spin-ice model . Therefore, long-range dipolar interactions cannot be neglected when modeling such arrays.

Second, the spin correlations we measure differ significantly from those of the true ground state of the DSI model. To further determine how far we are from the ground state manifold, we compared our experimental values for the spin correlators with those given by Monte Carlo simulations by employing a correlation scattering analysis . We find that the magnetic configurations we imaged within the charge domains correspond to spin configurations that are very close to the pseudo-ice manifold of the thermodynamic spin-ice 2 phase. In fact, many of our experimental correlations fall within the standard deviations of their corresponding Monte Carlo averages obtained for a temperature of $T/J_{\alpha\beta}$ $\simeq$ 0.034 (see green diamonds and purple spheres in Fig. 2).

Third, if we average the spin correlators computed globally (see black diamonds in Fig. 2), the resulting values significantly differ from those calculated locally, within the charge domains (green diamonds in Fig. 2). In terms of effective temperature, the best fit obtained with our correlation scattering analysis gives a $T/J_{\alpha\beta}$ value of about $0.05$ (note that the $C_{\alpha\nu}$ and $C_{\alpha\delta}$ spin correlators are clearly out of the standard deviations). This $T/J_{\alpha\beta}$ value corresponds to a charge correlator of about $-0.54$ (crossing between the charge correlation plot and the vertical line in the inset of Fig. 3a). This contrasts with the $-0.22$ value that we find experimentally when averaging at a global scale. We then have to solve an apparent contradiction: while pairwise spin correlations measured on the array-scale are those corresponding to configurations approaching the spin ice 2 phase ($T/J_{\alpha\beta} = 0.05$), the charge correlator severely deviates from the at-equilibrium Monte Carlo average found for this temperature (vertical line in the inset of Fig. 3a). In fact, the experimental $\langle Q_{i}Q_{i+1} \rangle$ value alone indicates that the system has barely reached the spin ice 1 manifold (first crossing between the charge correlation plot and the horizontal line in the inset of Fig 3.a), although our systems contain large magnetic charge crystallites. As we shall see further on, these differences are signatures of the kinetics of the thermalization process which leads to an out-of-equilibrium state.

To model the experimental procedure, we make the hypothesis that upon cooling below $T_{C}$, each nanomagnet is subject to thermal fluctuations, but once magnetized, magnetization reversal is no longer possible due to the relatively high energy-barriers at stake. Except for the first nanomagnet, the magnetization of a given nanomagnet is biased by the stray field of its environment . Therefore, we model the full sample magnetization by the following steps : 1 - choose randomly a nanomagnet and set the direction of its magnetization, 2 - calculate the resulting stray field over the entire array and pick the not-yet-magnetized nanomagnet that perceives the highest effective field, 3 - set the direction along which this nanomagnet will orient itself, according to a Boltzmann-like probability, 4 - go back to step 2 and keep repeating steps 2 to 4 until the full sample is magnetized . Once the array is fully magnetized, we calculate the resulting spin-spin and charge-charge correlators on the network scale and repeat these measurements for temperatures ($T/J_{\alpha\beta}$) ranging from $10^{2}$ to $10^{-3}$. The corresponding values of the spin correlators (defined in Fig. 3b) are reported in Fig. 3c with the same color code used for the Monte Carlo simulations.

Although there are some differences between the correlators deduced from the two numerical approaches, they still show striking similarities. This is mostly due to the fact that the spin interactions in both cases are described by the dipolar spin ice Hamiltonian, hence the energy landscape is the same . However, the Monte Carlo approach explores the different configurations at-equilibrium and in an ergodic manner, whereas the kinetic algorithm is a rather one-shot approach, sequentially magnetizing each spin according to a Boltzmann probability in its attempt to minimize the free energy. Therefore, although some correlations can differ in both magnitude and sign for certain temperature windows (see the behavior of the $C_{\alpha\delta}$ and $C_{\alpha\nu}$ correlators for $T/J_{\alpha\beta}$ ranging from $1$ to $0.1$ - Fig. 3), the overall matching is good, especially at low-temperature where our experimental correlations fall.

By performing a correlation-scattering analysis like in the Monte Carlo case, we can place our experimental values for both the pairwise spin and charge correlators on the temperature-dependent variations predicted by the kinetic model. Interestingly, they all agree upon the same effective temperature, $T/J_{\alpha\beta}$=0.56, and a vast majority of them fit within the standard deviations associated to this temperature (see Fig. 3c and inset). The kinetic algorithm thus describes well all our experimental findings and solves the apparent contradiction mentioned above. We emphasize that, if the dipolar interactions are suppressed in the kinetic algorithm, leaving only nearest neighbor interactions at play, the model fails to reproduce our experimental results.

We have finally compared the magnetic configurations we imaged with the ones obtained by applying ac demagnetization protocols, similar to those used in other works . We find out that, for the same arrays, the effective temperature deduced from the analysis of the pairwise spin correlators is about one order of magnitude lower if the thermal kinetic approach is used. This feature raises further questions on the underlying mechanisms that make one procedure more efficient than the other in driving the system to a low-temperature state and clearly deserves more in-depth analysis. However, we have noticed that the kinetic procedure allows local magnetic configurations that are difficult to obtain through the use of a field protocol. Figure 4 shows snapshots of a kinetic simulation performed at low temperatures and at different steps of the magnetization process. The numbers labeling the nanomagnets indicate in which order they remagnetize. In this temperature regime, the magnetization process favors the formation of full hexagons with flux closure magnetic configurations. These local spin arrangements are specific to low-energy manifolds obtained with at-equilibrium Monte Carlo simulations, in which loop-flip algorithms are used to overcome the critical slowing down behavior encountered by single spin flip protocols . Interestingly, when working at the superparamagnetic limit, single spin flips are able to drive building blocks of such artificial systems into their corresponding ground state configurations, but they quickly become inefficient as the system size increases . However, the kinetic process at work in this study spontaneously favors closed-loop configurations. We thus relate the efficiency of our thermalization procedure to the kinetic pathways that the system follows when crossing the Curie temperature of the constituent material.

In conclusion, we have studied the behavior of thermally active artificial kagome spin ice systems by cooling the sample from a high-temperature paramagnetic state down below the Curie point of the constituent material. The resulting magnetic configurations present large magnetic charge crystallites and a detailed correlation analysis shows that our arrays are brought, locally, into the spin ice 2 phase, characterized by the emergence of a magnetic charge order within a disordered spin network. The kinetic processes that take place during the cooling procedure appear as an efficient mean to drive the system into low-energy manifolds where exotic physics emerges, opening new avenues to investigate unconventional magnetism in artificial spin systems.

This work was partially supported by the Region Lorraine and the Agence Nationale de la Recherche through the project ANR‐12-BS04-009 'Frustrated'. I.A Chioar acknowledges financial support from the Laboratoire d'Excellence LANEF Grenoble.


abstract: This paper proposes a potential approach to deploy the Coordinative Real-time Sub-Transmission Volt-Var Control Tool (CReST-VCT) into the energy management system (EMS) and distribution system management system (DMS) under the form of a system integration framework. The framework provides an architectural basis to serve as a guideline for implementing CReST-VCT and other advanced grid support tools to co-optimize the use of DERs and assets in both transmission and distribution domains. As CReST-VCT operation heavily depends on communication infrastructure and data flow, potential communication protocols for different physical domains of a real application is also included. Performance and security issues are discussed to highlight considerations that should be evaluated when applying the CReST-VCT to a specific deployment. Finally, the paper presents an example showing how CReST-VCT and other advanced grid support tools can be integrated in an open-source standardized-based platform that supports distribution utilities.
author:
bibliography: QuanNguyen_References.bib
title: EMS and DMS Integration of the Coordinative Real-time Sub-Transmission Volt-Var Control Tool under High Photovoltaic Penetration

# Introduction

generation has been steadily increasing and reshaping the landscape of power generation. The generation capacity of utility-scale and small-scale photovoltaic (PV) in 2018 was 30 and 20 GW, respectively, which doubles the numbers in 2008 . Therefore, industry, government organizations, and researchers have been developing methods and tools that can increase system operational efficiency and reliability: substation automation, automated fault location identification, volt/var control, etc. Recently, Eaton developed the CYME Volt/var Optimization analysis module that helps find the optimal distribution network operation . In , ETAP introduced a volt/var optimization tool that optimally manages system-wide voltage levels and reactive power flow to achieve efficient distribution grid operation. ABB, Siemens, Schweitzer Engineering Laboratory, and Dominion Energy also developed volt/var optimization solutions to increase distribution system operational efficiency .

As described above, volt-var optimization has been widely implemented in distribution systems. Increasing DER penetration in distribution systems make them more and more dynamic, with different volt-var profiles from traditional centralized supply models. With sufficient distributed energy resources (DER) penetration, DERs can have both positive and negative effects on the upstream transmission system. By controlling the power injected into or absorbed by the distribution systems, transmission system operators in the energy management system (EMS) control center can take advantage of the DERs and also mitigate their negative effects on planning and operation of the transmission system. To achieve an optimal coordinative transmission and distribution (T&D) operation, T&D operators need to establish bidirectional communication and efficiently handle the computational effort required by the operation framework on each side. As the grid assets in transmission and distribution network are operated by different entities within a utility or across separate owners, these capabilities have yet to be implemented. Such a very limited information exchange and coordination exist between current EMS and distribution management system (DMS) control centers impedes the possibility of T&D co-optimization.

Research focusing on closing the gap between T&D has recently been emerging. In , the power flow model and solution of a combined T&D system are discussed. The work in discusses the benefits of improving long-term transmission voltage stability by incorporating a distribution model. In , a coordinated T&D co-optimization framework is developed. In , a coordinated restoration of a T&D system using a decentralized scheme is proposed. As an extension of the existing works on co-optimization for integrated T&D systems, the Coordinated Real-time Sub-Transmission Volt-var Control Tool (CReST-VCT) is proposed to coordinate and optimize the use of var resources in both T&D sides . While these works mainly focuses on the T&D modeling and numerical analysis, aspects such as the communication and practical integration in EMS and DMS are not yet addressed.

In this paper, a high-level integration framework, practical challenges, and deployment approaches for a successful and reliable implementation of CReST-VCT in EMS and DMS are proposed and discussed. This work also serves as a reference for implementing an T&D analytical tool in real system control centers. The general model provided by the proposed framework can be tailored to a specific environment based on its specific constraints and requirements.

The contribution of this paper includes:

-   A generalized design framework for CReST-VCT and other T&D tools in EMS and DMS for enabling high levels of DER penetration.

-   An analysis of potential performance and security issues that should be considered during a practical implementation.

-   An example of implementing the distribution part of CReST-VCT in a open-source platform used for designing, studying, and testing advanced operational grid functionalities.

# Overview of CReST-VCT

## Objective

CReST-VCT is an interactive tool that coordinates and optimizes the action of static and dynamic reactive power components, such as capacitor banks and PV inverters, respectively, from a sub-transmission system down to its distribution feeders . CReST-VCT aims to enable high PV penetration at the distribution feeders while alleviating or eliminating voltage violations in the sub-transmission system.

## Control Architecture

The control architecture of CReST-VCT is shown in Fig. . The CReST-VCT algorithm performs coordinated optimizations between the EMS and DMS at regular intervals. The frequency of operation for the algorithm is flexible, but it is optimal if the chosen interval is consistent with the highest data resolution of load and DER. In , a 5-minute interval is chosen in this work due to the available 5-min solar data.

At the beginning of each control interval, operation limits of the distribution systems are reported to the EMS residing at the sub-transmission control center. CReST-VCT, which is integrated in the EMS, runs an optimal volt-var-control (VVC) algorithm that considers all sub-transmission operational constraints to determine the operating points for var resources in the sub-transmission system as well as send the operating points for the associated distribution systems back to their DMS at the distribution control centers.

At the DMS, a voltage-load sensitivity matrix (VLSM)-based VVC algorithm was developed to achieve the optimal dispatch of DERs at the distribution level . As shown in Figure 6, after receiving the operation points from the energy management system (EMS), the VLSM-based VVC will dispatch DER resources, e.g. smart thermostats, energy storage devices and PV inverter units, to follow the operation point sent to it while satisfying distribution system operational requirements and constraints. The VLSM-based VCC algorithm also will calculate the upper and lower limits of the real and reactive powers of the distribution feeders and forecast the next real power operating point for the next time step. That information is sent back to EMS to be used in the next control interval.

## Required communication and time flow

Preliminary tests of CReST-VCT indicate the following time flow for one operating interval:

-   The CReST-VCT algorithm at a EMS takes about 45 seconds when applying for a large-scale sub-transmission system with 3246 buses .

-   About 15 seconds are expected for sending the request on substation voltage, PV curtailment, and demand response from EMS to DMS. Because an EMS might communicate with multiple DMSs, the communication may be staggered within those 15 seconds.

-   At a DMS, solving the VLSM-based volt-var control algorithm takes about 30 seconds.

-   The bidirectional messages about real-time dispatch and supporting reactive power capacity between DMS and downstream DERs may take up to 120 seconds over various elements of communication infrastructure.

In addition, the information about the limits of demand response, PV curtailment, and reactive power support from DMSs to EMS, which is used for the next interval, may take about 60 seconds.

# Proposed CReST-VCT Design Framework at EMS and DMS

In this section, a framework for an implementation of CReST-VCT in an operational system is presented. The framework outlines one possible system design approach that defines the architecture structure, functional allocation to components, and resulting interfaces. While only one approach is described, key considerations that may drive design choices are presented throughout to help tailor it to a particular operational environment.

## CReST-VCT system structure

### EMS Subsystem

CReST-VCT would reside as one of the advanced applications in an EMS. It would have access to the data monitored by supervisory control and data acquisition (SCADA) and be able to store data in a historian or database. CReST-VCT would also require the state of the system from the state estimator application.

The current implementation of CReST-VCT acts as a solver for optimal power flow (OPF) using the General Algebraic Modeling System (GAMS) modeling language and commercial solver KNITRO. Many EMS systems are also equipped with OPF tools. CReST-VCT can either leverage the built-in tools or maintain stand-alone functionalities, so as to easily port between various vendors.

It is not feasible to specify the application program interfaces (APIs) for the CReST-VCT advanced application to be included within the EMS system, because this is usually proprietary information of the EMS system providers. An option to minimize API issues is to containerize the CReST-VCT advanced application using a platform like Docker. This would provide a common implementation that could be ported to a variety of different EMS systems. However, the EMS integration would still require translation from the EMS's unique APIs to those established by the container. Though CReST-VCT is not available as a container at this writing, that can be implemented if desired.

### DMS Subsystem

Just as in an EMS, the VLSM-based VVC algorithm at the distribution side of CReST-VCT resides in the DMS. As shown in Fig. , this application takes into account requirements from CReST-VCT and allocates these requirements to its distributed resources so that the demand-response cost and the overall voltage deviation are minimized.

The CReST-VCT component will require telemetry data for DER and distribution assets through the SCADA system. It is assumed the telemetry data is available within the DMS and the CReST-VCT component will interface using an appropriate data access API. To issue new operating points to the distribution assets in accordance with the CReST-VCT output, the DMS may communicate directly to DERs or may interface through a distributed energy resource management system (DERMS). The protocols used for these interfaces to these devices are discussed in Section IV.

Similar to the EMS integration, the details of the internal interfaces between the CReST-VCT and the DMS system are expected to be dependent on the vendor-specific DMS software and are not specified here. It is assumed that the CReST-VCT software will need to conform to the appropriate vendor-specific API. Again, this does not preclude the existing CReST-VCT software implementation from being implemented in a container as a stand-alone element or the DMS from implementing appropriate interfaces to the container.

## CReST-VCT system interfaces

In this section, the communication interfaces between the components used by the CReST-VCT solution are identified. These communication flows cross several different communication network domains in the utility operational environment and use several different protocols. The diverse communication network landscape was discussed in . For reference, the high-level network diagram is shown in Fig. .

### EMS-DMS interface

Connection from EMS to DMS is often point-to-point, dedicated connection. The messages that would be passed through this channel include the reactive power required to maintain voltage stability in the sub-transmission network, substation voltage, aggregated PV curtailment, and demand-response requirement. No generally accepted standard protocol was identified for this interface in our research. This may be a sign of the relatively new nature of coordination between EMS and DMS. Because both EMS and DMS systems use Distributed Network Protocol 3 (DNP3) to communicate with power system components, our assumption for this paper is that the EMS-to-DMS interface to support CReST-VCT would be implemented in DNP3.

### DMS to inverters

Connection from DMS to DERs is an area of rapid evolution. Standards are being updated or developed to take into consideration the variability in the many DERs systems from different vendors. California Rule 21 and Hawaii Rule 14 have mandated communication capabilities and availability of advanced functions in solar inverters with remote controllability . These developments are timely, as the penetration of PVs has been increasing, especially in the states of California and Hawaii. The IEEE Standard 2030.5 defines open standards for communications and data format for DERs . IEEE Standard 1547-2018 specifies standard choices for communication at the DER as well as mandatory advanced capabilities . Such standards support realization of requirements such as those specified by California Rule 21. However, a variety of different protocol options are still in use today.

Utilities might face interoperability challenges for communication networks. Each utility standardizes internally, adopting the collection of standards that best suits their particular environment, and uses it to dictate procurement decisions. Therefore, the protocols used to integrate CReST-VCT into a utility's operations may differ according to the network architecture and technology of that given utility. It is expected that the chosen set of protocols and the communication networks used will need to be able to handle various configurations of DERs.

# Performance Analysis

## Potential communication variation

CReST-VCT operation relies on a bidirectional flow of operational parameters and constraints between the T&D sides, as mentioned in Section II.C. The performance of the system depends greatly on the latency of the communication process. The faster a message flows between the two sides, the better performance system achieves.

However, there is a significant potential variation in the performance of communication networks that can be adopted for CReST-VCT applications. First, the communication latency and bandwidth vary with the chosen protocol and the layer of the network where the control resources of CReST-VCT, such as an inverter-based PV, are connected. For example, a DER connected at the neighborhood area network layer will have significantly more latency than a device connected at the substation local area network. Feeder devices such as regulators, capacitor banks, or DERs connected at the field area network layer or the neighborhood area network layer are subject to the widest range of possible communication network performance. Advanced metering infrastructure networks also contribute to the variation because their data rates can range from 10 kbps to 1.2 Mbps. Second, the specific geographic location with reference to network infrastructure, the surrounding topology, and network utilization also introduce variations in performance in a communication network such as a wireless network .

## Effects of communication delays

As described above, CReST-VCT operation spans multiple tiers of a typical utility's communication network. Different communication system technologies may be used in each, with its own performance and reliability characteristics.

Based on control architecture in Section II.B and , if solving the optimization problems at both EMS and DMS and sending data through the communication network takes $x$ minutes, the system is only expected to achieve the minimum losses while maintaining all operational constraints in the remaining (5 - $x$) minutes. As a result, the CreST-VCT solution will be more meaningful in minimizing losses and satisfying all the operational constraints over a longer period if $x$ is significantly small compared to 5 minutes.

## Scalability

CReST-VCT assumes that the DMS is fully capable of establishing two-way communication with all its control resources. Therefore, centralized control in a distribution system will reach a scalability limit, because the computational effort and communicational message flow are constrained by the 5-minute interval in CReST-VCT input to guarantee a desired level of system performance. To make up for this insufficient communication, the distribution side might devise decentralized and distributed approaches to control DERs locally with significantly less stringent or no communication requirements, respectively. However, these approaches only consider local state measurements rather than comprehensive system network constraints. Therefore, it cannot match the CReST-VCT solution in terms of minimizing losses and guaranteeing all operational constraints, such as voltage stabilization. Another possible solution that is worth to considering is to make a compromise between the nominal 5 minutes specified in CReST-VCT design and the scale of hosted inverter-based PVs. With an extended interval such as 10 minutes or 15 minutes instead of 5 minutes, fewer limitations are put on communication and the number of DERs.

The CReST-VCT tool also provides a mechanism to parallelize the algorithm computation by dividing the system into substation groups. However, each of these groups is expected to be synchronized to the same five-minute operational interval to make sure the DER dispatch and supporting field measurements are coordinated across the grid. This means that the DER dispatch and field measurements will be executed in parallel. This could pose scalability issues for the communication network infrastructure.

SCADA systems at the substation level operate at the typical poll cycle of one to several seconds for all points. Because these points are maintained continuously, it is assumed that communication networks will be scaled appropriately for normal operations, independently of CReST-VCT. Connection to devices outside the substation to take advantage of CReST-VCT optimization capabilities may drive new connections to the SCADA system. Further, CReST-VCT increases the potential hosting capacity of DERs that might otherwise be curtailed, increasing the possible number of DERs a given feeder can support. The potential optimizations with a five-minute interval could again provide incentive to connect more devices and interact with those devices more often.

## Security considerations

CReST-VCT has features such as remote control and remote access to DERs, and those DERs are equipped with digital communication and control interfaces that pose as cyberattack surfaces. Also, larger numbers of third-party devices associated with DER deployment can also increase cybersecurity risk if the devices are not sufficiently secure. Communication protocols used for CReST-VCT are vulnerable to attacks by attempts to intercept, modify, and/or corrupt the control signal packets. The U.S. National Electric Sector Cybersecurity Organization Resource (NESCOR) Technical Working Group 1 describes some realistic cyber threats that concerning the DER domain . Each of the Open Systems Interconnection (OSI) communication layers in the DER domain is vulnerable to potential cyberattacks. Cybersecurity should be considered with any type of EMS/DMS applications that are related to DER devices and their respective communication channels. The current cybersecurity postures of IEEE 2030.5, IEC 61850, and IEEE 1815 are discussed below.

# Implementation of the Distribution Side of CReST-VCT in GridAPPS-D

This section describes at a high level an example of how to integrate the distribution side of CReST-VCT in the open-source, data integration, and application development platform GridAPPS-D.

## Overview about GridAPPS-D

As mentioned in the Introduction about the structure and characteristics of future power systems, it is important to enable a mature, efficient, and validated procedure for the development of future grid support functions for utility operators. Therefore, GridAPPS-D was introduced in , and it enables distribution system application developers including commercial vendors, utilities, and research institutions to develop innovative applications related to planning and operation:

-   Researchers can use GridAPPS-D as a standardized platform to develop and test new advanced functionalities.

-   Application developers can use GridAPPS-D to streamline a development of an application for a full range of distribution systems.

-   System vendors can adopt the architecture and interface definitions of GridAPPS-D to support portable applications.

-   As GridAPPS-D-compliant products enter the marketplace, utilities will be able to integrate more advanced tools that meet their specific requirements with significantly reduced cost and time into their operations.

In general, GridAPPS-D™ platform aims to shift utilities away from a closed, vertically integrated architecture to collaborating in an open, layered, and standards-based architecture.

## Integrated design view for CReST-VCT in GridAPPS-D

As mentioned in Section 4.1, the key components of the CReST-VCT is a VLSM-based VVC algorithm installed at the DMS and a CReST-VCT algorithm installed at EMS.

For the integration of the co-simulation CReST-VCT in GridAPPS-D, the VLSM-based algorithm will be integrated as a new distribution system application, among other existing applications. It is a high-level application that iteratively receives the operating point determined from the EMS through the Data Ingest Interface, Standard-based Data Bus, and I/O Interface components inside GridAPPS-D. The VLSM-based algorithm then determines the optimal dispatch for the DER resources, e.g. smart thermostats, energy storage devices, and PV inverter units in the distribution model for each 5-minute period. This dispatch information is delivered to the distribution simulator, either in GridLAB-D or OpenDSS, via the I/O interface and Standard-based Data Bus components.

On the other hand, the CReST-VCT algorithm at the EMS is considered as an external tool that constantly coordinates and communicates with the VLSM-based VVC application by exchanging bidirectional data through the Data Ingest component of the GridAPPS-D. In addition, the transmission simulators such as PSSE, PSLF, and PowerWorld can also link and exchange information to the distribution system simulators such as GridLAB-D and OpenDSS to replicate the physical connection in practice.

# Conclusion

In this paper, a system integration framework for integrating the CReST-VCT and other advanced grid support function tools into the real-time operational environment of EMS and DMS. The framework serves as a guideline on designing system structure and interfaces that is consistent with the active and data driven characteristics of future transmission and distribution system under high penetration of DERs. Detailed analysis about communication variation, negative effects of communication delay, scalability, and security that utilities should consider for a practical deployment. To further clarify the proposed framework and analysis, the paper presents a high-level description about how the CReST-VCT can be integrated in GridAPPS-D, which is an open-source and standardized platform supporting system application developers, EMS, and DMS in testing and implementing new innovative supporting tools.

# <span style="background-color: yellow">Acknowledgment</span>


author: G. Giovannini , L. Feretti , M. Girardi , F. Govoni , M. Murgia , V. Vacca, J. Bagchi
date:
title: A giant radio halo in the low luminosity X-ray cluster Abell 523

We observe and analyze the diffuse radio emission present in the complex merging structure Abell 523, classified as a low luminosity X-ray cluster, to discuss its properties in the context of the halo total radio power versus X-ray luminosity correlation. We reduced VLA archive observations at 1.4 GHz to derive a deep radio image of the diffuse emission, and compared radio, optical, and X-ray data. Low-resolution VLA images detect a giant radio halo associated with a complex merging region. The properties of this new halo agree with those of radio halos in general discussed in the literature, but its radio power is about a factor of ten higher than expected on the basis of the cluster X-ray luminosity. Our study of this giant radio source demonstrates that radio halos can also be present in clusters with a low X-ray luminosity. Only a few similar cases have so far been found . This result suggests that this source represent a new class of objects, that cannot be explained by classical radio halo models. We suggest that the particle reacceleration related to merging processes is very efficient and/or the X-ray luminosity is not a good indicator of the past merging activity of a cluster.

# Introduction

Clusters of galaxies are characterized by X-ray emission from a hot intra-cluster medium (ICM, $T \sim 2{-}10$ keV). Thermal emission is a common property of all clusters of galaxies and has been detected even in poor groups as well as in optical filaments connecting rich clusters. The ICM is believed to be shock heated by merging during the structure formation.

Diffuse non-thermal radio sources with steep spectra have been detected in the central regions of a large number of clusters of galaxies (see e.g. , (and refs therein) and review papers such as , and ). They are giant radio sources with a spatial extent similar to that of the hot ICM, which are called radio halos. In addition, diffuse radio sources classified as radio relics are detected in cluster peripheral regions (see e.g. and ). These sources are not directly associated with the activity of individual galaxies and are related to physical properties of the whole cluster. There is a substantial evidence that radio halos are found in clusters with significant substructures in X-ray images and complex gas temperature distributions, both signatures of cluster mergers (e.g. , , , ). Classical radio halos show a centrally located and regular radio morphology. However, in several cases the radio halos are irregular and elongated as in A209, A401, and A2034. A few radio halos such as A851 and A2218 are highly asymmetric, being located mostly to one side with respect to the cluster center. interpreted these structures as possibly reflecting the geometry of the merger.

discussed that the radio power of halos at 1.4 GHz increases with the X-ray luminosity of the parent cluster, implying a direct connection between the radio and X-ray plasmas. Since the cluster X-ray luminosity and mass are correlated, it follows that radio halo power correlates with cluster mass (; ). The strong correlation between the radio halo total power at 1.4 GHz and the X-ray cluster luminosity was recently discussed using higher quality statistics by and , who show that the correlation is also present for low power radio halos (P$_{1.4}$ $<$ 10$^{24}$ W/Hz and L$_x$ $\sim$ 10$^{44}$ erg/s).

were the first to discuss whether a radio halo was present in A1213, an under-luminous X-ray cluster i.e. well outside the above correlation. A similar case was found in the cluster 0217+70 by . We do not consider radio relics because present models connect radio halos to turbulence in the ICM, while relic sources are correlated to the presence of shock waves. discussed whether the low X-ray luminosity in 0217+70 could be due to an underestimate of the X-ray luminosity because of the hydrogen column density, which might be patchy at these low Galactic latitudes. Therefore, further X-ray observations are required to measure the absorption directly. We note that discussed the peculiarity of the diffuse emission in A1213.

We present new VLA images from archive data of diffuse radio emission, found by chance by inspecting NVSS images, in A523, one more radio halo in a low X-ray luminosity cluster. We present in Sect. 2 the radio data and the X-ray and optical information. Discussion and conclusions are given in Sect. 3.

The intrinsic parameters quoted in this paper are computed for a $\Lambda$CDM cosmology with $H_0$ = 71 km s$^{-1}$Mpc$^{-1}$, $\Omega_m$ = 0.27, and $\Omega_{\Lambda}$ = 0.73. At z =0.10, the luminosity distance is 455 Mpc, and the angular conversion factor is 1.82 kpc/arcsec.

# Radio images

A523 was observed at 1.4 GHz in line mode in 2005, with the VLA in the C (5 hrs) and D (2 hrs) configuration (project AB1180). Calibration and imaging were performed with the NRAO Astronomical Image Processing System (AIPS). After the editing of bad points, self-calibration was applied to both datasets, to remove residual phase and gain variations.

In Fig. , we show the final image obtained with C configuration data and natural weights (contours), superimposed on the X-ray image (colors) obtained with the Rosat satellite. The diffuse emission is also visible at this relatively high angular resolution and discrete sources are well separated. To obtain a deeper image of the diffuse source, we performed a full resolution (using uniform weighting) image with C array data, cutting out in addition the shortest baselines where some flux density from the extended halo was present. In this image, only discrete sources are present, and we used it to select the clean components of sources in the halo region with the CCEDT task. These components were subtracted from the D array uv-data, and we used resulting data to produce a low resolution image of the diffuse component (see Fig. ). We used this image to measure the halo flux density and size. The main problems in this procedure are the following: 1) The presence of extended emission from a head-tail source located near the center of the radio halo. We compared the flux density of this source with the clean components used to subtract it from D-array UV data. We are confident to have subtracted more than 95$\%$ of the flux density of this extended tail source. 2) The elongated halo extension in the extreme east region. In this position, two discrete sources are present, one coincident with a X-ray source. However, the residual flux density in the low resolution image is too high to be interpreted as the residual of the subtracted sources. In any case, including this region or not will change the halo flux density by only $\sim$ 2.7 mJy.

The estimated halo flux density is 59 $\pm$ 5 mJy, where the error includes the uncertainties in the source subtraction), corresponding to a total radio power = 1.47 $\times$ 10$^{24}$ W/Hz (Log P = 24.17). The morphology is irregular with a largest angular size of $\sim$ 12 arcmin in the EW direction, corresponding to $\sim$ 1.3 Mpc.

The radio emission detected at high resolution and shown in Fig. 1 permeates both merging clumps, as found e.g. in the Bullet cluster (). This is also similar to the case e.g. of A2255 () or A665 () where the radio halos are elongated in the direction of the cluster merger. The radio structure in A523 is therefore typical of radio halos. When degraded to lower resolution, owing to the improved signal-to-noise ratio, the low brightness radio emission is more extended in particular in the E-W direction. The elongated structure does not show any morphological feature typical of radio relics such as high brightness filamentary structures or a transverse flux asymmetry (see e.g.). We conclude that its classification as a radio relic is very unlikely. We also note that a fainter X-ray emission, not detected by ROSAT, could be present. We therefore classify the diffuse radio emission in A523 as a giant radio halo. The size and radio power of the halo are in good agreement with the correlation presented for radio halos by .

## Optical and X-ray data

Optically, A523 is a very rich cluster of Abell richness class $=2$ (Abell et al. ). It was classified to be of the Rood-Sastry morphological type "L", i.e. to have a linear configuration of galaxies (Struble & Rood ). To date, only the redshift for the brightest cluster galaxy \[BCG at R.A.=$04^{\mathrm{h}}59^{\mathrm{m}}12\,\rlap{\hbox{$\mathrm{^s}$}}{\hbox{$.$}}\,97$, Dec.=$+08\ensuremath{\mathrm{^\circ}}
49\ensuremath{\mathrm{^\prime}}41.3\ensuremath{\ensuremath{\mathrm{^\prime}}\hskip -0.1em\ensuremath{\mathrm{^\prime}}\;}$ (J2000.0)\] has been acquired ($z=0.1036$; Crawford et al. ). No radio emission is clearly associated with this galaxy. The visual inspection of DSS2-red image shows the presence of two clumps. To quantify this impression, we compiled a photometric catalog in a large region around the BCG ($\sim 30\ensuremath{\mathrm{^\prime}}$) from the SuperCOSMOS Sky Surveys (SSS[^1]) extracting the objects classified as "galaxies" that have both $B_j$ and $R$ magnitudes available. In this photometric catalog, we selected likely cluster members on the basis of the color–magnitude relation (hereafter CMR), which indicates the early-type galaxy locus. To determine the CMR, we considered the 75 galaxies within a radius of 4$\mathrm{^\prime}\;$from the BCG and applied a two-sigma-clipping fitting procedure, obtaining $B_j$–$R$= 2.044-0.024$\times R$ for a subsample of 54 galaxies. The color-magnitude diagram is shown in Fig. .

From the photometric catalog, we considered as "likely cluster members" those galaxies lying within 0.25 mag of the CMR. Fig.  shows a zoomed region of the contour map for the likely members (734 galaxies within the whole $\sim$ 30$\mathrm{^\prime}\;$ region). The galaxy distribution reveals two significant peaks $\sim
5\ensuremath{\mathrm{^\prime}\;}$ ($\sim 0.5$ Mpc) far: the NNE peak, coincident with the BCG, and the SSW peak, with no dominant galaxy, at R.A.=$04^{\mathrm{h}}59^{\mathrm{m}}05\,\rlap{\hbox{$\mathrm{^s}$}}{\hbox{$.$}}\,4$, Dec.=$+08\ensuremath{\mathrm{^\circ}}
45\ensuremath{\mathrm{^\prime}}06\ensuremath{\ensuremath{\mathrm{^\prime}}\hskip -0.1em\ensuremath{\mathrm{^\prime}}\;}$ (J2000.0). The SSW peak has the higher density and a larger population – by factors of 1.2 and 2.7 according to the 2D adaptive–kernel method (2D DEDICA, Pisani ).

Figure  also shows the 105 galaxies of the photometric sample within a radius of 4$\mathrm{^\prime}\;$from the SSW peak, whose properties can be accurately described by the above CMR, too, suggesting that both clumps are at a similar redshift.

report a cluster X-ray luminosity in the 0.1-2.4 keV band of 1.07 $\times$ 10$^{44}$ erg/sec (after cosmological corrections). report a slightly lower value of 0.9 $\times$ 10$^{44}$ erg/sec. The hot ICM (see Fig. ) is clearly bimodal in agreement with the optical galaxy distribution. However, a clear shift between the galaxies and gas distribution is present (see Fig. ) and the BCG is not coincident with the X-ray peak. The SSW structure is more extended and shows an irregular shape. The NNE clump is more compact and regular. The bright source, visible in the X-ray map at R.A.=$04^{\mathrm{h}}59^{\mathrm{m}}37\,\rlap{\hbox{$\mathrm{^s}$}}{\hbox{$.$}}\,7$, Dec.=$+08\ensuremath{\mathrm{^\circ}}
45\ensuremath{\mathrm{^\prime}}50\ensuremath{\ensuremath{\mathrm{^\prime}}\hskip -0.1em\ensuremath{\mathrm{^\prime}}\;}$ (J2000.0) is identified with a discrete unrelated radio source.

We conclude that optical and X-ray data that we have presented indicate a merging cluster structure, where the main cluster is identified with the SSW structure. This cluster as shown by its irregular shape and lack of a dominant galaxy, was not in a relaxed stage, and is now strongly interacting (major merging) with a more compact cluster (the NNE structure) dominated by a bright BCG.

# Discussion

We present in Fig. the plot of the total radio power and X-ray luminosity for radio halos as shown by , where we have included 0217+70 and A523. The dots refer to classical powerful radio halos in X-ray luminous clusters, which have been found to show a correlation between the radio power and the X-ray luminosity. The red triangles are outliers, i.e. they refer to the few known radio halos that are overluminous in radio with respect to the empirical radio - X-ray correlation.

We focus here on the objects that are overluminous in radio with respect to the empirical radio - X-ray correlation. We note that there might be some uncertainties in the computation of their parameters. In particular, A1213 displays a peculiar diffuse emission, and its size and radio power is lower than that of classical giant radio halos. Even if showed that small-size radio halos have the same properties as giant radio halos, data of higher statistical quality are necessary to confirm this result. In addition, the cluster 0217+70, as discussed by , is on the galactic plane therefore the X-ray luminosity could be affected by an unusually high absorption. On the other hand, the uncertainties in the determination of the radio power and the X-ray luminosity in the cluster A1351 () are smaller than the deviations from the correlation (however see ).

In the case of A523, we are aware that the redshift is not well known, since only the redshift of the BCG is available. However, to be consistent with the correlation between the radio power and the X-ray luminosity, the clusters should be at a redshift of $\sim$ 0.2. Fig.  shows the predicted values of the $B_j$–$R$=1.4 and 2.0 for the predicted early-type galaxies at $z=0.0$ and $z=0.2$, respectively (see Fig. 1 of Brown et al. and refs. therein) after a 0.2 mag decorrection for the Galactic absorption. Although these predictions cannot be used to provide a redshift estimate, they are useful to show to the reader the difference expected in the color–magnitude relation for a $z\sim0.2$ cluster, thus strongly favouring the location of A523 at much lower distance. Moreover, at z = 0.2, the halo size would be $\sim$ 2.4 Mpc, the largest radio halo known to date (see ).

The radio halo in A523 therefore represents a robust case of a giant radio halo clearly identified with an X-ray underluminous cluster. As discussed in Sect. 2, the X-ray and optical distribution is similar (bimodal) but not in perfect agreement. The thermal gas distribution is shifted with respect to the galaxy distribution; this is more evident for the smaller and more compact NNE clump, where the position of the BCG is not at the center of the gas distribution. A higher quality X-ray image is necessary to investigate whether the peculiarity of A523 could be related to the possibility that the compact NNE group just crossed the SSW cluster as in the case of the Bullet cluster (1E0657-56, ), and this strong and recent interaction affected in a different way the gas and galaxy distribution.

We conclude that there is firm observational evidence that giant radio halos can be associated with low luminosity X-ray clusters. These halos are overluminous in radio by at least an order of magnitude with respect to that expected from the extrapolation of the observed radio power - X-ray luminosity.

The radio - X-ray correlation is consistent with reacceleration models, which have been extensively discussed in the literature (Cassano 2010 and refs therein) and supported by the observations (e.g. ). In the framework of these models, the radio emitting particles in cluster radio halos gain energy from merger-driven turbulence in the ICM. Turbulence development on timescales of the order of 1 Gyr is necessary to reaccelerate electrons to the energies needed to emit the observed synchrotron radiation at GHz frequencies, thus giant radio halos are expected to be strictly connected to massive systems that have undergone strong merging processes. In the framework of turbulent reacceleration models, Cassano 2010 suggests that very steep spectrum radio halos should exist that are detectable only at low frequencies, and be less luminous than predicted by the radio - X-ray correlation. The radio halo in A523, and a few similar objects, are instead more luminous in radio than predicted by the radio - X-ray correlation. These powerful radio halos, associated with clusters of low X-ray luminosity, do not appear to be described well by current models, hence are likely to represent a new class of objects, raising new questions about the origin of radio halos. They could be either young halos or clusters at a special time of the merger event, when particle acceleration processes have a higher efficiency (see ). Another possibility is that the X-ray luminosity might not be in these cases a good indicator of the previous cluster merging activity.

We note that the non-thermal diffuse radio emission in A523 does not follow the optical and the X-ray emission, unlike most other radio halos. The brighter halo region visible at high resolution (Fig. 1) is in good agreement with the X-ray bimodal distribution, but the more extended low brightness halo emission, which is most clearly imaged at low resolution, is elongated in the E-W direction i.e. in the direction perpendicular to the merger. This morphology could be related to peculiar cluster conditions that give rise to these overluminous radio halos and needs to be studied in better detail with more sensitive X-ray and radio data and radio spectral information.

[^1]: http://www-wfau.roe.ac.uk/sss/index.html