Search

We address the problem of verifying neural networks against geometric transformations of the input image, including rotation, scaling, shearing, and translation. The proposed method computes provably sound piecewise linear constraints for the pixel values by using sampling and linear approximations in combination with branch-and-bound Lipschitz optimisation. The method obtains provably tighter over-approximations of the perturbation region than the present state-of-the-art. We report results from experiments on a comprehensive set of verification benchmarks on MNIST and CIFAR10. We show that our proposed implementation resolves up to 32% more verification cases than present approaches., Comment: ECAI 2024

The application of large-language models (LLMs) to digital hardware code generation is an emerging field. Most LLMs are primarily trained on natural language and software code. Hardware code, such as Verilog, represents only a small portion of the training data and few hardware benchmarks exist. To address this gap, the open-source VerilogEval benchmark was released in 2023, providing a consistent evaluation framework for LLMs on code completion tasks. It was tested on state-of-the-art models at the time including GPT-4. However, VerilogEval and other Verilog generation benchmarks lack failure analysis and, in present form, are not conducive to exploring prompting techniques. Also, since VerilogEval's release, both commercial and open-source models have seen continued development. In this work, we evaluate new commercial and open-source models of varying sizes against an improved VerilogEval benchmark suite. We enhance VerilogEval's infrastructure and dataset by automatically classifying failures, introduce new prompts for supporting in-context learning (ICL) examples, and extend the supported tasks to specification-to-RTL translation. We find a measurable improvement in commercial state-of-the-art models, with GPT-4 Turbo achieving a 59% pass rate on spec-to-RTL tasks. We also study the performance of open-source and domain-specific models that have emerged, and demonstrate that models can benefit substantially from ICL. We find that recently-released Llama 3.1 405B achieves a pass rate of 58%, effectively matching that of GPT-4 Turbo, and that the much smaller domain-specific RTL-Coder 6.7B models achieve an impressive 37% pass rate. However, prompt engineering is key to achieving good pass rates, and varies widely with model and task. A benchmark infrastructure that allows for prompt engineering and failure analysis is key to continued model development and deployment., Comment: This paper revisits and improves the benchmark first presented in arXiv:2309.07544. Seven pages, three figures

DNA sequence alignment is an important workload in computational genomics. Reference-guided DNA assembly involves aligning many read sequences against candidate locations in a long reference genome. To reduce the computational load of this alignment, candidate locations can be pre-filtered using simpler alignment algorithms like edit distance. Prior work has explored accelerating filtering on simulated compute-in-DRAM, due to the massive parallelism of compute-in-memory architectures. In this paper, we present work-in-progress on accelerating filtering using a commercial compute-in-SRAM accelerator. We leverage the recently released Gemini accelerator platform from GSI Technology, which is the first, to our knowledge, commercial-scale compute-in-SRAM system. We accelerate the Myers' bit-parallel edit distance algorithm, producing average speedups of 14.1x over single-core CPU performance. Individual query/candidate alignments produce speedups of up to 24.1x. These early results suggest this novel architecture is well-suited to accelerating the filtering step of sequence-to-sequence DNA alignment.

We introduce two algorithms for computing tight guarantees on the probabilistic robustness of Bayesian Neural Networks (BNNs). Computing robustness guarantees for BNNs is a significantly more challenging task than verifying the robustness of standard Neural Networks (NNs) because it requires searching the parameters' space for safe weights. Moreover, tight and complete approaches for the verification of standard NNs, such as those based on Mixed-Integer Linear Programming (MILP), cannot be directly used for the verification of BNNs because of the polynomial terms resulting from the consecutive multiplication of variables encoding the weights. Our algorithms efficiently and effectively search the parameters' space for safe weights by using iterative expansion and the network's gradient and can be used with any verification algorithm of choice for BNNs. In addition to proving that our algorithms compute tighter bounds than the SoA, we also evaluate our algorithms against the SoA on standard benchmarks, such as MNIST and CIFAR10, showing that our algorithms compute bounds up to 40% tighter than the SoA., Comment: Accepted at AISTATS 2024

As the pace of progress that has followed Moore's law continues to diminish, it is critical that the US support Integrated Circuit (IC or chip) education and research to maintain technological innovation. Furthermore, US economic independence, security, and future international standing rely on having on-shore IC design capabilities. New devices with disparate technologies, improved design software toolchains and methodologies, and technologies to integrate heterogeneous systems will be needed to advance IC design capabilities. This will require rethinking both how we teach design to address the new complexity and how we inspire student interest in a hardware systems career path. The main recommendation of this workshop is that accessibility is the key issue. To this end, a National Chip Design Center (NCDC) should be established to further research and education by partnering academics and industry to train our future workforce. This should not be limited to R1 universities, but should also include R2, community college, minority serving institutions (MSI), and K-12 institutions to have the broadest effect. The NCDC should support the access, development, and maintenance of open design tools, tool flows, design kits, design components, and educational materials. Open-source options should be emphasized wherever possible to maximize accessibility. The NCDC should also provide access and support for chip fabrication, packaging and testing for both research and educational purposes., Comment: This material is based upon work supported by the NSF under Grant No. 2137629

There is a growing interest in single-class modelling and out-of-distribution detection as fully supervised machine learning models cannot reliably identify classes not included in their training. The long tail of infinitely many out-of-distribution classes in real-world scenarios, e.g., for screening, triage, and quality control, means that it is often necessary to train single-class models that represent an expected feature distribution, e.g., from only strictly healthy volunteer data. Conventional supervised machine learning would require the collection of datasets that contain enough samples of all possible diseases in every imaging modality, which is not realistic. Self-supervised learning methods with synthetic anomalies are currently amongst the most promising approaches, alongside generative auto-encoders that analyse the residual reconstruction error. However, all methods suffer from a lack of structured validation, which makes calibration for deployment difficult and dataset-dependant. Our method alleviates this by making use of multiple visually-distinct synthetic anomaly learning tasks for both training and validation. This enables more robust training and generalisation. With our approach we can readily outperform state-of-the-art methods, which we demonstrate on exemplars in brain MRI and chest X-rays. Code is available at https://github.com/matt-baugh/many-tasks-make-light-work ., Comment: Early accepted to MICCAI 2023

This technical report outlines our submission to the zero-shot track of the Visual Anomaly and Novelty Detection (VAND) 2023 Challenge. Building on the performance of the WINCLIP framework, we aim to enhance the system's localization capabilities by integrating zero-shot segmentation models. In addition, we perform foreground instance segmentation which enables the model to focus on the relevant parts of the image, thus allowing the models to better identify small or subtle deviations. Our pipeline requires no external data or information, allowing for it to be directly applied to new datasets. Our team (Variance Vigilance Vanguard) ranked third in the zero-shot track of the VAND challenge, and achieve an average F1-max score of 81.5/24.2 at a sample/pixel level on the VisA dataset., Comment: Ranked 3rd in zero-shot track of the Visual Anomaly and Novelty Detection (VAND) 2023 Challenge

Immunohistochemistry (IHC) biomarkers are essential tools for reliable cancer diagnosis and subtyping. It requires cross-staining comparison among Whole Slide Images (WSIs) of IHCs and hematoxylin and eosin (H&E) slides. Currently, pathologists examine the visually co-localized areas across IHC and H&E glass slides for a final diagnosis, which is a tedious and challenging task. Moreover, visually inspecting different IHC slides back and forth to analyze local co-expressions is inherently subjective and prone to error, even when carried out by experienced pathologists. Relying on digital pathology, we propose Composite Biomarker Image (CBI) in this work. CBI is a single image that can be composed using different filtered IHC biomarker images for better visualization. We present a CBI image produced in two steps by the proposed solution for better visualization and hence more efficient clinical workflow. In the first step, IHC biomarker images are aligned with the H&E images using one coordinate system and orientation. In the second step, the positive or negative IHC regions from each biomarker image (based on the pathologists recommendation) are filtered and combined into one image using a fuzzy inference system. For evaluation, the resulting CBI images, from the proposed system, were evaluated qualitatively by the expert pathologists. The CBI concept helps the pathologists to identify the suspected target tissues more easily, which could be further assessed by examining the actual WSIs at the same suspected regions.

Extracting complex structures from grid-based data is a common key step in automated medical image analysis. The conventional solution to recovering tree-structured geometries typically involves computing the minimal cost path through intermediate representations derived from segmentation masks. However, this methodology has significant limitations in the context of projective imaging of tree-structured 3D anatomical data such as coronary arteries, since there are often overlapping branches in the 2D projection. In this work, we propose a novel approach to predicting tree connectivity structure which reformulates the task as an optimization problem over individual steps of a recursive process. We design and train a two-stage model which leverages the UNet and Transformer architectures and introduces an image-based prompting technique. Our proposed method achieves compelling results on a pair of synthetic datasets, and outperforms a shortest-path baseline., Comment: 12 pages, 5 figures

We study in detail the sensitivity of the Antarctic Impulsive Transient Antenna (ANITA) to possible $\nu_\tau$ point source fluxes detected via $\tau$-lepton-induced air showers. This investigation is framed around the observation of four upward-going extensive air shower events very close to the horizon seen in ANITA-IV. We find that these four upgoing events are not observationally inconsistent with $\tau$-induced EASs from Earth-skimming $\nu_\tau$, both in their spectral properties as well as in their observed locations on the sky. These four events, as well as the overall diffuse and point source exposure to Earth-skimming $\nu_\tau$, are also compared against published ultrahigh-energy neutrino limits from the Pierre Auger Observatory. While none of these four events occurred at sky locations simultaneously visible by Auger, the implied fluence necessary for ANITA to observe these events is in strong tension with limits set by Auger across a wide range of energies and is additionally in tension with ANITA's Askaryan in-ice neutrino channel above $10^{19}$ eV. We conclude by discussing some of the technical challenges with simulating and analyzing these near horizon events and the potential for future observatories to observe similar events., Comment: 19 pages, 22 figures, will be published in Physical Review D (PRD)

Fast Radio Bursts (FRBs) are a promising tool for studying the low-density universe as their dispersion measures (DM) are extremely sensitive probes of electron column density. Active Galactic Nuclei (AGN) inject energy into the intergalactic medium, affecting the DM and its scatter. To determine the effectiveness of FRBs as a probe of AGN feedback, we analysed three different AGN models from the EAGLE simulation series. We measured the mean DM-redshift relation, and the scatter around it, using $2.56 \times 10^8$ sightlines at 131 redshift ($z$) bins between $0 \leq z \leq 3$. While the DM-redshift relation itself is highly robust against different AGN feedback models, significant differences are detected in the scatter around the mean: weaker feedback leads to more scatter. We find $\sim 10^4$ localised FRBs are needed to discriminate between the scatter in standard feedback and stronger, more intermittent feedback models. The number of FRBs required is dependent on the redshift distribution of the detected population. A log-normal redshift distribution at $z=0.5$ requires approximately 50% fewer localised FRBs than a distribution centred at $z=1$. With the Square Kilometre Array expected to detect $>10^3$ FRBs per day, in the future, FRBs will be able to provide constraints on AGN feedback., Comment: 5 pages, 4 figures; Submitted to MNRAS Letters. Comments welcome

The dispersion measures (DM) of fast radio bursts (FRBs) encode the integrated electron density along the line-of-sight, which is typically dominated by the intergalactic medium (IGM) contribution in the case of extragalactic FRBs. In this paper, we show that incorporating wide-field spectroscopic galaxy survey data in the foreground of localized FRBs can significantly improve constraints on the partition of diffuse cosmic baryons. Using mock DMs and realistic lightcone galaxy catalogs derived from the Millennium simulation, we define spectroscopic surveys that can be carried out with 4m and 8m-class wide field spectroscopic facilities. On these simulated surveys, we carry out Bayesian density reconstructions in order to estimate the foreground matter density field. In comparison with the `true' matter density field, we show that these can help reduce the uncertainties in the foreground structures by $\sim 2-3\times$ compared to cosmic variance. We calculate the Fisher matrix to forecast that $N=30\: (96)$ localized FRBs should be able to constrain the diffuse cosmic baryon fraction to $\sim 10\%\: (\sim 5\%) $, and parameters governing the size and baryon fraction of galaxy circumgalactic halos to within $\sim 20-25\%\: (\sim 8-12\%)$. From the Fisher analysis, we show that the foreground data increases the sensitivity of localized FRBs toward our parameters of interest by $\sim 25\times$. We briefly introduce FLIMFLAM, an ongoing galaxy redshift survey that aims to obtain foreground data on $\sim 30$ localized FRB fields., Comment: 27 pages, 12 Figures. Accepted by ApJ

Joint analysis of multiple biomarker images and tissue morphology is important for disease diagnosis, treatment planning and drug development. It requires cross-staining comparison among Whole Slide Images (WSIs) of immuno-histochemical and hematoxylin and eosin (H&E) microscopic slides. However, automatic, and fast cross-staining alignment of enormous gigapixel WSIs at single-cell precision is challenging. In addition to morphological deformations introduced during slide preparation, there are large variations in cell appearance and tissue morphology across different staining. In this paper, we propose a two-step automatic feature-based cross-staining WSI alignment to assist localization of even tiny metastatic foci in the assessment of lymph node. Image pairs were aligned allowing for translation, rotation, and scaling. The registration was performed automatically by first detecting landmarks in both images, using the scale-invariant image transform (SIFT), followed by the fast sample consensus (FSC) protocol for finding point correspondences and finally aligned the images. The Registration results were evaluated using both visual and quantitative criteria using the Jaccard index. The average Jaccard similarity index of the results produced by the proposed system is 0.942 when compared with the manual registration., Comment: Accepted in EMBC 2021 : 43rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society

Deep learning models for semantic segmentation are able to learn powerful representations for pixel-wise predictions, but are sensitive to noise at test time and do not guarantee a plausible topology. Image registration models on the other hand are able to warp known topologies to target images as a means of segmentation, but typically require large amounts of training data, and have not widely been benchmarked against pixel-wise segmentation models. We propose Atlas-ISTN, a framework that jointly learns segmentation and registration on 2D and 3D image data, and constructs a population-derived atlas in the process. Atlas-ISTN learns to segment multiple structures of interest and to register the constructed, topologically consistent atlas labelmap to an intermediate pixel-wise segmentation. Additionally, Atlas-ISTN allows for test time refinement of the model's parameters to optimize the alignment of the atlas labelmap to an intermediate pixel-wise segmentation. This process both mitigates for noise in the target image that can result in spurious pixel-wise predictions, as well as improves upon the one-pass prediction of the model. Benefits of the Atlas-ISTN framework are demonstrated qualitatively and quantitatively on 2D synthetic data and 3D cardiac computed tomography and brain magnetic resonance image data, out-performing both segmentation and registration baseline models. Atlas-ISTN also provides inter-subject correspondence of the structures of interest, enabling population-level shape and motion analysis., Comment: 33 pages, 15 figures

The dispersion measure (DM) of fast radio bursts (FRBs) provides a unique way to probe ionised baryons in the intergalactic medium (IGM). Cosmological models with different parameters lead to different DM-redshift ($\mathrm{DM}-z$) relations. Additionally, the over/under-dense regions in the IGM and the circumgalactic medium of intervening galaxies lead to scatter around the mean $\mathrm{DM}-z$ relations. We have used the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations to measure the mean $\mathrm{DM}-z$ relation and the scatter around it using over one billion lines-of-sight between redshifts $0

In medical imaging, outliers can contain hypo/hyper-intensities, minor deformations, or completely altered anatomy. To detect these irregularities it is helpful to learn the features present in both normal and abnormal images. However this is difficult because of the wide range of possible abnormalities and also the number of ways that normal anatomy can vary naturally. As such, we leverage the natural variations in normal anatomy to create a range of synthetic abnormalities. Specifically, the same patch region is extracted from two independent samples and replaced with an interpolation between both patches. The interpolation factor, patch size, and patch location are randomly sampled from uniform distributions. A wide residual encoder decoder is trained to give a pixel-wise prediction of the patch and its interpolation factor. This encourages the network to learn what features to expect normally and to identify where foreign patterns have been introduced. The estimate of the interpolation factor lends itself nicely to the derivation of an outlier score. Meanwhile the pixel-wise output allows for pixel- and subject- level predictions using the same model., Comment: Accepted for publication at the Journal of Machine Learning for Biomedical Imaging (MELBA) https://www.melba-journal.org

The ANtarctic Impulsive Transient Antenna (ANITA) long-duration balloon experiment is sensitive to interactions of ultra high-energy (E > 10^{18} eV) neutrinos in the Antarctic ice sheet. The third flight of ANITA, lasting 22 days, began in December 2014. We develop a methodology to search for energetic neutrinos spatially and temporally coincident with potential source classes in ANITA data. This methodology is applied to several source classes: the TXS 0506+056 blazar and NGC 1068, the first potential TeV neutrino sources identified by IceCube, flaring high-energy blazars reported by the Fermi All-Sky Variability Analysis, gamma-ray bursts, and supernovae. Among searches within the five source classes, one candidate was identified as associated with SN 2015D, although not at a statistically significant level. We proceed to place upper limits on the source classes. We further comment on potential applications of this methodology to more sensitive future instruments., Comment: 23 pages, 7 figures, version accepted to JCAP

The balloon-borne ANITA experiment is designed to detect ultra-high energy neutrinos via radio emissions produced by an in-ice shower. Although initially purposed for interactions within the Antarctic ice sheet, ANITA also demonstrated the ability to self-trigger on radio emissions from ultra-high energy charged cosmic rays interacting in the Earth's atmosphere. For showers produced above the Antarctic ice sheet, reflection of the down-coming radio signals at the Antarctic surface should result in a polarity inversion prior to subsequent observation at the $\sim$35-40 km altitude ANITA gondola. ANITA has published two anomalous instances of upcoming cosmic-rays with measured polarity opposite the remaining sample of $\sim$50 UHECR signals. The steep observed upwards incidence angles (25--30 degrees relative to the horizontal) require non-Standard Model physics if these events are due to in-ice neutrino interactions, as the Standard Model cross-section would otherwise prohibit neutrinos from penetrating the long required chord of Earth. Shoemaker et al. posit that glaciological effects may explain the steep observed anomalous events. We herein consider the scenarios offered by Shoemaker et al. and find them to be disfavored by extant ANITA and HiCal experimental data. We note that the recent report of four additional near-horizon anomalous ANITA-4 events, at $>3\sigma$ significance, are incompatible with their model, which requires significant signal transmission into the ice.

ANITA's fourth long-duration balloon flight in late 2016 detected 29 cosmic-ray (CR)-like events on a background of $0.37^{+0.27}_{-0.17}$ anthropogenic events. CRs are mainly seen in reflection off the Antarctic ice sheets, creating a characteristic phase-inverted waveform polarity. However, four of the below-horizon CR-like events show anomalous non-inverted polarity, a $p = 5.3 \times 10^{-4}$ chance if due to background. All anomalous events are from locations near the horizon; ANITA-IV observed no steeply-upcoming anomalous events similar to the two such events seen in prior flights., Comment: 6 pages, 4 figures, to appear in Physical Review Letters. Supplemental material (reference 17) available from corresponding author

Information diffusion within financial markets plays a crucial role in the process of price formation and the propagation of sentiment and risk. We perform a comparative analysis of information transfer between industry sectors of the Chinese and the USA stock markets, using daily sector indices for the period from 2000 to 2017. The information flow from one sector to another is measured by the transfer entropy of the daily returns of the two sector indices. We find that the most active sector in information exchange (i.e., the largest total information inflow and outflow) is the {\textit{non-bank financial}} sector in the Chinese market and the {\textit{technology}} sector in the USA market. This is consistent with the role of the non-bank sector in corporate financing in China and the impact of technological innovation in the USA. In each market, the most active sector is also the largest information sink that has the largest information inflow (i.e., inflow minus outflow). In contrast, we identify that the main information source is the {\textit{bank}} sector in the Chinese market and the {\textit{energy}} sector in the USA market. In the case of China, this is due to the importance of net bank lending as a signal of corporate activity and the role of energy pricing in affecting corporate profitability. There are sectors such as the {\textit{real estate}} sector that could be an information sink in one market but an information source in the other, showing the complex behavior of different markets. Overall, these findings show that stock markets are more synchronized, or ordered, during periods of turmoil than during periods of stability., Comment: 12 pages including 8 figures

We combine polarization-resolved low frequency Raman microscopy and dispersion-corrected density-functional calculations (DFT-D3) to study polymorph-dependent lattice vibrations in coronene, a model molecular system for nanographenes and disc-like organic semiconductors that exhibits two crystalline structures with distinct electronic and optical properties. Changes in low energy Raman-active lattice phonons are followed across the {\gamma}- to \b{eta}-phase transition at 150 K. Raman frequencies calculated using DFT-D3 agree to within 4 cm-1, and on the basis of polarisation dependence of peak positions and intensities we achieve a clear mode assignment. Studies of the associated atomic motions show how the pure librational and rotational modes of {\gamma}-coronene change into mixed roto-librations in the \b{eta}-phase, explaining the remarkable differences in Raman spectra of the two phases.

Growth of $>$100 $\mu$m thick diamond layer adherent on aluminium nitride is presented in this work. While thick films failed to adhere on untreated AlN films, hydrogen/nitrogen plasma treated AlN films retained the thick diamond layers. Clear differences in zeta potential measurement confirms the surface modification due to hydrogen/nitrogen plasma treatment. Areal Raman maps showed an increase in non-diamond carbon in the initial layers of diamond grown on pre-treated AlN. The presence of non-diamond carbon has minimal effect on the interface between diamond and AlN. The surfaces studied with x-ray photoelectron spectroscopy (XPS) revealed a clear distinction between pre-treated and untreated samples. The surface aluminium goes from nitrogen rich environment to an oxygen rich environment after pre-treatment. Cross section transmission electron microscopy shows a clean interface between diamond and AlN. Thermal barrier resistance between diamond and AlN was found to be in the range of 16 m$^2$K/GW which is a large improvement on the current state-of-the-art.

This study reports on the plasmon-mediated remote Raman sensing promoted by specially designed coaxial nanowires. This unusual geometry for Raman study is based on the separation, by several micrometres, of the excitation laser spot, on one tip of the nanowire, and the Raman detection at the other tip. The very weak efficiency of Raman emission makes it challenging in a remote configuration. For the proof-of-concept, we designed coaxial nanowires consisting in a gold core to propagate the surface plasmon polaritons and a Raman-emitting shell of poly(3,4-ethylene-dioxythiophene). The success of the fabrication was demonstrated by correlating, for the same single nanowire, a morphological analysis by electron microscopy and Raman spectroscopy analysis. Importantly for probing remote-Raman effect, the original hard template-based process allows to control the location of the polymer shell all along the nanowire, or only close to one or the two nanowire tips. Such all-in-one single nanowires could have applications in the remote detection of photo-degradable substances and for exploring 1D nanosources for integrated photonic and plasmonic systems., Comment: 8 pages, 6 figures, Supplementary Information available

Fruitbat is an open source Python 2/3 package for estimating redshifts, energies and the galactic dispersion measure contributions of fast radio bursts (FRBs). Fruitbat combines various dispersion measure (DM) and redshift relations with the YMW16 galactic dispersion measure model into a single easy to use API., Comment: Accepted by the Journal of Open Source Software, 4 pages, 2 figures, package available via pip and source code available on Github (see text for links)

Understanding collective electronic states such as superconductivity and charge density waves is pivotal for fundamental science and applications. The layered transition metal dichalcogenide 1T-TiSe2 hosts a unique charge density wave (CDW) phase transition whose origins are still not fully understood. Here, we present ultrafast time- and angle-resolved photoemission spectroscopy (TR-ARPES) measurements complemented by time-resolved reflectivity (TRR) which allows us to establish the contribution of excitonic and electron-phonon interactions to the CDW. We monitor the energy shift of the valence band (VB) and coupling to coherent phonons as a function of laser fluence. The VB shift, directly related to the CDW gap closure, exhibits a markedly slower recovery dynamics at fluences above Fth = 60 microJ cm-2. This observation coincides with a shift in the relative weight of coherently coupled phonons to higher frequency modes in time-resolved reflectivity (TRR), suggesting a phonon bottleneck. Using a rate equation model, the emergence of a high-fluence bottleneck is attributed to an abrupt reduction in coupled phonon damping and an increase in exciton dissociation rate linked to the loss of CDW superlattice phonons. Thus, our work establishes the important role of both excitonic and phononic interactions in the CDW phase transition and the advantage of combining complementary femtosecond techniques to understand the complex interactions in quantum materials., Comment: 11 pages, 4 figures

A Monte Carlo simulation program for the radio detection of Ultra High Energy (UHE) neutrino interactions in the Antarctic ice as viewed by the Antarctic Impulsive Transient Antenna (ANITA) is described in this article. The program, icemc, provides an input spectrum of UHE neutrinos, the parametrization of the Askaryan radiation generated by their interaction in the ice, and the propagation of the radiation through ice and air to a simulated model of the third and fourth ANITA flights. This paper provides an overview of the icemc simulation, descriptions of the physics models used and of the ANITA electronics processing chain, data/simulation comparisons to validate the predicted performance, and a summary of the impact of published results.

Raman microscopy is a valuable tool for detecting physical and chemical properties of a sample material. When probing nanomaterials or nanocomposites the spatial resolution of Raman microscopy is not always adequate as it is limited by the optical diffraction limit. Numerical post-processing with super-resolution algorithms provides a means to enhance resolution and can be straightforwardly applied. The aim of this work is to present interior-point least squares (IPLS) as a powerful tool for super-resolution in Raman imaging through constrained optimisation. IPLS's potential for super-resolution is illustrated on numerically generated test images. Its resolving power is demonstrated on Raman spectroscopic data of a polymer nanowire sample. Comparison to AFM data of the same sample substantiates that the presented method is a promising technique for analysing nanomaterial samples., Comment: 21 pages, 8 figures

The ANtarctic Impulsive Transient Antenna (ANITA) NASA long-duration balloon payload completed its fourth flight in December 2016, after 28 days of flight time. ANITA is sensitive to impulsive broadband radio emission from interactions of ultra-high-energy neutrinos in polar ice (Askaryan emission). We present the results of two separate blind analyses searching for signals from Askaryan emission in the data from the fourth flight of ANITA. The more sensitive analysis, with a better expected limit, has a background estimate of $0.64^{+0.69}_{-0.45}$ and an analysis efficiency of $82\pm2\%$. The second analysis has a background estimate of $0.34^{+0.66}_{-0.16}$ and an analysis efficiency of $71\pm6\%$. Each analysis found one event in the signal region, consistent with the background estimate for each analysis. The resulting limit further tightens the constraints on the diffuse flux of ultra-high-energy neutrinos at energies above $10^{19.5}$ eV., Comment: 11 pages, 7 figures

Recently, the ANITA collaboration reported on two upward-going extensive air shower events consistent with a primary particle that emerges from the surface of the ice. These events may be of $\nu_\tau$ origin, in which the neutrino interacts within the Earth to produce a $\tau$ lepton that emerges from the Earth, decays in the atmosphere, and initiates an extensive air shower. In this paper we estimate an upper bound on the ANITA acceptance to a diffuse $\nu_\tau$ flux detected via $\tau$-lepton-induced air showers within the bounds of Standard Model (SM) uncertainties. By comparing this estimate with the acceptance of Pierre Auger Observatory and IceCube and assuming SM interactions, we conclude that a $\nu_\tau$ origin of these events would imply a neutrino flux at least two orders of magnitude above current bounds., Comment: 12 pages, 7 figures

These proceedings address a recent publication by the ANITA collaboration of four upward- pointing cosmic-ray-like events observed in the first flight of ANITA. Three of these events were consistent with stratospheric cosmic-ray air showers where the axis of propagation does not inter- sect the surface of the Earth. The fourth event was consistent with a primary particle that emerges from the surface of the ice suggesting a possible {\tau}-lepton decay as the origin of this event. These proceedings follow-up on the modeling and testing of the hypothesis that this event was of {\tau} neutrino origin., Comment: 8 pages, 3 figures, presented at the International Cosmic Ray Conference 2017, Busan, South Korea

We report on an upward traveling, radio-detected cosmic-ray-like impulsive event with characteristics closely matching an extensive air shower. This event, observed in the third flight of the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload, is consistent with a similar event reported in a previous flight. These events may be produced by the atmospheric decay of an upward-propagating $\tau$-lepton produced by a $\nu_{\tau}$ interaction, although their relatively steep arrival angles create tension with the standard model (SM) neutrino cross section. Each of the two events have $a~posteriori$ background estimates of $\lesssim 10^{-2}$ events. If these are generated by $\tau$-lepton decay, then either the charged-current $\nu_{\tau}$ cross section is suppressed at EeV energies, or the events arise at moments when the peak flux of a transient neutrino source was much larger than the typical expected cosmogenic background neutrinos., Comment: 5 pages, 4 figures. Supplemental material available from corresponding author by request

The Antarctic Impulsive Transient Antenna (ANITA), a NASA long-duration balloon payload, searches for radio emission from interactions of ultra-high-energy neutrinos in polar ice. The third flight of ANITA (ANITA-III) was launched in December 2014 and completed a 22-day flight. We present the results of three analyses searching for Askaryan radio emission of neutrino origin. In the most sensitive of the analyses, we find one event in the signal region on an expected a priori background of $0.7^{+0.5}_{-0.3}$. Though consistent with the background estimate, the candidate event remains compatible with a neutrino hypothesis even after additional post-unblinding scrutiny., Comment: 15 pages, 15 figures, Accepted to PRD

The balloon-borne HiCal radio-frequency (RF) transmitter, in concert with the ANITA radio-frequency receiver array, is designed to measure the Antarctic surface reflectivity in the RF wavelength regime. The amplitude of surface-reflected transmissions from HiCal, registered as triggered events by ANITA, can be compared with the direct transmissions preceding them by O(10) microseconds, to infer the surface power reflection coefficient $\cal{R}$. The first HiCal mission (HiCal-1, Jan. 2015) yielded a sample of 100 such pairs, resulting in estimates of $\cal{R}$ at highly-glancing angles (i.e., zenith angles approaching $90^\circ$), with measured reflectivity for those events which exceeded extant calculations. The HiCal-2 experiment, flying from Dec., 2016-Jan., 2017, provided an improvement by nearly two orders of magnitude in our event statistics, allowing a considerably more precise mapping of the reflectivity over a wider range of incidence angles. We find general agreement between the HiCal-2 reflectivity results and those obtained with the earlier HiCal-1 mission, as well as estimates from Solar reflections in the radio-frequency regime. In parallel, our calculations of expected reflectivity have matured; herein, we use a plane-wave expansion to estimate the reflectivity R from both a flat, smooth surface (and, in so doing, recover the Fresnel reflectivity equations) and also a curved surface. Multiplying our flat-smooth reflectivity by improved Earth curvature and surface roughness corrections now provides significantly better agreement between theory and the HiCal 2a/2b measurements., Comment: submitted to Astropart. Phys

We consider model order reduction of a nonlinear cable-mass system modeled by a 1D wave equation with interior damping and dynamic boundary conditions. The system is driven by a time dependent forcing input to a linear mass-spring system at one boundary. The goal of the model reduction is to produce a low order model that produces an accurate approximation to the displacement and velocity of the mass in the nonlinear mass-spring system at the opposite boundary. We first prove that the linearized and nonlinear unforced systems are well-posed and exponentially stable under certain conditions on the damping parameters, and then consider a balanced truncation method to generate the reduced order model (ROM) of the nonlinear input-output system. Little is known about model reduction of nonlinear input-output systems, and so we present detailed numerical experiments concerning the performance of the nonlinear ROM. We find that the ROM is accurate for many different combinations of model parameters.

The NASA supported High-Altitude Calibration (HiCal)-2 instrument flew as a companion balloon to the ANITA-4 experiment in December 2016. Based on a HV discharge pulser producing radio-frequency (RF) calibration pulses, HiCal-2 comprised two payloads, which flew for a combined 18 days, covering 1.5 revolutions of the Antarctic continent. ANITA-4 captured over 10,000 pulses from HiCal, both direct and reflected from the surface, at distances varying from 100-800 km, providing a large dataset for surface reflectivity measurements. Herein we present details on the design, construction and performance of HiCal-2., Comment: Published in NIM-A, final version

The Antarctic Impulsive Transient Antenna (ANITA) is a NASA long-duration balloon experiment with the primary goal of detecting ultra-high-energy ($>10^{18}\,\mbox{eV}$) neutrinos via the Askaryan Effect. The fourth ANITA mission, ANITA-IV, recently flew from Dec 2 to Dec 29, 2016. For the first time, the Tunable Universal Filter Frontend (TUFF) boards were deployed for mitigation of narrow-band, anthropogenic noise with tunable, switchable notch filters. The TUFF boards also performed second-stage amplification by approximately 45 dB to boost the $\sim\,\mu\mbox{V-level}$ radio frequency (RF) signals to $\sim$ mV-level for digitization, and supplied power via bias tees to the first-stage, antenna-mounted amplifiers. The other major change in signal processing in ANITA-IV is the resurrection of the $90^{\circ}$ hybrids deployed previously in ANITA-I, in the trigger system, although in this paper we focus on the TUFF boards. During the ANITA-IV mission, the TUFF boards were successfully operated throughout the flight. They contributed to a factor of 2.8 higher total instrument livetime on average in ANITA-IV compared to ANITA-III due to reduction of narrow-band, anthropogenic noise before a trigger decision is made., Comment: 27 pages, 19 figures

We report on high resolution measurements of resonances in the spectrum of coherent synchrotron radiation (CSR) at the Canadian Light Source (CLS). The resonances permeate the spectrum at wavenumber intervals of $0.074 ~\textrm{cm}^{-1}$, and are highly stable under changes in the machine setup (energy, bucket filling pattern, CSR in bursting or continuous mode). Analogous resonances were predicted long ago in an idealized theory as eigenmodes of a smooth toroidal vacuum chamber driven by a bunched beam moving on a circular orbit. A corollary of peaks in the spectrum is the presence of pulses in the wakefield of the bunch at well defined spatial intervals. Through experiments and further calculations we elucidate the resonance and wakefield mechanisms in the CLS vacuum chamber, which has a fluted form much different from a smooth torus. The wakefield is observed directly in the 30-110 GHz range by RF diodes, and indirectly by an interferometer in the THz range. The wake pulse sequence found by diodes is less regular than in the toroidal model, and depends on the point of observation, but is accounted for in a simulation of fields in the fluted chamber. Attention is paid to polarization of the observed fields, and possible coherence of fields produced in adjacent bending magnets. Low frequency wakefield production appears to be mainly local in a single bend, but multi-bend effects cannot be excluded entirely, and could play a role in high frequency resonances. New simulation techniques have been developed, which should be invaluable in further work., Comment: 5 figures. Differs from first posting by correction of typo and inclusion of responses to referee's questions

We extend results related to maximal subalgebras and ideals from Lie to Leibniz algebras. In particular, we classify minimal non-elementary Leibniz algebras and Leibniz algebras with a unique maximal ideal. In both cases, there are types of these algebras with no Lie algebra analogue. We also give a classification of E-Leibniz algebras which is very similiar to its Lie algebra counterpart., Comment: 10 pages

We extend results on finite dimensional nilpotent Lie algebras to Leibniz algebras and counterexamples to others are found. One generator algebras are used in these examples and are investigated further.

A Frattini theory for non-associative algebras was developed by Towers and results for particular classes of algebras have appeared in various articles. Especially plentiful are results on Lie algebras. It is the purpose of this paper to extend some of the Lie algebra results to Leibniz algebras.

We generate inherent structures, local potential-energy minima, of the "$k$-space overlap potential" in two-dimensional many-particle systems using a cooling and quenching simulation technique. The ground states associated with the $k$-space overlap potential are stealthy ({\it i.e.,} completely suppress single scattering of radiation for a range of wavelengths) and hyperuniform ({\it i.e.,} infinite wavelength density fluctuations vanish). However, we show via quantitative metrics that the inherent structures exhibit a range of stealthiness and hyperuniformity depending on the fraction of degrees of freedom that are constrained. Inherent structures in two dimensions typically contain five-particle rings, wavy grain boundaries, and vacancy-interstitial defects. The structural and thermodynamic properties of inherent structures are relatively insensitive to the temperature from which they are sampled, signifying that the energy landscape is relatively flat and devoid of deep wells. Using the nudged-elastic-band algorithm, we construct paths from ground-state configurations to inherent structures and identify the transition points between them. In addition, we use point patterns generated from a random sequential addition (RSA) of hard disks, which are nearly stealthy, and examine the particle rearrangements necessary to make the configurations absolutely stealthy. We introduce a configurational proximity metric to show that only small local, but collective, particle rearrangements are needed to drive initial RSA configurations to stealthy disordered ground states. These results lead to a more complete understanding of the unusual behaviors exhibited by the family of "collective-coordinate" potentials to which the $k$-space overlap potential belongs., Comment: 36 pages, 16 figures

We introduce a one-parameter family, $0 \leq H \leq 1$, of pair potential functions with a single relative energy minimum that stabilize a range of vacancy-riddled crystals as ground states. The "quintic potential" is a short-ranged, nonnegative pair potential with a single local minimum of height $H$ at unit distance and vanishes cubically at a distance of $\rt$. We have developed this potential to produce ground states with the symmetry of the triangular lattice while favoring the presence of vacancies. After an exhaustive search using various optimization and simulation methods, we believe that we have determined the ground states for all pressures, densities, and $0 \leq H \leq 1$. For specific areas below $3\rt/2$, the ground states of the "quintic potential" include high-density and low-density triangular lattices, kagom\'{e} and honeycomb crystals, and stripes. We find that these ground states are mechanically stable but are difficult to self-assemble in computer simulations without defects. For specific areas above $3\rt/2$, these systems have a ground-state phase diagram that corresponds to hard disks with radius $\rt$. For the special case of H=0, a broad range of ground states is available. Analysis of this case suggests that among many ground states, a high-density triangular lattice, low-density triangular lattice, and striped phases have the highest entropy for certain densities. The simplicity of this potential makes it an attractive candidate for experimental realization with application to the development of novel colloidal crystals or photonic materials., Comment: 25 pages, 11 figures

For any real sequence {c(n)} tending to infinity as n tends to infinity, this constructs a function f which is continuous and integrable, and such that for every nonzero x, limsup c(n) f(n x) is infinite., Comment: 8 pages, LaTeX. Corrected 2010Dec15. Minor corrections made 2011Jan18, especially E contains instead of E equals in the 6th line above section 2; some comments added

The phase behavior of hard superballs is examined using molecular dynamics within a deformable periodic simulation box. A superball's interior is defined by the inequality $|x|^{2q} + |y|^{2q} + |z|^{2q} \leq 1$, which provides a versatile family of convex particles ($q \geq 0.5$) with cube-like and octahedron-like shapes as well as concave particles ($q < 0.5$) with octahedron-like shapes. Here, we consider the convex case with a deformation parameter q between the sphere point (q = 1) and the cube (q = 1). We find that the asphericity plays a significant role in the extent of cubatic ordering of both the liquid and crystal phases. Calculation of the first few virial coefficients shows that superballs that are visually similar to cubes can have low-density equations of state closer to spheres than to cubes. Dense liquids of superballs display cubatic orientational order that extends over several particle lengths only for large q. Along the ordered, high-density equation of state, superballs with 1 < q < 3 exhibit clear evidence of a phase transition from a crystal state to a state with reduced long-ranged orientational order upon the reduction of density. For $q \geq 3$, long-ranged orientational order persists until the melting transition. The width of coexistence region between the liquid and ordered, high-density phase decreases with q up to q = 4.0. The structures of the high-density phases are examined using certain order parameters, distribution functions, and orientational correlation functions. We also find that a fixed simulation cell induces artificial phase transitions that are out of equilibrium. Current fabrication techniques allow for the synthesis of colloidal superballs, and thus the phase behavior of such systems can be investigated experimentally., Comment: 33 pages, 14 figures

We show that classical many-particle systems interacting with certain soft pair interactions in two dimensions exhibit novel low-temperature behaviors. Ground states span from disordered to crystalline. At some densities, a large fraction of normal-mode frequencies vanish. Lattice ground-state configurations have more vanishing frequencies than disordered ground states at the same density and exhibit vanishing shear moduli. For the melting transition from a crystal, the thermal expansion coefficient is negative. These unusual results are attributed to the topography of the energy landscape., Comment: 10 pages, 4 figures. This paper has been accepted for publication in Physical Review Letters

Using a collective coordinate numerical optimization procedure, we construct ground-state configurations of interacting particle systems in various space dimensions so that the scattering of radiation exactly matches a prescribed pattern for a set of wave vectors. We show that the constructed ground states are, counterintuitively, disordered (i.e., possess no long-range order) in the infinite-volume limit. We focus on three classes of configurations with unique radiation scattering characteristics: (i)``stealth'' materials, which are transparent to incident radiation at certain wavelengths; (ii)``super-ideal'' gases, which scatter radiation identically to that of an ensemble of ideal gas configurations for a selected set of wave vectors; and (iii)``equi-luminous'' materials, which scatter radiation equally intensely for a selected set of wave vectors. We find that ground-state configurations have an increased tendency to contain clusters of particles as one increases the prescribed luminosity. Limitations and consequences of this procedure are detailed., Comment: 44 pages, 16 figures, revtek4