Catalog footprint

What is connected

117works
42topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

117 published item(s)

preprint2026arXiv

Geometry over Density: Few-Shot Cross-Domain OOD Detection

Out-of-distribution (OOD) detection identifies test samples that fall outside a model's training distribution, a capability critical for safe deployment in high-stakes applications. Standard OOD detectors are trained on a specific in-distribution (ID) dataset and detect deviations from that single domain. In contrast, we study few-shot cross-domain OOD detection: given a \emph{single} pre-trained model, can we perform OOD detection on \emph{arbitrary} new ID-OOD task pairs using only a handful of ID samples at inference time, with no additional training? We propose \textbf{UFCOD}, a unified framework that achieves this goal through information-geometric analysis of diffusion trajectories. Our key insight is that diffusion noise predictions are score functions (gradients of log-density), and we extract two energy features: \emph{Path Energy} (integrated score magnitude) and \emph{Dynamics Energy} (score smoothness), that form a discrete Sobolev norm capturing how samples interact with the learned diffusion process. The central contribution is a \textbf{train-once, deploy-anywhere} paradigm: a diffusion model trained on a single dataset (e.g., CelebA) serves as a universal feature extractor for OOD detection across semantically unrelated domains (e.g., CIFAR-10, SVHN, Textures). At deployment, each new task requires only $\sim$100 unlabeled ID samples for inference: no retraining, no fine-tuning, no task-specific adaptation. Using 100 ID samples per task, UFCOD achieves 93.7\% average AUROC across 12 cross-domain benchmarks, competitive with methods trained on 50k--163k samples, demonstrating $\sim$500$\times$ improvement in sample efficiency. See our code in https://github.com/lili0415/UFCOD.

preprint2024arXiv

Exploration of faint X-ray and radio sources in the massive globular cluster M14: A UV-bright counterpart to Nova Ophiuchus 1938

Using a 12 ks archival Chandra X-ray Observatory ACIS-S observation on the massive globular cluster (GC) M14, we detect a total of 7 faint X-ray sources within its half-light radius at a 0.5-7 keV depth of $2.5\times 10^{31}\,\mathrm{erg~s^{-1}}$. We cross-match the X-ray source positions with a catalogue of the Very Large Array radio point sources and a Hubble Space Telescope (HST) UV/optical/near-IR photometry catalogue, revealing radio counterparts to 2 and HST counterparts to 6 of the X-ray sources. In addition, we also identify a radio source with the recently discovered millisecond pulsar PSR 1737-0314A. The brightest X-ray source, CX1, appears to be consistent with the nominal position of the classic nova Ophiuchi 1938 (Oph 1938), and both Oph 1938 and CX1 are consistent with a UV-bright variable HST counterpart, which we argue to be the source of the nova eruption in 1938. This makes Oph 1938 the second classic nova recovered in a Galactic GC since Nova T Scorpii in M80. CX2 is consistent with the steep-spectrum radio source VLA8, which unambiguously matches a faint blue source; the steepness of VLA8 is suggestive of a pulsar nature, possibly a transitional millisecond pulsar with a late K dwarf companion, though an active galactic nucleus (AGN) cannot be ruled out. The other counterparts to the X-ray sources are all suggestive of chromospherically active binaries or background AGNs, so their nature requires further membership information.

preprint2024arXiv

Inverse source problem of the biharmonic equation from multi-frequency phaseless data

This work deals with an inverse source problem for the biharmonic wave equation. A two-stage numerical method is proposed to identify the unknown source from the multi-frequency phaseless data. In the first stage, we introduce some artificially auxiliary point sources to the inverse source system and establish a phase retrieval formula. Theoretically, we point out that the phase can be uniquely determined and estimate the stability of this phase retrieval approach. Once the phase information is retrieved, the Fourier method is adopted to reconstruct the source function from the phased multi-frequency data. The proposed method is easy-to-implement and there is no forward solver involved in the reconstruction. Numerical experiments are conducted to verify the performance of the proposed method.

preprint2023arXiv

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more.

preprint2023arXiv

The average degree of edge chromatic critical graphs with maximum degree seven

In this paper, by developing several new adjacency lemmas about a path on $4$ or $5$ vertices, we show that the average degree of 7-critical graphs is at least 6. It implies Vizing's planar graph conjecture for planar graphs with maximum degree $7$ and its extension to graphs embeddable in a surface with nonnegative Euler characteristic due to Sanders and Zhao (J. Combin. Theory Ser. B 83 (2001) 201-212 and J. Combin. Theory Ser. B 87 (2003) 254-263) and Zhang (Graphs and Combinatorics 16 (2000) 467-495).

preprint2022arXiv

A Deep Reinforcement Learning Framework for Rapid Diagnosis of Whole Slide Pathological Images

The deep neural network is a research hotspot for histopathological image analysis, which can improve the efficiency and accuracy of diagnosis for pathologists or be used for disease screening. The whole slide pathological image can reach one gigapixel and contains abundant tissue feature information, which needs to be divided into a lot of patches in the training and inference stages. This will lead to a long convergence time and large memory consumption. Furthermore, well-annotated data sets are also in short supply in the field of digital pathology. Inspired by the pathologist's clinical diagnosis process, we propose a weakly supervised deep reinforcement learning framework, which can greatly reduce the time required for network inference. We use neural network to construct the search model and decision model of reinforcement learning agent respectively. The search model predicts the next action through the image features of different magnifications in the current field of view, and the decision model is used to return the predicted probability of the current field of view image. In addition, an expert-guided model is constructed by multi-instance learning, which not only provides rewards for search model, but also guides decision model learning by the knowledge distillation method. Experimental results show that our proposed method can achieve fast inference and accurate prediction of whole slide images without any pixel-level annotations.

preprint2022arXiv

All-sky search for gravitational wave emission from scalar boson clouds around spinning black holes in LIGO O3 data

This paper describes the first all-sky search for long-duration, quasi-monochromatic gravitational-wave signals emitted by ultralight scalar boson clouds around spinning black holes using data from the third observing run of Advanced LIGO. We analyze the frequency range from 20~Hz to 610~Hz, over a small frequency derivative range around zero, and use multiple frequency resolutions to be robust towards possible signal frequency wanderings. Outliers from this search are followed up using two different methods, one more suitable for nearly monochromatic signals, and the other more robust towards frequency fluctuations. We do not find any evidence for such signals and set upper limits on the signal strain amplitude, the most stringent being $\approx10^{-25}$ at around 130~Hz. We interpret these upper limits as both an "exclusion region" in the boson mass/black hole mass plane and the maximum detectable distance for a given boson mass, based on an assumption of the age of the black hole/boson cloud system.

preprint2022arXiv

BA-Net: Bridge Attention for Deep Convolutional Neural Networks

In recent years, channel attention mechanism has been widely investigated due to its great potential in improving the performance of deep convolutional neural networks (CNNs) in many vision tasks. However, in most of the existing methods, only the output of the adjacent convolution layer is fed into the attention layer for calculating the channel weights. Information from other convolution layers has been ignored. With these observations, a simple strategy, named Bridge Attention Net (BA-Net), is proposed in this paper for better performance with channel attention mechanisms. The core idea of this design is to bridge the outputs of the previous convolution layers through skip connections for channel weights generation. Based on our experiment and theory analysis, we find that features from previous layers also contribute to the weights significantly. The Comprehensive evaluation demonstrates that the proposed approach achieves state-of-the-art(SOTA) performance compared with the existing methods in accuracy and speed. which shows that Bridge Attention provides a new perspective on the design of neural network architectures with great potential in improving performance. The code is available at https://github.com/zhaoy376/Bridge-Attention.

preprint2022arXiv

Combining Individual and Joint Networking Behavior for Intelligent IoT Analytics

The IoT vision of a trillion connected devices over the next decade requires reliable end-to-end connectivity and automated device management platforms. While we have seen successful efforts for maintaining small IoT testbeds, there are multiple challenges for the efficient management of large-scale device deployments. With Industrial IoT, incorporating millions of devices, traditional management methods do not scale well. In this work, we address these challenges by designing a set of novel machine learning techniques, which form a foundation of a new tool, it IoTelligent, for IoT device management, using traffic characteristics obtained at the network level. The design of our tool is driven by the analysis of 1-year long networking data, collected from 350 companies with IoT deployments. The exploratory analysis of this data reveals that IoT environments follow the famous Pareto principle, such as: (i) 10% of the companies in the dataset contribute to 90% of the entire traffic; (ii) 7% of all the companies in the set own 90% of all the devices. We designed and evaluated CNN, LSTM, and Convolutional LSTM models for demand forecasting, with a conclusion of the Convolutional LSTM model being the best. However, maintaining and updating individual company models is expensive. In this work, we design a novel, scalable approach, where a general demand forecasting model is built using the combined data of all the companies with a normalization factor. Moreover, we introduce a novel technique for device management, based on autoencoders. They automatically extract relevant device features to identify device groups with similar behavior to flag anomalous devices.

preprint2022arXiv

Constraints from LIGO O3 data on gravitational-wave emission due to r-modes in the glitching pulsar PSR J0537-6910

We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537-6910 using data from the LIGO-Virgo Collaboration observing run O3. PSR J0537-6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role in the spin evolution of this pulsar. Theoretical models confirm this possibility and predict emission at a level that can be probed by ground-based detectors. In order to explore this scenario, we search for r-mode emission in the epochs between glitches by using a contemporaneous timing ephemeris obtained from NICER data. We do not detect any signals in the theoretically expected band of 86-97 Hz, and report upper limits on the amplitude of the gravitational waves. Our results improve on previous amplitude upper limits from r-modes in J0537-6910 by a factor of up to 3 and place stringent constraints on theoretical models for r-mode driven spin-down in PSR J0537-6910, especially for higher frequencies at which our results reach below the spin-down limit defined by energy conservation.

preprint2022arXiv

Deep Supervised Information Bottleneck Hashing for Cross-modal Retrieval based Computer-aided Diagnosis

Mapping X-ray images, radiology reports, and other medical data as binary codes in the common space, which can assist clinicians to retrieve pathology-related data from heterogeneous modalities (i.e., hashing-based cross-modal medical data retrieval), provides a new view to promot computeraided diagnosis. Nevertheless, there remains a barrier to boost medical retrieval accuracy: how to reveal the ambiguous semantics of medical data without the distraction of superfluous information. To circumvent this drawback, we propose Deep Supervised Information Bottleneck Hashing (DSIBH), which effectively strengthens the discriminability of hash codes. Specifically, the Deep Deterministic Information Bottleneck (Yu, Yu, and Principe 2021) for single modality is extended to the cross-modal scenario. Benefiting from this, the superfluous information is reduced, which facilitates the discriminability of hash codes. Experimental results demonstrate the superior accuracy of the proposed DSIBH compared with state-of-the-arts in cross-modal medical data retrieval tasks.

preprint2022arXiv

ECOD: Unsupervised Outlier Detection Using Empirical Cumulative Distribution Functions

Outlier detection refers to the identification of data points that deviate from a general data distribution. Existing unsupervised approaches often suffer from high computational cost, complex hyperparameter tuning, and limited interpretability, especially when working with large, high-dimensional datasets. To address these issues, we present a simple yet effective algorithm called ECOD (Empirical-Cumulative-distribution-based Outlier Detection), which is inspired by the fact that outliers are often the "rare events" that appear in the tails of a distribution. In a nutshell, ECOD first estimates the underlying distribution of the input data in a nonparametric fashion by computing the empirical cumulative distribution per dimension of the data. ECOD then uses these empirical distributions to estimate tail probabilities per dimension for each data point. Finally, ECOD computes an outlier score of each data point by aggregating estimated tail probabilities across dimensions. Our contributions are as follows: (1) we propose a novel outlier detection method called ECOD, which is both parameter-free and easy to interpret; (2) we perform extensive experiments on 30 benchmark datasets, where we find that ECOD outperforms 11 state-of-the-art baselines in terms of accuracy, efficiency, and scalability; and (3) we release an easy-to-use and scalable (with distributed support) Python implementation for accessibility and reproducibility.

preprint2022arXiv

Federated Learning with Non-IID Data

Federated learning enables resource-constrained edge compute devices, such as mobile phones and IoT devices, to learn a shared model for prediction, while keeping the training data local. This decentralized approach to train models provides privacy, security, regulatory and economic benefits. In this work, we focus on the statistical challenge of federated learning when local data is non-IID. We first show that the accuracy of federated learning reduces significantly, by up to 55% for neural networks trained for highly skewed non-IID data, where each client device trains only on a single class of data. We further show that this accuracy reduction can be explained by the weight divergence, which can be quantified by the earth mover's distance (EMD) between the distribution over classes on each device and the population distribution. As a solution, we propose a strategy to improve training on non-IID data by creating a small subset of data which is globally shared between all the edge devices. Experiments show that accuracy can be increased by 30% for the CIFAR-10 dataset with only 5% globally shared data.

preprint2022arXiv

First joint observation by the underground gravitational-wave detector, KAGRA, with GEO600

We report the results of the first joint observation of the KAGRA detector with GEO600. KAGRA is a cryogenic and underground gravitational-wave detector consisting of a laser interferometer with three-kilometer arms, and located in Kamioka, Gifu, Japan. GEO600 is a British--German laser interferometer with 600 m arms, and located near Hannover, Germany. GEO600 and KAGRA performed a joint observing run from April 7 to 20, 2020. We present the results of the joint analysis of the GEO--KAGRA data for transient gravitational-wave signals, including the coalescence of neutron-star binaries and generic unmodeled transients. We also perform dedicated searches for binary coalescence signals and generic transients associated with gamma-ray burst events observed during the joint run. No gravitational-wave events were identified. We evaluate the minimum detectable amplitude for various types of transient signals and the spacetime volume for which the network is sensitive to binary neutron-star coalescences. We also place lower limits on the distances to the gamma-ray bursts analysed based on the non-detection of an associated gravitational-wave signal for several signal models, including binary coalescences. These analyses demonstrate the feasibility and utility of KAGRA as a member of the global gravitational-wave detector network.

preprint2022arXiv

Gaussian Kernel Variance For an Adaptive Learning Method on Signals Over Graphs

This paper discusses a special kind of a simple yet possibly powerful algorithm, called single-kernel Gradraker (SKG), which is an adaptive learning method predicting unknown nodal values in a network using known nodal values and the network structure. We aim to find out how to configure the special kind of the model in applying the algorithm. To be more specific, we focus on SKG with a Gaussian kernel and specify how to find a suitable variance for the kernel. To do so, we introduce two variables with which we are able to set up requirements on the variance of the Gaussian kernel to achieve (near-) optimal performance and can better understand how SKG works. Our contribution is that we introduce two variables as analysis tools, illustrate how predictions will be affected under different Gaussian kernels, and provide an algorithm finding a suitable Gaussian kernel for SKG with knowledge about the training network. Simulation results on real datasets are provided.

preprint2022arXiv

GCN-MIF: Graph Convolutional Network with Multi-Information Fusion for Low-dose CT Denoising

Being low-level radiation exposure and less harmful to health, low-dose computed tomography (LDCT) has been widely adopted in the early screening of lung cancer and COVID-19. LDCT images inevitably suffer from the degradation problem caused by complex noises. It was reported that deep learning (DL)-based LDCT denoising methods using convolutional neural network (CNN) achieved impressive denoising performance. Although most existing DL-based methods (e.g., encoder-decoder framework) can implicitly utilize non-local and contextual information via downsampling operator and 3D CNN, the explicit multi-information (i.e., local, non-local, and contextual) integration may not be explored enough. To address this issue, we propose a novel graph convolutional network-based LDCT denoising model, namely GCN-MIF, to explicitly perform multi-information fusion for denoising purpose. Concretely, by constructing intra- and inter-slice graph, the graph convolutional network is introduced to leverage the non-local and contextual relationships among pixels. The traditional CNN is adopted for the extraction of local information. Finally, the proposed GCN-MIF model fuses all the extracted local, non-local, and contextual information. Extensive experiments show the effectiveness of our proposed GCN-MIF model by quantitative and visualized results. Furthermore, a double-blind reader study on a public clinical dataset is also performed to validate the usability of denoising results in terms of the structural fidelity, the noise suppression, and the overall score. Models and code are available at https://github.com/tonyckc/GCN-MIF_demo.

preprint2022arXiv

GWTC-2.1: Deep Extended Catalog of Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run

The second Gravitational-Wave Transient Catalog reported on 39 compact binary coalescences observed by the Advanced LIGO and Advanced Virgo detectors between 1 April 2019 15:00 UTC and 1 October 2019 15:00 UTC. We present GWTC-2.1, which reports on a deeper list of candidate events observed over the same period. We analyze the final version of the strain data over this period with improved calibration and better subtraction of excess noise, which has been publicly released. We employ three matched-filter search pipelines for candidate identification, and estimate the astrophysical probability for each candidate event. While GWTC-2 used a false alarm rate threshold of 2 per year, we include in GWTC-2.1, 1201 candidates that pass a false alarm rate threshold of 2 per day. We calculate the source properties of a subset of 44 high-significance candidates that have an astrophysical probability greater than 0.5. Of these candidates, 36 have been reported in GWTC-2. If the 8 additional high-significance candidates presented here are astrophysical, the mass range of events that are unambiguously identified as binary black holes (both objects $\geq 3M_\odot$) is increased compared to GWTC-2, with total masses from $\sim 14 M_\odot$ for GW190924_021846 to $\sim 182 M_\odot$ for GW190426_190642. The primary components of two new candidate events (GW190403_051519 and GW190426_190642) fall in the mass gap predicted by pair instability supernova theory. We also expand the population of binaries with significantly asymmetric mass ratios reported in GWTC-2 by an additional two events (the mass ratio is less than $0.65$ and $0.44$ at $90\%$ probability for GW190403_051519 and GW190917_114630 respectively), and find that 2 of the 8 new events have effective inspiral spins $χ_\mathrm{eff} > 0$ (at $90\%$ credibility), while no binary is consistent with $χ_\mathrm{eff} < 0$ at the same significance.

preprint2022arXiv

Improved random batch Ewald method in molecular dynamics simulations

The random batch Ewald (RBE) is an efficient and accurate method for molecular dynamics (MD) simulations of physical systems at the nano-/micro- scale. The method shows great potential to solve the computational bottleneck of long-range interactions, motivating a necessity to accelerating short-range components of the non-bonded interactions for a further speedup of MD simulations. In this work, we present an improved RBE method for the non-bonding interactions by introducing the random batch idea to constructing neighbor lists for the treatment of both the short-range part of the Ewald splitting and the Lennard-Jones potential. The efficiency of the novel neighbor list algorithm owes to the stochastic minibatch strategy which can significantly reduce the total number of neighbors. We obtan the error estimate and convergence by theoretical analysis and implement the improved RBE method in the LAMMPS package. Benchmark simulations are performed to demonstrate the accuracy and stability of the algorithm. Numerical tests on computer performance by conducting large-scaled MD simulations for systems including up to 0.1 billion water molecules, run on massive cluster with up to 50 thousand CPU cores, demonstrating the attractive features such as the high parallel scalability and memory-saving of the method in comparison to the existing methods.

preprint2022arXiv

Intrinsically Motivated Self-supervised Learning in Reinforcement Learning

In vision-based reinforcement learning (RL) tasks, it is prevalent to assign auxiliary tasks with a surrogate self-supervised loss so as to obtain more semantic representations and improve sample efficiency. However, abundant information in self-supervised auxiliary tasks has been disregarded, since the representation learning part and the decision-making part are separated. To sufficiently utilize information in auxiliary tasks, we present a simple yet effective idea to employ self-supervised loss as an intrinsic reward, called Intrinsically Motivated Self-Supervised learning in Reinforcement learning (IM-SSR). We formally show that the self-supervised loss can be decomposed as exploration for novel states and robustness improvement from nuisance elimination. IM-SSR can be effortlessly plugged into any reinforcement learning with self-supervised auxiliary objectives with nearly no additional cost. Combined with IM-SSR, the previous underlying algorithms achieve salient improvements on both sample efficiency and generalization in various vision-based robotics tasks from the DeepMind Control Suite, especially when the reward signal is sparse.

preprint2022arXiv

Learning Robust Representation through Graph Adversarial Contrastive Learning

Existing studies show that node representations generated by graph neural networks (GNNs) are vulnerable to adversarial attacks, such as unnoticeable perturbations of adjacent matrix and node features. Thus, it is requisite to learn robust representations in graph neural networks. To improve the robustness of graph representation learning, we propose a novel Graph Adversarial Contrastive Learning framework (GraphACL) by introducing adversarial augmentations into graph self-supervised learning. In this framework, we maximize the mutual information between local and global representations of a perturbed graph and its adversarial augmentations, where the adversarial graphs can be generated in either supervised or unsupervised approaches. Based on the Information Bottleneck Principle, we theoretically prove that our method could obtain a much tighter bound, thus improving the robustness of graph representation learning. Empirically, we evaluate several methods on a range of node classification benchmarks and the results demonstrate GraphACL could achieve comparable accuracy over previous supervised methods.

preprint2022arXiv

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets.

preprint2022arXiv

Probing Early Universe Supercooled Phase Transitions with Gravitational Wave Data

We investigate the reach of the LIGO/Virgo/KAGRA detectors in the search for signatures of first-order phase transitions in the early Universe. Utilising data from the first three observing runs, we derive constraints on the parameters of the underlying gravitational-wave background, focusing on transitions characterised by strong supercooling. As an application of our analysis, we determine bounds on the parameter space of two representative particle physics models. We also comment on the expected reach of third-generation detectors in probing supercooled phase transitions.

preprint2022arXiv

Producing and detecting long-lived particles at different experiments at the LHC

We propose a new strategy to look for long-lived particles (LLP) at the LHC. The LLPs are produced at one experiment, but its decay products are detected by a detector at another experiment. We use a confining Hidden Valley scenario as a benchmark. Through showering and hadronization, the multiplicity of hidden mesons can be large, and their decay products, dimuon as chosen in this study, are typically too soft to pass triggers in traditional LHC searches. We find the best acceptance is achieved if we produce LLPs at collision points at the LHCb and ALICE experiments, and use the muon chamber of ATLAS for detection. This new search is cost-efficient since it does not require a new detector to be built. Meanwhile, it can provide coverage of interesting parameter space, which is complementary to other proposed LLP searches.

preprint2022arXiv

Revisiting Skeleton-based Action Recognition

Human skeleton, as a compact representation of human action, has received increasing attention in recent years. Many skeleton-based action recognition methods adopt graph convolutional networks (GCN) to extract features on top of human skeletons. Despite the positive results shown in previous works, GCN-based methods are subject to limitations in robustness, interoperability, and scalability. In this work, we propose PoseC3D, a new approach to skeleton-based action recognition, which relies on a 3D heatmap stack instead of a graph sequence as the base representation of human skeletons. Compared to GCN-based methods, PoseC3D is more effective in learning spatiotemporal features, more robust against pose estimation noises, and generalizes better in cross-dataset settings. Also, PoseC3D can handle multiple-person scenarios without additional computation cost, and its features can be easily integrated with other modalities at early fusion stages, which provides a great design space to further boost the performance. On four challenging datasets, PoseC3D consistently obtains superior performance, when used alone on skeletons and in combination with the RGB modality.

preprint2022arXiv

Search for anisotropic gravitational-wave backgrounds using data from Advanced LIGO and Advanced Virgo's first three observing runs

We report results from searches for anisotropic stochastic gravitational-wave backgrounds using data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. For the first time, we include Virgo data in our analysis and run our search with a new efficient pipeline called {\tt PyStoch} on data folded over one sidereal day. We use gravitational-wave radiometry (broadband and narrow band) to produce sky maps of stochastic gravitational-wave backgrounds and to search for gravitational waves from point sources. A spherical harmonic decomposition method is employed to look for gravitational-wave emission from spatially-extended sources. Neither technique found evidence of gravitational-wave signals. Hence we derive 95\% confidence-level upper limit sky maps on the gravitational-wave energy flux from broadband point sources, ranging from $F_{α, Θ} < {\rm (0.013 - 7.6)} \times 10^{-8} {\rm erg \, cm^{-2} \, s^{-1} \, Hz^{-1}},$ and on the (normalized) gravitational-wave energy density spectrum from extended sources, ranging from $Ω_{α, Θ} < {\rm (0.57 - 9.3)} \times 10^{-9} \, {\rm sr^{-1}}$, depending on direction ($Θ$) and spectral index ($α$). These limits improve upon previous limits by factors of $2.9 - 3.5$. We also set 95\% confidence level upper limits on the frequency-dependent strain amplitudes of quasimonochromatic gravitational waves coming from three interesting targets, Scorpius X-1, SN 1987A and the Galactic Center, with best upper limits range from $h_0 < {\rm (1.7-2.1)} \times 10^{-25},$ a factor of $\geq 2.0$ improvement compared to previous stochastic radiometer searches.

preprint2022arXiv

Search for continuous gravitational wave emission from the Milky Way center in O3 LIGO--Virgo data

We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO--Virgo run in the detector frequency band $[10,2000]\rm~Hz$ have been used. No significant detection was found and 95$\%$ confidence level upper limits on the signal strain amplitude were computed, over the full search band, with the deepest limit of about $7.6\times 10^{-26}$ at $\simeq 142\rm~Hz$. These results are significantly more constraining than those reported in previous searches. We use these limits to put constraints on the fiducial neutron star ellipticity and r-mode amplitude. These limits can be also translated into constraints in the black hole mass -- boson mass plane for a hypothetical population of boson clouds around spinning black holes located in the GC.

preprint2022arXiv

Search for continuous gravitational waves from 20 accreting millisecond X-ray pulsars in O3 LIGO data

Results are presented of searches for continuous gravitational waves from 20 accreting millisecond X-ray pulsars with accurately measured spin frequencies and orbital parameters, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. The search algorithm uses a hidden Markov model, where the transition probabilities allow the frequency to wander according to an unbiased random walk, while the $\mathcal{J}$-statistic maximum-likelihood matched filter tracks the binary orbital phase. Three narrow sub-bands are searched for each target, centered on harmonics of the measured spin frequency. The search yields 16 candidates, consistent with a false alarm probability of 30% per sub-band and target searched. These candidates, along with one candidate from an additional target-of-opportunity search done for SAX J1808.4$-$3658, which was in outburst during one month of the observing run, cannot be confidently associated with a known noise source. Additional follow-up does not provide convincing evidence that any are a true astrophysical signal. When all candidates are assumed non-astrophysical, upper limits are set on the maximum wave strain detectable at 95% confidence, $h_0^{95\%}$. The strictest constraint is $h_0^{95\%} = 4.7\times 10^{-26}$ from IGR J17062$-$6143. Constraints on the detectable wave strain from each target lead to constraints on neutron star ellipticity and $r$-mode amplitude, the strictest of which are $ε^{95\%} = 3.1\times 10^{-7}$ and $α^{95\%} = 1.8\times 10^{-5}$ respectively. This analysis is the most comprehensive and sensitive search of continuous gravitational waves from accreting millisecond X-ray pulsars to date.

preprint2022arXiv

Search of the Early O3 LIGO Data for Continuous Gravitational Waves from the Cassiopeia A and Vela Jr. Supernova Remnants

We present directed searches for continuous gravitational waves from the neutron stars in the Cassiopeia A (Cas A) and Vela Jr. supernova remnants. We carry out the searches in the LIGO data from the first six months of the third Advanced LIGO and Virgo observing run, using the Weave semi-coherent method, which sums matched-filter detection-statistic values over many time segments spanning the observation period. No gravitational wave signal is detected in the search band of 20--976 Hz for assumed source ages greater than 300 years for Cas A and greater than 700 years for Vela Jr. Estimates from simulated continuous wave signals indicate we achieve the most sensitive results to date across the explored parameter space volume, probing to strain magnitudes as low as ~$6.3\times10^{-26}$ for Cas A and ~$5.6\times10^{-26}$ for Vela Jr. at frequencies near 166 Hz at 95% efficiency.

preprint2022arXiv

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.

preprint2022arXiv

Stability for the multi-dimensional Borg--Levinson theorem of the biharmonic operator

In this paper, we prove a conditional Hölder stability estimate for the inverse spectral problem of the biharmonic operator. The proof employs the resolvent estimate and a Weyl-type law for the biharmonic operator which were obtained by the authors in \cite{LYZ}. This work extends nontrivially the result in \cite{stefanov} from the second order Schrödinger operator to the fourth order biharmonic operator.

preprint2022arXiv

Stochastic Gravitational Wave Background from PBH-ABH Mergers

The measurement of gravitational waves produced by binary black-hole mergers at the Advanced LIGO has encouraged extensive studies on the stochastic gravitational wave background. Recent studies have focused on gravitational wave sources made of the same species, such as mergers from binary primordial black holes or those from binary astrophysical black holes. In this paper, we study a new possibility -- the stochastic gravitational wave background produced by mergers of one primordial black hole and one astrophysical black hole. Such systems are necessarily present if primordial black holes exist. We study the isotropic gravitational wave background produced through the history of the Universe. We find it is very challenging to detect such a signal. We also demonstrate that it is improper to treat the gravitational waves produced by such binaries in the Milky Way as a directional stochastic background, due to a very low binary formation rate.

preprint2022arXiv

Stringent axion constraints with Event Horizon Telescope polarimetric measurements of M87$^\star$

The hitherto unprecedented angular resolution of the Event Horizon Telescope (EHT) has created exciting opportunities in the search for new physics. Recently, the linear polarization of radiation emitted near the supermassive black hole M87$^\star$ was measured on four separate days, precisely enabling tests of the existence of a dense axion cloud produced by a spinning black hole. The presence of an axion cloud leads to a frequency-independent oscillation in the electric vector position angle (EVPA) of this linear polarization. For a nearly face-on M87$^\star$, this oscillation in the EVPA appears as a propagating wave along the photon ring. In this paper, we leverage the azimuthal distribution of EVPA measured by the EHT to study the axion-photon coupling. We propose a novel differential analysis procedure to reduce the astrophysical background, and derive stringent constraints on the existence of axions in the previously unexplored mass window $\sim (10^{-21}-10^{-20})$~eV.

preprint2022arXiv

Symmetry breaking and anomalous conductivity in a double moiré superlattice

A double moiré superlattice can be realized by stacking three layers of atomically thin two-dimensional materials with designer interlayer twisting or lattice mismatches. In this novel structure, atomic reconstruction of constituent layers could introduce significant modifications to the lattice symmetry and electronic structure at small twist angles. Here, we employ conductive atomic force microscopy (cAFM) to investigate symmetry breaking and local electrical properties in twisted trilayer graphene. We observe clear double moiré superlattices with two distinct moire periods all over the sample. At neighboring domains of the large moiré, the current exhibit either two- or six-fold rotational symmetry, indicating delicate symmetry breaking beyond the rigid model. Moreover, an anomalous current appears at the 'A-A' stacking site of the larger moiré, contradictory to previous observations on twisted bilayer graphene. Both behaviors can be understood by atomic reconstruction, and we also show that the cAFM signal of twisted graphene samples is dominated by the tip-graphene contact resistance that maps the local work function of twisted graphene and the metallic tip qualitatively. Our results unveil cAFM is an effective probe for visualizing atomic reconstruction and symmetry breaking in novel moiré superlattices, which could provide new insights for exploring and manipulating more exotic quantum states based on twisted van der Waals heterostructures.

preprint2022arXiv

The MAVERIC Survey: The first radio and X-ray limits on the detached black holes in NGC 3201

The Galactic globular cluster NGC 3201 is the first Galactic globular cluster observed to host dynamically-confirmed stellar-mass black holes, containing two confirmed and one candidate black hole. This result indicates that globular clusters can retain black holes, which has important implications for globular cluster evolution. NGC 3201 has been observed as part of the MAVERIC survey of Galactic globular clusters. We use these data to confirm that there is no radio or X-ray detection of the three black holes, and present the first radio and X-ray limits on these sources. These limits indicate that any accretion present is at an extremely low rate and may be extremely inefficient. In particular, for the system ACS ID #21859, by assuming the system is tidally locked and any accretion is through the capture of the companion's winds, we constrain the radiative efficiency of any accretion to $\leq1.5\times10^{-5}$. We also combine the radio and X-ray source catalogues from the MAVERIC survey with the existing MUSE spectroscopic surveys and the HUGS catalogue of NGC 3201 to provide a catalogue of 42 multiwavelength sources in this cluster. We identify a new red straggler source with X-ray emission, and investigate the multiwavelength properties of the sub-subgiant population in the cluster.

preprint2022arXiv

Two-Stream Graph Convolutional Network for Intra-oral Scanner Image Segmentation

Precise segmentation of teeth from intra-oral scanner images is an essential task in computer-aided orthodontic surgical planning. The state-of-the-art deep learning-based methods often simply concatenate the raw geometric attributes (i.e., coordinates and normal vectors) of mesh cells to train a single-stream network for automatic intra-oral scanner image segmentation. However, since different raw attributes reveal completely different geometric information, the naive concatenation of different raw attributes at the (low-level) input stage may bring unnecessary confusion in describing and differentiating between mesh cells, thus hampering the learning of high-level geometric representations for the segmentation task. To address this issue, we design a two-stream graph convolutional network (i.e., TSGCN), which can effectively handle inter-view confusion between different raw attributes to more effectively fuse their complementary information and learn discriminative multi-view geometric representations. Specifically, our TSGCN adopts two input-specific graph-learning streams to extract complementary high-level geometric representations from coordinates and normal vectors, respectively. Then, these single-view representations are further fused by a self-attention module to adaptively balance the contributions of different views in learning more discriminative multi-view representations for accurate and fully automatic tooth segmentation. We have evaluated our TSGCN on a real-patient dataset of dental (mesh) models acquired by 3D intraoral scanners. Experimental results show that our TSGCN significantly outperforms state-of-the-art methods in 3D tooth (surface) segmentation. Github: https://github.com/ZhangLingMing1/TSGCNet.

preprint2021arXiv

All-sky search for long-duration gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

After the detection of gravitational waves from compact binary coalescences, the search for transient gravitational-wave signals with less well-defined waveforms for which matched filtering is not well-suited is one of the frontiers for gravitational-wave astronomy. Broadly classified into "short" $ \lesssim 1~$\,s and "long" $ \gtrsim 1~$\,s duration signals, these signals are expected from a variety of astrophysical processes, including non-axisymmetric deformations in magnetars or eccentric binary black hole coalescences. In this work, we present a search for long-duration gravitational-wave transients from Advanced LIGO and Advanced Virgo's third observing run from April 2019 to March 2020. For this search, we use minimal assumptions for the sky location, event time, waveform morphology, and duration of the source. The search covers the range of $2~\text{--}~ 500$~s in duration and a frequency band of $24 - 2048$ Hz. We find no significant triggers within this parameter space; we report sensitivity limits on the signal strength of gravitational waves characterized by the root-sum-square amplitude $h_{\mathrm{rss}}$ as a function of waveform morphology. These $h_{\mathrm{rss}}$ limits improve upon the results from the second observing run by an average factor of 1.8.

preprint2021arXiv

All-sky search for short gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

This paper presents the results of a search for generic short-duration gravitational-wave transients in data from the third observing run of Advanced LIGO and Advanced Virgo. Transients with durations of milliseconds to a few seconds in the 24--4096 Hz frequency band are targeted by the search, with no assumptions made regarding the incoming signal direction, polarization or morphology. Gravitational waves from compact binary coalescences that have been identified by other targeted analyses are detected, but no statistically significant evidence for other gravitational wave bursts is found. Sensitivities to a variety of signals are presented. These include updated upper limits on the source rate-density as a function of the characteristic frequency of the signal, which are roughly an order of magnitude better than previous upper limits. This search is sensitive to sources radiating as little as $\sim$10$^{-10} M_{\odot} c^2$ in gravitational waves at $\sim$70 Hz from a distance of 10~kpc, with 50\% detection efficiency at a false alarm rate of one per century. The sensitivity of this search to two plausible astrophysical sources is estimated: neutron star f-modes, which may be excited by pulsar glitches, as well as selected core-collapse supernova models.

preprint2021arXiv

All-sky, all-frequency directional search for persistent gravitational-waves from Advanced LIGO's and Advanced Virgo's first three observing runs

We present the first results from an all-sky all-frequency (ASAF) search for an anisotropic stochastic gravitational-wave background using the data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. Upper limit maps on broadband anisotropies of a persistent stochastic background were published for all observing runs of the LIGO-Virgo detectors. However, a broadband analysis is likely to miss narrowband signals as the signal-to-noise ratio of a narrowband signal can be significantly reduced when combined with detector output from other frequencies. Data folding and the computationally efficient analysis pipeline, {\tt PyStoch}, enable us to perform the radiometer map-making at every frequency bin. We perform the search at 3072 {\tt{HEALPix}} equal area pixels uniformly tiling the sky and in every frequency bin of width $1/32$~Hz in the range $20-1726$~Hz, except for bins that are likely to contain instrumental artefacts and hence are notched. We do not find any statistically significant evidence for the existence of narrowband gravitational-wave signals in the analyzed frequency bins. Therefore, we place $95\%$ confidence upper limits on the gravitational-wave strain for each pixel-frequency pair, the limits are in the range $(0.030 - 9.6) \times10^{-24}$. In addition, we outline a method to identify candidate pixel-frequency pairs that could be followed up by a more sensitive (and potentially computationally expensive) search, e.g., a matched-filtering-based analysis, to look for fainter nearly monochromatic coherent signals. The ASAF analysis is inherently independent of models describing any spectral or spatial distribution of power. We demonstrate that the ASAF results can be appropriately combined over frequencies and sky directions to successfully recover the broadband directional and isotropic results.

preprint2021arXiv

Chandra and HST Studies of Six Millisecond Pulsars in the Globular Cluster M13

We analyse 55 ks of Chandra X-ray observations of the Galactic globular cluster M13. Using the latest radio timing positions of six known millisecond pulsars (MSPs) in M13 from Wang et al. (2020), we detect confident X-ray counterparts to five of the six MSPs at X-ray luminosities of $L_X$(0.3-8 keV)$\sim 3 \times 10^{30} - 10^{31}~{\rm erg~s^{-1}}$, including the newly discovered PSR J1641+3627F. There are limited X-ray counts at the position of PSR J1641+3627A, for which we obtain an upper limit $L_X<1.3 \times 10^{30}~{\rm erg~s^{-1}}$. We analyse X-ray spectra of all six MSPs, which are well-described by either a single blackbody or a single power-law model. We also incorporate optical/UV imaging observations from the Hubble Space Telescope (HST) and find optical counterparts to PSR J1641+3627D and J1641+3627F. Our colour-magnitude diagrams indicate the latter contains a white dwarf, consistent with the properties suggested by radio timing observations. The counterpart to J1641+3627D is only visible in the V band; however, we argue that the companion to J1641+3627D is also a white dwarf, since we see a blackbody-like X-ray spectrum, while MSPs with nondegenerate companions generally show non-thermal X-rays from shocks between the pulsar and companion winds. Our work increases the sample of known X-ray and optical counterparts of MSPs in globular clusters.

preprint2021arXiv

Inorganic component imaging of aggregate glue droplets on spider orb webs by TOF-SIMS

In this review, we discuss the use of time-of-flight secondary-ion mass spectrometry (TOF-SIMS) technology for analyzing the viscous glue (is called aggregate glue droplets) of spider orb webs and examine the results obtained. Element distribution images of the aggregate glue droplets were observed by TOF-SIMS. A uniform element distribution is seen for suspended pristine aggregate glue droplets, and a differential spreading of aggregate glue components is seen for attached aggregate glue droplets. We also observed TOF-SIMS images of water in aggregate glue droplets, where water was observed to be consistent with the distribution of oozing salt. We also found that the alkali metal in the aggregate glue droplets showed similar characteristics by feeding cesium carbonate to spiders.

preprint2021arXiv

PyHealth: A Python Library for Health Predictive Models

Despite the explosion of interest in healthcare AI research, the reproducibility and benchmarking of those research works are often limited due to the lack of standard benchmark datasets and diverse evaluation metrics. To address this reproducibility challenge, we develop PyHealth, an open-source Python toolbox for developing various predictive models on healthcare data. PyHealth consists of data preprocessing module, predictive modeling module, and evaluation module. The target users of PyHealth are both computer science researchers and healthcare data scientists. With PyHealth, they can conduct complex machine learning pipelines on healthcare datasets with fewer than ten lines of code. The data preprocessing module enables the transformation of complex healthcare datasets such as longitudinal electronic health records, medical images, continuous signals (e.g., electrocardiogram), and clinical notes into machine learning friendly formats. The predictive modeling module provides more than 30 machine learning models, including established ensemble trees and deep neural network-based approaches, via a unified but extendable API designed for both researchers and practitioners. The evaluation module provides various evaluation strategies (e.g., cross-validation and train-validation-test split) and predictive model metrics. With robustness and scalability in mind, best practices such as unit testing, continuous integration, code coverage, and interactive examples are introduced in the library's development. PyHealth can be installed through the Python Package Index (PyPI) or https://github.com/yzhao062/PyHealth .

preprint2021arXiv

RNN-Test: Towards Adversarial Testing for Recurrent Neural Network Systems

While massive efforts have been investigated in adversarial testing of convolutional neural networks (CNN), testing for recurrent neural networks (RNN) is still limited and leaves threats for vast sequential application domains. In this paper, we propose an adversarial testing framework RNN-Test for RNN systems, focusing on the main sequential domains, not only classification tasks. First, we design a novel search methodology customized for RNN models by maximizing the inconsistency of RNN states to produce adversarial inputs. Next, we introduce two state-based coverage metrics according to the distinctive structure of RNNs to explore more inference logics. Finally, RNN-Test solves the joint optimization problem to maximize state inconsistency and state coverage, and crafts adversarial inputs for various tasks of different kinds of inputs. For evaluations, we apply RNN-Test on three sequential models of common RNN structures. On the tested models, the RNN-Test approach is demonstrated to be competitive in generating adversarial inputs, outperforming FGSM-based and DLFuzz-based methods to reduce the model performance more sharply with 2.78% to 32.5% higher success (or generation) rate. RNN-Test could also achieve 52.65% to 66.45% higher adversary rate on MNIST-LSTM model than relevant work testRNN. Compared with the neuron coverage, the proposed state coverage metrics as guidance excel with 4.17% to 97.22% higher success (or generation) rate.

preprint2021arXiv

Search for Gravitational Waves Associated with Gamma-Ray Bursts Detected by Fermi and Swift During the LIGO-Virgo Run O3b

We search for gravitational-wave signals associated with gamma-ray bursts detected by the Fermi and Swift satellites during the second half of the third observing run of Advanced LIGO and Advanced Virgo (1 November 2019 15:00 UTC-27 March 2020 17:00 UTC).We conduct two independent searches: a generic gravitational-wave transients search to analyze 86 gamma-ray bursts and an analysis to target binary mergers with at least one neutron star as short gamma-ray burst progenitors for 17 events. We find no significant evidence for gravitational-wave signals associated with any of these gamma-ray bursts. A weighted binomial test of the combined results finds no evidence for sub-threshold gravitational wave signals associated with this GRB ensemble either. We use several source types and signal morphologies during the searches, resulting in lower bounds on the estimated distance to each gamma-ray burst. Finally, we constrain the population of low luminosity short gamma-ray bursts using results from the first to the third observing runs of Advanced LIGO and Advanced Virgo. The resulting population is in accordance with the local binary neutron star merger rate.

preprint2021arXiv

Search for subsolar-mass binaries in the first half of Advanced LIGO and Virgo's third observing run

We report on a search for compact binary coalescences where at least one binary component has a mass between 0.2 $M_\odot$ and 1.0 $M_\odot$ in Advanced LIGO and Advanced Virgo data collected between 1 April 2019 1500 UTC and 1 October 2019 1500 UTC. We extend previous analyses in two main ways: we include data from the Virgo detector and we allow for more unequal mass systems, with mass ratio $q \geq 0.1$. We do not report any gravitational-wave candidates. The most significant trigger has a false alarm rate of 0.14 $\mathrm{yr}^{-1}$. This implies an upper limit on the merger rate of subsolar binaries in the range $[220-24200] \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$, depending on the chirp mass of the binary. We use this upper limit to derive astrophysical constraints on two phenomenological models that could produce subsolar-mass compact objects. One is an isotropic distribution of equal-mass primordial black holes. Using this model, we find that the fraction of dark matter in primordial black holes is $f_\mathrm{PBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 6\%$. The other is a dissipative dark matter model, in which fermionic dark matter can collapse and form black holes. The upper limit on the fraction of dark matter black holes depends on the minimum mass of the black holes that can be formed: the most constraining result is obtained at $M_\mathrm{min}=1 M_\odot$, where $f_\mathrm{DBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 0.003\%$. These are the tightest limits on spinning subsolar-mass binaries to date.

preprint2021arXiv

Stability for an inverse source problem of the biharmonic operator

In this paper, we study for the first time the stability of the inverse source problem for the biharmonic operator with a compactly supported potential in $\mathbb R^3$. Firstly, to connect the boundary data with the unknown source, we shall consider an eigenvalue problem for the bi-Schr$\ddot{\rm o}$dinger operator $Δ^2 + V(x)$ on a ball which contains the support of the potential $V$. We prove a Weyl-type law for the upper bounds of spherical normal derivatives of both the eigenfunctions $ϕ$ and their Laplacian $Δϕ$ corresponding to the bi-Schr$\ddot{\rm o}$dinger operator. This type of upper bounds was proved by Hassell and Tao for the Schr$\ddot{\rm o}$dinger operator. Secondly, we investigate the meromorphic continuation of the resolvent of the bi-Schr$\ddot{\rm o}$dinger operator and prove the existence of a resonance-free region and an estimate of $L^2_{\rm comp} - L^2_{\rm loc}$ type for the resolvent. As an application, we prove a bound of the analytic continuation of the data from the given data to the higher frequency data. Finally, we derive the stability estimate which consists of the Lipschitz type data discrepancy and the high frequency tail of the source function, where the latter decreases as the upper bound of the frequency increases.

preprint2021arXiv

Stability for an inverse source problem of the damped biharmonic plate equation

This paper is concerned with the stability of the inverse source problem for the damped biharmonic plate equation in three dimensions. The stability estimate consists of the Lipschitz type data discrepancy and the high frequency tail of the source function, where the latter decreases as the upper bound of the frequency increases. The stability also shows exponential dependence on the constant damping coefficient. The analysis employs Carleman estimates and time decay estimates for the damped plate wave equation to obtain an exact observability bound and depends on the study of the resonance-free region and an upper bound of the resolvent of the biharmonic operator with respect to the complex wavenumber.

preprint2021arXiv

SUOD: Accelerating Large-Scale Unsupervised Heterogeneous Outlier Detection

Outlier detection (OD) is a key machine learning (ML) task for identifying abnormal objects from general samples with numerous high-stake applications including fraud detection and intrusion detection. Due to the lack of ground truth labels, practitioners often have to build a large number of unsupervised, heterogeneous models (i.e., different algorithms with varying hyperparameters) for further combination and analysis, rather than relying on a single model. How to accelerate the training and scoring on new-coming samples by outlyingness (referred as prediction throughout the paper) with a large number of unsupervised, heterogeneous OD models? In this study, we propose a modular acceleration system, called SUOD, to address it. The proposed system focuses on three complementary acceleration aspects (data reduction for high-dimensional data, approximation for costly models, and taskload imbalance optimization for distributed environment), while maintaining performance accuracy. Extensive experiments on more than 20 benchmark datasets demonstrate SUOD's effectiveness in heterogeneous OD acceleration, along with a real-world deployment case on fraudulent claim analysis at IQVIA, a leading healthcare firm. We open-source SUOD for reproducibility and accessibility.

preprint2021arXiv

Superscalability of the random batch Ewald method

Coulomb interaction, following an inverse-square force-law, quantifies the amount of force between two stationary and electrically charged particles. The long-range nature of Coulomb interactions poses a major challenge to molecular dynamics simulations which are major tools for problems at the nano-/micro- scale. Various algorithms are developed to calculate the pairwise Coulomb interactions to a linear scaling but the poor scalability limits the size of simulated systems. Here, we conduct an efficient molecular dynamics algorithm with the random batch Ewald method on all-atom systems where the complete Fourier components in the Coulomb interaction are replaced by randomly selected mini-batches. By simulating the $N$-body systems up to 100 million particles using $10$ thousand CPU cores, we show that this algorithm furnishes $O(N)$ complexity, almost perfect scalability and an order of magnitude faster computational speed when compared to the existing state-of-the-art algorithms. Further examinations of our algorithm on distinct systems, including pure water, micro-phase-separated electrolyte and protein solution demonstrate that the spatiotemporal information on all time and length scales investigated and thermodynamic quantities derived from our algorithm are in perfect agreement with those obtained from the existing algorithms. Therefore, our algorithm provides a breakthrough solution on scalability of computing the Coulomb interaction. It is particularly useful and cost-effective to simulate ultra-large systems, which was either impossible or very costing to conduct using existing algorithms, thus would benefit the broad community of sciences.

preprint2020arXiv

Analytical solution for the surface states of antiferromagnetic topological insulator MnBi$_2$Te$_4$

Recently, the intrinsic magnetic topological insulator MnBi$_2$Te$_4$ has attracted great attention. It has an out-of-plane antiferromagnetic order, which is believed to open a sizable energy gap in the surface states. This gap, however, was not always observable in the latest angle-resolved photoemission spectroscopy (ARPES) experiments. To address this issue, we analytically derive an effective model for the two-dimensional (2D) surface states by starting from a three-dimensional (3D) Hamiltonian for bulk MnBi$_2$Te$_4$ and taking into account the spatial profile of the bulk magnetization. Our calculations suggest that the diminished surface gap may be caused by a much smaller and more localized intralayer ferromagnetic order. In addition, we calculate the spatial distribution and penetration depth of the surface states, which indicates that the surface states are mainly embedded in the first two septuple layers from the terminating surface. From our analytical results, the influence of the bulk parameters on the surface states can be found explicitly. Furthermore, we derive a $\bf{k}\cdot \bf{p}$ model for MnBi$_2$Te$_4$ thin films and show the oscillation of the Chern number between odd and even septuple layers. Our results will be helpful for the ongoing explorations of the MnBi$_x$Te$_y$ family.

preprint2020arXiv

Baryonic and Leptonic GeV Dark Matter

We perform a systematic analysis of models with GeV-scale dark matter coupled to baryons and leptons. Such theories provide a natural framework to explain the matter-antimatter asymmetry of the universe. We find that only a few baryonic dark matter models are free from tree-level proton decay without explicitly imposing baryon number conservation. We enumerate those cases and provide a brief overview of their phenomenology. We then focus on a leptonic dark matter model for a more detailed discussion of the baryon asymmetry generation via leptogenesis, the symmetry restoration in the dark sector and the expected dark matter annihilation signals in indirect detection experiments.

preprint2020arXiv

COPOD: Copula-Based Outlier Detection

Outlier detection refers to the identification of rare items that are deviant from the general data distribution. Existing approaches suffer from high computational complexity, low predictive capability, and limited interpretability. As a remedy, we present a novel outlier detection algorithm called COPOD, which is inspired by copulas for modeling multivariate data distribution. COPOD first constructs an empirical copula, and then uses it to predict tail probabilities of each given data point to determine its level of "extremeness". Intuitively, we think of this as calculating an anomalous p-value. This makes COPOD both parameter-free, highly interpretable, and computationally efficient. In this work, we make three key contributions, 1) propose a novel, parameter-free outlier detection algorithm with both great performance and interpretability, 2) perform extensive experiments on 30 benchmark datasets to show that COPOD outperforms in most cases and is also one of the fastest algorithms, and 3) release an easy-to-use Python implementation for reproducibility.

preprint2020arXiv

Distinct Topological Surface States on the Two Terminations of MnBi$_4$Te$_7$

The recent discovered intrinsic magnetic topological insulator MnBi2Te4 have been met with unusual success in hosting emergent phenomena such as the quantum anomalous Hall effect and the axion insulator states. However, the surface-bulk correspondence of the Mn-Bi-Te family, composed by the superlattice-like MnBi2Te4/(Bi2Te3)n (n = 0, 1, 2, 3 ...) layered structure, remains intriguing but elusive. Here, by using scanning tunneling microscopy (STM) and angle-resolved photoemission spectroscopy (ARPES) techniques, we unambiguously assign the two distinct surface states of MnBi4Te7 (n = 1) to the quintuple-layer (QL) Bi2Te3 termination and the septuple-layer (SL) MnBi2Te4 termination, respectively. A comparison of the experimental observations with theoretical calculations reveals the diverging topological behaviors, especially the hybridization effect between magnetic and nonmagnetic layers, on the two terminations: a gap on the QL termination originating from the topological surface states of the QL hybridizing with the bands of the beneath SL, and a gapless Dirac-cone band structure on the SL termination with time-reversal symmetry. The quasi-particle interference patterns further confirm the topological nature of the surface states for both terminations, continuing far above the Fermi energy. The QL termination carries a spin-helical Dirac state with hexagonal warping, while at the SL termination, a strongly canted helical state from the surface lies between a pair of Rashba-split states from its neighboring layer. Our work elucidates an unprecedented hybridization effect between the building blocks of the topological surface states, and also reveals the termination-dependent time-reversal symmetry breaking in a magnetic topological insulator, rendering an ideal platform to realize the half-integer quantum Hall effect and relevant quantum phenomena.

preprint2020arXiv

Epitaxial Growth of $β$-Ga$_2$O$_3$ Coated Wide Bandgap Semiconductor Tape for Flexible UV Photodetector

The epitaxial growth of technically-important $β$-Ga$_2$O$_3$ semiconductor thin films have not been realized on flexible substrates due to limitations by the high-temperature crystallization conditions and the lattice-matching requirements. In this report, for the first time single crystal $β$-Ga$_2$O$_3$(-201) thin films is epitaxially grown on the flexible CeO2 (001)-buffered hastelloy tape. The results indicate that CeO$_2$ (001) has a small bi-axial lattice mismatch with $β$-Ga$_2$O$_3$ (-201), thus inducing a simultaneous double-domain epitaxial growth. Flexible photodetectors are fabricated based on the epitaxial $β$-Ga$_2$O$_3$ coated tapes. Measurements show that the obtained photodetectors have a responsivity of 40 mA/W, with an on/off ratio reaching 1000 under 250 nm incident light and 5 V bias voltage. Such photoelectrical performance is already within the mainstream level of the $β$-Ga$_2$O$_3$ based photodetectors by using the conventional rigid single crystal substrates; and more importantly remained robust against more than 1000 cycles of bending tests. In addition, the epitaxy technique described in the report also paves the way for the fabrication of a wide range of flexible epitaxial film devices that utilize the materials with lattice parameters similar to $β$-Ga$_2$O$_3$, including GaN, AlN and SiC.

preprint2020arXiv

FineGym: A Hierarchical Video Dataset for Fine-grained Action Understanding

On public benchmarks, current action recognition techniques have achieved great success. However, when used in real-world applications, e.g. sport analysis, which requires the capability of parsing an activity into phases and differentiating between subtly different actions, their performances remain far from being satisfactory. To take action recognition to a new level, we develop FineGym, a new dataset built on top of gymnastic videos. Compared to existing action recognition datasets, FineGym is distinguished in richness, quality, and diversity. In particular, it provides temporal annotations at both action and sub-action levels with a three-level semantic hierarchy. For example, a "balance beam" event will be annotated as a sequence of elementary sub-actions derived from five sets: "leap-jump-hop", "beam-turns", "flight-salto", "flight-handspring", and "dismount", where the sub-action in each set will be further annotated with finely defined class labels. This new level of granularity presents significant challenges for action recognition, e.g. how to parse the temporal structures from a coherent action, and how to distinguish between subtly different action classes. We systematically investigate representative methods on this dataset and obtain a number of interesting findings. We hope this dataset could advance research towards action understanding.

preprint2020arXiv

Half-Magnetic Topological Insulator

Topological magnets are a new family of quantum materials providing great potential to realize emergent phenomena, such as quantum anomalous Hall effect and axion-insulator state. Here we present our discovery that stoichiometric ferromagnet MnBi8Te13 with natural heterostructure MnBi2Te4-(Bi2Te3)3 is an unprecedented half-magnetic topological insulator, with the magnetization existing at the MnBi2Te4 surface but not at the opposite surface terminated by triple Bi2Te3 layers. Our angle-resolved photoemission spectroscopy measurements unveil a massive Dirac gap at the MnBi2Te4 surface, and gapless Dirac cone on the other side. Remarkably, the Dirac gap (~28 meV) at MnBi2Te4 surface decreases monotonically with increasing temperature and closes right at the Curie temperature, thereby representing the first smoking-gun spectroscopic evidence of magnetization-induced topological surface gap among all known magnetic topological materials. We further demonstrate theoretically that the half-magnetic topological insulator is desirable to realize the half-quantized surface anomalous Hall effect, which serves as a direct proof of the general concept of axion electrodynamics in condensed matter systems.

preprint2020arXiv

Intra- and Inter-Action Understanding via Temporal Action Parsing

Current methods for action recognition primarily rely on deep convolutional networks to derive feature embeddings of visual and motion features. While these methods have demonstrated remarkable performance on standard benchmarks, we are still in need of a better understanding as to how the videos, in particular their internal structures, relate to high-level semantics, which may lead to benefits in multiple aspects, e.g. interpretable predictions and even new methods that can take the recognition performances to a next level. Towards this goal, we construct TAPOS, a new dataset developed on sport videos with manual annotations of sub-actions, and conduct a study on temporal action parsing on top. Our study shows that a sport activity usually consists of multiple sub-actions and that the awareness of such temporal structures is beneficial to action recognition. We also investigate a number of temporal parsing methods, and thereon devise an improved method that is capable of mining sub-actions from training data without knowing the labels of them. On the constructed TAPOS, the proposed method is shown to reveal intra-action information, i.e. how action instances are made of sub-actions, and inter-action information, i.e. one specific sub-action may commonly appear in various actions.

preprint2020arXiv

Mechanical Quantum Sensing in the Search for Dark Matter

Numerous astrophysical and cosmological observations are best explained by the existence of dark matter, a mass density which interacts only very weakly with visible, baryonic matter. Searching for the extremely weak signals produced by this dark matter strongly motivate the development of new, ultra-sensitive detector technologies. Paradigmatic advances in the control and readout of massive mechanical systems, in both the classical and quantum regimes, have enabled unprecedented levels of sensitivity. In this white paper, we outline recent ideas in the potential use of a range of solid-state mechanical sensing technologies to aid in the search for dark matter in a number of energy scales and with a variety of coupling mechanisms.

preprint2020arXiv

Omni-sourced Webly-supervised Learning for Video Recognition

We introduce OmniSource, a novel framework for leveraging web data to train video recognition models. OmniSource overcomes the barriers between data formats, such as images, short videos, and long untrimmed videos for webly-supervised learning. First, data samples with multiple formats, curated by task-specific data collection and automatically filtered by a teacher model, are transformed into a unified form. Then a joint-training strategy is proposed to deal with the domain gaps between multiple data sources and formats in webly-supervised learning. Several good practices, including data balancing, resampling, and cross-dataset mixup are adopted in joint training. Experiments show that by utilizing data from multiple sources and formats, OmniSource is more data-efficient in training. With only 3.5M images and 800K minutes videos crawled from the internet without human labeling (less than 2% of prior works), our models learned with OmniSource improve Top-1 accuracy of 2D- and 3D-ConvNet baseline models by 3.0% and 3.9%, respectively, on the Kinetics-400 benchmark. With OmniSource, we establish new records with different pretraining strategies for video recognition. Our best models achieve 80.4%, 80.5%, and 83.6 Top-1 accuracies on the Kinetics-400 benchmark respectively for training-from-scratch, ImageNet pre-training and IG-65M pre-training.

preprint2020arXiv

On Isometry Robustness of Deep 3D Point Cloud Models under Adversarial Attacks

While deep learning in 3D domain has achieved revolutionary performance in many tasks, the robustness of these models has not been sufficiently studied or explored. Regarding the 3D adversarial samples, most existing works focus on manipulation of local points, which may fail to invoke the global geometry properties, like robustness under linear projection that preserves the Euclidean distance, i.e., isometry. In this work, we show that existing state-of-the-art deep 3D models are extremely vulnerable to isometry transformations. Armed with the Thompson Sampling, we develop a black-box attack with success rate over 95% on ModelNet40 data set. Incorporating with the Restricted Isometry Property, we propose a novel framework of white-box attack on top of spectral norm based perturbation. In contrast to previous works, our adversarial samples are experimentally shown to be strongly transferable. Evaluated on a sequence of prevailing 3D models, our white-box attack achieves success rates from 98.88% to 100%. It maintains a successful attack rate over 95% even within an imperceptible rotation range $[\pm 2.81^{\circ}]$.

preprint2020arXiv

Probing Axions with Event Horizon Telescope Polarimetric Measurements

With high spatial resolution, polarimetric imaging of a supermassive black hole, like M87$^\star$ or Sgr A$^\star$, by the Event Horizon Telescope can be used to probe the existence of ultralight bosonic particles, such as axions. Such particles can accumulate around a rotating black hole through the superradiance mechanism, forming an axion cloud. When linearly polarized photons are emitted from an accretion disk near the horizon, their position angles oscillate due to the birefringent effect when traveling through the axion background. In particular, the observations of supermassive black holes M87$^\star$ (Sgr A$^\star$) can probe the dimensionless axion-photon coupling $c = 2 πg_{a γ} f_a$ for axions with mass around $O(10^{-20})$~eV ($O( 10^{-17}$)~eV) and decay constant $f_a < O(10^{16})$ GeV, which is complimentary to other axion measurements.

preprint2020arXiv

Spectral analysis of the quiescent low-mass X-ray binary in the globular cluster M30

We present a recent Chandra observation of the quiescent low-mass X-ray binary containing a neutron star, located in the globular cluster M30. We fit the thermal emission from the neutron star to extract its mass and radius. We find no evidence of flux variability between the two observations taken in 2001 and 2017, nor between individual 2017 observations, so we analyse them together to increase the signal to noise. We perform simultaneous spectral fits using standard light-element composition atmosphere models (hydrogen or helium), including absorption by the interstellar medium, correction for pile-up of X-ray photons on the detector, and a power-law for count excesses at high photon energy. Using a Markov-chain Monte Carlo approach, we extract mass and radius credible intervals for both chemical compositions of the atmosphere: $R_{\textrm{NS}}=7.94^{+0.76}_{-1.21}$ km and $M_{\textrm{NS}}<1.19$ M$_{\odot}$ assuming pure hydrogen, and $R_{\textrm{NS}}=10.50^{+2.88}_{-2.03}$ km and $M_{\textrm{NS}}<1.78$ M$_{\odot}$ for helium, where the uncertainties represent the 90% credible regions. For H, the small radius is difficult to reconcile with most current nuclear physics models (especially for nucleonic equations of state) and with other measurements of neutron star radii, with recent preferred values generally in the 11-14 km range. Whereas for He, the measured radius is consistent with this range. We discuss possible sources of systematic uncertainty that may result in an underestimation of the radius, identifying the presence of surface temperature inhomogeneities as the most relevant bias. According to this, we conclude that either the atmosphere is composed of He, or it is a H atmosphere with a significant contribution of hot spots to the observed radiation.

preprint2020arXiv

Spontaneous Surface Collapse and Reconstruction in Antiferromagnetic Topological Insulator MnBi$_2$Te$_4$

MnBi$_2$Te$_4$ is an antiferromagnetic topological insulator which stimulates intense interests due to the exotic quantum phenomena and promising device applications. Surface structure is a determinant factor to understand the novel magnetic and topological behavior of MnBi2Te4, yet its precise atomic structure remains elusive. Here, we discovered a spontaneous surface collapse and reconstruction in few-layer MnBi2Te4 exfoliated under delicate protection. Instead of the ideal septuple-layer structure in the bulk, the collapsed surface is shown to reconstruct as Mn-doped Bi$_2$Te$_3$ quintuple-layer and Mn$_x$Bi$_y$Te double-layer with a clear van der Waals gap in between. Combining with first-principles calculations, such spontaneous surface collapse is attributed to the abundant intrinsic Mn-Bi antisite defects and tellurium vacancy in the exfoliated surface, which is further supported by in-situ annealing and electron irradiation experiments. Our results shed light on the understanding of the intricate surface-bulk correspondence of MnBi$_2$Te$_4$, and provide insightful perspective of the surface-related quantum measurements in MnBi$_2$Te$_4$ few-layer devices.

preprint2020arXiv

SUOD: Toward Scalable Unsupervised Outlier Detection

Outlier detection is a key field of machine learning for identifying abnormal data objects. Due to the high expense of acquiring ground truth, unsupervised models are often chosen in practice. To compensate for the unstable nature of unsupervised algorithms, practitioners from high-stakes fields like finance, health, and security, prefer to build a large number of models for further combination and analysis. However, this poses scalability challenges in high-dimensional large datasets. In this study, we propose a three-module acceleration framework called SUOD to expedite the training and prediction with a large number of unsupervised detection models. SUOD's Random Projection module can generate lower subspaces for high-dimensional datasets while reserving their distance relationship. Balanced Parallel Scheduling module can forecast the training and prediction cost of models with high confidence---so the task scheduler could assign nearly equal amount of taskload among workers for efficient parallelization. SUOD also comes with a Pseudo-supervised Approximation module, which can approximate fitted unsupervised models by lower time complexity supervised regressors for fast prediction on unseen data. It may be considered as an unsupervised model knowledge distillation process. Notably, all three modules are independent with great flexibility to "mix and match"; a combination of modules can be chosen based on use cases. Extensive experiments on more than 30 benchmark datasets have shown the efficacy of SUOD, and a comprehensive future development plan is also presented.

preprint2020arXiv

SYNC: A Copula based Framework for Generating Synthetic Data from Aggregated Sources

A synthetic dataset is a data object that is generated programmatically, and it may be valuable to creating a single dataset from multiple sources when direct collection is difficult or costly. Although it is a fundamental step for many data science tasks, an efficient and standard framework is absent. In this paper, we study a specific synthetic data generation task called downscaling, a procedure to infer high-resolution, harder-to-collect information (e.g., individual level records) from many low-resolution, easy-to-collect sources, and propose a multi-stage framework called SYNC (Synthetic Data Generation via Gaussian Copula). For given low-resolution datasets, the central idea of SYNC is to fit Gaussian copula models to each of the low-resolution datasets in order to correctly capture dependencies and marginal distributions, and then sample from the fitted models to obtain the desired high-resolution subsets. Predictive models are then used to merge sampled subsets into one, and finally, sampled datasets are scaled according to low-resolution marginal constraints. We make four key contributions in this work: 1) propose a novel framework for generating individual level data from aggregated data sources by combining state-of-the-art machine learning and statistical techniques, 2) perform simulation studies to validate SYNC's performance as a synthetic data generation algorithm, 3) demonstrate its value as a feature engineering tool, as well as an alternative to data collection in situations where gathering is difficult through two real-world datasets, 4) release an easy-to-use framework implementation for reproducibility and scalability at the production level that easily incorporates new data.

preprint2020arXiv

The MAVERIC survey: A hidden pulsar and a black hole candidate in ATCA radio imaging of the globular cluster NGC 6397

Using a 16.2 hr radio observation by the Australia Telescope Compact Array (ATCA) and archival Chandra data, we found $>5σ$ radio counterparts to 4 known and 3 new X-ray sources within the half-light radius ($r_\mathrm{h}$) of the Galactic globular cluster NGC 6397. The previously suggested millisecond pulsar (MSP) candidate, U18, is a steep-spectrum ($S_ν\propto ν^α$; $α=-2.0^{+0.4}_{-0.5}$) radio source with a 5.5 GHz flux density of $54.7\pm 4.3~\mathrm{μJy}$. We argue that U18 is most likely a "hidden" MSP that is continuously hidden by plasma shocked at the collision between the winds from the pulsar and companion star. The nondetection of radio pulsations so far is probably the result of enhanced scattering in this shocked wind. On the other hand, we observed 5.5 GHz flux of the known MSP PSR J1740-5340 (U12) to decrease by a factor of $>2.8$ during epochs of 1.4 GHz eclipse, indicating that the radio flux is absorbed in its shocked wind. If U18 is indeed a pulsar whose pulsations are scattered, we note the contrast with U12's flux decrease in eclipse, which argues for two different eclipse mechanisms at the same radio frequency. In addition to U12 and U18, we also found radio associations for 5 other Chandra X-ray sources, four of which are likely background galaxies. The last, U97, which shows strong H$α$ variability, is mysterious; it may be either a quiescent black hole low-mass X-ray binary, or something more unusual.

preprint2020arXiv

Unsaturated Single Atoms on Monolayer Transition Metal Dichalcogenides for Ultrafast Hydrogen Evolution

Large scale implementation of electrochemical water splitting for hydrogen evolution requires cheap and efficient catalysts to replace expensive platinum. Molybdenum disulfide is one of the most promising alternative catalysts but its intrinsic activity is still inferior to platinum. There is therefore a need to explore new active site origins in molybdenum disulfide with ultrafast reaction kinetics and to understand their mechanisms. Here, we report a universal cold hydrogen plasma reduction method for synthesizing different single atoms sitting on two-dimensional monolayers. In case of molybdenum disulfide, we design and identify a new type of active site, i.e., unsaturated Mo single atoms on cogenetic monolayer molybdenum disulfide. The catalyst shows exceptional intrinsic activity with a Tafel slope of 35.1 mV dec-1 and a turnover frequency of ~10^3 s-1 at 100 mV, based on single flake microcell measurements. Theoretical studies indicate that coordinately unsaturated Mo single atoms sitting on molybdenum disulfide increase the bond strength between adsorbed hydrogen atoms and the substrates through hybridization, leading to fast hydrogen adsorption/desorption kinetics and superior hydrogen evolution activity. This work shines fresh light on preparing highly-efficient electrocatalysts for water splitting and other electrochemical processes, as well as provides a general method to synthesize single atoms on two-dimensional monolayers.

preprint2019arXiv

Combining Machine Learning Models using combo Library

Model combination, often regarded as a key sub-field of ensemble learning, has been widely used in both academic research and industry applications. To facilitate this process, we propose and implement an easy-to-use Python toolkit, combo, to aggregate models and scores under various scenarios, including classification, clustering, and anomaly detection. In a nutshell, combo provides a unified and consistent way to combine both raw and pretrained models from popular machine learning libraries, e.g., scikit-learn, XGBoost, and LightGBM. With accessibility and robustness in mind, combo is designed with detailed documentation, interactive examples, continuous integration, code coverage, and maintainability check; it can be installed easily through Python Package Index (PyPI) or https://github.com/yzhao062/combo.

preprint2019arXiv

Detecting Axion-like Dark Matter with Linearly Polarized Pulsar Light

Non-relativistic QCD axions or axion-like particles are among the most popular candidates for cold Dark Matter (DM) in the universe. We proposed to detect axion-like DM, using linearly polarized pulsar light as a probe. Because of birefringence effect potentially caused by an oscillating galactic axion DM background, when pulsar light travels across the galaxy, its linear polarization angle may vary with time. With a soliton+NFW galactic DM density profile, we show that this strategy can potentially probe an axion-photon coupling as small as $\sim 10^{-13}$ GeV$^{-1}$ for axion mass $m_a \sim 10^{-22}-10^{-20}$ eV, given the current measurement accuracy. An exclusion limit stronger than CAST ($ \sim 10^{-10}$ GeV$^{-1}$) and SN1987A ($ \sim 10^{-11}$ GeV$^{-1}$) could be extended up to $m_a \sim 10^{-18}$ eV and $\sim 10^{-19}$ eV, respectively.

preprint2019arXiv

Inverse scattering for the one-dimensional Helmholtz equation with piecewise constant wave speed

This paper analyzes inverse scattering for the one-dimensional Helmholtz equation in the case where the wave speed is piecewise constant. Scattering data recorded for an arbitrarily small interval of frequencies is shown to determine the wave speed uniquely, and a direct reconstruction algorithm is presented. The algorithm is exact provided data is recorded for a sufficiently wide range of frequencies and the jump points of the wave speed are equally spaced with respect to travel time. Numerical examples show that the algorithm works also in the general case of arbitrary wave speed (either with jumps or continuously varying etc.) giving progressively more accurate approximations as the range of recorded frequencies increases. A key underlying theoretical insight is to associate scattering data to compositions of automorphisms of the unit disk, which are in turn related to orthogonal polynomials on the unit circle. The algorithm exploits the three-term recurrence of orthogonal polynomials to reduce the required computation.

preprint2019arXiv

Magnetic order induced polarization anomaly of Raman scattering in 2D magnet CrI$_3$

The recent discovery of 2D magnets has revealed various intriguing phenomena due to the coupling between spin and other degree of freedoms (such as helical photoluminescence, nonreciprocal SHG). Previous research on the spin-phonon coupling effect mainly focuses on the renormalization of phonon frequency. Here we demonstrate that the Raman polarization selection rules of optical phonons can be greatly modified by the magnetic ordering in 2D magnet CrI$_3$. For monolayer samples, the dominant A$\rm_{1g}$ peak shows abnormally high intensity in the cross polarization channel at low temperature, which is forbidden by the selection rule based on the lattice symmetry. While for bilayer, this peak is absent in the cross polarization channel for the layered antiferromagnetic (AFM) state and reappears when it is tuned to the ferromagnetic (FM) state by an external magnetic field. Our findings shed light on exploring the emergent magneto-optical effects in 2D magnets.

preprint2019arXiv

Uniqueness to some inverse source problems for the wave equation in unbounded domains

This paper is concerned with inverse acoustic source problems in an unbounded domain with dynamical boundary surface data of Dirichlet kind. The measurement data are taken at a surface far away from the source support. We prove uniqueness in recovering source terms of the form $f(x)g(t)$ and $f(x_1,x_2,t) h(x_3)$, where $g(t)$ and $h(x_3)$ are given and $x=(x_1, x_2, x_3)$ is the spatial variable in three dimensions. Without these a priori information, we prove that the boundary data of a family of solutions can be used to recover general source terms depending on both time and spatial variables. For moving point sources radiating periodic signals, the data recorded at four receivers are proven sufficient to uniquely recover the orbit function. Simultaneous determination of embedded obstacles and source terms was verified in an inhomogeneous background medium using the observation data of infinite time period. Our approach depends heavily on the Laplace transform.

preprint2019arXiv

XGBOD: Improving Supervised Outlier Detection with Unsupervised Representation Learning

A new semi-supervised ensemble algorithm called XGBOD (Extreme Gradient Boosting Outlier Detection) is proposed, described and demonstrated for the enhanced detection of outliers from normal observations in various practical datasets. The proposed framework combines the strengths of both supervised and unsupervised machine learning methods by creating a hybrid approach that exploits each of their individual performance capabilities in outlier detection. XGBOD uses multiple unsupervised outlier mining algorithms to extract useful representations from the underlying data that augment the predictive capabilities of an embedded supervised classifier on an improved feature space. The novel approach is shown to provide superior performance in comparison to competing individual detectors, the full ensemble and two existing representation learning based algorithms across seven outlier datasets.

preprint2017arXiv

Mechanism of Re precipitation in irradiated W-Re alloys from kinetic Monte Carlo simulations

High-temperature, high-dose, neutron irradiation of W results in the formation of Re-rich clusters at concentrations one order of magnitude lower than the thermodynamic solubility limit. These clusters may eventually transform into brittle W-Re intermetallic phases, which can lead to high levels of hardening and thermal conductivity losses. Standard theories of radiation enhanced diffusion and precipitation cannot explain the formation of these precipitates and so understanding the mechanism by which nonequilibrium clusters form under irradiation is crucial to predict materials degradation and devise mitigation strategies. Here we carry out a thermodynamic study of W-Re alloys and conduct kinetic Monte Carlo simulations of Re cluster formation in irradiated W-2Re alloys using a generalized Hamiltonian for crystals containing point defects parameterized entirely with electronic structure calculations. Our model incorporates recently-gained mechanistic information of mixed-interstitial solute transport, which is seen to control cluster nucleation and growth by forming quasi-spherical nuclei after an average incubation time of 20 s at 1800 K. These nuclei are seen to grow by attracting more mixed interstitials bringing solute atoms, which in turns attracts vacancies leading to recombination and solute agglomeration. The clusters grow to a maximum size of approximately 4-nm radius, and are not fully dense with Re, containing 50% or less near the center. Our simulations are in reasonable agreement with recent atom probe examinations of ion irradiated W-2Re systems at 773 K.

preprint2016arXiv

Adaptive estimation of the copula correlation matrix for semiparametric elliptical copulas

We study the adaptive estimation of copula correlation matrix $Σ$ for the semi-parametric elliptical copula model. In this context, the correlations are connected to Kendall's tau through a sine function transformation. Hence, a natural estimate for $Σ$ is the plug-in estimator $\hatΣ$ with Kendall's tau statistic. We first obtain a sharp bound on the operator norm of $\hatΣ-Σ$. Then we study a factor model of $Σ$, for which we propose a refined estimator $\widetildeΣ$ by fitting a low-rank matrix plus a diagonal matrix to $\hatΣ$ using least squares with a nuclear norm penalty on the low-rank matrix. The bound on the operator norm of $\hatΣ-Σ$ serves to scale the penalty term, and we obtain finite sample oracle inequalities for $\widetildeΣ$. We also consider an elementary factor copula model of $Σ$, for which we propose closed-form estimators. All of our estimation procedures are entirely data-driven.

preprint2016arXiv

Detecting Superlight Dark Matter with Fermi-Degenerate Materials

We examine in greater detail the recent proposal of using superconductors for detecting dark matter as light as the warm dark matter limit of O(keV). Detection of such light dark matter is possible if the entire kinetic energy of the dark matter is extracted in the scattering, and if the experiment is sensitive to O(meV) energy depositions. This is the case for Fermi-degenerate materials in which the Fermi velocity exceeds the dark matter velocity dispersion in the Milky Way of ~10^-3. We focus on a concrete experimental proposal using a superconducting target with a transition edge sensor in order to detect the small energy deposits from the dark matter scatterings. Considering a wide variety of constraints, from dark matter self-interactions to the cosmic microwave background, we show that models consistent with cosmological/astrophysical and terrestrial constraints are observable with such detectors. A wider range of viable models with dark matter mass below an MeV is available if dark matter or mediator properties (such as couplings or masses) differ at BBN epoch or in stellar interiors from those in superconductors. We also show that metal targets pay a strong in-medium suppression for kinetically mixed mediators; this suppression is alleviated with insulating targets.

preprint2016arXiv

Distribution System Outage Detection using Consumer Load and Line Flow Measurements

An outage detection framework for power distribution networks is proposed. Given the tree structure of the distribution system, a method is developed combining the use of real-time power flow measurements on edges of the tree with load forecasts at the nodes of the tree. A maximum a posteriori detector {\color{black} (MAP)} is formulated for arbitrary number and location of outages on trees which is shown to have an efficient detector. A framework relying on the maximum missed detection probability is used for optimal sensor placement and is solved for tree networks. Finally, a set of case studies is considered using feeder data from the Pacific Northwest National Laboratories. We show that a 10\% loss in mean detection reliability network wide reduces the required sensor density by 60 \% for a typical feeder if efficient use of measurements is performed.

preprint2016arXiv

Explicitly Broken Supersymmetry with Exactly Massless Moduli

There is an avatar of the little hierarchy problem of the MSSM in 3-dimensional supersymmetry. We propose a solution to this problem in AdS$_3$ based on the AdS/CFT correspondence. The bulk theory is a supergravity theory in which U(1) $\times$ U(1) R-symmetry is gauged by Chern-Simons fields. The bulk theory is deformed by a boundary term quadratic in the gauge fields. It breaks SUSY completely and sources an exactly marginal operator in the dual CFT. SUSY breaking is communicated by gauge interactions to bulk scalar fields and their spinor superpartners. Since the R-charges of scalar and spinor differ, this generates a SUSY breaking shift of their masses. The Ward identity facilitates the calculation of these mass shifts to any desired order in the strength of the deformation. Moduli fields are massless $R$-neutral bulk scalars with vanishing potential in the undeformed theory. These properties are maintained to all orders in the deformation despite the fact that moduli couple in the bulk to loops of R-charged fields.

preprint2016arXiv

Minimum Sparsity of Unobservable Power Network Attacks

Physical security of power networks under power injection attacks that alter generation and loads is studied. The system operator employs Phasor Measurement Units (PMUs) for detecting such attacks, while attackers devise attacks that are unobservable by such PMU networks. It is shown that, given the PMU locations, the solution to finding the sparsest unobservable attacks has a simple form with probability one, namely, $κ(G^M) + 1$, where $κ(G^M)$ is defined as the vulnerable vertex connectivity of an augmented graph. The constructive proof allows one to find the entire set of the sparsest unobservable attacks in polynomial time. Furthermore, a notion of the potential impact of unobservable attacks is introduced. With optimized PMU deployment, the sparsest unobservable attacks and their potential impact as functions of the number of PMUs are evaluated numerically for the IEEE 30, 57, 118 and 300-bus systems and the Polish 2383, 2737 and 3012-bus systems. It is observed that, as more PMUs are added, the maximum potential impact among all the sparsest unobservable attacks drops quickly until it reaches the minimum sparsity.

preprint2016arXiv

Phase-locked array of quantum cascade lasers with an intracavity spatial filter

Phase-locking an array of quantum cascade lasers is an effective way to achieve higher output power and beam shaping. In this article, based on Talbot effect, we show a new-type phase-locked array of mid-infrared quantum cascade lasers with an integrated spatial- filtering Talbot cavity. All the arrays show stable in-phase operation from the threshold current to full power current. The beam divergence of the array device is smaller than that of a single-ridge laser. We use the multi-slit Fraunhofer diffraction mode to interpret the far-field radiation profile and give a solution to get better beam quality. The maximum power is just about 5 times that of a single-ridge laser for eleven-laser array device and 3 times for seven-laser array device. Considering the great modal selection ability, simple fabricating process and the potential for achieving better beam quality and smaller cavity loss, this new-type phase-locked array may be a hopeful and elegant solution to get high power or beam shaping.

preprint2015arXiv

3D Automatic Segmentation Method for Retinal Optical Coherence Tomography Volume Data Using Boundary Surface Enhancement

With the introduction of spectral-domain optical coherence tomography (SDOCT), much larger image datasets are routinely acquired compared to what was possible using the previous generation of time-domain OCT. Thus, there is a critical need for the development of 3D segmentation methods for processing these data. We present here a novel 3D automatic segmentation method for retinal OCT volume data. Briefly, to segment a boundary surface, two OCT volume datasets are obtained by using a 3D smoothing filter and a 3D differential filter. Their linear combination is then calculated to generate new volume data with an enhanced boundary surface, where pixel intensity, boundary position information, and intensity changes on both sides of the boundary surface are used simultaneously. Next, preliminary discrete boundary points are detected from the A-Scans of the volume data. Finally, surface smoothness constraints and a dynamic threshold are applied to obtain a smoothed boundary surface by correcting a small number of error points. Our method can extract retinal layer boundary surfaces sequentially with a decreasing search region of volume data. We performed automatic segmentation on eight human OCT volume datasets acquired from a commercial Spectralis OCT system, where each volume of data consisted of 97 OCT images with a resolution of 496 512; experimental results show that this method can accurately segment seven layer boundary surfaces in normal as well as some abnormal eyes.

preprint2015arXiv

A Radio for Hidden-Photon Dark Matter Detection

We propose a resonant electromagnetic detector to search for hidden-photon dark matter over an extensive range of masses. Hidden-photon dark matter can be described as a weakly coupled "hidden electric field," oscillating at a frequency fixed by the mass, and able to penetrate any shielding. At low frequencies (compared to the inverse size of the shielding), we find that observable effect of the hidden photon inside any shielding is a real, oscillating magnetic field. We outline experimental setups designed to search for hidden-photon dark matter, using a tunable, resonant LC circuit designed to couple to this magnetic field. Our "straw man" setups take into consideration resonator design, readout architecture and noise estimates. At high frequencies,there is an upper limit to the useful size of a single resonator set by $1/ν$. However, many resonators may be multiplexed within a hidden-photon coherence length to increase the sensitivity in this regime. Hidden-photon dark matter has an enormous range of possible frequencies, but current experiments search only over a few narrow pieces of that range. We find the potential sensitivity of our proposal is many orders of magnitude beyond current limits over an extensive range of frequencies, from 100 Hz up to 700 GHz and potentially higher.

preprint2015arXiv

AdS/CFT and the Little Hierarchy Problem

The AdS/CFT correspondence is applied to a close analogue of the little hierarchy problem in $AdS_{d+1}$, $d \geq 3$. The new mechanism requires a Maxwell field that gauges a $U(1)_R$ symmetry in a bulk supergravity theory with a negative cosmological constant. Supersymmetry is explicitly broken by a non-local boundary term with dimensionless coupling $h$. Non-locality appears to cause no pathology, and the SUSY breaking deformation engendered is exactly marginal. SUSY breaking effects in the bulk are computed using the $U(1)$ Ward identity. Conformal dimensions and thus masses of scalar and spinor partners are split simply because they carry different R-charges. However SUSY breaking effects cancel to all orders for R-neutral fields, even in diagrams with internal R-charged loops. SUSY breaking corrections can be summed to all orders in $h$. Diagrams involving graviton loops do not introduce any further SUSY breaking corrections. A possible scenario for a flat spacetime limit is outlined.

preprint2015arXiv

Creating and Probing Electron Whispering Gallery Modes in Graphene

Designing high-finesse resonant cavities for electronic waves faces challenges due to short electron coherence lengths in solids. Previous approaches, e.g. the seminal nanometer-sized quantum corrals, depend on careful positioning of adatoms at clean surfaces. Here we demonstrate an entirely different approach, inspired by the peculiar acoustic phenomena in whispering galleries. Taking advantage of graphene's unique properties, namely gate-tunable light-like carriers, we create Whispering Gallery Mode (WGM) resonators defined by circular pn-junctions, induced by a scanning tunneling probe. We can tune the resonator size and the carrier concentration under the probe in a back-gated graphene device over a wide range, independently and in situ. The confined modes, revealed through characteristic resonances in the tunneling spectrum, originate from Klein scattering at pn junction boundaries. The WGM-type confinement and resonances are a new addition to the quantum electron-optics toolbox, paving the way to real-world electronic lenses and resonators.

preprint2015arXiv

Dark Matter Annihilation Decay at The LHC

Collider experiments provide an opportunity to shed light on dark matter (DM) self-interactions. In this work, we study the possibility of generating DM bound states -- the Darkonium -- at the LHC and discuss how the annihilation decay of the Darkonium produces force carriers. We focus on two popular scenarios that contain large DM self-couplings: the Higgsinos in the $λ$-SUSY model, and self-interacting DM (SIDM) framework. After forming bound states, the DM particles annihilate into force mediators, which decay into the standard model particles either through a prompt or displaced process. This generates interesting signals for the heavy resonance search. We calculate the production rate of bound states and study the projected future constraints from the existing heavy resonance searches.

preprint2015arXiv

Superconducting Detectors for Super Light Dark Matter

We propose and study a new class of superconducting detectors which are sensitive to O(meV) electron recoils from dark matter-electron scattering. Such devices could detect dark matter as light as the warm dark matter limit, mX > keV. We compute the rate of dark matter scattering off of free electrons in a (superconducting) metal, including the relevant Pauli blocking factors. We demonstrate that classes of dark matter consistent with terrestrial and cosmological/astrophysical constraints could be detected by such detectors with a moderate size exposure.

preprint2015arXiv

Weak convergence of empirical copula processes indexed by functions

Weak convergence of the empirical copula process indexed by a class of functions is established. Two scenarios are considered in which either some smoothness of these functions or smoothness of the underlying copula function is required. A novel integration by parts formula for multivariate, right continuous functions of bounded variation, which is perhaps of independent interest, is proved. It is a key ingredient in proving weak convergence of a general empirical process indexed by functions of bounded variation.

preprint2014arXiv

A Parametrically Enhanced Hidden Photon Search

Many theories beyond the Standard Model contain hidden photons. A light hidden photon will generically couple to the Standard Model through a kinetic mixing term, giving a powerful avenue for detection using "Light-Shining-Through-A-Wall"-type transmission experiments with resonant cavities. We demonstrate a parametric enhancement of the signal in such experiments, resulting from transmission of the longitudinal mode of the hidden photon. While previous literature has focused on the production and detection of transverse modes, the longitudinal mode allows a significant improvement in experimental sensitivity. Although optical experiments such as ALPS are unable to take useful advantage of this enhancement, the reach of existing microwave cavity experiments such as CROWS is significantly enhanced beyond their published results. Future microwave cavity experiments, designed with appropriate geometry to take full advantage of the longitudinal mode, will provide a powerful probe of hidden-photon parameter space extending many orders of magnitude beyond current limits, including significant regions where the hidden photon can be dark matter.

preprint2014arXiv

Dark Matter Induced Nucleon Decay: Model and Signatures

If dark matter (DM) carries anti-baryon number, a DM particle may annihilate with a nucleon by flipping to anti-DM. Inspired by Hylogenesis models, we introduce a single component DM model where DM is asymmetric and carries B and L as -1/2. It can annihilate with a nucleon to an anti-lepton and an anti-DM at leading order or with an additional meson at sub-leading order. Such signals may be observed in proton decay experiments. If DM is captured in the Sun, the DM induced nucleon decay can generate a large flux of anti-neutrinos, which could be observed in neutrino experiments. Furthermore, the anti-DM particle in the final state obtains a relatively large momentum (few hundred MeV), and escapes the Sun. These fast-moving anti-DM particles could also induce interesting signals in various underground experiments.

preprint2014arXiv

Detecting Boosted Dark Matter from the Sun with Large Volume Neutrino Detectors

We study novel scenarios where thermal dark matter (DM) can be efficiently captured in the Sun and annihilate into boosted dark matter. In models with semi-annihilating DM, where DM has a non-minimal stabilization symmetry, or in models with a multi-component DM sector, annihilations of DM can give rise to stable dark sector particles with moderate Lorentz boosts. We investigate both of these possibilities, presenting concrete models as proofs of concept. Both scenarios can yield viable thermal relic DM with masses O(1)-O(100) GeV. Taking advantage of the energetic proton recoils that arise when the boosted DM scatters off matter, we propose a detection strategy which uses large volume terrestrial detectors, such as those designed to detect neutrinos or proton decays. In particular, we propose a search for proton tracks pointing towards the Sun. We focus on signals at Cherenkov-radiation-based detectors such as Super-Kamiokande (SK) and its upgrade Hyper-Kamiokande (HK). We find that with spin-dependent scattering as the dominant DM-nucleus interaction at low energies, boosted DM can leave detectable signals at SK or HK, with sensitivity comparable to DM direct detection experiments while being consistent with current constraints. Our study provides a new search path for DM sectors with non-minimal structure.

preprint2014arXiv

Indirect Detection Signatures for the Origin of Asymmetric Dark Matter

We study the decay signatures of Asymmetric Dark Matter (ADM) via higher dimension operators which are responsible for generating the primordial dark matter (DM) asymmetry. Since the signatures are sensitive both to the nature of the higher dimension operator generating the DM asymmetry and to the sign of the baryon or lepton number that the DM carries, indirect detection may provide a window into the nature of the mechanism which generates the DM asymmetry. We consider in particular dimension-6 fermionic operators of the form ${\cal O}_{ADM} = X {\cal O}_{B-L}/M^2$, where ${\cal O}_{B-L} = u^c d^c d^c,~\ell \ell e^c,~q \ell d^c$ (or operators related through a Hermitian conjugate) with the scale $M$ around or just below the GUT scale. We derive constraints on ADM particles both in the natural mass range (around a few GeV), as well as in the range between 100 GeV to 10 TeV. For light ADM, we focus on constraints from both the low energy gamma ray data and proton/anti-proton fluxes. For heavy ADM, we consider $γ$-rays and proton/anti-proton fluxes, and we fit $e^+/e^-$ data from AMS-02 and H.E.S.S. (neglecting the Fermi charged particle fluxes which disagree with AMS-02 below 100 GeV). We show that, although the best fit regions from electron/positron measurement are still in tension with other channels on account of the H.E.S.S. measurement at high energies, compared to an ordinary symmetric dark matter scenario, the decay of DM with a primordial asymmetry reduces the tension. Better measurement of the flux at high energy will be necessary to draw a definite conclusion about the viability of decaying DM as source for the signals.

preprint2014arXiv

Semiparametric Gaussian copula classification

This paper studies the binary classification of two distributions with the same Gaussian copula in high dimensions. Under this semiparametric Gaussian copula setting, we derive an accurate semiparametric estimator of the log density ratio, which leads to our empirical decision rule and a bound on its associated excess risk. Our estimation procedure takes advantage of the potential sparsity as well as the low noise condition in the problem, which allows us to achieve faster convergence rate of the excess risk than is possible in the existing literature on semiparametric Gaussian copula classification. We demonstrate the efficiency of our empirical decision rule by showing that the bound on the excess risk nearly achieves a convergence rate of $n^{-1/2}$ in the simple setting of Gaussian distribution classification.

preprint2013arXiv

A 2.0 Gb/s Throughput Decoder for QC-LDPC Convolutional Codes

This paper propose a decoder architecture for low-density parity-check convolutional code (LDPCCC). Specifically, the LDPCCC is derived from a quasi-cyclic (QC) LDPC block code. By making use of the quasi-cyclic structure, the proposed LDPCCC decoder adopts a dynamic message storage in the memory and uses a simple address controller. The decoder efficiently combines the memories in the pipelining processors into a large memory block so as to take advantage of the data-width of the embedded memory in a modern field-programmable gate array (FPGA). A rate-5/6 QC-LDPCCC has been implemented on an Altera Stratix FPGA. It achieves up to 2.0 Gb/s throughput with a clock frequency of 100 MHz. Moreover, the decoder displays an excellent error performance of lower than $10^{-13}$ at a bit-energy-to-noise-power-spectral-density ratio ($E_b/N_0$) of 3.55 dB.

preprint2013arXiv

A Low-Fluorine Solution with the F/Ba Mole Ratio of 2 for the Fabrication of YBCO Films

In the reported low-fluorine MOD-YBCO studies, the lowest F/Ba mole ratio of the precursor solution was 4.5. However, further lowering the F/Ba ratio is important according to the researches of YBCO thick film. On the other hand, the F/Ba ratio is necessary to be at least 2 for the full conversion of the Ba precursor to BaF_2 to avoid the formation of BaCO_3, which is detrimental to the superconducting performance. In this study, a novel solution with the F/Ba mole ratio of 2 was developed, in which the fluorine content was only about 10.3% of that used in the conventional TFA-MOD method. Attenuated total reflectance-Fourier transformed-infrared spectra(ATR-FT-IR) revealed that BaCO_3 was remarkably suppressed in the as-pyrolyzed film and eliminated at 700 Celsius degree. Thus YBCO films with a critical current density (J_c) over 5 MA cm^{-2} (77 K, 0 T, 200 nm thickness) could be obtained on LAO single crystal substrates. In-situ FT-IR spectra showed that no obvious fluorinated gaseous by-products were detected in the pyrolysis step, which indicated that all of the F atoms might remain in the film as fluorides. X-ray diffraction (XRD) θ/2θ-scan showed that BaF_2, but neither YF_3 nor CuF_2, was detected in the films quenched at 400 - 800 Celsius degree. The formation priority of BaF_2 over YF_3 and CuF_2 was interpreted by the chemical equilibrium of the potential reactions. Our study could enlarge the synthesis window of the precursor solution for MOD-YBCO fabrication and open a gate to study the fluorine content in the precursor solution continuously and systematically.

preprint2013arXiv

Induced Electroweak Symmetry Breaking and Supersymmetric Naturalness

In this paper we study a new class of supersymmetric models that can explain a 125 GeV Higgs without fine-tuning. These models contain additional `auxiliary Higgs' fields with large tree-level quartic interaction terms but no Yukawa couplings. These have electroweak-breaking vacuum expectation values, and contribute to the VEVs of the MSSM Higgs fields either through an induced quartic or through an induced tadpole. The quartic interactions for the auxiliary Higgs fields can arise from either D-terms or F-terms. The tadpole mechanism has been previously studied in strongly-coupled models with large D-terms, referred to as `superconformal technicolor.' The perturbative models studied here preserve gauge coupling unification in the simplest possible way, namely that all new fields are in complete SU(5) multiplets. The models are consistent with the observed properties of the 125 GeV Higgs-like boson as well as precision electroweak constraints, and predict a rich phenomenology of new Higgs states at the weak scale. The tuning is less than 10% in almost all of the phenomenologically allowed parameter space. If electroweak symmetry is broken by an induced tadpole, the cubic and quartic Higgs self-couplings are significantly smaller than in the standard model.

preprint2013arXiv

The Atacama Cosmology Telescope: Cosmological parameters from three seasons of data

We present constraints on cosmological and astrophysical parameters from high-resolution microwave background maps at 148 GHz and 218 GHz made by the Atacama Cosmology Telescope (ACT) in three seasons of observations from 2008 to 2010. A model of primary cosmological and secondary foreground parameters is fit to the map power spectra and lensing deflection power spectrum, including contributions from both the thermal Sunyaev-Zeldovich (tSZ) effect and the kinematic Sunyaev-Zeldovich (kSZ) effect, Poisson and correlated anisotropy from unresolved infrared sources, radio sources, and the correlation between the tSZ effect and infrared sources. The power ell^2 C_ell/2pi of the thermal SZ power spectrum at 148 GHz is measured to be 3.4 +\- 1.4 muK^2 at ell=3000, while the corresponding amplitude of the kinematic SZ power spectrum has a 95% confidence level upper limit of 8.6 muK^2. Combining ACT power spectra with the WMAP 7-year temperature and polarization power spectra, we find excellent consistency with the LCDM model. We constrain the number of effective relativistic degrees of freedom in the early universe to be Neff=2.79 +\- 0.56, in agreement with the canonical value of Neff=3.046 for three massless neutrinos. We constrain the sum of the neutrino masses to be Sigma m_nu < 0.39 eV at 95% confidence when combining ACT and WMAP 7-year data with BAO and Hubble constant measurements. We constrain the amount of primordial helium to be Yp = 0.225 +\- 0.034, and measure no variation in the fine structure constant alpha since recombination, with alpha/alpha0 = 1.004 +/- 0.005. We also find no evidence for any running of the scalar spectral index, dns/dlnk = -0.004 +\- 0.012.

preprint2013arXiv

The Naked-eye Optical Transient OT 120926

A previously unknown optical transient (OT 120926) has been observed in the constellation Bootes. The transient flared to magnitude 4.7, which is comparable to the visual magnitudes of the nearby stars $π$ Boo and $ο$ Boo. Database searches do not yield an unambiguous identification of a quiescent counterpart of this transient but do identify several candidates. However, none of the candidate stellar counterparts have shown any credible evidence of previous variability in the All-Sky Automated Survey or the Catalina Real-time Transient Survey. A flare on the nearby high proper motion, probable M dwarf star LP 440-48 could have produced OT 120926, but the amplitude of the flare would be an unprecedented 11.3 magnitudes. The current record amplitude for such flares on M dwarfs is 9.5 magnitudes.

preprint2012arXiv

A Complete Model of Low-Scale Gauge Mediation

Recent signs of a Standard Model-like Higgs at 125 GeV point towards large A-terms in the MSSM. This presents special challenges for gauge mediation, which by itself predicts vanishing A-terms at the messenger scale. In this paper, we review the general problems that arise when extending gauge mediation to achieve large A-terms, and the mechanisms that exist to overcome them. Using these mechanisms, we construct weakly-coupled models of low-scale gauge mediation with extended Higgs-messenger couplings that generate large A-terms at the messenger scale and viable mu/B_mu-terms. Our models are simple, economical, and complete realizations of supersymmetry at the weak scale.

preprint2012arXiv

Capacity Factors of a Point-to-point Network

In this paper, we investigate some properties on capacity factors, which were proposed to investigate the link failure problem from network coding. A capacity factor (CF) of a network is an edge set, deleting which will cause the maximum flow to decrease while deleting any proper subset will not. Generally, a $k$-CF is a minimal (not minimum) edge set which will cause the network maximum flow decrease by $k$. Under point to point acyclic scenario, we characterize all the edges which are contained in some CF, and propose an efficient algorithm to classify. And we show that all edges on some $s$-$t$ path in an acyclic point-to-point acyclic network are contained in some 2-CF. We also study some other properties of CF of point to point network, and a simple relationship with CF in multicast network. On the other hand, some computational hardness results relating to capacity factors are obtained. We prove that deciding whether there is a capacity factor of a cyclic network with size not less a given number is NP-complete, and the time complexity of calculating the capacity rank is lowered bounded by solving the maximal flow. Besides that, we propose the analogous definition of CF on vertices and show it captures edge capacity factors as a special case.

preprint2012arXiv

Implementation Of Decoders for LDPC Block Codes and LDPC Convolutional Codes Based on GPUs

With the use of belief propagation (BP) decoding algorithm, low-density parity-check (LDPC) codes can achieve near-Shannon limit performance. In order to evaluate the error performance of LDPC codes, simulators running on CPUs are commonly used. However, the time taken to evaluate LDPC codes with very good error performance is excessive. In this paper, efficient LDPC block-code decoders/simulators which run on graphics processing units (GPUs) are proposed. We also implement the decoder for the LDPC convolutional code (LDPCCC). The LDPCCC is derived from a pre-designed quasi-cyclic LDPC block code with good error performance. Compared to the decoder based on the randomly constructed LDPCCC code, the complexity of the proposed LDPCCC decoder is reduced due to the periodicity of the derived LDPCCC and the properties of the quasi-cyclic structure. In our proposed decoder architecture, $Γ$ (a multiple of a warp) codewords are decoded together and hence the messages of $Γ$ codewords are also processed together. Since all the $Γ$ codewords share the same Tanner graph, messages of the $Γ$ distinct codewords corresponding to the same edge can be grouped into one package and stored linearly. By optimizing the data structures of the messages used in the decoding process, both the read and write processes can be performed in a highly parallel manner by the GPUs. In addition, a thread hierarchy minimizing the divergence of the threads is deployed, and it can maximize the efficiency of the parallel execution. With the use of a large number of cores in the GPU to perform the simple computations simultaneously, our GPU-based LDPC decoder can obtain hundreds of times speedup compared with a serial CPU-based simulator and over 40 times speedup compared with an 8-thread CPU-based simulator.

preprint2012arXiv

Magnetoresistance measurements of Graphene at the Charge Neutrality Point

We report on transport measurements of the insulating state that forms at the charge neutrality point of graphene in a magnetic field. Using both conventional two-terminal measurements, sensitive to bulk and edge conductance, and Corbino measurements, sensitive only to the bulk conductance, we observed a vanishing conductance with increasing magnetic fields. By examining the resistance changes of this insulating state with varying perpendicular and in-plane fields, we probe the spin-active components of the excitations in total fields of up to 45 Tesla. Our results indicate that ν=0 quantum Hall state in single layer graphene is not spin polarized.

preprint2012arXiv

On PMU Location Selection for Line Outage Detection in Wide-area Transmission Networks

The optimal PMU locations to collect voltage phase angle measurements for detecting line outages in wide-area transmission networks are investigated. The problem is established as one of maximizing the minimum distance among the voltage phase angle signatures of the outages, which can be equivalently formulated as an integer programming problem. Based on a greedy heuristic and a linear programming relaxation, a branch and bound algorithm is proposed to find the globally optimal PMU locations. Using this algorithm, the optimal tradeoff between the number of PMUs and the outage detection performance is characterized for IEEE 14, 24 and 30 bus systems. The algorithm is shown to find the globally optimal PMU locations in a small number of iterations. It is observed that it is sufficient to have roughly one third of the buses providing PMU measurements in order to achieve the same outage detection performance as with all the buses providing PMU measurements.

preprint2012arXiv

Signature splitting inversion and backbending in 80Rb

High spin states of 80Rb are studied via the fusion-evaporation reactions 65Cu+19F, 66Zn+18O and 68Zn+16O with the beam energies of 75 MeV, 76 MeV and 80 MeV, respectively. Twenty-three new states with twenty-eight new γtransitions were added to the previously proposed level scheme, where the second negative-parity band is significantly pushed up to spins of 22^{-} and 15^{-} and two new sidebands are built on the known first negative-parity band. Two successive band crossings with frequencies 0.51 MeV and 0.61 MeV in the α=0 branch as well as another one in the α=1 branch of the second negative-parity band are observed for the first time. Signature inversions occur in the positive- and first negative-parity bands at the spins of 11\hbar and 15\hbar, respectively. The signature splitting is seen obviously in the second negative-parity band, but the signature inversion is not observed. It is also found that the structure of the two negative-parity bands is similar to that of its isotone ^{82}Y. Signature inversion in the positive-parity yrast band with configuration πg_{9/2} \otimes νg_{9/2} in this nucleus is discussed using the projected shell model (PSM).

preprint2012arXiv

The Atacama Cosmology Telescope: Data Characterization and Map Making

We present a description of the data reduction and mapmaking pipeline used for the 2008 observing season of the Atacama Cosmology Telescope (ACT). The data presented here at 148 GHz represent 12% of the 90 TB collected by ACT from 2007 to 2010. In 2008 we observed for 136 days, producing a total of 1423 hours of data (11 TB for the 148 GHz band only), with a daily average of 10.5 hours of observation. From these, 1085 hours were devoted to a 850 deg^2 stripe (11.2 hours by 9.1 deg) centered on a declination of -52.7 deg, while 175 hours were devoted to a 280 deg^2 stripe (4.5 hours by 4.8 deg) centered at the celestial equator. We discuss sources of statistical and systematic noise, calibration, telescope pointing, and data selection. Out of 1260 survey hours and 1024 detectors per array, 816 hours and 593 effective detectors remain after data selection for this frequency band, yielding a 38% survey efficiency. The total sensitivity in 2008, determined from the noise level between 5 Hz and 20 Hz in the time-ordered data stream (TOD), is 32 micro-Kelvin sqrt{s} in CMB units. Atmospheric brightness fluctuations constitute the main contaminant in the data and dominate the detector noise covariance at low frequencies in the TOD. The maps were made by solving the least-squares problem using the Preconditioned Conjugate Gradient method, incorporating the details of the detector and noise correlations. Cross-correlation with WMAP sky maps, as well as analysis from simulations, reveal that our maps are unbiased at multipoles ell > 300. This paper accompanies the public release of the 148 GHz southern stripe maps from 2008. The techniques described here will be applied to future maps and data releases.

preprint2011arXiv

Diagnosing the top-quark angular asymmetry using LHC intrinsic charge asymmetries

Flavor-violating interactions involving new heavy particles are among proposed explanations for the $t\bar t$ forward-backward asymmetry observed at the Tevatron. Many of these models generate a $t\bar t$-plus-jet signal at the LHC. In this paper we identify several new charge asymmetric variables in $t\bar tj$ events that can contribute to the discovery of such models at the LHC. We propose a data-driven method for the background, largely eliminating the need for a Monte Carlo prediction of $t\bar t$-plus-jets, and thus reducing systematic errors. With a fast detector simulation, we estimate the statistical sensitivity of our variables for one of these models, finding that charge-asymmetric variables could materially assist in the exclusion of the Standard Model across much of the mass and coupling range, given 5 inverse fb of data. Should any signal appear, our variables will be useful in distinguishing classes of models from one another.

preprint2011arXiv

Measurement of the ν= 1/3 fractional quantum Hall energy gap in suspended graphene

We report on magnetotransport measurements of multi-terminal suspended graphene devices. Fully developed integer quantum Hall states appear in magnetic fields as low as 2 T. At higher fields the formation of longitudinal resistance minima and transverse resistance plateaus are seen corresponding to fractional quantum Hall states, most strongly for ν= 1/3. By measuring the temperature dependence of these resistance minima, the energy gap for the 1/3 fractional state in graphene is determined to be at ~20 K at 14 T.

preprint2011arXiv

On the Sum-Capacity with Successive Decoding in Interference Channels

In this paper, we investigate the sum-capacity of the two-user Gaussian interference channel with Gaussian superposition coding and successive decoding. We first examine an approximate deterministic formulation of the problem, and introduce the complementarity conditions that capture the use of Gaussian coding and successive decoding. In the deterministic channel problem, we find the constrained sum-capacity and its achievable schemes with the minimum number of messages, first in symmetric channels, and then in general asymmetric channels. We show that the constrained sum-capacity oscillates as a function of the cross link gain parameters between the information theoretic sum-capacity and the sum-capacity with interference treated as noise. Furthermore, we show that if the number of messages of either of the two users is fewer than the minimum number required to achieve the constrained sum-capacity, the maximum achievable sum-rate drops to that with interference treated as noise. We provide two algorithms (a simple one and a finer one) to translate the optimal schemes in the deterministic channel model to the Gaussian channel model. We also derive two upper bounds on the sum-capacity of the Gaussian Han-Kobayashi schemes, which automatically upper bound the sum-capacity using successive decoding of Gaussian codewords. Numerical evaluations show that, similar to the deterministic channel results, the constrained sum-capacity in the Gaussian channels oscillates between the sum-capacity with Han-Kobayashi schemes and that with single message schemes.

preprint2011arXiv

Optimal Spectrum Management in Multiuser Interference Channels

In this paper, we study the non-convex problem of continuous frequency optimal spectrum management in multiuser frequency selective interference channels. Firstly, a simple pairwise channel condition for FDMA schemes to achieve all Pareto optimal points of the rate region is derived. It enables fully distributed global optimal decision making on whether any two users should use orthogonal channels. Next, we present in detail an analytical solution to finding the global optimum of sum-rate maximization in two-user symmetric flat channels. Generalizing this solution to frequency selective channels, a convex optimization is established that solves the global optimum. Finally, we show that our method generalizes to K-user (K>=2) weighted sum-rate maximization in asymmetric frequency selective channels, and transform this classic non-convex optimization in the primal domain to an equivalent convex optimization. The complexity is shown to be separable in its dependence on the channel parameters and the power constraints.

preprint2011arXiv

Recovering Particle Masses from Missing Energy Signatures with Displaced Tracks

If physics beyond the Standard Model contains metastable new particle states, its signatures may manifest themselves at the LHC through the presence of displaced vertices or displaced tracks. Such signatures are common in low-scale supersymmetry breaking scenarios as well as a host of other well-studied new physics models. Here we analyze the kinematics of dual cascade decays with missing energy resulting from a stable non-interacting new particle state at the bottom of the decay chain. We find that if the final step of the decay chain involves a metastable particle, the resulting presence of displaced tracks provides strong constraints with which to solve for the unknown kinematic quantities lost through missing energy. In addition we develop techniques with which to recover all kinematic unknowns, even in situations where the unknown quantities outnumber the constraints. These techniques can all be performed using a very small number of events and can thus be applied to very early discovery level searches at the LHC.

preprint2011arXiv

Simplified Models for LHC New Physics Searches

This document proposes a collection of simplified models relevant to the design of new-physics searches at the LHC and the characterization of their results. Both ATLAS and CMS have already presented some results in terms of simplified models, and we encourage them to continue and expand this effort, which supplements both signature-based results and benchmark model interpretations. A simplified model is defined by an effective Lagrangian describing the interactions of a small number of new particles. Simplified models can equally well be described by a small number of masses and cross-sections. These parameters are directly related to collider physics observables, making simplified models a particularly effective framework for evaluating searches and a useful starting point for characterizing positive signals of new physics. This document serves as an official summary of the results from the "Topologies for Early LHC Searches" workshop, held at SLAC in September of 2010, the purpose of which was to develop a set of representative models that can be used to cover all relevant phase space in experimental searches. Particular emphasis is placed on searches relevant for the first ~50-500 pb-1 of data and those motivated by supersymmetric models. This note largely summarizes material posted at http://lhcnewphysics.org/, which includes simplified model definitions, Monte Carlo material, and supporting contacts within the theory community. We also comment on future developments that may be useful as more data is gathered and analyzed by the experiments.

preprint2010arXiv

Entropic force and its fluctuation in Euclidian quantum gravity

In this paper, we study the idea about gravity as entropic force proposed by Verlinde. By applying the identification between partition functions of gravity theory and the dual field theory, we find the gravitational force can be calculated from the changing of entropy induced from its dual field theory side. Unlike Verlinde's calculation, our method does not depend on other assumptions such as equipartition of energy, so it is a safe check on the entropic force idea. Also we show how entropic force can be derived by similar method from gravity side along, which will provide a clue to a better understanding on the gauge/gravity duality. At last, we study the fluctuation of the force, and we find the fluctuation acting on the point-like particle can never be larger than the expectation value of the force. For black hole in AdS space, by gauge/gravity duality, the fluctuation is induced by thermal fluctuation from field theory living on conformal boundary. And for the metric of black hole in flat space, the ratio between fluctuation and force goes to a constant$\frac{T}{m}$ at infinity. .

preprint2010arXiv

The Atacama Cosmology Telescope: A Measurement of the Cosmic Microwave Background Power Spectrum at 148 and 218 GHz from the 2008 Southern Survey

We present measurements of the cosmic microwave background (CMB) power spectrum made by the Atacama Cosmology Telescope at 148 GHz and 218 GHz, as well as the cross-frequency spectrum between the two channels. Our results clearly show the second through the seventh acoustic peaks in the CMB power spectrum. The measurements of these higher-order peaks provide an additional test of the ΛCDM cosmological model. At l > 3000, we detect power in excess of the primary anisotropy spectrum of the CMB. At lower multipoles 500 < l < 3000, we find evidence for gravitational lensing of the CMB in the power spectrum at the 2.8σ level. We also detect a low level of Galactic dust in our maps, which demonstrates that we can recover known faint, diffuse signals.

preprint2010arXiv

The Atacama Cosmology Telescope: Calibration with WMAP Using Cross-Correlations

We present a new calibration method based on cross-correlations with WMAP and apply it to data from the Atacama Cosmology Telescope (ACT). ACT's observing strategy and map making procedure allows an unbiased reconstruction of the modes in the maps over a wide range of multipoles. By directly matching the ACT maps to WMAP observations in the multipole range of 400 < ell < 1000, we determine the absolute calibration with an uncertainty of 2% in temperature. The precise measurement of the calibration error directly impacts the uncertainties in the cosmological parameters estimated from the ACT power spectra. We also present a combined map based on ACT and WMAP data that has high signal-to-noise over a wide range of multipoles.

preprint2010arXiv

The Atacama Cosmology Telescope: Cosmology from Galaxy Clusters Detected via the Sunyaev-Zel'dovich Effect

We present constraints on cosmological parameters based on a sample of Sunyaev-Zel'dovich-selected galaxy clusters detected in a millimeter-wave survey by the Atacama Cosmology Telescope. The cluster sample used in this analysis consists of 9 optically-confirmed high-mass clusters comprising the high-significance end of the total cluster sample identified in 455 square degrees of sky surveyed during 2008 at 148 GHz. We focus on the most massive systems to reduce the degeneracy between unknown cluster astrophysics and cosmology derived from SZ surveys. We describe the scaling relation between cluster mass and SZ signal with a 4-parameter fit. Marginalizing over the values of the parameters in this fit with conservative priors gives sigma_8 = 0.851 +/- 0.115 and w = -1.14 +/- 0.35 for a spatially-flat wCDM cosmological model with WMAP 7-year priors on cosmological parameters. This gives a modest improvement in statistical uncertainty over WMAP 7-year constraints alone. Fixing the scaling relation between cluster mass and SZ signal to a fiducial relation obtained from numerical simulations and calibrated by X-ray observations, we find sigma_8 = 0.821 +/- 0.044 and w = -1.05 +/- 0.20. These results are consistent with constraints from WMAP 7 plus baryon acoustic oscillations plus type Ia supernoava which give sigma_8 = 0.802 +/- 0.038 and w = -0.98 +/- 0.053. A stacking analysis of the clusters in this sample compared to clusters simulated assuming the fiducial model also shows good agreement. These results suggest that, given the sample of clusters used here, both the astrophysics of massive clusters and the cosmological parameters derived from them are broadly consistent with current models.

preprint2009arXiv

Tuning the graphene work function by electric field effect

We report variation of the work function for single and bi-layer graphene devices measured by scanning Kelvin probe microscopy (SKPM). Using the electric field effect, the work function of graphene can be adjusted as the gate voltage tunes the Fermi level across the charge neutrality point. Upon biasing the device, the surface potential map obtained by SKPM provides a reliable way to measure the contact resistance of individual electrodes contacting graphene.