Catalog footprint

What is connected

339works
77topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

339 published item(s)

preprint2026arXiv

PanoWorld: Towards Spatial Supersensing in 360$^\circ$ Panorama World

Multimodal large laboratory models (MLLMs) still struggle with spatial understanding under the dominant perspective-image paradigm, which inherits the narrow field of view of human-like perception. For navigation, robotic search, and 3D scene understanding, 360-degree panoramic sensing offers a form of supersensing by capturing the entire surrounding environment at once. However, existing MLLM pipelines typically decompose panoramas into multiple perspective views, leaving the spherical structure of equirectangular projection (ERP) largely implicit. In this paper, we study pano-native understanding, which requires an MLLM to reason over an ERP panorama as a continuous, observer-centered space. To this end, we first define the key abilities for pano-native understanding, including semantic anchoring, spherical localization, reference-frame transformation, and depth-aware 3D spatial reasoning. We then build a large-scale metadata construction pipeline that converts mixed-source ERP panoramas into geometry-aware, language-grounded, and depth-aware supervision, and instantiate these signals as capability-aligned instruction tuning data. On the model side, we introduce PanoWorld with Spherical Spatial Cross-Attention, which injects spherical geometry into the visual stream. We further construct PanoSpace-Bench, a diagnostic benchmark for evaluating ERP-native spatial reasoning. Experiments show that PanoWorld substantially outperforms both proprietary and open-source baselines on PanoSpace-Bench, H* Bench, and R2R-CE Val-Unseen benchmarks. These results demonstrate that robust panoramic reasoning requires dedicated pano-native supervision and geometry-aware model adaptation. All source code and proposed data will be publicly released.

preprint2025arXiv

Demonstration of Hardware Efficient Photonic Variational Quantum Algorithm

Quantum computing has brought a paradigm change in computer science, where non-classical technologies have promised to outperform their classical counterpart. Such an advantage was only demonstrated for tasks without practical applications, still out of reach for the state-of-art quantum technologies. In this context, a promising strategy to find practical use of quantum computers is to exploit hybrid quantum-classical models, where a quantum device estimates a hard-to-compute quantity, while a classical optimizer trains the parameters of the model. In this work, we demonstrate that single photons and linear optical networks are sufficient for implementing Variational Quantum Algorithms, when the problem specification, or ansatz, is tailored to this specific platform. We show this by a proof-of-principle demonstration of a variational approach to tackle an instance of a factorization task, whose solution is encoded in the ground state of a suitable Hamiltonian. This work which combines Variational Quantum Algorithms with hardware efficient ansatzes for linear-optics networks showcases a promising pathway towards practical applications for photonic quantum platforms.

preprint2025arXiv

OmniVCus: Feedforward Subject-driven Video Customization with Multimodal Control Conditions

Existing feedforward subject-driven video customization methods mainly study single-subject scenarios due to the difficulty of constructing multi-subject training data pairs. Another challenging problem that how to use the signals such as depth, mask, camera, and text prompts to control and edit the subject in the customized video is still less explored. In this paper, we first propose a data construction pipeline, VideoCus-Factory, to produce training data pairs for multi-subject customization from raw videos without labels and control signals such as depth-to-video and mask-to-video pairs. Based on our constructed data, we develop an Image-Video Transfer Mixed (IVTM) training with image editing data to enable instructive editing for the subject in the customized video. Then we propose a diffusion Transformer framework, OmniVCus, with two embedding mechanisms, Lottery Embedding (LE) and Temporally Aligned Embedding (TAE). LE enables inference with more subjects by using the training subjects to activate more frame embeddings. TAE encourages the generation process to extract guidance from temporally aligned control signals by assigning the same frame embeddings to the control and noise tokens. Experiments demonstrate that our method significantly surpasses state-of-the-art methods in both quantitative and qualitative evaluations. Video demos are at our project page: https://caiyuanhao1998.github.io/project/OmniVCus/. Our code, models, data are released at https://github.com/caiyuanhao1998/Open-OmniVCus

preprint2024arXiv

Depth Map Denoising Network and Lightweight Fusion Network for Enhanced 3D Face Recognition

With the increasing availability of consumer depth sensors, 3D face recognition (FR) has attracted more and more attention. However, the data acquired by these sensors are often coarse and noisy, making them impractical to use directly. In this paper, we introduce an innovative Depth map denoising network (DMDNet) based on the Denoising Implicit Image Function (DIIF) to reduce noise and enhance the quality of facial depth images for low-quality 3D FR. After generating clean depth faces using DMDNet, we further design a powerful recognition network called Lightweight Depth and Normal Fusion network (LDNFNet), which incorporates a multi-branch fusion block to learn unique and complementary features between different modalities such as depth and normal images. Comprehensive experiments conducted on four distinct low-quality databases demonstrate the effectiveness and robustness of our proposed methods. Furthermore, when combining DMDNet and LDNFNet, we achieve state-of-the-art results on the Lock3DFace database.

preprint2024arXiv

Transfer the linguistic representations from TTS to accent conversion with non-parallel data

Accent conversion aims to convert the accent of a source speech to a target accent, meanwhile preserving the speaker's identity. This paper introduces a novel non-autoregressive framework for accent conversion that learns accent-agnostic linguistic representations and employs them to convert the accent in the source speech. Specifically, the proposed system aligns speech representations with linguistic representations obtained from Text-to-Speech (TTS) systems, enabling training of the accent voice conversion model on non-parallel data. Furthermore, we investigate the effectiveness of a pretraining strategy on native data and different acoustic features within our proposed framework. We conduct a comprehensive evaluation using both subjective and objective metrics to assess the performance of our approach. The evaluation results highlight the benefits of the pretraining strategy and the incorporation of richer semantic features, resulting in significantly enhanced audio quality and intelligibility.

preprint2023arXiv

A Lightweight Clustering Framework for Unsupervised Semantic Segmentation

Unsupervised semantic segmentation aims to categorize each pixel in an image into a corresponding class without the use of annotated data. It is a widely researched area as obtaining labeled datasets is expensive. While previous works in the field have demonstrated a gradual improvement in model accuracy, most required neural network training. This made segmentation equally expensive, especially when dealing with large-scale datasets. We thus propose a lightweight clustering framework for unsupervised semantic segmentation. We discovered that attention features of the self-supervised Vision Transformer exhibit strong foreground-background differentiability. Therefore, clustering can be employed to effectively separate foreground and background image patches. In our framework, we first perform multilevel clustering across the Dataset-level, Category-level, and Image-level, and maintain consistency throughout. Then, the binary patch-level pseudo-masks extracted are upsampled, refined and finally labeled. Furthermore, we provide a comprehensive analysis of the self-supervised Vision Transformer features and a detailed comparison between DINO and DINOv2 to justify our claims. Our framework demonstrates great promise in unsupervised semantic segmentation and achieves state-of-the-art results on PASCAL VOC and MS COCO datasets.

preprint2023arXiv

High precision atom interferometer-based dynamic gravimeter measurement by eliminating the cross-coupling effect

A dynamic gravimeter with an atomic interferometer (AI) can perform absolute gravity measurements with high precision. AI-based dynamic gravity measurement is a type of joint measurement that uses AI sensors and a classical accelerometer. The coupling of the two sensors may degrade the measurement precision. In this study, we analyzed the cross-coupling effect and introduced a recovery vector to suppress this effect. We improved the phase noise of the interference fringe by a factor of 1.9 by performing marine gravity measurements using an AI-based gravimeter and optimizing the recovery vector. Marine gravity measurements were performed, and high gravity measurement precision was achieved. The external and inner coincidence accuracies of the gravity measurement are 0.42 mGal and 0.46 mGal, which were improved by factors of 4.18 and 4.21 by optimizing the cross-coupling effect.

preprint2023arXiv

LaFFi: Leveraging Hybrid Natural Language Feedback for Fine-tuning Language Models

Fine-tuning Large Language Models (LLMs) adapts a trained model to specific downstream tasks, significantly improving task-specific performance. Supervised Fine-Tuning (SFT) is a common approach, where an LLM is trained to produce desired answers. However, LLMs trained with SFT sometimes make simple mistakes and result in hallucinations on reasoning tasks such as question-answering. Without external feedback, it is difficult for SFT to learn a good mapping between the question and the desired answer, especially with a small dataset. This paper introduces an alternative to SFT called Natural Language Feedback for Finetuning LLMs (LaFFi). LaFFi has LLMs directly predict the feedback they will receive from an annotator. We find that requiring such reflection can significantly improve the accuracy in in-domain question-answering tasks, providing a promising direction for the application of natural language feedback in the realm of SFT LLMs. Additional ablation studies show that the portion of human-annotated data in the annotated datasets affects the fine-tuning performance.

preprint2023arXiv

Shortcuts to Adiabatic Soliton Compression in Active Nonlinear Kerr Media

We implement variational shortcuts to adiabaticity for optical pulse compression in an active nonlinear Kerr medium with distributed amplification and spatially varying dispersion and nonlinearity. Starting with the hyperbolic secant ansatz, we employ a variational approximation to systematically derive dynamical equations, establishing analytical relationships linking the amplitude, width, and chirp of the pulse. Through the inverse engineering approach, we manipulate the distributed gain/loss, nonlinearity and dispersion profiles to efficiently compress the optical pulse over a reduced distance with high fidelity. In addition, we explore the dynamical stability of the system to illustrate the advantage of our protocol over conventional adiabatic approaches. Finally, we analyze the impact of tailored higher-order dispersion on soliton self-compression and derive physical constraints on the final soliton width for the complementary case of soliton expansion. The broader implications of our findings extend beyond optical systems, encompassing areas such as cold-atom and magnetic systems highlighting the versatility and relevance of our approach in various physical contexts.

preprint2022arXiv

Accelerating Adaptive Cubic Regularization of Newton's Method via Random Sampling

In this paper, we consider an unconstrained optimization model where the objective is a sum of a large number of possibly nonconvex functions, though overall the objective is assumed to be smooth and convex. Our bid to solving such model uses the framework of cubic regularization of Newton's method. As well known, the crux in cubic regularization is its utilization of the Hessian information, which may be computationally expensive for large-scale problems. To tackle this, we resort to approximating the Hessian matrix via sub-sampling. In particular, we propose to compute an approximated Hessian matrix by either \textit{uniformly}\/ or \textit{non-uniformly}\/ sub-sampling the components of the objective. Based upon such sampling strategy, we develop accelerated adaptive cubic regularization approaches and provide theoretical guarantees on global iteration complexity of $O(ε^{-1/3})$ with high probability, which matches that of the original accelerated cubic regularization methods \cite{Jiang-2017-Unified} using the \textit{full}\/ Hessian information. Interestingly, we show that in the worst case scenario our algorithm still achieves an $O\left(\log(ε^{-1})ε^{-5/6}\right)$ iteration complexity bound. The performances of the proposed methods on the regularized logistic regression problems show a clear effect of acceleration in terms of the epoch counts on several real data sets.

preprint2022arXiv

Active Learning for Contextual Search with Binary Feedbacks

In this paper, we study the learning problem in contextual search, which is motivated by applications such as first-price auction, personalized medicine experiments, and feature-based pricing experiments. In particular, for a sequence of arriving context vectors, with each context associated with an underlying value, the decision-maker either makes a query at a certain point or skips the context. The decision-maker will only observe the binary feedback on the relationship between the query point and the value associated with the context. We study a PAC learning setting, where the goal is to learn the underlying mean value function in context with a minimum number of queries. To address this challenge, we propose a tri-section search approach combined with a margin-based active learning method. We show that the algorithm only needs to make $O(1/\varepsilon^2)$ queries to achieve an $ε$-estimation accuracy. This sample complexity significantly reduces the required sample complexity in the passive setting, at least $Ω(1/\varepsilon^4)$.

preprint2022arXiv

All You May Need for VQA are Image Captions

Visual Question Answering (VQA) has benefited from increasingly sophisticated models, but has not enjoyed the same level of engagement in terms of data creation. In this paper, we propose a method that automatically derives VQA examples at volume, by leveraging the abundance of existing image-caption annotations combined with neural models for textual question generation. We show that the resulting data is of high-quality. VQA models trained on our data improve state-of-the-art zero-shot accuracy by double digits and achieve a level of robustness that lacks in the same model trained on human-annotated VQA data.

preprint2022arXiv

AutoMO-Mixer: An automated multi-objective Mixer model for balanced, safe and robust prediction in medicine

Accurately identifying patient's status through medical images plays an important role in diagnosis and treatment. Artificial intelligence (AI), especially the deep learning, has achieved great success in many fields. However, more reliable AI model is needed in image guided diagnosis and therapy. To achieve this goal, developing a balanced, safe and robust model with a unified framework is desirable. In this study, a new unified model termed as automated multi-objective Mixer (AutoMO-Mixer) model was developed, which utilized a recent developed multiple layer perceptron Mixer (MLP-Mixer) as base. To build a balanced model, sensitivity and specificity were considered as the objective functions simultaneously in training stage. Meanwhile, a new evidential reasoning based on entropy was developed to achieve a safe and robust model in testing stage. The experiment on an optical coherence tomography dataset demonstrated that AutoMO-Mixer can obtain safer, more balanced, and robust results compared with MLP-Mixer and other available models.

preprint2022arXiv

Bayesian posterior repartitioning for nested sampling

Priors in Bayesian analyses often encode informative domain knowledge that can be useful in making the inference process more efficient. Occasionally, however, priors may be unrepresentative of the parameter values for a given dataset, which can result in inefficient parameter space exploration, or even incorrect inferences, particularly for nested sampling (NS) algorithms. Simply broadening the prior in such cases may be inappropriate or impossible in some applications. Hence our previous solution to this problem, known as posterior repartitioning (PR), redefines the prior and likelihood while keeping their product fixed, so that the posterior inferences and evidence estimates remain unchanged, but the efficiency of the NS process is significantly increased. In its most practical form, PR raises the prior to some power beta, which is introduced as an auxiliary variable that must be determined on a case-by-case basis, usually by lowering beta from unity according to some pre-defined `annealing schedule' until the resulting inferences converge to a consistent solution. Here we present a very simple yet powerful alternative Bayesian approach, in which beta is instead treated as a hyperparameter that is inferred from the data alongside the original parameters of the problem, and then marginalised over to obtain the final inference. We show through numerical examples that this Bayesian PR (BPR) method provides a very robust, self-adapting and computationally efficient `hands-off' solution to the problem of unrepresentative priors in Bayesian inference using NS. Moreover, unlike the original PR method, we show that even for representative priors BPR has a negligible computational overhead relative to standard nesting sampling, which suggests that it should be used as the default in all NS analyses.

preprint2022arXiv

Collaborative learning of images and geometrics for predicting isocitrate dehydrogenase status of glioma

The isocitrate dehydrogenase (IDH) gene mutation status is an important biomarker for glioma patients. The gold standard of IDH mutation detection requires tumour tissue obtained via invasive approaches and is usually expensive. Recent advancement in radiogenomics provides a non-invasive approach for predicting IDH mutation based on MRI. Meanwhile, tumor geometrics encompass crucial information for tumour phenotyping. Here we propose a collaborative learning framework that learns both tumor images and tumor geometrics using convolutional neural networks (CNN) and graph neural networks (GNN), respectively. Our results show that the proposed model outperforms the baseline model of 3D-DenseNet121. Further, the collaborative learning model achieves better performance than either the CNN or the GNN alone. The model interpretation shows that the CNN and GNN could identify common and unique regions of interest for IDH mutation prediction. In conclusion, collaborating image and geometric learners provides a novel approach for predicting genotype and characterising glioma.

preprint2022arXiv

Context-Based Dynamic Pricing with Online Clustering

We consider a context-based dynamic pricing problem of online products, which have low sales. Sales data from Alibaba, a major global online retailer, illustrate the prevalence of low-sale products. For these products, existing single-product dynamic pricing algorithms do not work well due to insufficient data samples. To address this challenge, we propose pricing policies that concurrently perform clustering over product demand and set individual pricing decisions on the fly. By clustering data and identifying products that have similar demand patterns, we utilize sales data from products within the same cluster to improve demand estimation for better pricing decisions. We evaluate the algorithms using regret, and the result shows that when product demand functions come from multiple clusters, our algorithms significantly outperform traditional single-product pricing policies. Numerical experiments using a real dataset from Alibaba demonstrate that the proposed policies, compared with several benchmark policies, increase the revenue. The results show that online clustering is an effective approach to tackling dynamic pricing problems associated with low-sale products.

preprint2022arXiv

Coverage Optimization of Camera Network for Continuous Deformable Object

In this paper, a deformable object is considered for cameras deployment with the aim of visual coverage. The object contour is discretized into sampled points as meshes, and the deformation is represented as continuous trajectories for the sampled points. To reduce the computational complexity, some feature points are carefully selected representing the continuous deformation process, and the visual coverage for the deformable object is transferred to cover the specific feature points. In particular, the vertexes of a rectangle that can contain the entire deformation trajectory of every sampled point on the object contour are chosen as the feature points. An improved wolf pack algorithm is then proposed to solve the optimization problem. Finally, simulation results are given to demonstrate the effectiveness of the proposed deployment method of camera network.

preprint2022arXiv

Development of a compact high-resolution absolute gravity gradiometer based on atom interferometers

We present a compact high-resolution gravity gradiometer based on dual Rb-85 atom interferometers using stimulated Raman transitions. A baseline L=44.5 cm and an interrogation time T=130 ms are realized in a sensor head with volume of less than 100 liters. Experimental parameters are optimized to improve the short-term sensitivity while a rejection algorithm relying on inversion of the Raman wave vector is implemented to improve the long-term stability. After an averaging time of 17000 s, a phase resolution of 104 μrad is achieved, which corresponds to a gravity gradient resolution of 0.86 E. As far as we know, this is the sub-E atom gravity gradiometer with the highest level of compactness to date. After the evaluation and correction of system errors induced by light shift, residual Zeeman shift, Coriolis effect and self-attraction effect, the instrument serves as an absolute gravity gradiometer and with it the local gravity gradient is measured to be 3114 (53) E.

preprint2022arXiv

Digitized-counterdiabatic quantum approximate optimization algorithm

The quantum approximate optimization algorithm (QAOA) has proved to be an effective classical-quantum algorithm serving multiple purposes, from solving combinatorial optimization problems to finding the ground state of many-body quantum systems. Since QAOA is an ansatz-dependent algorithm, there is always a need to design ansatz for better optimization. To this end, we propose a digitized version of QAOA enhanced via the use of shortcuts to adiabaticity. Specifically, we use a counterdiabatic (CD) driving term to design a better ansatz, along with the Hamiltonian and mixing terms, enhancing the global performance. We apply our digitized-counterdiabatic QAOA to Ising models, classical optimization problems, and the P-spin model, demonstrating that it outperforms standard QAOA in all cases we study.

preprint2022arXiv

Digitized-Counterdiabatic Quantum Optimization

We propose digitized-counterdiabatic quantum optimization (DCQO) to achieve polynomial enhancement over adiabatic quantum optimization for the general Ising spin-glass model, which includes the whole class of combinatorial optimization problems. This is accomplished via the digitization of adiabatic quantum algorithms that are catalysed by the addition of non-stoquastic counterdiabatic terms. The latter are suitably chosen, not only for escaping classical simulability, but also for speeding up the performance. Finding the ground state of a general Ising spin-glass Hamiltonian is used to illustrate that the inclusion of k-local non-stoquastic counterdiabatic terms can always outperform the traditional adiabatic quantum optimization with stoquastic Hamiltonians. In particular, we show that a polynomial enhancement in the ground-state success probability can be achieved for a finite-time evolution, even with the simplest 2-local counterdiabatic terms. Furthermore, the considered digitization process, within the gate-based quantum computing paradigm, provides the flexibility to introduce arbitrary non-stoquastic interactions. Along these lines, using our proposed paradigm on current NISQ computers, quantum speed-up may be reached to find approximate solutions for NP-complete and NP-hard optimization problems. We expect DCQO to become a fast-lane paradigm towards quantum advantage in the NISQ era.

preprint2022arXiv

Dynamic Car Dispatching and Pricing: Revenue and Fairness for Ridesharing Platforms

A major challenge for ridesharing platforms is to guarantee profit and fairness simultaneously, especially in the presence of misaligned incentives of drivers and riders. We focus on the dispatching-pricing problem to maximize the total revenue while keeping both drivers and riders satisfied. We study the computational complexity of the problem, provide a novel two-phased pricing solution with revenue and fairness guarantees, extend it to stochastic settings and develop a dynamic (a.k.a., learning-while-doing) algorithm that actively collects data to learn the demand distribution during the scheduling process. We also conduct extensive experiments to demonstrate the effectiveness of our algorithms.

preprint2022arXiv

Efficient broadband frequency conversion via shortcut to adiabaticity

The method of adiabatic frequency conversion, in analogy with the two level atomic system, has been put forward recently and verified experimentally to achieve robust frequency mixing processes such as sum and difference frequency generation. Here we present a comparative study of efficient frequency mixing using various techniques of shortcuts to adiabaticity (STA) such as counter-diabatic driving and invariant-based inverse engineering. We show that, it is possible to perform sum frequency generation by properly designing the poling structure of a periodically poled crystal and the coupling between the input lights and the crystal. The required crystal length for frequency conversion is significantly decreases beyond the adiabatic limit. Our approach significantly improves the robustness of the process against the variation in temperature as well as the signal frequency. By introducing a single parameter control technique with constant coupling and combining with the inverse engineering, perturbation theory and optimal control, we show that the phase mismatch can be further optimized with respect to the fluctuations of input wavelength and crystal temperature that results into a novel experimentally realizable mixing scheme.

preprint2022arXiv

Encoder-Decoder Architecture for Supervised Dynamic Graph Learning: A Survey

In recent years, the prevalent online services generate a sheer volume of user activity data. Service providers collect these data in order to perform client behavior analysis, and offer better and more customized services. Majority of these data can be modeled and stored as graph, such as the social graph in Facebook, user-video interaction graph in Youtube. These graphs need to evolve over time to capture the dynamics in the real world, leading to the invention of dynamic graphs. However, the temporal information embedded in the dynamic graphs brings new challenges in analyzing and deploying them. Events staleness, temporal information learning and explicit time dimension usage are some example challenges in dynamic graph learning. In order to offer a convenient reference to both the industry and academia, this survey presents the Three Stages Recurrent Temporal Learning Framework based on dynamic graph evolution theories, so as to interpret the learning of temporal information with a generalized framework. Under this framework, this survey categories and reviews different learnable encoder-decoder architectures for supervised dynamic graph learning. We believe that this survey could supply useful guidelines to researchers and engineers in finding suitable graph structures for their dynamic learning tasks.

preprint2022arXiv

Finding Influential Instances for Distantly Supervised Relation Extraction

Distant supervision (DS) is a strong way to expand the datasets for enhancing relation extraction (RE) models but often suffers from high label noise. Current works based on attention, reinforcement learning, or GAN are black-box models so they neither provide meaningful interpretation of sample selection in DS nor stability on different domains. On the contrary, this work proposes a novel model-agnostic instance sampling method for DS by influence function (IF), namely REIF. Our method identifies favorable/unfavorable instances in the bag based on IF, then does dynamic instance sampling. We design a fast influence sampling algorithm that reduces the computational complexity from $\mathcal{O}(mn)$ to $\mathcal{O}(1)$, with analyzing its robustness on the selected sampling function. Experiments show that by simply sampling the favorable instances during training, REIF is able to win over a series of baselines that have complicated architectures. We also demonstrate that REIF can support interpretable instance selection.

preprint2022arXiv

Fixed-Support Wasserstein Barycenters: Computational Hardness and Fast Algorithm

We study the fixed-support Wasserstein barycenter problem (FS-WBP), which consists in computing the Wasserstein barycenter of $m$ discrete probability measures supported on a finite metric space of size $n$. We show first that the constraint matrix arising from the standard linear programming (LP) representation of the FS-WBP is \textit{not totally unimodular} when $m \geq 3$ and $n \geq 3$. This result resolves an open question pertaining to the relationship between the FS-WBP and the minimum-cost flow (MCF) problem since it proves that the FS-WBP in the standard LP form is not an MCF problem when $m \geq 3$ and $n \geq 3$. We also develop a provably fast \textit{deterministic} variant of the celebrated iterative Bregman projection (IBP) algorithm, named \textsc{FastIBP}, with a complexity bound of $\tilde{O}(mn^{7/3}\varepsilon^{-4/3})$, where $\varepsilon \in (0, 1)$ is the desired tolerance. This complexity bound is better than the best known complexity bound of $\tilde{O}(mn^2\varepsilon^{-2})$ for the IBP algorithm in terms of $\varepsilon$, and that of $\tilde{O}(mn^{5/2}\varepsilon^{-1})$ from accelerated alternating minimization algorithm or accelerated primal-dual adaptive gradient algorithm in terms of $n$. Finally, we conduct extensive experiments with both synthetic data and real images and demonstrate the favorable performance of the \textsc{FastIBP} algorithm in practice.

preprint2022arXiv

FocalClick: Towards Practical Interactive Image Segmentation

Interactive segmentation allows users to extract target masks by making positive/negative clicks. Although explored by many previous works, there is still a gap between academic approaches and industrial needs: first, existing models are not efficient enough to work on low power devices; second, they perform poorly when used to refine preexisting masks as they could not avoid destroying the correct part. FocalClick solves both issues at once by predicting and updating the mask in localized areas. For higher efficiency, we decompose the slow prediction on the entire image into two fast inferences on small crops: a coarse segmentation on the Target Crop, and a local refinement on the Focus Crop. To make the model work with preexisting masks, we formulate a sub-task termed Interactive Mask Correction, and propose Progressive Merge as the solution. Progressive Merge exploits morphological information to decide where to preserve and where to update, enabling users to refine any preexisting mask effectively. FocalClick achieves competitive results against SOTA methods with significantly smaller FLOPs. It also shows significant superiority when making corrections on preexisting masks. Code and data will be released at github.com/XavierCHEN34/ClickSEG

preprint2022arXiv

Improved Upper Bounds for Finding Tarski Fixed Points

We study the query complexity of finding a Tarski fixed point over the $k$-dimensional grid $\{1,\ldots,n\}^k$. Improving on the previous best upper bound of $\smash{O(\log^{\lceil 2k/3\rceil} n)}$ [FPS20], we give a new algorithm with query complexity $\smash{O(\log^{\lceil (k+1)/2\rceil} n)}$. This is based on a novel decomposition theorem about a weaker variant of the Tarski fixed point problem, where the input consists of a monotone function $f:[n]^k\rightarrow [n]^k$ and a monotone sign function $b:[n]^k\rightarrow \{-1,0,1\}$ and the goal is to find an $x\in [n]^k$ that satisfies $either$ $f(x)\preceq x$ and $b(x)\le 0$ $or$ $f(x)\succeq x$ and $b(x)\ge 0$.

preprint2022arXiv

Investigation of Singing Voice Separation for Singing Voice Detection in Polyphonic Music

Singing voice detection (SVD), to recognize vocal parts in the song, is an essential task in music information retrieval (MIR). The task remains challenging since singing voice varies and intertwines with the accompaniment music, especially for some complicated polyphonic music such as choral music recordings. To address this problem, we investigate singing voice detection while discarding the interference from the accompaniment. The proposed SVD has two steps: i. The singing voice separation (SVS) technique is first utilized to filter out the singing voice's potential part coarsely. ii. Upon the continuity of vocal in the time domain, Long-term Recurrent Convolutional Networks (LRCN) is used to learn compositional features. Moreover, to eliminate the outliers, we choose to use a median filter for time-domain smoothing. Experimental results show that the proposed method outperforms the existing state-of-the-art works on two public datasets, the Jamendo Corpus and the RWC pop dataset.

preprint2022arXiv

Latent-Variable Advantage-Weighted Policy Optimization for Offline RL

Offline reinforcement learning methods hold the promise of learning policies from pre-collected datasets without the need to query the environment for new transitions. This setting is particularly well-suited for continuous control robotic applications for which online data collection based on trial-and-error is costly and potentially unsafe. In practice, offline datasets are often heterogeneous, i.e., collected in a variety of scenarios, such as data from several human demonstrators or from policies that act with different purposes. Unfortunately, such datasets can exacerbate the distribution shift between the behavior policy underlying the data and the optimal policy to be learned, leading to poor performance. To address this challenge, we propose to leverage latent-variable policies that can represent a broader class of policy distributions, leading to better adherence to the training data distribution while maximizing reward via a policy over the latent variable. As we empirically show on a range of simulated locomotion, navigation, and manipulation tasks, our method referred to as latent-variable advantage-weighted policy optimization (LAPO), improves the average performance of the next best-performing offline reinforcement learning methods by 49% on heterogeneous datasets, and by 8% on datasets with narrow and biased distributions.

preprint2022arXiv

Leveraging Intrinsic Gradient Information for Further Training of Differentiable Machine Learning Models

Designing models that produce accurate predictions is the fundamental objective of machine learning (ML). This work presents methods demonstrating that when the derivatives of target variables (outputs) with respect to inputs can be extracted from processes of interest, e.g., neural networks (NN) based surrogate models, they can be leveraged to further improve the accuracy of differentiable ML models. This paper generalises the idea and provides practical methodologies that can be used to leverage gradient information (GI) across a variety of applications including: (1) Improving the performance of generative adversarial networks (GANs); (2) efficiently tuning NN model complexity; (3) regularising linear regressions. Numerical results show that GI can effective enhance ML models with existing datasets, demonstrating its value for a variety of applications.

preprint2022arXiv

Machine-learning assisted quantum control in random environment

Disorder in condensed matter and atomic physics is responsible for a great variety of fascinating quantum phenomena, which are still challenging for understanding, not to mention the relevant dynamical control. Here we introduce proof of the concept and analyze neural network-based machine learning algorithm for achieving feasible high-fidelity quantum control of a particle in random environment. To explicitly demonstrate its capabilities, we show that convolutional neural networks are able to solve this problem as they can recognize the disorder and, by supervised learning, further produce the policy for the efficient low-energy cost control of a quantum particle in a time-dependent random potential. We have shown that the accuracy of the proposed algorithm is enhanced by a higher-dimensional mapping of the disorder pattern and using two neural networks, each properly trained for the given task. The designed method, being computationally more efficient than the gradient-descent optimization, can be applicable to identify and control various noisy quantum systems on a heuristic basis.

preprint2022arXiv

Memory Bounds for Continual Learning

Continual learning, or lifelong learning, is a formidable current challenge to machine learning. It requires the learner to solve a sequence of $k$ different learning tasks, one after the other, while retaining its aptitude for earlier tasks; the continual learner should scale better than the obvious solution of developing and maintaining a separate learner for each of the $k$ tasks. We embark on a complexity-theoretic study of continual learning in the PAC framework. We make novel uses of communication complexity to establish that any continual learner, even an improper one, needs memory that grows linearly with $k$, strongly suggesting that the problem is intractable. When logarithmically many passes over the learning tasks are allowed, we provide an algorithm based on multiplicative weights update whose memory requirement scales well; we also establish that improper learning is necessary for such performance. We conjecture that these results may lead to new promising approaches to continual learning.

preprint2022arXiv

Meta-Learning Digitized-Counterdiabatic Quantum Optimization

Solving optimization tasks using variational quantum algorithms has emerged as a crucial application of the current noisy intermediate-scale quantum devices. However, these algorithms face several difficulties like finding suitable ansatz and appropriate initial parameters, among others. In this work, we tackle the problem of finding suitable initial parameters for variational optimization by employing a meta-learning technique using recurrent neural networks. We investigate this technique with the recently proposed digitized-counterdiabatic quantum approximate optimization algorithm (DC-QAOA) that utilizes counterdiabatic protocols to improve the state-of-the-art QAOA. The combination of meta learning and DC-QAOA enables us to find optimal initial parameters for different models, such as MaxCut problem and the Sherrington-Kirkpatrick model. Decreasing the number of iterations of optimization as well as enhancing the performance, our protocol designs short depth circuit ansatz with optimal initial parameters by incorporating shortcuts-to-adiabaticity principles into machine learning methods for the near-term devices.

preprint2022arXiv

Mitigating Data Heterogeneity in Federated Learning with Data Augmentation

Federated Learning (FL) is a prominent framework that enables training a centralized model while securing user privacy by fusing local, decentralized models. In this setting, one major obstacle is data heterogeneity, i.e., each client having non-identically and independently distributed (non-IID) data. This is analogous to the context of Domain Generalization (DG), where each client can be treated as a different domain. However, while many approaches in DG tackle data heterogeneity from the algorithmic perspective, recent evidence suggests that data augmentation can induce equal or greater performance. Motivated by this connection, we present federated versions of popular DG algorithms, and show that by applying appropriate data augmentation, we can mitigate data heterogeneity in the federated setting, and obtain higher accuracy on unseen clients. Equipped with data augmentation, we can achieve state-of-the-art performance using even the most basic Federated Averaging algorithm, with much sparser communication.

preprint2022arXiv

Molecular conformer search with low-energy latent space

Identifying low-energy conformers with quantum mechanical accuracy for molecules with many degrees of freedom is challenging. In this work, we use the molecular dihedral angles as features and explore the possibility of performing molecular conformer search in a latent space with a generative model named variational auto-encoder (VAE). We bias the VAE towards low-energy molecular configurations to generate more informative data. In this way, we can effectively build a reliable energy model for the low-energy potential energy surface. After the energy model has been built, we extract local-minimum conformations and refine them with structure optimization. We have tested and benchmarked our low-energy latent-space (LOLS) structure search method on organic molecules with $5-9$ searching dimensions. Our results agree with previous studies.

preprint2022arXiv

Multi-modal learning for predicting the genotype of glioma

The isocitrate dehydrogenase (IDH) gene mutation is an essential biomarker for the diagnosis and prognosis of glioma. It is promising to better predict glioma genotype by integrating focal tumor image and geometric features with brain network features derived from MRI. Convolutions neural networks show reasonable performance in predicting IDH mutation, which, however, cannot learn from non-Euclidean data, e.g., geometric and network data. In this study, we propose a multi-modal learning framework using three separate encoders to extract features of focal tumor image, tumor geometrics and global brain networks. To mitigate the limited availability of diffusion MRI, we develop a self-supervised approach to generate brain networks from anatomical multi-sequence MRI. Moreover, to extract tumor-related features from the brain network, we design a hierarchical attention module for the brain network encoder. Further, we design a bi-level multi-modal contrastive loss to align the multi-modal features and tackle the domain gap at the focal tumor and global brain. Finally, we propose a weighted population graph to integrate the multi-modal features for genotype prediction. Experimental results on the testing set show that the proposed model outperforms the baseline deep learning models. The ablation experiments validate the performance of different components of the framework. The visualized interpretation corresponds to clinical knowledge with further validation. In conclusion, the proposed learning framework provides a novel approach for predicting the genotype of glioma.

preprint2022arXiv

Nested sampling for physical scientists

We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the application of NS in three scientific fields in which the algorithm has proved to be useful: cosmology, gravitational-wave astronomy, and materials science. We close by making recommendations for best practice when using NS and by summarizing potential limitations and optimizations of NS.

preprint2022arXiv

Neural 3D Reconstruction in the Wild

We are witnessing an explosion of neural implicit representations in computer vision and graphics. Their applicability has recently expanded beyond tasks such as shape generation and image-based rendering to the fundamental problem of image-based 3D reconstruction. However, existing methods typically assume constrained 3D environments with constant illumination captured by a small set of roughly uniformly distributed cameras. We introduce a new method that enables efficient and accurate surface reconstruction from Internet photo collections in the presence of varying illumination. To achieve this, we propose a hybrid voxel- and surface-guided sampling technique that allows for more efficient ray sampling around surfaces and leads to significant improvements in reconstruction quality. Further, we present a new benchmark and protocol for evaluating reconstruction performance on such in-the-wild scenes. We perform extensive experiments, demonstrating that our approach surpasses both classical and neural reconstruction methods on a wide variety of metrics.

preprint2022arXiv

No Weighted-Regret Learning in Adversarial Bandits with Delays

Consider a scenario where a player chooses an action in each round $t$ out of $T$ rounds and observes the incurred cost after a delay of $d_{t}$ rounds. The cost functions and the delay sequence are chosen by an adversary. We show that in a non-cooperative game, the expected weighted ergodic distribution of play converges to the set of coarse correlated equilibria if players use algorithms that have "no weighted-regret" in the above scenario, even if they have linear regret due to too large delays. For a two-player zero-sum game, we show that no weighted-regret is sufficient for the weighted ergodic average of play to converge to the set of Nash equilibria. We prove that the FKM algorithm with $n$ dimensions achieves an expected regret of $O\left(nT^{\frac{3}{4}}+\sqrt{n}T^{\frac{1}{3}}D^{\frac{1}{3}}\right)$ and the EXP3 algorithm with $K$ arms achieves an expected regret of $O\left(\sqrt{\log K\left(KT+D\right)}\right)$ even when $D=\sum_{t=1}^{T}d_{t}$ and $T$ are unknown. These bounds use a novel doubling trick that, under mild assumptions, provably retains the regret bound for when $D$ and $T$ are known. Using these bounds, we show that FKM and EXP3 have no weighted-regret even for $d_{t}=O\left(t\log t\right)$. Therefore, algorithms with no weighted-regret can be used to approximate a CCE of a finite or convex unknown game that can only be simulated with bandit feedback, even if the simulation involves significant delays.

preprint2022arXiv

On Digital Subcarrier Multiplexing under A Bandwidth Limitation and ASE Noise

We show that digital subcarrier multiplexing (DSM) systems require much greater complexity for Nyquist pulse shaping than single-carrier (SC) systems, and it is a misconception that both systems use the same bandwidth when using the same pulse shaping. Through back-to-back (B2B) experiments with realistic transmitter (TX) modules and amplified spontaneous emission (ASE) noise loading, we show that even with optimized waterfilling and entropy loading, DSM does not achieve a larger net data rate (NDR) compared to SC when only ASE noise exists in the channel in long-haul transmission scenarios.

preprint2022arXiv

On the Complexity of Dynamic Submodular Maximization

We study dynamic algorithms for the problem of maximizing a monotone submodular function over a stream of $n$ insertions and deletions. We show that any algorithm that maintains a $(0.5+ε)$-approximate solution under a cardinality constraint, for any constant $ε>0$, must have an amortized query complexity that is $\mathit{polynomial}$ in $n$. Moreover, a linear amortized query complexity is needed in order to maintain a $0.584$-approximate solution. This is in sharp contrast with recent dynamic algorithms of [LMNF+20, Mon20] that achieve $(0.5-ε)$-approximation with a $\mathsf{poly}\log(n)$ amortized query complexity. On the positive side, when the stream is insertion-only, we present efficient algorithms for the problem under a cardinality constraint and under a matroid constraint with approximation guarantee $1-1/e-ε$ and amortized query complexities $\smash{O(\log (k/ε)/ε^2)}$ and $\smash{k^{\tilde{O}(1/ε^2)}\log n}$, respectively, where $k$ denotes the cardinality parameter or the rank of the matroid.

preprint2022arXiv

Optimization of Directional Landmark Deployment for Visual Observer on SE(3)

An optimization method is proposed in this paper for novel deployment of given number of directional landmarks (location and pose) within a given region in the 3-D task space. This new deployment technique is built on the geometric models of both landmarks and the monocular camera. In particular, a new concept of Multiple Coverage Probability (MCP) is defined to characterize the probability of at least n landmarks being covered simultaneously by a camera at a fixed position. The optimization is conducted with respect to the position and pose of the given number of landmarks to maximize MCP through globally exploration of the given 3-D space. By adopting the elimination genetic algorithm, the global optimal solutions can be obtained, which are then applied to improve the convergent performance of the visual observer on SE(3) as a demonstration example. Both simulation and experimental results are presented to validate the effectiveness of the proposed landmark deployment optimization method.

preprint2022arXiv

PAC-Bayes Information Bottleneck

Understanding the source of the superior generalization ability of NNs remains one of the most important problems in ML research. There have been a series of theoretical works trying to derive non-vacuous bounds for NNs. Recently, the compression of information stored in weights (IIW) is proved to play a key role in NNs generalization based on the PAC-Bayes theorem. However, no solution of IIW has ever been provided, which builds a barrier for further investigation of the IIW's property and its potential in practical deep learning. In this paper, we propose an algorithm for the efficient approximation of IIW. Then, we build an IIW-based information bottleneck on the trade-off between accuracy and information complexity of NNs, namely PIB. From PIB, we can empirically identify the fitting to compressing phase transition during NNs' training and the concrete connection between the IIW compression and the generalization. Besides, we verify that IIW is able to explain NNs in broad cases, e.g., varying batch sizes, over-parameterization, and noisy labels. Moreover, we propose an MCMC-based algorithm to sample from the optimal weight posterior characterized by PIB, which fulfills the potential of IIW in enhancing NNs in practice.

preprint2022arXiv

PACTran: PAC-Bayesian Metrics for Estimating the Transferability of Pretrained Models to Classification Tasks

With the increasing abundance of pretrained models in recent years, the problem of selecting the best pretrained checkpoint for a particular downstream classification task has been gaining increased attention. Although several methods have recently been proposed to tackle the selection problem (e.g. LEEP, H-score), these methods resort to applying heuristics that are not well motivated by learning theory. In this paper we present PACTran, a theoretically grounded family of metrics for pretrained model selection and transferability measurement. We first show how to derive PACTran metrics from the optimal PAC-Bayesian bound under the transfer learning setting. We then empirically evaluate three metric instantiations of PACTran on a number of vision tasks (VTAB) as well as a language-and-vision (OKVQA) task. An analysis of the results shows PACTran is a more consistent and effective transferability measure compared to existing selection methods.

preprint2022arXiv

Radio properties of the OH megamaser galaxy IIZw 096

Based on the two epochs EVN archive data from OH line observations of IIZw 096, we confirm that the high-resolution OH emission in this source mainly comes from two spots (OH1 and OH2) of comp D1 of this merging system. We found no significant variations in the OH line emission. The OH 1665 MHz line emission is detected at about 6 $σ$ level in the OH1 region by combining two epoch EVN observations. We found that the comp D1 shows the brightest CO, HCO+ line emission, as well as multi-band radio continuum emission. The environment around D1 shows no clear velocity structure associated with circular motions, making it different from most other OHMs in the literature, which might have been caused by an effect during the merger stage. Meanwhile, we found that the CO emission shows three velocity structures around D1, including the central broad FWHM region, the double peak region where the CO line profile shows two separated peaks, and the region of the high-velocity clouds where the CO line peaks at a high velocity ($\sim$ 11000 \kms). \HI in absorption also show high-velocity clouds around the D1 region, which might be due to inflows caused by the merging of two or more galaxy components. Based on the high-resolution K-band VLA and L-band VLBA observations of the radio continuum emission, we derived the brightness temperature in the range $10^{5}$ K to $10^{6}$ K, which is consistent with other starburst dominant OHM sources in the literature. The multi-band VLA observations show that the radio continuum emission of comp D might also have contributions from free-free emission, besides synchrotron emission. As a concenquence, these results support a starburst origin for the OHMs, without the presence of an AGN.

preprint2022arXiv

Robust Dynamic Assortment Optimization in the Presence of Outlier Customers

We consider the dynamic assortment optimization problem under the multinomial logit model (MNL) with unknown utility parameters. The main question investigated in this paper is model mis-specification under the $\varepsilon$-contamination model, which is a fundamental model in robust statistics and machine learning. In particular, throughout a selling horizon of length $T$, we assume that customers make purchases according to a well specified underlying multinomial logit choice model in a $(1-\varepsilon)$-fraction of the time periods, and make arbitrary purchasing decisions instead in the remaining $\varepsilon$-fraction of the time periods. In this model, we develop a new robust online assortment optimization policy via an active elimination strategy. We establish both upper and lower bounds on the regret, and show that our policy is optimal up to logarithmic factor in $T$ when the assortment capacity is constant. %% capacity of assortments has a constant upper limit. We further develop a fully adaptive policy that does not require any prior knowledge of the contamination parameter $\varepsilon$. In the case of the existence a sub-optimality gap between optimal and sub-optimal products, we also established gap-dependent logarithmic regret upper bounds and lower bounds in both the known-$\varepsilon$ and unknown-$\varepsilon$ cases. Our simulation study shows that our policy outperforms the existing policies based on upper confidence bounds (UCB) and Thompson sampling.

preprint2022arXiv

Robust One Round Federated Learning with Predictive Space Bayesian Inference

Making predictions robust is an important challenge. A separate challenge in federated learning (FL) is to reduce the number of communication rounds, particularly since doing so reduces performance in heterogeneous data settings. To tackle both issues, we take a Bayesian perspective on the problem of learning a global model. We show how the global predictive posterior can be approximated using client predictive posteriors. This is unlike other works which aggregate the local model space posteriors into the global model space posterior, and are susceptible to high approximation errors due to the posterior's high dimensional multimodal nature. In contrast, our method performs the aggregation on the predictive posteriors, which are typically easier to approximate owing to the low-dimensionality of the output space. We present an algorithm based on this idea, which performs MCMC sampling at each client to obtain an estimate of the local posterior, and then aggregates these in one round to obtain a global ensemble model. Through empirical evaluation on several classification and regression tasks, we show that despite using one round of communication, the method is competitive with other FL techniques, and outperforms them on heterogeneous settings. The code is publicly available at https://github.com/hasanmohsin/FedPredSpace_1Round.

preprint2022arXiv

S2MS: Self-Supervised Learning Driven Multi-Spectral CT Image Enhancement

Photon counting spectral CT (PCCT) can produce reconstructed attenuation maps in different energy channels, reflecting energy properties of the scanned object. Due to the limited photon numbers and the non-ideal detector response of each energy channel, the reconstructed images usually contain much noise. With the development of Deep Learning (DL) technique, different kinds of DL-based models have been proposed for noise reduction. However, most of the models require clean data set as the training labels, which are not always available in medical imaging field. Inspiring by the similarities of each channel's reconstructed image, we proposed a self-supervised learning based PCCT image enhancement framework via multi-spectral channels (S2MS). In S2MS framework, both the input and output labels are noisy images. Specifically, one single channel image was used as output while images of other single channels and channel-sum image were used as input to train the network, which can fully use the spectral data information without extra cost. The simulation results based on the AAPM Low-dose CT Challenge database showed that the proposed S2MS model can suppress the noise and preserve details more effectively in comparison with the traditional DL models, which has potential to improve the image quality of PCCT in clinical applications.

preprint2022arXiv

Shortcuts to Adiabaticity for Fast Qubit Readout in Circuit Quantum Electrodynamics

We propose how to engineer the longitudinal coupling to accelerate the measurement of a qubit longitudinally coupled to a cavity, motivated by the concept of shortcuts to adiabaticity. Different modulations are inversely designed from two methods of inverse engineering and counter-diabatic driving, for achieving larger values of the signal-to-noise ratio (SNR) at nanosecond scale. By comparison, we demonstrate that our protocols outperform the usual periodic modulations on the pointer state separation and SNR. Finally, we show a possible implementation considering state-of-the-art circuit quantum electrodynamics architecture, estimating the minimal time allowed for the measurement process.

preprint2022arXiv

Two-stage Hypothesis Tests for Variable Interactions with FDR Control

In many scenarios such as genome-wide association studies where dependences between variables commonly exist, it is often of interest to infer the interaction effects in the model. However, testing pairwise interactions among millions of variables in complex and high-dimensional data suffers from low statistical power and huge computational cost. To address these challenges, we propose a two-stage testing procedure with false discovery rate (FDR) control, which is known as a less conservative multiple-testing correction. Theoretically, the difficulty in the FDR control dues to the data dependence among test statistics in two stages, and the fact that the number of hypothesis tests conducted in the second stage depends on the screening result in the first stage. By using the Cramér type moderate deviation technique, we show that our procedure controls FDR at the desired level asymptotically in the generalized linear model (GLM), where the model is allowed to be misspecified. In addition, the asymptotic power of the FDR control procedure is rigorously established. We demonstrate via comprehensive simulation studies that our two-stage procedure is computationally more efficient than the classical BH procedure, with a comparable or improved statistical power. Finally, we apply the proposed method to a bladder cancer data from dbGaP where the scientific goal is to identify genetic susceptibility loci for bladder cancer.

preprint2022arXiv

Unbiased Implicit Feedback via Bi-level Optimization

Implicit feedback is widely leveraged in recommender systems since it is easy to collect and provides weak supervision signals. Recent works reveal a huge gap between the implicit feedback and user-item relevance due to the fact that implicit feedback is also closely related to the item exposure. To bridge this gap, existing approaches explicitly model the exposure and propose unbiased estimators to improve the relevance. Unfortunately, these unbiased estimators suffer from the high gradient variance, especially for long-tail items, leading to inaccurate gradient updates and degraded model performance. To tackle this challenge, we propose a low-variance unbiased estimator from a probabilistic perspective, which effectively bounds the variance of the gradient. Unlike previous works which either estimate the exposure via heuristic-based strategies or use a large biased training set, we propose to estimate the exposure via an unbiased small-scale validation set. Specifically, we first parameterize the user-item exposure by incorporating both user and item information, and then construct an unbiased validation set from the biased training set. By leveraging the unbiased validation set, we adopt bi-level optimization to automatically update exposure-related parameters along with recommendation model parameters during the learning. Experiments on two real-world datasets and two semi-synthetic datasets verify the effectiveness of our method.

preprint2021arXiv

Adversarial Combinatorial Bandits with General Non-linear Reward Functions

In this paper we study the adversarial combinatorial bandit with a known non-linear reward function, extending existing work on adversarial linear combinatorial bandit. {The adversarial combinatorial bandit with general non-linear reward is an important open problem in bandit literature, and it is still unclear whether there is a significant gap from the case of linear reward, stochastic bandit, or semi-bandit feedback.} We show that, with $N$ arms and subsets of $K$ arms being chosen at each of $T$ time periods, the minimax optimal regret is $\widetildeΘ_{d}(\sqrt{N^d T})$ if the reward function is a $d$-degree polynomial with $d< K$, and $Θ_K(\sqrt{N^K T})$ if the reward function is not a low-degree polynomial. {Both bounds are significantly different from the bound $O(\sqrt{\mathrm{poly}(N,K)T})$ for the linear case, which suggests that there is a fundamental gap between the linear and non-linear reward structures.} Our result also finds applications to adversarial assortment optimization problem in online recommendation. We show that in the worst-case of adversarial assortment problem, the optimal algorithm must treat each individual $\binom{N}{K}$ assortment as independent.

preprint2021arXiv

Bayesian Meta-Learning for Few-Shot Policy Adaptation Across Robotic Platforms

Reinforcement learning methods can achieve significant performance but require a large amount of training data collected on the same robotic platform. A policy trained with expensive data is rendered useless after making even a minor change to the robot hardware. In this paper, we address the challenging problem of adapting a policy, trained to perform a task, to a novel robotic hardware platform given only few demonstrations of robot motion trajectories on the target robot. We formulate it as a few-shot meta-learning problem where the goal is to find a meta-model that captures the common structure shared across different robotic platforms such that data-efficient adaptation can be performed. We achieve such adaptation by introducing a learning framework consisting of a probabilistic gradient-based meta-learning algorithm that models the uncertainty arising from the few-shot setting with a low-dimensional latent variable. We experimentally evaluate our framework on a simulated reaching and a real-robot picking task using 400 simulated robots generated by varying the physical parameters of an existing set of robotic platforms. Our results show that the proposed method can successfully adapt a trained policy to different robotic platforms with novel physical parameters and the superiority of our meta-learning algorithm compared to state-of-the-art methods for the introduced few-shot policy adaptation problem.

preprint2021arXiv

Confronting the Carbon-footprint Challenge of Blockchain

The distributed consensus mechanism is the backbone of the rapidly developing blockchain network. Blockchain platforms consume vast amounts of electricity based on the current consensus mechanism of Proof of Work. Here, we point out an advanced consensus mechanism named Proof of Stake that can eliminate the extensive energy consumption of the current PoW-based blockchain. We comprehensively elucidate the current and projected energy consumption and carbon footprint of the PoW and PoS based Bitcoin and Ethereum blockchain platforms.

preprint2021arXiv

Connecting The Dots To Combat Collective Fraud

Modern fraudsters write malicious programs to coordinate a group of accounts to commit collective fraud for illegal profits in online platforms. These programs have access to a set of finite resources - a set of IPs, devices, and accounts etc. and sometime manipulate fake accounts to collaboratively attack the target system. Inspired by these observations, we share our experience in building two real-time risk control systems to detect collective fraud. We show that with TigerGraph, a powerful graph database, and its innovative query language - GSQL, data scientists and fraud experts can conveniently implement and deploy an end-to-end risk control system as a graph database application.

preprint2021arXiv

Connection between inverse engineering and optimal control in shortcuts to adiabaticity

We consider fast high-fidelity quantum control by using a shortcut to adiabaticity (STA) technique and optimal control theory (OCT). Three specific examples, including expansion of cold atoms from the harmonic trap, atomic transport by moving harmonic trap, and spin dynamics in the presence of dissipation, are explicitly detailed. Using OCT as a qualitative guide, we demonstrate how STA protocols designed from inverse engineering method, can approach with very high precision optimal solutions built about physical constraints, by a proper choice of the interpolation function and with a very reduced number of adjustable parameters.

preprint2021arXiv

Dimension Independent Generalization Error by Stochastic Gradient Descent

One classical canon of statistics is that large models are prone to overfitting, and model selection procedures are necessary for high dimensional data. However, many overparameterized models, such as neural networks, perform very well in practice, although they are often trained with simple online methods and regularization. The empirical success of overparameterized models, which is often known as benign overfitting, motivates us to have a new look at the statistical generalization theory for online optimization. In particular, we present a general theory on the generalization error of stochastic gradient descent (SGD) solutions for both convex and locally convex loss functions. We further discuss data and model conditions that lead to a ``low effective dimension". Under these conditions, we show that the generalization error either does not depend on the ambient dimension $p$ or depends on $p$ via a poly-logarithmic factor. We also demonstrate that in several widely used statistical models, the ``low effective dimension'' arises naturally in overparameterized settings. The studied statistical applications include both convex models such as linear regression and logistic regression and non-convex models such as $M$-estimator and two-layer neural networks.

preprint2021arXiv

Distributed Estimation for Principal Component Analysis: an Enlarged Eigenspace Analysis

The growing size of modern data sets brings many challenges to the existing statistical estimation approaches, which calls for new distributed methodologies. This paper studies distributed estimation for a fundamental statistical machine learning problem, principal component analysis (PCA). Despite the massive literature on top eigenvector estimation, much less is presented for the top-$L$-dim ($L>1$) eigenspace estimation, especially in a distributed manner. We propose a novel multi-round algorithm for constructing top-$L$-dim eigenspace for distributed data. Our algorithm takes advantage of shift-and-invert preconditioning and convex optimization. Our estimator is communication-efficient and achieves a fast convergence rate. In contrast to the existing divide-and-conquer algorithm, our approach has no restriction on the number of machines. Theoretically, the traditional Davis-Kahan theorem requires the explicit eigengap assumption to estimate the top-$L$-dim eigenspace. To abandon this eigengap assumption, we consider a new route in our analysis: instead of exactly identifying the top-$L$-dim eigenspace, we show that our estimator is able to cover the targeted top-$L$-dim population eigenspace. Our distributed algorithm can be applied to a wide range of statistical problems based on PCA, such as principal component regression and single index model. Finally, We provide simulation studies to demonstrate the performance of the proposed distributed estimator.

preprint2021arXiv

DoubleEnsemble: A New Ensemble Method Based on Sample Reweighting and Feature Selection for Financial Data Analysis

Modern machine learning models (such as deep neural networks and boosting decision tree models) have become increasingly popular in financial market prediction, due to their superior capacity to extract complex non-linear patterns. However, since financial datasets have very low signal-to-noise ratio and are non-stationary, complex models are often very prone to overfitting and suffer from instability issues. Moreover, as various machine learning and data mining tools become more widely used in quantitative trading, many trading firms have been producing an increasing number of features (aka factors). Therefore, how to automatically select effective features becomes an imminent problem. To address these issues, we propose DoubleEnsemble, an ensemble framework leveraging learning trajectory based sample reweighting and shuffling based feature selection. Specifically, we identify the key samples based on the training dynamics on each sample and elicit key features based on the ablation impact of each feature via shuffling. Our model is applicable to a wide range of base models, capable of extracting complex patterns, while mitigating the overfitting and instability issues for financial market prediction. We conduct extensive experiments, including price prediction for cryptocurrencies and stock trading, using both DNN and gradient boosting decision tree as base models. Our experiment results demonstrate that DoubleEnsemble achieves a superior performance compared with several baseline methods.

preprint2021arXiv

Dynamic Assortment Selection under the Nested Logit Models

We study a stylized dynamic assortment planning problem during a selling season of finite length $T$. At each time period, the seller offers an arriving customer an assortment of substitutable products and the customer makes the purchase among offered products according to a discrete choice model. The goal of the seller is to maximize the expected revenue, or equivalently, to minimize the worst-case expected regret. One key challenge is that utilities of products are unknown to the seller and need to be learned. Although the dynamic assortment planning problem has received increasing attention in revenue management, most existing work is based on the multinomial logit choice models (MNL). In this paper, we study the problem of dynamic assortment planning under a more general choice model -- the nested logit model, which models hierarchical choice behavior and is ``the most widely used member of the GEV (generalized extreme value) family''. By leveraging the revenue-ordered structure of the optimal assortment within each nest, we develop a novel upper confidence bound (UCB) policy with an aggregated estimation scheme. Our policy simultaneously learns customers' choice behavior and makes dynamic decisions on assortments based on the current knowledge. It achieves the accumulated regret at the order of $\tilde{O}(\sqrt{MNT})$, where $M$ is the number of nests and $N$ is the number of products in each nest. We further provide a lower bound result of $Ω(\sqrt{MT})$, which shows the near optimality of the upper bound when $T$ is much larger than $M$ and $N$. When the number of items per nest $N$ is large, we further provide a discretization heuristic for better performance of our algorithm. Numerical results are presented to demonstrate the empirical performance of our proposed algorithms.

preprint2021arXiv

Echo state graph neural networks with analogue random resistor arrays

Recent years have witnessed an unprecedented surge of interest, from social networks to drug discovery, in learning representations of graph-structured data. However, graph neural networks, the machine learning models for handling graph-structured data, face significant challenges when running on conventional digital hardware, including von Neumann bottleneck incurred by physically separated memory and processing units, slowdown of Moore's law due to transistor scaling limit, and expensive training cost. Here we present a novel hardware-software co-design, the random resistor array-based echo state graph neural network, which addresses these challenges. The random resistor arrays not only harness low-cost, nanoscale and stackable resistors for highly efficient in-memory computing using simple physical laws, but also leverage the intrinsic stochasticity of dielectric breakdown to implement random projections in hardware for an echo state network that effectively minimizes the training cost thanks to its fixed and random weights. The system demonstrates state-of-the-art performance on both graph classification using the MUTAG and COLLAB datasets and node classification using the CORA dataset, achieving 34.2x, 93.2x, and 570.4x improvement of energy efficiency and 98.27%, 99.46%, and 95.12% reduction of training cost compared to conventional graph learning on digital hardware, respectively, which may pave the way for the next generation AI system for graph learning.

preprint2021arXiv

Enquire One's Parent and Child Before Decision: Fully Exploit Hierarchical Structure for Self-Supervised Taxonomy Expansion

Taxonomy is a hierarchically structured knowledge graph that plays a crucial role in machine intelligence. The taxonomy expansion task aims to find a position for a new term in an existing taxonomy to capture the emerging knowledge in the world and keep the taxonomy dynamically updated. Previous taxonomy expansion solutions neglect valuable information brought by the hierarchical structure and evaluate the correctness of merely an added edge, which downgrade the problem to node-pair scoring or mini-path classification. In this paper, we propose the Hierarchy Expansion Framework (HEF), which fully exploits the hierarchical structure's properties to maximize the coherence of expanded taxonomy. HEF makes use of taxonomy's hierarchical structure in multiple aspects: i) HEF utilizes subtrees containing most relevant nodes as self-supervision data for a complete comparison of parental and sibling relations; ii) HEF adopts a coherence modeling module to evaluate the coherence of a taxonomy's subtree by integrating hypernymy relation detection and several tree-exclusive features; iii) HEF introduces the Fitting Score for position selection, which explicitly evaluates both path and level selections and takes full advantage of parental relations to interchange information for disambiguation and self-correction. Extensive experiments show that by better exploiting the hierarchical structure and optimizing taxonomy's coherence, HEF vastly surpasses the prior state-of-the-art on three benchmark datasets by an average improvement of 46.7% in accuracy and 32.3% in mean reciprocal rank.

preprint2021arXiv

Experimental quantum simulation of superradiant phase transition beyond no-go theorem via antisqueezing

Superradiant phase transition (SPT) in thermal equilibrium, as a fundamental concept bridging the statistical physics and electrodynamics, can offer the key resources for quantum information science. Notwithstanding its fundamental and practical significances, equilibrium SPT has never been observed in experiments since the first proposal in the 1970s. Furthermore, the existence of equilibrium SPT in the cavity quantum electrodynamics (QED) systems is still subject of ongoing debates, due to the no-go theorem induced by the so-called A2 term. Based on the platform of nuclear magnetic resonance (NMR), here we experimentally demonstrate the occurrence of equilibrium SPT beyond no-go theorem by introducing the antisqueezing effect. The mechanism relies on the antisqueezing that recovers the singularity of the ground state via exponentially enhancing the zero point fluctuation (ZPF) of system. The strong entanglement and squeezed Schrodinger cat states of spins are achieved experimentally in the superradiant phase, which may play an important role in fundamental tests of quantum theory, implementing quantum metrology and high-efficient quantum information processing. Our experiment also shows the antisqueezing-enhanced signal-to-noise rate (SNR) of NMR spectrum, providing a general method for implementing high-precision NMR experiments.

preprint2021arXiv

Experimentally Realizing Efficient Quantum Control with Reinforcement Learning

Robust and high-precision quantum control is crucial but challenging for scalable quantum computation and quantum information processing. Traditional adiabatic control suffers severe limitations on gate performance imposed by environmentally induced noise because of a quantum system's limited coherence time. In this work, we experimentally demonstrate an alternative approach {to quantum control} based on deep reinforcement learning (DRL) on a trapped $^{171}\mathrm{Yb}^{+}$ ion. In particular, we find that DRL leads to fast and robust {digital quantum operations with running time bounded by shortcuts to adiabaticity} (STA). Besides, we demonstrate that DRL's robustness against both Rabi and detuning errors can be achieved simultaneously without any input from STA. Our experiments reveal a general framework of digital quantum control, leading to a promising enhancement in quantum information processing.

preprint2021arXiv

Fast-forward scaling of atom-molecule conversion in Bose-Einstein condensates

Robust stimulated Raman exact passages are requisite for controlling nonlinear quantum systems, with the wide applications ranging from ultracold molecules, non-linear optics to superchemistry. Inspired by shortcuts to adiabaticity, we propose the fast-forward scaling of stimulated Raman adiabatic processes with the nonlinearity involved, describing the transfer from an atomic Bose-Einstein condensate to a molecular one by controllable external fields. The fidelity and robustness of atom-molecule conversion are shown to surpass those of conventional adiabatic passages, assisted by fast-forward driving field. Finally, our results are extended to the fractional stimulated Raman adiabatic processes for the coherent superposition of atomic and molecular states.

preprint2021arXiv

First-order Newton-type Estimator for Distributed Estimation and Inference

This paper studies distributed estimation and inference for a general statistical problem with a convex loss that could be non-differentiable. For the purpose of efficient computation, we restrict ourselves to stochastic first-order optimization, which enjoys low per-iteration complexity. To motivate the proposed method, we first investigate the theoretical properties of a straightforward Divide-and-Conquer Stochastic Gradient Descent (DC-SGD) approach. Our theory shows that there is a restriction on the number of machines and this restriction becomes more stringent when the dimension $p$ is large. To overcome this limitation, this paper proposes a new multi-round distributed estimation procedure that approximates the Newton step only using stochastic subgradient. The key component in our method is the proposal of a computationally efficient estimator of $Σ^{-1} w$, where $Σ$ is the population Hessian matrix and $w$ is any given vector. Instead of estimating $Σ$ (or $Σ^{-1}$) that usually requires the second-order differentiability of the loss, the proposed First-Order Newton-type Estimator (FONE) directly estimates the vector of interest $Σ^{-1} w$ as a whole and is applicable to non-differentiable losses. Our estimator also facilitates the inference for the empirical risk minimizer. It turns out that the key term in the limiting covariance has the form of $Σ^{-1} w$, which can be estimated by FONE.

preprint2021arXiv

Hybrid quantum-classical approach to enhanced quantum metrology

Quantum metrology plays a fundamental role in many scientific areas. However, the complexity of engineering entangled probes and the external noise raise technological barriers for realizing the expected precision of the to-be-estimated parameter with given resources. Here, we address this problem by introducing adjustable controls into the encoding process and then utilizing a hybrid quantum-classical approach to automatically optimize the controls online. Our scheme does not require any complex or intractable off-line design, and it can inherently correct certain unitary errors during the learning procedure. We also report the first experimental demonstration of this promising scheme for the task of finding optimal probes for frequency estimation on a nuclear magnetic resonance (NMR) processor. The proposed scheme paves the way to experimentally auto-search optimal protocol for improving the metrology precision.

preprint2021arXiv

Implementation of a Hybrid Classical-Quantum Annealing Algorithm for Logistic Network Design

The logistic network design is an abstract optimization problem that, under the assumption of minimal cost, seeks the optimal configuration of the supply chain's infrastructures and facilities based on customer demand. Key economic decisions are taken about the location, number, and size of manufacturing facilities and warehouses based on the optimal solution. Therefore, improvements in the methods to address this question, which is known to be in the NP-hard complexity class, would have relevant financial consequences. Here, we implement in the D-Wave quantum annealer a hybrid classical-quantum annealing algorithm. The cost function with constraints is translated to a spin Hamiltonian, whose ground state encodes the searched result. As a benchmark, we measure the accuracy of results for a set of paradigmatic problems against the optimal published solutions (the error is on average below $1\%$), and the performance is compared against the classical algorithm, showing a remarkable reduction in the number of iterations. This work shows that state-of-the-art quantum annealers may codify and solve relevant supply-chain problems even still far from useful quantum supremacy.

preprint2021arXiv

Learning a Product Relevance Model from Click-Through Data in E-Commerce

The search engine plays a fundamental role in online e-commerce systems, to help users find the products they want from the massive product collections. Relevance is an essential requirement for e-commerce search, since showing products that do not match search query intent will degrade user experience. With the existence of vocabulary gap between user language of queries and seller language of products, measuring semantic relevance is necessary and neural networks are engaged to address this task. However, semantic relevance is different from click-through rate prediction in that no direct training signal is available. Most previous attempts learn relevance models from user click-through data that are cheap and abundant. Unfortunately, click behavior is noisy and misleading, which is affected by not only relevance but also factors including price, image and attractive titles. Therefore, it is challenging but valuable to learn relevance models from click-through data. In this paper, we propose a new relevance learning framework that concentrates on how to train a relevance model from the weak supervision of click-through data. Different from previous efforts that treat samples as either relevant or irrelevant, we construct more fine-grained samples for training. We propose a novel way to consider samples of different relevance confidence, and come up with a new training objective to learn a robust relevance model with desirable score distribution. The proposed model is evaluated on offline annotated data and online A/B testing, and it achieves both promising performance and high computational efficiency. The model has already been deployed online, serving the search traffic of Taobao for over a year.

preprint2021arXiv

Lifelong Learning based Disease Diagnosis on Clinical Notes

Current deep learning based disease diagnosis systems usually fall short in catastrophic forgetting, i.e., directly fine-tuning the disease diagnosis model on new tasks usually leads to abrupt decay of performance on previous tasks. What is worse, the trained diagnosis system would be fixed once deployed but collecting training data that covers enough diseases is infeasible, which inspires us to develop a lifelong learning diagnosis system. In this work, we propose to adopt attention to combine medical entities and context, embedding episodic memory and consolidation to retain knowledge, such that the learned model is capable of adapting to sequential disease-diagnosis tasks. Moreover, we establish a new benchmark, named Jarvis-40, which contains clinical notes collected from various hospitals. Our experiments show that the proposed method can achieve state-of-the-art performance on the proposed benchmark.

preprint2021arXiv

Online Disease Self-diagnosis with Inductive Heterogeneous Graph Convolutional Networks

We propose a Healthcare Graph Convolutional Network (HealGCN) to offer disease self-diagnosis service for online users based on Electronic Healthcare Records (EHRs). Two main challenges are focused in this paper for online disease diagnosis: (1) serving cold-start users via graph convolutional networks and (2) handling scarce clinical description via a symptom retrieval system. To this end, we first organize the EHR data into a heterogeneous graph that is capable of modeling complex interactions among users, symptoms and diseases, and tailor the graph representation learning towards disease diagnosis with an inductive learning paradigm. Then, we build a disease self-diagnosis system with a corresponding EHR Graph-based Symptom Retrieval System (GraphRet) that can search and provide a list of relevant alternative symptoms by tracing the predefined meta-paths. GraphRet helps enrich the seed symptom set through the EHR graph when confronting users with scarce descriptions, hence yield better diagnosis accuracy. At last, we validate the superiority of our model on a large-scale EHR dataset.

preprint2021arXiv

Random Restrictions of High-Dimensional Distributions and Uniformity Testing with Subcube Conditioning

We give a nearly-optimal algorithm for testing uniformity of distributions supported on $\{-1,1\}^n$, which makes $\tilde O (\sqrt{n}/\varepsilon^2)$ queries to a subcube conditional sampling oracle (Bhattacharyya and Chakraborty (2018)). The key technical component is a natural notion of random restriction for distributions on $\{-1,1\}^n$, and a quantitative analysis of how such a restriction affects the mean vector of the distribution. Along the way, we consider the problem of mean testing with independent samples and provide a nearly-optimal algorithm.

preprint2021arXiv

Shape-Enforcing Operators for Point and Interval Estimators

A common problem in econometrics, statistics, and machine learning is to estimate and make inference on functions that satisfy shape restrictions. For example, distribution functions are nondecreasing and range between zero and one, height growth charts are nondecreasing in age, and production functions are nondecreasing and quasi-concave in input quantities. We propose a method to enforce these restrictions ex post on point and interval estimates of the target function by applying functional operators. If an operator satisfies certain properties that we make precise, the shape-enforced point estimates are closer to the target function than the original point estimates and the shape-enforced interval estimates have greater coverage and shorter length than the original interval estimates. We show that these properties hold for six different operators that cover commonly used shape restrictions in practice: range, convexity, monotonicity, monotone convexity, quasi-convexity, and monotone quasi-convexity. We illustrate the results with two empirical applications to the estimation of a height growth chart for infants in India and a production function for chemical firms in China.

preprint2021arXiv

Some continuous and discontinuous Galerkin methods and structure preservation for incompressible flows

In this paper, we present consistent and inconsistent discontinuous Galerkin methods for incompressible Euler and Navier-Stokes equations with the kinematic pressure, Bernoulli function and EMAC function. Semi- and fully discrete energy stability of the proposed dG methods are proved in a unified fashion. Conservation of total energy, linear and angular momentum is discussed with both central and upwind fluxes. Numerical experiments are presented to demonstrate our findings and compare our schemes with conventional schemes in the literature in both unsteady and steady problems. Numerical results show that global conservation of the physical quantities may not be enough to demonstrate the performance of the schemes, and our schemes are competitive and able to capture essential physical features in several benchmark problems.

preprint2021arXiv

Tight Regret Bounds for Infinite-armed Linear Contextual Bandits

Linear contextual bandit is an important class of sequential decision making problems with a wide range of applications to recommender systems, online advertising, healthcare, and many other machine learning related tasks. While there is a lot of prior research, tight regret bounds of linear contextual bandit with infinite action sets remain open. In this paper, we address this open problem by considering the linear contextual bandit with (changing) infinite action sets. We prove a regret upper bound on the order of $O(\sqrt{d^2T\log T})\times \text{poly}(\log\log T)$ where $d$ is the domain dimension and $T$ is the time horizon. Our upper bound matches the previous lower bound of $Ω(\sqrt{d^2 T\log T})$ in [Li et al., 2019] up to iterated logarithmic terms.

preprint2021arXiv

Variance Reduced Median-of-Means Estimator for Byzantine-Robust Distributed Inference

This paper develops an efficient distributed inference algorithm, which is robust against a moderate fraction of Byzantine nodes, namely arbitrary and possibly adversarial machines in a distributed learning system. In robust statistics, the median-of-means (MOM) has been a popular approach to hedge against Byzantine failures due to its ease of implementation and computational efficiency. However, the MOM estimator has the shortcoming in terms of statistical efficiency. The first main contribution of the paper is to propose a variance reduced median-of-means (VRMOM) estimator, which improves the statistical efficiency over the vanilla MOM estimator and is computationally as efficient as the MOM. Based on the proposed VRMOM estimator, we develop a general distributed inference algorithm that is robust against Byzantine failures. Theoretically, our distributed algorithm achieves a fast convergence rate with only a constant number of rounds of communications. We also provide the asymptotic normality result for the purpose of statistical inference. To the best of our knowledge, this is the first normality result in the setting of Byzantine-robust distributed learning. The simulation results are also presented to illustrate the effectiveness of our method.

preprint2020arXiv

44 GHz methanol masers: Observations toward 95 GHz methanol masers

We report a simultaneous 44 and 95 GHz class I methanol maser survey toward 144 sources from the 95 GHz class I methanol maser catalog. The observations were made with the three telescopes of the Korean very long baseline interferometry network operating in single-dish mode. The detection rates are 89% at 44 GHz and 77% at 95 GHz. There are 106 new discoveries at 44 GHz. Comparing the previous 95 GHz detections with new observations of the same transitions made using the Purple Mountain Observatory 13.7 m radio telescope shows no clear evidence of variability on a timescale of six years. Emission from the 44 and 95 GHz transitions shows strong correlations in peak velocity, peak flux density, and integrated flux density, indicating that they are likely cospatial. We found that the peak flux density ratio Spk,95/Spk,44 decreases as the 44 GHz peak flux density increases. We found that some class I methanol masers in our sample could be associated with infrared dark clouds, while others are associated with H II regions, indicating that some sources occur at an early stage of high-mass star formation, while others are located toward more evolved sources.

preprint2020arXiv

A 4-6 GHz Radio Recombination Line Survey in the Milky Way

We performed a radio recombination line (RRL) survey to construct a high-mass star-forming region (HMSFR) sample in the Milky Way based on the all-sky Wide-Field Infrared Survey Explorer ($\textit{All-WISE}$) point source catalog. The survey was observed with the Shanghai 65m Tianma radio telescope (TMRT) covering 10 hydrogen RRL transitions ranging from H98$α$ to H113$α$ (corresponding to the rest frequencies of 4.5$-$6.9 GHz) simultaneously. Out of 3348 selected targets, we identified an HMSFR sample consisting of 517 sources traced by RRLs, a large fraction of this sample (486) locate near the Galactic plane ($|$$\textit{b}$$|$ $<$ 2 deg). In addition to the hydrogen RRLs, we also detected helium and carbon RRLs towards 49 and 23 sources respectively. We cross-match the RRL detections with the 6.7 methanol maser sources built up in previous works for the same target sample, as a result, 103 HMSFR sources were found to harbor both emissions. In this paper, we present the HMSFR catalog accompanied by the measured RRL line properties and a correlation with our methanol maser sample, which is believed to tracer massive stars at earlier stages. The construction of an HMSFR sample consisting of sources in various evolutionary stages indicated by different tracers is fundamental for future studies of high-mass star formation in such regions.

preprint2020arXiv

A picture of pseudogap phase related to charge fluxes

Recently, charge density fluctuations or charge fluxes attract strong interests in understanding the unconventional superconductivity. In this paper, a new emergent configuration in cuprates is identified by density functional theory simulations, called the charge pseudoplane, which exhibits the property of confining the dynamic charge fluxes for higher superconducting transition temperatures. It further redefines the fundamental collective excitation in cuprates as pQon with the momentum-dependent and ultrafast localization-delocalization duality. It is shown that both pseudogap and superconducting phases can be born from and intertwined through the charge flux confinement property of the charge pseudoplane region. Our experimental simulations based on the new picture provide good agreements with previous angle resolved photoemission spectroscopy and scanning tunneling microscopy results. Our work thus opens a new perspective on the origin of pseudogap phase and other related phases in cuprates, and further provides a critical descriptor to search and design higher temperature superconductors.

preprint2020arXiv

Adaptive Discrete Smoothing for High-Dimensional and Nonlinear Panel Data

In this paper we develop a data-driven smoothing technique for high-dimensional and non-linear panel data models. We allow for individual specific (non-linear) functions and estimation with econometric or machine learning methods by using weighted observations from other individuals. The weights are determined by a data-driven way and depend on the similarity between the corresponding functions and are measured based on initial estimates. The key feature of such a procedure is that it clusters individuals based on the distance / similarity between them, estimated in a first stage. Our estimation method can be combined with various statistical estimation procedures, in particular modern machine learning methods which are in particular fruitful in the high-dimensional case and with complex, heterogeneous data. The approach can be interpreted as a \textquotedblleft soft-clustering\textquotedblright\ in comparison to traditional\textquotedblleft\ hard clustering\textquotedblright that assigns each individual to exactly one group. We conduct a simulation study which shows that the prediction can be greatly improved by using our estimator. Finally, we analyze a big data set from didichuxing.com, a leading company in transportation industry, to analyze and predict the gap between supply and demand based on a large set of covariates. Our estimator clearly performs much better in out-of-sample prediction compared to existing linear panel data estimators.

preprint2020arXiv

ALPINE: Active Link Prediction using Network Embedding

Many real-world problems can be formalized as predicting links in a partially observed network. Examples include Facebook friendship suggestions, consumer-product recommendations, and the identification of hidden interactions between actors in a crime network. Several link prediction algorithms, notably those recently introduced using network embedding, are capable of doing this by just relying on the observed part of the network. Often, the link status of a node pair can be queried, which can be used as additional information by the link prediction algorithm. Unfortunately, such queries can be expensive or time-consuming, mandating the careful consideration of which node pairs to query. In this paper we estimate the improvement in link prediction accuracy after querying any particular node pair, to use in an active learning setup. Specifically, we propose ALPINE (Active Link Prediction usIng Network Embedding), the first method to achieve this for link prediction based on network embedding. To this end, we generalized the notion of V-optimality from experimental design to this setting, as well as more basic active learning heuristics originally developed in standard classification settings. Empirical results on real data show that ALPINE is scalable, and boosts link prediction accuracy with far fewer queries.

preprint2020arXiv

Bright solitons in a spin-tensor-momentum-coupled Bose-Einstein condensate

Synthetic spin-tensor-momentum coupling has recently been proposed to realize in atomic Bose-Einstein condensates. Here we study bright solitons in Bose-Einstein condensates with spin-tensor-momentum coupling and spin-orbit coupling. The properties and dynamics of spin-tensor-momentum-coupled and spin-orbit-coupled bright solitons are identified to be different. We contribute the difference to the different symmetries.

preprint2020arXiv

Composition and Configuration Patterns in Multiple-View Visualizations

Multiple-view visualization (MV) is a layout design technique often employed to help users see a large number of data attributes and values in a single cohesive representation. Because of its generalizability, the MV design has been widely adopted by the visualization community to help users examine and interact with large, complex, and high-dimensional data. However, although ubiquitous, there has been little work to categorize and analyze MVs in order to better understand its design space. As a result, there has been little to no guideline in how to use the MV design effectively. In this paper, we present an in-depth study of how MVs are designed in practice. We focus on two fundamental measures of multiple-view patterns: composition, which quantifies what view types and how many are there; and configuration, which characterizes spatial arrangement of view layouts in the display space. We build a new dataset containing 360 images of MVs collected from IEEE VIS, EuroVis, and PacificVis publications 2011 to 2019, and make fine-grained annotations of view types and layouts for these visualization images. From this data we conduct composition and configuration analyses using quantitative metrics of term frequency and layout topology. We identify common practices around MVs, including relationship of view types, popular view layouts, and correlation between view types and layouts. We combine the findings into a MV recommendation system, providing interactive tools to explore the design space, and support example-based design.

preprint2020arXiv

Correcting Knowledge Base Assertions

The usefulness and usability of knowledge bases (KBs) is often limited by quality issues. One common issue is the presence of erroneous assertions, often caused by lexical or semantic confusion. We study the problem of correcting such assertions, and present a general correction framework which combines lexical matching, semantic embedding, soft constraint mining and semantic consistency checking. The framework is evaluated using DBpedia and an enterprise medical KB.

preprint2020arXiv

COVID-19 Public Opinion and Emotion Monitoring System Based on Time Series Thermal New Word Mining

With the spread and development of new epidemics, it is of great reference value to identify the changing trends of epidemics in public emotions. We designed and implemented the COVID-19 public opinion monitoring system based on time series thermal new word mining. A new word structure discovery scheme based on the timing explosion of network topics and a Chinese sentiment analysis method for the COVID-19 public opinion environment is proposed. Establish a "Scrapy-Redis-Bloomfilter" distributed crawler framework to collect data. The system can judge the positive and negative emotions of the reviewer based on the comments, and can also reflect the depth of the seven emotions such as Hopeful, Happy, and Depressed. Finally, we improved the sentiment discriminant model of this system and compared the sentiment discriminant error of COVID-19 related comments with the Jiagu deep learning model. The results show that our model has better generalization ability and smaller discriminant error. We designed a large data visualization screen, which can clearly show the trend of public emotions, the proportion of various emotion categories, keywords, hot topics, etc., and fully and intuitively reflect the development of public opinion.

preprint2020arXiv

DAN: A Deformation-Aware Network for Consecutive Biomedical Image Interpolation

The continuity of biological tissue between consecutive biomedical images makes it possible for the video interpolation algorithm, to recover large area defects and tears that are common in biomedical images. However, noise and blur differences, large deformation, and drift between biomedical images, make the task challenging. To address the problem, this paper introduces a deformation-aware network to synthesize each pixel in accordance with the continuity of biological tissue. First, we develop a deformation-aware layer for consecutive biomedical images interpolation that implicitly adopting global perceptual deformation. Second, we present an adaptive style-balance loss to take the style differences of consecutive biomedical images such as blur and noise into consideration. Guided by the deformation-aware module, we synthesize each pixel from a global domain adaptively which further improves the performance of pixel synthesis. Quantitative and qualitative experiments on the benchmark dataset show that the proposed method is superior to the state-of-the-art approaches.

preprint2020arXiv

Distributed High-dimensional Regression Under a Quantile Loss Function

This paper studies distributed estimation and support recovery for high-dimensional linear regression model with heavy-tailed noise. To deal with heavy-tailed noise whose variance can be infinite, we adopt the quantile regression loss function instead of the commonly used squared loss. However, the non-smooth quantile loss poses new challenges to high-dimensional distributed estimation in both computation and theoretical development. To address the challenge, we transform the response variable and establish a new connection between quantile regression and ordinary linear regression. Then, we provide a distributed estimator that is both computationally and communicationally efficient, where only the gradient information is communicated at each iteration. Theoretically, we show that, after a constant number of iterations, the proposed estimator achieves a near-oracle convergence rate without any restriction on the number of machines. Moreover, we establish the theoretical guarantee for the support recovery. The simulation analysis is provided to demonstrate the effectiveness of our method.

preprint2020arXiv

EdgeDRNN: Enabling Low-latency Recurrent Neural Network Edge Inference

This paper presents a Gated Recurrent Unit (GRU) based recurrent neural network (RNN) accelerator called EdgeDRNN designed for portable edge computing. EdgeDRNN adopts the spiking neural network inspired delta network algorithm to exploit temporal sparsity in RNNs. It reduces off-chip memory access by a factor of up to 10x with tolerable accuracy loss. Experimental results on a 10 million parameter 2-layer GRU-RNN, with weights stored in DRAM, show that EdgeDRNN computes them in under 0.5 ms. With 2.42 W wall plug power on an entry level USB powered FPGA board, it achieves latency comparable with a 92 W Nvidia 1080 GPU. It outperforms NVIDIA Jetson Nano, Jetson TX2 and Intel Neural Compute Stick 2 in latency by 6X. For a batch size of 1, EdgeDRNN achieves a mean effective throughput of 20.2 GOp/s and a wall plug power efficiency that is over 4X higher than all other platforms.

preprint2020arXiv

Effects of coherence on quantum speed limits and shortcuts to adiabaticity in many-particle systems

We discuss the effects of many-body coherence on the speed of evolution of ultracold atomic gases and the relation to quantum speed limits. Our approach is focused on two related systems, spinless fermions and the bosonic Tonks-Girardeau gas, which possess equivalent density dynamics but very different coherence properties. To illustrate the effect of the coherence on the dynamics we consider squeezing an anharmonic potential which confines the particles and find that the speed of the evolution exhibits subtle, but fundamental differences between the two systems. Furthermore, we explore the difference in the driven dynamics by implementing a shortcut to adiabaticity designed to reduce spurious excitations. We show that collisions between the strongly interacting bosons can lead to changes in the coherence which results in different evolution speeds and therefore different fidelities of the final states.

preprint2020arXiv

Efficient Cysteine Conformer Search with Bayesian Optimization

Finding low-energy molecular conformers is challenging due to the high dimensionality of the search space and the computational cost of accurate quantum chemical methods for determining conformer structures and energies. Here, we combine active-learning Bayesian optimization (BO) algorithms with quantum chemistry methods to address this challenge. Using cysteine as an example, we show that our procedure is both efficient and accurate. After only one thousand single-point calculations and approximately thirty structure relaxations, which is less than 10% computational cost of the current fastest method, we have found the low-energy conformers in good agreement with experimental measurements and reference calculations.

preprint2020arXiv

Epitaxial Growth and Band Structure of Antiferromagnetic Mott Insulator CeOI

The van der Waals material CeOI is predicted to be a layered antiferromagnetic Mott insulator by DFT+U calculation. We successfully grow the CeOI films down to monolayer on graphene/6H-SiC(0001) substrate by using molecular beam epitaxy. Films are studied by {\it in-situ} scanning tunneling microscopy and spectroscopy, which shows a band gap of 4.4 eV. A metallic phase with composition unidentified also exists. This rare earth oxyhalide adds a new member to the two-dimensional magnetic materials.

preprint2020arXiv

Formation of a fast "lane" for positional transition in a microparticle-suspended nematic liquid crystal cell

In this paper, based on the numerical calculation of total energy utilizing the Green's function method, we found that the external electric field applied to a microparticle-suspended nematic liquid crystal cell, if reaching a critical value, combined with its direction, surface anchoring feature and molecular dielectric anisotropy, is possible to create an anisotropic "bubble" around the microparticle with a vertical fast "lane", in which the microparticle can, driven by the asymmetric buoyant force, vertically move swiftly from the cell's midplane to a new equilibrium position, triggering a positional transition discovered by the author previously. Such a new equilibrium position is decided via a competition between the buoyant force and the effective force built upon the microparticle by the elastic energy gradient along the "lane". The threshold value of external field, depends on thickness $L$ and Frank elastic constant $K$ and slightly on the microparticle size and density, in a Fréedericksz-like manner, but by a factor. For a nematic liquid crystal cell with planar surface alignment, a bistable equilibrium structure for the transition is found when the direction of the applied electric field is (a) perpendicular to the two plates of the cell with positive molecular dielectric anisotropy, or (b) parallel to the two plates and the anchoring direction of the cell with negative molecular dielectric anisotropy. Except for the formation of a vertical fast "lane", when the electric field applied is parallel to both the two plates and perpendicular to the anchoring direction, the microparticle suspended in the nematic liquid crystal tends to be trapped in the midplane, regardless of the sign of the molecular dielectric anisotropy. Such phenomenon also occurs for negative molecular dielectric anisotropy while the external is applied perpendicular to the two plates.

preprint2020arXiv

Human-centered collaborative robots with deep reinforcement learning

We present a reinforcement learning based framework for human-centered collaborative systems. The framework is proactive and balances the benefits of timely actions with the risk of taking improper actions by minimizing the total time spent to complete the task. The framework is learned end-to-end in an unsupervised fashion addressing the perception uncertainties and decision making in an integrated manner. The framework is shown to provide more fluent coordination between human and robot partners on an example task of packaging compared to alternatives for which perception and decision-making systems are learned independently, using supervised learning. The foremost benefit of the proposed approach is that it allows for fast adaptation to new human partners and tasks since tedious annotation of motion data is avoided and the learning is performed on-line.

preprint2020arXiv

Inertia indices and eigenvalue inequalities for Hermitian matrices

We present a characterization of eigenvalue inequalities between two Hermitian matrices by means of inertia indices. As applications, we deal with some classical eigenvalue inequalities for Hermitian matrices, including the Cauchy interlacing theorem and the Weyl inequality, in a simple and unified approach. We also give a common generalization of eigenvalue inequalities for (Hermitian) normalized Laplacian matrices of simple (signed, weighted, directed) graphs. Our approach is also suitable for Hermitian matrices of the second kind of digraphs recently introduced by Mohar.

preprint2020arXiv

Learning and Testing Junta Distributions with Subcube Conditioning

We study the problems of learning and testing junta distributions on $\{-1,1\}^n$ with respect to the uniform distribution, where a distribution $p$ is a $k$-junta if its probability mass function $p(x)$ depends on a subset of at most $k$ variables. The main contribution is an algorithm for finding relevant coordinates in a $k$-junta distribution with subcube conditioning [BC18, CCKLW20]. We give two applications: 1. An algorithm for learning $k$-junta distributions with $\tilde{O}(k/ε^2) \log n + O(2^k/ε^2)$ subcube conditioning queries, and 2. An algorithm for testing $k$-junta distributions with $\tilde{O}((k + \sqrt{n})/ε^2)$ subcube conditioning queries. All our algorithms are optimal up to poly-logarithmic factors. Our results show that subcube conditioning, as a natural model for accessing high-dimensional distributions, enables significant savings in learning and testing junta distributions compared to the standard sampling model. This addresses an open question posed by Aliakbarpour, Blais, and Rubinfeld [ABR17].

preprint2020arXiv

Non-iterative Simultaneous Rigid Registration Method for Serial Sections of Biological Tissue

In this paper, we propose a novel non-iterative algorithm to simultaneously estimate optimal rigid transformation for serial section images, which is a key component in volume reconstruction of serial sections of biological tissue. In order to avoid error accumulation and propagation caused by current algorithms, we add extra condition that the position of the first and the last section images should remain unchanged. This constrained simultaneous registration problem has not been solved before. Our algorithm method is non-iterative, it can simultaneously compute rigid transformation for a large number of serial section images in a short time. We prove that our algorithm gets optimal solution under ideal condition. And we test our algorithm with synthetic data and real data to verify our algorithm's effectiveness.

preprint2020arXiv

On Stationary-Point Hitting Time and Ergodicity of Stochastic Gradient Langevin Dynamics

Stochastic gradient Langevin dynamics (SGLD) is a fundamental algorithm in stochastic optimization. Recent work by Zhang et al. [2017] presents an analysis for the hitting time of SGLD for the first and second order stationary points. The proof in Zhang et al. [2017] is a two-stage procedure through bounding the Cheeger's constant, which is rather complicated and leads to loose bounds. In this paper, using intuitions from stochastic differential equations, we provide a direct analysis for the hitting times of SGLD to the first and second order stationary points. Our analysis is straightforward. It only relies on basic linear algebra and probability theory tools. Our direct analysis also leads to tighter bounds comparing to Zhang et al. [2017] and shows the explicit dependence of the hitting time on different factors, including dimensionality, smoothness, noise strength, and step size effects. Under suitable conditions, we show that the hitting time of SGLD to first-order stationary points can be dimension-independent. Moreover, we apply our analysis to study several important online estimation problems in machine learning, including linear regression, matrix factorization, and online PCA.

preprint2020arXiv

On the Sample Complexity of Reinforcement Learning with Policy Space Generalization

We study the optimal sample complexity in large-scale Reinforcement Learning (RL) problems with policy space generalization, i.e. the agent has a prior knowledge that the optimal policy lies in a known policy space. Existing results show that without a generalization model, the sample complexity of an RL algorithm will inevitably depend on the cardinalities of state space and action space, which are intractably large in many practical problems. To avoid such undesirable dependence on the state and action space sizes, this paper proposes a new notion of eluder dimension for the policy space, which characterizes the intrinsic complexity of policy learning in an arbitrary Markov Decision Process (MDP). Using a simulator oracle, we prove a near-optimal sample complexity upper bound that only depends linearly on the eluder dimension. We further prove a similar regret bound in deterministic systems without the simulator.

preprint2020arXiv

Organized Self-Emulsification toward Structural Color

The formation of water-in-oil-in-water (W/O/W) double emulsions can be well-controlled through an organized self-emulsification mechanism in the presence of rigid bottlebrush amphiphilic block copolymers. Nanoscale water droplets with well-controlled diameters form ordered spatial arrangements within the micron-scale oil droplets. Upon solvent evaporation, solid microspheres with hexagonal close packed nanopore arrays are obtained resulting in bright structural colors. The reflected color is precisely tunable across the whole visible light range through tailoring contour length of the bottlebrush molecule. In-situ observation of the W/O interface using confocal laser scanning microscopy provides insights into the mechanism of the organized self-emulsification. This work provides a powerful strategy for the fabrication of structural colored materials in an easy and scalable manner.

preprint2020arXiv

Polynomial-time trace reconstruction in the smoothed complexity model

In the \emph{trace reconstruction problem}, an unknown source string $x \in \{0,1\}^n$ is sent through a probabilistic \emph{deletion channel} which independently deletes each bit with probability $δ$ and concatenates the surviving bits, yielding a \emph{trace} of $x$. The problem is to reconstruct $x$ given independent traces. This problem has received much attention in recent years both in the worst-case setting where $x$ may be an arbitrary string in $\{0,1\}^n$ \cite{DOS17,NazarovPeres17,HHP18,HL18,Chase19} and in the average-case setting where $x$ is drawn uniformly at random from $\{0,1\}^n$ \cite{PeresZhai17,HPP18,HL18,Chase19}. This paper studies trace reconstruction in the \emph{smoothed analysis} setting, in which a ``worst-case'' string $x^{\worst}$ is chosen arbitrarily from $\{0,1\}^n$, and then a perturbed version $\bx$ of $x^{\worst}$ is formed by independently replacing each coordinate by a uniform random bit with probability $σ$. The problem is to reconstruct $\bx$ given independent traces from it. Our main result is an algorithm which, for any constant perturbation rate $0<σ< 1$ and any constant deletion rate $0 < δ< 1$, uses $\poly(n)$ running time and traces and succeeds with high probability in reconstructing the string $\bx$. This stands in contrast with the worst-case version of the problem, for which $\text{exp}(O(n^{1/3}))$ is the best known time and sample complexity \cite{DOS17,NazarovPeres17}. Our approach is based on reconstructing $\bx$ from the multiset of its short subwords and is quite different from previous algorithms for either the worst-case or average-case versions of the problem. The heart of our work is a new $\poly(n)$-time procedure for reconstructing the multiset of all $O(\log n)$-length subwords of any source string $x\in \{0,1\}^n$ given access to traces of $x$.

preprint2020arXiv

Quantum Advantage in Cryptography with a Low-Connectivity Quantum Annealer

The application in cryptography of quantum algorithms for prime factorization fostered the interest in quantum computing. However, quantum computers, and particularly quantum annealers, can also be helpful to construct secure cryptographic keys. Indeed, finding robust Boolean functions for cryptography is an important problem in sequence ciphers, block ciphers, and hash functions, among others. Due to the super-exponential size $\mathcal{O}(2^{2^n})$ of the associated space, finding $n$-variable Boolean functions with global cryptographic constraints is computationally hard. This problem has already been addressed employing generic low-connected incoherent D-Wave quantum annealers. However, the limited connectivity of the Chimera graph, together with the exponential growth in the complexity of the Boolean function design problem, limit the problem scalability. Here, we propose a special-purpose coherent quantum annealing architecture with three couplers per qubit, designed to optimally encode the bent function design problem. A coherent quantum annealer with this tree-type architecture has the potential to solve the $8$-variable bent function design problem, which is classically unsolved, with only $127$ physical qubits and $126$ couplers. This paves the way to reach useful quantum supremacy within the framework of quantum annealing for cryptographic purposes.

preprint2020arXiv

Quantum computing cryptography: Finding cryptographic Boolean functions with quantum annealing by a 2000 qubit D-wave quantum computer

As the building block in symmetric cryptography, designing Boolean functions satisfying multiple properties is an important problem in sequence ciphers, block ciphers, and hash functions. However, the search of $n$-variable Boolean functions fulfilling global cryptographic constraints is computationally hard due to the super-exponential size $\mathcal{O}(2^{2^n})$ of the space. Here, we introduce a codification of the cryptographically relevant constraints in the ground state of an Ising Hamiltonian, allowing us to naturally encode it in a quantum annealer, which seems to provide a quantum speedup. Additionally, we benchmark small $n$ cases in a D-Wave machine, showing its capacity of devising bent functions, the most relevant set of cryptographic Boolean functions. We have complemented it with local search and chain repair to improve the D-Wave quantum annealer performance related to the low connectivity. This work shows how to codify super-exponential cryptographic problems into quantum annealers and paves the way for reaching quantum supremacy with an adequately designed chip.

preprint2020arXiv

Radial Coverage Strength for Optimization of Multi-Camera Deployment

In this paper, a new concept, radial coverage strength, is first proposed to characterize the visual sensing performance when the orientation of the target pose is considered. In particular, the elevation angle of the optical pose of the visual sensor is taken to decompose the visual coverage strength into effective and ineffective components, motivated by the imaging intuition. An optimization problem is then formulated for a multi-camera network to maximize the coverage of the object area based on the strength information fusion along the effective coverage strength direction through the deployment of the angle between radial coverage vector of the camera optical pose. Both simulation and experiments are conducted to validate the proposed approach and comparison with existing methods is also provided.

preprint2020arXiv

Reply to "Comment on 'Apical charge flux-modulated in-plane transport properties of cuprate superconductors'"

This Reply to preceding Comment of arXiv:1909.09867 shows why the statements in the Comment are misleading. We point out that our physical picture and theirs are fundamentally different, therefore the claim of using their correlation to include ours shows a very limited physical relevance, which instead may impede the precise understanding of either of the two pictures.

preprint2020arXiv

Resolution Adaptive Networks for Efficient Inference

Adaptive inference is an effective mechanism to achieve a dynamic tradeoff between accuracy and computational cost in deep networks. Existing works mainly exploit architecture redundancy in network depth or width. In this paper, we focus on spatial redundancy of input samples and propose a novel Resolution Adaptive Network (RANet), which is inspired by the intuition that low-resolution representations are sufficient for classifying "easy" inputs containing large objects with prototypical features, while only some "hard" samples need spatially detailed information. In RANet, the input images are first routed to a lightweight sub-network that efficiently extracts low-resolution representations, and those samples with high prediction confidence will exit early from the network without being further processed. Meanwhile, high-resolution paths in the network maintain the capability to recognize the "hard" samples. Therefore, RANet can effectively reduce the spatial redundancy involved in inferring high-resolution inputs. Empirically, we demonstrate the effectiveness of the proposed RANet on the CIFAR-10, CIFAR-100 and ImageNet datasets in both the anytime prediction setting and the budgeted batch classification setting.

preprint2020arXiv

Retrieving Quantum Information with Active Learning

Active learning is a machine learning method aiming at optimal design for model training. At variance with supervised learning, which labels all samples, active learning provides an improved model by labeling samples with maximal uncertainty according to the estimation model. Here, we propose the use of active learning for efficient quantum information retrieval, which is a crucial task in the design of quantum experiments. Meanwhile, when dealing with large data output, we employ active learning for the sake of classification with minimal cost in fidelity loss. Indeed, labeling only 5% samples, we achieve almost 90% rate estimation. The introduction of active learning methods in the data analysis of quantum experiments will enhance applications of quantum technologies.

preprint2020arXiv

Screening and understanding Li adsorption on 2-dimensional metallic materials by learning physics

Two-dimensional (2D) materials have received considerable attention as possible electrodes in Li-ion batteries (LIBs), although a deeper understanding of the Li adsorption behavior as well as broad screening of the materials space is still needed. In this work, we build a high-throughput screening scheme that incorporates a learned interaction. First, density functional theory and graph convolution networks are utilized to calculate minimum Li adsorption energies for a small set of 2D metallic materials. The data is then used to find a dependence of the minimum Li adsorption energies on the sum of ionization potential, work function of the 2D metal, and coupling energy between Li+ and substrate. Our results show that variances of elemental properties and density are the most correlated features with coupling. To illustrate the applicability of this approach, the model is employed to show that some fluorides and chromium oxides are potential high-voltage materials with adsorption energies < -7 eV, and the found physics is used as the design principle to enhance the Li adsorption ability of graphene. This physics-driven approach shows higher accuracy and transferability compared with purely data-driven models.

preprint2020arXiv

Shortcuts to adiabaticity for an interacting Bose-Einstein condensate via exact solutions of the generalized Ermakov equation

Shortcuts to adiabatic expansion of the effectively one-dimensional Bose-Einstein condensate (BEC) loaded in the harmonic-oscillator (HO) trap is investigated by combining techniques of the variational approximation and inverse engineering. Piecewise-constant (discontinuous) intermediate trap frequencies, similar to the known bang-bang forms in the optimal-control theory, are derived from an exact solution of a generalized Ermakov equation. Control schemes considered in the paper include imaginary trap frequencies at short time scales, i.e., the HO potential replaced by the quadratic repulsive one. Taking into regard the BEC's intrinsic nonlinearity, results are reported for the minimal transfer time, excitation energy (which measures deviation from the effective adiabaticity), and stability for the shortcut-to-adiabaticity protocols. These results are not only useful for the realization of fast frictionless cooling, but also help to address fundamental problems of the quantum speed limit and thermodynamics.

preprint2020arXiv

Shortcuts to Adiabaticity in Digitized Adiabatic Quantum Computing

Shortcuts to adiabaticity are well-known methods for controlling the quantum dynamics beyond the adiabatic criteria, where counter-diabatic (CD) driving provides a promising means to speed up quantum many-body systems. In this work, we show the applicability of CD driving to enhance the digitized adiabatic quantum computing paradigm in terms of fidelity and total simulation time. We study the state evolution of an Ising spin chain using the digitized version of the standard CD driving and its variants derived from the variational approach. We apply this technique in the preparation of Bell and Greenberger-Horne-Zeilinger states with high fidelity using a very shallow quantum circuit. We implement this proposal in the IBM quantum computer, proving its usefulness for the speed up of adiabatic quantum computing in noisy intermediate-scale quantum devices.

preprint2020arXiv

Smooth bang-bang shortcuts to adiabaticity for atomic transport in a moving harmonic trap

Bang-bang control is often used to implement a minimal-time shortcut to adiabaticity for efficient transport of atoms in a moving harmonic trap. However, drastic changes of the on-off controller, leading to high transport-mode excitation and energy consumption, become infeasible under realistic experimental conditions. To circumvent these problems, we propose smooth bang-bang protocols with near-minimal time, by setting the physical constraints on the relative displacement, speed, and acceleration between the mass center of the atom and the trap center. We adopt Pontryagin's maximum principle to obtain the analytical solutions of smooth bang-bang protocol for near-time-minimal control. More importantly, it is found that the energy excitation and sloshing amplitude are significantly reduced at the expense of operation time. We also present a multiple shooting method for the self-consistent numerical analysis. Finally, this method is applied to other tasks, e.g., energy minimization, where obtaining smooth analytical form is complicated.

preprint2020arXiv

State-Aware Tracker for Real-Time Video Object Segmentation

In this work, we address the task of semi-supervised video object segmentation(VOS) and explore how to make efficient use of video property to tackle the challenge of semi-supervision. We propose a novel pipeline called State-Aware Tracker(SAT), which can produce accurate segmentation results with real-time speed. For higher efficiency, SAT takes advantage of the inter-frame consistency and deals with each target object as a tracklet. For more stable and robust performance over video sequences, SAT gets awareness for each state and makes self-adaptation via two feedback loops. One loop assists SAT in generating more stable tracklets. The other loop helps to construct a more robust and holistic target representation. SAT achieves a promising result of 72.3% J&F mean with 39 FPS on DAVIS2017-Val dataset, which shows a decent trade-off between efficiency and accuracy. Code will be released at github.com/MegviiDetection/video_analyst.

preprint2020arXiv

Suppression of Coriolis error in weak equivalence principle test using ^[85]Rb-^[87]Rb dual-species atom interferometer

Coriolis effect is an important error source in the weak equivalence principle (WEP) test using atom interferometer. In this paper, the problem of Coriolis error in WEP test is studied theoretically and experimentally. In theoretical simulation, Coriolis effect is analyzed by establishing an error model. The measurement errors of Eotvos coefficient (eta) in WEP test related to experimental parameters, such as horizontal-velocity difference and horizontal-position difference of atomic clouds, horizontal-position difference of detectors and rotation compensation of Raman laser's mirror are calculated. In experimental investigation, the position difference between Rb-85 and Rb-87 atomic clouds is reduced to 0.1 mm by optimizing the experimental parameters, an alternating detection method is used to suppress the error caused by detection position difference, thus the Coriolis error related to atomic clouds and detectors is eliminated to 1.1E-9. This Coriolis error is further corrected by compensating the rotation of Raman laser's mirror, and the total uncertainty of eta measurement related to Coriolis effect is reduced as 4.4E-11.

preprint2020arXiv

The radio properties of the OH megamaser galaxy IRAS 02524+2046

We present results from VLBI observations of continuum and OH line emission in IRAS 02524+2046 and also arcsecond-scale radio properties of this galaxy using VLA archive data. We found that there is no significant detection of radio continuum emission from VLBI observations. The arcsecond-scale radio images of this source show no clear extended emission, the total radio flux density at L and C band are around 2.9 mJy and 1.0 mJy respectively, which indicate a steep radio spectral index between the two band. Steep spectral index, low brightness temperature and high $q$-ratio (the FIR to the radio flux density), which are three critical indicators in classification of radio activity in the nuclei of galaxies, are all consistent with the classification of this source as a starburst galaxy from its optical spectrum. The high-resolution line profile show that both of \textbf{the 1665 and 1667 MHz OH maser} line have been detected which show three and two clear components respectively. The channel maps show that the maser emission are distributed in a region $\sim$ 210 pc $\times$ 90 pc, the detected maser components at different region show similar double spectral feature, which might be an evidence that this galaxy is at a stage of major merger as seen from the optical morphology.

preprint2020arXiv

Thresholding Bandit Problem with Both Duels and Pulls

The Thresholding Bandit Problem (TBP) aims to find the set of arms with mean rewards greater than a given threshold. We consider a new setting of TBP, where in addition to pulling arms, one can also \emph{duel} two arms and get the arm with a greater mean. In our motivating application from crowdsourcing, dueling two arms can be more cost-effective and time-efficient than direct pulls. We refer to this problem as TBP with Dueling Choices (TBP-DC). This paper provides an algorithm called Rank-Search (RS) for solving TBP-DC by alternating between ranking and binary search. We prove theoretical guarantees for RS, and also give lower bounds to show the optimality of it. Experiments show that RS outperforms previous baseline algorithms that only use pulls or duels.

preprint2020arXiv

Transformer with Bidirectional Decoder for Speech Recognition

Attention-based models have made tremendous progress on end-to-end automatic speech recognition(ASR) recently. However, the conventional transformer-based approaches usually generate the sequence results token by token from left to right, leaving the right-to-left contexts unexploited. In this work, we introduce a bidirectional speech transformer to utilize the different directional contexts simultaneously. Specifically, the outputs of our proposed transformer include a left-to-right target, and a right-to-left target. In inference stage, we use the introduced bidirectional beam search method, which can not only generate left-to-right candidates but also generate right-to-left candidates, and determine the best hypothesis by the score. To demonstrate our proposed speech transformer with a bidirectional decoder(STBD), we conduct extensive experiments on the AISHELL-1 dataset. The results of experiments show that STBD achieves a 3.6\% relative CER reduction(CERR) over the unidirectional speech transformer baseline. Besides, the strongest model in this paper called STBD-Big can achieve 6.64\% CER on the test set, without language model rescoring and any extra data augmentation strategies.

preprint2020arXiv

Uncertainty Quantification for Demand Prediction in Contextual Dynamic Pricing

Data-driven sequential decision has found a wide range of applications in modern operations management, such as dynamic pricing, inventory control, and assortment optimization. Most existing research on data-driven sequential decision focuses on designing an online policy to maximize the revenue. However, the research on uncertainty quantification on the underlying true model function (e.g., demand function), a critical problem for practitioners, has not been well explored. In this paper, using the problem of demand function prediction in dynamic pricing as the motivating example, we study the problem of constructing accurate confidence intervals for the demand function. The main challenge is that sequentially collected data leads to significant distributional bias in the maximum likelihood estimator or the empirical risk minimization estimate, making classical statistics approaches such as the Wald's test no longer valid. We address this challenge by developing a debiased approach and provide the asymptotic normality guarantee of the debiased estimator. Based this the debiased estimator, we provide both point-wise and uniform confidence intervals of the demand function.

preprint2020arXiv

Variable Skipping for Autoregressive Range Density Estimation

Deep autoregressive models compute point likelihood estimates of individual data points. However, many applications (i.e., database cardinality estimation) require estimating range densities, a capability that is under-explored by current neural density estimation literature. In these applications, fast and accurate range density estimates over high-dimensional data directly impact user-perceived performance. In this paper, we explore a technique, variable skipping, for accelerating range density estimation over deep autoregressive models. This technique exploits the sparse structure of range density queries to avoid sampling unnecessary variables during approximate inference. We show that variable skipping provides 10-100$\times$ efficiency improvements when targeting challenging high-quantile error metrics, enables complex applications such as text pattern matching, and can be realized via a simple data augmentation procedure without changing the usual maximum likelihood objective.

preprint2020arXiv

Versatile Mixed Methods for Non-Isothermal Incompressible Flows

The purpose of this paper is to extend the versatile mixed methods originally developed by Chen and Williams for isothermal flows in "Versatile Mixed Methods for the Incompressible Navier-Stokes Equations," Computers & Mathematics with Applications, 2020, (under review), to simulate non-isothermal incompressible flows. These new mixed methods are particularly interesting, as with only minor modifications they can be applied to a much broader range of flows, including non-isothermal weakly-compressible flows, and fully-compressible flows. In the main body of this paper, we carefully develop these mixed methods for solving the Boussinesq model equations. Thereafter, we prove the L2-stability of the discrete temperature field, and assess the practical behavior of the methods by applying them to a set of well-known convection problems.

preprint2020arXiv

Versatile Mixed Methods for the Incompressible Navier-Stokes Equations

In the spirit of the "Principle of Equipresence" introduced by Truesdell & Toupin, The Classical Field Theories (1960), we use the full version of the viscous stress tensor which was originally derived for compressible flows, instead of the classical incompressible stress tensor. In our approach, the divergence-free constraint for the viscous stress term is not enforced ahead of discretization. Instead, our formulation allows the scheme itself to "choose" a consistent way to interpret the divergence-free constraint: i.e., the divergence-free constraint is interpreted (or enforced) in a consistent fashion in both the mass conservation equation and the stress tensor term (in the momentum equation). Furthermore, our approach preserves the original symmetrical properties of the stress tensor, e.g. its rotational invariance, and it remains physically correct in the context of compressible flows. As a result, our approach facilitates versatility and code reuse. In this paper, we introduce our approach and establish some important mathematical properties for the resulting class of finite element schemes. More precisely, for general mixed methods, which are not necessarily pointwise divergence-free, we establish the existence of a new norm induced by the full, viscous bilinear form. Thereafter, we prove the coercivity of the viscous bilinear form and the semi-coercivity of a convective trilinear form. In addition, we demonstrate L2-stability of the discrete velocity fields for the general class of methods and (by deduction) the H(div)-conforming methods. Finally, we run some numerical experiments to illustrate the behavior of the versatile mixed methods, and we make careful comparisons with a conventional H(div)-conforming scheme.

preprint2020arXiv

Wolf phase tomography (WPT) of transparent structures using partially coherent illumination

Diffraction tomography using coherent holographic imaging has been proposed by Emil Wolf in 1969 to extract 3D information from transparent, inhomogeneous objects. At the same time, the Wolf equations describe the propagation correlations associated with partially coherent fields. Combining these two concepts, here we present Wolf phase tomography (WPT), which is a method for performing diffraction tomography using partially coherent fields. The WPT reconstruction works in the direct space-time domain, without the need of Fourier transformation, and decouples the refractive index distribution from the thickness of the sample. We demonstrate the WPT principle using data acquired by spatial light interference microscopy (SLIM). SLIM is a quantitative phase imaging method that upgrades an existing phase contrast microscope by introducing controlled phase shifts between the incident and scattered fields. The illumination field in SLIM is spatially partially coherent (emerging from a ring-shaped pupil function) and of low temporal coherence (white light), thus, suitable for the Wolf equations. From three intensity measurements corresponding to different phase-contrast frames in SLIM, the 3D refractive index distribution is obtained right away by computing the Laplacian and second time derivative of the measured complex correlation function. The high-throughput and simplicity of this method enables the study of 3D, dynamic events in living cells across entire multi-well plates, with RI sensitivity on the order of 10-5. We validate WPT with measurements of standard samples (microbeads), spermatozoa, and live neural networks.

preprint2020arXiv

Zariski's conjecture and Euler-Chow series

We study the relations between the finite generation of Cox ring, the rationality of Euler-Chow series and Poincaré series and Zariski's conjecture on dimensions of linear systems. We prove that if the Cox ring of a smooth projective variety is finitely generated, then all Poincaré series of the variety are rational. We also prove that the multi-variable Poincaré series associated to big divisors on a smooth projective surface are rational, assuming the rationality of multi-variable Poincare series on curves.

preprint2019arXiv

All-optical frequency resolved optical gating for isolated attosecond pulse reconstruction

We demonstrate an all-optical approach for precise characterization of attosecond extreme ultraviolet pulses. Isolated attosecond pulse is produced from high order harmonics using intense driving pulse with proper gating technique. When a weak field is synchronized with the driver, it perturbs the harmonics generation process via altering the accumulated phase of the electron trajectories. The perturbed harmonic spectrum can be formulated as a convolution of the unperturbed dipole and a phase gate, implying the validity of complete reconstruction of isolated attosecond pulses using conventional frequency resolved optical gating method. This in situ measurement avoids the central momentum approximation assumed in the widely used attosecond streaking measurement, providing a simple and reliable metrology for isolated attosecond pulse.

preprint2019arXiv

Bayesian Dynamic Modeling and Monitoring of Network Flows

In the context of a motivating study of dynamic network flow data on a large-scale e-commerce web site, we develop Bayesian models for on-line/sequential analysis for monitoring and adapting to changes reflected in node-node traffic. For large-scale networks, we customize core Bayesian time series analysis methods using dynamic generalized linear models (DGLMs). These are integrated into the context of multivariate networks using the concept of decouple/recouple that was recently introduced in multivariate time series. This method enables flexible dynamic modeling of flows on large-scale networks and exploitation of partial parallelization of analysis while maintaining coherence with an over-arching multivariate dynamic flow model. This approach is anchored in a case-study on internet data, with flows of visitors to a commercial news web site defining a long time series of node-node counts on over 56,000 node pairs. Central questions include characterizing inherent stochasticity in traffic patterns, understanding node-node interactions, adapting to dynamic changes in flows and allowing for sensitive monitoring to flag anomalies. The methodology of dynamic network DGLMs applies to many dynamic network flow studies.

preprint2019arXiv

Coupled density-spin Bose-Einstein condensates dynamics and collapse in systems with quintic nonlinearity

We investigate the effects of spin-orbit coupling and Zeeman splitting on the coupled density-spin dynamics and collapse of the Bose-Einstein condensate driven by the quintic self-attraction in the same- and cross-spin channels. The characteristic feature of the collapse is the decrease in the width as given by the participation ratio of the density rather than by the expectation values of the coordinate. Qualitative arguments and numerical simulations reveal the existence of a critical spin-orbit coupling strength which either prohibits or leads to the collapse, and its dependence on other parameters, such as the condensates norm, spin-dependent nonlinear coupling, and the Zeeman splitting. The entire nonlinear dynamics critically depend on the initial spin sate.

preprint2019arXiv

Experimental Observation of Equilibrium and Dynamical Quantum Phase Transitions via Out-of-Time-Ordered Correlators

The out-of-time-ordered correlators (OTOC) have been established as a fundamental concept for quantifying quantum information scrambling and diagnosing quantum chaotic behavior. Recently, it was theoretically proposed that the OTOC can be used as an order parameter to dynamically detect both equilibrium quantum phase transitions (EQPTs) and dynamical quantum phase transitions (DQPTs) in one-dimensional many-body systems. Here we report the first experimental observation of EQPTs and DQPTs in a quantum spin chain via quench dynamics of OTOC on a nuclear magnetic resonance quantum simulator. We observe that the quench dynamics of both the order parameter and the two-body correlation function cannot detect the DQPTs, but the OTOC can unambiguously detect the DQPTs. Moreover, we demonstrate that the long-time average value of the OTOC in quantum quench signals the equilibrium quantum critical point and ordered quantum phases, thus one can measure the EQPTs from the non-equilibrium quantum quench dynamics. Our experiment paves a way for experimentally investigating DQPTs through OTOCs and for studying the EQPTs through the non-equilibrium quantum quench dynamics with quantum simulators.

preprint2019arXiv

Ferromagnetic spin correlations in the two-dimensional Hubbard model

We analyze the dynamical nearest-neighbor and next-nearest-neighbor spin correlations in the 4-site and 8-site dynamical cluster approximation to the two-dimensional Hubbard model. Focusing on the robustness of these correlations at long imaginary times, we reveal enhanced ferromagnetic correlations on the lattice diagonal, consistent with the emergence of composite spin-1 moments at a temperature scale that essentially coincides with the pseudo-gap temperature $T^*$. We discuss these results in the context of the spin-freezing theory of unconventional superconductivity.

preprint2019arXiv

Lecture Notes of Tensor Network Contractions

Tensor network (TN), a young mathematical tool of high vitality and great potential, has been undergoing extremely rapid developments in the last two decades, gaining tremendous success in condensed matter physics, atomic physics, quantum information science, statistical physics, and so on. In this lecture notes, we focus on the contraction algorithms of TN as well as some of the applications to the simulations of quantum many-body systems. Starting from basic concepts and definitions, we first explain the relations between TN and physical problems, including the TN representations of classical partition functions, quantum many-body states (by matrix product state, tree TN, and projected entangled pair state), time evolution simulations, etc. These problems, which are challenging to solve, can be transformed to TN contraction problems. We present then several paradigm algorithms based on the ideas of the numerical renormalization group and/or boundary states, including density matrix renormalization group, time-evolving block decimation, coarse-graining/corner tensor renormalization group, and several distinguished variational algorithms. Finally, we revisit the TN approaches from the perspective of multi-linear algebra (also known as tensor algebra or tensor decompositions) and quantum simulation. Despite the apparent differences in the ideas and strategies of different TN algorithms, we aim at revealing the underlying relations and resemblances in order to present a systematic picture to understand the TN contraction approaches.

preprint2019arXiv

Reentrance of Topological Phase in Spin-1 Frustrated Heisenberg Chain

For the Haldane phase, the magnetic field usually tends to break the symmetry and drives the system into a topologically trivial phase. Here, we report a novel reentrance of the Haldane phase at zero temperature in the spin-1 antiferromagnetic Heisenberg model on sawtooth chain. A partial Haldane phase is induced by the magnetic field, which is the combination of the Haldane state in one sublattice and a ferromagnetically ordered state in the other sublattice. Such a partial topological order is a result of the zero-temperature entropy due to quantum fluctuations caused by geometrical frustration.

preprint2019arXiv

Towards Pricing Financial Derivatives with an IBM Quantum Computer

Pricing interest-rate financial derivatives is a major problem in finance, in which it is crucial to accurately reproduce the time-evolution of interest rates. Several stochastic dynamics have been proposed in the literature to model either the instantaneous interest rate or the instantaneous forward rate. A successful approach to model the latter is the celebrated Heath-Jarrow-Morton framework, in which its dynamics is entirely specified by volatility factors. On its multifactor version, this model considers several noisy components to capture at best the dynamics of several time-maturing forward rates. However, as no general analytical solution is available, there is a trade-off between the number of noisy factors considered and the computational time to perform a numerical simulation. Here, we employ the quantum principal component analysis to reduce the number of noisy factors required to accurately simulate the time evolution of several time-maturing forward rates. The principal components are experimentally estimated with the $5$-qubit IBMQX2 quantum computer for $2\times 2$ and $3\times 3$ cross-correlation matrices, which are based on historical data for two and three time-maturing forward rates. This manuscript is a first step towards the design of a general quantum algorithm to fully simulate on quantum computers the Heath-Jarrow-Morton model for pricing interest-rate financial derivatives. It shows indeed that practical applications of quantum computers in finance will be achievable in the near future.

preprint2019arXiv

ZAIGA: Zhaoshan Long-baseline Atom Interferometer Gravitation Antenna

The Zhaoshan long-baseline Atom Interferometer Gravitation Antenna (ZAIGA) is a new type of underground laser-linked interferometer facility, and is currently under construction. It is in the 200-meter-on-average underground of a mountain named Zhaoshan which is about 80 km southeast to Wuhan. ZAIGA will be equipped with long-baseline atom interferometers, high-precision atom clocks, and large-scale gyros. ZAIGA facility will take an equilateral triangle configuration with two 1-km-apart atom interferometers in each arm, a 300-meter vertical tunnel with atom fountain and atom clocks mounted, and a tracking-and-ranging 1-km-arm-length prototype with lattice optical clocks linked by locked lasers. The ZAIGA facility will be used for experimental research on gravitation and related problems including gravitational wave detection, high-precision test of the equivalence principle of micro-particles, clock based gravitational red-shift measurement, rotation measurement and gravito-magnetic effect.

preprint2016arXiv

Aerodynamic Heating in Hypersonic Boundary Layers:\ Role of Dilatational Waves

The evolution of multi-mode instabilities in a hypersonic boundary layer and their effects on aerodynamic heating are investigated. Experiments are conducted in a Mach 6 wind tunnel using Rayleigh-scattering flow visualization, fast-response pressure sensors, fluorescent temperature-sensitive paint (TSP), and particle image velocimetry (PIV). Calculations are also performed based on both parabolized stability equations (PSE) and direct numerical simulations (DNS). It is found that second-mode dilatational waves, accompanied by high-frequency alternating fluid compression and expansion, produce intense aerodynamic heating in a small region that rapidly heats the fluid passing through it. As a result, the surface temperature rapidly increases and results in an overshoot over the nominal transitional value. When the dilatation waves decay downstream, the surface temperature decreases gradually until transition is completed. A theoretical analysis is provided to interpret the temperature distribution affected by the aerodynamic heating.

preprint2016arXiv

Analytic prediction for planar turbulent boundary layers

Analytic predictions of mean velocity profile (MVP) and streamwise ($x$) development of related integral quantities are presented for flows in channel and turbulent boundary layer (TBL), based on a symmetry analysis of eddy length and total stress. Specific predictions are the friction velocity $u_τ$: ${ U_e/u_τ}\approx 2.22\ln Re_x+2.86-3.83\ln(\ln Re_x)$; the boundary layer thickness $δ_e$: $x/δ_e \approx 7.27\ln Re_x-5.18-12.52\ln(\ln Re_x)$; the momentum thickness Reynolds number: $Re_x/Re_θ=4.94[{(\ln {{\mathop{\rm Re}\nolimits} _θ } + 1.88)^2} + 1]$, all in good agreement with empirical data.

preprint2016arXiv

Architectures of Soft Robotic Locomotion Enabled by Simple Mechanical Principles

In nature, a variety of limbless locomotion patterns flourish from the small or basic life form (Escherichia coli, the amoeba, etc.) to the large or intelligent creatures (e.g., slugs, starfishes, earthworms, octopuses, jellyfishes, and snakes). Many bioinspired soft robots based on locomotion have been developed in the past decades. In this work, based on the kinematics and dynamics of two representative locomotion modes (i.e., worm-like crawling and snake-like slithering), we propose a broad set of innovative designs for soft mobile robots through simple mechanical principles. Inspired by and go beyond existing biological systems, these designs include 1-D (dimensional), 2-D, and 3-D robotic locomotion patterns enabled by simple actuation of continuous beams. We report herein over 20 locomotion modes achieving various locomotion functions, including crawling, rising, running, creeping, squirming, slithering, swimming, jumping, turning, turning over, helix rolling, wheeling, etc. Some of them are able to reach high speed, high efficiency, and overcome obstacles. All these locomotion strategies and functions can be integrated into a simple beam model. The proposed simple and robust models are adaptive for severe and complex environments. These elegant designs for diverse robotic locomotion patterns are expected to underpin future deployments of soft robots and to inspire series of advanced designs.

preprint2016arXiv

Bayesian Decision Process for Cost-Efficient Dynamic Ranking via Crowdsourcing

Rank aggregation based on pairwise comparisons over a set of items has a wide range of applications. Although considerable research has been devoted to the development of rank aggregation algorithms, one basic question is how to efficiently collect a large amount of high-quality pairwise comparisons for the ranking purpose. Because of the advent of many crowdsourcing services, a crowd of workers are often hired to conduct pairwise comparisons with a small monetary reward for each pair they compare. Since different workers have different levels of reliability and different pairs have different levels of ambiguity, it is desirable to wisely allocate the limited budget for comparisons among the pairs of items and workers so that the global ranking can be accurately inferred from the comparison results. To this end, we model the active sampling problem in crowdsourced ranking as a Bayesian Markov decision process, which dynamically selects item pairs and workers to improve the ranking accuracy under a budget constraint. We further develop a computationally efficient sampling policy based on knowledge gradient as well as a moment matching technique for posterior approximation. Experimental evaluations on both synthetic and real data show that the proposed policy achieves high ranking accuracy with a lower labeling cost.

preprint2016arXiv

Benchmarking Deep Reinforcement Learning for Continuous Control

Recently, researchers have made significant progress combining the advances in deep learning for learning feature representations with reinforcement learning. Some notable examples include training agents to play Atari games based on raw pixel data and to acquire advanced manipulation skills using raw sensory inputs. However, it has been difficult to quantify progress in the domain of continuous control due to the lack of a commonly adopted benchmark. In this work, we present a benchmark suite of continuous control tasks, including classic tasks like cart-pole swing-up, tasks with very high state and action dimensionality such as 3D humanoid locomotion, tasks with partial observations, and tasks with hierarchical structure. We report novel findings based on the systematic evaluation of a range of implemented reinforcement learning algorithms. Both the benchmark and reference implementations are released at https://github.com/rllab/rllab in order to facilitate experimental reproducibility and to encourage adoption by other researchers.

preprint2016arXiv

Bound exciton and free exciton states in GaSe thin slab

The photoluminescence (PL) and absorption experiments have been performed in GaSe slab with incident light polarized perpendicular to c-axis of sample at 10K. An obvious energy difference of about 34meV between exciton absorption peak and PL peak (the highest energy peak) is observed. By studying the temperature dependence of PL spectra, we attribute it to energy difference between free exciton and bound exciton states, where main exciton absorption peak comes from free exciton absorption, and PL peak are attributed to recombination of bound exciton at 10K. This strong bound exciton effect is stable up to 50K. Moreover, the temperature dependence of integrated PL intensity and PL lifetime reveals that a non-radiative process, with active energy extracted as 0.5meV, dominates PL emission.

preprint2016arXiv

Bulk flow scaling for turbulent channel and pipe flows

We report a theory deriving bulk flow scaling for canonical wall-bounded flows. The theory accounts for the symmetries of boundary geometry (flat plate channel versus circular pipe) by a variational calculation for a large-scale energy length, which characterizes its bulk flow scaling by a simple exponent, i.e. $m=4$ for channel and 5 for pipe. The predicted mean velocity shows excellent agreement with several dozen sets of quality empirical data for a wide range of the Reynolds number (Re), with a universal bulk flow constant $κ\approx0.45$. Predictions for dissipation and turbulent transport in the bulk flow are also given, awaiting data verification.

preprint2016arXiv

Comparison of three Statistical Classification Techniques for Maser Identification

We applied three statistical classification techniques - linear discriminant analysis (LDA), logistic regression and random forests - to three astronomical datasets associated with searches for interstellar masers. We compared the performance of these methods in identifying whether specific mid-infrared or millimetre continuum sources are likely to have associated interstellar masers. We also discuss the ease, or otherwise, with which the results of each classification technique can be interpreted. Non-parametric methods have the potential to make accurate predictions when there are complex relationships between critical parameters. We found that for the small datasets the parametric methods logistic regression and LDA performed best, for the largest dataset the non-parametric method of random forests performed with comparable accuracy to parametric techniques, rather than any significant improvement. This suggests that at least for the specific examples investigated here accuracy of the predictions obtained is not being limited by the use of parametric models. We also found that for LDA, transformation of the data to match a normal distribution in the input parameters led to big improvements in accuracy. The different classification techniques had significant overlap in their predictions, further astronomical observations will enable the accuracy of these predictions to be tested.

preprint2016arXiv

Competitive analysis of the top-K ranking problem

Motivated by applications in recommender systems, web search, social choice and crowdsourcing, we consider the problem of identifying the set of top $K$ items from noisy pairwise comparisons. In our setting, we are non-actively given $r$ pairwise comparisons between each pair of $n$ items, where each comparison has noise constrained by a very general noise model called the strong stochastic transitivity (SST) model. We analyze the competitive ratio of algorithms for the top-$K$ problem. In particular, we present a linear time algorithm for the top-$K$ problem which has a competitive ratio of $\tilde{O}(\sqrt{n})$; i.e. to solve any instance of top-$K$, our algorithm needs at most $\tilde{O}(\sqrt{n})$ times as many samples needed as the best possible algorithm for that instance (in contrast, all previous known algorithms for the top-$K$ problem have competitive ratios of $\tildeΩ(n)$ or worse). We further show that this is tight: any algorithm for the top-$K$ problem has competitive ratio at least $\tildeΩ(\sqrt{n})$.

preprint2016arXiv

Effects of intrinsic strain on the structural stability and mechanical properties of phosphorene nanotubes

Using molecular dynamics (MD) simulations, we explore the structural stability and mechanical integrity of phosphorene nanotubes (PNTs), where the intrinsic strain in the tubular PNT structure plays an important role. It is proposed that the atomic structure of larger-diameter armchair PNTs (armPNTs) can remain stable at higher temperature, but the high intrinsic strain in the hoop direction renders zigzag PNTs (zigPNTs) less favorable. The mechanical properties of PNTs, including the Young's modulus and fracture strength, are sensitive to the diameter, showing a size dependence. A simple model is proposed to express the Young's modulus as a function of the intrinsic axial strain which in turns depends on the diameter of PNTs. In addition, the compressive buckling of armPNTs is length-dependent, whose instability modes transit from column buckling to shell buckling are observed as the ratio of diameter/length increases.

preprint2016arXiv

Experimental Observation of the Quantum Anomalous Hall Effect in a Magnetic Topological Insulator

The quantized version of the anomalous Hall effect has been predicted to occur in magnetic topological insulators, but the experimental realization has been challenging. Here, we report the observation of the quantum anomalous Hall (QAH) effect in thin films of Cr-doped (Bi,Sb)2Te3, a magnetic topological insulator. At zero magnetic field, the gate-tuned anomalous Hall resistance reaches the predicted quantized value of h/e^2,accompanied by a considerable drop of the longitudinal resistance. Under a strong magnetic field, the longitudinal resistance vanishes whereas the Hall resistance remains at the quantized value. The realization of the QAH effect may lead to the development of low-power-consumption electronics.

preprint2016arXiv

Experimental observation of topological Fermi arcs in type-II Weyl semimetal MoTe2

Weyl semimetal is a new quantum state of matter [1-12] hosting the condensed matter physics counterpart of relativisticWeyl fermion [13] originally introduced in high energy physics. The Weyl semimetal realized in the TaAs class features multiple Fermi arcs arising from topological surface states [10, 11, 14-16] and exhibits novel quantum phenomena, e.g., chiral anomaly induced negative mag-netoresistance [17-19] and possibly emergent supersymmetry [20]. Recently it was proposed theoretically that a new type (type-II) of Weyl fermion [21], which does not have counterpart in high energy physics due to the breaking of Lorentz invariance, can emerge as topologically-protected touching between electron and hole pockets. Here, we report direct spectroscopic evidence of topological Fermi arcs in the predicted type-II Weyl semimetal MoTe2 [22-24]. The topological surface states are confirmed by directly observing the surface states using bulk-and surface-sensitive angle-resolved photoemission spectroscopy (ARPES), and the quasi-particle interference (QPI) pattern between the two putative Fermi arcs in scanning tunneling microscopy (STM). Our work establishes MoTe2 as the first experimental realization of type-II Weyl semimetal, and opens up new opportunities for probing novel phenomena such as exotic magneto-transport [21] in type-II Weyl semimetals.

preprint2016arXiv

Experimental realization of stimulated Raman shortcut-to-adiabatic passage with cold atoms

Accurate control of a quantum system is a fundamental requirement in many areas of modern science ranging from quantum information processing to high-precision measurements. A significantly important goal in quantum control is to prepare a desired state as fast as possible with sufficiently high fidelity allowed by available resources and experimental constraints. Stimulated Raman adiabatic passage (STIRAP) is a robust way to realize high-fidelity state transfer but it requires a sufficiently long operation time to satisfy the adiabatic criteria. We here theoretically propose and then experimentally demonstrate a shortcut-to-adiabatic protocol to speed up the STIRAP. By modifying the shapes of the Raman pulses, we experimentally realize a fast and high-fidelity stimulated Raman shortcut-to-adiabatic passage that is robust against control parameter variations. The all-optical, robust, and fast protocol demonstrated here provides an efficient and practical way to control quantum systems.

preprint2016arXiv

Extracting the differential phase in dual atom interferometers by modulating magnetic fields

We present a new scheme for measuring the differential phase in dual atom interferometers. The magnetic field is modulated in one interferometer, and the differential phase can be extracted without measuring the amplitude of the magnetic field by combining the ellipse and linear fitting methods. The gravity gradient measurements are discussed based on dual atom interferometers. Numerical simulation shows that the systematic error of the differential phase measurement is largely decreased when the duration of the magnetic field is symmetrically modulated. This combined fitting scheme has a high accuracy for measuring an arbitrary differential phase in dual atom interferometers.

preprint2016arXiv

Generating Discriminative Object Proposals via Submodular Ranking

A multi-scale greedy-based object proposal generation approach is presented. Based on the multi-scale nature of objects in images, our approach is built on top of a hierarchical segmentation. We first identify the representative and diverse exemplar clusters within each scale by using a diversity ranking algorithm. Object proposals are obtained by selecting a subset from the multi-scale segment pool via maximizing a submodular objective function, which consists of a weighted coverage term, a single-scale diversity term and a multi-scale reward term. The weighted coverage term forces the selected set of object proposals to be representative and compact; the single-scale diversity term encourages choosing segments from different exemplar clusters so that they will cover as many object patterns as possible; the multi-scale reward term encourages the selected proposals to be discriminative and selected from multiple layers generated by the hierarchical image segmentation. The experimental results on the Berkeley Segmentation Dataset and PASCAL VOC2012 segmentation dataset demonstrate the accuracy and efficiency of our object proposal model. Additionally, we validate our object proposals in simultaneous segmentation and detection and outperform the state-of-art performance.

preprint2016arXiv

High-temperature superconductivity in single-unit-cell FeSe films on anatase TiO2(001)

We report on the observation of high-temperature ($T_\textrm{c}$) superconductivity and magnetic vortices in single-unit-cell FeSe films on anatase TiO$_2$(001) substrate by using scanning tunneling microscopy. A systematic study and engineering of interfacial properties has clarified the essential roles of substrate in realizing the high-$T_\textrm{c}$ superconductivity, probably via interface-induced electron-phonon coupling enhancement and charge transfer. By visualizing and tuning the oxygen vacancies at the interface, we find their very limited effect on the superconductivity, which excludes interfacial oxygen vacancies as the primary source for charge transfer between the substrate and FeSe films. Our findings have placed severe constraints on any microscopic model for the high-$T_\textrm{c}$ superconductivity in FeSe-related heterostructures.

preprint2016arXiv

Magnons and Phonons Optically Driven Out of Local Equilibrium in a Magnetic Insulator

Magnons are the energy quanta of fundamental spin excitations, namely spin waves, and they can make a considerable contribution to energy transport in some magnetic materials in a similar manner as lattice vibration waves or phonons. The coupling and possible non-equilibrium between magnons and other energy carriers have been used to explain several recently discovered thermally driven spin transport and energy conversion phenomena. Here, we report experiments in which local non-equilibrium between magnons and phonons in a single crystalline bulk magnetic insulator, Y3Fe5O12 (yttrium iron garnet, or YIG), has been created optically within a focused laser spot and probed directly with the use of micro-Brillouin light scattering (BLS). By analyzing the experimental results with a thermally induced magnon diffusion model, we obtain the magnon diffusion length of thermal magnons. By explicitly establishing non-equilibrium between magnons and phonons, our studies represent an important step toward a quantitative understanding of various spin-heat coupling phenomena.

preprint2016arXiv

Molecular Lines of 13 Galactic Infrared Bubble Regions

We investigated the physical properties of molecular clouds and star formation processes around infrared bubbles which are essentially expanding HII regions. We performed observations of 13 galactic infrared bubble fields containing 18 bubbles. Five molecular lines, 12CO (J=1-0), 13CO (J=1-0), C18O(J=1-0), HCN (J=1-0), and HCO+ (J=1-0), were observed, and several publicly available surveys, GLIMPSE, MIPSGAL, ATLASGAL, BGPS, VGPS, MAGPIS, and NVSS, were used for comparison. We find that these bubbles are generally connected with molecular clouds, most of which are giant. Several bubble regions display velocity gradients and broad shifted profiles, which could be due to the expansion of bubbles. The masses of molecular clouds within bubbles range from 100 to 19,000 solar mass, and their dynamic ages are about 0.3-3.7 Myr, which takes into account the internal turbulence pressure of surrounding molecular clouds. Clumps are found in the vicinity of all 18 bubbles, and molecular clouds near four of these bubbles with larger angular sizes show shell-like morphologies, indicating that either collect-and-collapse or radiation-driven implosion processes may have occurred. Due to the contamination of adjacent molecular clouds, only six bubble regions are appropriate to search for outflows, and we find that four of them have outflow activities. Three bubbles display ultra-compact HII regions at their borders, and one of them is probably responsible for its outflow. In total, only six bubbles show star formation activities in the vicinity, and we suggest that star formation processes might have been triggered.

preprint2016arXiv

New method of averaging diffeomorphisms based on Jacobian determinant and curl vector

Averaging diffeomorphisms is a challenging problem, and it has great applications in areas like medical image atlases. The simple Euclidean average can neither guarantee the averaged transformation is a diffeomorphism, nor get reasonable result when there is a local rotation. The goal of this paper is to propose a new approach to averaging diffeomorphisms based on the Jacobian determinant and the curl vector of the diffeomorphisms. Instead of averaging the diffeomorphisms directly, we average the Jacobian determinants and the curl vectors, and then construct a diffeomorphism based on the averaged Jacobian determinant and averaged curl vector as the average of diffeomorphisms. Numerical examples with convincible results are presented to demonstrate the method.

preprint2016arXiv

On Bayes Risk Lower Bounds

This paper provides a general technique for lower bounding the Bayes risk of statistical estimation, applicable to arbitrary loss functions and arbitrary prior distributions. A lower bound on the Bayes risk not only serves as a lower bound on the minimax risk, but also characterizes the fundamental limit of any estimator given the prior knowledge. Our bounds are based on the notion of $f$-informativity, which is a function of the underlying class of probability measures and the prior. Application of our bounds requires upper bounds on the $f$-informativity, thus we derive new upper bounds on $f$-informativity which often lead to tight Bayes risk lower bounds. Our technique leads to generalizations of a variety of classical minimax bounds (e.g., generalized Fano's inequality). Our Bayes risk lower bounds can be directly applied to several concrete estimation problems, including Gaussian location models, generalized linear models, and principal component analysis for spiked covariance models. To further demonstrate the applications of our Bayes risk lower bounds to machine learning problems, we present two new theoretical results: (1) a precise characterization of the minimax risk of learning spherical Gaussian mixture models under the smoothed analysis framework, and (2) lower bounds for the Bayes risk under a natural prior for both the prediction and estimation errors for high-dimensional sparse linear regression under an improper learning setting.

preprint2016arXiv

Prediction of a Two-dimensional Phosphorus Nitride Monolayer

Today, 2D semiconductor materials have been extended into the nitrogen group: phosphorene, arsenene, antimonene and even nitrogene. Motivated by them, based upon first-principles density functional calculations, we propose a new two-dimensional phosphorus nitride (PN) structure that is stable well above the room temperature, due to its extremely high cohesive energy. Unlike phosphorene, PN structure is resistant to high temperature oxidation. The structure is predicted to be a semiconductor with a wide, indirect band gap of 2.64 eV. More interestingly, the phosphorus nitride monolayer experiences an indirect-to-direct band-gap transition at a relatively small tensile strain. Such dramatic transformation in the electronic structure combined with structural stability and oxidation resistance at high temperature could pave the way for exciting innovations in high-speed ultrathin transistors, power electronic modules, ultra-high efficiency LEDs and semiconductor lasers.

preprint2016arXiv

Predictions of canonical wall bounded turbulent flows via a modified $k-ω$ equation

A major challenge in computation of engineering flows is to derive and improve turbulence models built on turbulence physics. Here, we present a physics-based modified $k-ω$ equation for canonical wall bounded turbulent flows (boundary layer, channel and pipe), predicting both mean velocity profile (MVP) and streamwise mean kinetic energy profile (SMKP) with high accuracy over a wide range of Reynolds number ($Re$). The result builds on a multi-layer quantification of wall flows, which allows a significant modification of the $k-ω$ equation. Three innovations are introduced: First, an adjustment of the Karman constant to 0.45 is set for the overlap region with a logarithmic MVP. Second, a wake parameter models the turbulent transport near the centerline. Third, an anomalous dissipation factor represents the effect of a meso layer in the overlap region. Then, a highly accurate (above 99\%) prediction of MVPs is obtained in Princeton pipes, improving the original model prediction by up to 10\%. Moreover, the entire SMKP, including the newly observed outer peak, is predicted. With a slight change of the wake parameter, the model also yields accurate predictions for channels and boundary layers.

preprint2016arXiv

Quantum phase transition, universality and scaling behaviors in the spin-1/2 Heisenberg model with ferromagnetic and antiferromagnetic competing interactions on honeycomb lattice

The quantum phase transition, scaling behaviors, and thermodynamics in the spin-1/2 quantum Heisenberg model with antiferromagnetic coupling $J>0$ in armchair direction and ferromagnetic interaction $J'<0$ in zigzag direction on a honeycomb lattice are systematically studied using the continuous-time quantum Monte Carlo method. By calculating the Binder ratio $Q_{2}$ and spin stiffness $ρ$ in two directions for various coupling ratio $α=J'/J$ under different lattice sizes, we found that a quantum phase transition from the dimerized phase to the stripe phase occurs at the quantum critical point $α_c=-0.93$. Through the finite-size scaling analysis on $Q_{2}$, $ρ_{x}$ and $ρ_{y}$, we determined the critical exponent related to the correlation length $ν$ to be 0.7212(8), implying that this transition falls into a classical Heisenberg O(3) universality. A zero magnetization plateau is observed in the dimerized phase, whose width decreases with increasing $α$. A phase diagram in the coupling ratio $α$-magnetic field $h$ plane is obtained, where four phases, including dimerized, stripe, canted stripe and polarized phases are identified. It is also unveiled that the temperature dependence of the specific heat $C(T)$ for different $α$'s intersects precisely at one point, similar to that of liquid $^{3}$He under different pressures and several magnetic compounds under various magnetic fields. The scaling behaviors of $Q_{2}$, $ρ$ and $C(T)$ are carefully analyzed. The susceptibility is well compared with the experimental data to give the magnetic parameters of both compounds.

preprint2016arXiv

RL$^2$: Fast Reinforcement Learning via Slow Reinforcement Learning

Deep reinforcement learning (deep RL) has been successful in learning sophisticated behaviors automatically; however, the learning process requires a huge number of trials. In contrast, animals can learn new tasks in just a few trials, benefiting from their prior knowledge about the world. This paper seeks to bridge this gap. Rather than designing a "fast" reinforcement learning algorithm, we propose to represent it as a recurrent neural network (RNN) and learn it from data. In our proposed method, RL$^2$, the algorithm is encoded in the weights of the RNN, which are learned slowly through a general-purpose ("slow") RL algorithm. The RNN receives all information a typical RL algorithm would receive, including observations, actions, rewards, and termination flags; and it retains its state across episodes in a given Markov Decision Process (MDP). The activations of the RNN store the state of the "fast" RL algorithm on the current (previously unseen) MDP. We evaluate RL$^2$ experimentally on both small-scale and large-scale problems. On the small-scale side, we train it to solve randomly generated multi-arm bandit problems and finite MDPs. After RL$^2$ is trained, its performance on new MDPs is close to human-designed algorithms with optimality guarantees. On the large-scale side, we test RL$^2$ on a vision-based navigation task and show that it scales up to high-dimensional problems.

preprint2016arXiv

Scalable Bayesian modeling, monitoring and analysis of dynamic network flow data

Traffic flow count data in networks arise in many applications, such as automobile or aviation transportation, certain directed social network contexts, and Internet studies. Using an example of Internet browser traffic flow through site-segments of an international news website, we present Bayesian analyses of two linked classes of models which, in tandem, allow fast, scalable and interpretable Bayesian inference. We first develop flexible state-space models for streaming count data, able to adaptively characterize and quantify network dynamics efficiently in real-time. We then use these models as emulators of more structured, time-varying gravity models that allow formal dissection of network dynamics. This yields interpretable inferences on traffic flow characteristics, and on dynamics in interactions among network nodes. Bayesian monitoring theory defines a strategy for sequential model assessment and adaptation in cases when network flow data deviates from model-based predictions. Exploratory and sequential monitoring analyses of evolving traffic on a network of web site-segments in e-commerce demonstrate the utility of this coupled Bayesian emulation approach to analysis of streaming network count data.

preprint2016arXiv

Shortcut to adiabatic population transfer in quantum three-level systems: effective two-level problems and feasible counter-diabatic driving

Shortcut to adiabaticity in various quantum systems has attracted much attention with the wide applications in quantum information processing and quantum control. In this paper, we concentrate on stimulated Raman shortcut-to-adiabatic passage in quantum three-level systems. To implement counter-diabatic driving but without additional coupling, we first reduce the quantum three-level systems to effective two-level problems at large intermediate-level detuning, or on resonance, apply counter-diabatic driving along with the unitary transformation, and eventually modify the pump and Stokes pulses for achieving fast and high-fidelity population transfer. The required laser intensity and stability against parameter variation are further discussed, to demonstrate the advantage of shortcuts to adiabaticity.

preprint2016arXiv

The bandmerged Planck Early Release Compact Source Catalogue: Probing sub-structure in the molecular gas at high Galactic latitude

The Planck Early Release Compact Source Catalogue (ERCSC) includes nine lists of highly reliable sources, individually extracted at each of the nine Planck frequency channels. To facilitate the study of the Planck sources, especially their spectral behaviour across the radio/infrared frequencies, we provide a "bandmerged" catalogue of the ERCSC sources. This catalogue consists of 15191 entries, with 79 sources detected in all nine frequency channels of Planck and 6818 sources detected in only one channel. We describe the bandmerging algorithm, including the various steps used to disentangle sources in confused regions. The multi-frequency matching allows us to develop spectral energy distributions of sources between 30 and 857 GHz, in particular across the 100 GHz band, where the energetically important CO J=1->0 line enters the Planck bandpass. We find ~3-5sigma evidence for contribution to the 100 GHz intensity from foreground CO along the line of sight to 147 sources with |b|>30 deg. The median excess contribution is 4.5+/-0.9 percent of their measured 100 GHz flux density which cannot be explained by calibration or beam uncertainties. This translates to 0.5+/-0.1 K km s^{-1} of CO which must be clumped on the scale of the Planck 100 GHz beam, i.e., ~10 arcmin. If this is due to a population of low mass (~15 Msun) molecular gas clumps, the total mass in these clumps may be more than 2000 Msun. Further, high-spatial-resolution, ground-based observations of the high-latitude sky will help shed light on the origin of this diffuse, clumpy CO emission.

preprint2016arXiv

The Role of Context Selection in Object Detection

We investigate the reasons why context in object detection has limited utility by isolating and evaluating the predictive power of different context cues under ideal conditions in which context provided by an oracle. Based on this study, we propose a region-based context re-scoring method with dynamic context selection to remove noise and emphasize informative context. We introduce latent indicator variables to select (or ignore) potential contextual regions, and learn the selection strategy with latent-SVM. We conduct experiments to evaluate the performance of the proposed context selection method on the SUN RGB-D dataset. The method achieves a significant improvement in terms of mean average precision (mAP), compared with both appearance based detectors and a conventional context model without the selection scheme.

preprint2016arXiv

TMRT Observations of Carbon-chain molecules in Serpens South 1A

We report Shanghai Tian Ma Radio Telescope detections of several long carbon-chain molecules at C and Ku band, including HC3N, HC5N, HC7N, HC9N, C3S, C6H and C8H toward the starless cloud Serpens South 1a. We detected some transitions (HC9N J=13-12 F=12-11 and F=14-13, H13CCCN J=2-1 F=1-0 and F=1-1, HC13CCN J=2-1 F=2-2, F=1-0 and F=1-1, HCC13CN J=2-1 F=1-0 and F=1-1) and resolved some hyperfine components (HC5N J=6-5 F=5-4, H13CCCN J=2-1 F=2-1) for the first time in the interstellar medium. The column densities of these carbon-chain molecules in a range of 10^{12}-10^{13} cm^{-2} are comparable to two carbon-chain molecule rich sources, TMC-1 and Lupus-1A. The abundance ratios are 1.00:(1.11\pm0.15):(1.47\pm0.18) for [H13CCCN]:[HC13CCN]:[HCC13CN]. This result implies that the 13C isotope is also concentrated in the carbon atom adjacent to the nitrogen atom in HC3N in Serpens south 1a, which is similar to TMC-1. The [HC3N]/[H13CCCN] ratio of 78\pm9, the [HC3N]/[HC13CCN] ratio of 70\pm8, and the [HC3N]/[HCC13CN] ratio of 53\pm4 are also comparable to those in TMC-1. In any case, Serpens South 1a proves a testing ground for understanding carbon-chain chemistry.

preprint2016arXiv

Total positivity of recursive matrices

Let $A=[a_{n,k}]_{n,k\ge 0}$ be an infinite lower triangular matrix defined by the recurrence $$a_{0,0}=1,\quad a_{n+1,k}=r_{k}a_{n,k-1}+s_{k}a_{n,k}+t_{k+1}a_{n,k+1},$$ where $a_{n,k}=0$ unless $n\ge k\ge 0$ and $r_k,s_k,t_k$ are all nonnegative. Many well-known combinatorial triangles are such matrices, including the Pascal triangle, the Stirling triangle (of the second kind), the Bell triangle, the Catalan triangles of Aigner and Shapiro. We present some sufficient conditions such that the recursive matrix $A$ is totally positive. As applications we give the total positivity of the above mentioned combinatorial triangles in a unified approach.

preprint2016arXiv

Unsupervised Network Pretraining via Encoding Human Design

Over the years, computer vision researchers have spent an immense amount of effort on designing image features for the visual object recognition task. We propose to incorporate this valuable experience to guide the task of training deep neural networks. Our idea is to pretrain the network through the task of replicating the process of hand-designed feature extraction. By learning to replicate the process, the neural network integrates previous research knowledge and learns to model visual objects in a way similar to the hand-designed features. In the succeeding finetuning step, it further learns object-specific representations from labeled data and this boosts its classification power. We pretrain two convolutional neural networks where one replicates the process of histogram of oriented gradients feature extraction, and the other replicates the process of region covariance feature extraction. After finetuning, we achieve substantially better performance than the baseline methods.

preprint2016arXiv

Visualizing the elongated vortices in $γ$-Ga nanostrips

We study the magnetic response of superconducting $γ$-Ga via low temperature scanning tunneling microscopy and spectroscopy. The magnetic vortex cores rely substantially on the Ga geometry, and exhibit an unexpectedly-large axial elongation with aspect ratio up to 40 in rectangular Ga nano-strips (width $l$ $<$ 100 nm). This is in stark contrast with the isotropic circular vortex core in a larger round-shaped Ga island. We suggest that the unusual elongated vortices in Ga nanostrips originate from geometric confinement effect probably via the strong repulsive interaction between the vortices and Meissner screening currents at the sample edge. Our finding provides novel conceptual insights into the geometrical confinement effect on magnetic vortices and forms the basis for the technological applications of superconductors.

preprint2016arXiv

Weak coupling of pseudoacoustic phonons and magnon dynamics in incommensurate spin ladder compound Sr14Cu24O41

Intriguing lattice dynamics has been predicted for aperiodic crystals that contain incommensurate substructures. Here we report inelastic neutron scattering measurements of phonon and magnon dispersions in Sr14Cu24O41, which contains incommensurate one-dimensional (1D) chain and two-dimensional (2D) ladder substructures. Two distinct pseudoacoustic phonon modes, corresponding to the sliding motion of one sublattice against the other, are observed for atomic motions polarized along the incommensurate axis. In the long wavelength limit, it is found that the sliding mode shows a remarkably small energy gap of 1.7-1.9 meV, indicating very weak interactions between the two incommensurate sublattices. The measurements also reveal a gapped and steep linear magnon dispersion of the ladder sublattice. The high group velocity of this magnon branch and weak coupling with acoustic and pseudoacoustic phonons can explain the large magnon thermal conductivity in Sr14Cu24O41 crystals. In addition, the magnon specific heat is determined from the measured total specific heat and phonon density of states, and exhibits a Schottky anomaly due to gapped magnon modes of the spin chains. These findings offer new insights into the phonon and magnon dynamics and thermal transport properties of incommensurate magnetic crystals that contain low-dimensional substructures.

preprint2016arXiv

Zero-energy modes and valley asymmetry in the Hofstadter spectrum of bilayer graphene van der Waals heterostructures with hBN

We investigate the magnetic minibands of a heterostructure consisting of bilayer graphene (BLG) and hexagonal boron nitride (hBN) by numerically diagonalizing a two-band Hamiltonian that describes electrons in BLG in the presence of a moire potential. Due to inversion-symmetry breaking characteristic for the moire potential, the valley symmetry of the spectrum is broken, but despite this, the zero-energy Landau level in BLG survives, albeit with reduced degeneracy. In addition, we derive effective models for the low-energy features in the magnetic minibands and demonstrate the appearance of secondary Dirac points in the valence band, which we confirm by numerical simulations. Then, we analyze how single-particle gaps in the fractal energy spectrum produce a sequence of incompressible states observable under a variation of carrier density and magnetic field.

preprint2015arXiv

A Bayesian Particle Filtering Method For Brain Source Localisation

In this paper, we explore the multiple source localisation problem in the cerebral cortex using magnetoencephalography (MEG) data. We model neural currents as point-wise dipolar sources which dynamically evolve over time, then model dipole dynamics using a probabilistic state space model in which dipole locations are strictly constrained to lie within the cortex. Based on the proposed models, we develop a Bayesian particle filtering algorithm for localisation of both known and unknown numbers of dipoles. The algorithm consists of a region of interest (ROI) estimation step for initial dipole number estimation, a Gibbs multiple particle filter (GMPF) step for individual dipole state estimation, and a selection criterion step for selecting the final estimates. The estimated results from the ROI estimation are used to adaptively adjust particle filter's sample size to reduce the overall computational cost. The proposed models and the algorithm are tested in numerical experiments. Results are compared with existing particle filtering methods. The numerical results show that the proposed methods can achieve improved performance metrics in terms of dipole number estimation and dipole localisation.

preprint2015arXiv

A New Method for Triangular Mesh Generation

Computational mathematics plays an increasingly important role in computational fluid dynamics (CFD). The aeronautics and aerospace re- search community is working on next generation of CFD capacity that is accurate, automatic, and fast. A key component of the next generation of CFD is a greatly enhanced capacity for mesh generation and adaptivity of the mesh according to solution and geometry. In this paper, we propose a new method that generates triangular meshes on domains of curved boundary. The method deforms a Cartesian mesh that covers the domain to generate a mesh with prescribed boundary nodes. The deformation fields are generated by a system of divergence and curl equations which are solved effectively by the least square finite element method.

preprint2015arXiv

Addition is exponentially harder than counting for shallow monotone circuits

Let $U_{k,N}$ denote the Boolean function which takes as input $k$ strings of $N$ bits each, representing $k$ numbers $a^{(1)},\dots,a^{(k)}$ in $\{0,1,\dots,2^{N}-1\}$, and outputs 1 if and only if $a^{(1)} + \cdots + a^{(k)} \geq 2^N.$ Let THR$_{t,n}$ denote a monotone unweighted threshold gate, i.e., the Boolean function which takes as input a single string $x \in \{0,1\}^n$ and outputs $1$ if and only if $x_1 + \cdots + x_n \geq t$. We refer to circuits that are composed of THR gates as monotone majority circuits. The main result of this paper is an exponential lower bound on the size of bounded-depth monotone majority circuits that compute $U_{k,N}$. More precisely, we show that for any constant $d \geq 2$, any depth-$d$ monotone majority circuit computing $U_{d,N}$ must have size $\smash{2^{Ω(N^{1/d})}}$. Since $U_{k,N}$ can be computed by a single monotone weighted threshold gate (that uses exponentially large weights), our lower bound implies that constant-depth monotone majority circuits require exponential size to simulate monotone weighted threshold gates. This answers a question posed by Goldmann and Karpinski (STOC'93) and recently restated by Hastad (2010, 2014). We also show that our lower bound is essentially best possible, by constructing a depth-$d$, size-$2^{O(N^{1/d})}$ monotone majority circuit for $U_{d,N}$. As a corollary of our lower bound, we significantly strengthen a classical theorem in circuit complexity due to Ajtai and Gurevich (JACM'87). They exhibited a monotone function that is in AC$^0$ but requires super-polynomial size for any constant-depth monotone circuit composed of unbounded fan-in AND and OR gates. We describe a monotone function that is in depth-$3$ AC$^0$ but requires exponential size monotone circuits of any constant depth, even if the circuits are composed of THR gates.

preprint2015arXiv

Characterizing Information Spreading in Online Social Networks

Online social networks (OSNs) are changing the way in which the information spreads throughout the Internet. A deep understanding of the information spreading in OSNs leads to both social and commercial benefits. In this paper, we characterize the dynamic of information spreading (e.g., how fast and widely the information spreads against time) in OSNs by developing a general and accurate model based on the Interactive Markov Chains (IMCs) and mean-field theory. This model explicitly reveals the impacts of the network topology on information spreading in OSNs. Further, we extend our model to feature the time-varying user behaviors and the ever-changing information popularity. The complicated dynamic patterns of information spreading are captured by our model using six key parameters. Extensive tests based on Renren's dataset validate the accuracy of our model, which demonstrate that it can characterize the dynamic patterns of video sharing in Renren precisely and predict future spreading tendency successfully.

preprint2015arXiv

Collapse of spin-orbit coupled Bose-Einstein condensates

A finite-size quasi two-dimensional Bose-Einstein condensate collapses if the attraction between atoms is sufficiently strong. Here we present a theory of collapse for condensates with the interatomic attraction and spin-orbit coupling. We consider two realizations of spin-orbit coupling: the axial Rashba coupling and balanced, effectively one-dimensional, Rashba-Dresselhaus one. In both cases spin-dependent "anomalous" velocity, proportional to the spin-orbit coupling strength, plays a crucial role. For the Rashba coupling, this velocity forms a centrifugal component in the density flux opposite to that arising due to the attraction between particles and prevents the collapse at a sufficiently strong coupling. For the balanced Rashba-Dresselhaus coupling, the spin-dependent velocity can spatially split the initial state in one dimension and form spin-projected wavepackets, reducing the total condensate density. Depending on the spin-orbit coupling strength, interatomic attraction, and the initial state, this splitting either prevents the collapse or modifies the collapse process. These results show that the collapse can be controlled by a spin-orbit coupling, thus, extending the domain of existence of condensates of attracting atoms.

preprint2015arXiv

Detection of a methanol megamaser in a major-merger galaxy

We have detected emission from both the 4_{-1}-3_{0} E (36.2~GHz) class I and 7_{-2}-8_{-1} E (37.7~GHz) class II methanol transitions towards the centre of the closest ultra-luminous infrared galaxy Arp 220. The emission in both the methanol transitions show narrow spectral features and have luminosities approximately 8 orders of magnitude stronger than that observed from typical class I methanol masers observed in Galactic star formation regions. The emission is also orders of magnitude stronger than the expected intensity of thermal emission from these transitions and based on these findings we suggest that the emission from the two transitions are masers. These observations provides the first detection of a methanol megamaser in the 36.2 and 37.7 GHz transitions and represents only the second detection of a methanol megamaser, following the recent report of an 84 GHz methanol megamaser in NGC1068. We find the methanol megamasers are significantly offset from the nuclear region and arise towards regions where there is Ha emission, suggesting that it is associated with starburst activity. The high degree of correlation between the spatial distribution of the 36.2 GHz methanol and X-ray plume emission suggests that the production of strong extragalactic class I methanol masers is related to galactic outflow driven shocks and perhaps cosmic rays. In contrast to OH and H2O megamasers which originate close to the nucleus, methanol megamasers provide a new probe of feedback (e.g. outflows) processes on larger-scales and of star formation beyond the circumnuclear starburst regions of active galaxies.

preprint2015arXiv

Detection of a superconducting phase in a two-atom layer of hexagonal Ga film grown on semiconducting GaN(0001)

The recent observation of superconducting state at atomic scale has motivated the pursuit of exotic condensed phases in two-dimensional (2D) systems. Here we report on a superconducting phase in two-monolayer crystalline Ga films epitaxially grown on wide band-gap semiconductor GaN(0001). This phase exhibits a hexagonal structure and only 0.552 nm in thickness, nevertheless, brings about a superconducting transition temperature Tc as high as 5.4 K, confirmed by in situ scanning tunneling spectroscopy, and ex situ electrical magneto-transport and magnetization measurements. The anisotropy of critical magnetic field and Berezinski-Kosterlitz-Thouless-like transition are observed, typical for the 2D superconductivity. Our results demonstrate a novel platform for exploring atomic-scale 2D superconductor, with great potential for understanding of the interface superconductivity.

preprint2015arXiv

Disentangling the magnetoelectric and thermoelectric transport in topological insulator thin films

We report transport studies on (Bi,Sb)2Te3 topological insulator thin films with tunable electronic band structure. We find a doping and temperature regime in which the Hall coefficient is negative indicative of electron-type carriers, whereas the Seebeck coefficient is positive indicative of hole-type carriers. This sign anomaly is due to the distinct transport behaviors of the bulk and surface states: the surface Dirac fermions dominate magnetoelectric transport while the thermoelectric effect is mainly determined by the bulk states. These findings may inspire new ideas for designing topological insulator-based high efficiency thermoelectric devices.

preprint2015arXiv

Fast and optimal transport of atoms with non-harmonic traps

We investigate the fast transport of an atom or a packet of atoms by different kinds of non-harmonic traps including power-law traps. The study is based on the reverse engineering method. Exact results are obtained and applied to design robust transport protocols. The optimization of the transport trajectory is performed with classical trajectories, and remains valid for the transport of a wave packet.

preprint2015arXiv

Fast control of topological vortex formation in BEC by counter-diabatic driving

Topological vortex formation has been known as the simplest method for vortex formation in BEC of alkali atoms. This scheme requires inversion of the bias magnetic field along the axis of the condensate, which leads to atom loss when the bias field crosses zero. In this Letter, we propose a scheme with which the atom loss is greatly suppressed by adding counter-diabatic magnetic field. A naive counter-diabatic field violates the Maxwell equations and we need to introduce an approximation to make it physically feasible. The resulting field requires an extra currents, which is experimentally challenging. Finally we solve this problem by applying a gauge transformation so that the counter-diabatic field is generated by controlling the original trap field with the additional control of the bias field.

preprint2015arXiv

Hierarchy of gaps and magnetic minibands in graphene in the presence of the Abrikosov vortex lattice

We determine the structure of band and gaps in graphene encapsulated in hexagonal boron nitride and subjected to magnetic field of Abrikosov lattice of vortices in the underlying superconducting film. The spectrum features one non-dispersive magnetic miniband at zero energy, separated by the largest gaps in the miniband spectrum from a pair of minibands resembling slightly broadened first Landau levels in graphene, suggesting the persistence of $ν= \pm 2$ quantum Hall effect states. Also, we identify occasional merging point of magnetic minibands which feature Dirac-type dispersion at the consecutive miniband edges.

preprint2015arXiv

Interface enhanced electron-phonon coupling and high temperature superconductivity in potassium-coated ultra-thin FeSe films on SrTiO3

Alkali-metal (potassium) adsorption on FeSe thin films with thickness from two unit cells (UC) to 4-UC on SrTiO3 grown by molecular beam epitaxy is investigated with a low-temperature scanning tunneling microscope. At appropriate potassium coverage (0.2-0.3 monolayer), the tunneling spectra of the films all exhibit a superconducting-like gap larger than 11 meV (five times the gap value of bulk FeSe), and two distinct features of characteristic phonon modes at 11 meV and 21 meV. The results reveal the critical role of the interface enhanced electron-phonon coupling for possible high temperature superconductivity in the system and is consistent with recent theories. Our study provides compelling evidence for the conventional pairing mechanism for this type of heterostructure superconducting systems.

preprint2015arXiv

Molecular Beam Epitaxy Growth of Superconducting LiFeAs Film on SrTiO3(001) Substrate

The stoichiometric "111" iron-based superconductor, LiFeAs, has attacted great research interest in recent years. For the first time, we have successfully grown LiFeAs thin film by molecular beam epitaxy (MBE) on SrTiO3(001) substrate, and studied the interfacial growth behavior by reflection high energy electron diffraction (RHEED) and low-temperature scanning tunneling microscope (LT-STM). The effects of substrate temperature and Li/Fe flux ratio were investigated. Uniform LiFeAs film as thin as 3 quintuple-layer (QL) is formed. Superconducting gap appears in LiFeAs films thicker than 4 QL at 4.7 K. When the film is thicker than 13 QL, the superconducting gap determined by the distance between coherence peaks is about 7 meV, close to the value of bulk material. The ex situ transport measurement of thick LiFeAs film shows a sharp superconducting transition around 16 K. The upper critical field, Hc2(0)=13.0 T, is estimated from the temperature dependent magnetoresistance. The precise thickness and quality control of LiFeAs film paves the road of growing similar ultrathin iron arsenide films.

preprint2015arXiv

Near-optimal small-depth lower bounds for small distance connectivity

We show that any depth-$d$ circuit for determining whether an $n$-node graph has an $s$-to-$t$ path of length at most $k$ must have size $n^{Ω(k^{1/d}/d)}$. The previous best circuit size lower bounds for this problem were $n^{k^{\exp(-O(d))}}$ (due to Beame, Impagliazzo, and Pitassi [BIP98]) and $n^{Ω((\log k)/d)}$ (following from a recent formula size lower bound of Rossman [Ros14]). Our lower bound is quite close to optimal, since a simple construction gives depth-$d$ circuits of size $n^{O(k^{2/d})}$ for this problem (and strengthening our bound even to $n^{k^{Ω(1/d)}}$ would require proving that undirected connectivity is not in $\mathsf{NC^1}.$) Our proof is by reduction to a new lower bound on the size of small-depth circuits computing a skewed variant of the "Sipser functions" that have played an important role in classical circuit lower bounds [Sip83, Yao85, Hås86]. A key ingredient in our proof of the required lower bound for these Sipser-like functions is the use of \emph{random projections}, an extension of random restrictions which were recently employed in [RST15]. Random projections allow us to obtain sharper quantitative bounds while employing simpler arguments, both conceptually and technically, than in the previous works [Ajt89, BPU92, BIP98, Ros14].

preprint2015arXiv

New Classes of Permutation Binomials and Permutation Trinomials over Finite Fields

Permutation polynomials over finite fields play important roles in finite fields theory. They also have wide applications in many areas of science and engineering such as coding theory, cryptography, combinatorial design, communication theory and so on. Permutation binomials and trinomials attract people's interest due to their simple algebraic form and additional extraordinary properties. In this paper, several new classes of permutation binomials and permutation trinomials are constructed. Some of these permutation polynomials are generalizations of known ones.

preprint2015arXiv

New Variational Method of Grid Generation with Prescribed Jacobian Determinant and Prescribed Curl

Adaptive grid generation is an active research topic for numer- ical solution of differential equations. In this paper, we propose a variational method which generates transformations with prescribed Jacobian determinant and curl. Then we use this transformation to achieve adaptive grid generation task, and show the importance of curl in a transformation.

preprint2015arXiv

Optimal transport of two ions under slow spring-constant drifts

We investigate the effect of slow spring-constant drifts of the trap used to shuttle two ions of different mass. We design transport protocols to suppress or mitigate the final excitation energy by applying invariant-based inverse engineering, perturbation theory, and a harmonic dynamical normal-mode approximation. A simple, explicit trigonometric protocol for the trap trajectory is found to be robust with respect to the spring-constant drifts.

preprint2015arXiv

Probing Dirac Fermion Dynamics in Topological Insulator Bi$_2$Se$_3$ Films with Scanning Tunneling Microscope

Scanning tunneling microscopy and spectroscopy have been used to investigate the femtosecond dynamics of Dirac fermions in the topological insulator Bi$_2$Se$_3$ ultrathin films. At two-dimensional limit, bulk electrons becomes quantized and the quantization can be controlled by film thickness at single quintuple layer level. By studying the spatial decay of standing waves (quasiparticle interference patterns) off steps, we measure directly the energy and film thickness dependence of phase relaxation length $l_ϕ$ and inelastic scattering lifetime $τ$ of topological surface-state electrons. We find that $τ$ exhibits a remarkable $(E-E_F)^{-2}$ energy dependence and increases with film thickness. We show that the features revealed are typical for electron-electron scattering between surface and bulk states.

preprint2015arXiv

Pulse design without rotating wave approximation

We design realizable time-dependent semiclassical pulses to invert the population of a two-level system faster than adiabatically when the rotating-wave approximation cannot be applied. Different approaches, based on the counterdiabatic method or on invariants, may lead to singularities in the pulse functions. Ways to avoid or cancel the singularities are put forward when the pulse spans few oscillations. For many oscillations an alternative numerical minimization method is proposed and demonstrated.

preprint2015arXiv

Quantifying wall turbulence via a symmetry approach. Part I. A Lie group theory

First principle based prediction of mean flow quantities of wall-bounded turbulent flows (channel, pipe, and turbulent boundary layer - TBL) is of great importance from both physics and engineering standpoints. Here (Part I), we present a symmetry-based approach which derives analytic expressions governing the mean velocity profile (MVP) from an innovative Lie-group analysis. The new approach begins by identifying a set of order functions (e.g. stress length, shear-induced eddy length), in analogy with the order parameter in Landau's mean-field theory, which aims at capturing symmetry aspects of the fluctuations (e.g. Reynolds stress, dissipation). The order functions are assumed to satisfy a dilation group invariance - representing the effects of the wall on fluctuations - which allows us to postulate three new kinds of invariant solutions of the mean momentum equation (MME), focusing on group invariants of the order functions (rather than those of the mean velocity as done in previous studies). The first - a power law solution - gives functional forms for the viscous sublayer, the buffer layer, the log-layer, and a newly identified central `core' (for channel and pipe, but non-existent for TBL). The second - a defect power law of form $1-r^{m}$ ($r$ being the distance from the center line) - describes the `bulk zone' (the region of balance between production and dissipation). The third - a relation between the group invariants of the stress length function and its first derivative - describes scaling transition between adjacent layers. A combination of these three expressions yields a multi-layer formula covering the entire flow domain, identifying three important parameters: scaling exponent, layer thickness, and transition sharpness. All three kinds of invariant solutions are validated, individually and in combination, by data from direct numerical simulations (DNS).

preprint2015arXiv

Rationality of Euler-Chow series and finite generation of Cox rings

In this paper we work with a series whose coefficients are the Euler characteristic of Chow varieties of a given projective variety. For varieties where the Cox ring is defined, it is easy to see that in this case the ring associated to the series is the Cox ring. If this ring is noetherian then the series is rational. It is an open question whether the converse holds. In this paper we give an example showing the converse fails. However we conjecture that it holds when the variety is rationally connected. As an evidence of this conjecture, It is proved that the series is not rational, and in a sense defined, not algebraic, in the case of the blowup of the projective plane at nine or more points in general position. Furthermore, we also treat some other examples of varieties with infinitely generated Cox ring, studied by Mukai and Hassett-Tschinkel. These are the first examples known where the series is not rational. We also compute the series for Del Pezzo surfaces.

preprint2015arXiv

Resolvent and spectral measure on non-trapping asymptotically hyperbolic manifolds I: Resolvent construction at high energy

This is the first in a series of papers in which we investigate the resolvent and spectral measure on non-trapping asymptotically hyperbolic manifolds with applications to the restriction theorem, spectral multiplier results and Strichartz estimates. In this first paper, we use semiclassical Lagrangian distributions and semiclassical intersecting Lagrangian distributions, along with Mazzeo-Melrose 0-calculus, to construct the high energy resolvent on general non- trapping asymptotically hyperbolic manifolds, generalizing the work due to Melrose, Sa Barreto and Vasy. We note that there is an independent work by Y. Wang which also constructs the high-energy resolvent.

preprint2015arXiv

Resolvent and spectral measure on non-trapping asymptotically hyperbolic manifolds II: spectral measure, restriction theorem, spectral multiplier

We consider the Laplacian $Δ$ on an asymptotically hyperbolic manifold $X$, as defined by Mazzeo and Melrose. We give pointwise bounds on the spectral measure for the operator $(Δ- n^2/4)_+^{1/2}$ on such manifolds, under the assumptions that $X$ is nontrapping and there is no resonance at the bottom of the spectrum. This uses the construction of the resolvent given by Mazzeo and Melrose (valid when the spectral parameter lies in a compact set), Melrose, Sá Barreto and Vasy (high energy estimates for a perturbation of the hyperbolic metric) and the present authors (see also the work by Yiran Wang) in the general high-energy case. Under the additional assumption that $X$ has negative curvature everywhere, we give two applications of the spectral measure estimates. The first, following work due to Guillarmou and Sikora with the second author in the asymptotically conic case, is a restriction theorem, that is, a $L^p(X) \to L^{p'}(X)$ operator norm bound on the spectral measure. The second is a spectral multiplier result, that is, a bound on functions $F((Δ- n^2/4)_+^{1/2})$ of the square root of the Laplacian, in terms of norms of the function $F$. Compared to the asymptotically conic case, our spectral multiplier result is weaker, but the restriction estimate is stronger. In both cases, the difference can be traced to the exponential volume growth at infinity for asymptotically hyperbolic manifolds, as opposed to polynomial growth in the asymptotically conic setting. The pointwise bounds on the spectral measure established here will also be applied to Strichartz estimates in another work by the first author.

preprint2015arXiv

Shortcuts to adiabaticity

Quantum adiabatic processes -that keep constant the populations in the instantaneous eigenbasis of a time-dependent Hamiltonian- are very useful to prepare and manipulate states, but take typically a long time. This is often problematic because decoherence and noise may spoil the desired final state, or because some applications require many repetitions. "Shortcuts to adiabaticity" are alternative fast processes which reproduce the same final populations, or even the same final state, as the adiabatic process in a finite, shorter time. Since adiabatic processes are ubiquitous, the shortcuts span a broad range of applications in atomic, molecular and optical physics, such as fast transport of ions or neutral atoms, internal population control and state preparation (for nuclear magnetic resonance or quantum information), cold atom expansions and other manipulations, cooling cycles, wavepacket splitting, and many-body state engineering or correlations microscopy. Shortcuts are also relevant to clarify fundamental questions such as a precise quantification of the third principle of thermodynamics and quantum speed limits. We review different theoretical techniques proposed to engineer the shortcuts, the experimental results, and the prospects.

preprint2015arXiv

Shortcuts to adiabaticity for non-Hermitian systems

Adiabatic processes driven by non-Hermitian, time-dependent Hamiltonians may be sped up by generalizing inverse engineering techniques based on Berry's transitionless driving algorithm or on dynamical invariants. We work out the basic theory and examples described by two-level Hamiltonians: the acceleration of rapid adiabatic passage with a decaying excited level and of the dynamics of a classical particle on an expanding harmonic oscillator.

preprint2015arXiv

Simultaneous electrical-field-effect modulation of both top and bottom Dirac surface states of epitaxial thin films of three-dimensional topological insulators

It is crucial for the studies of the transport properties and quantum effects related to Dirac surface states of three-dimensional topological insulators (3D TIs) to be able to simultaneously tune the chemical potentials of both top and bottom surfaces of a 3D TI thin film. We have realized this in molecular beam epitaxy-grown thin films of 3D TIs, as well as magnetic 3D TIs, by fabricating dual-gate structures on them. The films could be tuned between n-type and p-type by each gate alone. Combined application of two gates can reduce the carrier density of a TI film to a much lower level than with only one of them and enhance the film resistance by 10000 %, implying that Fermi level is tuned very close to the Dirac points of both top and bottom surface states without crossing any bulk band. The result promises applications of 3D TIs in field effect devices.

preprint2015arXiv

Stein-Tomas Restriction Theorem via Spectral Measure on Metric Measure Spaces

The Stein-Tomas restriction theorem on Euclidean space says one can meaningfully restrict $\hat{f}$ to the unit sphere of $\mathbb{R}^n$ provided $f \in L^p(\mathbb{R}^n)$ with $1 < p < 2$. This result can be rewritten in terms of the estimates for the spectral measure of Laplacian. Guillarmou, Hassell and Sikora formulated a sufficient condition of the restriction theorem, via spectral measure, on abstract metric measure spaces. But they only proved the result in a special case. The present note aims to give a complete proof. In the end, it will be applied to the restriction theorem on asymptotically conic manifolds.

preprint2015arXiv

Superconducting Fluctuations in the Normal State of the Two-Dimensional Hubbard Model

We compute the two-particle quantities relevant for superconducting correlations in the two-dimensional Hubbard model within the dynamical cluster approximation. In the normal state we identify the parameter regime in density, interaction, and second-nearest-neighbor hopping strength that maximizes the $d_{x^2-y^2}$ superconducting transition temperature. We find in all cases that the optimal transition temperature occurs at intermediate coupling strength, and is suppressed at strong and weak interaction strengths. Similarly, superconducting fluctuations are strongest at intermediate doping and suppressed towards large doping and half-filling. We find a change in sign of the vertex contributions to $d_{xy}$ superconductivity from repulsive near half filling to attractive at large doping. $p$-wave superconductivity is not found at the parameters we study, and $s$-wave contributions are always repulsive. For negative second-nearest-neighbor hopping the optimal transition temperature shifts towards the electron-doped side in opposition to the van Hove singularity which moves towards hole doping. We surmise that an increase of the local interaction of the electron-doped compounds would increase $T_c$.

preprint2015arXiv

Superconductivity dichotomy in K-coated single and double unit cell FeSe films on SrTiO3

We report the superconductivity evolution of one unit cell (1-UC) and 2-UC FeSe films on SrTiO3(001) substrates with potassium (K) adsorption. By in situ scanning tunneling spectroscopy measurement, we find that the superconductivity in 1-UC FeSe films is continuously suppressed with increasing K coverage, whereas non-superconducting 2-UC FeSe films become superconducting with a gap of ~17 meV or ~11 meV depending on whether the underlying 1-UC films are superconducting or not. This work explicitly reveals that the interface electron-phonon coupling is strongly related to the charge transfer at FeSe/STO interface and plays vital role in enhancing Cooper pairing in both 1-UC and 2-UC FeSe films.

preprint2015arXiv

Test of Equivalence Principle at $10^{-8}$ Level by a Dual-species Double-diffraction Raman Atom Interferometer

We report an improved test of the weak equivalence principle by using a simultaneous $^{85}$Rb-$^{87}$Rb dual-species atom interferometer. We propose and implement a four-wave double-diffraction Raman transition scheme for the interferometer, and demonstrate its ability in suppressing common-mode phase noise of Raman lasers after their frequencies and intensity ratios are optimized. The statistical uncertainty of the experimental data for Eötvös parameter $η$ is $0.8\times10^{-8}$ at 3200 s. With various systematic errors corrected the final value is $η=(2.8\pm3.0)\times10^{-8}$. The major uncertainty is attributed to the Coriolis effect.

preprint2015arXiv

Tight Bounds for the Distribution-Free Testing of Monotone Conjunctions

We improve both upper and lower bounds for the distribution-free testing of monotone conjunctions. Given oracle access to an unknown Boolean function $f:\{0,1\}^n \rightarrow \{0,1\}$ and sampling oracle access to an unknown distribution $\mathcal{D}$ over $\{0,1\}^n$, we present an $\tilde{O}(n^{1/3}/ε^5)$-query algorithm that tests whether $f$ is a monotone conjunction versus $ε$-far from any monotone conjunction with respect to $\mathcal{D}$. This improves the previous best upper bound of $\tilde{O}(n^{1/2}/ε)$ by Dolev and Ron when $1/ε$ is small compared to $n$. For some constant $ε_0>0$, we also prove a lower bound of $\tildeΩ(n^{1/3})$ for the query complexity, improving the previous best lower bound of $\tildeΩ(n^{1/5})$ by Glasner and Servedio. Our upper and lower bounds are tight, up to a poly-logarithmic factor, when the distance parameter $ε$ is a constant. Furthermore, the same upper and lower bounds can be extended to the distribution-free testing of general conjunctions, and the lower bound can be extended to that of decision lists and linear threshold functions.

preprint2015arXiv

Time and spatial parity operations with trapped ions

We propose a physical implementation of time and spatial parity transformations, as well as Galilean boosts, in a trapped-ion quantum simulator. By embedding the simulated model into an enlarged simulating Hilbert space, these fundamental symmetry operations can be fully realized and measured with ion traps. We illustrate our proposal with analytical and numerical techniques of prototypical examples with state-of-the-art trapped-ion platforms. These results pave the way for the realization of time and spatial parity transformations in other models and quantum platforms.

preprint2015arXiv

Topological p+ip superconductivity in doped graphene-like single-sheet materials BC$_3$

We theoretically study exotic superconducting phases in graphenelike single-sheet material BC$_3$ doped to its type-II Van Hove singularity whose saddle point momenta are {\it not} time-reversal-invariant. From combined renormalization group analysis and RPA calculations, we show that the dominant superconducting instability induced by weak repulsive interactions is in the time-reversal-invariant $p+ip$ pairing channel because of the interplay among dominant ferromagnetic fluctuations, subleading spin fluctuations at finite momentum, and spin-orbit coupling. Such time-reversal-invariant $p+ip$ superconductivity has nontrivial $\mathbb{Z}_2$ topological invariant. Our results show that doped BC$_3$ provides a promising route to realize a genuine 2D helical $p+ip$ superconductor.

preprint2015arXiv

Transient particle energies in shortcuts to adiabatic expansions of harmonic traps

The expansion of a harmonic potential that holds a quantum particle may be realized without any final particle excitation but much faster than adiabatically via "shortcuts to adiabaticity" (STA). While ideally the process time can be reduced to zero, practical limitations and constraints impose minimal finite times for the externally controlled time-dependent frequency protocols. We examine the role of different time-averaged energies (total, kinetic, potential, non-adiabatic) and of the instantaneous power in characterizing or selecting different protocols.Specifically, we prove a virial theorem for STA processes, set minimal energies for specific times or viceversa, and discuss their realizability by means of Dirac impulses or otherwise.

preprint2015arXiv

Two Information Entropies of Square-Lattice $\pm {J}$ Ising Models

Two information entropies, $S_{IB}$ based on the local bonding energy configurations and $S_{IS}$ based on the local spin configurations, are applied in the square-lattice $\pm {J}$ Ising systems with Monte Carlo simulations. With $S_{IB}$ and $S_{IS}$, the spin glass states can be distinguished from paramagnetic states clearly. The results reveal that the nature of spin glass states is of ordering in bonding energy and disordering in spin orientations. The consistence of $S_{IB}$ and the thermodynamic entropy is found.

preprint2015arXiv

Vibrational mode multiplexing of ultracold atoms

Sending multiple messages on qubits encoded in different vibrational modes of cold atoms or ions along a transmission waveguide requires to merge first and then separate the modes at input and output ends. Similarly, different qubits can be stored in the modes of a trap and be separated later. We design the fast splitting of a harmonic trap into an asymmetric double well so that the initial ground vibrational state becomes the ground state of one of two final wells, and the initial first excited state becomes the ground state of the other well. This might be done adiabatically by slowly deforming the trap. We speed up the process by inverse engineering a double-function trap using dynamical invariants. The separation (demultiplexing) followed by an inversion of the asymmetric bias and then by the reverse process (multiplexing) provides a population inversion protocol based solely on trap reshaping.

preprint2015arXiv

Well-Supported versus Approximate Nash Equilibria: Query Complexity of Large Games

We study the randomized query complexity of approximate Nash equilibria (ANE) in large games. We prove that, for some constant $ε>0$, any randomized oracle algorithm that computes an $ε$-ANE in a binary-action, $n$-player game must make $2^{Ω(n/\log n)}$ payoff queries. For the stronger solution concept of well-supported Nash equilibria (WSNE), Babichenko previously gave an exponential $2^{Ω(n)}$ lower bound for the randomized query complexity of $ε$-WSNE, for some constant $ε>0$; the same lower bound was shown to hold for $ε$-ANE, but only when $ε=O(1/n)$. Our result answers an open problem posed by Hart and Nisan and Babichenko and is very close to the trivial upper bound of $2^n$. Our proof relies on a generic reduction from the problem of finding an $ε$-WSNE to the problem of finding an $ε/(4α)$-ANE, in large games with $α$ actions, which might be of independent interest.

preprint2015arXiv

What is the mass of alpha Cen B b?

We investigate the possibility of constraining the sin i degeneracy of alpha Cen B b -- with orbital period P=3.24 d; a = 0.042 AU; msini = 1.1 M_earth -- to estimate the true mass of the newly reported terrestrial exoplanet in the nearest stellar system to our Sun. We present detailed numerical simulations of the dynamical stability of the exoplanet in the alpha Cen AB binary system for a range of initial inclinations, eccentricities, and semi-major axes. The system represents a benchmark case for the interplay of the Kozai mechanism, general relativistic and tidal forces. From our simulations, there is only a small boundary in initial inclinations and initial semi-major axes that result in the migration via the Kozai mechanism of alpha Cen B b to its present location. Inside this boundary, the planet orbit is stable for up to 1 Gyr against the Kozai mechanism, and outside this boundary the planet collides with alpha Cen B or is ejected. In our three simulations where the planet migrates in towards the star via the Kozai mechanism, the final inclination is 46-53 degrees relative to the AB orbital plane, lower than the initial inclination of 75 degrees in each case. We discuss inclination constraints from the formation of alpha Cen B b in situ at its present location, migration in a proto-planetary disk, or migration in resonance with additional planets. We conclude that alpha Cen B b probably has a mass of less than 2.7 M_earth, implying a likely terrestrial composition warranting future confirmation.

preprint2014arXiv

A Scale-Free Topology Construction Model for Wireless Sensor Networks

A local-area and energy-efficient (LAEE) evolution model for wireless sensor networks is proposed. The process of topology evolution is divided into two phases. In the first phase, nodes are distributed randomly in a fixed region. In the second phase, according to the spatial structure of wireless sensor networks, topology evolution starts from the sink, grows with an energy-efficient preferential attachment rule in the new node's local-area, and stops until all nodes are connected into network. Both analysis and simulation results show that the degree distribution of LAEE follows the power law. This topology construction model has better tolerance against energy depletion or random failure than other non-scale-free WSN topologies.

preprint2014arXiv

Boolean function monotonicity testing requires (almost) $n^{1/2}$ non-adaptive queries

We prove a lower bound of $Ω(n^{1/2 - c})$, for all $c>0$, on the query complexity of (two-sided error) non-adaptive algorithms for testing whether an $n$-variable Boolean function is monotone versus constant-far from monotone. This improves a $\tildeΩ(n^{1/5})$ lower bound for the same problem that was recently given in [CST14] and is very close to $Ω(n^{1/2})$, which we conjecture is the optimal lower bound for this model.

preprint2014arXiv

Charge-Stripe Order in a Parent Compound of Iron-based Superconductors

Charge ordering is one of the most intriguing and extensively studied phenomena in correlated electronic materials because of its strong impact on electron transport properties including superconductivity. Despite its ubiquitousness in correlated systems, the occurrence of charge ordering in iron-based superconductors is still unresolved. Here we use scanning tunneling microscopy to reveal a long-range charge-stripe order and a highly anisotropic dispersion of electronic states in the ground state of stoichiometric FeTe, the parent compound of the Fe(Te, Se, S) superconductor family. The formation of charge order in a strongly correlated electron system with integer nominal valence (here Fe$^{2+}$) is unexpected and suggests that the iron-based superconductors may exhibit more complex charge dynamics than originally expected. We show that the present observations can be attributed to the surpassing of the role of local Coulomb interaction by the poorly screened longer-range Coulomb interactions, facilitated by large Hund's rule coupling.

preprint2014arXiv

Chemical potential dependent gap-opening at the Dirac surface states of Bi2Se3 induced by aggregated substitutional Cr atoms

With angle-resolved photoemission spectroscopy, gap-opening is resolved at up to room temperature in the Dirac surface states of molecular beam epitaxy grown Cr-doped Bi2Se3 topological insulator films, which however show no long-range ferromagnetic order down to 1.5 K. The gap size is found decreasing with increasing electron doping level. Scanning tunneling microscopy and first principles calculations demonstrate that substitutional Cr atoms aggregate into superparamagnetic multimers in Bi2Se3 matrix, which contribute to the observed chemical potential dependent gap-opening in the Dirac surface states without long-range ferromagnetic order.

preprint2014arXiv

Classifying and Visualizing Motion Capture Sequences using Deep Neural Networks

The gesture recognition using motion capture data and depth sensors has recently drawn more attention in vision recognition. Currently most systems only classify dataset with a couple of dozens different actions. Moreover, feature extraction from the data is often computational complex. In this paper, we propose a novel system to recognize the actions from skeleton data with simple, but effective, features using deep neural networks. Features are extracted for each frame based on the relative positions of joints (PO), temporal differences (TD), and normalized trajectories of motion (NT). Given these features a hybrid multi-layer perceptron is trained, which simultaneously classifies and reconstructs input data. We use deep autoencoder to visualize learnt features, and the experiments show that deep neural networks can capture more discriminative information than, for instance, principal component analysis can. We test our system on a public database with 65 classes and more than 2,000 motion sequences. We obtain an accuracy above 95% which is, to our knowledge, the state of the art result for such a large dataset.

preprint2014arXiv

Detection of 36 GHz class I methanol maser emission towards NGCS253

We have used the Australia Telescope Compact Array (ATCA) to search for emission from the $4_{-1} \rightarrow 3_{0}E$ transition of methanol (36.2 GHz) towards the center of the nearby starburst galaxy NGC253. Two regions of emission were detected, offset from the nucleus along the same position angle as the inner spiral arms. The emission is largely unresolved on a scale of 5 arcsec, has a full-width half maximum (FWHM) line width of < 30 km s$^{-1}$, and an isotropic luminosity orders of magnitude larger than that observed in any Galactic star formation regions. These characteristics suggest that the 36.2 GHz methanol emission is most likely a maser, although observations with higher angular and spectral resolution are required to confirm this. If it is a maser this represents the first detection of a class I methanol maser outside the Milky Way. The 36.2 GHz methanol emission in NGC253 has more than an order of magnitude higher isotropic luminosity than the widespread emission recently detected towards the center of the Milky Way. If emission from this transition scales with nuclear star formation rate then it may be detectable in the central regions of many starburst galaxies. Detection of methanol emission in ultra-luminous infra-red galaxies (ULIRGs) would open up a new tool for testing for variations in fundamental constants (in particular the proton-to-electron mass ratio) on cosmological scales.

preprint2014arXiv

Dichotomy of Electronic Structure and Superconductivity between Single-Layer and Double-Layer FeSe/SrTiO3 Films

The latest discovery of possible high temperature superconductivity in the single-layer FeSe film grown on a SrTiO3 substrate, together with the observation of its unique electronic structure and nodeless superconducting gap, has generated much attention. Initial work also found that, while the single-layer FeSe/SrTiO3 film exhibits a clear signature of superconductivity, the double-layer FeSe/SrTiO3 film shows an insulating behavior. Such a dramatic difference between the single-layer and double-layer FeSe/SrTiO3 films is surprising and the underlying origin remains unclear. Here we report our comparative study between the single-layer and double-layer FeSe/SrTiO3 films by performing a systematic angle-resolved photoemission study on the samples annealed in vacuum. We find that, like the single-layer FeSe/SrTiO3 film, the as-prepared double-layer FeSe/SrTiO3 film is insulating and possibly magnetic, thus establishing a universal existence of the magnetic phase in the FeSe/SrTiO3 films. In particular, the double-layer FeSe/SrTiO3 film shows a quite different doping behavior from the single-layer film in that it is hard to get doped and remains in the insulating state under an extensive annealing condition. The difference originates from the much reduced doping efficiency in the bottom FeSe layer of the double-layer FeSe/SrTiO3 film from the FeSe-SrTiO3 interface. These observations provide key insights in understanding the origin of superconductivity and the doping mechanism in the FeSe/SrTiO3 films. The property disparity between the single-layer and double-layer FeSe/SrTiO3 films may facilitate to fabricate electronic devices by making superconducting and insulating components on the same substrate under the same condition.

preprint2014arXiv

Dirac edges of fractal magnetic minibands in graphene with hexagonal moire superlattices

We find a systematic reappearance of massive Dirac features at the edges of consecutive minibands formed at magnetic fields B_{p/q}= pϕ_0/(qS) providing rational magnetic flux through a unit cell of the moire superlattice created by a hexagonal substrate for electrons in graphene. The Dirac-type features in the minibands at B=B_{p/q} determine a hierarchy of gaps in the surrounding fractal spectrum, and show that these minibands have topological insulator properties. Using the additional $q$-fold degeneracy of magnetic minibands at B_{p/q}, we trace the hierarchy of the gaps to their manifestation in the form of incompressible states upon variation of the carrier density and magnetic field.

preprint2014arXiv

Early-stage star forming cloud cores in GLIMPSE Extended Green Objects (EGOs) as traced by organic species

In order to investigate the physical and chemical properties of massive star forming cores in early stages, we analyse the excitation and abundance of four organic species, CH3OH, CH3OCH3, HCOOCH3 and CH3CH2CN, toward 29 Extended Green Object (EGO) cloud cores that were observed by our previous single dish spectral line survey. The EGO cloud cores are found to have similar methanol J_3-J_2 rotation temperatures of ~44 K, a typical linear size of ~0.036 pc, and a typical beam averaged methanol abundance of several 10^(-9) (the beam corrected value could reach several 10^(-7)). The abundances of the latter three species, normalized by that of methanol, are found to be correlated also across a large variety of clouds such as EGO cloud cores, hot corinos, massive hot cores and Galactic Center clouds. The chemical properties of the EGO cloud cores lie between that of hot cores and hot corinos. However, the abundances and abundance ratios of the four species can not be satisfactorily explained by recent chemical models either among the EGO cloud cores or among the various types of cloud cores from literature.

preprint2014arXiv

Electrically tuned magnetic order and magnetoresistance in a topological insulator

The Dirac-like surface states of the topological insulators (TIs) are protected by time reversal symmetry (TRS) and exhibit a host of novel properties. Introducing magnetism into TI, which breaks the TRS, is expected to create exotic topological magnetoelectric effects. A particularly intriguing phenomenon in this case is the magnetic field dependence of electrical resistance, or magnetoresistance (MR). The intricate interplay between topological protection and broken-TRS may lead to highly unconventional MR behaviour that can find unique applications in magnetic sensing and data storage. However, so far the MR of TI with spontaneously broken TRS is still poorly understood, mainly due to the lack of well-controlled experiments. In this work, we investigate the magneto transport properties of a ferromagnetic TI thin film fabricated into a field effect transistor device. We observe an unusually complex evolution of MR when the Fermi level (EF) is tuned across the Dirac point (DP) by gate voltage. In particular, MR tends to be positive when EF lies close to the DP but becomes negative at higher energies. This trend is opposite to that expected from the Berry phase picture for localization, but is intimately correlated with the gate-tuned magnetic order. We show that the underlying physics is the competition between the topology-induced weak antilocalization and magnetism-induced negative MR. The simultaneous electrical control of magnetic order and magneto transport facilitates future TI-based spintronic devices.

preprint2014arXiv

Electronic Evidence of an Insulator-Superconductor Transition in Single-Layer FeSe/SrTiO3 Films

In high temperature cuprate superconductors, it is now generally agreed that the parent compound is a Mott insulator and superconductivity is realized by doping the antiferromagnetic Mott insulator. In the iron-based superconductors, however, the parent compound is mostly antiferromagnetic metal, raising a debate on whether an appropriate starting point should go with an itinerant picture or a localized picture. It has been proposed theoretically that the parent compound of the iron-based superconductors may be on the verge of a Mott insulator, but so far no clear experimental evidence of doping-induced Mott transition has been available. Here we report an electronic evidence of an insulator-superconductor transition observed in the single-layer FeSe films grown on the SrTiO3 substrate. By taking angle-resolved photoemission measurements on the electronic structure and energy gap, we have identified a clear evolution of an insulator to a superconductor with the increasing doping. This observation represents the first example of an insulator-superconductor transition via doping observed in the iron-based superconductors. It indicates that the parent compound of the iron-based superconductors is in proximity of a Mott insulator and strong electron correlation should be considered in describing the iron-based superconductors.

preprint2014arXiv

Electronic Raman Scattering On Individual Single Walled Carbon Nanotubes

We report experimental measurements of electronic Raman scattering under resonant conditions by electrons in individual single-walled carbon nanotubes (SWNTs). The inelastic Raman scattering at low frequency range reveals a single particle excitation feature and the dispersion of electronic structure around the center of Brillouin zone of a semiconducting SWNT (14, 13) is extracted.

preprint2014arXiv

Exciton Binding Energy of Monolayer WS2

The optical properties of monolayer transition metal dichalcogenides (TMDC) feature prominent excitonic natures. Here we report an experimental approach toward measuring the exciton binding energy of monolayer WS2 with linear differential transmission spectroscopy and two-photon photoluminescence excitation spectroscopy (TP-PLE). TP-PLE measurements show the exciton binding energy of 0.71eV around K valley in the Brillouin zone. The trion binding energy of 34meV, two-photon absorption cross section 4X10^{4}cm^{2}W^{-2}S^{-1} at 780nm and exciton-exciton annihilation rate around 0.5cm^{2}/s are experimentally obtained.

preprint2014arXiv

Exclusive quantum channels in quantum networks

Quantum state can be teleported to a remote site by only local measurement and classical communication if the prior Einstein-Podolsky-Rosen quantum channel is available between the sender and the receiver. Those quantum channels shared by multiple nodes can constitute a quantum network. Yet, studies on the efficiency of quantum communication between nodes of quantum networks remain limited, which differs from classical case in that the quantum channel will be consumed if teleportation is performed. Here, we introduce the exclusive quantum channels (EQC) as the measure of efficiency of quantum information transmission. It quantifies the amount of quantum information which can be teleported between nodes in a quantum network. We show that different types of EQC are local quantities with effective circles. Significantly, capacity of quantum communication of quantum networks quantified by EQC is independent of distance for the communicating nodes. Thus, the quantum network can be dealt as the isotropic medium where quantum communication is no-decaying. EQC are studied by both analytical and numerical methods. The EQC can be enhanced by transformations of lattices of quantum network via entanglement swapping. Our result opens the avenue in studying the quantum communication of the quantum networks.

preprint2014arXiv

Experimental observation of Dirac-like surface states and topological phase transition in Pb$_{1-x}$Sn$_x$Te(111) films

The surface of a topological crystalline insulator (TCI) carries an even number of Dirac cones protected by crystalline symmetry. We epitaxially grew high quality Pb$_{1-x}$Sn$_x$Te(111) films and investigated the TCI phase by in-situ angle-resolved photoemission spectroscopy. Pb$_{1-x}$Sn$_x$Te(111) films undergo a topological phase transition from trivial insulator to TCI via increasing the Sn/Pb ratio, accompanied by a crossover from n-type to p-type doping. In addition, a hybridization gap is opened in the surface states when the thickness of film is reduced to the two-dimensional limit. The work demonstrates an approach to manipulating the topological properties of TCI, which is of importance for future fundamental research and applications based on TCI.

preprint2014arXiv

Fast shuttling of a trapped ion in the presence of noise

We theoretically investigate the motional excitation of a single ion caused by spring-constant and position uctuations of a harmonic trap during trap shuttling processes. A detailed study of the sensitivity on noise for several transport protocols and noise spectra is provided. The effect of slow spring-constant drifts is also analyzed. Trap trajectories that minimize the excitation are designed combining invariant-based inverse engineering, perturbation theory, and optimal control.

preprint2014arXiv

Fast transitionless expansions of Gaussian anharmonic traps for cold atoms: bang-singular-bang control

Combining invariant-based inverse engineering, perturbation theory, and Optimal Control Theory, we design fast, transitionless expansions of cold neutral atoms or ions in Gaussian anharmonic traps. Bounding the possible trap frequencies and using a "bang-singular-bang" control we find fast processes for a continuum of durations up to a minimum time that corresponds to a purely bang-bang (stepwise frequency constant) control.

preprint2014arXiv

HERMES: A Hierarchical Broadcast-Based Silicon Photonic Interconnect for Scalable Many-Core Systems

Optical interconnection networks, as enabled by recent advances in silicon photonic device and fabrication technology, have the potential to address on-chip and off-chip communication bottlenecks in many-core systems. Although several designs have shown superior power efficiency and performance compared to electrical alternatives, these networks will not scale to the thousands of cores required in the future. In this paper, we introduce Hermes, a hybrid network composed of an optimized broadcast for power-efficient low-latency global-scale coordination and circuit-switch sub-networks for high-throughput data delivery. This network will scale for use in thousand core chip systems. At the physical level, SoI-based adiabatic coupler has been designed to provide low-loss and compact optical power splitting. Based on the adiabatic coupler, a topology based on 2-ary folded butterfly is designed to provide linear power division in a thousand core layout with minimal cross-overs. To address the network agility and provide for efficient use of optical bandwidth, a flow control and routing mechanism is introduced to dynamically allocate bandwidth and provide fairness usage of network resources. At the system level, bloom filter-based filtering for localization of communication are designed for reducing global traffic. In addition, a novel greedy-based data and workload migration are leveraged to increase the locality of communication in a NUCA (non-uniform cache access) architecture. First order analytic evaluation results have indicated that Hermes is scalable to at least 1024 cores and offers significant performance improvement and power savings over prior silicon photonic designs.

preprint2014arXiv

Imaging the Electron-Boson Coupling in Superconducting FeSe

Scanning tunneling spectroscopy has been used to reveal signatures of a bosonic mode in the local quasiparticle density of states of superconducting FeSe films. The mode appears below Tc as a 'dip-hump' feature at energy Ω~ 4.7 kBTc beyond the superconducting gap Δ. Spectra on strained regions of the FeSe films reveal simultaneous decreases in Δand Ω. This contrasts with all previous reports on other high-Tc superconductors, where Δlocally anti-correlates with Ω. A local strong coupling model is found to reconcile the discrepancy well, and to provide a unified picture of the electron-boson coupling in unconventional superconductors.

preprint2014arXiv

New algorithms and lower bounds for monotonicity testing

We consider the problem of testing whether an unknown Boolean function $f$ is monotone versus $ε$-far from every monotone function. The two main results of this paper are a new lower bound and a new algorithm for this well-studied problem. Lower bound: We prove an $\tildeΩ(n^{1/5})$ lower bound on the query complexity of any non-adaptive two-sided error algorithm for testing whether an unknown Boolean function $f$ is monotone versus constant-far from monotone. This gives an exponential improvement on the previous lower bound of $Ω(\log n)$ due to Fischer et al. [FLN+02]. We show that the same lower bound holds for monotonicity testing of Boolean-valued functions over hypergrid domains $\{1,\ldots,m\}^n$ for all $m\ge 2$. Upper bound: We give an $\tilde{O}(n^{5/6})\text{poly}(1/ε)$-query algorithm that tests whether an unknown Boolean function $f$ is monotone versus $ε$-far from monotone. Our algorithm, which is non-adaptive and makes one-sided error, is a modified version of the algorithm of Chakrabarty and Seshadhri [CS13a], which makes $\tilde{O}(n^{7/8})\text{poly}(1/ε)$ queries.

preprint2014arXiv

Power-law scaling of asymmetries in viscous flow with PT-symmetric inflow and outflow

In recent years, open systems with balanced loss and gain, that are invariant under the combined parity and time-reversal ($\mathcal{PT}$) operations, have been studied via asymmetries of their solutions. They represent systems as diverse as coupled optical waveguides and electrical or mechanical oscillators. We numerically investigate the asymmetries of incompressible viscous flow in two and three dimensions with "balanced" inflow-outflow ($\mathcal{PT}$-symmetric) configurations. By introducing configuration-dependent classes of asymmetry functions in velocity, kinetic energy density, and vorticity fields, we find that the flow asymmetries exhibit power-law scaling with a single exponent in the laminar regime with the Reynolds number ranging over four decades. We show that such single-exponent scaling is expected for small Reynolds numbers, although its robustness at large values of Reynolds numbers is unexpected. Our results imply that $\mathcal{PT}$-symmetric inflow-outflow configurations provide a hitherto unexplored avenue to tune flow properties.

preprint2014arXiv

Prediction of temperature distribution in turbulent Rayleigh-Benard convection

A quantitative theory is developed for the vertical mean temperature profile (MTP) in turbulent Rayleigh-Benard convection (RBC), which explains the recent experimental and numerical observations of a logarithmic law by Ahlers et al.(Phys. Rev. Lett., 2012). A multi-layer model is formulated and quantified, whose predictions agree with DNS and experimental data for the Rayleigh-number (Ra) over seven decades. In particular, a thermal buffer layer follows a 1/7 scaling like the previously postulated mixing zone (Procaccia et al, Phys. Rev. A,1991), and yields a Ra-dependent log law constant. A new parameterization of Nu(Ra) dependence is proposed, based on the present multi-layer quantification of the bulk MTP.

preprint2014arXiv

Shortcuts to adiabaticity in three-level systems using Lie transforms

Sped-up protocols (shortcuts to adiabaticity) that drive a system quickly to the same populations than a slow adiabatic process may involve Hamiltonian terms difficult to realize in practice. We use the dynamical symmetry of the Hamiltonian to find, by means of Lie transforms, alternative Hamiltonians that achieve the same goals without the problematic terms. We apply this technique to three-level systems (two interacting bosons in a double well, and beam splitters with two and three output channels) driven by Hamiltonians that belong to the four-dimensional algebra U3S3.

preprint2014arXiv

Spectral Methods meet EM: A Provably Optimal Algorithm for Crowdsourcing

Crowdsourcing is a popular paradigm for effectively collecting labels at low cost. The Dawid-Skene estimator has been widely used for inferring the true labels from the noisy labels provided by non-expert crowdsourcing workers. However, since the estimator maximizes a non-convex log-likelihood function, it is hard to theoretically justify its performance. In this paper, we propose a two-stage efficient algorithm for multi-class crowd labeling problems. The first stage uses the spectral method to obtain an initial estimate of parameters. Then the second stage refines the estimation by optimizing the objective function of the Dawid-Skene estimator via the EM algorithm. We show that our algorithm achieves the optimal convergence rate up to a logarithmic factor. We conduct extensive experiments on synthetic and real datasets. Experimental results demonstrate that the proposed algorithm is comparable to the most accurate empirical approach, while outperforming several other recently proposed methods.

preprint2014arXiv

Spin-polarized currents generated by magnetic Fe atomic chains

Fe-based devices are widely used in spintronics because of high spin-polarization and magnetism. In this work, free-standing Fe atomic chains were proposed to be used as the thinnest wires to generate spin-polarized currents due to the spin-polarized energy bands. By ab initio calculations, the zigzag structure was found more stable than the wide-angle zigzag structure and has higher ratio of spin-up and spin-down currents. By our theoretical prediction, Fe atomic chains have sufficiently long thermal lifetime only at T<=150 K, while C atomic chains are very stable even at T=1000 K. This result means that the spintronic devices based on Fe chains could only work at low temperature. A system constructed by a short Fe chain sandwiched between two graphene electrodes was proposed as a spin-polarized current generator, while a C chain does not have such property. The present work may be instructive and meaningful to further practical applications based on recent technical development on the preparation of metal atomic chains [Proc. Natl. Acad. Sci. U.S.A. 107, 9055 (2010)].

preprint2014arXiv

Statistical Decision Making for Optimal Budget Allocation in Crowd Labeling

In crowd labeling, a large amount of unlabeled data instances are outsourced to a crowd of workers. Workers will be paid for each label they provide, but the labeling requester usually has only a limited amount of the budget. Since data instances have different levels of labeling difficulty and workers have different reliability, it is desirable to have an optimal policy to allocate the budget among all instance-worker pairs such that the overall labeling accuracy is maximized. We consider categorical labeling tasks and formulate the budget allocation problem as a Bayesian Markov decision process (MDP), which simultaneously conducts learning and decision making. Using the dynamic programming (DP) recurrence, one can obtain the optimal allocation policy. However, DP quickly becomes computationally intractable when the size of the problem increases. To solve this challenge, we propose a computationally efficient approximate policy, called optimistic knowledge gradient policy. Our MDP is a quite general framework, which applies to both pull crowdsourcing marketplaces with homogeneous workers and push marketplaces with heterogeneous workers. It can also incorporate the contextual information of instances when they are available. The experiments on both simulated and real data show that the proposed policy achieves a higher labeling accuracy than other existing policies at the same budget level.

preprint2014arXiv

The catalogues and mid-infrared environment of Interstellar OH Masers

Data for a number of OH maser lines have been collected from surveys. The posi- tions are compared to recent mid-infrared (MIR) surveys such as Spitzer-GLIMPSE and WISE, restricting the comparison to point sources. The colors and intensities of the IR sources are compared. There are many 18 cm OH masers, but far fewer in lines arising from higher energy levels. We also make a comparison with the 5 cm Class II methanol masers. We have divided the results into 3 subsamples: those associated with OH masers only, those associated with OH masers and Class II methanol masers, and those only associated with Class II methanol masers. There are no obvious dif- ferences in the color-color or color-magnitude results for the GLIMPSE point sources. However, according to the results from the WISE 22 μm survey, the sources associ- ated with OH masers are brighter than those associated with methanol masers. We interpret the presence of OH and methanol masers mark the locations of regions where stars are forming. The OH masers are located on the borders of sharp features found in the IR. These are referred to as bubbles. If the OH masers mark the positions of protostars, the result provides indirect evidence for triggered star formation caused by the expansion of the bubbles.

preprint2014arXiv

The complexity of approximating conservative counting CSPs

We study the complexity of approximately solving the weighted counting constraint satisfaction problem #CSP(F). In the conservative case, where F contains all unary functions, there is a classification known for the case in which the domain of functions in F is Boolean. In this paper, we give a classification for the more general problem where functions in F have an arbitrary finite domain. We define the notions of weak log-modularity and weak log-supermodularity. We show that if F is weakly log-modular, then #CSP(F)is in FP. Otherwise, it is at least as difficult to approximate as #BIS, the problem of counting independent sets in bipartite graphs. #BIS is complete with respect to approximation-preserving reductions for a logically-defined complexity class #RHPi1, and is believed to be intractable. We further sub-divide the #BIS-hard case. If F is weakly log-supermodular, then we show that #CSP(F) is as easy as a (Boolean) log-supermodular weighted #CSP. Otherwise, we show that it is NP-hard to approximate. Finally, we give a full trichotomy for the arity-2 case, where #CSP(F) is in FP, or is #BIS-equivalent, or is equivalent in difficulty to #SAT, the problem of approximately counting the satisfying assignments of a Boolean formula in conjunctive normal form. We also discuss the algorithmic aspects of our classification.

preprint2014arXiv

Ultrafast dynamics and fragmentation of C60 in intense laser pulses

The radiation-induced fragmentation of the C60 fullerene was investigated by the tight-binding electron-ion dynamics simulations. In intense laser field, the breathing vibrational mode is much more strongly excited than the pentagonal-pinch mode. The fragmentation effect was found more remarkable at long wavelength lambda>800 nm rather than the resonant wavelengths due to the internal laser-induced dipole force, and the production ratio of C and C2 rapidly grows with increasing wavelength. By such fragmentation law, C atoms, C2 dimers or large Cn fragments could be selectively obtained by changing the laser wavelength. And the fragmentation of C60 by two laser pulses like the multi-step atomic photoionization was investigated.

preprint2013arXiv

A scheme for realizing continuously tunable spectrum in the visible light region based on monatomic carbon chains

We propose a scheme for realizing the continuously tunable spectrum based on monatomic carbon chains. By hybrid density functional calculations, we first show that the direct band gap of monatomic carbon chains change continuously from 1.58 to 3.8 eV as strain is applied from -5 to 10% to the chain, with separated Van Hove singularity peaks enhanced. To realize this tunability, a realistic stretching device is proposed by contacting the chain with graphene sheets, which can apply up to 9% elongation to the chain, yielding tunable light-emitting wavelengths from 345 to 561 nm.

preprint2013arXiv

An age-of-allele test of neutrality for transposable element insertions

How natural selection acts to limit the proliferation of transposable elements (TEs) in genomes has been of interest to evolutionary biologists for many years. To describe TE dynamics in populations, many previous studies have used models of transposition-selection equilibrium that rely on the assumption of a constant rate of transposition. However, since TE invasions are known to happen in bursts through time, this assumption may not be reasonable in natural populations. Here we propose a test of neutrality for TE insertions that does not rely on the assumption of a constant transposition rate. We consider the case of TE insertions that have been ascertained from a single haploid reference genome sequence and have subsequently had their allele frequency estimated in a population sample. By conditioning on the age of an individual TE insertion (using information contained in the number of substitutions that have occurred within the TE sequence since insertion), we determine the probability distribution for the insertion allele frequency in a population sample under neutrality. Taking models of varying population size into account, we then evaluate predictions of our model against allele frequency data from 190 retrotransposon insertions sampled from North American and African populations of Drosophila melanogaster. Using this non-equilibrium model, we are able to explain about 80% of the variance in TE insertion allele frequencies based on age alone. Controlling both for nonequilibrium dynamics of transposition and host demography, we provide evidence for negative selection acting against most TEs as well as for positive selection acting on a small subset of TEs. Our work establishes a new framework for the analysis of the evolutionary forces governing large insertion mutations like TEs, gene duplications or other copy number variants.

preprint2013arXiv

Direct observation of high temperature superconductivity in one-unit-cell FeSe films

Heterostructure based interface engineering has been proved an effective method for finding new superconducting systems and raising superconductivity transition temperature (TC). In previous work on one unit-cell (UC) thick FeSe films on SrTiO3 (STO) substrate, a superconducting-like energy gap as large as 20 meV, was revealed by in situ scanning tunneling microscopy/spectroscopy (STM/STS). Angle resolved photoemission spectroscopy (ARPES) further revealed a nearly isotropic gap of above 15 meV, which closes at a temperature of ~ 65 K. If this transition is indeed the superconducting transition, then the 1-UC FeSe represents the thinnest high TC superconductor discovered so far. However, up to date direct transport measurement of the 1-UC FeSe films has not been reported, mainly because growth of large scale 1-UC FeSe films is challenging and the 1-UC FeSe films are too thin to survive in atmosphere. In this work, we successfully prepared 1-UC FeSe films on insulating STO substrates with non-superconducting FeTe protection layers. By direct transport and magnetic measurements, we provide definitive evidence for high temperature superconductivity in the 1-UC FeSe films with an onset TC above 40 K and a extremely large critical current density JC ~ 1.7*106 A/cm2 at 2 K. Our work may pave the way to enhancing and tailoring superconductivity by interface engineering.

preprint2013arXiv

Electronic analogy of Goos-Hänchen effect: a review

The analogies between optical and electronic Goos-Hänchen effects are established based on electron wave optics in semiconductor or graphene-based nanostructures. In this paper, we give a brief overview of the progress achieved so far in the field of electronic Goos-Hänchen shifts, and show the relevant optical analogies. In particular, we present several theoretical results on the giant positive and negative Goos-Hänchen shifts in various semiconductor or graphen-based nanostructures, their controllability, and potential applications in electronic devices, e.g. spin (or valley) beam splitters.

preprint2013arXiv

Evaluating computational models of explanation using human judgments

We evaluate four computational models of explanation in Bayesian networks by comparing model predictions to human judgments. In two experiments, we present human participants with causal structures for which the models make divergent predictions and either solicit the best explanation for an observed event (Experiment 1) or have participants rate provided explanations for an observed event (Experiment 2). Across two versions of two causal structures and across both experiments, we find that the Causal Explanation Tree and Most Relevant Explanation models provide better fits to human data than either Most Probable Explanation or Explanation Tree models. We identify strengths and shortcomings of these models and what they can reveal about human explanation. We conclude by suggesting the value of pursuing computational and psychological investigations of explanation in parallel.

preprint2013arXiv

Fast and robust population transfer in two-level quantum systems with dephasing noise and/or systematic frequency errors

We design, by invariant-based inverse engineering, driving fields that invert the population of a two-level atom in a given time, robustly with respect to dephasing noise and/or systematic frequency shifts. Without imposing constraints, optimal protocols are insensitive to the perturbations but need an infinite energy. For a constrained value of the Rabi frequency, a flat $π$ pulse is the least sensitive protocol to phase noise but not to systematic frequency shifts, for which we describe and optimize a family of protocols.

preprint2013arXiv

Fast spin control in a two-electron double quantum dot by dynamical invariants

Inverse engineering of electric fields has been recently proposed to achieve fast and robust spin control in a single-electron quantum dot with spin-orbit coupling. In this paper we design, by inverse engineering based on Lewis-Riesenfeld invariants, time-dependent electric fields to realize fast transitions in the selected singlet-triplet subspace of a two-electron double quantum dot. We apply two-mode driving schemes, directly employing the Lewis-Riesenfeld phases, to minimize the electric field necessary to design flexible protocols and perform spin manipulation on the chosen timescale.

preprint2013arXiv

Generalized Fibonacci polynomials and Fibonomial coefficients

The focus of this paper is the study of generalized Fibonacci polynomials and Fibonomial coefficients. The former are polynomials {n} in variables s and t given by {0} = 0, {1} = 1, and {n} = s{n-1}+t{n-2} for n ge 2. The latter are defined by {n choose k} = {n}!/({k}!{n-k}!) where {n}! = {1}{2}...{n}. These quotients are also polynomials in s and t, and specializations give the ordinary binomial coefficients, the Fibonomial coefficients, and the q-binomial coefficients. We present some of their fundamental properties, including a more general recursion for {n}, an analogue of the binomial theorem, a new proof of the Euler-Cassini identity in this setting with applications to estimation of tails of series, and valuations when s and t take on integral values. We also study a corresponding analogue of the Catalan numbers. Conjectures and open problems are scattered throughout the paper.

preprint2013arXiv

Mass acquisition of Dirac fermions in Cr-doped topological insulator Sb2Te3 films

We introduce time-reversal-symmetry breaking by doping Cr atoms into the topmost quintuple layer or into the bulk of Sb2Te3 thin films. We find that even at a high Cr-doping level the Landau level spectrum keeps a good quality, enabling the first demonstration of deviation of the zero-mode Landau level, induced by the acquisition of a mass term in the surface states in the presence of surface or bulk magnetic doping. The magnitude of the mass term in the surface states increases with increasing Cr-doping level. Our observation suggests Cr-doped Sb2Te3 is a promising candidate for the realization of the proposed novel magnetoelectric effects.

preprint2013arXiv

Observations of 6.7 GHz Methanol Masers with EAVN I: VLBI Images of the first Epoch of Observations

Very long baseline interferometry (VLBI) monitoring of the 6.7 GHz methanol maser allows us to measure the internal proper motions of the maser spots and therefore study the gas motion around high-mass young stellar objects. To this end, we have begun monitoring observations with the East-Asian VLBI Network. In this paper we present the results of the first epoch observation for 36 sources, including 35 VLBI images of the methanol maser. Since two independent sources were found in three images, respectively, images of 38 sources were obtained. In 34 sources, more than or equal to 10 spots were detected. The observed spatial scale of the maser distribution was from 9 to 4900 astronomical units, and the following morphological categories were observed: elliptical, arched, linear, paired, and complex. The position of the maser spot was determined to an accuracy of approximately 0.1 mas, sufficiently high to measure the internal proper motion from two years of monitoring observations. The VLBI observation, however, detected only approximately 20% of all maser emission, suggesting that the remaining 80% of the total flux was spread into an undetectable extended distribution. Therefore, in addition to high-resolution observations, it is important to observe the whole structure of the maser emission including extended low-brightness structures, to reveal the associated site of the maser and gas motion.

preprint2013arXiv

Predicting the chemical stability of monatomic chains

A simple model for evaluating the thermal atomic transfer rates in nanosystems [EPL 94, 40002 (2011)] was developed to predict the chemical reaction rates of nanosystems with small gas molecules. The accuracy of the model was verified by MD simulations for molecular adsorption and desorption on a monatomic chain. By the prediction, a monatomic carbon chain should survive for 120 years in the ambient of 1 atm O2 at room temperature, and it is very invulnerable to N2, H2O, NO2, CO and CO2, while a monatomic gold chain quickly ruptures in vacuum. It is worth noting that since the model can be easily applied via common ab initio calculations, it could be widely used in the prediction of chemical stability of nanosystems.

preprint2013arXiv

Predicting the growth rate of helium bubbles in metal tritide

Helium bubbles nucleation and growth in metals or metal tritide is a long-standing problem attracting considerable attention in nuclear industry but the mechanism remains indistinct and predicting the growth rate of helium bubble is inexistence still up to new. Here, the rate of helium bubbles nucleation and growth in metal tritide is developed based on a dynamical model, which describes the diameter of helium bubbles increasing linearly as t**(1/3) in titanium tritide at room temperature, agreeing quite well with the experimental phenomenon. The way of reducing storage temperature from 300 to 225 K or increasing the helium atoms diffusion barrier from 0.81 to 1.1 eV can effectively restrain bubbles growth and prolong lifetime of titanium tritide more than 4 times, which provides a useful reference to relevant experiment exploration and applications. This model also can be used to predict lifetime of new tritium-storage materials and plasma facing materials in nuclear industry.

preprint2013arXiv

Single molecule capturer by doped monatomic carbon chain

The B-doped monatomic carbon chain has fine molecular capture ability to H2O and especially to NO2, better than other doped monatomic carbon chains.At 300 K and 1 atm, the capture probability of the B-doped monatomic carbon chain is appreciable even in a NO2 concentration of 1 p.p.m., and the adsorbates' influence on the quantum transport is notable for the detection. In contrast, the pure monatomic carbon chain shows its invulnerability to N2, O2, H2O, NO2, CO and CO2, and it is incapable for molecule capturing due to too low adsorption ability and weak response on the quantum conductance. In the investigation of these issues, a statistic mechanical model [EPL 94, 40002 (2011); Chin. Phys. Lett. 29, 2012 (080501)] was extended to predict the adsorption and desorption rates of molecules on nanodevices. The theoretical foundation of this model was further discussed and its accuracy was verified by molecular dynamics simulations.

preprint2013arXiv

The fractal nature of the Fibonomial triangle

It is well known that Pascal's triangle exhibits fractal behavior when reduced modulo a prime. We show that the triangle of Fibonomial coefficients has a similar nature modulo two. Specifically, for any $m \ge 0$, the subtriangle consisting of the first $3 \cdot 2^m$ rows is duplicated on the left and right sides of the next $3 \cdot 2^m$ rows, with an inverted triangle of zeros in between. We give three proofs of this fact. The first uses a combinatorial interpretation of the Fibonomials due to Sagan and Savage. The second employs an analogue of Lucas' congruence for the parity of binomial coefficients. The final one is inductive. We also use induction to show that the Fibonomial triangle has a similar structure modulo three. We end with some open questions.

preprint2013arXiv

Tuning the Conductance of Monatomic Carbon Chain

Ab initio calculations show that the conductance of short monatomic carbon chain can be dramatically modified by adhering a single H, N, or O atom to the chain. For example, the conductance of the pristine chain gets about two orders of magnitude smaller if an H atom is adhered to the chain. By a statistical model, the structure of the carbon chain with the single atom adhered is found to be quite stable at room temperature, indicating that the method can be used to tune the conductance of monatomic carbon chain.

preprint2012arXiv

A 95 GHz Class I Methanol Maser Survey Toward A Sample of GLIMPSE Point Sources Associated with BGPS Clumps

We report a survey with the Purple Mountain Observatory (PMO) 13.7-m radio telescope for class I methanol masers from the 95 GHz (8_0 - 7_1 A^+) transition. The 214 target sources were selected by combining information from both the Spitzer GLIMPSE and 1.1 mm BGPS survey catalogs. The observed sources satisfy both the GLIMPSE mid-IR criteria of [3.6]-[4.5]>1.3, [3.6]-[5.8]>2.5, [3.6]-[8.0]>2.5 and 8.0 um magnitude less than 10, and also have an associated 1.1 mm BGPS source. Class I methanol maser emission was detected in 63 sources, corresponding to a detection rate of 29% for this survey. For the majority of detections (43), this is the first identification of a class I methanol maser associated with these sources. We show that the intensity of the class I methanol maser emission is not closely related to mid-IR intensity or the colors of the GLIMPSE point sources, however, it is closely correlated with properties (mass and beam-averaged column density) of the BGPS sources. Comparison of measures of star formation activity for the BGPS sources with and without class I methanol masers indicate that the sources with class I methanol masers usually have higher column density and larger flux density than those without them. Our results predict that the criteria log(S_{int})<-38.0+1.72log(N_{H_{2}}^{beam}) and log(N_{H_{2}}^{beam})>22.1, which utilizes both the integrated flux density (S_{int}) and beam-averaged column density (N_{H_{2}}^{beam}) of the BGPS sources, are very efficient for selecting sources likely to have an associated class I methanol maser. Our expectation is that searches using these criteria will detect 90% of the predicted number of class I methanol masers from the full BGPS catalog (~ 1000), and do so with a high detection efficiency (~75%).

preprint2012arXiv

A multi-layer model for turbulent kinetic energy in pipe flows

A multi-layer model of an energy length function is developed by employing recent results of the authors. The theory predicts the complete, mean streamwise turbulent kinetic-energy profile (MKP), in good agreement with empirical data for a wide range of Reynolds numbers (Re). In particular, a critical $Re_τ$ is predicted, beyond which a scaling anomaly appears and MKP develops a second peak.

preprint2012arXiv

A search for 95 GHz class I methanol masers in molecular outflows

We have observed a sample of 288 molecular outflow sources including 123 high-mass and 165 low-mass sources to search for class I methanol masers at 95 GHz transition and to investigate relationship between outflow characteristics and class I methanol maser emission with the PMO-13.7m radio telescope. Our survey detected 62 sources with 95 GHz methanol masers above 3$σ$ detection limit, which include 47 high-mass sources and 15 low-mass sources. Therefore the detection rate is 38% for high-mass outflow sources and 9% for low-mass outflow sources, suggesting that class I methanol maser is relatively easily excited in high-mass sources. There are 37 newly detected 95 GHz methanol masers (including 27 high-mass and 10 low-mass sources), 19 of which are newly identified (i.e. first identification) class I methanol masers (including 13 high-mass and 6 low-mass sources). Statistical analysis for the distributions of maser detections with the outflow parameters reveals that the maser detection efficiency increases with outflow properties (e.g. mass, momentum, kinetic energy and mechanical luminosity of outflows etc.). Systematic investigations of relationships between the intrinsic luminosity of methanol maser and the outflow properties (including mass, momentum, kinetic energy, bolometric luminosity and mass loss rate of central stellar sources) indicate a positive correlations. This further supports that class I methanol masers are collisionally pumped and associated with shocks, where outflows interact with the surrounding ambient medium.

preprint2012arXiv

Anisotropic vortex lattice structures in the FeSe superconductor

In the recent work by Song et al. [Science 332, 1410 (2011)], the scanning tunneling spectroscopy experiment in the stoichiometric FeSe reveals evidence for nodal superconductivity and strong anisotropy. The nodal structure can be explained with the extended s-wave pairing structure with the mixture of the $s_{x^2+y^2}$ and $s_{x^2y^2}$ pairing symmetries. We calculate the anisotropic vortex structure by using the self-consistent Bogoliubov-de Gennes mean-field theory. In considering the absence of magnetic ordering in the FeSe at the ambient pressure, orbital ordering is introduced, which breaks the $C_4$ lattice symmetry down to $C_2$, to explain the anisotropy in the vortex tunneling spectra.

preprint2012arXiv

Crystal structure and epitaxy of Bi2Te3 films grown on Si

We report comprehensive x-ray diffraction studies of the crystal structure and epitaxy of thin films of the topological insulator Bi2Te3 grown on Si (1 1 1). The films are single crystals of high crystalline quality, which strongly depends on that of their substrates, with in-plane epitaxial relationships of Bi2Te3 [2 1 -3 0] || Si [1 -1 0] and Bi2Te3 [0 1 -1 0] || Si [1 1 -2] along which the lattices of 1x3 Bi2Te3 and 2x2 Si supercells are well matched. As the samples age, we observe loss of crystalline Bi2Te3 film thickness accompanied with roughening of the crystalline interfaces, formation of new crystalline phases as well as compositional and structural modification of the Si substrate, consistent with the diffusion of Te into the Si substrate.

preprint2012arXiv

Double-periodic quasi-periodic graphene superlattice: non-Bragg band gap and electronic transport

Electronic band gap and transport in quasi-periodic graphene superlattice of double-periodic sequence have been investigated. It is found that such quasi-periodic structure can possess a zero-averaged wave number (zero-$\bar{k}$) gap which associated with an unusual Dirac point. Different from Bragg gap, the zero-$\bar{k}$ gap is less sensitive to the incidence angle, and robust against the lattice constants. The locations of Dirac point and multi-Dirac-points in the graphene superlattices of various sequences are also compared. The control of electron transport over the zero-$\bar{k}$ band gap in graphene superlattice may facilitate the development of many graphene-based electronics.

preprint2012arXiv

Electronic Origin of High Temperature Superconductivity in Single-Layer FeSe Superconductor

The latest discovery of high temperature superconductivity signature in single-layer FeSe is significant because it is possible to break the superconducting critical temperature ceiling (maximum Tc~55 K) that has been stagnant since the discovery of Fe-based superconductivity in 2008. It also blows the superconductivity community by surprise because such a high Tc is unexpected in FeSe system with the bulk FeSe exhibiting a Tc at only 8 K at ambient pressure which can be enhanced to 38 K under high pressure. The Tc is still unusually high even considering the newly-discovered intercalated FeSe system A_xFe_{2-y}Se_2 (A=K, Cs, Rb and Tl) with a Tc at 32 K at ambient pressure and possible Tc near 48 K under high pressure. Particularly interesting is that such a high temperature superconductivity occurs in a single-layer FeSe system that is considered as a key building block of the Fe-based superconductors. Understanding the origin of high temperature superconductivity in such a strictly two-dimensional FeSe system is crucial to understanding the superconductivity mechanism in Fe-based superconductors in particular, and providing key insights on how to achieve high temperature superconductivity in general. Here we report distinct electronic structure associated with the single-layer FeSe superconductor. Its Fermi surface topology is different from other Fe-based superconductors; it consists only of electron pockets near the zone corner without indication of any Fermi surface around the zone center. Our observation of large and nearly isotropic superconducting gap in this strictly two-dimensional system rules out existence of node in the superconducting gap. These results have provided an unambiguous case that such a unique electronic structure is favorable for realizing high temperature superconductivity.

preprint2012arXiv

Engineering fast and stable splitting of matter waves

When attempting to split coherent cold atom clouds or a Bose-Einstein condensate (BEC) by bifurcation of the trap into a double well, slow adiabatic following is unstable with respect to any slight asymmetry, and the wave "collapses" to the lower well, whereas a generic fast chopping splits the wave but it also excites it. Shortcuts to adiabaticity engineered to speed up the adiabatic process through non-adiabatic transients, provide instead quiet and robust fast splitting. The non-linearity of the BEC makes the proposed shortcut even more stable.

preprint2012arXiv

Fast and robust spin manipulation in a quantum dot by electric fields

We apply an invariant-based inverse engineering method to control by time-dependent electric fields electron spin dynamics in a quantum dot with spin-orbit coupling in a weak magnetic field. The designed electric fields provide a shortcut to adiabatic processes that flips the spin rapidly, thus avoiding decoherence effects. This approach, being robust with respect to the device-dependent noise, can open new possibilities for the spin-based quantum information processing.

preprint2012arXiv

Group Sparse Additive Models

We consider the problem of sparse variable selection in nonparametric additive models, with the prior knowledge of the structure among the covariates to encourage those variables within a group to be selected jointly. Previous works either study the group sparsity in the parametric setting (e.g., group lasso), or address the problem in the non-parametric setting without exploiting the structural information (e.g., sparse additive models). In this paper, we present a new method, called group sparse additive models (GroupSpAM), which can handle group sparsity in additive models. We generalize the l1/l2 norm to Hilbert spaces as the sparsity-inducing penalty in GroupSpAM. Moreover, we derive a novel thresholding condition for identifying the functional sparsity at the group level, and propose an efficient block coordinate descent algorithm for constructing the estimate. We demonstrate by simulation that GroupSpAM substantially outperforms the competing methods in terms of support recovery and prediction accuracy in additive models, and also conduct a comparative experiment on a real breast cancer dataset.

preprint2012arXiv

Inapproximability After Uniqueness Phase Transition in Two-Spin Systems

A two-state spin system is specified by a 2 x 2 matrix A = {A_{0,0} A_{0,1}, A_{1,0} A_{1,1}} = {β1, 1 γ} where β, γ\ge 0. Given an input graph G=(V,E), the partition function Z_A(G) of a system is defined as Z_A(G) = \sum_{σ: V -> {0,1}} \prod_{(u,v) \in E} A_{σ(u), σ(v)} We prove inapproximability results for the partition function in the region specified by the non-uniqueness condition from phase transition for the Gibbs measure. More specifically, assuming NP \ne RP, for any fixed β, γin the unit square, there is no randomized polynomial-time algorithm that approximates Z_A(G) for d-regular graphs G with relative error ε= 10^{-4}, if d = Ω(Δ(β,γ)), where Δ(β,γ) > 1/(1-βγ) is the uniqueness threshold. Up to a constant factor, this hardness result confirms the conjecture that the uniqueness phase transition coincides with the transition from computational tractability to intractability for Z_A(G). We also show a matching inapproximability result for a region of parameters β, γoutside the unit square, and all our results generalize to partition functions with an external field.

preprint2012arXiv

Interface induced high temperature superconductivity in single unit-cell FeSe films on SrTiO3

Searching for superconducting materials with high transition temperature (TC) is one of the most exciting and challenging fields in physics and materials science. Although superconductivity has been discovered for more than 100 years, the copper oxides are so far the only materials with TC above 77 K, the liquid nitrogen boiling point. Here we report an interface engineering method for dramatically raising the TC of superconducting films. We find that one unit-cell (UC) thick films of FeSe grown on SrTiO3 (STO) substrates by molecular beam epitaxy (MBE) show signatures of superconducting transition above 50 K by transport measurement. A superconducting gap as large as 20 meV of the 1 UC films observed by scanning tunneling microcopy (STM) suggests that the superconductivity could occur above 77 K. The occurrence of superconductivity is further supported by the presence of superconducting vortices under magnetic field. Our work not only demonstrates a powerful way for finding new superconductors and for raising TC, but also provides a well-defined platform for systematic study of the mechanism of unconventional superconductivity by using different superconducting materials and substrates.

preprint2012arXiv

KFe_2Se_2 is the parent compound of K-doped iron selenide superconductors

We elucidate the existing controversies in the newly discovered K-doped iron selenide (KxFe2-ySe2-z) superconductors. The stoichiometric KFe2Se2 with \surd2\times\surd2 charge ordering was identified as the parent compound of KxFe2-ySe2-z superconductor using scanning tunneling microscopy and spectroscopy. The superconductivity is induced in KFe2Se2 by either Se vacancies or interacting with the anti-ferromagnetic K2Fe4Se5 compound. Totally four phases were found to exist in KxFe2-ySe2-z: parent compound KFe2Se2, superconducting KFe2Se2 with \surd2\times\surd5 charge ordering, superconducting KFe2Se2-z with Se vacancies and insulating K2Fe4Se5 with \surd5\times\surd5 Fe vacancy order. The phase separation takes place at the mesoscopic scale under standard molecular beam epitaxy condition.

preprint2012arXiv

Logarithmic distribution of mean velocity and turbulent kinetic energy in a pipe flow

A Lie-group based similarity theory is developed for both momentum and energy distributions in a turbulent pipe flow, leading to asymptotic logarithmic profiles of mean velocity and turbulent kinetic energy. Both channel and pipe data over a wide range of Re yield 0.45 to be the universal Karman constant. A new spatial invariant characterizing outer dynamics is discovered and validated by reliable experimental data. The theory predicts the mean velocity profile (MVP) with 99% accuracy for high Re experimental data (up to 40 millions), and offers a quantitative explanation for recent observation of logarithmic kinetic energy distribution by Hullmak et al. (Phys. Rev. Lett. 108, 094501).

preprint2012arXiv

Non-magnetic and magnetic thiolate-protected Au25 superatoms on Cu(111), Ag(111) and Au(111) surfaces

Geometry, electronic structure, and magnetic properties of methylthiolate-stabilized Au$_{25}$L$_{18}$ and MnAu$_{24}$L$_{18}$ (L = SCH$_3$) clusters adsorbed on noble-metal (111) surfaces have been investigated by using spin-polarized density functional theory computations. The interaction between the cluster and the surface is found to be mediated by charge transfer mainly from or into the ligand monolayer. The electronic properties of the 13-atom metal core remain in all cases rather undisturbed as compared to the isolated clusters in gas phase. The Au$_{25}$L$_{18}$ cluster retains a clear HOMO - LUMO energy gap in the range of 0.7 eV to 1.0 eV depending on the surface. The ligand layer is able to decouple the electronic structure of the magnetic MnAu$_{24}$L$_{18}$ cluster from Au(111) surface, retaning a high local spin moment of close to 5 $μ_{B}$ arising from the spin-polarized Mn(3d) electrons. These computations imply that the thiolate monolayer-protected gold clusters may be used as promising building blocks with tunable energy gaps, tunneling barriers, and magnetic moments for applications in the area of electron and/or spin transport.

preprint2012arXiv

Normal functions, Picard-Fuchs equations, and elliptic fibrations on K3 surfaces

Using Gauss-Manin derivatives of normal functions, we arrive at some remarkable results on the non-triviality of the transcendental regulator for $K_m$ of a very general projective algebraic manifold. Our strongest results are for the transcendental regulator for $K_1$ of a very general $K3$ surface. We also construct an explicit family of $K_1$ cycles on $H \oplus E_8 \oplus E_8$-polarized $K3$ surfaces, and show they are indecomposable by a direct evaluation of the real regulator. Critical use is made of natural elliptic fibrations, hypersurface normal forms, and an explicit parametrization by modular functions.

preprint2012arXiv

Phase Diagram and High Temperature Superconductivity at 65 K in Tuning Carrier Concentration of Single-Layer FeSe Films

Superconductivity in the cuprate superconductors and the Fe-based superconductors is realized by doping the parent compound with charge carriers, or by application of high pressure, to suppress the antiferromagnetic state. Such a rich phase diagram is important in understanding superconductivity mechanism and other physics in the Cu- and Fe-based high temperature superconductors. In this paper, we report a phase diagram in the single-layer FeSe films grown on SrTiO3 substrate by an annealing procedure to tune the charge carrier concentration over a wide range. A dramatic change of the band structure and Fermi surface is observed, with two distinct phases identified that are competing during the annealing process. Superconductivity with a record high transition temperature (Tc) at ~65 K is realized by optimizing the annealing process. The wide tunability of the system across different phases, and its high-Tc, make the single-layer FeSe film ideal not only to investigate the superconductivity physics and mechanism, but also to study novel quantum phenomena and for potential applications.

preprint2012arXiv

Shortcuts to adiabaticity by superadiabatic iterations

Different techniques to speed up quantum adiabatic processes are currently being explored for applications in atomic, molecular and optical physics, such as transport, cooling and expansions, wavepacket splitting, or internal state control. Here we examine the capabilities of superadiabatic iterations to produce a sequence of shortcuts to adiabaticity. The general formalism is worked out as well as examples for population inversion in a two-level system.

preprint2012arXiv

Smoothing Proximal Gradient Method for General Structured Sparse Learning

We study the problem of learning high dimensional regression models regularized by a structured-sparsity-inducing penalty that encodes prior structural information on either input or output sides. We consider two widely adopted types of such penalties as our motivating examples: 1) overlapping group lasso penalty, based on the l1/l2 mixed-norm penalty, and 2) graph-guided fusion penalty. For both types of penalties, due to their non-separability, developing an efficient optimization method has remained a challenging problem. In this paper, we propose a general optimization approach, called smoothing proximal gradient method, which can solve the structured sparse regression problems with a smooth convex loss and a wide spectrum of structured-sparsity-inducing penalties. Our approach is based on a general smoothing technique of Nesterov. It achieves a convergence rate faster than the standard first-order method, subgradient method, and is much more scalable than the most widely used interior-point method. Numerical results are reported to demonstrate the efficiency and scalability of the proposed method.

preprint2012arXiv

Smoothing proximal gradient method for general structured sparse regression

We study the problem of estimating high-dimensional regression models regularized by a structured sparsity-inducing penalty that encodes prior structural information on either the input or output variables. We consider two widely adopted types of penalties of this kind as motivating examples: (1) the general overlapping-group-lasso penalty, generalized from the group-lasso penalty; and (2) the graph-guided-fused-lasso penalty, generalized from the fused-lasso penalty. For both types of penalties, due to their nonseparability and nonsmoothness, developing an efficient optimization method remains a challenging problem. In this paper we propose a general optimization approach, the smoothing proximal gradient (SPG) method, which can solve structured sparse regression problems with any smooth convex loss under a wide spectrum of structured sparsity-inducing penalties. Our approach combines a smoothing technique with an effective proximal gradient method. It achieves a convergence rate significantly faster than the standard first-order methods, subgradient methods, and is much more scalable than the most widely used interior-point methods. The efficiency and scalability of our method are demonstrated on both simulation experiments and real genetic data sets.

preprint2012arXiv

Superconductivity in a single layer alkali-doped FeSe: a weakly coupled two-leg ladder system

We prepare single layer potassium-doped iron selenide (110) film by molecular beam expitaxy. Such a single layer film can be viewed as a two-dimensional system composed of weakly coupled two-leg iron ladders. Scanning tunneling spectroscopy reveals that superconductivity is developed in this two-leg ladder system. The superconducting gap is similar to that of the multi-layer films. However, the Fermi surface topology given by first-principles calculation is remarkably different from that of the bulk materials. Our results suggest that superconducting pairing is very short-ranged or takes place rather locally in iron-chalcogenides. The superconductivity is most likely driven by electron-electron correlation effect and is insensitive to the change of Fermi surfaces.

preprint2012arXiv

Suppression of Superconductivity by Twin Boundaries in FeSe

Low-temperature scanning tunneling microscopy and spectroscopy are employed to investigate twin boundaries in stoichiometric FeSe films grown by molecular beam epitaxy. Twin boundaries can be unambiguously identified by imaging the 90° change in the orientation of local electronic dimers from Fe site impurities on either side. Twin boundaries run at approximately 45° to the Fe-Fe bond directions, and noticeably suppress the superconducting gap, in contrast with the recent experimental and theoretical findings in other iron pnictides. Furthermore, vortices appear to accumulate on twin boundaries, consistent with the degraded superconductivity there. The variation in superconductivity is likely caused by the increased Se height in the vicinity of twin boundaries, providing the first local evidence for the importance of this height to the mechanism of superconductivity.

preprint2012arXiv

The Complexity of Non-Monotone Markets

We introduce the notion of non-monotone utilities, which covers a wide variety of utility functions in economic theory. We then prove that it is PPAD-hard to compute an approximate Arrow-Debreu market equilibrium in markets with linear and non-monotone utilities. Building on this result, we settle the long-standing open problem regarding the computation of an approximate Arrow-Debreu market equilibrium in markets with CES utility functions, by proving that it is PPAD-complete when the Constant Elasticity of Substitution parameter ρis any constant less than -1.

preprint2012arXiv

Topological p-n Junction

We consider a junction between surface $p$-type and surface $n$-type on an ideal topological insulator in which carrier type and density in two adjacent regions are locally controlled by composition graded doping or electrical gating. Such junction setting on topological insulators are fundamental for possible device application. A single gapless chiral edge state localized along the junction interface appears in the presence of an external magnetic field, and it can be probed by scanning tunneling microscopy and transport measurements. We propose to realize this topological \emph{p-n} junction in (Bi$_{1-x}$Sb$_x$)$_2$Te$_3$, which has insulating bulk properties and a tunable surface state across the Dirac cone.

preprint2012arXiv

VLA/JVLA Monitoring of Bright Northern Radio Sources

We report multiple epoch VLA/JVLA observations of 89 northern hemisphere sources, most with 37\,GHz flux density > 1 Jy, observed at 4.8, 8.5, 33.5, and 43.3 GHz. The high frequency selection leads to a predominantly flat spectrum sample, with 85% of our sources being in the Planck Early Release Compact Source Catalog (ERCSC). These observations allow us to: 1) validate Planck's 30 and 44 GHz flux density scale, 2) extend the radio SEDs of Planck sources to lower frequencies allowing for the full 5-857GHz regime to be studied, and 3) characterize the variability of these sources. At 30 GHz and 44 GHz, the JVLA and Planck flux densities agree to within 3%. On timescales of less than two months the median variability of our sources is 2%. On timescales of about a year the median variability increases to 14%. Using the WMAP 7-year data, the 30 GHz median variability on a 1-6 years timescale is 16%.

preprint2011arXiv

A 95 GHz Class I Methanol Maser Survey Toward GLIMPSE Extended Green Objects (EGOs)

We report the results of a systematic survey for 95 GHz class I methanol masers towards a new sample of 192 massive young stellar object (MYSO) candidates associated with ongoing outflows (known as extended green objects or EGOs) identified from the Spitzer GLIMPSE survey. The observations were made with the Australia Telescope National Facility (ATNF) Mopra 22-m radio telescope and resulted in the detection of 105 new 95 GHz class I methanol masers. For 92 of the sources our observations provide the first identification of a class I maser transition associated with these objects (i.e. they are new class I methanol maser sources). Our survey proves that there is indeed a high detection rate (55%) of class I methanol masers towards EGOs. Comparison of the GLIMPSE point sources associated with EGOs with and without class I methanol maser detections shows they have similar mid-IR colors, with the majority meeting the color selection criteria -0.6<[5.8]-[8.0]<1.4 and 0.5<[3.6]-[4.5]<4.0. Investigations of the IRAC and MIPS 24 um colors and the associated millimeter dust clump properties (mass and density) of the EGOs for the sub-samples based on which class of methanol masers they are associated with suggests that the stellar mass range associated with class I methanol masers extends to lower masses than for class II methanol masers, or alternatively class I methanol masers may be associated with more than one evolutionary phase during the formation of a high-mass star.

preprint2011arXiv

A Smoothing Stochastic Gradient Method for Composite Optimization

We consider the unconstrained optimization problem whose objective function is composed of a smooth and a non-smooth conponents where the smooth component is the expectation a random function. This type of problem arises in some interesting applications in machine learning. We propose a stochastic gradient descent algorithm for this class of optimization problem. When the non-smooth component has a particular structure, we propose another stochastic gradient descent algorithm by incorporating a smoothing method into our first algorithm. The proofs of the convergence rates of these two algorithms are given and we show the numerical performance of our algorithm by applying them to regularized linear regression problems with different sets of synthetic data.

preprint2011arXiv

A Turing model of pattern formation in atmospheric pressure gas discharges

In this letter we propose a Turing model of the formation of patterns of visible light emission intensity in atmospheric pressure gas discharges. The electron density and the electron temperature take the roles of activator and inhibitor respectively in a two-reactant Turing model, with the activator diffusion coeffcient being much smaller than that of the inhibitor, and ionization and excitation from excited state atoms considered as the dominant reaction processes. The model predicts striations in a 1D system, which quantitatively agree with experimental results in terms of the spatial variation scale lengths. Additionally, the model also predicts changes in the discharge structure observed experimentally when input power and gas pressure are varied.

preprint2011arXiv

Band structure engineering in (Bi1-xSbx)2Te3 ternary topological insulators

Three-dimensional (3D) topological insulators (TI) are novel quantum materials with insulating bulk and topologically protected metallic surfaces with Dirac-like band structure. The spin-helical Dirac surface states are expected to host exotic topological quantum effects and find applications in spintronics and quantum computation. The experimental realization of these ideas requires fabrication of versatile devices based on bulk-insulating TIs with tunable surface states. The main challenge facing the current TI materials exemplified by Bi2Se3 and Bi2Te3 is the significant bulk conduction, which remains unsolved despite extensive efforts involving nanostructuring, chemical doping and electrical gating. Here we report a novel approach for engineering the band structure of TIs by molecular beam epitaxy (MBE) growth of (Bi1-xSbx)2Te3 ternary compounds. Angle-resolved photoemission spectroscopy (ARPES) and transport measurements show that the topological surface states exist over the entire composition range of (Bi1-xSbx)2Te3 (x = 0 to 1), indicating the robustness of bulk Z2 topology. Most remarkably, the systematic band engineering leads to ideal TIs with truly insulating bulk and tunable surface state across the Dirac point that behave like one quarter of graphene. This work demonstrates a new route to achieving intrinsic quantum transport of the topological surface states and designing conceptually new TI devices with well-established semiconductor technology.

preprint2011arXiv

Complexity of Counting CSP with Complex Weights

We give a complexity dichotomy theorem for the counting Constraint Satisfaction Problem (#CSP in short) with complex weights. To this end, we give three conditions for its tractability. Let F be any finite set of complex-valued functions, then we prove that #CSP(F) is solvable in polynomial time if all three conditions are satisfied; and is #P-hard otherwise. Our complexity dichotomy generalizes a long series of important results on counting problems: (a) the problem of counting graph homomorphisms is the special case when there is a single symmetric binary function in F; (b) the problem of counting directed graph homomorphisms is the special case when there is a single not-necessarily-symmetric binary function in F; and (c) the standard form of #CSP is when all functions in F take values in {0,1}.

preprint2011arXiv

Controllable Goos-Hänchen shifts and spin beam splitter for ballistic electrons in a parabolic quantum well under a uniform magnetic field

The quantum Goos-Hänchen shift for ballistic electrons is investigated in a parabolic potential well under a uniform vertical magnetic field. It is found that the Goos-Hänchen shift can be negative as well as positive, and becomes zero at transmission resonances. The beam shift depends not only on the incident energy and incidence angle, but also on the magnetic field and Landau quantum number. Based on these phenomena, we propose an alternative way to realize the spin beam splitter in the proposed spintronic device, which can completely separate spin-up and spin-down electron beams by negative and positive Goos-Hänchen shifts.

preprint2011arXiv

Crossover between weak localization and weak antilocalization in magnetically doped topological insulator

Topological insulators (TI) are a new class of quantum materials with insulating bulk enclosed by topologically protected metallic boundaries. The surface states of three-dimensional TIs have spin helical Dirac structure, and are robust against time reversal invariant perturbations. This extraordinary property is notably exemplified by the absence of backscattering by nonmagnetic impurities and the weak antilocalization (WAL) of Dirac fermions. Breaking the time reversal symmetry (TRS) by magnetic element doping is predicted to create a variety of exotic topological magnetoelectric effects. Here we report transport studies on magnetically doped TI Cr-Bi2Se3. With increasing Cr concentration, the low temperature electrical conduction exhibits a characteristic crossover from WAL to weak localization (WL). In the heavily doped regime where WL dominates at the ground state, WAL reenters as temperature rises, but can be driven back to WL by strong magnetic field. These complex phenomena can be explained by a unified picture involving the evolution of Berry phase with the energy gap opened by magnetic impurities. This work demonstrates an effective way to manipulate the topological transport properties of the TI surface states by TRS-breaking perturbations.

preprint2011arXiv

Domain wall partition function of the eight-vertex model with a non-diagonal reflecting end

With the help of the Drinfeld twist or factorizing F-matrix for the eight-vertex SOS model, we obtain the explicit determinant expression of the partition function of the eight-vertex model with a generic non-diagonal reflecting end and domain wall boundary condition. Our result shows that, contrary to the eight-vertex model without a reflection end, the partition function can be expressed as a single determinant.

preprint2011arXiv

Drinfeld twist and symmetric Bethe vectors of the open XYZ chain with non-diagonal boundary terms

With the help of the Drinfeld twist or factorizing F-matrix for the eight-vertex solid-on-solid (SOS) model, we find that in the F-basis provided by the twist the two sets of pseudo-particle creation operators simultaneously take completely symmetric and polarization free form. This allows us to obtain the explicit and completely symmetric expressions of the two sets of Bethe states of the model.

preprint2011arXiv

Electronic band gap and transport in Fibonacci quasi-periodic graphene superlattice

We investigate electronic band gap and transport in Fibonacci quasi-periodic graphene superlattice. It is found that such structure can possess a zero-$\bar{k}$ gap which exists in all Fibonacci sequences. Different from Bragg gap, zero-$\bar{k}$ gap associated with Dirac point is less sensitive to the incidence angle and lattice constants. The defect mode appeared inside the zero-$\bar{k}$ gap has a great effect on transmission, conductance and shot noise, which can be applicable to control the electron transport.

preprint2011arXiv

Electronic beam shifts in monolayer graphene superlattice

Electronic analogue of generalized Goos-Hänchen shifts is investigated in the monolayer graphene superlattice with one-dimensional periodic potentials of square barriers. It is found that the lateral shifts for the electron beam transmitted through the monolayer graphene superlattice can be negative as well as positive near the band edges of zero-$\bar{k}$ gap, which are different from those near the band edges of Bragg gap. These negative and positive beam shifts have close relation to the Dirac point. When the condition $q_A d_A= -q_B d_B= m π$ ($m=1,2,3...$) is satisfied, the beam shifts can be controlled from negative to positive when the incident energy is above the Dirac point, and vice versa. In addition, the beam shifts can be greatly enhanced by the defect mode inside the zero-$\bar{k}$ gap. These intriguing phenomena can be verified in a relatively simple optical setup, and have potential applications in the graphene-based electron wave devices.

preprint2011arXiv

Fast transitionless expansions of cold atoms in optical Gaussian beam traps

We study fast expansions of cold atoms in a three-dimensional Gaussian-beam optical trap. Three different methods to avoid final motional excitation are compared: inverse engineering using Lewis-Riesenfeld invariants, which provides the best overall performance, a bang-bang approach, and a fast adiabatic approach. We analyze the excitation effect of anharmonic terms, radial-longitudinal coupling, and radial-frequency mismatch. In the inverse engineering approach these perturbations can be suppressed or mitigated by increasing the laser beam waist.

preprint2011arXiv

Fast transport of Bose-Einstein condensates

We propose an inverse method to accelerate without final excitation the adiabatic transport of a Bose Einstein condensate. The method, applicable to arbitrary potential traps, is based on a partial extension of the Lewis-Riesenfeld invariants, and provides transport protocols that satisfy exactly the no-excitation conditions without constraints or approximations. This inverse method is complemented by optimizing the trap trajectory with respect to different physical criteria and by studying the effect of noise.

preprint2011arXiv

Fermi Level Tuning of Epitaxial Sb2Te3 Thin Films on Graphene by Regulating Intrinsic Defects and Substrate Transfer Doping

High-quality Sb2Te3 films are obtained by molecular beam epitaxy on graphene substrate and investigated by in situ scanning tunneling microscopy/spectroscopy. Intrinsic defects responsible for the natural p-type conductivity of Sb2Te3 are identified to be the Sb vacancies and SbTe antisites in agreement with first-principles calculations. By minimizing defect densities, coupled with a transfer doping by the graphene substrate, the Fermi level of Sb2Te3 thin films can be tuned over the entire range of the bulk band gap. This establishes the necessary condition to explore topological insulator behaviors near the Dirac point.

preprint2011arXiv

Graph Homomorphisms with Complex Values: A Dichotomy Theorem

Graph homomorphism has been studied intensively. Given an m x m symmetric matrix A, the graph homomorphism function is defined as \[Z_A (G) = \sum_{f:V->[m]} \prod_{(u,v)\in E} A_{f(u),f(v)}, \] where G = (V,E) is any undirected graph. The function Z_A can encode many interesting graph properties, including counting vertex covers and k-colorings. We study the computational complexity of Z_A for arbitrary symmetric matrices A with algebraic complex values. Building on work by Dyer and Greenhill, Bulatov and Grohe, and especially the recent beautiful work by Goldberg, Grohe, Jerrum and Thurley, we prove a complete dichotomy theorem for this problem. We show that Z_A is either computable in polynomial-time or #P-hard, depending explicitly on the matrix A. We further prove that the tractability criterion on A is polynomial-time decidable.

preprint2011arXiv

Karman constant and accurate mean flow prediction in a turbulent pipe

The Karman constant κ- widely used in atmospheric science and engineering turbulence modelling, and proposed by Prandtl in 1925 and von Karman in 1930 to describe the mean velocity of a turbulent wall-bounded flow - leads to a logarithmic profile in an overlap region near the wall. For over eighty years, its value was believed to be ~0.41. But more recently, many argue that it is not a constant, because of measured variations in different flows and at different Reynolds numbers (Re). Here, a multi-layer analytic theory is shown to lead to a re-interpretation of κas a global constant for both the overlap region and outer flow, and to yield a new method for its measurement. The newly determined value is 0.45 for both channel and pipe. It is shown that this new κ, together with other wall constants, yields a 99% accuracy in the prediction of mean velocity data at all points in high Re (up to 40 million) pipe flow. The theory also describes finite Re effect, and discovers a transition at the friction Re (i.e. Re_τ)=5000. An accurate model for the prediction of turbulent transport in canonical pipe and channel flows is achieved here, and we propose the model to be valid for a wide class of turbulent flows.

preprint2011arXiv

Landau quantization and the thickness limit of topological insulator thin films of Sb2Te3

We report the experimental observation of Landau quantization of molecular beam epitaxy grown Sb2Te3 thin films by a low-temperature scanning tunneling microscope. Different from all the reported systems, the Landau quantization in Sb2Te3 topological insulator is not sensitive to the intrinsic substitutional defects in the films. As a result, a nearly perfect linear energy dispersion of surface states as 2D massless Dirac fermion system is achieved. We demonstrate that 4 quintuple layers are the thickness limit for Sb2Te3 thin film being a 3D topological insulator. The mechanism of the Landau level broadening is discussed in terms of enhanced quasiparticle lifetime.

preprint2011arXiv

Lewis-Riesenfeld invariants and transitionless tracking algorithm

Different methods have been recently put forward and implemented experimentally to inverse engineer the time dependent Hamiltonian of a quantum system and accelerate slow adiabatic processes via non-adiabatic shortcuts. In the "transitionless tracking algorithm" proposed by Berry, shortcut Hamiltonians are designed so that the system follows exactly, in an arbitrarily short time, the approximate adiabatic path defined by a reference Hamiltonian. A different approach is based on designing first a Lewis-Riesenfeld invariant to carry the eigenstates of a Hamiltonian from specified initial to final configurations, again in an arbitrary time, and then constructing from the invariant the transient Hamiltonian connecting these boundary configurations. We show that the two approaches, apparently quite different in form and so far in results, are in fact strongly related and potentially equivalent, so that the inverse-engineering operations in one of them can be reinterpreted and understood in terms of the concepts and operations of the other one. We study as explicit examples the expansions of time-dependent harmonic traps and state preparation of two level systems.

preprint2011arXiv

Magnetic phase diagram of interacting nanoparticle systems under the mean field model

The disordered random-anisotropy magnetic nanoparticle systems with competing dipolar interactions and ferromagnetic exchange couplings are investigated by Monte Carlo simulations. Superspin glass (SSG) and superferromagnetic (SFM) behaviors are found at low temperatures depending on the interactions. Based on the mean field approximation, the Curie-Weiss temperature TCW = 0 is suggested as the phase boundary between the SSG systems and the SFM systems, which is convinced by the spontaneous magnetizations and relaxations. The magnetic phase diagram is plotted.

preprint2011arXiv

Multiple Schrödinger pictures and dynamics in shortcuts to adiabaticity

A Schrödinger equation may be transformed by unitary operators into dynamical equations in different interaction pictures which share with it a common physical frame, i.e., the same underlying interactions, processes and dynamics. In contrast to this standard scenario, other relations are also possible, such as a common interaction-picture dynamical equation corresponding to several Schrödinger equations that represent different physics. This may enable us to design alternative and feasible experimental routes for operations that are a priori difficult or impossible to perform. The power of this concept is exemplified by engineering Hamiltonians that improve the performance or make realizable several shortcuts to adiabaticity.

preprint2011arXiv

Optimal trajectories for efficient atomic transport without final excitation

We design optimal harmonic-trap trajectories to transport cold atoms without final excitation, combining an inverse engineering techniqe based on Lewis-Riesenfeld invariants with optimal control theory. Since actual traps are not really harmonic, we keep the relative displacement between the center of mass and the trap center bounded. Under this constraint, optimal protocols are found according to different physical criteria. The minimum time solution has a "bang-bang" form, and the minimum displacement solution is of "bang-off-bang" form. The optimal trajectories for minimizing the transient energy are also discussed.

preprint2011arXiv

Phase Separation and Magnetic Order in K-doped Iron Selenide Superconductor

Alkali-doped iron selenide is the latest member of high Tc superconductor family, and its peculiar characters have immediately attracted extensive attention. We prepared high-quality potassium-doped iron selenide (KxFe2-ySe2) thin films by molecular beam epitaxy and unambiguously demonstrated the existence of phase separation, which is currently under debate, in this material using scanning tunneling microscopy and spectroscopy. The stoichiometric superconducting phase KFe2Se2 contains no iron vacancies, while the insulating phase has a \surd5\times\surd5 vacancy order. The iron vacancies are shown always destructive to superconductivity in KFe2Se2. Our study on the subgap bound states induced by the iron vacancies further reveals a magnetically-related bipartite order in the superconducting phase. These findings not only solve the existing controversies in the atomic and electronic structures in KxFe2-ySe2, but also provide valuable information on understanding the superconductivity and its interplay with magnetism in iron-based superconductors.

preprint2011arXiv

Power-Law Decay of Standing Waves on the Surface of Topological Insulators

We propose a general theory on the standing waves (quasiparticle interference pattern) caused by the scattering of surface states off step edges in topological insulators, in which the extremal points on the constant energy contour of surface band play the dominant role. Experimentally we image the interference patterns on both Bi$_2$Te$_3$ and Bi$_2$Se$_3$ films by measuring the local density of states using a scanning tunneling microscope. The observed decay indices of the standing waves agree excellently with the theoretical prediction: In Bi$_2$Se$_3$, only a single decay index of -3/2 exists; while in Bi$_2$Te$_3$ with strongly warped surface band, it varies from -3/2 to -1/2 and finally to -1 as the energy increases. The -1/2 decay indicates that the suppression of backscattering due to time-reversal symmetry does not necessarily lead to a spatial decay rate faster than that in the conventional two-dimensional electron system. Our formalism can also explain the characteristic scattering wave vectors of the standing wave caused by non-magnetic impurities on Bi$_2$Te$_3$.

preprint2011arXiv

Scalar products of the open XYZ chain with non-diagonal boundary terms

With the help of the F-basis provided by the Drinfeld twist or factorizing F-matrix of the eight-vertex solid-on-solid (SOS) model, we obtain the determinant representations of the scalar products of Bethe states for the open XYZ chain with non-diagonal boundary terms. By taking the on shell limit, we obtain the determinant representations (or Gaudin formula) of the norms of the Bethe states.

preprint2011arXiv

VLBI observations of 10 CSO candidates: expansion velocities of hot spots

Observations of ten Compact Symmetric Objects ({\rm CSO}) candidates have been made with the Very Long Baseline Array at 8.4 GHz in 2005 and with a combined Chinese and European VLBI array at 8.4 GHz in 2009. The 2009 observations incorporate for the first time the two new Chinese telescopes at Miyun and Kunming for international astrophysical observations. The observational data, in combination with archival VLBA data from previous epochs, have been used to derive the proper motions of the VLBI components. Because of the long time baseline of $\sim$16 years of the VLBI data sets, the expansion velocities of the hot spots can be measured at an accuracy as high as $\sim$1.3 $μ$as yr$^{-1}$. Six of the ten sources are identified as CSOs with a typical double or triple morphology on the basis of both spectral index maps and their mirror-symmetry of proper motions of the terminal hot spots. The compact double source J1324+4048 is also identified as a CSO candidate. Among the three remaining sources, J1756+5748 and J2312+3847 are identified as core-jet sources with proper motions of their jet components relating to systemic source expansion. The third source J0017+5312 is likely also a core-jet source, but a robust detection of a core is needed for an unambiguous identification. The kinematic ages of the CSOs derived from proper motions range from 300 to 2500 years. The kinematic age distribution of the CSOs confirm an overabundance of compact young CSOs with ages less than 500 years. CSOs with known kinematic ages may be used to study the dynamical evolution of extragalactic radio sources at early stages.

preprint2010arXiv

A Complexity View of Markets with Social Influence

In this paper, inspired by the work of Megiddo on the formation of preferences and strategic analysis, we consider an early market model studied in the field of economic theory, in which each trader's utility may be influenced by the bundles of goods obtained by her social neighbors. The goal of this paper is to understand and characterize the impact of social influence on the complexity of computing and approximating market equilibria. We present complexity-theoretic and algorithmic results for approximating market equilibria in this model with focus on two concrete influence models based on the traditional linear utility functions. Recall that an Arrow-Debreu market equilibrium in a conventional exchange market with linear utility functions can be computed in polynomial time by convex programming. Our complexity results show that even a bounded-degree, planar influence network can significantly increase the difficulty of equilibrium computation even in markets with only a constant number of goods. Our algorithmic results suggest that finding an approximate equilibrium in markets with hierarchical influence networks might be easier than that in markets with arbitrary neighborhood structures. By demonstrating a simple market with a constant number of goods and a bounded-degree, planar influence graph whose equilibrium is PPAD-hard to approximate, we also provide a counterexample to a common belief, which we refer to as the myth of a constant number of goods, that equilibria in markets with a constant number of goods are easy to compute or easy to approximate.

preprint2010arXiv

A Decidable Dichotomy Theorem on Directed Graph Homomorphisms with Non-negative Weights

The complexity of graph homomorphism problems has been the subject of intense study. It is a long standing open problem to give a (decidable) complexity dichotomy theorem for the partition function of directed graph homomorphisms. In this paper, we prove a decidable complexity dichotomy theorem for this problem and our theorem applies to all non-negative weighted form of the problem: given any fixed matrix A with non-negative algebraic entries, the partition function Z_A(G) of directed graph homomorphisms from any directed graph G is either tractable in polynomial time or #P-hard, depending on the matrix A. The proof of the dichotomy theorem is combinatorial, but involves the definition of an infinite family of graph homomorphism problems. The proof of its decidability is algebraic using properties of polynomials.

preprint2010arXiv

Electron interaction-driven insulating ground state in Bi2Se3 topological insulators in the two dimensional limit

We report a transport study of ultrathin Bi2Se3 topological insulators with thickness from one quintuple layer to six quintuple layers grown by molecular beam epitaxy. At low temperatures, the film resistance increases logarithmically with decreasing temperature, revealing an insulating ground state. The sharp increase of resistance with magnetic field, however, indicates the existence of weak antilocalization, which should reduce the resistance as temperature decreases. We show that these apparently contradictory behaviors can be understood by considering the electron interaction effect, which plays a crucial role in determining the electronic ground state of topological insulators in the two dimensional limit.

preprint2010arXiv

Fast atomic transport without vibrational heating

We use the dynamical invariants associated with the Hamiltonian of an atom in a one dimensional moving trap to inverse engineer the trap motion and perform fast atomic transport without final vibrational heating. The atom is driven non-adiabatically through a shortcut to the result of adiabatic, slow trap motion. For harmonic potentials this only requires designing appropriate trap trajectories, whereas perfect transport in anharmonic traps may be achieved by applying an extra field to compensate the forces in the rest frame of the trap. The results can be extended to atom stopping or launching. The limitations due to geometrical constraints, energies and accelerations involved are analyzed, as well as the relation to previous approaches (based on classical trajectories or "fast-forward" and "bang-bang" methods) which can be integrated in the invariant-based framework.

preprint2010arXiv

Goos-Hänchen-like shifts for Dirac fermions in monolayer graphene barrier

We investigate the Goos-Hänchen-like shifts for Dirac fermions in transmission through a monolayer graphene barrier. The lateral shifts, as the functions of the barrier's width and the incidence angle, can be negative and positive in Klein tunneling and classical motion, respectively. Due to their relations to the transmission gap, the lateral shifts can be enhanced by the transmission resonances when the incidence angle is less than the critical angle for total reflection, while their magnitudes become only the order of Fermi wavelength when the incidence angle is larger than the critical angle. These tunable beam shifts can also be modulated by the height of potential barrier and the induced gap, which gives rise to the applications in graphene-based devices.

preprint2010arXiv

Graph-Structured Multi-task Regression and an Efficient Optimization Method for General Fused Lasso

We consider the problem of learning a structured multi-task regression, where the output consists of multiple responses that are related by a graph and the correlated response variables are dependent on the common inputs in a sparse but synergistic manner. Previous methods such as l1/l2-regularized multi-task regression assume that all of the output variables are equally related to the inputs, although in many real-world problems, outputs are related in a complex manner. In this paper, we propose graph-guided fused lasso (GFlasso) for structured multi-task regression that exploits the graph structure over the output variables. We introduce a novel penalty function based on fusion penalty to encourage highly correlated outputs to share a common set of relevant inputs. In addition, we propose a simple yet efficient proximal-gradient method for optimizing GFlasso that can also be applied to any optimization problems with a convex smooth loss and the general class of fusion penalty defined on arbitrary graph structures. By exploiting the structure of the non-smooth ''fusion penalty'', our method achieves a faster convergence rate than the standard first-order method, sub-gradient method, and is significantly more scalable than the widely adopted second-order cone-programming and quadratic-programming formulations. In addition, we provide an analysis of the consistency property of the GFlasso model. Experimental results not only demonstrate the superiority of GFlasso over the standard lasso but also show the efficiency and scalability of our proximal-gradient method.

preprint2010arXiv

Graph-Valued Regression

Undirected graphical models encode in a graph $G$ the dependency structure of a random vector $Y$. In many applications, it is of interest to model $Y$ given another random vector $X$ as input. We refer to the problem of estimating the graph $G(x)$ of $Y$ conditioned on $X=x$ as ``graph-valued regression.'' In this paper, we propose a semiparametric method for estimating $G(x)$ that builds a tree on the $X$ space just as in CART (classification and regression trees), but at each leaf of the tree estimates a graph. We call the method ``Graph-optimized CART,'' or Go-CART. We study the theoretical properties of Go-CART using dyadic partitioning trees, establishing oracle inequalities on risk minimization and tree partition consistency. We also demonstrate the application of Go-CART to a meteorological dataset, showing how graph-valued regression can provide a useful tool for analyzing complex data.

preprint2010arXiv

Growth of quantum well films of topological insulator Bi2Se3 on insulating substrate

Insulating substrates are crucial for electrical transport study and room temperature application of topological insulator films at thickness of only several nanometers. High quality quantum well films of Bi2Se3, a typical three-dimensional topological insulator, have been grown on α-Al2O3 (sapphire) (0001) by molecular beam epitaxy. The films exhibit well-defined quantum well states and surface states, suggesting the uniform thickness over macroscopic area. The Bi2Se3 thin films on sapphire (0001) provide a good system to study low-dimensional physics of topological insulators since conduction contribution from the substrate is negligibly small.

preprint2010arXiv

Landau Quantization of Massless Dirac Fermions in Topological Insulator

The recent theoretical prediction and experimental realization of topological insulators (TI) has generated intense interest in this new state of quantum matter. The surface states of a three-dimensional (3D) TI such as Bi_2Te_3, Bi_2Se_3 and Sb_2Te_3 consist of a single massless Dirac cones. Crossing of the two surface state branches with opposite spins in the materials is fully protected by the time reversal (TR) symmetry at the Dirac points, which cannot be destroyed by any TR invariant perturbation. Recent advances in thin-film growth have permitted this unique two-dimensional electron system (2DES) to be probed by scanning tunneling microscopy (STM) and spectroscopy (STS). The intriguing TR symmetry protected topological states were revealed in STM experiments where the backscattering induced by non-magnetic impurities was forbidden. Here we report the Landau quantization of the topological surface states in Bi_2Se_3 in magnetic field by using STM/STS. The direct observation of the discrete Landau levels (LLs) strongly supports the 2D nature of the topological states and gives direct proof of the nondegenerate structure of LLs in TI. We demonstrate the linear dispersion of the massless Dirac fermions by the square-root dependence of LLs on magnetic field. The formation of LLs implies the high mobility of the 2DES, which has been predicted to lead to topological magneto-electric effect of the TI.

preprint2010arXiv

Non-negative Weighted #CSPs: An Effective Complexity Dichotomy

We prove a complexity dichotomy theorem for all non-negative weighted counting Constraint Satisfaction Problems (CSP). This caps a long series of important results on counting problems including unweighted and weighted graph homomorphisms and the celebrated dichotomy theorem for unweighted #CSP. Our dichotomy theorem gives a succinct criterion for tractability. If a set F of constraint functions satisfies the criterion, then the counting CSP problem defined by F is solvable in polynomial time; if it does not satisfy the criterion, then the problem is #P-hard. We furthermore show that the question of whether F satisfies the criterion is decidable in NP. Surprisingly, our tractability criterion is simpler than the previous criteria for the more restricted classes of problems, although when specialized to those cases, they are logically equivalent. Our proof mainly uses Linear Algebra, and represents a departure from Universal Algebra, the dominant methodology in recent years.

preprint2010arXiv

On Tractable Exponential Sums

We consider the problem of evaluating certain exponential sums. These sums take the form $\sum_{x_1,...,x_n \in Z_N} e^{f(x_1,...,x_n) {2 πi / N}} $, where each x_i is summed over a ring Z_N, and f(x_1,...,x_n) is a multivariate polynomial with integer coefficients. We show that the sum can be evaluated in polynomial time in n and log N when f is a quadratic polynomial. This is true even when the factorization of N is unknown. Previously, this was known for a prime modulus N. On the other hand, for very specific families of polynomials of degree \ge 3, we show the problem is #P-hard, even for any fixed prime or prime power modulus. This leads to a complexity dichotomy theorem - a complete classification of each problem to be either computable in polynomial time or #P-hard - for a class of exponential sums. These sums arise in the classifications of graph homomorphisms and some other counting CSP type problems, and these results lead to complexity dichotomy theorems. For the polynomial-time algorithm, Gauss sums form the basic building blocks. For the hardness results, we prove group-theoretic necessary conditions for tractability. These tests imply that the problem is #P-hard for even very restricted families of simple cubic polynomials over fixed modulus N.

preprint2010arXiv

Photon spin operator and Pauli matrix

Any polarization vector of a plane wave can be decomposed into a pair of mutually orthogonal base vectors, known as a polarization basis. Regarding this decomposition as a quasi-unitary transformation from a three-component vector to a corresponding two-component spinor, one is led to a representation formalism for the photon spin. The spin operator $\hat{\boldsymbol γ}$ defined on the space of unit spinors, referred to as the Jones space, has only component along the wave vector and is represented by one of the Pauli matrices in the commonly used polarization basis. It is deformed by the quasi-unitary transformation from the spin operator that is defined on the space of unit polarization vectors, referred to as the Pancharatnam space. On the basis of this theory, it is shown that the Cartesian components of spin operator $\hat{\boldsymbol γ}$ are mutually commutative and the spin angular momentum in units of $\hbar$ is exactly the component of the Stokes vector along the wave vector.

preprint2010arXiv

Revising the kinematics of 12GHz CH3OH masers in W3(OH)

We derive accurate proper motions of the CH3OH 12 GHz masers towards the W3(OH) UCHII region, employing seven epochs of VLBA observations spanning a time interval of about 10 yr. The achieved velocity accuracy is of the order of 0.1 km/s, adequate to precisely measure the relative velocities of most of the 12 GHz masers in W3(OH), with amplitude varying in the range 0.3 - 3 km/s. Towards W3(OH), the most intense 12 GHz masers concentrate in a small area towards the north (the northern clump) of the UCHII region. We have compared the proper motions of the CH3OH 12 GHz masers with those (derived from literature data) of the OH 6035 MHz masers, emitting from the same region of the methanol masers. In the northern clump, the two maser emissions emerge from nearby (but likely distinct) cloudlets of masing gas with, in general, a rather smooth variation of line-of-sight and sky-projected velocities, which suggests some connection of the environments and kinematics traced by both maser types. The conical outflow model, previously proposed to account for the 12 GHz maser kinematics in the northern clump, does not reproduce the new, accurate measurements of 12 GHz maser proper motions and has to be rejected. We focus on the subset of 12 GHz masers of the northern clump belonging to the "linear structure at P.A. = 130-140 degree", whose regular variation of LSR velocities with position presents evidence for some ordered motion. We show that the 3-dimensional velocities of this "linear distribution" of 12GHz masers can be well fitted considering a flat, rotating disk, seen almost edge-on.

preprint2010arXiv

Robust zero-averaged wave-number gap inside gapped graphene superlattices

In this paper, the electronic band structures and its transport properties in the gapped graphene superlattices, with one-dimensional (1D) periodic potentials of square barriers, are systematically investigated. It is found that a zero averaged wave-number (zero-$\overline{k}$ ) gap is formed inside the gapped graphene-based superlattices, and the condition for obtaining such a zero-$\overline{k}$ gap is analytically presented. The properties of this zero-$\overline{k}$ gap including its transmission, conductance and Fano factor are studied in detail. Finally it is revealed that the properties of the electronic transmission, conductance and Fano factor near the zero-$\overline{k}$ gap are very insensitive to the structural disorder for the finite graphene-based periodic-barrier systems.

preprint2010arXiv

Shortcut to adiabatic passage in two and three level atoms

We propose a method to transfer the population and control the state of two-level and three-level atoms speeding-up Adiabatic Passage techniques while keeping their robustness versus parameter variations. The method is based on supplementing the standard laser beam setup of Adiabatic Passage methods with auxiliary steering laser pulses of orthogonal polarization. This provides a shortcut to adiabaticity driving the system along the adiabatic path defined by the standard setup.

preprint2010arXiv

Spin transfer in a ferromagnet-quantum dot and tunnel barrier coupled Aharonov-Bohm ring system with Rashba spin-orbit interactions

The spin transfer effect in ferromagnet-quantum dot (insulator)-ferromagnet Aharonov-Bohm (AB) ring system with Rashba spin-orbit (SO) interactions is investigated by means of Keldysh nonequilibrium Green function method. It is found that both the magnitude and direction of the spin transfer torque (STT) acting on the right ferromagnet electrode can be effectively controlled by changing the magnetic flux threading the AB ring or the gate voltage on the quantum dot. The STT can be greatly augmented by matching a proper magnetic flux and an SO interaction at a cost of low electrical current. The STT, electrical current, and spin current are uncovered to oscillate with the magnetic flux. The present results are expected to be useful for information storage in nanospintronics.

preprint2010arXiv

Transient energy excitation in shortcuts to adiabaticity for the time dependent harmonic oscillator

There is recently a surge of interest to cut down the time it takes to change the state of a quantum system adiabatically. We study for the time-dependent harmonic oscillator the transient energy excitation in speed-up processes designed to reproduce the initial populations at some predetermined final frequency and time, providing lower bounds and examples. Implications for the limits imposed to the process times and for the principle of unattainability of the absolute zero, in a single expansion or in quantum refrigerator cycles, are drawn.

preprint2009arXiv

A Search for Infall Evidence in EGOs I: the Northern Sample

We report the first systematic survey of molecular lines (including HCO+ (1-0) and 12CO, 13CO, C18O (1-0) lines at 3 mm band) towards a new sample of 88 massive young stellar object (MYSO) candidates associated with ongoing outflows (known as extended green objects or EGOs) identified from the Spitzer GLIMPSE survey in the northern hemisphere with the PMO-13.7 m radio telescope. By analyzing the asymmetries of the optically thick line HCO+ for 69 of 72 EGOs with HCO+ detection, we found 29 sources with blue asymmetric profiles and 19 sources with red asymmetric profiles. This results in a blue excess of 0.14, seen as a signature of collapsing cores in the observed EGO sample. The relatively small blue excess measured in our full sample due to that the observed EGOs are mostly dominated by outflows and at an earlier evolutionary phase associated with IRDCs and 6.7 GHz methanol masers. The physical properties of clouds surrounding EGOs derived from CO lines are similar to those of massive clumps wherein the massive star forming cores associated with EGOs possibly embedded. The infall velocities and mass infall rates derived for 20 infall candidates are also consistent with the typical values found in MYSOs. Thus our observations further support the speculation of Cyganowski et al. (2008) that EGOs trace a population with ongoing outflow activity and active rapid accretion stage of massive protostellar evolution from a statistical view, although there maybe have limitations due to single-pointing survey with a large beam.

preprint2009arXiv

Experimental demonstration of the topological surface states protected by the time-reversal symmetry

We report direct imaging of standing waves of the nontrivial surface states of topological insulator Bi$_2$Te$_3$ by using a low temperature scanning tunneling microscope. The interference fringes are caused by the scattering of the topological states off Ag impurities and step edges on the Bi$_2$Te$_3$(111) surface. By studying the voltage-dependent standing wave patterns, we determine the energy dispersion $E(k)$, which confirms the Dirac cone structure of the topological states. We further show that, very different from the conventional surface states, the backscattering of the topological states by nonmagnetic impurities is completely suppressed. The absence of backscattering is a spectacular manifestation of the time-reversal symmetry, which offers a direct proof of the topological nature of the surface states.

preprint2009arXiv

Extragalactic Point Source Search in Five-year WMAP 41, 61 and 94 GHz Maps

We present the results of an extragalactic point source search using the five-year WMAP 41, 61 and 94 GHz (Q-, V- and W-band) temperature maps. This work is an extension of our point source search in the WMAP maps applying a CMB-free technique. An internal linear combination (ILC) map has been formed from the three-band maps, with the weights chosen to remove the CMB anisotropy signal as well as to favor a selection of flat-spectrum sources. We find 381 sources at the > 5 sigma level outside the WMAP point source detection mask in the ILC map, among which 89 are "new" (i.e., not present in the WMAP catalog). Source fluxes have been calculated and corrected for the Eddington bias. We have solidly identified 367 (96.3%) of our sources. The 1 sigma positional uncertainty is estimated to be 2'. The 14 unidentified sources could be either extended radio structure or obscured by Galactic emission. We have also applied the same detection process on simulated maps and found 364+/-21 detections on average. The recovered source distribution N(>S) agrees well with the simulation input, which proves the reliability of this method.

preprint2009arXiv

Fast optimal frictionless atom cooling in harmonic traps

A method is proposed to cool down atoms in a harmonic trap without phase-space compression as in a perfectly slow adiabatic expansion, i.e., keeping the populations of the instantaneous initial and final levels invariant, but in a much shorter time. This may require that the harmonic trap becomes an expulsive parabolic potential in some time interval. The cooling times achieved are also shorter than previous minimal times using optimal-control bang-bang methods and real frequencies.

preprint2009arXiv

Frictionless dynamics of Bose-Einstein condensates under fast trap variations

A method is proposed to design the time dependence of the trap frequency and achieve in a short time an adiabatic-like (frictionless) evolution of Bose-Einstein condensates governed by the Gross-Pitaevskii equation. Different cases depending on the effective dimension of the trap and the interaction regimes are considered. 2D traps are particularly suitable as the method can be applied without the need to impose any additional time-dependent change in the strength of the interatomic interaction or a Thomas-Fermi regime as it occurs for 1D and 3D traps.

preprint2008arXiv

Unusual giant magnetoresistance effect in heterojunction structure of ultra-thin single-crystal Pb film on silicon substrate

Superconductor films on semiconductor substrates draw much attention recently since the derived superconductor-based electronics have been shown promising for future data process and storage technologies. By growing atomically uniform single-crystal epitaxial Pb films of several nanometers thick on Si wafers to form a sharp superconductor-semiconductor heterojunction, we have obtained an unusual giant magnetoresistance effect when the Pb film is superconducting. In addition to the great fundamental interest of this effect, the simple structure and compatibility and scalability with current Si-based semiconductor technology offer a great opportunity for integrating superconducting circuits and detectors in a single chip.

preprint2007arXiv

Assertion-Based Design Exploration of DVS in Network Processor Architectures

With the scaling of technology and higher requirements on performance and functionality, power dissipation is becoming one of the major design considerations in the development of network processors. In this paper, we use an assertion-based methodology for system-level power/performance analysis to study two dynamic voltage scaling (DVS) techniques, traffic-based DVS and execution-based DVS, in a network processor model. Using the automatically generated distribution analyzers, we analyze the power and performance distributions and study their trade-offs for the two DVS policies with different parameter settings such as threshold values and window sizes. We discuss the optimal configurations of the two DVS policies under different design requirements. By a set of experiments, we show that the assertion-based trace analysis methodology is an efficient tool that can help a designer easily compare and study optimal architectural configurations in a large design space.