abstract: The transition from Asymptotic Giant Branch star to Planetary Nebula is short-lived and mysterious. Though it lasts only a few thousand years, it is thought to be the time when the asymmetries observed in subsequent phases arise. During this epoch, the star is shrouded in thick clouds of dust and molecular gas; infrared observations are needed to reveal these objects at their most pivotal moment. I present preliminary results of a Spitzer study of targets spanning the range from post-AGB stars to Planetary Nebulae with the goal of determining the genesis of asymmetry in these objects.
author: Alexa Hart$^{1,2}$, Joe Hora$^2$, Luciano Cerrigone$^3$, Grazia Umana$^4$, Corrado Trigilio$^4$, Martin Cohen$^5$ and Massimo Marengo$^6$
bibliography: apn5v1.bib
title: Proto-Planetary Nebulae with the Spitzer Space Telescope

# The Evolution of Intermediate-Mass Stars

The life cycle of a star and the manner of its ultimate demise depends on its mass. Intermediate-mass stars gradually fling their outer layers into space, forming planetary nebulae.

At the end of the Giant phase, the envelope of intermediate mass stars (defined as $0.8 \mbox{--}8 M_{\odot}$) expands and the core contracts, resulting in a degenerate C/O core. This is the beginning of the Asymptotic Giant Branch (AGB) phase; the star will start to burn nuclear fuel in shells surrounding the spent core. During this stage, there is a slow, dense and (presumably) spherically symmetric stream of material leaving the star called the AGB wind. This wind is interspersed with shells of enhanced density due to pulsation. These shells have been observed in many Planetary Nebulae such as the Cat's Eye Nebula (Balick et al., 2001).

As the mass loss process continues, the thick cloud of dust and gas it creates begins to obscure the star itself in visible wavelengths (Garcı́a-Lario, 2006). Eventually, the central star will have lost so much mass that the AGB wind will cease; this marks the beginning of the post-AGB phase (Kwok et al. 1993). Though we can conjecture about what happens next, we have precious little observational evidence for it due to the obscuration of the central star. However, post-AGB stars observed in the mid-infrared with IRAS apparently fit into an evolutionary sequence with effective temperatures steadily increasing until the star is hot enough to ionize its envelope (spectral types marching from K, G, F, A to B). Near the earliest onset of ionization, the star becomes a Transition Object, sometimes called a proto-Planetary Nebula. The transition phase lasts mere thousands of years.

The planetary nebula phase commences when the star ionizes its envelope. The ionized material is optically thin, so the central star re-emerges in visible wavelengths once the ionization front has traveled sufficiently far out into the envelope. Morphologies of these young PNe are diverse: bipolar, elliptical, point-symmetric and spherical.

Though few AGB stars have been spatially resolved, they are generally assumed to have spherically symmetric outflows (as long as they are single). Yet when the star emerges just a short time later, it is suddenly the central star of a spectacular, and often highly asymmetric, Planetary Nebulae. Somehow, during this phase of mere thousands of years, the structure of the winds, and perhaps the star itself, change fundamentally. What happens behind that dark cloud of dust and gas is revealed in infrared observations with the Spitzer Space Telescope.

# The Sample

A sample of targets was chosen to span the entire temporal sequence from star to nebula. This sequence was divided into three constituent epochs: post-AGB stars, true Transition Objects, and Planetary Nebulae. Each of these epochs were the subject of an observational program on *Spitzer*, divided as follows: 26 post-AGB stars, 36 Transition Objects, and 18 Planetary Nebulae.

The post-AGB sample was designed to catch stars that are nearing the end of the post-AGB phase, that is, stars that are on the brink of ionizing their envelopes. To find such stars, we started from targets that were characterized as post-AGB candidates based on their IRAS colors (Suarez et al 2006) and then further constrained them to have a central star of spectral type A and a far-IR excess typical of PNe (Parthasarathy 2000). The spectral type condition should limit the sample to objects with very hot central stars, which should be those that are at the very end of the post-AGB phase. Type A stars should not generally be hot enough to ionize surrounding material.

To isolate true Transition Objects, i.e. stars that have just begun to ionize their circumstellar envelopes, a sample of post-AGB candidate (again selected based on IRAS colors and infrared excess) with slightly hotter spectral type B central stars were chosen. These stars should be just hot enough to emit ionizing radiation. As an independent measure of ionization, these objects were observed with the VLA; 16 of the 36 were detected, confirming the presence of ionized shells (Umana et al. 2004; Cerrigone et al. 2008).

The chosen Planetary Nebulae all have H-deficient central stars of late, carbon-rich Wolf-Rayet spectral type (\[WC 8–12\]). These PNe may also have in common a peculiar spectroscopic property: simultaneous evidence for both oxygen-based crystalline silicates and carbon-based PAHs in their dust.

# Dust Chemistry

An important observational aspect of these stars is their dust chemistry, that is, if the dust around them is Oxygen-rich or Carbon-rich. This distinction is a consequence of the high binding energy of CO; since the CO molecule will form first, only the more abundant of C or O will be left to form other dust grains. All AGB stars start out O-rich; stars with roughly $1.5 \mbox{--}4 M_{\odot}$ can become C-rich via the Third Dredge Up.

In Carbon-rich sources, we see Polycyclic Aromatic Hydrocarbon (PAH) features at 3.3, 6.2, 7.7, 8.6 and 11.3 $\mu$m$\,$ (Hrivnak et al. 2007), as well as an SiC feature at 11.3 $\mu$m$\,$ (Speck et al. 2009). Oxygen-rich sources have silicate features at $\sim$`<!-- -->`{=html}10-18 $\mu$m$\,$ and show evidence for water ice. A small minority of stars show evidence of both PAHs and silicates in their dust (called dual-dust objects), a phenomena which has yet to be understood. These dual dust spectra could be indicative of either an O-rich environment preserved in a disk while the rest of the envelope transformed to C-rich or simply the formation of PAHs in a C-rich environment as a result of the photo-dissociation of CO molecules.

Interestingly, Cerrigone et al. 2009 found that approximately $40\%$ of the Transition Object sample were dual-dust objects, a striking contrast to the global occurence of $\sim10\%$ among PNe. The authors noted that this might be due to a selection effect, in particular the choice of objects with a far-IR excess which may be indicative of the presence of a disk. Additionally, $70\%$ of those dual-dust objects were radio-detected in free-free emission, establishing an apparent correlation between ionization and the dual-dust phenomenon. Modeling of the PAH features resulted in the conclusion that the PAH molecules are located in the outflows, far from the central star. For the detailed analysis, see Cerrigone et al. 2009.

The PNe sample was chosen with this peculiarity in mind: Cohen et al. 2002 found that four out of six PNe with central stars of type \[WCL\] in their study were dual-dust objects; the PNe sample was designed to further explore this apparent connection.

# Observations & Results

The sample of 80 targets were first observed with *Spitzer*; follow-up observations of 30 targets were performed using the Magellan I & II Telescopes at Las Campanas Observatory in Chile.

The three programs that were taken with *Spitzer* (Werner et al. 2004) that form the basis of this project were completed in 2009. The sample comprises about 75 targets in total. All have infrared spectra covering 5–40 $\mu$m$\,$ measured on the Infrared Spectrograph (IRS) (Houck et al. 2004) and were also observed with the Infrared Array Camera (IRAC) imager (Fazio et al. 2004), which images each object in four broadband channels centered at 3.6, 4.5, 5.8, and 8.0 $\mu$m$\,$ .

The Magellan Clay telescope with the MMIRS (McLeod et al. 2004) was used to obtain near-infrared (H and K band) spectra of the observable subset of the *Spitzer* sample. In addition, the Magellan Baade telescope was used to image the subset of sources in the *Spitzer* sample that have extended emission using near-Infrared narrowband filters that isolate $H_{2}$ and $Br \gamma$ emission with PANIC (Martini et al. 2004).

Results of the Transition Object program have been published by Cerrigone et al. 2009. The imaging search (with IRAC and PANIC) for faint, extended emission around the post-AGB sample yielded a null result; SED and spectral analysis of the data is underway. A sample IRAC image from the Planetary Nebula program is shown in Figure 1. SED and spectral analysis is underway.


abstract: Network testing plays an important role in the iterative process of developing new communication protocols and algorithms.
  .
  However, test environments have to keep up with the evolution of technology and require continuous update and redesign. In this paper, we propose COINS, a framework that can be used by wireless technology developers to enable continuous integration (CI) practices in their testbed infrastructure. As a proof-of-concept, we provide a reference architecture and implementation of COINS for controlled testing of multi-technology 5G Machine Type Communication (MTC) networks. The implementation upgrades an existing wireless experimentation testbed with new software and hardware functionalities. It blends web service technology and operating system virtualization technologies with emerging Internet of Things technologies enabling CI for wireless networks. Moreover, we also extend an existing qualitative methodology for comparing similar frameworks and identify and discuss open challenges for wider use of CI practices in wireless technology development.
author:
bibliography: IEEEabrv.bib; testlab-testbed-infrastructure.bib
title: COINS: ContinuOus IntegratioN in wirelesS technology development

# Introduction

The next generation of mobile networks will have to support industrial and traffic safety applications, enhanced multimedia applications on smartphones and a plethora of connected meters, sensors and other devices that will comprise living and working environments in the future. In order to meet the requirements of all these applications, mobile networks will have to accommodate both human-type communications (HTC) and machine-type communications (MTC) . MTC capabilities of existing predominantly HTC-oriented cellular technologies will be enhanced through the inclusion of non-cellular short-range and low-power wide area technologies under the 5G umbrella in the form of capillary networks. This will result in a multi-technology, high-interference radio environment such as can be found in some frequency bands in dense urban centers. In Figure , spectrum measurements in the shared 200 kHz unlicensed sub-1 GHz band over the course of 24 hours show signals from IEEE 802.15.4-based, LoRa and Sigfox networks as well as a number of unidentifiable proprietary technologies.

Given the increasing complexity of wireless technologies, evaluation and testing of new devices, protocols or solutions in realistic environments is desirable but tedious and expensive. Although complicated and expensive to build and maintain, distributed and heterogeneous wireless testbeds offer the most realistic experimental conditions . While the cost of adding additional hardware functionality to support emerging technologies is difficult to overcome, the cost of developing software for the increased variety and number of wireless devices can be significantly reduced using continuous integration (CI) principles and tools . In particular, open source software projects have benefited significantly from adopting CI principles that significantly reduce the workload of the core maintainers. Contributed code can be automatically tested for syntax, style and errors using sophisticated toolchains that automate previously manual development practices. Studies have shown that the CI methodology significantly increases productivity and code quality .

A recently proposed software-centric systematic testing methodology for mobile networks incorporates unit testing, integration testing and system testing, which are central to software development and represent an important part of the CI process. Testbed support enabling such a methodology can be economically implemented using modern tools, such as software containers, that can be automatically deployed to test new code.

The main contribution of this paper is COINS, a framework for ContinuOus IntegratioN in wirelesS technology development that can be used by developers to enable conventional CI practices on wireless infrastructure. To illustrate and facilitate discussion on the framework capabilities, we also provide a reference architecture and implementation. In particular, the existing LOG-a-TEC testbed for controlled testing of emerging multi-technology 5G capillary MTC networks was upgraded to enable CI. LOG-a-TEC is a combined outdoor and indoor heterogeneous wireless testbed for experimentation with sensor networks and machine-type communication[^1].

Two further contributions include (i) extending an existing qualitative methodology for comparing similar frameworks; and (ii) identification of and discussion on open challenges for wider adoption of CI practices in the development of wireless technologies.

The paper is structured as follows. Section identifies the main functionalities of a wireless experimentation testbed required by a software developer and explains how the CI process can be utilized in handling them. Section details the reference architecture and implementation of COINS. Section discusses the deployment of the new wireless experimentation technology as an extension of the existing LOG-a-TEC testbed and lists examples of supported experiments, while Section compares COINS with related work. Section identifies and discusses open challenges. Finally, Section summarizes the paper.

# Continuous integration process on wireless testbeds

There are four main functionalities that wireless network developers and testers need from a testbed environment:

-   Simple, low-effort integration of software and hardware components that require testing; the priority is to minimize the overhead that testbed usage introduces into the technology development cycle.

-   Controllable generation of traffic in the network under test, such as on-demand packet transmission according to custom patterns and distributions, to test new network protocols under realistic traffic conditions.

-   Controlled emulation of a radio operating environment using various technologies in order to identify potential coexistence issues between new and legacy technologies as early in the development cycle as possible.

-   Monitoring and analysis of network and radio spectrum metrics.

The introduction of the CI process into the wireless technology development proposed in this paper significantly lowers the barriers to the first and last of these functionalities while relying on the other two to successfully carry out the CI tests. According to existing CI practices , the proposed process should properly interface three main entities depicted in Figure a:

1.  *Repository* that holds the code, configuration files, tests and documentation, i.e., everything needed for the automated build. Generally, the repository is implemented using a version control system such as Git or Mercurial.

2.  *Automation system* that performs all the instantiations and builds and transfers software components. Conceptually, it does for CI what the production line does for a factory. Its implementation typically includes a large set of tools and components. In modern software development, the most efficient approach is to use web and containerization technologies.

3.  *Testbed device* is used to run the tests and experiments. It can be off-the-shelf equipment, a custom built device or a hybrid. The first approach is more suitable for improving existing technologies, while the second is for developing new ones. For the sake of generality, we assume in the following that a testbed device comprises separate functional blocks represented by the infrastructure node for device management and the target node for test execution.

As depicted in Figure a, the CI process begins in step 1 when a new testbed device automatically registers itself in the automation system and becomes available for testing. In step 2, the developer sends new code (pushes a commit) to the repository, initiating step 3, which consists of triggering the automation system. In step 4, the automation system pulls the changes from the repository and executes the deployment scripts specifying the testbed devices required by the test case. Step 5 deploys the changes from the repository to the testbed device and instructs it to start the build process. In step 6, the testbed device initiates the build process on the infrastructure node. The build system is contained inside a container; therefore, the infrastructure part of the testbed device stays unaffected by any build failures or run-away processes. If the building process succeeds, the infrastructure node in step 7 flashes the target node with freshly built firmware and in step 8 executes the test process. In step 9, the infrastructure node obtains the results from the target node and examines them. In step 10, the automation system updates the repository with the result and communicates the result to the developer, typically via email.

# COINS architecture and implementation

The architecture of COINS is illustrated in Figure b and further explained in this section. It consists of the three main architectural components identified and described in Section , which are mapped on the elements of the existing LOG-a-TEC infrastructure . The COINS architecture is implementation independent with well defined modules which interact through specified interfaces for which we provide a reference implementation although another implementation using a set of different tools is certainly possible.

For the reference implementation of COINS[^2], we chose widely used open source tools. In this way, COINS benefits from and contributes to the work of established communities and is more likely to be adopted. In particular, the container system is based on Docker[^3], which we use to package software for distribution to nodes in the testbed. The CI hook service adds support for the GitHub WebHook[^4] client. We use Ansible[^5] to describe and automate tasks. Ansible is a popular automation engine for configuration management and software deployment. The networked resource monitoring tool Munin[^6] serves for testbed device monitoring. We also implemented a custom node registry system and released it as free software under the AGPL license.

## Testbed devices

The LOG-a-TEC testbed devices can use several wireless technologies, as depicted on the left side of Figure b. Architecturally, we consider hybrid testbed devices with two separate functional blocks represented by the infrastructure node and the target node. The main reason for the separation of the testbed device into two functional blocks is to have a generic infrastructure node that can be combined with various target nodes to support a range of heterogeneous experiments with a unified way of management and reconfiguration. Another reason is to achieve hardware separation between robust and stable testbed management functionality on the infrastructure node and the experimental code, which is prone to errors and crashes on the target node. The infrastructure node can usually also support containerization, while the target node is often based on an embedded system with limited hardware resources.

In our implementation, the infrastructure node is a custom-designed single board computer based on the BeagleCore module running the Debian GNU/Linux operating system. It features wired Ethernet and Wi-Fi for infrastructure connectivity, an SD card slot for storage expansion, VESNA and USB interfaces for modular extensibility and a development/debug process. The target node is a custom VESNA (VErsatile platform for Sensor Network Applications) device with an ARM Cortex-M3 microcontroller application module and dedicated experimentation transceivers. VESNA can run a dedicated OS (e.g., Contiki) or custom firmware. The integration of the two is done via custom hardware design and the adoption of selected interfaces and protocols, thus going beyond similar approaches such as . Both nodes are interconnected via application and development interfaces, providing the exchange of application data as well as remote low-level application debugging. The reliable development interface is based on JTAG and the application interface on protocols such as LCSP (Light-weight Client Server Protocol; an HTTP-inspired protocol ) running on top of a serial interface.

## Repository

The repository support was added to the system to abstract the underlying complexity and enable the CI process described in Section . It includes (i) the target test that has to be executed on the testbed device; (ii) the test control that runs on the management server and communicates with the testbed device; (iii) the developer's code that has to be tested; (iv) dependencies such as software libraries; (v) the container file, which describes how to package the target test, test control, dependencies and developer's code into a container; and (vi) the deployment configuration, which specifies where and how to build and run the container.

In our implementation, the target test and test control are realized using C and Python, but in general, the proposed framework is programming language independent. The container support is implemented through Docker, while the deployment is automated by Ansible. The repository that contains everything needed for the completely automated build is realized through a GitHub repository. Examples of more advanced tests and experiments are summarized in Section .

## Infrastructure management and build automation system

The infrastructure management and build automation system is a complex setup composed of self-sufficient systems packaged in a Docker container following the micro-service architecture approach. It contains (i) the node registry, which has a list with details about all testbed devices; (ii) the node monitoring service, which is in charge of monitoring the activity and health of the testbed devices; (iii) the orchestration and configuration, which performs regular maintenance and reconfiguration of the testbed; (iv) the CI hook service, which waits for signals of changes in the repository and initiates the build process; and (v) the container file, which describes how to package all the previous components in a container.

There are several periodic background jobs running on the management system that perform house-keeping tasks:

-   building the Ansible host file from the node registry database,

-   updating the information on the availability of individual nodes using the Ansible ping command,

-   stand-alone testbed device monitoring based on Munin, and

-   CI Hook service triggering the build process.

Externally, the management server exposes a user-facing web interface and an HTTP REST API . The infrastructure management and build automation system allows browsing through clusters of individual devices and nodes in the database, positioning devices on a map, visualizing reports, monitoring device activity, etc. Thus, it gives an instant overview of the state of the testbed at any time. Warnings can be defined to be triggered when specific monitored variables exceed a given threshold, such as running out of system memory or storage space. In such cases, the system sends an email to the system administrator. Depicting the testbed devices on a map is particularly useful for the selection of devices for different wireless test cases. Programmatic access to the management system is performed through the REST API; i.e., the nodes use the REST API to perform automatic registration and update changes in the configuration.

# Testbed deployment and supported functionality

In order to enable testing and experimentation with the COINS framework in a dense and heterogeneous 5G capillary radio environment, we upgraded the combined outdoor/indoor LOG-a-TEC testbed at the JSI campus with the necessary tools and services described in Section . Table summarizes the number and type of testbed devices in an indoor and outdoor environment, whereas the layout of devices in the outdoor environment is depicted in Figure .

```latex
\begin{tabular}{llll}
        \bfseries Node type & \bfseries Radio & \bfseries Frequency & \bfseries No. of dev. \\
        [1ex] \hline \\ [-1ex]
        \multirow{2}*{Target node SRD A}    & Atmel AT86RF212   & 868 MHz   & \multirow{2}*{21 outdoor} \\
                    & TI CC2500     & 2.4 GHz   & \\
        [1ex] \hline \\ [-1ex]
        \multirow{2}*{Target node SRD B}    & TI CC1101     & 868 MHz   & \multirow{2}*{21 outdoor} \\
                    & Atmel AT86RF231   & 2.4 GHz   & \\
        [1ex] \hline \\ [-1ex]
        \multirow{2}*{Target node LPWA} & \multirow{2}*{Semtech SX1272} & \multirow{2}*{868 MHz}    & 3 outdoor \\
        & & & 1 indoor \\
        [1ex] \hline \\ [-1ex]
        \multirow{2}*{Target node UWB}  & \multirow{2}*{DecaWave DW1000}    & 3.5 - 6.5 GHz     & 11 outdoor \\
                    &                   & (6 channels)      & 20 indoor \\
        [1ex] \hline \\ [-1ex]
        Target node     & \multirow{2}*{NXP TDA18219HN}         & \multirow{2}*{470 - 866 MHz}      & \multirow{2}*{2 outdoor} \\
        UHF spec. sens. \\
        [1ex] \hline \\ [-1ex]
        Infrastructure      & \multirow{2}*{TI WL1837}      & \multirow{2}*{5 GHz}  & 58 outdoor \\
        node            &                   &           & 21 indoor \\
    \end{tabular}
```

The outdoor part of the wireless testbed is located in and around a park area of 55 m by 60 m. Testbed devices are placed on light poles 3.5 m above the ground and on the surrounding buildings at heights from 2.0 to 9.3 m. To enable experiments in indoor as well as in mixed indoor/outdoor scenarios, the testbed is extended in an indoor environment with additional 20 UWB devices and one LPWA device. These are deployed on the second and third floors of the 28.4 m by 16.6 m building. In addition, several portable devices of all types described in Section are at disposal for analyzing the use cases incorporating device mobility or further network densification.

The diversity of types, the relatively high number and density of testbed devices, and the possibility of hosting multiple radio technologies in most devices to transmit/receive/sense at the same time at the same location enables experimentation scenarios in a very dense and heterogeneous wireless environment. The various types of experiments include: (i) investigation of assisted and non-assisted indoor/outdoor radio localization algorithms; (ii) spectrum sensing for dynamic spectrum access and cognitive radio networking; (iii) link and/or network performance evaluation in a controlled interference and/or heterogeneous radio environment; (iv) design and evaluation of new custom-developed wireless protocols and their benchmarking against standardized protocols; and (v) performance evaluation and controlled piloting of prototype IoT devices.

In order to support such diverse testing capabilities, the repository comprises a set of testbed libraries and tools. They enable: out of the box localization; Line-of-Sight/non-Line-of-Sight link classification; link quality estimation (LQE), classification and prediction; custom protocol stack development using modular components; spectrum sensing by energy detection, cyclostationary detection and eigenvalue- and covariance-based methods; signal generation including wireless microphone profiles; and game theoretic power allocation. These software components and modules can be incorporated in custom developed algorithms or adapted to the particular use case and subsequently used as a testing benchmark. They can also be used for setting the required testing environment. In particular, for evaluation and comparison of spectrum sensing algorithms and other wireless technologies or algorithms, the testbed enables generating controlled interference. Furthermore, localization and distributed spectrum sensing can be used to build geolocation spectrum occupancy databases, which in turn can be used in investigation of primary/secondary spectrum access regimes.

A minimum test case for a wireless communication system includes software for three entities: the transmitter, the receiver and the test controller. The first two will run on the target node within the testbed device, while the third runs on the management server. The test controller instructs one testbed device to go in receive mode and another testbed device to transmit the test sequence. When the first testbed device receives the test sequence, it reports the results back to the test controller, which asserts that the received sequence matches the transmitted one. If it does, the test succeeds; otherwise, it fails i.e., `assertEqual(nodeTX.data, nodeRX.data)`. The described testing approach was implemented in our testbed for the purpose of LoRa firmware development. Especially in industry, there are situations where multiple developers work on the same code-base simultaneously. In such environments it can happen that a fix from one developer breaks a feature of another. By introducing automated testing on a real testbed, the integration problems get discovered and fixed early in the development process.

# Comparison with related work

While we did not find any other framework addressing wireless experimentation with CI support, we have identified three existing frameworks that are to some extent similar to COINS, i.e., the control and Management Framework (OMF) , the Network Implementation Testbed Laboratory (NITOS) and the Berlin Open Wireless Lab (BOWL) . We performed a feature-wise comparison of these frameworks and summarized it in Table .

In the comparison, we considered an extensive list of features that can be used to compare the experimentation systems . The list of features, however, lacks the CI-specific properties that our work focuses on; therefore, we extended the list with four core CI concepts :

-   Existence of a single source repository that contains everything needed for the completely automated build process.

-   Support for completely automated build process executed on each commit.

-   Self-tests included inside the repository, which run on each build.

-   Fast build process so that each commit can be built and tested.

| **Features**         | **OMF** | **NITOS** |      **BOWL**      |   **COINS**    |
|:---------------------|:-------:|:---------:|:------------------:|:--------------:|
| Description Language |  OEDL   |   REST    |        UCI         | Ansible/Docker |
| Experiment Types     |   Any   |    Any    |      Wireless      |    Wireless    |
| Operation            | Testbed |  Testbed  | Production/Testbed |    Testbed     |
| Interoperability     |         |           |                    |                |
| Reproducibility      |         |           |                    |                |
| Fault Tolerance      |         |           |                    |                |
| Debugging            |         |           |                    |                |
| Monitoring           |         |           |                    |                |
| Data Management      |         |           |                    |                |
| Source Repository    |         |           |                    |                |
| Automated Build      |         |           |                    |                |
| Self Testing         |         |           |                    |                |
| Fast Build           |         |           |                    |                |

Framework comparison.

For each feature presented in Table , we explain how we addressed the problem compared to the other frameworks. The last four features are CI specific and are realized only in COINS.

#### Description Language

The description language is the language used to describe how and where to build and run the experiment. In OMF, the experiment is described using OMF Experiment Description Language. The NITOS testbed is abstracted through a REST API. The BOWL system is based on the OpenWRT Unified Configuration Interface (UCI) scripts. In COINS we did not want to introduce a new language specifically for this purpose. Hence we used well-known tools that many developers already use in their development practices. In order to define what tests to run on which node and how to obtain test results we use the Ansible Playbook, a language for describing software deployment steps. For building our containers we use Dockerfile, a text format used by Docker to define container build procedures.

#### Experiment Types

The type of experiment describes the platform where the experiment is going to run. OMF and NITOS can support wired or wireless experiments, whereas BOWL is limited to WiFi-based wireless systems. COINS focuses on wireless systems and is suitable also for restricted capability devices. While OMF, NITOS and COINS are suitable for developing completely new as well as optimizing existing wireless solutions, BOWL is primarily suitable for optimizing existing solutions.

#### Operation

The operation refers to the ability of the testbed to also serve as a production system. OMF and NITOS are exclusively testbed frameworks, while BOWL is designed to work simultaneously as a testbed and production WiFi access network. In this respect, COINS is also exclusively a testbed framework.

#### Interoperability

Framework interoperability covers an infrastructure independent design, support for existing infrastructure services, resource discovery and software interoperability. OMF and NITOS address interoperability issues, while BOWL UCI configuration scripts are implementation specific. OMF can be extended, and NITOS specifies a testbed API. COINS addresses interoperability by offering a modular and implementation independent design.

#### Reproducibility

Reproducibility refers to the ability to automatically run the same experiments multiple times. OMF enables complete reproducibility, while NITOS offers limited reproducibility since it implements automated node selection based on availability. BOWL also addresses the problem of making experiments reproducible. COINS supports reproducibility of test cases by creating tags inside the repository and manually choosing the devices used in the test. While declaratively all frameworks address reproducibility, we argue that it may be severely affected in case of missing repository support on the framework level and challenged by the ability to control the conditions in the physical environment.

#### Fault Tolerance

Fault tolerance addresses recovery capabilities when the experiment does not execute as expected. OMF and NITOS frameworks are designed to be fault tolerant. In OMF, there is a Management Plane responsible for rebooting a resource and setting up particular disk images on PC-based resources. NITOS uses an automatic monitoring tool to diagnose a malfunction in one of the resources. BOWL is not fully fault tolerant as it also serves as a production WiFi network for actual users. COINS is inherently fault tolerant due to its architecture design having separated infrastructure and target node as well as separate development and application interfaces.

#### Debugging

The debugging process addresses the ability to identify problems and their causes during the experiment execution. All compared frameworks are designed to produce extensive log files during the build process. COINS also collects extensive debug log during the execution. Furthermore the debug log is automatically emailed to the developer in case the execution fails.

#### Monitoring

Monitoring is a feature that enables the developer or a maintainer to oversee the status of the system at any time. All compared frameworks have a system for monitoring. COINS uses the Ansible ping method, which works on top of SSH to identify which nodes are accessible. Furthermore, it integrates a standalone open-source monitoring tool Munin as part of the management system.

#### Data Management

Data management is responsible for the distribution and collection of data. OMF uses the Measurement Collection system based on a relational database. NITOS and BOWL do not prescribe any dedicated data collection mechanism on the framework level. COINS uses the same system for build data management as for debugging; i.e., each build produces the result files along with the debug files in the build repository.

#### Source Repository

The COINS framework is designed to keep everything needed for a completely automated build inside a single version-controlled source repository. This is in line with one of the fundamental CI concepts.

#### Automated Build

The automated build is another core CI concept. It ensures that a new code committed by the developer does not break the process, which is particularly important in CI. We showed that COINS defines the automated build system, which deploys the modified source code to the appropriate targets where it is built and tested.

#### Self Testing

Tests of each committed change are extremely important for the CI process. In COINS, all the tests are included inside the repository. The test process is executed automatically on each commit, and the result is communicated to the developer.

####  Fast Build

The last core CI concept requires a fast build process so that the developer does not have to wait long for the test results. To achieve fast builds, COINS relies on the container technology and incremental changes delivery to the testbed device. Therefore, only the changes are rebuilt within the container, thus making the build process fast.

# Discussion on open Challenges

For wider adoption of CI practices in the development and testing of wireless technologies, we identified three groups of challenges. The first group revolves around the constraints posed by the physical environment; the second one is concerned with process automation; and the third group is related to the lifecycle management of the wireless product.

## Challenges related to the physical environment

A CI system for wireless networks comprised of distributed hardware has to take into account that there may be unpredictable radio interference on the physical transmission medium or failures in hardware, and it must be ensured that the test reports to the developer are immune to them. In other words, from the test reports, it should be clear if the test failed because of software issues or because of conditions in the physical environment.

#### Reproducibility

A possible approach for ensuring reproducibility of the test results in the unpredictable wireless environment is for the CI framework to use spectrum sensing to identify a suitable test channel and thus avoid interference. Additionally, using controlled interference in test cases can sometimes prove useful to verify the interference specific performance of a particular wireless algorithm.

#### Reliability

A possible approach for ensuring reliability is to have redundancy in terms of distributed hardware and perform the same tests on more disjoint subsets of testbed devices.

## Challenges related to the process automation

With the implementation of the proposed framework, it is possible to realize automated system tests. However, while the CI system and process are fully automated, the detection of errors caused by the physical environment is currently manual. This is an issue specific to wireless environments and has yet to be rigorously addressed. Additionally, the developer of the test currently needs to manually select one or more testbed devices to be used in the test.

#### Automated tests

A possible approach for enabling fully automated tests is to design and develop a mechanism that would analyze spectrum sensing data in correlation with the tests and repeat the test if it failed because of a temporally jammed channel. This process could be made transparent to the developer.

#### Automated resource allocation while retaining test results reproducibility

A possible approach for enabling automated resource allocation is to develop a mechanism that would automatize the resource (i.e., testbed device) assignment for individual test execution. It would be interesting to see if the concurrent tests could be somehow categorized and how this information could be utilized by an algorithm for automatic selection of testbed devices to optimize test execution and testbed utilization while keeping test results reproducible as well.

## Challenges related to the lifecycle management of the wireless product

Continuous software development should also consider continuous delivery of the builds in production. The proposed framework is focused on developing quality software also for restricted capability devices. However, additional consideration could be given to extending the framework with continuous delivery capabilities to solve the issue of wireless devices, which are not maintained after entering the market. Continuous delivery could also solve security issues in devices on the market. Currently, every electronics producer addresses the problem of continuous delivery separately; however, it would be beneficial to find some standard method.

# Summary

In this paper, we introduced continuous integration (CI) practices for wireless network development. We proposed the COINS framework, which simplifies the code development process on real wireless testbeds and enables frequent code changes. We then showed that the proposed framework can be implemented on a testbed composed of state-of-the art hardware modules suitable for 5G capillary wireless MTC experimentation. We realized a reference implementation of COINS by integrating existing open-source components and cloud services into a fully functional system supporting CI for wireless networks. We performed a qualitative comparison of the proposed framework with three related existing frameworks. The comparison uses a list of features already used in the community for comparing experiment management tools. However, to fully illustrate the advantage of the proposed framework, we had to extend the available list of features with four CI-related ones. Finally, we identified and discussed open challenges.

# Acknowledgment

This work was partly funded by the Slovenian Research Agency (Grant no. P2-0016) and the European Community under the H2020 eWINE (Grant no. 688 116), WiSHFUL (Grant no. 645 274) and Fed4FIRE+ (Grant no. 732 638) projects.

# Biographies

MATEVŽ VUČNIK received his B.Sc. in electrical engineering from the University Of Ljubljana, Faculty of Electrical Engineering in 2010. He is currently a PhD candidate at the Jožef Stefan International Postgraduate School and works as a senior research assistant in the Department for Communication Systems at the Jožef Stefan Institute. His main interests include distributed software development and Web technologies, focusing on communication protocols and IoT integration. He has actively collaborated on several EU projects.

TOMAŽ ŠOLC studied electronics at the Faculty of Electrical Engineering, University of Ljubljana and received his B.Sc. in 2007. He is currently a senior research assistant at the Department of Communication Systems, Jožef Stefan Institute and a Ph.D. student at the Jožef Stefan International Postgraduate School. His work involves hardware design, measurements and embedded software development as well as research in the field of spectrum sensing. He participated in various national and EU projects.

URBAN GREGORC received a Master degree in Electrical Engineering from the University of Ljubljana, Faculty of Electrical Engineering in 2016. After graduation, he worked at the Department of Communication Systems at the Jožef Stefan Institute, Ljubljana as a research assistant and in 2018 moved to a startup ComSensus. His research is focused on Low Power Wide Area Network technologies for the Internet of Things.

ANDREJ HROVAT received his B.Sc. in 2004. He obtained a Ph.D. degree in Electrical Engineering from the Jožef Stefan International Postgraduate School, Slovenia in 2011. He is currently a researcher in the Department of Communication Systems of the Jožef Stefan Institute and an assistant professor at the Jožef Stefan International Postgraduate School. His research interests include radio signal propagation, channel modelling, terrestrial and satellite fixed and mobile wireless communications, radio signal measurements and emergency communications.

KLEMEN BREGAR received a degree in electrical engineering from the University Of Ljubljana, Faculty of Electrical Engineering in 2013. Since 2013, he has worked as a young researcher and Ph.D. candidate at the Department of Communication Systems at Jožef Stefan Institute, Ljubljana. His research interests include indoor localization, ultra-wide band communication systems, biomedical wireless sensor networks and machine learning.

MIHA SMOLNIKAR received his BSc from the University of Ljubljana and has more than 10 years of experience in networked embedded systems and wireless communications. His main interests include wireless sensor networks, energy harvesting, sensor technologies and the optimization of the protocol stack. In the area of IoT, he focuses on the development of products and services for the smart home, smart city and smart grid segments. He actively participates in several EU projects.

MIHAEL MOHORČIČ is head of the Department of Communication Systems and Scientific Counselor at the Jožef Stefan Institute. He is also associate professor at the Jožef Stefan International Postgraduate School. His research experience include development and performance evaluation of network protocols for mobile and wireless communication systems, cognitive radio networks, wireless sensor networks, dynamic composition of communication services and wireless experimental testbeds. He has participated in many EU projects. He is a Senior Member of IEEE.

CAROLINA FORTUNA received her B.Sc. in 2006, Ph.D. in 2013 and was a postdoctoral research associate at IBCN, Ghent University, 2014-2015. Currently, she is research fellow at the Department of Communication Systems, Jožef Stefan Institute and an assistant at the Jožef Stefan International Postgraduate School. Her research is interdisciplinary, focusing on data and knowledge driven modelling of communication and sensor systems. She has participated in EU and industrial projects with Bloomberg LP and Siemens PSE.

[^1]: LOG-a-TEC, <http://log-a-tec.eu/>

[^2]: COINS, <a href=" https://github.com/sensorlab/SensorManagementSystem" class="uri"> https://github.com/sensorlab/SensorManagementSystem</a>

[^3]: Docker, <https://www.docker.com>

[^4]: GitHub, <https://developer.github.com/webhooks>

[^5]: Ansible, <https://www.ansible.com>

[^6]: Munin, <http://munin-monitoring.org>


abstract: Hidden U(1) gauge symmetries are common to many extensions of the Standard Model proposed to explain dark matter. The hidden gauge vector bosons of such extensions may mix kinetically with Standard Model photons, providing a means for electromagnetic power to pass through conducting barriers. The ADMX detector was used to search for hidden vector bosons originating in an emitter cavity driven with microwave power. We exclude hidden vector bosons with kinetic couplings $\chi>3.48\times 10^{-8}$ for masses less than $3~\mueV$. This limit represents an improvement of more than two orders of magnitude in sensitivity relative to previous cavity experiments.
author: A. Wagner, G. Rybka, M. Hotz, L. J Rosenberg; S.J. Asztalos[^1], G. Carosi, C. Hagmann, D. Kinion and K. van Bibber[^2]; J. Hoskins, C. Martin, P. Sikivie and D.B. Tanner; R. Bradley; J. Clarke
date: 2024-10-02
title: A Search for Hidden Sector Photons with ADMX

The Standard Model (SM) of particle physics explains the electromagnetic, weak and strong interactions of the leptons and quarks that comprise the visible universe. Extensions of the SM allow for a so-called hidden sector of particles; these particles may be so extremely weakly interacting with the known SM particles that they have little effect on observations. This hidden sector has been suggested as a model for dark matter in the universe  where the particles of the hidden sector are charged under a hidden U(1) gauge symmetry. The hidden U(1) vector boson may mix kinetically  with the SM photon providing a weak coupling of the hidden sector to the SM. In those models where the kinetic mixing provides the only coupling of the hidden sector to the SM, the hidden vector boson is referred to as a paraphoton. Paraphotons with a mass between $50~\MeV$ and $1~\GeV$ could provide a mechanism to explain the excess positron fraction observed by PAMELA  while light paraphotons in the mass range $1~\mueV - 1~\meV$ have been shown to produce a hidden cosmic microwave background (hCMB) .

Paraphoton-photon mixing may be detected directly via the measurement of transmitted power between shielded microwave cavities . We use the Axion Dark Matter eXperiment (ADMX) detector in this way to search for evidence of hidden sector kinetic mixing. This experiment is sensitive to paraphotons with $\Mpara\approx 3\mueV$, the low end of the paraphoton mass range that may produce a hCMB.

In general the Lagrangian describing the SM and a hidden sector U(1) with kinetic coupling $\chi$ is

$$\begin{gathered}
\label{eq:L}
{\cal L} = -\frac{1}{4}F^{\mu\nu}F_{\mu\nu} -{1\over 4}B^{\mu\nu}B_{\mu\nu} \\
-{1\over 2}\chi F^{\mu\nu}B_{\mu\nu} +{1\over 2}\Mpara^2 B^{\mu}B_{\mu},
\end{gathered}$$

where $B^\mu$ is the vector potential of the paraphoton, $\Mpara$ is the mass of the paraphoton, and $F^{\mu\nu}$ and $B^{\mu\nu}$ are the electromagnetic and hidden sector field strength tensors, respectively. The kinetic coupling described by the $F^{\mu\nu}B_{\mu\nu}$ term in  results in a mixing of photons with paraphotons. The origin of this mixing is made explicit with the change of basis $B^\mu\rightarrow \tilde{B}^\mu - \chi A^\mu$ to diagonalize the kinetic terms in ,

$$\begin{gathered}
\label{eq:Lint}
{\cal L} = -\frac{1}{4}F^{\mu\nu}F_{\mu\nu} -{1\over 4}\tilde{B}^{\mu\nu}\tilde{B}_{\mu\nu} \\
+{1\over 2}\Mpara^2(\tilde{B}^{\mu}\tilde{B}_{\mu} - 2\chi\tilde{B}^\mu A_\mu + \chi^2 A^\mu A_\mu).
\end{gathered}$$

Here, we have absorbed the renormalization of the electromagnetic gauge coupling $e^2 \rightarrow {e^2/(1 - \chi^2)}$. The mass term in  coupling the electromagnetic and hidden sector vector potentials gives rise to a mixing angle between hidden sector interaction and mass eigenstates. This provides the mechanism for paraphoton-photon oscillations  analogous to that of neutrino flavors.

The coupling of the paraphoton to the photon modifies the electromagnetic potential with an additional Yukawa term . Constraints from the non-observation of deviations from Coulombs law  and the black-body spectrum of the cosmic microwave background  exclude couplings $\chi>10^{-7}$ in the paraphoton mass range $10^{-9}<\Mpara<10^{-4}~\rm{eV}$.

Microwave cavities can be used as sensitive searches for paraphotons via the resonant enhancement of photons mixing into paraphotons in an emitter cavity and the observation of transmitted power in a detector cavity resulting from the paraphotons oscillating back into photons . Such cavities typically cover the paraphoton mass range $10^{-6}<\Mpara<10^{-3}~\rm{eV}$. The power observed in the detector cavity $\Pdet$, for an emitter cavity supplied with power $\Pem$, is given by

$$\label{eq:pow}
\Pdet=\chi^4Q_{\rm det}Q_{\rm em}G^2\bigg(\frac{\Mpara}{\omega}\bigg)^8\Pem,$$

where the emitter cavity is driven with frequency $\nu=\omega/2\pi$, and $Q_{\rm det},~Q_{\rm em}$ are the quality factors of the detector and emitter cavities respectively. These factors account for the resonant enhancement of photon paraphoton mixing probabilities. They are typically of order $10^4$ for cavities used by ADMX. The dimensionless coefficient $G$ accounts for the interference of photons in the detector cavity as paraphotons mix back in different locations with different relative phases and is given by

$$\label{eq:G}
G=\omega^2 \int_{V_{\rm em}}\int_{V_{\rm det}}{e^{ik|\bold{x-y}|}\over {4\pi|\bold{x-y}|}}\vec{A}_{\rm em}(\bold{y})\cdot\vec{A}_{\rm det}(\bold{x})d^3\bold{y}d^3\bold{x}.$$

Here, $k$ is the wave number of the paraphoton, $k^2 = \omega^2 - \Mpara^2$, and the vector potentials of the EM normal modes, $\vec{A}_{\rm em},~\vec{A}_{\rm det}$, obey the normalization condition

$$\label{eq:A_norm}
\int_{V}|\vec{A}|^2(\bold{x})d^3\bold{x} = 1,$$

where we use the Coulomb gauge $\vec{\nabla}\cdot\vec{A}=0$. The magnitude of $G$ is shown in  as a function of the relative paraphoton wave number, $k/\omega$, for the arrangement of cavities in this experiment, shown in .

ADMX employs a tunable microwave cavity immersed in a 7.5 T magnetic field. It was originally built to detect the conversion of dark matter axions into photons. The ADMX detector is described in detail in Ref. , and recent results of the axion search are given in Ref. .

Microwave power from paraphoton mixing can be detected in an identical manner to that from axion conversions. Electromagnetic power in the ADMX cavity excites the  normal mode which is read out by an antenna, amplified, mixed, and Fourier transformed with a frequency resolution of 125 Hz. The sensitivity of power measurements is limited by black body noise from the cavity and the intrinsic noise of the amplifier. It is thus important to operate at as low a temperature as possible. The present experiment maintains the cavity at a physical temperature of 2 K with pumped helium while the High electron mobility Field Effect Transistor (HFET) amplifiers  have a noise temperature of 4 K. This results in a total noise temperature of 6 K.

A significantly lower cavity temperature would reduce the cavity black body noise, but not the HFET noise. The use of Superconducting Quantum Interference Device (SQUID) amplifiers, however, with a noise temperature approaching the quantum limit at low temperatures , reduces the system noise temperature substantially, for example, to a projected 200 mK or less for an operating temperature of 100 mK. The SQUID amplifiers used during normal ADMX operations were not operational for this experiment.

The signal-to-noise ratio for power detected in the ADMX cavity is set by the radiometer equation

$$\label{eq:SNR}
{S\over N}={\Pdet\over k_BT}\sqrt{t\over b},$$

where $T$ is the total noise temperature, $t$ is the time over which the power is measured and $b$ is the frequency bandwidth resolution. A paraphoton signal would be observed as a peak above the noise spectrum in the detector cavity at the frequency of the driven emitter cavity. We find that the noise temperature of the system dominates any systematic effects of the receiver chain for integration times less than 30 days .

The paraphoton search procedure is to operate the ADMX detector cavity in the  mode, tuned via the placement of two 0.0508 m diameter rods, to a resonant frequency of 722.725 MHz. A second cavity is supplied with 150 mW of power in the same mode and frequency. Both the cavities and tuning rods are copper plated stainless steel cylinders. The cavities each measure 0.40 m in diameter and 0.926 m in height and are separated 3.30 m vertically and 1.35 m horizontally due to the presence of the detector cryostat. The quality factors of the cavities, $Q_{\rm det},Q_{\rm em}$ are $7\times 10^4$ and $1.4\times 10^4$ respectively at 722 MHz. A schematic of the ADMX paraphoton search is shown in .

We average the power in a 125 Hz wide bin around 722.725 MHz for a period of $7.2\times10^3$ s. We limit the excess power to $\Pdet=0.69\pm1.08\times10^{-23}~{\rm W}$ (). Since this is consistent with no excess power we set a limit on the paraphoton coupling $\chi$ as function of mass, $\Mpara$ using  and . We exclude paraphotons over a broad mass range. For masses near $3~\mueV$, we exclude kinetic couplings $\chi > 3.48\times 10^{-8}$ at the $95\%$ confidence level as shown in . We exclude couplings to lower paraphoton masses with reduced sensitivity to the coupling.

Feed-through of microwave power from the emitter to the receiver antennas is a potentially serious systematic effect that dramatically limited the sensitivity of microwave cavities to paraphoton mixing in a previous experiment . However; ADMX has been optimized to eliminate external microwave signals which could limit the sensitivity of the axion search. The ADMX cryostat functions as an excellent Faraday cage where all signal and sensor lines have been tested for leakage with an external microwave source. We have terminated all unnecessary lines and supplied the remainder with additional shielding. The HFET amplifiers supply 30 dB of amplification and are well shielded within the cyrostat. This greatly suppresses the effect of any feed-through from outside the cryostat. Consequently, the ADMX paraphoton search was not limited by the feed-through of microwave power.

The ADMX receiver will be upgraded to improve the experiment's sensitivity. This upgrade will employ SQUID amplifiers with noise temperatures close to the quantum limit and a dilution refrigerator cooled to 100 mK to reduce the total noise temperature to 200 mK. Improved frequency stability will allow bandwidth resolution of better than 0.02 Hz. The reinstallation of ADMX at the University of Washington will also enable the emitter and detector cavities to be positioned within 1.5 m so that $|G|\approx 1$. The power driving the emitter cavity was supplied by the source used for nominal ADMX operations; power levels in excess of 1 W will be available in the receiver upgrade. We estimate this upgrade will provide sensitivity to paraphotons with a kinetic coupling $\chi< 10^{-9}$ ().

In summary, we have used the ADMX detector in a search for hidden sector physics. We exclude paraphotons over a broad mass range. For paraphoton masses near $3~\mueV$ we exclude kinetic couplings $\chi > 3.48\times 10^{-8}$ at the $95\%$ confidence level. Lower paraphoton masses are excluded with reduced sensitivity to the kinetic coupling. Upgrades to ADMX will permit an improvement in this limit by more than an order of magnitude. This measurement is a significant improvement in microwave cavity limits and is comparable to the sensitivity of the limit set in Refs.  by tests of Coulomb's law. Finally, this experiment, together with the axion  and chameleon  searches performed by ADMX, demonstrates not only the extraordinary sensitivity of microwave cavity experiments but also the breadth of physics such experiments can expolore.

The ADMX collaboration gratefully acknowledges support by the U.S. Department of Energy, Office of High Energy Physics under contract numbers DE-FG02-96ER40956 (University of Washington), DE-AC52-07NA27344 (Lawrence Livermore National Laboratory), and DE-FG02-97ER41029 (University of Florida). Additional support was provided by Lawrence Livermore National Laboratory under the LDRD program. Development of the SQUID amplifier (JC) was supported by the Director, Office of Science, Office of Basic Energy Sciences, Materials Sciences and Engineering Division, of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.

[^1]: Currently at XIA LLC, 31057 Genstar Rd., Hayward CA, 94544.

[^2]: Currently at Naval Postgraduate School, Monterey CA, 93943.


**Step by Step to Peace in Syria**

Raphael Parens and Yaneer Bar Yam

*New England Complex Systems Institute*

210 Broadway Suite 101 Cambridge MA 02139, USA

(Dated February 9, 2016)

**The revolution and Civil War in Syria has led to substantial death and suffering, a massive refugee crisis, and growth of ISIS extremism and its terror attacks globally. Conflict between disparate groups is ongoing. Here we propose that interventions should be pursued to stop specific local conflicts, creating safe zones, that can be expanded gradually and serve as examples for achieving a comprehensive solution for safety, peace and stable local governance in Syria.**

Nearly five years since the deterioration of government control and order, stemming from the "Arab Spring" protests of 2011, Syria has become a "failed state" with widespread violence, the source of ISIS ideological violent extremism and terrorism, and the origin of massive population displacement and exodus. Often called a civil war, violence in Syria includes conflict between many parties: the government forces, ISIS, several distinct rebel groups with various affiliations, and many external parties including national actors (Saudi Arabia, Turkey, Iran, Russia and the US) and non-state actors (Hezbollah, Druze and Kurdish groups). Despite the existence of alliances on the ground, specifically in opposition to ISIS, these groups have different agendas and support different local groups and may interfere with, or even fight and kill, members of other groups.

Because of the complexity of the conflict there are various framings of the situation as a battle between the opposition and Assad, a battle against ISIS, a battle of proxies, etc. Here we consider a strategy of achieving peace by addressing local conflicts with local solutions prior to building a comprehensive solution to the entire set of conflicts at the national scale. This is a grass roots approach to building safety and security from the ground up. Such efforts already exist in some places. We provide a framework for this approach in terms of a validated scientific analysis of essential drivers of conflict rooted in ethnic geography. While some of those involved will consider our discussion narrowly in terms of whether we provide support for their cause, our objective is to address directly the suffering of populations through establishing robust local safety and governance that will be an immediate relief to local populations severely affected by the conflict. The safety of the populace should not and need not be a hostage for the national solution which may follow.

Complex ethnic geography is a central reason for the large number of different groups in conflict as this geography leads to local allegiances that do not aggregate at the national scale. The makeup, and even existence of several of these groups is a matter of debate, complicated by the possibility that some groups choose to hide their differences (e.g. Alawites and Nusaries) . These debates do not affect the overall framework of our discussion, although they may play a role in the subsequent consideration of local conditions. A well researched map of ethnic geography is shown in Fig. showing 13 ethnic groups that have local geographic areas in which they are the majority population. The largest of these ethnic groups are, in descending order by population size, Arab Sunni Muslims (the majority, 60% of the population), Arabic speaking Alawites, Arabic speaking Christians (Levantines), Kurds, Druze, Ismailis, Nusaries, and Imamis .

We have previously demonstrated that, where ethnic groups exist in geographic patches of 20-60 km in diameter, there is a high probability of conflict . We have also shown that establishing local autonomy through subnational boundaries (as is found in Switzerland) is a means of alleviating conflict . Our analysis suggests that providing some level of local autonomy to the ethnic groups would reduce the impetus for local conflict and could serve as a basis for peace and stability .

The objective of local autonomy is to enable local decision making to reflect values in determining ordinances, institutions and services to the local population. As is found in Switzerland, where there are multiple cantons of each of catholic and protestant dominated areas, there is no need to have all of the members of a group in a single governance structure. The objective rather, is to reduce friction that leads to conflict by enabling those with widely different values to impose their values on public spaces. We note that the establishment of local partially autonomous regions does not itself determine the structure of federal national governance, which would be the framework for national security and interactions among the local groups as in, for example, the federal government of Switzerland. Other countries have federal governance systems with various ways of balancing local and national governance roles, including the United States, where four levels of governance are present in municipal, county, state and federal systems.

In order to advance the understanding of the conflict and potential solutions in Syria, we have constructed a visualization identifying areas susceptible to violence. This study used a simplified methodology compared to previous studies . We superimposed circles of 20 km diameter on the ethnic map of Syria by visual inspection. Circles are placed where one ethnic population is the majority within the circle and others are the majority around the periphery. Fig. manifests why and where ethnic group affiliation, whether defined by language, religion, or both, plays a major role in the conflict in Syria.

Our analytic results for the locations of likely violence are consistent with the current state of conflict shown in Fig. . We see that ethnic violence in Syria can be divided into four main regions and an isolated patch: northeast, north, west, southwest, and Deir az-Zur.

*Northeast*—The northeast region is dominated by Kurd and Sunni Arab groups. Violence is particularly likely between them along the northern border with Turkey, where the two groups are interspersed in 20 km patches. In addition, some combination of Armenians, Assyrians/Syriacs/Chaldians, and Chechens are surrounded by Kurds and Sunnis in several patches in this area, perhaps also contributing to violence.

*North*—The north region contains areas of likely conflict between Kurds, Sunnis, Circassians, Imamis, Levantines, and Turkomans.

*West*—The west region contains the largest number of likely ethnic violence areas, as a collection of Kurds, Druze, Imamis, Circassians, Assyrians/Chaldians/Syriacs, Ismailis, Nusairis, Turkomans, and Arabs all live in ethnic pockets of the critical size.

*South*—The south region in the area of Damascus also contains a significant number of ethnic patches, with possible conflict between Levantines, Druze, Aramaeans, Assyrians/Chaldians/Syriacs, Circassians, and Kabardas.

*Deir az-Zur*—Finally, the city of Deir az-Zur contains an isolated ethnic patch, where a population of Armenians is surrounded by Arab populations.

According to our previous analysis, the establishment of boundaries between ethnic populations to provide partial local autonomy would increase stability and inhibit current and future conflict in Syria. Some suitable natural barriers already exist due to the topography of mountains, lakes, and rivers. However, many boundaries will have to be established through political borders or artificial barriers.

Under the current conditions of fragmentary control and multiple competing groups, implementing a national process to resolve all conflicts is difficult. The complexity of local conflicts, and the many parties involved, will be a barrier to any such comprehensive plan. We propose that a step-by-step bottom up strategy provides a useful and robust alternative to a national plan. Indeed, in this context, the natural scale of intervention is at the community level. Our analysis should serve as a motivation for local governance creation and maintenance in Syria, rather than a blueprint. The complexity of governance creation on the ground will require adaptation due to the specifics of local conditions.

In such a step-by-step strategy "safe zones" should be established. Efforts should be made to identify specific local areas of conflict for intervention to establish the safe zones, including local governance and sub-national political borders, potentially redrawing governorates (muhafazat), districts (manatiq), or subdistricts. Such interventions should recognize specific conditions of villages and urban neighborhoods, their values and capacity to provide for their own security, as well as the relationships they have with nearby groups and national or international groups. The complexity of these local conditions must be addressed through direct engagement with the local population as an integral part of the process of achieving security. In recent years new local governance structures have emerged in response to the disorder . These may serve as a basis for the robust and longer-term structures that are needed. Once established, a safe zone will allow the local population to rebuild and reestablish normal lives.

Among the challenges to be faced in achieving safe zones is establishing reliable international support for security, and where necessary, assigning ownership of geographically associated economic resources and religious/cultural site control. Identifying locations where the framework for peace, even if complex, can be more readily achieved is key to early progress and establishing precedent for later more difficult areas.

The greatest challenge in implementing this approach may be the psychological shift from seeing power as an absolute quality on a nation-state basis, and shifting to a perspective in which local power is balanced against the role of allegiances and against larger scales of power, including national. The current conflict is often seen as an irreconcilable power struggle between local groups and the nation-state, as embodied by its government. Here the objective is to show how local groups can coexist with national power. This coexistence is not about full autonomy but rather about a balance between local autonomy within relationships among the local, national and international groups.

By focusing on the local nature of interactions, the basis of community governance can be established. The structure of national power, in whatever form then occurs, becomes secondary as individual day-to-day existence can be primarily determined by the local authority and only secondarily with the national governance.

**Acknowledgements:**

We would like to thank Francisco Prieto-Castrillo, Nima Dehghani and Matthew Hardcastle for helpful comments on the manuscript and Nadim Shehadi, Nassim Taleb and Tom Garvey for helpful discussions.


abstract: We numerically investigate the intriguing effects produced by random percolative disorder in two-dimensional Josephson-junction arrays. By dynamic scaling analysis, we evaluate critical temperatures and critical exponents with high accuracy. It is observed that, with the introduction of site-diluted disorder, the Kosterlitz-Thouless phase transition is eliminated and evolves into a continuous transition with power-law divergent correlation length. Moreover, genuine depinning transition and creep motion are studied, evidence for distinct creep motion types is provided. Our results not only are in good agreement with the recent experimental findings, but also shed some light on the relevant phase transitions.
address:  $^{1}$ Department of Physics, Zhejiang University, Hangzhou 310027, P. R. China
  $^{2}$ Center for Statistical and Theoretical Condensed Matter Physics, Zhejiang Normal University, Jinhua 321004, P. R. China
author: Jian-Ping Lv$^{1}$, Huan Liu$^{1}$, and Qing-Hu Chen$^{2,1,\dag}$
date: 2024-10-01
title: Phase transition in site-diluted Josephson-junction arrays: A numerical study

# Introduction

Understanding the critical behavior of Josephson-junction arrays (JJA's) with various disorders is always a challenging issue and has been intensely studied in recent years –. However, the properties of different phases and various phase transitions are not well understood. Josephson-junction arrays gives an excellent realization to both two-dimensional (2D) $XY$ model and granular High-$T_{c}$ superconductors . As we know, the pure JJA's undergoes the celebrated Kosterlitz-Thouless (KT) phase transition from the superconducting state to the normal one, this transition is driven by the unbinding of thermally activated topological defects . When the disorder is introduced, the interplays among the repulsive vortex-vortex interaction, the periodic pinning potential caused by the discreteness of the arrays, and the defects produced by the disorder provide a rich physical picture.

In site-diluted JJA's, the crosses around the randomly selected sites are removed from the square lattice. Since it is a representative model for realizing the irregular JJA's systems, how the percolation influences the physical properties of JJA's has attracted considerable attention . Harris $et$ $al$ introduced random percolative disorder into Nb-Au-Nb proximity-coupled junctions, the current-voltage ($I$-$V$) characteristics were measured and the results demonstrated that the only difference of the phase transition compared with that in ideal JJA's system is the decrease of critical temperature, while the transition type still belongs to the KT one with the disorder strength spanning from $p=0.7$ to $p=1.0$ (here $1-p$ is the fraction of diluted sites) . However, in a recent experiment, Yun et al showed that the phase transition changes into a non-KT-type one when the disorder strength increases to a moderate value ($p=0.86$) . Therefore, the existence of the KT-type phase transition in site-diluted JJA's remains a topic of controversy, the nature of these phase transitions and the various phases is not clear.

On the other hand, much effort has been devoted to the zero-temperature depinning transition (ZTDT) and the related low-temperature creep motion (LTCM) both theoretically and numerically in a large variety of physical problems, such as charge density waves , random-field Ising model , and flux lines in type-II superconductors . Since the non-linear dynamic response is a striking problem, there is increasing interest in its properties and characteristics, especially in the flux lines of type-II superconductors . In a recent numerical study on the three-dimensional glass states of flux lines, Arrhenius creep motion was observed at a strong collective pinning, while the non-Arrhenius creep motion was demonstrated at a weak collective pinning .

In this work, we numerically study the finite-temperature phase transition (FTPT) in site-diluted JJA's at different percolative disorder strengths, the ZTDT and the LTCM are also investigated. The outline of this paper is as follows. Sec. describes the model and the numerical method briefly. In Sec. , we present the main results, where some discussions are also made. Sec. gives a short summary of the main conclusions.

# Model and simulation method

JJA's can be described by the 2D $XY$ model on a simple square lattice, the Hamiltonian reads $$H=-\sum\limits_{ <i,j>}J_{ij}\cos(\phi_{i}-\phi_{j}-A_{ij}),$$ where the sum is over all nearest neighboring pairs on a 2D square lattice, $J_{ij}$ denotes the strength of Josephson coupling between site i and site j, $\phi_{i}$ specifies the phase of the superconducting order parameter on site i, and $A_{ij}=(2\pi/\Phi_{0})\int \mathbf{A} \cdot d \mathbf{l}$ is the integral of magnetic vector potential from site i to site j, $\Phi_{0}$ denotes the flux quantum. The direct sum of $A_{ij}$ around an elementary plaquette is $2\pi f$, with $f$ the magnetic flux penetrating each plaquette produced by the uniformly applied field, which is measured in unit of $\Phi_{0}$. $f=0$ and $f=2/5$ are the focuses of this paper. The system sizes are selected as $128\times128$ for $f=0$ and $100\times100$ for $f=2/5$, where the finite size effects are negligible. We introduce the site-diluted disorder similar to the previous experiments . We first select the diluted sites randomly with the probability $1-p$, then remove the nearest four bonds around the selected sites from the lattice. The distributions of the diluted sites are the same for all the samples considered. The percolative threshold concentration $p_{c}$ is about $0.592$ .

The resistivity-shunted-junction (RSJ) dynamics is incorporated in the simulations, which can be described as $${\frac{\sigma \hbar }{2e}}\sum_{j}(\dot{\phi _{i}}-\dot{\phi _{j}})=-{\frac{%
\partial H}{\partial \phi _{i}}}+J_{{\rm ex},i}-\sum_{j}\eta
_{ij},$$ where $\sigma$ is the normal conductivity, $J_{ex,i}$ refers to the external current, $\eta_{ij}$ denotes the thermal noise current with $<\eta_{ij}(t)>=0$ and $<\eta_{ij}(t)\eta_{ij}(t')>=2\sigma
k_{B} T\delta(t-t^{'})$.

The fluctuating twist boundary condition is applied in the $xy$ plane to maintain the current, thus the new phase angle $\theta_{i}=\phi_{i}+r_{i}\cdot\Delta$ $(
\Delta=(\triangle_{x},\triangle_{y})$ is the twist variable) is periodic in each direction. In this way, supercurrent between site i and site j is given by $J^{s}_{i\rightarrow
j}=J_{ij}sin(\theta_{i}-\theta_{j}-A_{ij}-r_{ij}\cdot \Delta),$ and the dynamics of $\Delta_{\alpha}$ can be written as

$$\dot{\Delta}_{\alpha}=\frac{1}{L^{2}}\sum\limits_{<i,j>\alpha}[J_{i\rightarrow
j}+\eta_{ij}]-I_{\alpha},$$ where $\alpha$ denotes the $x$ or $y$ direction, the voltage drop in $\alpha$ direction is $V=-L \dot{\Delta}_{\alpha}$. For convenience, units are taken as $2e=\hbar=J_{0}=\sigma =k_{B}=1$ in the following. Above equations can be solved efficiently by a pseudo-spectral algorithm due to the periodicity of phase in all directions. The time stepping is done using a second-order Runge-Kutta scheme with $\Delta t=0.05$. Our runs are typically $(4-8)\times 10^{7}$ time steps and the latter half time steps are for the measurements. The detailed procedure in the simulations was described in Ref. . In this work, a uniform external current $I$ along $x$ direction is fed into the system.

Since RSJ simulations with direct numerical integrations of stochastic equations of motion are very time-consuming, it is practically difficult to perform any serious disorder averaging in the present rather large systems. Our results are based on one realization of disorder. For these very large samples, it is expected to exist a good self-averaging effect, which is confirmed by two additional simulations with different realizations of disorder. This point is also supported by a recent study of JJA's by Um $et$ $al$ . In addition, simulations with different initial states are performed and the results are nearly the same. Actually, the hysteric phenomenon is usually negligible in previous RSJ dynamical simulations on JJA's . For these reasons, the results from simulations with a unique initial state (random phases in this work) are accurate and convincing.

# RESULTS AND DISCUSSIONS

## Finite temperature phase transition

The $I$-$V$ characteristics are measured at different disorder strengths and temperatures. At each temperature, we try to probe the system at a current as low as possible. To check the method used in this work, we investigate the $I$-$V$ characteristics for $f=0,p=1.0$. As shown in Fig. 1(a), the slope of the $I$-$V$ curve in log-log plot at the transition temperature $T_{c}$ ($\approx
0.894$) is equal to 3, demonstrating that the $I$-$V$ index jumps from 3 to 1, consistent with the well-known fact that the pure JJA's experiences a KT-type phase transition at $T_{c}\approx 0.894$. Figs. 1(b) and (c) show the $I$-$V$ traces at different disorder strengths in unfrustrated JJA's, while Fig. 1(d) is for $f=2/5,p=0.65$. It is clear that, at lower temperatures, $R=V/I$ tends to zero as the current decreases, which follows that there is a true superconducting phase with zero linear resistivity.

It is crucial to use a powerful scaling method to analyze the $I$-$V$ characteristics. In this paper, we adopt the Fisher-Fisher-Huse (FFH) dynamic scaling method, which provides an excellent approach to analyze the superconducting phase transition . If the properly scaled $I$-$V$ curves collapse onto two scaling curves above and below the transition temperature, a continuous superconducting phase transition is ensured. Such a method is widely used , the scaling form of which in 2D is $$\label{eq:ffh}
V=I\xi^{-z}\psi_{\pm}(I\xi),$$ where $\psi_{+(-)}(x)$ is the scaling function above (below) $T_{c}$, $z$ is the dynamic exponent, $\xi$ is the correlation length, and $V \sim I^{z+1}$ at $T=T_{c}$.

Assuming that the transition is continuous and characterized by the divergence of the characteristic length $\xi\sim|T-T_{c}|^{-\nu}$ and time scale $t\sim\xi^{z}$, FFH dynamic scaling takes the following form $$\label{eq:scv}
(V/I)|T-T_{c}|^{-z\nu}=\psi_{\pm}(I|T-T_{c}|^{-\nu}).$$

On the other hand, to certify a KT-type phase transition in JJA's, a new scaling form is proposed as follows $$\label{eq:sc2}
(I/T)(I/V)^{1/z}=P_{\pm}(I\xi/T),$$ which can be derived directly from Eq. () after some simple algebra. The correlation length of KT-type phase transition above $T_{c}$ is well defined as $\xi\sim e^{({c/|T-T_{c}|})^{1/2}}$ and Eq. () is rewritten as $$\label{eq:sckt}
(I/T)(I/V)^{1/z}=P_{+}(I e^{({c/|T-T_{c}|})^{1/2}}/T).$$

We perform the dynamic scaling analysis at a strong disorder ($p=0.65$) in unfrustrated system ($f=0$). Using $T_{c}=0.24\pm0.01$, $z=1.2\pm0.02$ and $\nu=1.0\pm0.02$, an excellent collapse is achieved according to Eq. (), which is shown in Fig. 2. In addition, all the low-temperature $I$-$V$ curves can be fitted to $V\sim I$exp$(-(\alpha/I)^\mu)$ with $\mu=0.9\sim1.1$. These results certify a continuous superconducting phase with long-rang phase coherence. The critical temperature for such a strongly disordered system is very close to that in 2D gauge glass model ($T_{c}=0.22$) .

For $f=0,p=0.86$, we first still adopt the scaling form in Eq. () to investigate the $I$-$V$ characteristics. As displayed in Fig. 3, we get a good collapse for $T<T_{c}$ with $T_{c}=0.58\pm0.01$, $z=2.0\pm0.01$ and $\nu=1.4\pm0.02$, demonstrating a superconducting phase with power-law divergent correlation for $T<T_{c}$. Note that the collapse is poor for $T>T_{c}$, implying that the phase transition is not a completely non-KT-type one. Next, we use the scaling form in Eq. () to analyze the $I$-$V$ data above $T_{c}$. Interestingly, using $T_{c}=0.58$ and $z=2.0$ determined above, a good collapse for $T>T_{c}$ is achieved, which is shown in Fig. 4. That is to say, the $I$-$V$ characteristics at $T<T_{c}$ are like those of a continuous phase transition with power-law divergent correlation length while at $T>T_{c}$ are like those of KT-type phase transition, which are well consistent with the recent experimental observations .

To make a comprehensive comparison with the experimental findings in Ref. , we also investigate the FTPT in frustrated JJA's ($f=2/5$) at a strong site-diluted disorder ($p=0.65$). As shown in Fig. 5, a superconducting phase transition with power-law divergent correlation is clearly observed.

As is well known, non-KT-type phase transition in frustrated systems is a natural result. However, it is intriguing to see that in unfrustrated systems, one may ask what our results really imply and what is the mechanism. It has been revealed that in the presence of a strong random pinning which is produced by random site dilutions, a breaking of ergodicity due to large energy barrier against vortex motion may allow enough vortices to experience a non-KT-type continuous transition .

Interestingly, we recover the phenomena in experiments by the present model and give some insights into the FTPT. More information on the low-temperature phase calls for further equilibrium Monte Carlo simulations as in Refs. . Table I summarizes the critical temperatures at different frustrations and disorder strengths. One can find that, the critical temperature in unfrustrated system decreases with increasing diluted sites.

|        | f=0     | f=2/5   |
|:-------|:--------|:--------|
| p=0.95 | 0.85(2) | 0.16(2) |
| p=0.86 | 0.58(1) | 0.13(1) |
| p=0.7  | 0.27(2) | 0.12(1) |
| p=0.65 | 0.24(1) | 0.14(1) |

Summary of $T_{c}$.

The systems considered in our work are site-diluted JJA's, which are not the same as bond-diluted JJA's in Ref. . In bond-diluted systems the diluted bonds are randomly removed, while in the site-diluted systems, the diluted sites are randomly selected, then the nearest four bonds around the selected sites are removed. Although the JJA's in Ref. and the present work are diluted in different ways, it is interesting to note that the obtained exponents in FTPT are very close, possibly due to the similar disorder effect produced.

## Depinning transition and creep motion

Next, we turn to the ZTDT and the LTCM for the typical site-diluted JJA's systems mentioned above. Depinning can be described as a critical phenomenon with scaling law $V \sim (I-I_{c})^{\beta}$, demonstrating a transition from a pinned state below critical driving force $I_{c}$ to a sliding state above $I_{c}$. The $(I-I_{c}).vs.V$ traces at $T=0$ for $f=0,p=0.86$; $f=0,p=0.65$ and $f=2/5,p=0.65$ are displayed in Fig. 6, linear-fittings of Log$(I-I_{c}).vs.$Log$V$ curves are also shown as solid lines. As for $f=0,p=0.86$, the depinning exponent $\beta$ is determined to be $2.62\pm0.1$ and the critical current $I_{c}$ is $0.302\pm0.005$, while for the cases $f=0,p=0.65$ and $f=2/5,p=0.65$, the depinning exponents are evaluated to be $2.37\pm0.1$ and $2.27\pm0.05$ with the critical currents $I_{c}=0.039\pm0.001$ and $I_{c}=0.035\pm0.002$, respectively.

When the temperature increases slightly, creep motions can be observed. In the low-temperature regime, the $I$-$V$ traces are rounded near the zero-temperature critical current due to thermal fluctuations. Fisher first suggested to map such a phenomenon for the ferromagnet in magnetic field where the second-order phase transition occurs  . This mapping was then extended to the random-field Ising model   and the flux lines in type-II superconductors  . For the flux lines in type-II superconductors, if the voltage is identified as the order parameter, the current and the temperature are taken as the inverse temperature and the field respectively, analogous to the second-order phase transition in the ferromagnet, the voltage, current and the temperature will satisfy the following scaling ansatz   $$\label{eq:sca}
V(T,I)=T^{1/\delta}S[(1-I_{c}/I)T^{-1/\beta\delta}].$$ Where $S(x)$ is a scaling function. The relation $V(T,I=I_{c})=S(0)T^{1/\delta}$ can be easily derived at $I=I_c$, by which the critical current $I_{c}$ and the critical exponent $\delta$ can be determined through the linear fitting of the Log$T$-Log$V$ curve at $I_{c}$.

The Log$T$-Log$V$ curves are plotted in Fig. 7(a) for $f=0,p=0.86$. We can observe that the critical current is between 0.3 and 0.32. In order to locate the critical current precisely, we calculate other values of voltage at current within (0.3,0.32) with a current step 0.01 by quadratic interpolation . Deviation of the $T$-$V$ curves from the power law is calculated as the square deviations $SD=\sum[V(T)-y(T)]^{2}$ between the temperature range we calculated, here the functions $y(T)=C1T^{-C2}$ are obtained by linear fitting of the Log$T$-Log$V$ curves. The current at which the $SD$ is minimum is defined as the critical current. The critical current is then determined to be $0.302\pm0.001$. Simultaneously, we obtain the exponent $1/\delta =1.688\pm0.001$ from the slope of Log$T$-Log$V$ curve at $I_{c}=0.302$. The similar method is applied to investigate the cases $f=0,p=0.65$ and $f=2/5,p=0.65$. As shown in Figs. 7 (b) and (c), the critical current $I_{c}$ and critical exponent $1/\delta$ for $f=0,p=0.65$ are determined to be $0.03875\pm0.0005$, $2.24\pm0.02$ respectively, for $f=2/5,p=0.65$, the result is $I_{c}=0.034\pm0.001$, $1/\delta=2.29\pm0.01$.

We then draw the scaling plots according to Eq. (). By one-parameter tuning (only $\beta$), we get the best collapses of data in the regime $I \leq I_{c}$ with $\beta=2.61\pm0.02$ and $2.28\pm0.02$ for $f=0,p=0.86$ and $f=0,p=0.65$ respectively, which are shown in Figs. 8(a) and (b). For $f=0,p=0.86$, this curve can be fitted by $S(x)=0.0994$exp$(1.9x)$, combined with the relation $\beta\delta=1.55$, suggesting a non-Arrhenius creep motion. However, for the strongly site-diluted system with $f=0,p=0.65$, the scaling curve can be fitted by $S(x)=0.037$exp$(0.5x)$, combined with the relation $\beta\delta \approx 1.0$, indicative of an Arrhenius creep motion. Interestingly, as displayed in Fig. 8(c) for $f=2/5,p=0.65$, the exponent $\beta$ is found to be $2.30\pm0.02$, which yields $\beta \delta \approx 1.0$. The scaling curve in the regime $I \leq I_{c}$ can be fitted by $S(x)=0.105$exp$(0.25x)$. These two combined facts suggest an Arrhenius creep motion in this case.

It is worthwhile to note that both the FTPT and the LTCM for strongly disordered JJA's ($p=0.65$) with and without frustration are very similar. The $I$-$V$ curves in low temperature for all three cases can be described by $V\propto T^{1/\delta}$ exp$[
A(1-I_{c}/I)/T^{\beta\delta}]$, this is one of the main characteristics of glass phases , while the $I$-$V$ traces for KT-type phases can be fitted to $V \propto I^{a}$. Hence, by the scaling ansatze in Eq. (), we have provided another evidence for the existence of non-KT-type phases in the low-temperature regime for these three cases ($f=0,p=0.86;
f=0,p=0.65; f=2/5,p=0.65$).

# SUMMARY

To explore the properties of various critical phenomena in site-diluted JJA's, we performed large-scale simulations for two typical percolative strengths $p=0.86$ and $p=0.65$ as in a recent experimental work . We investigated the FTPT, the ZTDT and the LTCM in these systems. The RSJ dynamics was applied in our work, from which we measured the $I$-$V$ characteristics at different temperatures.

The results obtained in this work about FTPT are well consistent with the recent experimental findings in Ref. and are inconsistent with the earlier experimental study in Ref. , possibly due to the large noise in the measurement of voltage in Ref. (larger than 0.2nv), which was considerably reduced in the experiments by Yun et al . The evidence for non-KT-type phase transition was revealed by two different scaling ansatzes (Eq. () and Eq. ()). Our results also shed some light on the various phases and the phase transitions where the different divergent correlations at various disorder strengths were suggested, and the critical exponents were evaluated. Furthermore, the results in this paper are useful for understanding not only the site-diluted systems, but also the whole class of disordered JJA's, for instance, the combination of two different phase transitions may exist in other disordered JJA's systems.

In addition, the ZTDT and the LTCM were also touched. It was demonstrated by the scaling analysis that the creep law for $f=0,p=0.86$ is non-Arrhenius type while those for $f=0,p=0.65$ and $f=2/5,p=0.65$ belong to the Arrhenius type. It is interesting to note that the non-Arrhenius type creep law for $f=0,p=0.86$ is similar to that in three-dimensional flux lines with a weak collective pinning . The product of the two exponents $1.55$ is also very close to $3/2$ determined in Ref. . For $f=0,p=0.65$ and $f=2/5,p=0.65$, the observed Arrhenius type creep law is also similar to that in the glass states of flux lines with a strong collective pinning as in Ref. . Future experimental work is needed to clarify this observation.

# ACKNOWLEDGEMENTS

This work was supported by National Natural Science Foundation of China under Grant Nos. 10774128, PCSIRT (Grant No. IRT0754) in University in China, National Basic Research Program of China (Grant Nos. 2006CB601003 and 2009CB929104),Zhejiang Provincial Natural Science Foundation under Grant No. Z7080203, and Program for Innovative Research Team in Zhejiang Normal University.

$^{\dag}$ Corresponding author. Email:email@example.com