Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
75works
0followers
31topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

75 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&#39;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&#39;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&#39; 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 &#39;A-A&#39; 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&#39;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 &#34;short&#34; $ \lesssim 1~$\,s and &#34;long&#34; $ \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&#39;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&#39;s and Advanced Virgo&#39;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&#39;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&#39;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&#39;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 &#34;extremeness&#34;. 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 &#34;balance beam&#34; event will be annotated as a sequence of elementary sub-actions derived from five sets: &#34;leap-jump-hop&#34;, &#34;beam-turns&#34;, &#34;flight-salto&#34;, &#34;flight-handspring&#34;, and &#34;dismount&#34;, 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&#39;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 &#34;mix and match&#34;; 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&#39;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 &#34;hidden&#34; 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&#39;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.