abstract: The feasibility of a next generation underground water Cherenkov detector is examined and a conceptual design (UNO) is presented. The design has a linear detector configuration with a total volume of 650 kton which is 13 times the total volume of the Super-Kamiokande detector. It corresponds to a 20 times increase in fiducial volume for physics analyses. The physics goals of UNO are to increase the sensitivity of the searches for nucleon decays about a factor of ten and to make precision measurements of the solar and atmospheric neutrino properties. In addition, the detection sensitivity for Supernova neutrinos will reach as far as the Andromeda galaxy.
author: Chang Kee Jung
title: Feasibility of a Next Generation Underground Water Cherenkov Detector: UNO [^1]

# Introduction

Large scale underground water Cherenkov detector experiments, Super-Kamiokande and its predecessors (IMB and Kamiokande), have been extremely successful in producing crucial physics results during the last two decades. Their accomplishments includes: the first real time measurement of the solar neutrinos, confirmation of the solar neutrino flux deficit, observation of the neutrino oscillations in the atmospheric neutrinos, observation of neutrinos from the Supernova 1987A, and setting the world best limits on the nucleon decays.

These detectors were originally conceived for searches for nucleon decays predicted by various Grand Unification Theories (GUTs). While no positive observation of nucleon decays have been made to date with these detectors, the evidence for neutrino oscillations now firmly established by the Super-Kamiokande atmospheric neutrino analysis provides us with a breakthrough in particle physics beyond the Standard Model. This finding indicates that the neutrino masses are indeed very small if we assume no degeneracy in mass eigenstates, which in turn indicates that there may be a new very high energy physics scale that facilitates small neutrino masses via "See-saw" mechanism and allows protons to decay.

There are many theoretical models that predict proton decays and some of them were presented in this workshop (NNN99). A specific example of such models can be found in Ref.  which presents a complete and detailed description of interplay between the neutrino masses, proton decays and other Standard Model observables in G(2,2,4) and SO(10) frameworks. The model predicts proton decay lifetimes within the reach of the Super-Kamiokande, especially in the SUSY favored decay modes. These predictions along with many other predictions from other models encourage us to turn our attention back to the proton decay searches with higher sensitivity.

If discovered, proton decay provides only and unambiguous direct evidence of existence of GUT scale physics at low energies. It certainly will be remembered as one of the most revolutionary discoveries in particle physics history. Some say the discovery is around the corner.

In order to further our effort to search for nucleon decays which I believe is in the category of "must-do" physics and to do high statistics studies of neutrino physics including solar neutrinos, atmospheric neutrinos and supernova neutrinos, I propose a next generation large underground water Cherenkov detector which is named UNO (Ultra underground Nucleon decay and neutrino Observatory).

In the following sections, the design considerations, baseline configuration, physics capabilities and cost estimation of the UNO detector is described in some detail.

# Ultra underground Nucleon decay and neutrino Observatory Detector

The design philosophy of UNO is to make a relatively simple extension of the well established water Cherenkov detector technology beyond Super-Kamiokande to achieve an order of magnitude better sensitivity in nucleon decay searches and to study various neutrino physics with higher precision. With water Cherenkov detector technology, we can utilize the tremendous amount of experience and expertise gained from the IMB, Kamioka and Super-Kamiokande detectors.

In order to establish reasonable detector parameters for physics capability studies and cost estimates, I set the benchmark fiducial volume of the UNO detector to be 20 times that of Super-Kamiokande. The design of the detector is kept to be simple and robust, and broad physics capabilities are required.

Several design options are considered keeping in mind the practical limitations in conventional water Cherenkov detector technique: namely, (1) the largest depth of water of a detector is limited by the current PMT pressure stress limit ( 8 atm for current 20$^{\prime\prime}$ Hamamatsu PMTs) which can be overcome if one is willing use an optical high pressure water tight container for each PMT and compromising the PMT efficiency; 2) The maximum dimension of the detector without active detection element is limited by the finite light attenuation length in pure water ($\sim$`<!-- -->`{=html}80 m @400 nm in Super-Kamiokande ).

Three different detector shapes are considered: Big, Torus and Linear. According to D. L. Petersen who is a rock engineer at the CNA Consulting Engineering at Chicago, Illinois, the excavation costs are more or less same independent of the shapes of the detector. Then, if we want to keep the detector cost at minimum for a fixed fiducial volume (445 kton; $\sim$`<!-- -->`{=html}20 times the Super-Kamiokande), we need to keep the ratio $rv$=(fiducial volume/total volume) large and keep $rs$=(PMT surface area/total volume) small. These requirements obviously prefer fat detectors.

The Big detector design option considered has a cube shape with 86x86x86m outer dimensions. It has a potential problem with the PMT pressure limit for the PMTs at the bottom of the detector where the water pressure will be about 9 atmospheric pressure. It also has a diagonal length of 150 m which is much longer than the light attenuation length currently measured in the pure water.

The Torus detector design option considered turns out to be very inefficient in term of the $rv$ value and is physically not possible if the cross-section is 60x60 m for a benchmark detector size. Even with a 50x50 m cross-section a torus would be too tight, *i.e.* the diameter of the central rock column will be too small to support the structure.

Thus, a Torus design will have to have a small cross-section which results in a small $rv$ and large $rs$ values. For example, a 40x40m cross-section torus has $rv$=0.6 compared to 0.7 for a Linear option considered below.

The Linear detector design option considered has a 60mx60mx180m outer dimensions and it appears to be the most optimal one. When compartmentalized as three 60mx60mx60m cubes in terms of active detection elements, it satisfies all of the requirements mentioned above and results in a reasonable $rv$ and $rs$ values. The compartmentalization option naturally makes the detector cost more expensive but provides several significant benefits compared to an open geometry option. First, it minimizes so-called flasher background events, which commonly occur in all water Cherenkov detectors that use PMTs, by confining them in an isolated compartment. Second, it minimizes inefficiency and variation in efficiency due to finite light attenuation length in water. Third, it keeps the detector operational live time to close to 100%. The current detector operational live time for Super-Kamiokande is about 90%. The 10% inefficiency is mainly due to detector calibrations. With compartmentalization we can ensure that at least one compartment will be alive during the detector calibration times. Keeping the detector live time very high (virtually at 100%) is of course very important for physics programs such as a supernova watch.

Another issue that needs to be dealt with concerning the compartmentalization is whether we should make the divisions among the compartments rigid so that the water can be filled and drained independently. Certainly making the divisions rigid walls will make any future repairing of the detector easier. (It should be noted that it takes two months to drain the current Super-Kamiokande tank.) However, it is not clear whether this is an important enough reason for us to consider such an option which will incur substantial additional cost.

One could also argue about whether the next generation detector should be built underground or underwater. Indeed there are a few ideas presented in recent meetings that consider very large underwater detectors. In my opinion, however, one of the most serious disadvantage of a underwater detector will be its inaccessibility for calibration and repair. The experience gained from the Super-Kamiokande tells us that a well selected and maintained underground detector environment is a wonderful environment for us to work and to operate the detector. It provides us an easy access to the detector for various calibrations and for trying out new ideas. On the other hand, I would assume that precision calibrations of a deep underwater detector will be extremely difficult if not impossible, especially if one is aiming to do a calibration for solar neutrino measurements using electron linac or DT generator etc. Thus, for a underwater detector design to become a viable option, a new sets of elaborate remote calibration systems has to be developed and tested. In addition, location and accessibility of other service facilities such as water purification system needs to be considered, and a practical and realistic solution must be found. Overall, I expect many more technical challenges for a underwater detector than for a underground detector.

Considering all of the above issues, I believe that, a large underground water Cherenkov detector with a linear configuration is the best option for a next generation nucleon decay and neutrino detector which can be built within next 10 years without requiring major R&D. A conceptual design of such a detector (UNO detector) is shown in Figure .

The detector has three compartmentalized sections with 60mx60mx60m dimensions resulting in a total length of 180 m and a total volume of 648 kton. The outer detector region of the detector has 2.5 m depth of veto shield and is instrumented with 14,901 8$^{\prime\prime}$ PMTs with a PMT density of 0.33 PMTs/m$^2$. The inner detector region has a total fiducial volume of 445 kton which is defined as the water volume 2 m inside of the inner PMT planes. The inner detector region is instrumented with 20$^{\prime\prime}$ PMTs with a PMT density of 1.96 PMTs/m$^2$ (40% photocathode coverage a la Super-Kamiokande) for the central section and 0.49 PMTs/m$^2$ (10% photocathode coverage) for the two sections at the wings. The total number of 20$^{\prime\prime}$ PMTs is 56,650 for this configuration.

## UNO Physics Goals and Capabilities

The detector configuration proposed above presents excellent physics capabilities of UNO for broad ranges of nucleon decay and neutrino physics. The two wings of the detector with 10% photocathode coverage will have an energy threshold below 10 MeV which should be sufficient for nucleon decay searches, supernova neutrino detections, high statistics atmospheric neutrino studies and high statistics solar neutrino studies at the high energy end of the energy spectrum. The central compartment with a 40% photocathode coverage will serve as a dedicated low energy solar neutrino detector with a 5 MeV analysis threshold which will also provide UNO with a detection capability for 6 MeV $\gamma$'s from $p\rightarrow\nu K^+$ decays in the oxygen nucleus. It also provides the capability to observe low energy supernova neutrino events ($\sim 5$ MeV) which are important for core collapse modeling.

The primary physics goal of UNO is to obtain an order of magnitude better sensitivity in nucleon decay searches (especially in $p \rightarrow e^+\pi^0$ mode) than Super-Kamiokande. The exact amount of improvements in the sensitivity vary depending on the decay modes, particularly on their expected backgrounds. Obviously for the background free decay modes we expect more than a order of magnitude improvements but for the background limited decay modes we would expect a factor of 4 or 5 improvements.

The amount of improvement also strongly depends on the optimization of the analyses for high statistics sample. For example, if we continue current standard Super-Kamiokande analysis for $p \rightarrow e^+\pi^0$, one will encounter a background limited analysis situation in ten years (0.5 years for UNO). However, if we make a tighter cut on total momentum in an event, say 150 MeV (rather than 250 MeV for current Super-Kamiokande standard analysis), the signal efficiency will go down to about 25-30% from current 44%, but we will be virtually background free for many UNO years. Since we gain a factor of 20 in total fiducial volume, by keeping the efficiency better than 22% we can insure a factor of 10 or more gain in the search sensitivity as long as we can keep the background to be zero.

In order to estimate the potential improvements reliably, we need to do much more work. First of all, we need to have better understanding of the characteristics of the atmospheric neutrino background events. This can be accomplished by studying the K2K 1 kton neutrino events which are similar in many aspects to the atmospheric neutrino background events relevant to nucleon decay searches. Second, we need to reduce the systematic uncertainties associated with the atmospheric neutrino flux and the neutrino interaction cross-section for water target. Some of this information can come from the K2K fine-grained detector data. Finally, we need to improve the search analyses by developing more sophisticated algorithms and optimize the analyses for maximum sensitivity. I personally believe that there is much room for improvement in these areas, especially for the multi-ring event analysis. We also need to generate much larger MC events samples with detailed simulations for these purposes. While it is difficult to predict the exact sensitivity of the UNO detector for nucleon decay searches at this time, I am convinced that an order of magnitude improvements in the search sensitivity for various decay modes are attainable.

Another equally important physics goal of UNO is supernova neutrino detection. With the proposed configuration UNO will be able to record 100k neutrino interaction events from a supernova explosion at 10 kpc away. This will allow us to make a detailed mapping of the time structure of the neutrino flux providing valuable information for the theoretical modeling of the supernovae. The detection reach for neutrinos from supernova explosions will be about 1 Mpc from the earth which includes most of the local group of galaxies including Andromeda. We expect (optimistically) about one supernova explosion in every ten years within this range.

Although current Super-Kamiokande measurements of the solar neutrino properties provide crucial information on solar neutrino problems with unprecedented statistics, the measurements on energy spectrum, day-night flux ratio, and seasonal variations are limited by the lack of statistics in providing an unambiguous resolution to the solar neutrino problems. Thus, it is strongly desirable to build a large scale detector that has a capability of precision studies of the solar neutrinos. The proposed configuration of UNO will provide about 20 times more statistics for high energy end of the solar neutrino events and at least 7 times more statistics for the low energy end events. Even if we are able to establish a solution to the solar neutrino problem with the results from Super-Kamiokande, SNO, KAMLAND and Borexino within this decade, we will need more precision measurements to firmly establish the SSM (Standard Solar Model) and neutrino oscillation parameters. I believe that we will be entering a precision measurement era for SSM in this decade as the two last decades have been for the Standard Model.

If we assume the neutrino oscillations observed in the Super-Kamiokande atmospheric neutrino events are due to $\nu_{\mu}$ to $\nu_{\tau}$ oscillations, there should be about 20 tau appearance events produced per year from the $\nu_{\tau}$ interactions in the Super-Kamiokande detector. For UNO the rate will be about 400 tau appearance events per year. While extracting a tau signal from background in water Cherenkov detector is not expected to be easy, it will be a worth while project to pursue, especially with a reasonably high statistics event sample. (Currently, we are pursuing such an analysis with the Super-Kamiokande data and the initial results appear to be promising. We are hoping to get about three sigma effect.) If the $\Delta m^2$ turns out to be relative low, say $\sim 3\times 10^{-3}$ eV$^2$ or lower, the long baseline neutrino experiments will have difficult time in observing unambiguous signal of tau appearance. And thus, establishing a direct observation of tau appearance coming from neutrino oscillations may take much longer time than we currently hope for.

Finally, a possibility of using UNO as a far detector of a Muon Storage Ring produced neutrino beam experiment should be considered. A water Cherenkov detector is a good candidate for a far detector. It provides an relatively inexpensive large target volume, particle ID for electrons and muons, and reasonable calorimetry. It could also provide good charge separation if it is aided by an external spectrometer as shown in Figure .

Thus, if UNO is built, it will serve as a truly multi-purpose general underground detector for variety of great physics. Especially for nucleon decay searches, it will serve as a robust multi-decay mode search detector.

## UNO Cost Estimates

In order to obtain a ball park figure cost estimate for UNO detector, I use the construction cost data for the Super-Kamiokande detector provided by K. Nakamura as a reference. I divide the major expense items into two categories according to their nature of scaling, i.e. volume-like scaling or surface-like scaling. And I apply my reasonable guesses to ultimately determine the scaling factor from Super-Kamiokande to UNO. Thus, these figures should not be taken too seriously. Table  shows my initial estimate of itemize costs for UNO along with the actual costs for Super-Kamiokande. The total cost for UNO with the baseline configuration proposed in this paper is then about \$520M and the cost for a detector with full 40% photocathode coverage for all inner surfaces (a la Super-Kamiokande) is about \$680M. And all of these figures do not include contingencies.

| Item                                       |  SuperK |     |     UNO |
|:-------------------------------------------|--------:|:---:|--------:|
| Cavity Excavation                          |  27,640 |  v  | 168,000 |
| Water piping and pumps                     |     630 |  v  |   4,082 |
| Water Purification Sys                     |   1,850 |  v  |  11,988 |
| Power Station                              |     720 |  v  |   2,160 |
| Crane                                      |     760 |  v  |   2,280 |
| Water Tank                                 |  18,400 |  s  |  92,480 |
| PMT support structure                      |   4,580 |  s  |  23,019 |
| Counting Room                              |     330 |  s  |     990 |
| Computer Building                          |   1,860 |  s  |   2,232 |
| Main Building                              |   3,000 |  s  |   3,600 |
| 20$^{\prime\prime}$ PMT (including cables) |  34,670 |  s  | 173,664 |
| Electronics                                |   6,330 |  s  |   9,495 |
| DAQ                                        |   1,090 |  s  |   1,635 |
| Air Conditioning                           |     210 |  s  |     315 |
| Veto instrumentation                       |   3,000 |  s  |   9,000 |
| 8$^{\prime\prime}$ PMT (including cables)  |   2,262 |  s  |  17,881 |
| Total                                      | 102,070 |     | 522,822 |

UNO Cost Estimates (in thousands of US dollars): This is a quite complete list of expense items for all aspect of detector construction. The "v" and "s" symbols note the volume-like and surface-like nature of cost scalings used when the costs are extrapolated from the Super-Kamiokande (SuperK) to UNO. A conversion rate of 1\$ = 100 yen is used.

For the major cost items of the detector, I use more reliable quotes and estimates for unit costs from vendors and experts. As shown in Table , the unit cost for excavation is taken from the estimates by D. L. Petersen, and the unit costs for PMTs are taken from Hamamatsu corporation. I take liberty of reducing the per channel electronic cost by a factor of 5 and the water tank construction cost per unit surface by a factor of 2 from the Super-Kamiokande cost. The cost of these five items comprise the majority of the UNO detector cost.

| Item                     | Unit Cost     | Source      |
|:-------------------------|:--------------|:------------|
| Excavation               | \$260/m$^3$   | L. Petersen |
| 20$^{\prime\prime}$ PMTs | \$3,100       | Hamamatsu   |
| 8$^{\prime\prime}$ PMTs  | \$1,200       | Hamamatsu   |
| Electronics              | \$170/channel |             |
| Water Tank               | \$2,076/m$^2$ |             |

Estimated Unit Costs: The excavation cost is assuming a horizontal access tunnel and rock quality (Q value) of 100. The PMT unit cost including cable cost is based on a 50k PMT order. It is \$2,850 if 100k PMTs are ordered. A conversion rate of 1\$ = 100 yen is used.

## Cost Reduction Possibilities

I believe that there is quite a bit of room to reduce the total cost for UNO. First of all we can reduce the excavation cost by finding an existing cavity or underground facility. Or we could simply get a better quote. It is said that the actual excavation cost charged to the mining company in Kamioka, Japan is only \$50/m$^3$. Second, we could use Mineguard type of liner rather than a stainless-steel container for the container tank structure. Third, we should optimize the PMT size and granuality for our physics goals and look into a way of developing new cheaper photo-detectors. For example, we could use two 8$^{\prime\prime}$ PMTs instead of a 20$^{\prime\prime}$ PMT, which may result in overall better performance and cheaper cost. We could also look for other vendors (from Russia?) than Hamamatsu.

## Site Selection

Obviously, the issue of site selection is an important issue for the success of a project like UNO with this magnitude. In principle, a site near the Equator is desired if we are aiming for maximum information from UNO for the solar neutrino physics. However, the consensus among solar neutrino theorists indicate that it is not crucial unless the MSW small mixing angle solution looks very promising. Thus, I concentrate my attention to possible sites in the United States.

An extensive site search has been conducted by R.R. Sharp, Jr. and A. Mann during the early 1980s for a site for potential National Underground Science Facility (NUSF). They surveyed mostly in the western states for existing tunnels, inactive mines, active mines and new locations. They found a couple of sites in Nevada particularly attractive. Although the size and depth requirements of the potential sites they searched for are quite different from the UNO requirements, the information gathered is found to be quite useful.

Another site search was conducted by W. Kropp *et. al.* concentrating on the site at San Jocinto, California. The results from this search also provides useful information for our current search.

The most attractive site so far presented to us is the WIPP (Waste Isolation Pilot Project) site at Carlsbad, New Mexico. The facility which is managed by DOE provides a hard salt rock geology, an already existing laboratory facility and tunneling machinery. With strong government interest to utilize this facility for scientific research, the site is indeed ideal for UNO or for any other large scale next generation underground detectors.

# Conclusions

Feasibility of a next generation very large underground water Cherenkov detector, UNO, is considered. The detector configuration of UNO is a linear compartmentalized detector (60mx60mx60mx3) with a middle section with a 40% photocathode coverage and two wing sections with 10% photocathode coverage. The middle section is dedicated for low energy solar neutrino and supernova neutrino studies and for detecting 6 MeV prompt $\gamma$'s from $p \rightarrow \nu K^+$ decays.

The detector which has a total volume of 650 kton and a corresponding fiducial volume of 445 kton (20 times larger than the fiducial volume of the Super-Kamiokande detector) can be built today at about \$500M of construction cost. There are no known serious technical challenges in building such a detector. Although rigorous R&D is desired for cost reduction and improvement in detector performance, there are no critical path R&D items. If funding is available, we should be able to build this detector within next ten years.

Such a detector, if built, will provide us with a bonanza of exciting physics programs: Nucleon decay searches, precision solar and atmospheric neutrino studies, and supernova neutrino observations. The detector can also serve as a far detector for a possible future neutrino factory. It will also compliment any other major accelerator based initiatives such as NLC or neutrino factories in terms of providing diverse but crucial physics programs to the US HEP community.

In my opinion, physics program of searches for nucleon decays undoubtedly belongs to "MUST-DO" physics category and it should be vigorously pursued. I also believe if such a large scale detector is built it should be a multi-purpose general detector with a robust multi-decay mode search capability for nucleon decay searches. It should not be specialized single purpose detector with specific decay mode search capability.

# Acknowledgements

This work is supported by the funding from the U.S. Department of Energy under the contract No. DEFG0292ER40697. I wish to express my gratitude to Kenzo Nakamura who provided vital information on the Super-Kamiokande construction costs, D. Lee Petersen who provided information on the technical details concerning underground excavation, Al Mann who provided information on the underground site searches. I also like to thank John Bahcall, Bill Marciano, Jogesh Pati and Frank Wilczek for their theoretical inputs and encouragements, and Maurice Goldhaber and other workshop organizing committee members for their tireless work in making the workshop successful and support for this work.

[^1]: To be published in the Proceedings of the Next generation Nucleon decay and Neutrino detector (NNN99) Workshop, Sep. 23-25, 1999, Stony Brook, New York


abstract: The two charge-density wave (CDW) transitions in NbSe$_3$ were investigated by scanning tunneling microscopy (STM) on *in situ* cleaved $(\bm{b},\bm{c})$ plane. The temperature dependence of first-order CDW satellite spots, obtained from the Fourier transform of the STM images, was measured between 5-140 K to extract the surface critical temperatures (T$_s$). The low T CDW transition occurs at T$_{2s}$=70-75 K, more than 15 K above the bulk T$_{2b}=59$K while at exactly the same wave number. Plausible mechanism for such an unusually high surface enhancement is a softening of transverse phonon modes involved in the CDW formation. The regime of 2D fluctuations is analyzed according to a Berezinskii-Kosterlitz-Thouless type of surface transition, expected for this incommensurate 2D CDW, by extracting the temperature dependence of the order parameter correlation functions.
author: Christophe Brun$^{1,*}$, Zhao-Zhong Wang$^1$, Pierre Monceau$^2$, and Serguei Brazovskii$^3$
bibliography: biblioNbSe3.bib
date: 13-01-2010
title: Surface Charge Density Wave Transition in NbSe$_3$

Powerful techniques have been recently developed, such as grazing incidence inelastic x-ray scattering, angle resolved photoemission, or scanning tunneling microscopy (STM), to probe the surface layer(s) of correlated electron systems. An important question is currently raised regarding the surface electronic states of these systems: are they identical to the bulk ones? This question becomes even more crucial when the material undergoes a phase transition to a broken symmetry state, for instance a charge-density wave (CDW) transition; the free surface must reflect this broken symmetry. The "extraordinary transition" corresponds to the very interesting case where the surface orders at higher temperature than the bulk . It opens the opportunity to study a truly two-dimensional (2D) system which is more perfect than usual ones prepared on a substrate.

Up to now experimental systems known to reveal such surface transitions are extremely scarce . The antiferromagnetic system NiO(100) is one of these . Recently, at the surface of the quasi-2D compound NbSe$_2$, the CDW was shown to order at about $1.5$ K above the bulk transition temperature . This was subsequently corroborated by emphasizing a modified behavior of the surface phonon modes . On the other hand "the commensurate checkerboard order" discovered by STM in the cuprate system Ca$_{2-x}$Na$_x$CuO$_2$Cl$_{2}$ (NaCCOC) , furnishes a example of charge ordering stabilized at the surface whereas the bulk doesn't show such counterpart at lower temperature. Hereafter, we report high-resolution *in situ* STM measurements on the quasi-one dimensional (1D) compound NbSe$_3$ which undergoes in the bulk two independent CDW transitions at $T_{1b}=144$ K and at $T_{2b}=59$ K, involving two CDW vectors, $\bm{q_1} = 0.24$ $\bm{b^*}$ for the high-temperature (HT) CDW and $\bm{q_2}= 0.5$ $\bm{a^*} + 0.26$ $\bm{b^*} + 0.5$ $\bm{c^*}$ for the low-temperature (LT) CDW. By measuring the temperature dependence of the amplitude of the $\bm{q_1}$ and $\bm{q_2}$ first-order satellites we found that at the $(\bm{b},\bm{c})$ surface of NbSe$_3$, the $\bm{q_2}$ transition occurs at 15 K higher than in the bulk, demonstrating the huge effect of the surface on the LT CDW phase transition. Owing to the phase degeneracy of this incommensurate CDW and to the anisotropic electronic properties of NbSe$_3$, the Berezinskii-Kosterlitz-Thouless (BKT) transition is investigated for the first time in real space.

NbSe$_3$, one of the most widely studied CDW system, was the first inorganic low-dimensional conductor to exhibit a sliding CDW state when an electric field above a threshold value is applied . It has a linear structure consisting of three pairs of metallic chains per unit cell running along the $\bm{b}$ axis, which are denoted as type I, II and III according to the strength of the chalcogen-chalcogen bond in the triangular basis of the chain. By STM we could unambiguously identify for the first time the three types of chains lying in the $(\bm{b},\bm{c})$ surface plane in the temperature range 5-140 K . The surface CDW wave-vectors are $\bm{q_1} = 0.24$ $\bm{b^*}$ and $\bm{{q_{2p}}}= 0.26$ $\bm{b^*}$ $+$ $0.5$ $\bm{c^*}$, in excellent agreement with the bulk reported values projected on the $(\bm{b},\bm{c})$ plane. A detailed analysis of the spatial distribution of the $\bm{q_1}$ and $\bm{q_2}$ CDWs on the various types of chain was performed and compared to bulk x-ray diffraction results . A long range new modulation defined by the wave-vector $\bm{u} \simeq 2 \times (0.26 - 0.24)$ $\bm{b^*}$ was observed, resulting from the interaction of both CDWs being present on chains III .

The present work was performed with an Omicron LT ultrahigh vacuum (UHV) STM system equipped with two UHV separated chambers. Well characterized NbSe$_3$ single crystals with typical dimensions of $0.01 \times 10 \times 0.05$ mm$^3$ were selected, cleaved *in situ* at room temperature along the $(\bm{b},\bm{c})$ planes and cooled down to 5, 63 or 78 K. Samples were further thermalized at temperatures between 5 and 140 K. Our thermometer is a silicon diode in good thermal contact with the sample. Intermediate temperatures were stabilized by a feedback system controlling a heater element while the STM head is continuously cooled by the cryostat through copper braids. Experiments were conducted when the thermal drift became acceptable. Both mechanically sharpened Pt/Ir and electrochemically etched W tips were used leading to similar results. All the STM images shown in the following are measured with constant current.

Fig.  shows three STM images measured at 78, 63 and 5 K, on large atomically flat terraces, with their 2D Fourier transform (2D FT). Strikingly at 78 K, almost 20 K above $T_{2b}$, one can easily detect the presence of diffuse $\bm{{q_{2p}}}$ superlattice spots, having a much lower amplitude than the $\bm{q_1}$ ones. At 63 K, i.e. 4 K above $T_{2b}$, the $\bm{{q_{2p}}}$ superlattice is well developed leading to sharp spots in the FT of the STM image, their amplitude being larger than the one of the $\bm{q_1}$ satellites. This indicates that the $\bm{q_2}$ CDW ordering occurs at higher temperature at the surface than in the bulk. It contrasts with the conventional behavior of the $\bm{q_1}$ CDW for which we do not find a noticeable increment of the surface transition temperature over the bulk one.

The $\bm{q_2}$ CDW surface transition temperature surface ($T_{2s}$) was determined by analyzing the temperature dependence of the amplitude of the first-order $\bm{{q_{2p}}}$ satellites, extracted from 2D FT of STM images (areas: $20\times20$ to $40\times40$ nm$^2$) . A selection of the measured STM images is presented in Fig. a for $5<T<78$ K. As we have shown in , only chains I and III are clearly visible in these tunneling conditions, while the location of chains II is shown for clarity, at 5 and 78 K. The corrugation of the $\bm{{q_{2p}}}$ CDW is about 15 pm at 5 K on chains I, has very weak variations up to 62 K, and decreases strongly between 62 K and 78 K. This is accompanied by a simultaneous diminution of the beating effect between $\bm{q_1}$ and $\bm{{q_{2p}}}$ on chains III (see ). At 78 K, the $\bm{{q_{2p}}}$ modulation nearly vanishes and the $\bm{q_1}$ wavefronts on chains III becomes quasi-parallel and regular.

Quantitatively, the extracted $\bm{{q_{2p}}}$ Fourier amplitude as a function of temperature is shown in Fig. b, with x-ray diffraction results and CDW-CDW tunneling data . A sharp transition is observed by STM in the range 70-75 K, much above $T_{2b}=59$ K . From the results of the $T_{2s}$ of an unfree surface can be extrapolated to $\approx60$ K. Hence there is evidence for an "extraordinary transition" occurring at the $(\bm{b},\bm{c})$ surface of NbSe$_3$ for the $\bm{q_2}$ CDW, leading to CDW ordering at a temperature exceeding the bulk one by about 15-20%. The temperature dependence is more abrupt in the vicinity of $T_{2s}$ at the surface than in the bulk. The dispersion of the experimental data reflects that at different positions along the surface the measured CDW corrugation is not identical, particularly in the vicinity of the transition around 70 K.

To explain the increase of $T_{2b}$ at the surface plane, it is appealing to consider the differences existing between phonons propagating at the surface and in the bulk. At a crystal surface, atom displacements normal to the surface are larger than those encountered in the bulk . Energies of the transverse phonon modes propagating in the surface plane are then smaller than those of the bulk phonons. The charge ordering occurring at the surface of NaCCOC has been recently explained in this way . Also the inelastic x-ray diffraction measurements on NbSe$_2$ have shown that the Kohn anomaly at $2\bm{k_F}$ is more pronounced at the surface and occurs at a lower energy than in the bulk underlying the changes in the phonon spectrum to explain the increase of $T_c$ at the surface . We believe that a similar situation is likely to happen in NbSe$_3$. Above $T_{2b}$, when the $\bm{q_2}$ CDW fluctuations become 2D, these were shown to be transverse in the $(\bm{a},\bm{b})$ plane and not in the $(\bm{b},\bm{c})$ plane . This is consistent with x-ray refinements of the $\bm{q_2}$ superlattice modulation at 4 K, showing that larger lattice displacements occur along $\bm{a}$ than along $\bm{c}$ . Thus, the phonon modes involved in the $\bm{q_2}$ distortion possess a transverse component at the surface, supporting the idea of softer transverse phonon modes for the top NbSe$_3$ layer, which would lead to a CDW ordering at a higher critical temperature through an increased electron-phonon coupling.

We get a deeper insight on the surface transition nature by measuring the in-plane correlation functions and extracting the inverse correlation lengths (see Figs. c and d). We observe a continuous evolution with $T$ showing that the system is in a 2D critical regime between 88 and 62 K. At 88 K, the correlations lengths $\xi_{b}$ (along $\bm{b}$) and $\xi_{c}$ (along $\bm{c}$) are equal to 25$b$ and 3$c$ respectively and increase to 70$b$ and 10$c$ at T$\approx$`<!-- -->`{=html}62K. Below 62 K both correlation lengths saturate, whereas the $\bm{{q_{2p}}}$ Fourier amplitude is stabilized according to Fig. b. Here, the 2D regime connects with the bulk pre-transitional region $T_{2b}=59$ K$<T<63$ K of 3D ordering  - which will lead to a growing interaction with the bulk.

To analyze the experimental results further we should take into account some limitations regarding the length scale, a finite window in the STM experiment, as well as regarding the time scale, which is not instantaneous like e.g. in X-ray diffraction experiment. A window of size $L\times L$ with $L=30$ nm, comprises 85 lattice units - $20$ CDW periods - along $\bm{b}$ and $20$ lattice units - $10$ CDW periods - along $\bm{c}$. We see from Fig. d that $L=30$ nm determines the smallest measurable FWHM which are reached at $\approx62$ K. A useful scale to compare with is the expected correlation length at $T=70$ K due to the thermal fluctuations of an isolated chain, $\xi_{T}=\hbar v_{F}/(\pi k_BT)\approx10$Å$\approx3b$, with $v_{F}=3\cdot10^{7}cm/\sec$. This value is an order of magnitude smaller than the observed one, supporting the presence of the observed 2D regime. We infer from this analysis of the correlation lengths that at $62$ K appears a ghost transition because of i) the finite size of the observation window, ii) the beginning of 3D correlations in the bulk which suppresses the surface 2D fluctuations.

The BKT transition refers to 2D systems like the He4, a superconductor or an incommensurate CDW, described by a complex order parameter. The BKT scheme contains two features; i) below $T_{BKT}$ the order parameter correlation falls off as a power law $\sim r^{-2\eta}$; ii) approaching $T_{BKT}$ the correlation length grows very fast, the nature of the transition being the confinement of vortex/anti-vortex pairs which may become unbound at $T>T_{BKT}$. In principle this description applies to the regime where the instantaneous local order parameter amplitude $A_{t,1D}$ (unlike the coherent component $A_{2D} \ll A_{t,1D}$ measured by the 2D FT) is already fixed, i.e. at low $T$ starting from a close vicinity of a mean-field (MF) transition at $T_{MF}$, which usually limits the observations. This region may be extended to higher $T$ in anisotropic systems like the CDWs, where the planes consist of chains. The 1D-2D-3D hierarchy promotes the regime with phase-only fluctuations: the 1D regime coalesces directly into the BKT state . In our case this facilitates the observation of the wide 2D fluctuating regime.

The lack of true long range order in space is ultimately related to its lack also in time. In the 1D regime, i.e. for $T>T_{BKT}$ (above $\approx90$ K), the estimated fluctuation time scale is orders of magnitude smaller than the typical STM acquisition time of 1 sec per line ; hence the temporal fluctuations wash out the CDW from observations. For $T<T_{BKT}$, the tremendous slowing down of the fluctuations allows for the order observation (the local order time decay law $\sim t^{-\eta}$ is characterized by the very small $\eta\le 1/8$). This discussion leads us to a usually overlooked fundamental distinction between the slow STM measurement and the instantaneous one by X-rays. Grazing incidence x-ray diffraction would be suitable tools to compare with our conclusions by tracking the $\bm{q_2}$ transition at the surface and in the bulk. Finally we mention that some of our scans allow to visualize directly the vortices (dislocations in the CDW language) and their pairs (phase solitons) which strengthens the advocated BKT scenario.

In conclusion, high-resolution STM images of NbSe$_3$ show that the $\bm{q_2}$ CDW, unlike the $\bm{q_1}$ one, undergoes a sharp "extraordinary" transition at the surface. Our analysis of the correlation lengths shows at $\approx88$ K a transition from the 1D regime to the 2D BKT state. This seems to be the first real space observation of the BKT regime. Below $\approx63$ K, 2D correlations are stabilized by the critical increase of 3D correlations in the bulk. When probing the fluctuations of the local CDW order by STM, our analysis emphasizes that the time decay of the local order dominates over the space one. This conclusion may have impact on other STM studies in strongly fluctuating systems like stripes in doped oxides.

The high-quality NbSe$_3$ samples were synthesized by H. Berger and F. Lévy. We thank J.-C. Girard, E. Machado and E. Canadell for many helpful discussions and C. David for technical assistance. This work was partly supported by the ANR grant BLAN07-03-192276.


abstract: Blockchain and Cryptocurrencies are gaining unprecedented popularity and understanding. Meanwhile, Ethereum is gaining a significant popularity in the blockchain community, mainly due to the fact that it is designed in a way that enables developers to write smart contract and decentralized applications (Dapps). This new paradigm of applications opens the door to many possibilities and opportunities. However, the security of Ethereum smart contracts has not received much attention; several Ethereum smart contracts malfunctioning have recently been reported. Unlike many previous works that have applied static and dynamic analyses to find bugs in smart contracts, we do not attempt to define and extract any features; instead we focus on reducing the expert's labor costs. We first present a new in-depth analysis of potential attacks methodology and then translate the bytecode of solidity into RGB color code. After that, we transform them to a fixed-sized encoded image. Finally, the encoded image is fed to convolutional neural network (CNN) for automatic feature extraction and learning, detecting compiler bugs of Ethereum smart contract.
author: TonTon Hsien-De Huang
Leopard Mobile Inc. (Cheetah Mobile Taiwan Agency)
Department of Computer Science and Information Engineering, National Cheng Kung University
firstname.lastname@example.com
title: Hunting the Ethereum Smart Contract: Color-inspired Inspection of Potential Attacks

# Introduction

Since Satoshi Nakamoto published the article "Bitcoin: A Peer-to-Peer Electronic Cash System" in 2008 , and after the official launch of Bitcoin in 2009, technologies such as blockchain and cryptocurrency have attracted attentions from academia and industry. At present, the technologies have been applied to many fields such as medical science, economics, Internet of Things . Since the launch of Ethereum (Next Generation Encryption Platform) with Vital contract function proposed by Vitalik Buterin in 2015, lots of attention has been obtained on its dedicated cryptocurrency Ether, smart contract, blockchain and its decentralized Ethereum Virtual Machine (EVM). The main reason is that its design method provides developers with the ability to develop Decentralized apps (Dapps), and thus obtain wider applications. A new application paradigm opens the door to many possibilities and opportunities.

Smart contracts are stored in blocks in order to assist in the negotiation and implementation of the verification contract. The New York Times stated that it distributes and pays the use of public computers that constitute the Internet by many users through the use of Ethereum ; Bloomberg Businessweek claims that it is a software that is shared to everyone and cannot be tampered . Smart contracts can interact with the database at a low cost, resulting in a decentralized autonomous organization made up of Ethereum; at the same time, because it is transparent, the smart contract can prove that its claimed function is real, eg. virtual casinos can prove its fairness . However, the security of smart contracts has not received the same attention, and because of its openness, if there are vulnerability in the contract, anyone can immediately see it, but it cannot be stopped immediately. The DAO vulnerabilities that worth more than \$50 million US dollars is an example . And it raises a similar error through static analysis. There are many vulnerabilities in the Ethereum smart contract, including "Call Stack Attack Vulnerability" , "Time Dependance Vulnerability" and various other known vulnerabilities . These development defects caused thousands of dollars in losses. Compared to traditional programming which can be patched if an error is detected, the smart contract on the blockchain is irreversible and immutable; regardless of its popularity or how much it is worth, the vulnerability cannot be repaired. Hence it is important to further determine the correctness of smart contracts before deployment.

At present, the most widely used security analysis technology is static analysis. However, static analysis needs to be preceded by expert manpower through the analysis of the source code and the acquisition of key code execution sections as a feature to check whether they may result in malicious trading behavior. The static analysis method is less applicable because of the irreversible and immutable nature of the development of smart contracts. Therefore, we have developed a color representation for translating the bytecode of solidity (Ethereum smart contract develop language) into RGB color code and transform them to a fixed-sized encoded image. After that, we also implement a Convolutional Neural Network (CNN) model to perform single-label classification of solidity bytecode images (without extracting features from the solidity source code manually in advance). The models we have tried (AlexNet, GoogleNet, and Inception-v3) have achieved good results in our experiments in addition to the good results achieved with the ImageNet Large-Scale Visual Recognition Challenge (ILSVRC); however, the solidity compilation vulnerability usually does not happen alone, so we further optimize it by using multi-label classification of solidity bytecode images through Transfer Learning, which takes only 1.5 seconds per analysis and obtains more accurate multi-label classification results. Because of the inability of smart contracts to be repaired after deployment, we hope that with the tools we have developed, we can provide developers with early detection of relevant vulnerabilities to repair and reduce smart contracts without deploying source code before deployment to reduce the security issues.

# Related Work

A very recent study conducted analysis on nearly one million contracts; among them, 34,200 contracts are vulnerable and nearly 4,000 contracts are practically exploitable. Smart contracts can hardly be adjusted once deployed. To achieve secure contracts, the key step is to have a thorough security examination before deployment . Loi Luu et al, considered that it is more accurate than the dynamic analysis if the results of the path-by-path method are inferred by static analysis, especially for Ethereum. The uncertainty and complexity of the blockchain behavior makes it difficult to simulate the execution environment. It is the biggest issue for dynamic analysis. Therefore, they developed an execution tool for symbols, called Oyente, to discover potential vulnerability in Ethereum's smart contracts.

Each symbolic path has a path condition which is a formula over the symbolic inputs built by accumulating constraints which those inputs must satisfy in order for execution to follow that path. A path is infeasible if its path condition is unsatiable. Otherwise, the path is feasible. Among 19366 existing Ethereum smart contracts, Oyente flags 8833 of them as vulnerable, including the The DAO which led to a \$50 million US dollar loss in June 2016 . IBM Research presents a framework (ZEUS) to verify the correctness and validate the fairness of smart contracts. ZEUS abstract interpretation and symbolic model checking, along with the power of constrained horn clauses to quickly verify contracts for safety . At the same time, they have analyzed almost 22.4K smart contracts,and discover 94.6% of contracts (containing cryptocurrency worth more than \$0.5 billion) are vulnerable. Microsoft Research, Inria, and Harvard University also have developed a dependent types and monadic effects framework F\* (a functional programming language aimed at program verification), and automated queries to statically verify properties on EVM bytecode and Solidity sources .

In summary, currently, most of the Ethereum security research still rely on labor-intensive examination and focus only on analyzing the control flow graph (CFG) or symbolic execution of the smart contract. so as to determine whether the program under test is causing malicious transaction behavior. This is insufficient in identifying the security flaws in codes. However, Machine Learning/Deep Learning (ML/DL) already has a wide range of applications, especially in security problems, such as spam filtering, botnet detection, and malware classification. Mean in while our proposed system is designed particularly for the security examination of contracts with the minimum labor cost.

# Our Proposed Methodology

We have developed a color representation for translating the bytecode of solidity into RGB color code and transform them to a fixed-sized encoded image. After that, the encoded image is fed to convolutional neural network for automatic feature extraction and learning (without extracting features from the solidity source code manually in advance.), reducing the expert's labor costs. Such translation is also featured by the fact that more complex information in the solidity source code can be preserved in the color image with 16777216 colors (each sampling with 24 bit pixels) compared to the gray scale image with only 256 colors (each sampling with 8 bit pixels). With the fully connected network infrastructure of DNN, though it can deal with its large amount of parameters, however, the local receptive fields and shared weights of CNN make it more suited for more complex structure. It not only decreases the amount of parameters, but also reflects the complexity of smart contract, saving the time for huge computation with current method.

## The Core Technology of Our Methodology

We now explain further in details. Firstly, taking the smart contract of The DAO (0xBB9bc244D798123fDe783fCc1C72d3Bb8C189413) as an example, its bytecode can be obtained in "https://etherscan.io" (as shown in Fig. ). After that, we perform translating the bytecode of solidity into RGB color code (eg 606060 =(R:96, G:96, B:96), 405260 = (R:64, G:82, B:96), 00357c = ( R: 64, G: 81, B: 96)) and then convert bytecode to rgb color code. After this, we can get color images. Then we input these images to CNN for training the compiler bugs detection model for smart contracts. Fig. is an example of converting the bytecode of solidity to jpg.

According to our experience, we also found that two approaches might be used to escape our smart contract detection:

-   Since the traditional filter size of CNN is 3\*3 or 5\*5, the uncorrelated bytecode might become correlated when we transform smart contract into images. The compiler bugs of smart contract may evade the detection by taking advantage of such a mismatch.

-   Meanwhile, smart contract color images are not natural images; instead, they are formed from solidity source code. Thus, the pooling inevitably destroys the contexts and semantics of the malware code, causing the detection inaccuracy.

To address the above two issues, we did many experiments with CNN models (includes AlexNet, GoogleNet, and Inception-v3). We found the characteristics of 1x1 convolution in "Network in Network " . 1x1 convolution is equivalent to cross-channel parametric pooling layer, and this cascaded cross channel parametric pooling structure allows complex and learnable interactions of cross channel information.

## The Architecture of Our Methodology

Fig. is a screenshot of our internal proof-of-concept UI, and Fig. is the result of the passback of our RESTful API (multi-label). A brief description of our system flow chart is shown as follows (steps 1-4 are off-line phases for our internal development and steps 5-6 are online phases, where the user can interact with the system. More specifically, after uploading the bytecode and invokes our RESTful API, the user can obtain the analysis result.):

-   Step 1. Crawling the bytecode of smart contract from etherscan by pyspider (including benign and malicious ones) as sample;

-   Step 2. Transforming the bytecode of smart contract into RGB color code and transforming them to a fixed-sized encoded image;

-   Step 3. After that, the encoded image is fed to convolutional neural network for automatic feature extraction and learning (without extracting features from the solidity source code manually in advance.)

-   Step 4. Finally, once the model has been trained and validated, we deploy it on the backend server.

-   Step 5. Only provide the bytecode of smart contract, we will transform it into smart contract color image.

-   Step 6. After the bytecode of smart contract are all identified, the scanned results will be provided to the users through our visualization tool and public RESTful API.

However, we also found more key problems. Because the gap between the number of normal contract and the number of vulnerable contract samples is very large, and currently there is only 17 Solidity known bugs such as "optimizerStateKnowledgeNotResetForJumpdest", "ArrayAccessCleanHigherOrderBits", and " AncientCompiler" etc., and their gap is also very large. We need to collect those imbalanced data and we need sufficient training data in this domain. It means that a picture will be defined in multiple categories. Therefore, based on our previous research , we re-implement Inception-v3 to multi-label classification through Transfer Learning. For example, 0xbBCf10D6bc180172d8d352BE5bBCfB814E8f3474 is at the same time with Solidity bugs including SolidFunSelectSelector, DelegateCallReturnValue, ECRecoverMalformedInput and SkipEmptyStringLiteral.

# Experiment Result

We run it on a 64-bit Ubuntu 14.04, and hardware setting are 128 GB DDR4 2400 RAM and Intel(R) Xeon(R) E5-2620 v4 CPU, NVIDIA TITAN V, TITAN XP and GTX 1080 GPUs; more specifically, the software setting is the nvidia-docker tensorflow:18.04-py2 on NVIDIA cloud. The research results and the data will be found on our website http://R2D2.TWMAN.ORG.

We first crawled from the etherscan.io for the verified contract information from Jan. 2018 to Apr. 2018. The smart contract message on this website represents that the ethereum mainnet has been verified. We further translate them into color images and feed them to CNN for training.

-   AlexNet

    -   100 epoch, LR: 0.01, accuracy: 83.85%, loss: 0.34

    -   250 epoch, LR: 0.01, accuracy: 86%, loss: 0.39

-   GoogleNet

    -   100 epoch, LR: 0.01, accuracy: 86%, loss: 0.32

    -   250 epoch, LR: 0.01, accuracy: 90%, loss: 0.53

-   Inception-v3

    -   100 epoch, LR: 0.001, accuracy: 97.10%, loss: 0.1345

    -   100 epoch, LR: 0.0001, accuracy: 96.52%, loss: 0.0935

    -   500 epoch, LR: 0.001, accuracy: 97.39%, loss: 0.1876

    -   500 epoch, LR: 0.0001, accuracy: 83.76%, loss: 1.08445

After that, we collected from May. 2018 to Jun 2018 from etherscan.io as our test dataset with Inception-v3. We calculate metrics such as accuracy, precision and recall and show the results as follows.

-   100 epoch, LR 0.001, accuracy: 95.44%, precision: 95.42%, recall: 98.5%

-   100 epoch, LR 0.0001, accuracy: 95.61%, precision: 95.26%, recall: 98.9%

-   500 epoch, LR 0.001, accuracy: 95.85%, precision: 95.44%, recall: 99.04%

-   500 epoch, LR 0.0001, accuracy: 83.65%, precision: 86.48%, recall: 91.46%

On the other hand, since the solidity compile bugs of smart contracts have the characteristics of multi-label, we initially identify if there are compile bugs. If the analysis results are malicious, we use transfer learning for multi-label. After training and inference, the results obtained are shown in Figure 4. Fig. is a screenshot of our public proof-of-concept scan result.

# Conclusion

According to our data, there are 1800 new smart contract produced on ethereum main net per day. Amongst less than 30% is verified by etherscan. Meanwhile according to the third party blockchain evaluation system, Rating Token (https://ratingtoken.io) and coinschedule, in 2018, there has been 11.75 billion US dollars financed by ICO projects. One of the key factors for the success of a smart contract for each project is the existence of a loophole; Our goal is to optimize the amount of parameters, network structure, and release automated verification tools and public RESTful API. The above experiment results demonstrate that our proposed system can have accurate security analysis on contracts with very limited labor cost.