Catalog footprint

What is connected

89works
53topics
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

89 published item(s)

preprint2026arXiv

SkillOS: Learning Skill Curation for Self-Evolving Agents

LLM-based agents are increasingly deployed to handle streaming tasks, yet they often remain one-off problem solvers that fail to learn from past interactions. Reusable skills distilled from experience provide a natural substrate for self-evolution, where high-quality skill curation serves as the key bottleneck. Existing approaches either rely on manual skill curation, prescribe heuristic skill operations, or train for short-horizon skill operations. However, they still struggle to learn complex long-term curation policies from indirect and delayed feedback. To tackle this challenge, we propose SkillOS, an experience-driven RL training recipe for learning skill curation in self-evolving agents. SkillOS pairs a frozen agent executor that retrieves and applies skills with a trainable skill curator that updates an external SkillRepo from accumulated experience. To provide learning signals for curation, we design composite rewards and train on grouped task streams based on skill-relevant task dependencies, where earlier trajectories update the SkillRepo, and later related tasks evaluate these updates. Across multi-turn agentic tasks and single-turn reasoning tasks, SkillOS consistently outperforms memory-free and strong memory-based baselines in both effectiveness and efficiency, with the learned skill curator generalizing across different executor backbones and task domains. Further analyses show that the learned curator produces more targeted skill use, while the skills in SkillRepo evolve into more richly structured Markdown files that encode higher-level meta-skills over time.

preprint2024arXiv

Exception-aware Lifecycle Model Construction for Framework APIs

The implementation of complex software systems usually depends on low-level frameworks or third-party libraries. During their evolution, the APIs adding and removing behaviors may cause unexpected compatibility problems. So, precisely analyzing and constructing the framework/ library's API lifecycle model is of great importance. Existing works have proposed the API existence-changing model for defect detection, while not considering the influence of semantic changes in APIs. In some cases, developers will not remove or deprecate APIs but modify their semantics by adding, removing, or modifying their exception-thrown code, which may bring potential defects to upper-level code. Therefore, besides the API existence model, it is also necessary for developers to be concerned with the exception-related code evolution in APIs, which requires the construction of exception-aware API lifecycle models for framework/library projects. To achieve automatic exception-aware API lifecycle model construction, this paper adopts a static analysis technique to extract exception summary information in the framework API code and adopts a multi-step matching strategy to obtain the changing process of exceptions. Then, it generates exception-aware API lifecycle models for the given framework/library project. With this approach, the API lifecycle extraction tool, JavaExP, is implemented, which is based on Java bytecode analysis. Compared to the state-of-the-art tool, JavaExP achieves both a higher F1 score (+60%) and efficiency (+7x), whose precision of exception matching and changing results is 98%. Compared to the exception-unaware API lifecycle modeling on 60 versions, JavaExp can identify 18% times more API changes. Among the 75,433 APIs under analysis, 20% of APIs have changed their exception-throwing behavior at least once after API introduction, which may bring many hidden compatibility issues.

preprint2024arXiv

Multi-agent Reinforcement Learning for Cooperative Lane Changing of Connected and Autonomous Vehicles in Mixed Traffic

Autonomous driving has attracted significant research interests in the past two decades as it offers many potential benefits, including releasing drivers from exhausting driving and mitigating traffic congestion, among others. Despite promising progress, lane-changing remains a great challenge for autonomous vehicles (AV), especially in mixed and dynamic traffic scenarios. Recently, reinforcement learning (RL), a powerful data-driven control method, has been widely explored for lane-changing decision makings in AVs with encouraging results demonstrated. However, the majority of those studies are focused on a single-vehicle setting, and lane-changing in the context of multiple AVs coexisting with human-driven vehicles (HDVs) have received scarce attention. In this paper, we formulate the lane-changing decision making of multiple AVs in a mixed-traffic highway environment as a multi-agent reinforcement learning (MARL) problem, where each AV makes lane-changing decisions based on the motions of both neighboring AVs and HDVs. Specifically, a multi-agent advantage actor-critic network (MA2C) is developed with a novel local reward design and a parameter sharing scheme. In particular, a multi-objective reward function is proposed to incorporate fuel efficiency, driving comfort, and safety of autonomous driving. Comprehensive experimental results, conducted under three different traffic densities and various levels of human driver aggressiveness, show that our proposed MARL framework consistently outperforms several state-of-the-art benchmarks in terms of efficiency, safety and driver comfort.

preprint2022arXiv

A Comprehensive Evaluation of Android ICC Resolution Techniques

Inter-component communication (ICC) is a widely used mechanism in mobile apps, which enables message-based control flow transferring and data passing between Android components. Effective ICC resolution requires precisely identifying entry points, analyzing data values of ICC fields, modeling related framework APIs, etc. Due to various control-flow- and data-flow-related characteristics involved and the lack of oracles for real-world apps, the comprehensive evaluation of ICC resolution techniques is challenging. To fill this gap, we collect multiple-type benchmark suites with 4,104 apps, covering hand-made apps, open-source, and commercial ones. Considering their differences, various evaluation metrics, e.g., number count, graph structure, and reliable oracle based metrics, are adopted on-demand. As the oracle for real-world apps is unavailable, we design a dynamic analysis approach to extract the real ICC links triggered during GUI exploration. By auditing the code implementations, we carefully check the extracted ICCs and confirm 1,680 ones to form a reliable oracle set, in which each ICC is labeled with 25 code characteristic tags. The evaluation performed on six state-of-the-art ICC resolution tools shows that 1) the completeness of static ICC resolution results on real-world apps is not satisfactory, as up to 38%-85% ICCs are missed by tools; 2) many wrongly reported ICCs are sent from or received by only a few components and the graph structure information can help the identification; 3) the efficiency of fundamental tools, like ICC resolution ones, should be optimized in both engineering and research aspects. By investigating both the missed and wrongly reported ICCs, we discuss the strengths of different tools for users and summarize eight common FN/FP patterns in ICC resolution for tool developers.

preprint2022arXiv

An Efficient Algorithm for Generating Directed Networks with Predetermined Assortativity Measures

Assortativity coefficients are important metrics to analyze both directed and undirected networks. In general, it is not guaranteed that the fitted model will always agree with the assortativity coefficients in the given network, and the structure of directed networks is more complicated than the undirected ones. Therefore, we provide a remedy by proposing a degree-preserving rewiring algorithm, called DiDPR, for generating directed networks with given directed assortativity coefficients. We construct the joint edge distribution of the target network by accounting for the four directed assortativity coefficients simultaneously, provided that they are attainable, and obtain the desired network by solving a convex optimization problem.Our algorithm also helps check the attainability of the given assortativity coefficients. We assess the performance of the proposed algorithm by simulation studies with focus on two different network models, namely Erdös--Rényi and preferential attachment random networks. We then apply the algorithm to a Facebook wall post network as a real data example. The codes for implementing our algorithm are publicly available in R package wdnet.

preprint2022arXiv

AUGER: Automatically Generating Review Comments with Pre-training Models

Code review is one of the best practices as a powerful safeguard for software quality. In practice, senior or highly skilled reviewers inspect source code and provide constructive comments, considering what authors may ignore, for example, some special cases. The collaborative validation between contributors results in code being highly qualified and less chance of bugs. However, since personal knowledge is limited and varies, the efficiency and effectiveness of code review practice are worthy of further improvement. In fact, it still takes a colossal and time-consuming effort to deliver useful review comments. This paper explores a synergy of multiple practical review comments to enhance code review and proposes AUGER (AUtomatically GEnerating Review comments): a review comments generator with pre-training models. We first collect empirical review data from 11 notable Java projects and construct a dataset of 10,882 code changes. By leveraging Text-to-Text Transfer Transformer (T5) models, the framework synthesizes valuable knowledge in the training stage and effectively outperforms baselines by 37.38% in ROUGE-L. 29% of our automatic review comments are considered useful according to prior studies. The inference generates just in 20 seconds and is also open to training further. Moreover, the performance also gets improved when thoroughly analyzed in case study.

preprint2022arXiv

Density-Aware Hyper-Graph Neural Networks for Graph-based Semi-supervised Node Classification

Graph-based semi-supervised learning, which can exploit the connectivity relationship between labeled and unlabeled data, has been shown to outperform the state-of-the-art in many artificial intelligence applications. One of the most challenging problems for graph-based semi-supervised node classification is how to use the implicit information among various data to improve the performance of classifying. Traditional studies on graph-based semi-supervised learning have focused on the pairwise connections among data. However, the data correlation in real applications could be beyond pairwise and more complicated. The density information has been demonstrated to be an important clue, but it is rarely explored in depth among existing graph-based semi-supervised node classification methods. To develop a flexible and effective model for graph-based semi-supervised node classification, we propose a novel Density-Aware Hyper-Graph Neural Networks (DA-HGNN). In our proposed approach, hyper-graph is provided to explore the high-order semantic correlation among data, and a density-aware hyper-graph attention network is presented to explore the high-order connection relationship. Extensive experiments are conducted in various benchmark datasets, and the results demonstrate the effectiveness of the proposed approach.

preprint2022arXiv

Ergodic problems for contact Hamilton-Jacobi equations

This paper deals with the generalized ergodic problem \[ H(x,u(x),Du(x))=c, \quad x\in M, \] where the unknown is a pair $(c,u)$ of a constant $c \in \mathbb{R}$ and a function $u$ on $M$ for which $u$ is a viscosity solution. We assume $H=H(x,u,p)$ satisfies Tonelli conditions in the argument $p\in T^*_xM$ and the Lipschitz condition in the argument $u\in\R$. For a given $c\in \R$, we first discuss necessary and sufficient conditions for the existence of viscosity solutions. Let $\mathfrak{C}$ denote the set of all real numbers $c$'s for which the above equation admits viscosity solutions. Then we show $\mathfrak{C}$ is an interval, whose endpoints $\x$, $\y$ with $\x\leqslant\y$ can be characterized by a min-max formula and a max-min formula, respectively. The most significant finding is that we figure out the structure of $\mathfrak{C}$ without monotonicity assumptions on $u$.

preprint2022arXiv

FedSSO: A Federated Server-Side Second-Order Optimization Algorithm

In this work, we propose FedSSO, a server-side second-order optimization method for federated learning (FL). In contrast to previous works in this direction, we employ a server-side approximation for the Quasi-Newton method without requiring any training data from the clients. In this way, we not only shift the computation burden from clients to server, but also eliminate the additional communication for second-order updates between clients and server entirely. We provide theoretical guarantee for convergence of our novel method, and empirically demonstrate our fast convergence and communication savings in both convex and non-convex settings.

preprint2022arXiv

Identifying structure-absorption relationships and predicting absorption strength of non-fullerene acceptors for organic photovoltaics

Non-fullerene acceptors (NFAs) are excellent light harvesters, yet the origin of such high optical extinction is not well understood. In this work, we investigate the absorption strength of NFAs by building a database of time-dependent density functional theory (TDDFT) calculations of ~500 pi-conjugated molecules. The calculations are first validated by comparison with experimental measurements on liquid and solid state using common fullerene and non-fullerene acceptors. We find that the molar extinction coefficient (ε_(d,max)) shows reasonable agreement between calculation in vacuum and experiment for molecules in solution, highlighting the effectiveness of TDDFT for predicting optical properties of organic pi-conjugated molecules. We then perform a statistical analysis based on molecular descriptors to identify which features are important in defining the absorption strength. This allows us to identify structural features that are correlated with high absorption strength in NFAs and could be used to guide molecular design: highly absorbing NFAs should possess a planar, linear, and fully conjugated molecular backbone with highly polarisable heteroatoms. We then exploit a random decision forest to draw predictions for ε_(d,max) using a computational framework based on extended tight-binding Hamiltonians, which shows reasonable predicting accuracy with lower computational cost than TDDFT. This work provides a general understanding of the relationship between molecular structure and absorption strength in pi-conjugated organic molecules, including NFAs, while introducing predictive machine-learning models of low computational cost.

preprint2022arXiv

Improving Tese Case Generation for Python Native Libraries Through Constraints on Input Data Structures

Modern Python projects execute computational functions using native libraries and give Python interfaces to boost execution speed; hence, testing these libraries becomes critical to the project's robustness. One challenge is that existing approaches use coverage to guide generation, but native libraries run as black boxes to Python code with no execution information. Another is that dynamic binary instrumentation reduces testing performance as it needs to monitor both native libraries and the Python virtual machine. To address these challenges, in this paper, we propose an automated test case generation approach that works at the Python code layer. Our insight is that many path conditions in native libraries are for processing input data structures through interacting with the VM. In our approach, we instrument the Python Interpreter to monitor the interactions between native libraries and VM, derive constraints on the structures, and then use the constraints to guide test case generation. We implement our approach in a tool named PyCing and apply it to six widely-used Python projects. The experimental results reveal that with the structure constraint guidance, PyCing can cover more execution paths than existing test cases and state-of-the-art tools. Also, with the checkers in the testing framework Pytest, PyCing can identify segmentation faults in 10 Python interfaces and memory leaks in 9. Our instrumentation strategy also has an acceptable influence on testing efficiency.

preprint2022arXiv

Knowledge Management for Cloud Computing Field

Migration legacy systems to cloud platforms is a knowledge intensive process. There is an ever increasing body of knowledge reporting empirical scenarios of successful and problematic cloud migration. Reusing this body of knowledge, dispersed and fragmented over the academic/multi-vocal literature, has practical values to mitigate costly risks and pitfalls in further projects of legacy to-cloud and cloud-to-cloud migration. In line with this, knowledge management systems/platforms pertinent to cloud migration are a prime prerequisite and a strategic imperative for an organization. We have conducted a qualitative exploratory study to understand the benefits and challenges of developing Knowledge Management Systems (KMS) for cloud migration in real trials. Whilst our prototype system demonstration supported the importance and bene-fits of developing Cloud Migration KMS (CM-KMS), our semi-structured industry interview study with 11 participants highlighted challenging impediments against developing this class of KMS. As a result, this study proposes nine significant challenges that cause the abandon of the design and maintenance of CM-KMS, including continuous changes and updates, integration of knowledge, knowledge granularity, preservation of context, automation, deconstruction of traditional knowledge, dependency on experts, hybrid knowledge of both vendor-specific and vendor-neutral cloud platforms, and parsimony. Our results inform cloud architects to pay attention to adopt CM-KMS for the legacy-to-cloud migration in their organizations.

preprint2022arXiv

Nonlinear semigroup approach to Hamilton-Jacobi equations -- A toy model

In this paper, we discuss the existence and multiplicity problem of viscosity solution to the Hamilton-Jacobi equation $$h(x,d_x u)+λ(x)u=c,\quad x\in M,$$ where $M$ is a closed manifold and $λ:M\rightarrow\mathbb{R}$ changes signs on $M$, via nonlinear semigroup method. It turns out that a bifurcation phenomenon occurs when parameter $c$ strides over the critical value. As an application of the main result, we analyse the structure of the set of viscosity solutions of an one-dimensional example in detail.

preprint2022arXiv

On the Minimax Spherical Designs

Distributing points on a (possibly high-dimensional) sphere with minimal energy is a long-standing problem in and outside the field of mathematics. This paper considers a novel energy function that arises naturally from statistics and combinatorial optimization, and studies its theoretical properties. Our result solves both the exact optimal spherical point configurations in certain cases and the minimal energy asymptotics under general assumptions. Connections between our results and the L1-Principal Component analysis and Quasi-Monte Carlo methods are also discussed.

preprint2022arXiv

On the Robustness of Reading Comprehension Models to Entity Renaming

We study the robustness of machine reading comprehension (MRC) models to entity renaming -- do models make more wrong predictions when the same questions are asked about an entity whose name has been changed? Such failures imply that models overly rely on entity information to answer questions, and thus may generalize poorly when facts about the world change or questions are asked about novel entities. To systematically audit this issue, we present a pipeline to automatically generate test examples at scale, by replacing entity names in the original test sample with names from a variety of sources, ranging from names in the same test set, to common names in life, to arbitrary strings. Across five datasets and three pretrained model architectures, MRC models consistently perform worse when entities are renamed, with particularly large accuracy drops on datasets constructed via distant supervision. We also find large differences between models: SpanBERT, which is pretrained with span-level masking, is more robust than RoBERTa, despite having similar accuracy on unperturbed test data. We further experiment with different masking strategies as the continual pretraining objective and find that entity-based masking can improve the robustness of MRC models.

preprint2022arXiv

Regression Modeling for Recurrent Events Using R Package reReg

Recurrent event analyses have found a wide range of applications in biomedicine, public health, and engineering, among others, where study subjects may experience a sequence of event of interest during follow-up. The R package reReg (Chiou and Huang 2021) offers a comprehensive collection of practical and easy-to-use tools for regression analysis of recurrent events, possibly with the presence of an informative terminal event. The regression framework is a general scale-change model which encompasses the popular Cox-type model, the accelerated rate model, and the accelerated mean model as special cases. Informative censoring is accommodated through a subject-specific frailty without no need for parametric specification. Different regression models are allowed for the recurrent event process and the terminal event. Also included are visualization and simulation tools.

preprint2022arXiv

Wavelet Regularization Benefits Adversarial Training

Adversarial training methods are state-of-the-art (SOTA) empirical defense methods against adversarial examples. Many regularization methods have been proven to be effective with the combination of adversarial training. Nevertheless, such regularization methods are implemented in the time domain. Since adversarial vulnerability can be regarded as a high-frequency phenomenon, it is essential to regulate the adversarially-trained neural network models in the frequency domain. Faced with these challenges, we make a theoretical analysis on the regularization property of wavelets which can enhance adversarial training. We propose a wavelet regularization method based on the Haar wavelet decomposition which is named Wavelet Average Pooling. This wavelet regularization module is integrated into the wide residual neural network so that a new WideWaveletResNet model is formed. On the datasets of CIFAR-10 and CIFAR-100, our proposed Adversarial Wavelet Training method realizes considerable robustness under different types of attacks. It verifies the assumption that our wavelet regularization method can enhance adversarial robustness especially in the deep wide neural networks. The visualization experiments of the Frequency Principle (F-Principle) and interpretability are implemented to show the effectiveness of our method. A detailed comparison based on different wavelet base functions is presented. The code is available at the repository: \url{https://github.com/momo1986/AdversarialWaveletTraining}.

preprint2021arXiv

Assortativity measures for weighted and directed networks

Assortativity measures the tendency of a vertex in a network being connected by other vertexes with respect to some vertex-specific features. Classical assortativity coefficients are defined for unweighted and undirected networks with respect to vertex degree. We propose a class of assortativity coefficients that capture the assortative characteristics and structure of weighted and directed networks more precisely. The vertex-to-vertex strength correlation is used as an example, but the proposed measure can be applied to any pair of vertex-specific features. The effectiveness of the proposed measure is assessed through extensive simulations based on prevalent random network models in comparison with existing assortativity measures. In application World Input-Ouput Networks,the new measures reveal interesting insights that would not be obtained by using existing ones. An implementation is publicly available in a R package "wdnet".

preprint2021arXiv

Atomic Norm Based Localization of Far-Field and Near-Field Signals with Generalized Symmetric Arrays

Most localization methods for mixed far-field (FF) and near-field (NF) sources are based on uniform linear array (ULA) rather than sparse linear array (SLA). In this paper, we propose a localization method for mixed FF and NF sources based on the generalized symmetric linear arrays, which include ULAs, Cantor array, Fractal array and many other SLAs. Our method consists of two steps. In the first step, the high-order statistics of the array output is exploited to increase the degree of freedom. Then the direction-of-arrivals (DOAs) of the FF and NF sources are jointly estimated by using the recently proposed atomic norm minimization (ANM), which belongs to the gridless super-resolution method since the discretization of the parameter space is not required. In the second step, the ranges are given by MUSIC-like one-dimensional searching. Simulations results are provided to demonstrate the advantages of our method.

preprint2021arXiv

Axial Higgs Mode Detected by Quantum Pathway Interference in RTe3

The observation of the Higgs boson solidified the standard model of particle physics. However, explanations of anomalies (e.g. dark matter) rely on further symmetry breaking calling for an undiscovered axial Higgs mode. In condensed matter the Higgs was seen in magnetic, superconducting and charge density wave(CDW) systems. Uncovering a low energy mode's vector properties is challenging, requiring going beyond typical spectroscopic or scattering techniques. Here, we discover an axial Higgs mode in the CDW system RTe3 using the interference of quantum pathways. In RTe3 (R=La,Gd), the electronic ordering couples bands of equal or different angular momenta. As such, the Raman scattering tensor associated to the Higgs mode contains both symmetric and antisymmetric components, which can be excited via two distinct, but degenerate pathways. This leads to constructive or destructive interference of these pathways, depending on the choice of the incident and Raman scattered light polarization. The qualitative behavior of the Raman spectra is well-captured by an appropriate tight-binding model including an axial Higgs mode. The elucidation of the antisymmetric component provides direct evidence that the Higgs mode contains an axial vector representation (i.e. a pseudo-angular momentum) and hints the CDW in RTe3 is unconventional. Thus we provide a means for measuring collective modes quantum properties without resorting to extreme experimental conditions.

preprint2021arXiv

Clustering US States by Time Series of COVID-19 New Case Counts with Non-negative Matrix Factorization

The spreading pattern of COVID-19 differ a lot across the US states under different quarantine measures and reopening policies. We proposed to cluster the US states into distinct communities based on the daily new confirmed case counts via a nonnegative matrix factorization (NMF) followed by a k-means clustering procedure on the coefficients of the NMF basis. A cross-validation method was employed to select the rank of the NMF. Applying the method to the entire study period from March 22 to July 25, we clustered the 49 continental states (including District of Columbia) into 7 groups, two of which contained a single state. To investigate the dynamics of the clustering results over time, the same method was successively applied to the time periods with increment of one week, starting from the period of March 22 to March 28. The results suggested a change point in the clustering in the week starting on May 30, which might be explained by a combined impact of both quarantine measures and reopening policies.

preprint2021arXiv

Extended calculations of energy levels, radiative properties, and lifetimes for P-like Ge XVIII

Using the multiconfiguration Dirac-Hartree-Fock (MCDHF) and the relativistic configuration interaction (RCI) methods, a consistent set of transition energies and radiative transition data for the lowest 150 states of the $3s^2 3p^3$, $3s 3p^4$, $3s^2 3p^2 3d$, $3s 3p^3 3d$, $3p^5$, and $3s^2 3p 3d^2$ configurations in P-like Ge XVIII is provided. To assess the accuracy of the MCDHF transition energies, we have also performed calculations using the many-body perturbation theory (MBPT). Comparisons are made between the present MCDHF and MBPT data sets, as well as with other available experimental and theoretical values. The resulting accurate and consistent MCDHF data set will be useful for line identification and modeling purposes. These data can be considered as a benchmark for other calculations.

preprint2021arXiv

Time-periodic solutions of contact Hamilton-Jacobi equations on the circle

We are concerned with the existence and multiplicity of nontrivial time-periodic viscosity solutions to \[ \partial_t w(x,t) + H( x,\partial_x w(x,t),w(x,t) )=0,\quad (x,t)\in \mathbb{S} \times [0,+\infty). \] We find that there are infinitely many nontrivial time-periodic viscosity solutions with different periods when $\frac{\partial H}{\partial u}(x,p,u)\leqslant-δ<0$ by analyzing the asymptotic behavior of the dynamical system $(C(\mathbb{S} ,\mathbb{R}),\{T_t\}_{t\geqslant 0})$, where $\{T_t\}_{t\geqslant 0}$ was introduced in \cite{WWY1}. Moreover, in view of the convergence of $T_{t_n}φ$, we get the existence of nontrivial periodic points of $T_t$, where $φ$ are initial data satisfying certain properties. This is a long-time behavior result for the solution to the above equation with initial data $φ$. At last, as an application, we describe to readers a bifurcation phenomenon for \[ \partial_t w(x,t) + H( x,\partial_x w(x,t),λw(x,t) )=0,\quad (x,t)\in \mathbb{S} \times [0,+\infty), \] when the sign of the parameter $λ$ varies. The structure of the unit circle $\mathbb{S}$ plays an essential role here. The most important novelty is the discovery of the nontrivial recurrence of $(C(\mathbb{S} ,\mathbb{R}),\{T_t\}_{t\geqslant 0})$.

preprint2020arXiv

A Fast Radio Burst discovered in FAST drift scan survey

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

preprint2020arXiv

Computational Performance of a Germline Variant Calling Pipeline for Next Generation Sequencing

With the booming of next generation sequencing technology and its implementation in clinical practice and life science research, the need for faster and more efficient data analysis methods becomes pressing in the field of sequencing. Here we report on the evaluation of an optimized germline mutation calling pipeline, HummingBird, by assessing its performance against the widely accepted BWA-GATK pipeline. We found that the HummingBird pipeline can significantly reduce the running time of the primary data analysis for whole genome sequencing and whole exome sequencing while without significantly sacrificing the variant calling accuracy. Thus, we conclude that expansion of such software usage will help to improve the primary data analysis efficiency for next generation sequencing.

preprint2020arXiv

Discovery and timing of pulsars in the globular cluster M13 with FAST

We report the discovery of a binary millisecond pulsar (namely PSR J1641+3627F or M13F) in the globular cluster M13 (NGC 6205) and timing solutions of M13A to F using observations made with the Five-hundred-metre Aperture Spherical radio Telescope (FAST). PSR J1641+3627F has a spin period of 3.00 ms and an orbital period of 1.4 days. The most likely companion mass is 0.16 M$_{\odot}$. M13A to E all have short spin periods and small period derivatives. We also confirm that the binary millisecond pulsar PSR J1641$+$3627E (also M13E) is a black widow with a companion mass around 0.02 M$_{\odot}$. We find that all the binary systems have low eccentricities compared to those typical for globular cluster pulsars and that they decrease with distance from the cluster core. This is consistent with what is expected as this cluster has a very low encounter rate per binary.

preprint2020arXiv

Enabling Cyberattack-Resilient Load Forecasting through Adversarial Machine Learning

In the face of an increasingly broad cyberattack surface, cyberattack-resilient load forecasting for electric utilities is both more necessary and more challenging than ever. In this paper, we propose an adversarial machine learning (AML) approach, which can respond to a wide range of attack behaviors without detecting outliers. It strikes a balance between enhancing a system's robustness against cyberattacks and maintaining a reasonable degree of forecasting accuracy when there is no attack. Attack models and configurations for the adversarial training were selected and evaluated to achieve the desired level of performance in a simulation study. The results validate the effectiveness and excellent performance of the proposed method.

preprint2020arXiv

Ground and Excited Exciton Polarons in Monolayer MoSe2

Monolayer transition metal dichalcogenide semiconductors, with versatile experimentally accessible exciton species, offer an interesting platform for investigating the interaction between excitons and a Fermi sea of charges. Using hexagonal boron nitride encapsulated monolayer MoSe2, we study the impact of charge density tuning on the ground and excited Rydberg states in the atomic layer. Consistent excitonpolaron behavior is revealed in both photoluminescence and reflection spectra of the A exciton 1s (A:1s) Rydberg state, in contrast to previous studies. The A:2s Rydberg state provides an opportunity to understand such interactions with greatly reduced exciton binding energy. We found that the impact of the Fermi sea becomes much more dramatic. With a photoluminescence upconversion technique, we further verify the 2s polaron-like behavior for the repulsive branch of B:2s exciton whose energy is well above the bare bandgap. Our studies show that the polaron-like interaction features are quite generic and highly robust, offering key insights into the dressed manybody state in a Fermi sea.

preprint2020arXiv

Learning from Explanations with Neural Execution Tree

While deep neural networks have achieved impressive performance on a range of NLP tasks, these data-hungry models heavily rely on labeled data, which restricts their applications in scenarios where data annotation is expensive. Natural language (NL) explanations have been demonstrated very useful additional supervision, which can provide sufficient domain knowledge for generating more labeled data over new instances, while the annotation time only doubles. However, directly applying them for augmenting model learning encounters two challenges: (1) NL explanations are unstructured and inherently compositional, which asks for a modularized model to represent their semantics, (2) NL explanations often have large numbers of linguistic variants, resulting in low recall and limited generalization ability. In this paper, we propose a novel Neural Execution Tree (NExT) framework to augment training data for text classification using NL explanations. After transforming NL explanations into executable logical forms by semantic parsing, NExT generalizes different types of actions specified by the logical forms for labeling data instances, which substantially increases the coverage of each NL explanation. Experiments on two NLP tasks (relation extraction and sentiment analysis) demonstrate its superiority over baseline methods. Its extension to multi-hop question answering achieves performance gain with light annotation effort.

preprint2020arXiv

Long-lived populations of momentum- and spin-indirect excitons in monolayer WSe$_2$

Monolayer transition metal dichalcogenides are a promising platform to investigate many-body interactions of excitonic complexes. In monolayer tungsten diselenide, the ground-state exciton is dark (spin-indirect), and the valley degeneracy allows low-energy dark momentum-indirect excitons to form. Interactions between the dark exciton species and the optically accessible bright exciton (X) are likely to play significant roles in determining the optical properties of X at high power, as well as limiting the ultimate exciton densities that can be achieved, yet so far little is known about these interactions. Here, we demonstrate long-lived dense populations of momentum-indirect intervalley ($X_K$) and spin-indirect intravalley (D) dark excitons by time-resolved photoluminescence measurements (Tr-PL). Our results uncover an efficient inter-state conversion between X to D excitons through the spin-flip process and the one between D and $X_K$ excitons mediated by the exchange interaction (D + D to $X_K$ + $X_K$). Moreover, we observe a persistent redshift of the X exciton due to strong excitonic screening by $X_K$ exciton with a response time in the timescale of sub-ns, revealing a non-trivial inter-state exciton-exciton interaction. Our results provide a new insight into the interaction between bright and dark excitons, and point to a possibility to employ dark excitons for investigating exciton condensation and the valleytronics.

preprint2020arXiv

Multivariate Distributionally Robust Convex Regression under Absolute Error Loss

This paper proposes a novel non-parametric multidimensional convex regression estimator which is designed to be robust to adversarial perturbations in the empirical measure. We minimize over convex functions the maximum (over Wasserstein perturbations of the empirical measure) of the absolute regression errors. The inner maximization is solved in closed form resulting in a regularization penalty involves the norm of the gradient. We show consistency of our estimator and a rate of convergence of order $ \widetilde{O}\left( n^{-1/d}\right) $, matching the bounds of alternative estimators based on square-loss minimization. Contrary to all of the existing results, our convergence rates hold without imposing compactness on the underlying domain and with no a priori bounds on the underlying convex function or its gradient norm.

preprint2020arXiv

Orbital angular momentum beam generation using a free-electron laser oscillator

With wavelength tunability, free-electron lasers (FELs) are well-suited for generating orbital angular momentum (OAM) beams in a wide photon energy range. We report the first experimental demonstration of OAM beam generation using an oscillator FEL. Lasing around 458 nm, we have produced the four lowest orders of coherently mixed OAM beams with good beam quality, excellent stability, and substantial intracavity power. We have also developed a pulsed mode operation of the OAM beam with a highly reproducible temporal structure for a range of modulation frequencies from 1 to 30 Hz. This development can be extended to short wavelengths, for example to x-rays using a future x-ray FEL oscillator. The operation of such an OAM FEL also paves the way for the generation of OAM gamma-ray beams via Compton scattering.

preprint2020arXiv

Probing new physics with multi-vacua quantum tunnelings beyond standard model through gravitational waves

We report on a novel phenomenon of particle cosmology, which features specific cosmological phase transitions via quantum tunnelings through multiple vacua. The latter is inspired by the axiverse ideas and enables us to probe the associated new physics models through a potential observation of specific patterns in the stochastic gravitational waves background. Multiple vacua may induce the nucleation of co-existing bubbles over the phase transition epoch, hence enhancing the overall process of bubbles' nucleation. Our detailed analysis of semi-analytical and numerical solutions to the bounce equations of the path integral in three vacua case has enabled us to determine the existence of three instanton solutions. This new mechanism of cosmological phase transitions clearly predicts a possibly sizeable new source of gravitational waves, with its energy spectrum being featured with particular patterns, which could be probed by the future gravitational wave interferometers.

preprint2020arXiv

Radial distribution of charm quarks in jets in high-energy heavy-ion collisions

Heavy flavor physics in high-energy heavy-ion collisions is a promising and active area to study the mass dependence of the "jet quenching" effects both at the RHIC and the LHC. In this talk, we present the first theoretical study on the $D^0$ meson radial distributions relative to the jet axis both in p+p and Pb+Pb collisions at $\sqrt{s_{NN}}=5.02$ TeV, where a nice agreement of our results with experimental data is observed. The in-medium parton propagations are described by a Monte Carlo transport model which uses the next-to-leading order (NLO) plus parton shower (PS) event generator SHERPA as input and includes elastic (collisional) and inelastic (radiative) in-medium interaction of heavy flavor jet. We find that, at low $D^0$ meson $p_T$, the radial distribution significantly shifts to larger radius indicating a strong diffusion effect, and the diffusion effects decrease quickly with $p_T$ ,which is consistent with the recent CMS measurements. We demonstrate that the angular deviation of charm quarks is sensitive to $D_s$ but not $\hat{q}$, which may provide new constrains on the collisional and radiative heavy quark energy loss.

preprint2020arXiv

Remarks on the distribution of colors in Gallai colorings

A Gallai coloring of a complete graph $K_n$ is an edge coloring without triangles colored with three different colors. A sequence $e_1\ge \dots \ge e_k$ of positive integers is an $(n,k)$-sequence if $\sum_{i=1}^k e_i=\binom{n}{2}$. An $(n,k)$-sequence is a G-sequence if there is a Gallai coloring of $K_n$ with $k$ colors such that there are $e_i$ edges of color $i$ for all $i,1\le i \le k$. Gyárfás, Pálvölgyi, Patkós and Wales proved that for any integer $k\ge 3$ there exists an integer $g(k)$ such that every $(n,k)$-sequence is a G-sequence if and only if $n\ge g(k)$. They showed that $g(3)=5, g(4)=8$ and $2k-2\le g(k)\le 8k^2+1$. We show that $g(5)=10$ and give almost matching lower and upper bounds for $g(k)$ by showing that with suitable constants $α,β>0$, $\frac{αk^{1.5}}{\ln k}\le g(k) \le βk^{1.5}$ for all sufficiently large $k$.

preprint2020arXiv

Scalable Multi-Hop Relational Reasoning for Knowledge-Aware Question Answering

Existing work on augmenting question answering (QA) models with external knowledge (e.g., knowledge graphs) either struggle to model multi-hop relations efficiently, or lack transparency into the model's prediction rationale. In this paper, we propose a novel knowledge-aware approach that equips pre-trained language models (PTLMs) with a multi-hop relational reasoning module, named multi-hop graph relation network (MHGRN). It performs multi-hop, multi-relational reasoning over subgraphs extracted from external knowledge graphs. The proposed reasoning module unifies path-based reasoning methods and graph neural networks to achieve better interpretability and scalability. We also empirically show its effectiveness and scalability on CommonsenseQA and OpenbookQA datasets, and interpret its behaviors with case studies.

preprint2020arXiv

The Fundamental Performance of FAST with 19-beam Receiver at L Band

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) passed national acceptance and is taking pilot cycle of 'Shared-Risk' observations. The 19-beam receiver covering 1.05-1.45 GHz was used for most of these observations. The electronics gain fluctuation of the system is better than 1\% over 3.5 hours, enabling enough stability for observations. Pointing accuracy, aperture efficiency and system temperature are three key parameters of FAST. The measured standard deviation of pointing accuracy is 7.9$''$, which satisfies the initial design of FAST. When zenith angle is less than 26.4$^\circ$, the aperture efficiency and system temperature around 1.4 GHz are $\sim$ 0.63 and less than 24 K for central beam, respectively. The measured value of these two parameters are better than designed value of 0.6 and 25 K, respectively. The sensitivity and stability of the 19-beam backend are confirmed to satisfy expectation by spectral HI observations toward N672 and polarization observations toward 3C286. The performance allows FAST to take sensitive observations in various scientific goals, from studies of pulsar to galaxy evolution.

preprint2020arXiv

The range of non-Kitaev terms and fractional particles in RuCl$_3$

Significant efforts have focused on the magnetic excitations of relativistic Mott insulators, predicted to realize the Kitaev quantum spin liquid (QSL). This exactly solvable model involves a highly entangled state resulting from bond-dependent Ising interactions that produce excitations which are non-local in terms of spin flips. A key challenge in real materials is identifying the relative size of the non-Kitaev terms and their role in the emergence or suppression of fractional excitations. Here, we identify the energy and temperature boundaries of non-Kitaev interactions by direct comparison of the Raman susceptibility of RuCl3 with quantum Monte Carlo (QMC) results for the Kitaev QSLs. Moreover, we further confirm the fractional nature of the magnetic excitations, which is given by creating a pair of fermionic quasiparticles. Interestingly, this fermionic response remains valid in the non-Kitaev range. Our results and focus on the use of the Raman susceptibility provide a stringent new test for future theoretical and experimental studies of QSLs.

preprint2020arXiv

Up- and Down-Conversion between Intra- and Inter-Valley Excitons in Waveguide Coupled Monolayer WSe2

The presence of two spin-split valleys in monolayer (1L) transition metal dichalcogenide (TMD) semiconductors supports versatile exciton species classified by their spin and valley quantum numbers. While the spin-0 intravalley exciton, known as the "bright" exciton, is readily observable, other types of excitons, such as the spin-1 intravalley (spin-dark) and spin-0 intervalley (momentum-dark) excitons, are more difficult to access. Here we develop a waveguide coupled 1L tungsten diselenide (WSe2) device to probe these exciton species. In particular, TM coupling to the atomic layer's out-of-plane dipole moments enabled us to not only efficiently collect, but also resonantly populate the spin-1 dark excitons, promising for developing devices with long valley lifetimes. Our work reveals several upconversion processes that bring out an intricate coupling network linking spin-0 and spin-1 intra- and inter-valley excitons, demonstrating that intervalley scattering and spin-flip are very common processes in the atomic layer. These experimental results deepen our understanding of tungsten diselenide exciton physics and illustrate that planar photonic devices are capable of harnessing versatile exciton species in TMD semiconductors.

preprint2020arXiv

Weak KAM solutions of Hamilton-Jacobi equations with decreasing dependence on unknown functions

We consider the Hamilton-Jacobi equation \[{H}(x,u,Du)=0,\quad x\in M, \] where $M$ is a connected, closed and smooth Riemannian manifold, ${H}(x,u,p)$ satisfies Tonelli conditions with respect to $p$ and certain decreasing condition with respect to $u$. Based on a dynamical approach developed in \cite{WWY,WWY1,WWY2}, we obtain a series of properties for weak KAM solutions (equivalently, viscosity solutions) of the stationary equation and the long time behavior of viscosity solutions of the evolutionary equation on the Cauchy problem \begin{equation*} \begin{cases} w_t+{H}(x,w,w_x)=0,\quad (x,t)\in M\times (0,+\infty),\\ w(x,0)=φ(x), \quad x\in M. \end{cases} \end{equation*}

preprint2018arXiv

On estimation for Brownian motion governed by telegraph process with multiple off states

Brownian motion whose infinitesimal variance changes according to a three-state continuous time Markov Chain is studied. This Markov Chain can be viewed as a telegraph process with one on state and two off states. We first derive the distribution of occupation time of the on state. Then the result is used to develop a likelihood estimation procedure when the stochastic process at hand is observed at discrete, possibly irregularly spaced time points. The likelihood function is evaluated with the forward algorithm in the general framework of hidden Markov models. The analytic results are confirmed with simulation studies. The estimation procedure is applied to analyze the position data from a mountain lion.

preprint2016arXiv

Activation of new Raman modes by inversion symmetry breaking in type II Weyl semimetal candidate T'-MoTe2

We synthesized distorted octahedral(T') molybdenum ditelluride (MoTe2) and investigated its vibrational properties with Raman spectroscopy, density functional theory and symmetry analysis. Compared to the results from high temperature centrosymmetric monoclinic (T'mo) phase, four new Raman bands emerge in the low temperature orthorhombic (T'or) phase, which was recently predicted to be a type II Weyl semimetal. Crystal-angle-dependent, light-polarization-resolved measurements indicate that all the observed Raman peaks belong to two categories: those vibrating along the zigzag Mo atomic chain (z-mode) and those vibrating in the mirror plane (m-mode) perpendicular to the zigzag chain. Interestingly the low energy shear z-mode and shear m-mode, absent from the T'mo spectra, become activated when sample cooling induces a phase transition to the T'or crystal structure. We interpret this observation as a consequence of inversion-symmetry breaking, which is crucial for the existence of Weyl fermions in the layered crystal. Our temperature dependent Raman measurements further show that both the high energy m-mode at 130 cm-1 and the low energy shear m-mode at 12 cm-1 provide useful gauges for monitoring the broken inversion symmetry in the crystal.

preprint2016arXiv

Detection of False Data Injection Attacks in Smart Grid under Colored Gaussian Noise

In this paper, we consider the problems of state estimation and false data injection detection in smart grid when the measurements are corrupted by colored Gaussian noise. By modeling the noise with the autoregressive process, we estimate the state of the power transmission networks and develop a generalized likelihood ratio test (GLRT) detector for the detection of false data injection attacks. We show that the conventional approach with the assumption of Gaussian noise is a special case of the proposed method, and thus the new approach has more applicability. {The proposed detector is also tested on an independent component analysis (ICA) based unobservable false data attack scheme that utilizes similar assumptions of sample observation.} We evaluate the performance of the proposed state estimator and attack detector on the IEEE 30-bus power system with comparison to conventional Gaussian noise based detector. The superior performance of {both observable and unobservable false data attacks} demonstrates the effectiveness of the proposed approach and indicates a wide application on the power signal processing.

preprint2016arXiv

Extended relativistic configuration interaction and many-body perturbation calculations of spectroscopic data for the $n \leq 6$ configurationsin ne-like ions between Cr XV and Kr XXVII

Level energies, wavelengths, electric dipole, magnetic dipole, electric quadrupole, and magnetic quadrupole transition rates, oscillator strengths, and line strengths from combined relativistic configuration interaction and many-body perturbation calculations are reported for the 201 fine-structure states of the $2s^2 2p^6$, $2s^2 2p^5 3l$, $2s 2p^6 3l$, $2s^2 2p^5 4l$, $2s 2p^6 4l$, $2s^2 2p^5 5l$, and $2s^2 2p^5 6l$ configurations in all Ne-like ions between Cr XV and Kr XXVII. Calculated level energies and transition data are compared with experiments from the NIST and CHIANTI databases, and other recent benchmark calculations. The mean energy difference with the NIST experiments is only 0.05%. The present calculations significantly increase the amount of accurate spectroscopic data for the $n >3$ states in a number of Ne-like ions of astrophysics interest. A complete dataset should be helpful in analyzing new observations from the solar and other astrophysical sources, and is also likely to be useful for modeling and diagnosing a variety of plasmas including astronomical and fusion plasma.

preprint2016arXiv

Intrinsic Phonon Bands in High-Quality Monolayer T' Molybdenum Ditelluride

The topologically nontrivial and chemically functional distorted octahedral (T') transition metal dichalcogenides (TMDCs) are a type of layered semimetal that has attracted significant recent attention. However, the properties of monolayer (1L) T'-TMDC, a fundamental unit of the system, is still largely unknown due to rapid sample degradation in air. Here we report that well-protected 1L CVD T'-MoTe2 exhibits sharp and robust intrinsic Raman bands, with intensities about one order of magnitude stronger than those from bulk T'-MoTe2. The high-quality samples enabled us to reveal the set of all nine even-parity zone-center optical phonons, providing reliable fingerprints for the previously elusive crystal. By performing light polarization and crystal orientation resolved scattering analysis, we can effectively distinguish the intrinsic modes from Te-metalloid-like modes A (~122 cm-1) and B (~141 cm-1) which are related to the sample degradation. Our studies offer a powerful non-destructive method for assessing sample quality and for monitoring sample degradation in situ, representing a solid advance in understanding the fundamental properties of 1L-T'-TMDCs.

preprint2016arXiv

Inversion-symmetry-breaking-activated shear Raman bands in $T'$-MoTe$_2$

Type-II Weyl fermion nodes, located at the touching points between electron and hole pockets, have been recently predicted to occur in distorted octahedral ($T'$) transition metal dichalcogenide semimetals, contingent upon the condition that the layered crystal has the noncentrosymmetric orthorhombic ($T'_{or}$) stacking. Here, we report on the emergence of two shear Raman bands activated by inversion symmetry breaking in $T'$-MoTe$_2$ due to sample cooling. Polarization and crystal orientation resolved measurements further point to a phase transition from the monoclinic ($T'_{mo}$) structure to the desired $T'_{or}$ lattice. These results provide spectroscopic evidence that low-temperature $T'$-MoTe$_2$ is suitable for probing type-II Weyl physics.

preprint2016arXiv

Raman scattering and anomalous Stokes anti-Stokes ratio in MoTe2 atomic layers

Stokes and anti-Stokes Raman scattering are performed on atomic layers of hexagonal molybdenum ditelluride (MoTe2), a prototypical transition metal dichalcogenide (TMDC) semiconductor. The data reveal all six types of zone center optical phonons, along with their corresponding Davydov splittings, which have been challenging to see in other TMDCs. We discover that the anti-Stokes Raman intensity of the low energy layer-breathing mode becomes more intense than the Stokes peak under certain experimental conditions, and find the effect to be tunable by excitation frequency and number of atomic layers. These observations are interpreted as a result of resonance effects arising from the C excitons in the vicinity of the Brillouin zone center in the photon-electron-phonon interaction process.

preprint2015arXiv

Antenna enhanced graphene THz emitter and detector

Recent intense electrical and optical studies of graphene have pushed the material to the forefront of optoelectronic research. Of particular interest is the few terahertz (THz) frequency regime where efficient light sources and highly sensitive detectors are very challenging to make. Here we present THz sources and detectors made with graphene field effect transistors (GFETs) enhanced by a double-patch antenna and an on-chip silicon lens. We report the first experimental observation of 1-3 THz radiation from graphene, as well as four orders of magnitude performance improvements in a GFET thermoelectric detector operating at ~2 THz. The quantitative analysis of the emitting power and its unusual charge density dependence indicate significant non-thermal contribution from the GFET. The polarization resolved detection measurements with different illumination geometries allow for detailed and quantitative analysis of various factors that contribute to the overall detector performance. Our experimental results represent a significant advance towards practically useful graphene THz devices.

preprint2015arXiv

Eclipsing Binaries From the CSTAR Project at Dome A, Antarctica

The Chinese Small Telescope ARray (CSTAR) has observed an area around the Celestial South Pole at Dome A since 2008. About $20,000$ light curves in the i band were obtained lasting from March to July, 2008. The photometric precision achieves about 4 mmag at i = 7.5 and 20 mmag at i = 12 within a 30 s exposure time. These light curves are analyzed using Lomb--Scargle, Phase Dispersion Minimization, and Box Least Squares methods to search for periodic signals. False positives may appear as a variable signature caused by contaminating stars and the observation mode of CSTAR. Therefore the period and position of each variable candidate are checked to eliminate false positives. Eclipsing binaries are removed by visual inspection, frequency spectrum analysis and locally linear embedding technique. We identify 53 eclipsing binaries in the field of view of CSTAR, containing 24 detached binaries, 8 semi-detached binaries, 18 contact binaries, and 3 ellipsoidal variables. To derive the parameters of these binaries, we use the Eclipsing Binaries via Artificial Intelligence (EBAI) method. The primary and the secondary eclipse timing variations (ETVs) for semi-detached and contact systems are analyzed. Correlated primary and secondary ETVs confirmed by false alarm tests may indicate an unseen perturbing companion. Through ETV analysis, we identify two triple systems (CSTAR J084612.64-883342.9 and CSTAR J220502.55-895206.7). The orbital parameters of the third body in CSTAR J220502.55-895206.7 are derived using a simple dynamical model.

preprint2015arXiv

Helicity resolved Raman scattering of MoS2, MoSe2, WS2 and WSe2 atomic layers

The two-fold valley degeneracy in two-dimensional (2D) semiconducting transition metal dichalcogenides (TMDCs) (Mo,W)(S,Se)2 is suitable for "valleytronics", the storage and manupulation of information utilizing the valley degree of freedom. The conservation of luminescent photon helicity in these 2D crystal monolayers has been widely regarded as a benchmark indicator for charge carrier valley polarization. Here we perform helicity-resolved Raman scattering of the TMDC atomic layers. In drastic contrast to luminescence, the dominant first-order zone-center Raman bands, including the low energy breathing and shear modes as well as the higher energy optical phonons, are found to either maintain or completely switch the helicity of incident photons. These experimental observations, in addition to providing a useful tool for characterization of TMDC atomic layers, shed new light on the connection between photon helicity and valley polarization.

preprint2015arXiv

LINE: Large-scale Information Network Embedding

This paper studies the problem of embedding very large information networks into low-dimensional vector spaces, which is useful in many tasks such as visualization, node classification, and link prediction. Most existing graph embedding methods do not scale for real world information networks which usually contain millions of nodes. In this paper, we propose a novel network embedding method called the "LINE," which is suitable for arbitrary types of information networks: undirected, directed, and/or weighted. The method optimizes a carefully designed objective function that preserves both the local and global network structures. An edge-sampling algorithm is proposed that addresses the limitation of the classical stochastic gradient descent and improves both the effectiveness and the efficiency of the inference. Empirical experiments prove the effectiveness of the LINE on a variety of real-world information networks, including language networks, social networks, and citation networks. The algorithm is very efficient, which is able to learn the embedding of a network with millions of vertices and billions of edges in a few hours on a typical single machine. The source code of the LINE is available online.

preprint2015arXiv

Photometric Variability in the CSTAR Field: Results From the 2008 Data Set

The Chinese Small Telescope ARray (CSTAR) is the first telescope facility built at Dome A, Antarctica. During the 2008 observing season, the installation provided long-baseline and high-cadence photometric observations in the i-band for 18,145 targets within 20 deg2 CSTAR field around the South Celestial Pole for the purpose of monitoring the astronomical observing quality of Dome A and detecting various types of photometric variability. Using sensitive and robust detection methods, we discover 274 potential variables from this data set, 83 of which are new discoveries. We characterize most of them, providing the periods, amplitudes and classes of variability. The catalog of all these variables is presented along with the discussion of their statistical properties.

preprint2014arXiv

A class of goodness-of-fit tests for spatial extremes models based on max-stable processes

Parametric max-stable processes are increasingly used to model spatial extremes. Starting from the fact that the dependence structure of a max-stable process is completely characterized by an extreme-value copula, a class of goodness-of-fit tests is proposed based on the comparison between a nonparametric and a parametric estimator of the corresponding unknown multivariate Pickands dependence function. Because of the high-dimensional setting under consideration, these functional estimators are only compared at a specific set of points at which they coincide, up to a multiplicative constant, with estimators of the extremal coefficients. The nonparametric estimators of the Pickands dependence function used in this work are those recently studied by Gudendorf and Segers. The parametric estimators rely on the use of the {\em pairwise pseudo-likelihood} which extends the concept of pairwise (composite) likelihood to a rank-based context. Approximate $p$-values for the resulting margin-free tests are obtained by means of a {\em one- or two-level parametric bootstrap}. Conditions for the asymptotic validity of these resampling procedures are given based on the work of Genest and Rémillard. The finite-sample performance of the tests is investigated in dimension 10 under the Smith, Schlather and geometric Gaussian models. An application of the tests to rainfall data is finally presented.

preprint2014arXiv

A Dynamical Approach to Viscosity Solutions of Hamilton-Jacobi Equations

In this paper, we consider the following Hamilton-Jacobi equation with initial condition: \begin{equation*} \begin{cases} \partial_tu(x,t)+H(x,t,u(x,t),\partial_xu(x,t))=0, u(x,0)=ϕ(x). \end{cases} \end{equation*} Under some assumptions on the convexity of $H(x,t,u,p)$ w.r.t. $p$, we develop a dynamical approach to viscosity solutions and show that there exists an intrinsic connection between viscosity solutions and certain minimal characteristics.

preprint2014arXiv

A Novel Compressed Sensing Based Model for Reconstructing Sparse Signals Using Phase Sparse Character

Phase modulation is a commonly used modulation mode in digital communication, which usually brings phase sparsity to digital signals. It is naturally to connect the sparsity with the newly emerged theory of compressed sensing (CS), which enables sub-Nyquist sampling of high-bandwidth to sparse signals. For the present, applications of CS theory in communication field mainly focus on spectrum sensing, sparse channel estimation etc. Few of current researches take the phase sparse character into consideration. In this paper, we establish the novel model of phase modulation signals based on phase sparsity, and introduce CS theory to the phase domain. According to CS theory, rather than the bandwidth, the sampling rate required here is scaling with the symbol rate, which is usually much lower than the Nyquist rate. In this paper, we provide analytical support for the model, and simulations verify its validity.

preprint2014arXiv

Automatic Modulation Recognition of PSK Signals Using Nonuniform Compressive Samples Based on High Order Statistics

Phase modulation is a commonly used modulation mode in digital communication, which usually brings phase sparsity to digital signals. It is naturally to connect the sparsity with the newly emerged theory of compressed sensing (CS), which enables sub-Nyquist sampling of high-bandwidth to sparse signals. For the present, applications of CS theory in communication field mainly focus on spectrum sensing, sparse channel estimation etc. Few of current researches take the phase sparse character into consideration. In this paper, we establish the novel model of phase modulation signals based on phase sparsity, and introduce CS theory to the phase domain. According to CS theory, rather than the bandwidth, the sampling rate required here is scaling with the symbol rate, which is usually much lower than the Nyquist rate. In this paper, we provide analytical support for the model, and simulations verify its validity.

preprint2014arXiv

Automatic Modulation Recognition of PSK Signals with Sub-Nyquist Sampling Based on High Order Statistics

Sampling rate required in the Nth Power Nonlinear Transformation (NPT) method is typically much greater than Nyquist rate, which causes heavy burden for the Analog to Digital Converter (ADC). Taking advantage of the sparse property of PSK signals' spectrum under NPT, we develop the NPT method for PSK signals with Sub-Nyquist rate samples. In this paper, combined the NPT method with Compressive Sensing (CS) theory, frequency spectrum reconstruction of the Nth power nonlinear transformation of PSK signals is presented, which can be further used for AMR and rough estimations of unknown carrier frequency and symbol rate.

preprint2014arXiv

Planetary Transit Candidates in the CSTAR Field: Analysis of the 2008 Data

The Chinese Small Telescope ARray (CSTAR) is a group of four identical, fully automated, static 14.5 cm telescopes. CSTAR is located at Dome A, Antarctica and covers 20 square degree of sky around the South Celestial Pole. The installation is designed to provide high-cadence photometry for the purpose of monitoring the quality of the astronomical observing conditions at Dome A and detecting transiting exoplanets. CSTAR has been operational since 2008, and has taken a rich and high-precision photometric data set of 10,690 stars. In the first observing season, we obtained 291,911 qualified science frames with 20-second integrations in the i-band. Photometric precision reaches about 4 mmag at 20-second cadence at i=7.5, and is about 20 mmag at i=12. Using robust detection methods, ten promising exoplanet candidates were found. Four of these were found to be giants using spectroscopic follow-up. All of these transit candidates are presented here along with the discussion of their detailed properties as well as the follow-up observations.

preprint2014arXiv

Sensitive Room-Temperature Terahertz Detection via Photothermoelectric Effect in Graphene

Terahertz (THz) radiation has uses from security to medicine; however, sensitive room-temperature detection of THz is notoriously difficult. The hot-electron photothermoelectric effect in graphene is a promising detection mechanism: photoexcited carriers rapidly thermalize due to strong electron-electron interactions, but lose energy to the lattice more slowly. The electron temperature gradient drives electron diffusion, and asymmetry due to local gating or dissimilar contact metals produces a net current via the thermoelectric effect. Here we demonstrate a graphene thermoelectric THz photodetector with sensitivity exceeding 10 V/W (700 V/W) at room temperature and noise equivalent power less than 1100 pW/Hz^1/2 (20 pW/Hz^1/2), referenced to the incident (absorbed) power. This implies a performance which is competitive with the best room-temperature THz detectors for an optimally coupled device, while time-resolved measurements indicate that our graphene detector is eight to nine orders of magnitude faster than those. A simple model of the response, including contact asymmetries (resistance, work function and Fermi-energy pinning) reproduces the qualitative features of the data, and indicates that orders-of-magnitude sensitivity improvements are possible.

preprint2014arXiv

Weak KAM theory for general Hamilton-Jacobi equations I: the solution semigroup under proper conditions

We consider the following evolutionary Hamilton-Jacobi equation with initial condition: \begin{equation*} \begin{cases} \partial_tu(x,t)+H(x,u(x,t),\partial_xu(x,t))=0,\\ u(x,0)=ϕ(x). \end{cases} \end{equation*} Under some assumptions on $H(x,u,p)$ with respect to $p$ and $u$, we provide a variational principle on the evolutionary Hamilton-Jacobi equation. By introducing an implicitly defined solution semigroup, we extend Fathi's weak KAM theory to certain more general cases, in which $H$ explicitly depends on the unknown function $u$. As an application, we show the viscosity solution of the evolutionary Hamilton-Jacobi equation with initial condition tends asymptotically to the weak KAM solution of the following stationary Hamilton-Jacobi equation: \begin{equation*} H(x,u(x),\partial_xu(x))=0. \end{equation*}.

preprint2014arXiv

Weak KAM theory for general Hamilton-Jacobi equations II: the fundamental solution under Lipschitz conditions

We consider the following evolutionary Hamilton-Jacobi equation with initial condition: \begin{equation*} \begin{cases} \partial_tu(x,t)+H(x,u(x,t),\partial_xu(x,t))=0,\\ u(x,0)=ϕ(x), \end{cases} \end{equation*} where $ϕ(x)\in C(M,\mathbb{R})$. Under some assumptions on the convexity of $H(x,u,p)$ with respect to $p$ and the uniform Lipschitz of $H(x,u,p)$ with respect to $u$, we establish a variational principle and provide an intrinsic relation between viscosity solutions and certain minimal characteristics. By introducing an implicitly defined {\it fundamental solution}, we obtain a variational representation formula of the viscosity solution of the evolutionary Hamilton-Jacobi equation. Moreover, we discuss the large time behavior of the viscosity solution of the evolutionary Hamilton-Jacobi equation and provide a dynamical representation formula of the viscosity solution of the stationary Hamilton-Jacobi equation with strictly increasing $H(x,u,p)$ with respect to $u$.

preprint2014arXiv

Weak KAM theory for general Hamilton-Jacobi equations III: the variational principle under Osgood conditions

We consider the following evolutionary Hamilton-Jacobi equation with initial condition: \begin{equation*} \begin{cases} \partial_tu(x,t)+H(x,u(x,t),\partial_xu(x,t))=0,\\ u(x,0)=ϕ(x), \end{cases} \end{equation*} where $ϕ(x)\in C(M,\mathbb{R})$. Under some assumptions on the convexity of $H(x,u,p)$ with respect to $p$ and the Osgood growth of $H(x,u,p)$ with respect to $u$, we establish an implicitly variational principle and provide an intrinsic relation between viscosity solutions and certain minimal characteristics. Moreover, we obtain a representation formula of the viscosity solution of the evolutionary Hamilton-Jacobi equation.

preprint2013arXiv

Anomalous melting behavior of solid hydrogen at high pressures

Hydrogen is the most abundant element in the universe, and its properties under conditions of high temperature and pressure are crucial to understand the interior of of large gaseous planets and other astrophysical bodies. At ultra high pressures solid hydrogen has been predicted to transform into a quantum fluid, because of its high zero point motion. Here we report first principles two phase coexistence and Z method determinations of the melting line of solid hydrogen in a pressure range spanning from 30 to 600 GPa. Our results suggest that the melting line of solid hydrogen, as derived from classical molecular dynamics simulations, reaches a minimum of 367 K at about 430 GPa, at higher pressures the melting line of the atomics Cs IV phase regain a positive slope. In view of the possible importance of quantum effects in hydrogen at such low temperatures, we also determined the melting temperature of the atomic CsIV phase at pressures of 400, 500, 600 GPa, employing Feynman path integral simulations. These result in a downward shift of the classical melting line by about 100 K, and hint at a possible secondary maximum in the melting line in the region between 500 and 600 GPa, testifying to the importance of quantum effects in this system. Combined, our results imply that the stability field of the zero temperature quantum liquid phase, if it exists at all, would only occur at higher pressures than previously thought.

preprint2013arXiv

Graphics Processing Unit acceleration of the Random Phase Approximation in the projector augmented wave method

The Random Phase Approximation (RPA) for correlation energy in the grid-based projector augmented wave (gpaw) code is accelerated by porting to the Graphics Processing Unit (GPU) architecture. The acceleration is achieved by grouping independent vectors/matrices and transforming the implementation from being memory bound to being computation/latency bound. With this approach, both the CPU and GPU implementations have been enhanced. We tested the GPU implementation on a few representative systems: molecules (O2), bulk solids (Li2O and MoO3) and molecules adsorbed on metal surfaces (N2/Ru(0001) and CO/Ni(111)). Improvements from 10+ to 40+ have been achieved (8-GPUs versus 8-CPUs). A realistic RPA calculation for CO/Ni(111) surface can be finished in 5.5 h using 8 GPUs. It is thus promising to employ non-self-consistent RPA for routine surface chemistry simulations.

preprint2013arXiv

Quantum molecular dynamics simulations of the thermophysical properties of shocked liquid ammonia for pressures up to 1.3 TPa

We investigate via quantum molecular-dynamics simulations the thermophysical properties of shocked liquid ammonia up to the pressure 1.3 TPa and temperature 120000 K. The principal Hugoniot is predicted from wide-range equation of state, which agrees well with available experimental measurements up to 64 GPa. Our systematic study of the structural properties demonstrates that liquid ammonia undergoes a gradual phase transition along the Hugoniot. At about 4800 K, the system transforms into a metallic, complex mixture state consisting of $\textnormal{N}\textnormal{H}_{3}$, $\textnormal{N}_{2}$, $\textnormal{H}_{2}$, N, and H. Furthermore, we discuss the implications for the interiors of Uranus and Neptune.

preprint2013arXiv

Weak KAM theorem for Hamilton-Jacobi equations

In this paper, we generalize weak KAM theorem from positive Lagrangian systems to "proper" Hamilton-Jacobi equations. We introduce an implicitly defined solution semigroup of evolutionary Hamilton-Jacobi equations. By exploring the properties of the solution semigroup, we prove the convergence of solution semigroup and existence of weak KAM solutions for stationary equations: \begin{equation*} H(x, u, d_x u)=0. \end{equation*}

preprint2012arXiv

Conventional and Acoustic Surface Plasmons on Noble Metal Surfaces: A Time-dependent Density Functional Theory Study

First-principles calculations of the conventional and acoustic surface plasmons (CSPs and ASPs) on the (111) surfaces of Cu, Ag, and Au are presented. The effect of $s-d$ interband transitions on both types of plasmons is investigated by comparing results from the local density approximation and an orbital dependent exchange-correlation (xc) potential that improves the position and width of the $d$ bands. The plasmon dispersions calculated with the latter xc-potential agree well with electron energy loss spectroscopy (EELS) experiments. For both the CSP and ASP, the same trend of Cu$<$Au$<$Ag is found for the plasmon energies and is attributed to the reduced screening by interband transitions from Cu, to Au and Ag. This trend for the ASP, however, contradicts a previous model prediction. While the ASP is seen as a weak feature in the EELS, it can be clearly identified in the static and dynamic dielectric band structure.

preprint2012arXiv

Goodness-of-fit testing based on a weighted bootstrap: A fast large-sample alternative to the parametric bootstrap

The process comparing the empirical cumulative distribution function of the sample with a parametric estimate of the cumulative distribution function is known as the empirical process with estimated parameters and has been extensively employed in the literature for goodness-of-fit testing. The simplest way to carry out such goodness-of-fit tests, especially in a multivariate setting, is to use a parametric bootstrap. Although very easy to implement, the parametric bootstrap can become very computationally expensive as the sample size, the number of parameters, or the dimension of the data increase. An alternative resampling technique based on a fast weighted bootstrap is proposed in this paper, and is studied both theoretically and empirically. The outcome of this work is a generic and computationally efficient multiplier goodness-of-fit procedure that can be used as a large-sample alternative to the parametric bootstrap. In order to approximately determine how large the sample size needs to be for the parametric and weighted bootstraps to have roughly equivalent powers, extensive Monte Carlo experiments are carried out in dimension one, two and three, and for models containing up to nine parameters. The computational gains resulting from the use of the proposed multiplier goodness-of-fit procedure are illustrated on trivariate financial data. A by-product of this work is a fast large-sample goodness-of-fit procedure for the bivariate and trivariate t distribution whose degrees of freedom are fixed.

preprint2012arXiv

Optical properties of bulk semiconductors and graphene/boron-nitride: The Bethe-Salpeter equation with derivative discontinuity-corrected DFT energies

We present an efficient implementation of the Bethe-Salpeter equation (BSE) for optical properties of materials in the projector augmented wave method GPAW. Single-particle energies and wave functions are obtained from the GLLBSC functional which explicitly includes the derivative discontinuity, is computationally inexpensive, and yields excellent fundamental gaps. Electron-hole interactions are included through the BSE using the statically screened interaction evaluated in the random phase approximation. For a representative set of semiconductors and insulators we find excellent agreement with experiments for the dielectric functions, onset of absorption, and lowest excitonic features. For the two-dimensional systems of graphene and hexagonal boron-nitride (h-BN) we find good agreement with previous many-body calculations. For the graphene/h-BN interface, we find that the fundamental and optical gaps of the h-BN layer are reduced by 2.0 eV and 0.7 eV, respectively, compared to freestanding h-BN. This reduction is due to image charge screening which shows up in the GLLBSC calculation as a reduction (vanishing) of the derivative discontinuity.

preprint2012arXiv

Raman Spectroscopy of magneto-phonon resonances in Graphene and Graphite

The magneto-phonon resonance or MPR occurs in semiconductor materials when the energy spacing between Landau levels is continuously tuned to cross the energy of an optical phonon mode. MPRs have been largely explored in bulk semiconductors, in two-dimensional systems and in quantum dots. Recently there has been significant interest in the MPR interactions of the Dirac fermion magnetoexcitons in graphene, and a rich splitting and anti-crossing phenomena of the even parity E2g long wavelength optical phonon mode have been theoretically proposed and experimentally observed. The MPR has been found to crucially depend on disorder in the graphene layer. This is a feature that creates new venues for the study of interplays between disorder and interactions in the atomic layers. We review here the fundamentals of MRP in graphene and the experimental Raman scattering works that have led to the observation of these phenomena in graphene and graphite.

preprint2012arXiv

Semiparametric Multivariate Accelerated Failure Time Model with Generalized Estimating Equations

The semiparametric accelerated failure time model is not as widely used as the Cox relative risk model mainly due to computational difficulties. Recent developments in least squares estimation and induced smoothing estimating equations provide promising tools to make the accelerate failure time models more attractive in practice. For semiparametric multivariate accelerated failure time models, we propose a generalized estimating equation approach to account for the multivariate dependence through working correlation structures. The marginal error distributions can be either identical as in sequential event settings or different as in parallel event settings. Some regression coefficients can be shared across margins as needed. The initial estimator is a rank-based estimator with Gehan's weight, but obtained from an induced smoothing approach with computation ease. The resulting estimator is consistent and asymptotically normal, with a variance estimated through a multiplier resampling method. In a simulation study, our estimator was up to three times as efficient as the initial estimator, especially with stronger multivariate dependence and heavier censoring percentage. Two real examples demonstrate the utility of the proposed method.

preprint2011arXiv

A goodness-of-fit test for bivariate extreme-value copulas

It is often reasonable to assume that the dependence structure of a bivariate continuous distribution belongs to the class of extreme-value copulas. The latter are characterized by their Pickands dependence function. In this paper, a procedure is proposed for testing whether this function belongs to a given parametric family. The test is based on a Cramér--von Mises statistic measuring the distance between an estimate of the parametric Pickands dependence function and either one of two nonparametric estimators thereof studied by Genest and Segers [Ann. Statist. 37 (2009) 2990--3022]. As the limiting distribution of the test statistic depends on unknown parameters, it must be estimated via a parametric bootstrap procedure, the validity of which is established. Monte Carlo simulations are used to assess the power of the test and an extension to dependence structures that are left-tail decreasing in both variables is considered.

preprint2011arXiv

Correlated charged impurity scattering in graphene

Understanding disorder in graphene is essential for electronic applications; in contrast to conventional materials, the extraordinarily low electron-phonon scattering1, 2 in graphene implies that disorder3-7 dominates its resistivity even at room temperature. Charged impurities5, 8-10 have been identified as an important disorder type in graphene on SiO2 substrates11, 12, giving a nearly linear carrier-density-dependent conductivity σ(n), and producing electron and hole puddles13-15 which determine the magnitude of graphene's minimum conductivity σmin10. Correlations of charged impurities are known to be essential in achieving the highest mobilities in remotely-doped semiconductor heterostructures16-18, and are present to some degree in any impurity system at finite temperature. Here we show that even modest correlations in the position of charged impurities, realized by annealing potassium on graphene, can increase the mobility by more than a factor of four. The results are well understood theoretically19 considering an impurity correlation length which is temperature dependent but independent of impurity density. Impurity correlations also naturally explain the sub-linear σ(n) commonly observed in substrate-bound graphene devices2, 11, 12, 20.

preprint2011arXiv

Large-sample tests of extreme-value dependence for multivariate copulas

Starting from the characterization of extreme-value copulas based on max-stability, large-sample tests of extreme-value dependence for multivariate copulas are studied. The two key ingredients of the proposed tests are the empirical copula of the data and a multiplier technique for obtaining approximate p-values for the derived statistics. The asymptotic validity of the multiplier approach is established, and the finite-sample performance of a large number of candidate test statistics is studied through extensive Monte Carlo experiments for data sets of dimension two to five. In the bivariate case, the rejection rates of the best versions of the tests are compared with those of the test of Ghoudi, Khoudraji and Rivest (1998) recently revisited by Ben Ghorbal, Genest and Neslehova (2009). The proposed procedures are illustrated on bivariate financial data and trivariate geological data.

preprint2011arXiv

Linear density response function in the projector-augmented wave method: Applications to solids, surfaces, and interfaces

We present an implementation of the linear density response function within the projector-augmented wave (PAW) method with applications to the linear optical and dielectric properties of both solids, surfaces, and interfaces. The response function is represented in plane waves while the single-particle eigenstates can be expanded on a real space grid or in atomic orbital basis for increased efficiency. The exchange-correlation kernel is treated at the level of the adiabatic local density approximation (ALDA) and crystal local field effects are included. The calculated static and dynamical dielectric functions of Si, C, SiC, AlP and GaAs compare well with previous calculations. While optical properties of semiconductors, in particular excitonic effects, are generally not well described by ALDA, we obtain excellent agreement with experiments for the surface loss function of the Mg(0001) surface with plasmon energies deviating by less than 0.2 eV. Finally, we apply the method to study the influence of substrates on the plasmon excitations in graphene. On SiC(0001), the long wavelength $π$ plasmons are significantly damped although their energies remain almost unaltered. On Al(111) the $π$ plasmon is completely quenched due to the coupling to the metal surface plasmon.

preprint2011arXiv

Quantum molecular dynamics simulations for the nonmetal-metal transition in shocked methane

We have performed quantum molecular-dynamics simulations for methane under shock compressions up to 80 GPa. We obtain good agreement with available experimental data for the principal Hugoniot, derived from the equation of state. A systematic study of the optical conductivity spectra, one-particle density of states, and the distributions of the electronic charge over supercell at Hugoniot points shows that the transition of shocked methane to a metallic state takes place close to the density at which methane dissociates significantly into molecular hydrogen and some long alkane chains. Through analyzing the pair correlation function, we predict the chemical picture of the shocked methane. In contrast to usual assumptions used for high pressure modeling of methane, we find that no diamond-like configurations occurs for the whole density-temperature range studied.

preprint2011arXiv

Rapid collapse of spin waves in non-uniform phases of the second Landau level

The spin degree of freedom in quantum phases of the second Landau level is probed by resonant light scattering. The long wavelength spin wave, which monitors the degree of spin polarization, is at the Zeeman energy in the fully spin-polarized state at $ν$=3. At lower filling factors the intensity of the Zeeman mode collapses indicating loss of polarization. A novel continuum of low-lying excitations emerges that dominates near $ν$=8/3 and $ν$=5/2. Resonant Rayleigh scattering reveals that quantum fluids for $ν<3$ break up into robust domain structures. While the state at $ν$=5/2 is considered to be fully polarized, these results reveal unprecedented roles for spin degrees of freedom.

preprint2011arXiv

The rate of convergence of new Lax-Oleinik type operators for time-periodic positive definite Lagrangian systems

Assume that the Aubry set of the time-periodic positive definite Lagrangian $L$ consists of one hyperbolic 1-periodic orbit. We provide an upper bound estimate of the rate of convergence of the family of new Lax-Oleinik type operators associated with $L$ introduced by the authors in \cite{W-Y}. In addition, we construct an example where the Aubry set of a time-independent positive definite Lagrangian system consists of one hyperbolic periodic orbit and the rate of convergence of the Lax-Oleinik semigroup cannot be better than $O(\frac{1}{t})$.

preprint2011arXiv

The sky brightness and transparency in i-band at Dome A, Antarctica

The i-band observing conditions at Dome A on the Antarctic plateau have been investigated using data acquired during 2008 with the Chinese Small Telescope ARray. The sky brightness, variations in atmospheric transparency, cloud cover, and the presence of aurorae are obtained from these images. The median sky brightness of moonless clear nights is 20.5 mag arcsec^{-2} in the SDSS $i$ band at the South Celestial Pole (which includes a contribution of about 0.06 mag from diffuse Galactic light). The median over all Moon phases in the Antarctic winter is about 19.8 mag arcsec^{-2}. There were no thick clouds in 2008. We model contributions of the Sun and the Moon to the sky background to obtain the relationship between the sky brightness and transparency. Aurorae are identified by comparing the observed sky brightness to the sky brightness expected from this model. About 2% of the images are affected by relatively strong aurorae.

preprint2010arXiv

A new kind of Lax-Oleinik type operator with parameters for time-periodic positive definite Lagrangian systems

In this paper we introduce a new kind of Lax-Oleinik type operator with parameters associated with positive definite Lagrangian systems for both the time-periodic case and the time-independent case. On one hand, the new family of Lax-Oleinik type operators with an arbitrary $u\in C(M,\mathbb{R}^1)$ as initial condition converges to a backward weak KAM solution in the time-periodic case, while it was shown by Fathi and Mather that there is no such convergence of the Lax-Oleinik semigroup. On the other hand, the new family of Lax-Oleinik type operators with an arbitrary $u\in C(M,\mathbb{R}^1)$ as initial condition converges to a backward weak KAM solution faster than the Lax-Oleinik semigroup in the time-independent case.

preprint2010arXiv

Anomalous optical and electronic properties of dense sodium

Based on ab initio density-functional-theory using generalized gradient approximation, we systematically study the optical and electronic properties of the insulating dense sodium phase (Na-hp4) reported recently [Ma \textit{et al.}, Nature \textbf{458}, 182 (2009)]. The structure is found optically anisotropic and transparent to visible light, which can be well interpreted using its electronic band structure and angular moment decomposed density of states. Through the bader analysis of Na-hp4 at different pressures, we conclude that ionicity exists in the structure and becomes stronger with increasing pressure. In addition, the absorption spectra in the energy range from 1.4 to 2.4 eV are compared with recent experimental results and found good agreement. It is found that the deep-lying valence electrons participate in the interband transition.

preprint2010arXiv

Charge transport in dual gated bilayer graphene with Corbino geometry

The resistance of dual-gated bilayer graphene is measured as a function of temperature and gating electric fields in the Corbino geometry which precludes edge transport. The temperature-dependent resistance is quantitatively described by a two channel conductance model including parallel thermal activation and variable range hopping channels, which gives the electric-field-dependent band gap whose magnitude is found to be in good agreement with infrared absorption experiments. Low temperature transport is similar to previous studies of dual-gated bilayer graphene with edges, suggesting that edge transport does not play an important role.

preprint2010arXiv

Quantum molecular dynamics simulations of lithium melting using Z-method

We performed first-principles molecular dynamics calculations for lithium using the projector augmented waves method and the generalized gradient approximation as exchange-correlation energy. The melting curve of lithium was computed using the \textit{Z}-method technique for pressures up to 30 GPa, which agrees well with the experimental and two-phase simulated results. The change of the melting line slope from positive to negative was predicted by the characteristic shape inversion of the \textit{Z} curve at about 8.2 GPa. Through analyzing the static properties, we conclude that no liquid-liquid phase transition accompanies the occurrence of the melting line maximum, which is caused by the higher compressibility of the liquid phase compared to the solid phase. In addition, we systematically studied the dynamic and optical properties of lithium near melting curve at critical superheating and melting temperatures. It was suggested that spectra difference at critical superheating and melting temperature may be able to diagnose the homogeneous melting.

preprint2010arXiv

Testing and Data Reduction of the Chinese Small Telescope Array (CSTAR) for Dome A, Antarctica

The Chinese Small Telescope ARray (hereinafter CSTAR) is the first Chinese astronomical instrument on the Antarctic ice cap. The low temperature and low pressure testing of the data acquisition system was carried out in a laboratory refrigerator and on the 4500m Pamirs high plateau, respectively. The results from the final four nights of test observations demonstrated that CSTAR was ready for operation at Dome A, Antarctica. In this paper we present a description of CSTAR and the performance derived from the test observations.

preprint2010arXiv

The First Release of the CSTAR Point Source Catalog from Dome A, Antarctica

In 2008 January the 24th Chinese expedition team successfully deployed the Chinese Small Telescope ARray (CSTAR) to DomeA, the highest point on the Antarctic plateau. CSTAR consists of four 14.5cm optical telescopes, each with a different filter (g, r, i and open) and has a 4.5degree x 4.5degree field of view (FOV). It operates robotically as part of the Plateau Observatory, PLATO, with each telescope taking an image every 30 seconds throughout the year whenever it is dark. During 2008, CSTAR #1 performed almost flawlessly, acquiring more than 0.3 million i-band images for a total integration time of 1728 hours during 158 days of observations. For each image taken under good sky conditions, more than 10,000 sources down to 16 mag could be detected. We performed aperture photometry on all the sources in the field to create the catalog described herein. Since CSTAR has a fixed pointing centered on the South Celestial Pole (Dec =-90 degree), all the sources within the FOV of CSTAR were monitored continuously for several months. The photometric catalog can be used for studying any variability in these sources, and for the discovery of transient sources such as supernovae, gamma-ray bursts and minor planets.